UHPLC-MS Chemical Fingerprinting and Antioxidant, Enzyme Inhibition, Anti-Inflammatory In Silico and Cytoprotective Activities of Cladonia chlorophaea and C. gracilis (Cladoniaceae) from Antarctica

The lichen species Cladonia chlorophaea and C. gracilis (Cladoniaceae) are widely distributed in the island archipelago of maritime Antarctica and represent a natural resource of scientific interest. In this work, the metabolomic characterization of the ethanolic extracts of these species and the determination of the antioxidant activity, enzymatic inhibition and anti-inflammatory potential of selected compounds on the 5-lipoxygenase enzyme by molecular docking and cytoprotective activity in the SH-SY5Y cell line were carried out. Nineteen compounds were identified by liquid chromatography coupled with quadrupole-time-of-flight mass spectrometry (UHPLC-ESI-QTOF-MS) in each of the species. The contents of phenolic compounds, antioxidant activity, the inhibition of cholinesterases (acetylcholinesterase and butyrylcholinesterase) and digestive enzymes (α-glucosidase and pancreatic lipase) were variable among species, with better results in C. chlorophaea. Molecular docking evidenced significant binding affinities of some compounds for the 5-lipoxygenase enzyme, together with outstanding pharmacokinetic properties. Both extracts were shown to promote cell viability and a reduction in reactive oxygen species production in an H2O2-induced oxidative stress model. This study contributes to the chemical knowledge of the Cladonia species and demonstrates the biological potential for the prevention and promising treatment of central nervous system pathologies, inflammatory disorders and metabolic alterations.


Introduction
Lichens are defined as a mutualistic relationship between a heterotrophic mycobiont and an autotrophic photobiont (algae and/or cyanobacteria) associated with an abundant and structured microbiome [1,2]. Antarctic lichens have aroused deep interest for their abundance, their sensitivity to environmental variables [3], their natural products and the abundance of particular bioactive compounds derived from the mycobiont that are organized in different chemical classes, such as depsides, depsidones, dibenzofurans, xanthones, terpenic derivatives, anthraquinones, etc., and synthesized through the acetylpolymalonate (APP), shikimic acid (SAP) and mevalonic acid (MAP) pathways in the abundance, their sensitivity to environmental variables [3], their natural products and the abundance of particular bioactive compounds derived from the mycobiont that are organized in different chemical classes, such as depsides, depsidones, dibenzofurans, xanthones, terpenic derivatives, anthraquinones, etc., and synthesized through the acetylpolymalonate (APP), shikimic acid (SAP) and mevalonic acid (MAP) pathways in the fungal component with multiple biological activities [4]. Species of the genus Cladonia (Clado-niaceae, Lecanorales) are characterized by an erect thallus consisting of podecia; primary thallus, persistent or absent; apothecia, absent or present and red or brown; and biatorine. They are usually found growing on mosses [5] (Figure 1).
(a) (b) Lichen secondary metabolites are widely studied by robust techniques such as ultrahigh-resolution chromatography (UHPLC) coupled with a mass spectrometer, and more than 1000 bioactive compounds are currently registered [6]. As for the genus Cladonia, there is evidence of the potential of extracts with major compounds, such as atranorin and oxalic, propionic, tartaric, butyric, malonic, lactic, citric, maleic, fumaric, stytic, succinic, norstytic, usnic, lecanoric, orsellinic and furmarprotocetraric acids, that show significant antioxidant activity, as in the case of the species Cladonia foliacea, C. fimbriata, C. digitata, C. verticillaris, C. furcata, C. pyxidata, C. rangiferina, C. uncialis, C. rappi, C. pocillum and C. clathrate [7]; likewise, recent studies show the positive effect of extracts of Antarctic lichens for the enzymatic inhibition of cholinesterases (acetylcholinesterase and butyrylcholineste-rase) that influence the interruption of nerve cell communication [8]. As a complement, this work shows the effect of these natural products on the inhibition of other enzymes involved in metabolic alterations associated with the development of chronic non-communicable diseases.
The study of the effect of extracts has also been enhanced on precursor targets of alterations of inflammatory origin; one of them, arachidonate 5-lipoxygenase (5-LOX), is a soluble monomeric enzyme that contains 673 amino acids and has two domains in its structure. The smaller domain is the N-terminal domain, while the larger domain is the area comprising the catalytic C-terminus with a non-heme iron atom. 5-LOX catalyzes the formation of 5(S)-hydroperoxyeicosatetraenoic acid (5-HPETE) from arachidonic acid (AA) by incorporating molecular oxygen, thus generating the dehydration of 5-HPETE to leukotriene A4 (LTA4); the conversion of this leukotriene (LTA4) causes the synthesis of other leukotrienes that, together with LTA4, are important mediators in the inflammatory response [9]. Therefore, 5-LOX is considered one of the important enzymes in the inflammatory response, it being the main cause of pathologies such as asthma and other inflammatory disorders, such as allergic rhinitis and rheumatoid arthritis [9,10]. On the other hand, recent studies with extracts of lichens such as Dactylina arctica, Nephromopsis stracheyi, Tuckermannopsis americana and Vulpicida pinastri and compounds such as evernic acid confirm their potential as neuroprotective agents in cellular models of oxidative stress by H2O2 and murine models of Parkinson's disease induced by MPTP (1-methyl-4phenyl-1,2,3,6-tetrahydropyridine) [11,12]. Lichen secondary metabolites are widely studied by robust techniques such as ultrahigh-resolution chromatography (UHPLC) coupled with a mass spectrometer, and more than 1000 bioactive compounds are currently registered [6]. As for the genus Cladonia, there is evidence of the potential of extracts with major compounds, such as atranorin and oxalic, propionic, tartaric, butyric, malonic, lactic, citric, maleic, fumaric, stytic, succinic, norstytic, usnic, lecanoric, orsellinic and furmarprotocetraric acids, that show significant antioxidant activity, as in the case of the species Cladonia foliacea, C. fimbriata, C. digitata, C. verticillaris, C. furcata, C. pyxidata, C. rangiferina, C. uncialis, C. rappi, C. pocillum and C. clathrate [7]; likewise, recent studies show the positive effect of extracts of Antarctic lichens for the enzymatic inhibition of cholinesterases (acetylcholinesterase and butyrylcholineste-rase) that influence the interruption of nerve cell communication [8]. As a complement, this work shows the effect of these natural products on the inhibition of other enzymes involved in metabolic alterations associated with the development of chronic non-communicable diseases.
The study of the effect of extracts has also been enhanced on precursor targets of alterations of inflammatory origin; one of them, arachidonate 5-lipoxygenase (5-LOX), is a soluble monomeric enzyme that contains 673 amino acids and has two domains in its structure. The smaller domain is the N-terminal domain, while the larger domain is the area comprising the catalytic C-terminus with a non-heme iron atom. 5-LOX catalyzes the formation of 5(S)-hydroperoxyeicosatetraenoic acid (5-HPETE) from arachidonic acid (AA) by incorporating molecular oxygen, thus generating the dehydration of 5-HPETE to leukotriene A4 (LTA4); the conversion of this leukotriene (LTA4) causes the synthesis of other leukotrienes that, together with LTA4, are important mediators in the inflammatory response [9]. Therefore, 5-LOX is considered one of the important enzymes in the inflammatory response, it being the main cause of pathologies such as asthma and other inflammatory disorders, such as allergic rhinitis and rheumatoid arthritis [9,10]. On the other hand, recent studies with extracts of lichens such as Dactylina arctica, Nephromopsis stracheyi, Tuckermannopsis americana and Vulpicida pinastri and compounds such as evernic acid confirm their potential as neuroprotective agents in cellular models of oxidative stress by H 2 O 2 and murine models of Parkinson's disease induced by MPTP (1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine) [11,12].
In this work, we report the identification of the compounds present in the ethanolic extract of C. chlorophaea and C. gracilis species, their content of phenolic compounds, antioxidant activity and enzymatic inhibition on cholinesterases, α-glucosidase and pancreatic lipase, an in silico analysis of the anti-inflammatory potential on the 5-LOX enzyme and their cytoprotection activity against H 2 O 2 -induced oxidative stress.

Preparation of Ethanolic Extracts
A total of 5 g of each lichen species was macerated with ethanol (three times with 30 mL each) via ultrasound at room temperature. Each extract was filtered, and the solutions were concentrated under reduced pressure at 38 • C to obtain two final gummy extracts.

LC Parameters and MS Parameters
The separation and identification of the compounds present in the lichen extracts were performed on a UHPLC-ESI-QTOF-MS system equipped with UHPLC Ultimate 3000 RS with Chromeleon 6.8 software (Dionex GmbH, Idstein, Germany) and Bruker maXis ESI-QTOF-MS with the software Data Analysis 4.0 (all Bruker Daltonik GmbH, Bremen, Germany). A total of 5 mg of each extract was dissolved in 2 mL of methanol for analysis and filtered with a polytetrafluoroethylene (PTFE) filter, and 10 µL was injected into the equipment. The chromatographic equipment consisted of a quaternary pump, an autosampler, a thermostated column compartment and a photodiode array detector. Elution was performed with a binary gradient system with eluent (A) 0.1% formic acid in the water, eluent (B) 0.1% formic acid in the acetonitrile and the gradient: 1% B isocratic (0-2 min), 1-5% B (2-3 min), 5% B isocratic (3-5 min), 5-10% B (5-8 min), 10-30% B (8-30 min), 30-95% B (319-38 min) and 1% B isocratic (39-50 min). Separation was carried out with a Thermo 5 µm C18 80 Å column (150 mm × 4.6 mm) at a flow rate of 1.0 mL/min. ESI-QTOF-MS experiments were recorded in negative ion mode, and the scan range was between 100 and 1200 m/z. Electrospray ionization (ESI) conditions included a capillary temperature of 200 • C, a capillary voltage of 2.0 kV, a dry gas flow rate of 8 L/min and a nebulizer pressure of 2 bar. The experiments were performed in automatic MS/MS mode. The structural characterization of secondary metabolites was based on HR full MS, fragmentation patterns and comparisons with the literature data.

Total Phenolic (TP) Content
The total phenolic content (TPC) of the extracts was measured by the Folin-Ciocalteu (FC) method [13,14]. Briefly, the test mixture was FC reagent (10%), the extract sample, Na 2 CO 3 (7%), was kept for 90 min in darkness and absorbance was measured at 750 nm in a microplate reader (BioTek Instrument, Inc., Winooski, VT, USA). The results were expressed as mg of gallic acid per gram of dried lichen.
2.6. Antioxidant Activity 2.6.1. Ferric-Reducing Antioxidant Power (FRAP) Assay The FRAP assay was performed following a method previously described in Karadag et al. and Parra et al. [15,16]. Briefly, the assay mixture was FRAP reagent, the extract sample was kept for 30 min in darkness and absorbance was measured at 593 nm in a microplate reader (BioTek Instrument, Inc., Winooski, VT, USA). Trolox was used as the reference compound. The results were expressed as micromoles of Trolox equivalents per gram of dried lichen.

Oxygen Radical Absorbance Capacity (ORAC) Assay
The assay was performed following a method previously described by Huang et al. [17]. Briefly, the assay mixture was fluorescein, Trolox and phosphate buffer, the extract sample was incubated for 30 min at 37 • C and AAPH solution was added as the peroxyl generator. Fluorescence intensity was measured at 520 nm on a microplate reader (BioTek Instrument, Inc., Winooski, VT, USA) every 2 min for 1 h and 30 min at an excitation/emission wavelength of 485/520 nm. Trolox was used as the reference compound. The results were expressed as micromoles of Trolox equivalents per gram of dried lichen.

DPPH Scavenging Activity
The assay was determined as previously reported [18]. Briefly, the assay mixture was a solution of DPPH (400 mM), the standard or extract sample, as appropriate, was kept for 30 min in the dark and absorbance was measured at 515 nm on a microplate reader (BioTek Instrument, Inc., Winooski, VT, USA). Gallic acid was used as a reference standard. The results were expressed as IC 50 values (µg of lichen/mL).

Cholinesterase Inhibition
According to Ellman's method [19], the inhibition activity was determined. DTNB solution (3mM), AChE (0.26 U/mL) and BChE (0.26 U/mL) enzyme solution were added, as appropriate, Tris-HCl buffer (50 mM, pH 8.0) and solutions of the extracts at different concentrations were incubated for 20 min at 25 • C, the substrates acetyl thiocholine iodide (15 mM) and butyryl thiocholine chloride (15 mM) were added, as appropriate, and absorbance was measured at 412 nm on a microplate reader (BioTek Instrument, Inc., Winooski, VT, USA). The results were expressed as IC 50 values (µg of lichen/mL). Galantamine was used as a positive control.

α-Glucosidase Inhibition Assay
α-glucosidase inhibition assay was performed according to the method reported by Costamagna et al. and Burgos-Edwards et al. [20,21]. Briefly, the reaction mixture contained sodium phosphate buffer (200 mM, pH 6.9), solutions of the extracts at different concentrations prepared in the same buffer and α-glucosidase (0.1 U/L). After 15 min of preincubation at 37 • C, the reaction was started by adding p-nitrophenyl-α-D-glucopyranoside (5 mM) into the wells and incubated for 30 min at 37 • C, and absorbance was measured at 415 nm in a microplate reader (BioTek Instrument, Inc., Winooski, VT, USA). The results were expressed as IC 50 values (µg of lichen/mL). Acarbose was used as a positive control.

Lipase Inhibition Assay
This assay was performed according to the method reported by McDougall et al. and Picot et al. [22,23]. Porcine pancreatic lipase type II was resuspended in ultrapure water at 20 mg/mL. The substrate 4-nitrophenyl-dodecanoate was prepared with DMSO and ethanol. Briefly, the assay mixture was Tris-HCl buffer (0.1 M, pH 8.5), lipase (10 mg/mL) and solutions of the extracts at different concentrations, the substrate solution (5 mM) was incubated for 15 min at 37 • C and absorbance was measured at 410 nm in a microplate reader (BioTek Instrument, Inc., Winooski, VT, USA). The results were expressed as IC 50 values (µg of lichen/mL). Orlistat was used as a positive control.

Calculation of ADME Parameters
Osiris Data Warrior toolkits (v 5.5.0) were used to verify the pharmacokinetic properties of compounds present in extracts of C. gracilis and C. chlorophaea species (usnic acid, thamnolic acid, fumarprotocetraric acid (FP acid), orsellinic acid, atranorin, cetraric acid, squamatic acid and methylorsellinate), compounds reported for the genus Cladonia (13-(beta-D-glucosyloxy)docosanoic acid (13D acid)), squaric acid, perlatolic acid and psoromic acid) and the commercial inhibitor of 5-LOX (Zileuton). The molecular descriptors that were calculated are the logarithm of the partition coefficient (cLogP), the number of hydrogen bond donors, the number of hydrogen bond acceptors, the molecular mass of compounds, the topological polar surface area (TPSA), the number of rotable bonds and violations of Lipinski's rule of five. The absorption percentage (% ABS) was calculated using the following equation [24][25][26]:

Calculation of Risk Toxicity
To calculate the toxicological properties of compounds present in extracts of C. gracilis and C. chlorophaea species (usnic acid, thamnolic acid, FP acid, orsellinic acid, atranorin, cetraric acid, squamatic acid and methylorsellinate), compounds reported for the genus Cladonia (13D acid, squaric acid, perlatolic acid and psoromic acid) and the commercial inhibitor of 5-LOX (Zileuton), the Osiris Data Warrior computational tool was used. The toxicity risks evaluated in each of the molecules were mutagenicity, tumorigenicity, irritation and reproductive effects [25,26].

In Silico Analysis
The crystallographic structure of the enzyme arachidonate 5-lipoxygenase (5-LOX) was obtained from the RCSB PDB (Research Collaboratory for Structural Bioinformatics Protein Data Bank) with the PDB ID code 6N2W. This crystallographic structure is crystallized with the nordihydroguaiaretic acid (NDGA) ligand which was used as a reference to carry out a directed coupling in the catalytic site of the NDGA/6N2W complex [9,25].
Enzyme optimization was carried out using UCSF Chimera software (v1.16, San Francisco, California, USA), where water molecules and ligands were removed from the active site of the crystallographic protein. In the same way, all polar hydrogen atoms were added at pH = 7.4. The appropriate ionization states for basic and acidic amino acid residues were considered [25,27].

Docking Simulation
The standard procedure for molecular docking was developed using the rigid crystallographic structure of the enzyme arachidonate 5-lipoxygenase (6N2W), nine flexible ligands that did not present any risk of toxicity (thamnolic acid, 13D acid, perlatolic acid, psoromic acid, orsellinic acid, atranorin, cetraric acid, squamatic acid and methylorsel-linate) and a commercial 5-LOX inhibitor ligand (zileuton) whose torsion angles were identified (during 10 independent runs per ligand). Targeted molecular docking was performed using the UCSF Chimera program [25,27,28], where the reference inhibitor catalytic pocket (NDGA) for arachidonate 5-lipoxygenase (5-LOX) was used. Polar hydrogens and Gasteiger partial charges were added, and a bounding box was created, whose size was fixed to a cube with sides that were 20 Å in length, using the Autodock Vina tools at UCSF Chimera. The centroid of the selected residue was chosen based on the enzyme's putative catalytic site (5-LOX), considering its known catalytic amino acids: His367, His372, His550 and Leu673. The coupling and analysis results were visualized using the Discovery Studio Visualizer [25,26]. After docking was complete, the best conformation for hydrogen bonding or π interactions was analyzed, including the free ligand binding energy (delta G, kcal/mol) [25,26,29].

Cell Treatments
A stock solution (1 mg/mL) of each lichen extract was dissolved in DMSO and phosphate-buffered saline (PBS). SH-SY5Y cells were pretreated with different concentrations of the extracts for 24 h, followed by H 2 O 2 (1 mM, 1 h). The highest concentration treatment had a final DMSO concentration of less than 0.1%.

Metabolic Activity Measurement
Cell survival and the cytoprotection rate were determined according to the MTT colorimetric assay [30], with some modifications. After exposure to the extract treatments, MTT stock solution (2 mg/mL) was added to the wells and incubated for 1 h; subsequently, DMSO (100 µL) was added, and absorbance was measured at 550 nm with a Spectrostar BMG microplate reader.

Intracellular ROS Production
According to the DCFH-DA assay [31] the intracellular production of ROS was determined. In the 96-well plate, DCFH-DA (24 µL) was added in the DMEM medium (1%) without phenol red (12 mL) for 30 min, and then the medium was aspirated. The most effective concentrations of the lichen extracts were added; after 1 h, hydrogen peroxide (250 mM) was added. Fluorescence was measured with a microplate reader (FLUOstar OPTIMA, BMG Labtech, Ortenberg, Germany) at an excitation/emission wavelength of 485/528 nm.

Statistical Analysis
Data were performed in triplicate and expressed as the mean data with standard deviations (SD) using Microsoft Excel 2019 software (Microsoft Office, Microsoft Corporation, Redmond, WA, USA). Statistical significance between groups was set at p < 0.05 and determined by one-way ANOVA with Tukey's post hoc test using the commercial software GraphPad Prism 8.

UHPLC Chromatographic Analysis of Lichens Extracts
The fingerprinting of the ethanolic extracts of C. gracilis and C. chlorophaea were obtained by high-resolution mass spectrometric analysis (UHPLC-MS). The negative mode was used for compound identification. In total, nineteen metabolites were detected in each of the species (Figure 2), and their tentative identification included aromatics, organic acids, lipids, dibenzofurans, depsides and depsidones (Table 1).  Table 1.
The group of compounds identified in the ethanolic extracts of the lichens species C. chlorophaea and C. gracilis form a chemical fingerprint related to the recent metabolomic reports of other Antarctic lichens such as C. metacorallifera with 26 compounds [32], Usnea antarctica with 21 compounds [33], Himantormia lugubris with 28 compounds [8], Lecania brailmontti with 18 compounds, Pseudephebe pubescens with 18 compounds and Sphaerophorus globosus with 14 compounds [26]; some of them, as fumarprotocetraric, lecanoric, squamatic and usnic acids, represent typical compounds in foliose lichens that participate in the biochemical mechanisms for the adaptation of lichens to the habitat [34]. Likewise, the chemical study of extracts is the basis for the detection of new compounds with biological potential, as in the case of the isolation of three p-terphenyls and two fumarprotocetraric acid lactones in the Antarctic lichens Stereocaulon alpinum [35] and C. metacorallifera [32], respectively.

Total Phenolic Contents and Antioxidant Activity
The results corresponding to phenolic content and antioxidant activity measured by the different colorimetric assays are detailed in Table 2. The ethanolic extract of C. chlorophaea showed six-times-higher phenolic content compared to C. gracilis, with 330.276 ± 0.006 and 53.563 ± 0.004 mg gallic acid/100 g of extract, respectively. For the antioxidant activity, in the FRAP and ORAC assays, the C. chlorophaea extract presented better results compared to the C. gracilis extract; as for the DPPH assay, both extracts presented a high IC50, it being higher for the C. chlorophaea extract.
The ethanolic extracts of C. chlorophaea and C. gracilis showed phenolic contents and antioxidant activities similar to those reported in the methanolic, acetonic, chloroformic and hydroalcoholic extracts of other species of Cladonia; however, there are variable cases such as the sample of C. foliaceae, where the methanolic extract did not show activity in the DPPH assay, and the percentage of phenolic compounds was 1.7% [36]. In C. furcata,  Table 1. Aromatic Derivatives For C. gracilis and C. chlorophaea: peak 1 corresponded to Na formiate (C 4 H 2 O 4 ) as the internal standard. Peak 6 was identified as squamatic acid (C 19 H 17 O 9 ), peak 15, with a molecular anion at m/z 355.0454, was identified as a criptostictic acid derivate (C 18 H 11 O 8 ) and peak 16, with a molecular anion at m/z 471.0569 and diagnostic peaks at m/z 375.0668, 355.0385, 167.0300 and 943.0975, was identified as fumarprotocetraric acid (C 22 H 15 O 12 ). In C. gracilis, peak 4 was identified as methylorsellinate (C 9 H 9 O 4 ). In C. chlorophaea, peak 3 was identified as orsellinic acid (C 8  Organic Acids For C. gracilis and C. chlorophaea: peak 2 was identified as citric acid (C6H7O7). Lipids For C. gracilis and C. chlorophaea: peak 7 was identified as hexahidroxioxohexacosanoic acid (C 26 H 49 O 9 ), peak 12, with a molecular anion at m/z 489.3432, was identified as pentahydroxyoxohexacosanoic acid (C 26 H 49 O 8 ), with diagnostic peaks at m/z 403.3001 and 979.6848, peak 13, with an [M-H]ion at m/z 417.3236, was identified as 9,10,12,13-tetrahidroxytricosanoic acid (C 23 H 45 O 6 ), with diagnostic peaks at m/z 235.0538 and 195.0616, and peak 24 was identified as octadeca-9,12,15-trienoic acid (C 18 H 29 O 2 ). In C. gracilis, peak 8 was identified as 9,10,12,13-tetrahydroxyheneicosanoic acid (C 21 H 41 O 6 ), peak 9 was identified as pentahydroxyhexacosanoic acid (C 26 H 51 O 7 ), with diagnostic peaks at m/z 448.3405 and 273.0163, peak 10 was identified as 9,10,12,13-tetrahydroxydocosanoic acid (C 22 H 43 O 6 ) and peak 11 was identified as stearin (C 23 H 45 O 7 ). In C. cholorophaea, peak 21, with a molecular anion at m/z 277.2173, was identified as linolenic acid (C 18 H 29 O 2 ), and peak 25 was detected but not identified. Dibenzofurans For C. gracilis and C. chlorophaea: peak 14 was identified as pseudoplacodiolic acid (C19H19O8), with diagnostic peaks at m/z 343.0807, 259.0598 and 231.0648, and peak 19 was identified as usnic acid (C18H15O7).
Depsidones For C. gracilis and C. chlorophaea: peak 18, with a molecular anion at m/z 401.0593, was identified as constictic acid or siphulellic acid (C19H14O10).
The group of compounds identified in the ethanolic extracts of the lichens species C. chlorophaea and C. gracilis form a chemical fingerprint related to the recent metabolomic reports of other Antarctic lichens such as C. metacorallifera with 26 compounds [32], Usnea antarctica with 21 compounds [33], Himantormia lugubris with 28 compounds [8], Lecania brailmontti with 18 compounds, Pseudephebe pubescens with 18 compounds and Sphaerophorus globosus with 14 compounds [26]; some of them, as fumarprotocetraric, lecanoric, squamatic and usnic acids, represent typical compounds in foliose lichens that participate in the biochemical mechanisms for the adaptation of lichens to the habitat [34]. Likewise, the chemical study of extracts is the basis for the detection of new compounds with biological potential, as in the case of the isolation of three p-terphenyls and two fumarprotocetraric acid lactones in the Antarctic lichens Stereocaulon alpinum [35] and C. metacorallifera [32], respectively.

Total Phenolic Contents and Antioxidant Activity
The results corresponding to phenolic content and antioxidant activity measured by the different colorimetric assays are detailed in Table 2. The ethanolic extract of C. chlorophaea showed six-times-higher phenolic content compared to C. gracilis, with 330.276 ± 0.006 and 53.563 ± 0.004 mg gallic acid/100 g of extract, respectively. For the antioxidant activity, in the FRAP and ORAC assays, the C. chlorophaea extract presented better results compared to the C. gracilis extract; as for the DPPH assay, both extracts presented a high IC 50 , it being higher for the C. chlorophaea extract.
The ethanolic extracts of C. chlorophaea and C. gracilis showed phenolic contents and antioxidant activities similar to those reported in the methanolic, acetonic, chloroformic and hydroalcoholic extracts of other species of Cladonia; however, there are variable cases such as the sample of C. foliaceae, where the methanolic extract did not show activity in the DPPH assay, and the percentage of phenolic compounds was 1.7% [36]. In C. furcata, an IC 50 >1000 and 471.3 ± 35.4 µg/mL were reported for DPPH in the methanolic and acetonic extracts, respectively [37]; in other works, an IC 50 in DPPH of 967.97 ± 2.35 µg/mL and a phenol content of 28.00 ± 1.065 mg GAE/g [38] were reported, as well as a 44.83% effect on DPPH activity and a phenolic content of 12.95 ± 1.065 µg pyrocatechol/g [39].
In other cases, such as C. pyxidata and C. rangiferina, for the acetone extract, the radical scavenging activity was 461.72 ± 2.93 and 987.64 ± 2.02 µg/mL, respectively, and the phenol content was 35.35 ± 1.013 and 22.00 ± 1.065 mg GAE/g, respectively [38]. In C. amaurocraea and C. mediterranea, the content of phenolic compounds in acetonic, ethanolic and methanolic extracts was considerably low, with values between 0.005 to 0.174 and 0.006 to 0.009 mg GAE/g, respectively, compared to the amount of proteins and carbohydrates [40]. In C. clathrata with the hydroalcoholic extract, a total phenol concentration of 371.8 ± 57.3 mg GAE/g and a DPPH of 69.25 ± 0.65 µg/mL were found [41], similar to C. verticillata, with phenol contents of 101. 37 ± 3.31 and 48.61 ±0.87 mg GAE/g in the acetonic and methanolic extract, respectively, and the DPPH assay, with 33.09 ± 0.84 and 17.25 ± 2.98 µg/mL for the acetonic and methanolic extract, respectively [42].
In other species such as C. rangiformis, the reported phenolic contents are low, with 0.528 ± 0.007 and 0.103 ± 0.012 mg GAE/g in the chloroform and methanol extracts, respectively [43]. In a work with a majority of compounds from 15 species of the genus Cladonia [44], they report, for DPPH radicals, an IC 50 [7]. For specimens of Cladonia sp., a DPPH of 38.5 ± 3.8 µg/mL has been reported [45]. In general, Antarctic lichen extracts show significant antioxidant activities that enhance their adaptation to extreme environmental conditions [46].  In comparison with other Antarctic lichens, the ethanolic extracts of the two Cladonia species performed better than the ethanolic extracts of Himantormia lugubris [8] and Lecania brialmontii, Pseudephebe pubescens and Sphaerophorus globosus [26], confirming the variable antioxidant potential of species from the polar territories.

Prediction of Pharmacokinetic Properties-ADME
The analysis of the pharmacokinetic properties of the compounds present in the species of the genus Cladonia (usnic acid, thamnolic acid, 13D acid, squaric acid, perlatolic acid, FP acid, psoromic acid, orsellinic acid, atranorin, cetraric acid, squamatic acid and methylorsellinate) was calculated with the Osiris Data Warrior computational tool, where the results of these compounds were compared with the pharmacokinetic properties of a commercial 5-LOX inhibitor (Zileuton) ( Table 4). Drugs that are administered orally must meet certain parameters, which are known as the Lipinski rules [25,26]. These rules allow for the evaluation and monitoring of drugs based on biological and pharmacological function. The parameters encompassed by Lipinski's rules are that the drug must have a molecular weight of less than 500 g/mol; cLogP must have a value less than 5; in the molecular structure there must be a maximum of five hydrogen donor sites; the compound must have 10 or fewer hydrogen bond receptor sites; finally, the compound must have less than 10 rotatable bonds in its chemical structure [25,26]. The molecules that had more than one violation of the Lipinski rules were the compound 13D acid (three violations) and perlatolic acid (two violations). These analyses show that these compounds may have a low bioavailability and pharmacokinetic properties, which make it difficult for them to be administered orally. All compounds were compared with the commercial 5-LOX inhibitor (zileuton), showing that usnic acid, squaric acid, psoromic acid, orsellinic acid, atranorin, cetraric acid, squamatic acid and methylorsellinate presented behaviors similar to zileuton, which did not present any violation of the Lipinski rules. Additionally, we observed that the compounds thamnolic acid and FP acid only presented a violation of the Lipinski rules, which are in the limit to be accepted as possible inhibitors of 5-LOX and to be administered orally; however, in silico and in vitro analyses are required to enrich and further support the in silico analysis. The bioavailability of compounds presents in Cladonia is also assessed by topological polar surface area analysis (TPSA) [24][25][26]. This parameter is highly related to the passive molecular transport that a drug can undergo through cell membranes and can predict the transport properties that the compounds can present, as well as their bioavailability. With the TPSA parameter, it was possible to predict the absorption percentage of each of the compounds, comparing them with the commercial inhibitor Zileuton [24]. The compounds that presented a higher percentage of absorption compared to zileuton (76.30%) were squaric acid (83.26%), orselinic acid (82.17%) and methylorselinate (85.97%).

Toxicity Prediction
To calculate the toxicity risk parameters that the compounds in the species of the genus Cladonia could present, the Osiris Data Warrior computational tool was used [25,26]. This computational tool allows for the evaluation of the molecules according to their chemical structure and the possible molecular fragments that can cause some toxicological risk. The toxicity risk parameters that were calculated are mutagenicity, tumorigenicity, irritation and reproductive toxicity [25,26]. The results are displayed using a blue color scale, where the most intense blue color indicates that the compound exhibits a high risk of toxicity, an intermediate blue color indicates that the compound exhibits an intermediate toxicity risk and a weaker blue color indicates that the compound does not present any toxicity risk (Figure 3). The compounds that presented risks of toxicity were usnic acid, squaric acid and FP acid. Usnic acid presented a high risk of toxicity in the reproductive effect because it presented a fragment in its chemical structure that confers this toxicity effect. The squaric acid exhibited a high probability of presenting tumorigenic properties because the entire chemical structure of squaric acid gives it these risks and toxicological properties. The FP acid compound presented various toxicological risks because its chemical fragmentation forms unstable chemical structures that give it these toxicological properties. FP acid exhibited five fragments with a high risk of an irritant effect; likewise, low risks of mutagenicity, tumorigenicity and toxicity in the reproductive effect were observed (Table 5).  further support the in silico analysis. The bioavailability of compounds presents in Cladonia is also assessed by topological polar surface area analysis (TPSA) [24][25][26]. This parameter is highly related to the passive molecular transport that a drug can undergo through cell membranes and can predict the transport properties that the compounds can present, as well as their bioavailability. With the TPSA parameter, it was possible to predict the absorption percentage of each of the compounds, comparing them with the commercial inhibitor Zileuton [24]. The compounds that presented a higher percentage of absorption compared to zileuton (76.30%) were squaric acid (83.26%), orselinic acid (82.17%) and methylorselinate (85.97%).

Toxicity Prediction
To calculate the toxicity risk parameters that the compounds in the species of the genus Cladonia could present, the Osiris Data Warrior computational tool was used [25,26]. This computational tool allows for the evaluation of the molecules according to their chemical structure and the possible molecular fragments that can cause some toxicological risk. The toxicity risk parameters that were calculated are mutagenicity, tumorigenicity, irritation and reproductive toxicity [25,26]. The results are displayed using a blue color scale, where the most intense blue color indicates that the compound exhibits a high risk of toxicity, an intermediate blue color indicates that the compound exhibits an intermediate toxicity risk and a weaker blue color indicates that the compound does not present any toxicity risk (Figure 3). The compounds that presented risks of toxicity were usnic acid, squaric acid and FP acid. Usnic acid presented a high risk of toxicity in the reproductive effect because it presented a fragment in its chemical structure that confers this toxicity effect. The squaric acid exhibited a high probability of presenting tumorigenic properties because the entire chemical structure of squaric acid gives it these risks and toxicological properties. The FP acid compound presented various toxicological risks because its chemical fragmentation forms unstable chemical structures that give it these toxicological properties. FP acid exhibited five fragments with a high risk of an irritant effect; likewise, low risks of mutagenicity, tumorigenicity and toxicity in the reproductive effect were observed (Table 5).

Evaluation of the Docking 5-LOX Inhibition
The directed coupling of 9 of the 12 compounds present in the species of the genus Cladonia was performed with the enzyme 5-LOX (PDB ID: 6N2W) to see the possible interactions with the catalytic site. The 5-LOX catalytic site is composed of a tetrad of catalytic residues (His367, His372, His550 and Leu673) that, together with the iron atom, coordinate for its catalysis. Ten conformations were tested in each of the compounds to be evaluated, selecting the best conformation that presented a better binding affinity (∆G, kcal/mol) and comparing them with a commercial 5-LOX inhibitor (Zileuton) [25,26].
The results in Table 6 showed that the one with the best binding affinity for the 5-LOX enzyme was psoromic acid (−8.30 kcal/mol); furthermore, the binding energy shown by this compound was very close compared to the commercial inhibitor zileuton (−8.70 kcal/mol). The compounds thamnolic acid, perlatolic acid, atranorin, cetraric acid and squamatic acid exhibited similar behaviors in terms of binding affinity (−7.40, −7.10, −7.20, −7.60 and −7.60 kcal/mol, respectively); however, they were inferior to zileuton and psoromic acid. The compounds that presented a lower binding affinity were 13D acid (−5.60 kcal/mol), orsellinic acid (−5.60 kcal/mol) and methylorsellinate (−5.40 kcal/mol); however, these compounds showed specific interactions with the catalytic site of 5-LOX. Note: a Percentage of absorption (%ABS); b topological polar surface area (TPSA); c molecular weight (MW); d logarithm of partition coefficient between n-octanol and water (cLogP); e number of hydrogen bond donors (HBD); f number of hydrogen bond acceptors (HBA); g number of rotable bonds (n-ROTB). * Commercial redox inhibitors of 5-LOX.

Evaluation of the Docking 5-LOX Inhibition
The directed coupling of 9 of the 12 compounds present in the species of the genus Cladonia was performed with the enzyme 5-LOX (PDB ID: 6N2W) to see the possible interactions with the catalytic site. The 5-LOX catalytic site is composed of a tetrad of catalytic residues (His367, His372, His550 and Leu673) that, together with the iron atom, Note: a Percentage of absorption (%ABS); b topological polar surface area (TPSA); c molecular weight (MW); d logarithm of partition coefficient between n-octanol and water (cLogP); e number of hydrogen bond donors (HBD); f number of hydrogen bond acceptors (HBA); g number of rotable bonds (n-ROTB). * Commercial redox inhibitors of 5-LOX.

Evaluation of the Docking 5-LOX Inhibition
The directed coupling of 9 of the 12 compounds present in the species of the genus Cladonia was performed with the enzyme 5-LOX (PDB ID: 6N2W) to see the possible interactions with the catalytic site. The 5-LOX catalytic site is composed of a tetrad of catalytic residues (His367, His372, His550 and Leu673) that, together with the iron atom, Note: a Percentage of absorption (%ABS); b topological polar surface area (TPSA); c molecular weight (MW); d logarithm of partition coefficient between n-octanol and water (cLogP); e number of hydrogen bond donors (HBD); f number of hydrogen bond acceptors (HBA); g number of rotable bonds (n-ROTB). * Commercial redox inhibitors of 5-LOX.

Evaluation of the Docking 5-LOX Inhibition
The directed coupling of 9 of the 12 compounds present in the species of the genus Cladonia was performed with the enzyme 5-LOX (PDB ID: 6N2W) to see the possible interactions with the catalytic site. The 5-LOX catalytic site is composed of a tetrad of catalytic residues (His367, His372, His550 and Leu673) that, together with the iron atom,

High-risk fragment indicating an irritant effect
Note: FP acid (furmarprotocetraric acid). Table 6. Binding energies resulting from molecular docking experiments of the selected compounds in the extracts of Cladonia, together with the standard inhibitor zileuton on 5-lipoxygenase (5-LOX).   (Figures 4A and 5A). The carboxylate present in one of the aromatic rings of the thamnolic acid presented strong interactions of the charge attraction type with the following residues: an Arg596 distance of 1.82 Å and a His372 distance of 3.92 Å (Figures 4A and 5A); the latter is highly involved in the catalytic site for the inhibition of the 5-LOX enzyme.

Compound Binding Energy (Kcal/mol) 5-Lipoxygenase (5-LOX)
Carbon hydrogen bond type interactions were observed, with a His367 distance of 2.83 Å, a Leu368 distance of 2.91 Å and a His600 distance of 2.87 Å, where the His367 residue is directly involved in the inhibition of the 5-LOX enzyme. The psoromic acid had hydrogen bonding interactions with the residues Gln363 (distance 1.99 Å), His367 (distance 2.88 Å) and Arg596 (distance 2.46 Å), which maintained the stability of the psoromic acid at the 5-LOX catalytic site (Figures 4D and 5D). Carbon hydrogen interactions occurred between the Pro569 residue (distance 3.42 Å) and the methoxyl group of one of the aromatic rings of psoromic acid; the other carbon hydrogen interaction was shown with the residue Thr364 (distance 2.83 Å) and the carboxylate group of the aromatic ring of psoromic acid (Figures 4D and 5D). As for Figures 6 and 7, Figures 6E and 7E show the molecular interactions between the catalytic site of 5-LOX and orsellinic acid; a hydrogen bridge type bond and a charge attraction interaction were shown with the Arg596 residue (distances of 1.87 Å and 3.13 Å, respectively). Two π-π T-shaped interactions were presented between the pi electrons of the aromatic ring of the orsellinic acid and the π electrons of the residues Phe359 (distance 5.23 Å) and Trp599 (distance 5.22 Å). The Ala603 residue had a π-alkyl type interaction between the π electrons of the aromatic ring of the orsellinic acid and the methyl group of the Ala603 residue (distance 4.70 Å).
Atranorin (Figures 6F and 7F) and cetraric acid ( Figures 6G and 7G) were the compounds that presented unfavorable donor-donor interactions and were precisely present in the hydroxyl group present in the aromatic ring in the structures of the compounds and between the residue Arg596 (distances of 1.17 Å and 1.18 Å); however, it was observed that cetraric acid presented strong hydrogen bonding interactions with the residues His432 (distance 2.01 Å), His367 (distance 2.51 Å) and Pro569 (distance 2.97 Å).
In Figures 6H and 7H, two hydrogen bond interactions are observed, which allowed for the stability of the squamatic acid in the catalytic site; these interactions occurred with His367 (distance 2.66 Å) and Gln363 (distance 2.12 Å) residues. Charge attraction interactions were observed between the carboxylate functional group of the squamatic acid and the amine group of the Arg596 residue (distance 3.31 Å). Figures 6I and 7I show the behavior of methylorsellinate at the catalytic site of the 5-LOX enzyme. In these images, we can observe the interactions that this compound presented with the catalytic site. The methylorsellinate presented hydrogen bonding interactions between the hydroxyl group of the aromatic ring and the residues Arg596 (distance 1.94 Å) and His600 (distance 2.08 Å). Additionally, π-π T-shaped interactions were presented between the π electrons of the aromatic ring of methylorsellinate and the π electrons of the aromatic rings of the residues Phe359 (distance 5.11 Å) and Trp599 (distance 5.38 Å).  The 13D acid presented two strong hydrogen bond type interactions ( Figures 4B and 5B) with the residues Arg596 (distance 1.87 Å) and His367 (distance 2.78 Å), which allowed it to stabilize within the catalytic pocket. Pi-sigma interactions ( Figures 4B and 5B) with His372 (distance 3.47 Å) were observed between the imidazole chain of His372 and the carbon of the tetrahydropyran group of 13D acid. Alkyl interactions (Figures 4B and 5B) were observed between the methylenes of the alkyl chain of 13D acid and the residues Ala410 (distance 4.37 Å), Leu414 (distance 4.46 Å), Trp599 (distance 5.00 Å) and Ala603 (distance 4.89 Å). Figures 4C and 5C show the molecular interactions between the catalytic site of 5-LOX and perlatolic acid. In this figure, three hydrogen bonding interactions were presented with the residues Gln363 (distance 1.93 Å), His367 (distance 2.23 Å) and His432 (distance 2.28 Å). Charge attraction interactions were exhibited between the residue Arg596 (distance 3.13 Å), and it was also observed that this residue also presented hydrogen bond type interactions (distance 2.11 and 2.48 Å) with the carboxylate group of perlatolic acid (Figures 4C and 5C).
The psoromic acid had hydrogen bonding interactions with the residues Gln363 (distance 1.99 Å), His367 (distance 2.88 Å) and Arg596 (distance 2.46 Å), which maintained the stability of the psoromic acid at the 5-LOX catalytic site ( Figures 4D and 5D). Carbon hydrogen interactions occurred between the Pro569 residue (distance 3.42 Å) and the methoxyl group of one of the aromatic rings of psoromic acid; the other carbon hydrogen interaction was shown with the residue Thr364 (distance 2.83 Å) and the carboxylate group of the aromatic ring of psoromic acid ( Figures 4D and 5D).
As for Figures 6 and 7, Figures 6E and 7E show the molecular interactions between the catalytic site of 5-LOX and orsellinic acid; a hydrogen bridge type bond and a charge attraction interaction were shown with the Arg596 residue (distances of 1.87 Å and 3.13 Å, respectively). Two π-π T-shaped interactions were presented between the pi electrons of the aromatic ring of the orsellinic acid and the π electrons of the residues Phe359 (distance 5.23 Å) and Trp599 (distance 5.22 Å). The Ala603 residue had a π-alkyl type interaction between the π electrons of the aromatic ring of the orsellinic acid and the methyl group of the Ala603 residue (distance 4.70 Å).
Atranorin (Figures 6F and 7F) and cetraric acid ( Figures 6G and 7G) were the compounds that presented unfavorable donor-donor interactions and were precisely present in the hydroxyl group present in the aromatic ring in the structures of the compounds and between the residue Arg596 (distances of 1.17 Å and 1.18 Å); however, it was observed that cetraric acid presented strong hydrogen bonding interactions with the residues His432 (distance 2.01 Å), His367 (distance 2.51 Å) and Pro569 (distance 2.97 Å).
In Figures 6H and 7H, two hydrogen bond interactions are observed, which allowed for the stability of the squamatic acid in the catalytic site; these interactions occurred with His367 (distance 2.66 Å) and Gln363 (distance 2.12 Å) residues. Charge attraction interactions were observed between the carboxylate functional group of the squamatic acid and the amine group of the Arg596 residue (distance 3.31 Å). Figures 6I and 7I show the behavior of methylorsellinate at the catalytic site of the 5-LOX enzyme. In these images, we can observe the interactions that this compound presented with the catalytic site. The methylorsellinate presented hydrogen bonding interactions between the hydroxyl group of the aromatic ring and the residues Arg596 (distance 1.94 Å) and His600 (distance 2.08 Å). Additionally, π-π T-shaped interactions were presented between the π electrons of the aromatic ring of methylorsellinate and the π electrons of the aromatic rings of the residues Phe359 (distance 5.11 Å) and Trp599 (distance 5.38 Å). First, the effect of the ethanolic extracts of C. chlorophaea and C. gracilis on the viability of the SH-SY5Y cell line was determined by quantification with the MTT (3-[4,5-dimethylthiazol-2-yl]-2,5 diphenyl tetrazolium bromide) assay. According to previ-ous cytotoxicity assays, five independent concentrations were determined for each extract: C. gracilis with 0.5, 1.0, 5.0, 10 and 25 µg/mL and C. chlorophaea with 0.5, 0.7, 1.0, 3.0 and 5.0 µg/mL. The C. gracilis extract promoted over 90% cell viability at all the concentrations evaluated; on the other hand, the C. chlorophaea extract showed optimal cell viability, but at 1.0 and 5.0 µg/mL, a slight reduction in viability was found (85% and 84%, respectively) ( Figure 8).   Subsequently, the most optimal concentrations of cell viability in the two extracts were selected to determine the protective effector against H 2 O 2 -induced oxidative stress: C. gracilis at 5.0, 10 and 25 µg/mL and C. chlorophaea at 0.5 and 0.7 µg/mL. According to Figure 9, H 2 O 2 treatment decreased the viability of SH-SY5Y cells by 75% compared to control cells. On the other hand, the two extracts promoted cell survival compared to H 2 O 2 treatment; thus, for C. gracilis, the increase in cell viability was 98.7%, 104.7% and 90.2% at 5.0, 10 and 25 µg/mL, respectively, and for C. chlorophaea, the increase was 100% and 88% at 0.5 and 0.7 µg/mL, respectively.

Effect of Cladonia Extracts on the Reduction in ROS Production in an H 2 O 2 -Induced Oxidative Stress Model
In this experiment, the concentrations of the extracts used in the previous assay were used to deepen their cytoprotective potential. According to Figure 10, a significant increase in intracellular ROS production of 143.9% was evidenced in the treatment of SH-SY5Y cells with H 2 O 2 versus the control treatment (100%). In parallel, treatments with ethanolic extracts of Cladonia species significantly reduced ROS production; thus, C. gracilis showed a reduction of 46.5%, 48.3%, and 38% at 5.0, 10 and 25 µg/mL, respectively, versus treatment with H 2 O 2 , and C. chlorophaea showed a reduction of 26.8% and 43.8% at 0.5 and 0.7 µg/mL, respectively, versus treatment with H 2 O 2 .

Effect of Cladonia Extracts on Cell Viability and Survival under Oxidative Str Induced by H2O2
First, the effect of the ethanolic extracts of C. chlorophaea and C. gracilis on the v of the SH-SY5Y cell line was determined by quantification with the MTT (3-[4,5-dim thiazol-2-yl]-2,5 diphenyl tetrazolium bromide) assay. According to previous cytot assays, five independent concentrations were determined for each extract: C. gracil 0.5, 1.0, 5.0, 10 and 25 µg/mL and C. chlorophaea with 0.5, 0.7, 1.0, 3.0 and 5.0 µg/m C. gracilis extract promoted over 90% cell viability at all the concentrations evalua the other hand, the C. chlorophaea extract showed optimal cell viability, but at 1.0 a µg/mL, a slight reduction in viability was found (85% and 84%, respectively) (Figu Subsequently, the most optimal concentrations of cell viability in the two e were selected to determine the protective effector against H2O2-induced oxidative C. gracilis at 5.0, 10 and 25 µg/mL and C. chlorophaea at 0.5 and 0.7 µg/mL. Accord Figure 9, H2O2 treatment decreased the viability of SH-SY5Y cells by 75% compa control cells. On the other hand, the two extracts promoted cell survival compared t treatment; thus, for C. gracilis, the increase in cell viability was 98.7%, 104.7% and

Effect of Cladonia Extracts on the Reduction in ROS Production in an H2O2-Induced Oxidative Stress Model
In this experiment, the concentrations of the extracts used in the previous assay were used to deepen their cytoprotective potential. According to Figure 10, a significant increase in intracellular ROS production of 143.9% was evidenced in the treatment of SH-SY5Y cells with H2O2 versus the control treatment (100%). In parallel, treatments with ethanolic extracts of Cladonia species significantly reduced ROS production; thus, C. gracilis showed a reduction of 46.5%, 48.3%, and 38% at 5.0, 10 and 25 µg/mL, respectively, versus treatment with H2O2, and C. chlorophaea showed a reduction of 26.8% and 43.8% at 0.5 and 0.7 µg/mL, respectively, versus treatment with H2O2. In other lichen species such as Dactylina arctica, Nephromopsis stracheyi, Tuckermannopsis americana and Vulpicida pinastri from cetrarioid clade, ideal viability ranges have been reported with methanolic extracts at concentrations of 5 to 10 µg/mL, along with a reduction of the same between 25 and 50 µg/mL; likewise, cell protection against H2O2 was evidenced, especially in D. artica and V. pinastri, with a viability percentage of 76.8% at 10 µg/mL and 78.9% at 5 µg/mL, respectively. In addition, for D. arctica, a reduction in ROS production of 42.4% was reported at 10 µg/mL compared to H2O2 treatment [12], which is similar to that obtained with C. gracilis species at concentrations of 5 and 10 µg/mL. This effective cytoprotective activity is triggered by major compounds in the extracts such as usnic acid and fumarprotetratric acid, which report high antioxidant properties that promote the reduction of oxidative damage in different cell lines [52]. In other studies, isolated compounds such as evernic acid, which are abundantly present in lichens of the genera Ramalina, Evernia and Hypogymnia, report a cell viability >50% using primary culture neurons [11]; biruloquinone, isolated from C. macilenta [47], showed a protective effect against H2O2-induced oxidative stress >80% at a concentration of 25 µg/mL and a reduction in ROS production at 25 and 12.5 µg/mL, which are similar to the values recorded for C. gracilis at working concentrations. In parallel, the viability assays can be contrasted with cytotoxicity analyses on cancer cell lines that show the promising use of com- In other lichen species such as Dactylina arctica, Nephromopsis stracheyi, Tuckermannopsis americana and Vulpicida pinastri from cetrarioid clade, ideal viability ranges have been reported with methanolic extracts at concentrations of 5 to 10 µg/mL, along with a reduction of the same between 25 and 50 µg/mL; likewise, cell protection against H 2 O 2 was e-videnced, especially in D. artica and V. pinastri, with a viability percentage of 76.8% at 10 µg/mL and 78.9% at 5 µg/mL, respectively. In addition, for D. arctica, a reduction in ROS production of 42.4% was reported at 10 µg/mL compared to H 2 O 2 treatment [12], which is similar to that obtained with C. gracilis species at concentrations of 5 and 10 µg/mL. This effective cytoprotective activity is triggered by major compounds in the extracts such as usnic acid and fumarprotetratric acid, which report high antioxidant properties that promote the reduction of oxidative damage in different cell lines [52]. In other studies, isolated compounds such as evernic acid, which are abundantly present in lichens of the genera Ramalina, Evernia and Hypogymnia, report a cell viability >50% using primary culture neurons [11]; biruloquinone, isolated from C. macilenta [47], showed a protective e-ffect against H 2 O 2 -induced oxidative stress >80% at a concentration of 25 µg/mL and a reduction in ROS production at 25 and 12.5 µg/mL, which are similar to the values recorded for C. gracilis at working concentrations. In parallel, the viability assays can be contrasted with cytotoxicity analyses on cancer cell lines that show the promising use of compounds isolated from lichens for pharmacological treatment, such as the case of (-)-usnic, fumarprotocetraric and 9 (O-methyl)protocetraric acids from C. convoluta [53] and di-ffractaic acid, vicanicin, lobaric acid, variolaric acid, protolichesterinic acid and usnic acid from Protousnea magellanica, Psoroma pallidum, Stereocaulon alpinum, Ochrolechia deceptionis, Cornicularia aculeata and C. lepidophora, respectively [54].

Conclusions
This work contributed to the understanding of the metabolomics of the species of the genus Cladonia, C. chlorophaea and C. gracilis, which are present in the Antarctic territory, reporting the identification of 19 compounds for each one by UHCPL-ESI-QTOF-MS analysis, including aromatics, organic acids, lipids, dibenzofurans, depsides and depsidones, which are typical of lichenized fungi. C. chlorophaea and C. gracilis species presented variable but significant results in terms of the phenolic compound content, FRAP ORAC antioxidant assays, cholinesterases (AChE and BChE) and α-glucosidase inhibition and moderate values in terms of the DPPH antioxidant assay and pancreatic lipase inhibition. On the other hand, in an anti-inflammatory potential approach, molecular docking evidenced significant binding affinities of some compounds for the 5-LOX enzyme, together with outstanding pharmacokinetic properties that compromise their good bioavailability. Likewise, the two extracts of the Cladonia species promoted protective activity in the human neuroblastoma cell line against cytotoxicity effects generated by hydrogen peroxide. These species represent a biological resource of scientific interest for the prevention and treatment of central nervous system pathologies, inflammatory disorders and metabolic alterations. Future research could address the isolation and study of bioactive compounds of both species and their molecular mechanisms of action.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors.