Traditionally Used Lathyrus Species: Phytochemical Composition, Antioxidant Activity, Enzyme Inhibitory Properties, Cytotoxic Effects, and in silico Studies of L. czeczottianus and L. nissolia

Members of the genus Lathyrus are used as food and as traditional medicines. In order to find new sources of biologically-active compounds, chemical and biological profiles of two Lathyrus species (L. czeczottianus and L. nissolia) were investigated. Chemical profiles were evaluated by HPLC-ESI-MSn, as well as by their total phenolic and flavonoid contents. In addition, antioxidant, enzyme inhibitory, and cytotoxic effects were also investigated. Antioxidant properties were tested by using different assays (DPPH, ABTS, CUPRAC, FRAP, phosphomolybdenum, and metal chelation). Cholinesterases (AChE and BChE), tyrosinase, α-amylase, and α-glucosidase were used to evaluate enzyme inhibitory effects. Moreover, vitexin (apigenin-8-C-glucoside) and 5-O-caffeoylquinic acid were further subjected to molecular docking experiments to provide insights about their interactions at molecular level with the tested enzymes. In vitro cytotoxic effects were examined against human embryonic kidney cells (HEK293) by using iCELLigence real time cell analysis system. Generally, L. czeczottianus exhibited stronger antioxidant properties than L. nissolia. However, L. nissolia had remarkable enzyme inhibitory effects against cholinesterase, amylase and glucosidase. HPLC-ESI-MSn analysis revealed that flavonoids were major components in these extracts. On the basis of these results, Lathyrus extracts were rich in biologically active components; thus, these species could be utilized to design new phytopharmaceutical and nutraceutical formulations.


INTRODUCTION
With recent advances in science, the therapeutic properties of plants have gained importance all over the world, due to their pharmacological activities (antioxidant, anti-cancer, antiinflammatory, antimicrobials, etc.) and nutraceuticals properties (Krishnaiah et al., 2011). Many reports have focused on benefits of phytochemicals from plants and their health-promoting effects. Among the various phytochemicals with biological properties, phenolic compounds are responsible for the abovementioned biological activities (Shahidi and Ambigaipalan, 2015). Actually, in global communities increasingly concerned on health and nutrition, these products are emerging as valuable alternatives for replacing synthetic drugs. On the basis of the aforementioned concerns, new plant-derived products or phytomedicines are considered as prospective materials.
The genus Lathyrus, which belongs to Fabaceae, is represented by more than 200 species worldwide (Alkin et al., 1986). This genus consists of 75 taxa, of which 18 species are endemic in Turkish Flora (Genç, 2009). Like most members of Fabaceae, Lathyrus genus presents different uses, such as alimentary, agricultural, industrial, ornamental, and in traditional medicine (Patto and Rubiales, 2014). For instance, some Lathyrus species such as L. hirsutus L., L. cicera L., L. odoratus L., L. ochrus L. (DC.), L. sylvestris L., and L. palustris L., which have high nutrient quality, are consumed as food (for both humans and animals) (Chavan et al., 2001). Especially, roots of L. odoratus are used as food in The East Anatolian Region of Turkey. Again, L. odoratus and L. sativus L. are used for agricultural processes in Turkey. Moreover, Lathyrus species have also importance in folk medicine for several purposes such as analgesic (seed of L. sativus), anti-inflammatory (aerial parts of L. cicera), and anti-rheumatism (leaf of L. rotundifolius Willd. subsp. miniatus (Bieb. ex Stev.) Davis) in Turkey (Altundag and Ozturk, 2011). However, up-to-date information about L. czeczottianus Bässler and L. nissolia L. chemical composition or biological characterization is scarce. Considering the importance of different Lathyrus species, as previously mentioned, the general aim of this research is to improve the knowledge about lessknown Lathyrus plants. In this direction, this study was designed to investigate: (i) the antioxidant activity of the methanolic extracts of L. czeczottianus and L. nissolia using various in vitro biochemical assays, (ii) the enzyme inhibitory potentials (anti-cholinesterase, anti-amylase, anti-glucosidase, and antityrosinase), (iii) the cytotoxicity in human embryonic kidney cells using the real time cellular impedance technology, (iv) the characterization of phytochemical compounds in the methanolic extracts by HPLC-ESI-MS, (v) and the components-enzymes interactions with in-silico techniques. Thus, our present study may contribute to offer new perspectives on the biological properties and phytochemical profile of the genus Lathyrus.

Plant Material and Extraction Procedure
Lathyrus species were collected at flowering stage and corresponding information and localities are explained below. Taxonomic identification of the plant materials was confirmed by senior taxonomist Dr. Murad Aydın Sanda [Selcuk University, Science Faculty, Department of Biology (Botany)] based on Flora of Turkey (Davis, 1970) and A checklist of the Flora of Turkey (Vascular Plants) (Güner et al., 2012). Voucher specimens were deposited at the KNYA Herbarium of Department of Biology, Selcuk University, Konya-Turkey. The plant materials were dried at room temperature. The dried aerial parts were ground to a fine powder using a laboratory mill.
To obtain methanolic extracts, the air-dried aerial parts (10 g) were macerated with 200 mL of methanol at room temperature (25 • C ± 1 • C) for 24 h. The extracts were concentrated under vacuum at 40 • C by using a rotary evaporator and stored at + 4 • C in dark until use.

Chromatographic Conditions
The HPLC system consisted of a vacuum degasser, an autosampler and a binary pump (Agilent Series 1100, Agilent Technologies, Santa Clara, CA, USA) equipped with a reversed phase Kinetex core-shell C 18 analytical column of 50 × 2.1 mm and 2.6 µm particle size (Phenomenex, Torrance, CA, USA). A C 18 Security Guard Ultra cartridge (Phenomenex) of 2.1 mm i.d. was placed before the analytical column. The best separation was achieved by using a mobile phase consisting of acetonitrile (A) and water-formic acid (100:0.1, v/v) (B). The following gradient program was used: 10% A (0 min), 25% A (10-20 min), 50% A (40 min), 100% A (42-47 min), and 10% A (49 min). The mobile phase flow rate was 0.4 mL min −1 . After filtration through 0.45 µm PTFE membrane filters, 10 µL of each extract was injected.
The HPLC system was connected to an ion trap mass spectrometer (Esquire 6,000, Bruker Daltonics, Billerica, MA, USA) equipped with an electrospray interface operating in negative ion mode. The scan range was set at m/z 100-1200 with a speed of 13,000 Da/s. The ESI conditions were as follows: drying gas (N 2 ) flow rate and temperature, 10 mL/min and 365 • C; nebulizer gas (N 2 ) pressure, 50 psi; capillary voltage, 4,500 V; capillary exit voltage, −117.3 V. The acquisition of MS n data was made in auto MS n mode, with isolation width of 4.0 m/z, and fragmentation amplitude of 0.6 V (MS n up to MS 4 ). Esquire control software was used for the data acquisition and Data Analysis for processing.

Receptors Preparation
All the crystallographic enzyme structures have been downloaded from the Protein Databank RCSB PDB (Berman et al., 2000): acetylcholinesterase (pdb:4X3C) (Pesaresi and Lamba) in complex with tacrine-nicotinamide hybrid inhibitor, butyrilcholinesterase (pdb:4BDS) (Nachon et al., 2013) in complex with tacrine, amylase (pdb:1VAH) (Zhuo et al., 2004) in complex with r-nitrophenyl-α-D-maltoside, glucosidase (pdb:3AXI) (Yamamoto et al., 2011) in complex with maltose, and tyrosinase (pdb:2Y9X) (Ismaya et al., 2011) in complex with tropolone. Waters, inhibitors and all the other molecules present in the pdb files were removed by using Pymol (DeLano, 2002) and the proteins alone were neutralized at pH 7.4 by Epik suit implemented in Maestro 9.2 suite (Maestro, 2011). Selenocysteines and seleno-methionines, if present, were converted respectively to cysteines and methionines. All the missing fragments and other errors present in the crystal structures were automatically solved by the Wizard Protein Preparation implemented in Maestro 9.2 suite (Maestro, 2011).

Ligands Preparation
5-O-caffeoylquinic acid (5-O-CQA, compound 4), isoschaftoside (compound 22) and vitexin (compound 33) were selected as representative compounds to carry out molecular docking studies, as these compounds were present in L. czeczottianus and L. nissolia extracts. The chemical structures, reported in Scheme 1, have been downloaded from Zinc databases (Irwin et al., 2012) and used for molecular modeling experiments. The ligands were prepared by the LigPrep tool embedded in Maestro 9.2, neutralized at pH 7.4 by Epik and minimized (Shelley et al., 2007).

Molecular Docking
Dockings of vitexin (33), isoschaftoside (22) and 5-O-CQA (4) have been performed for each selected enzyme employed for the in vitro enzymatic inhibition tests in this work. Glide (Friesner et al., 2006) has been employed for the docking calculations by using all the eXtra Precision scoring function for all the enzymes, with the exception of the docking to tyrosinase; in this case, Gold 6.0 (Jones et al., 1997) was used with the scoring function Gold Score, which has been previously found to be more suitable to produce reliable poses on metal containing enzymes (Mocan et al., 2016b). In both cases, the binding pocket was determined automatically by centering the grid on the crystallographic inhibitor, extended in a radius of 10 Angstroms from the center. The best pose for each compound docked to the selected enzymes was the best ranked among the 200 generated.

Cytotoxicity Assay with the iCELLigence System
Cytotoxic activity was performed by using the iCELLigence real time cell analysis technology (RTCA) (ACEA Biosciences, San Diego, CA, USA) as described previously (Lazarova et al., 2015). In brief, HEK-293 cells were seeded in E-Plate L8 at a density of 5 × 10 4 cells/well in 300 µL of growth medium. The system was set to take measurements every 30 min for 72 h. After 24 h, when the cells were in exponential phase, they were treated with 300 µL medium containing 2,000, 1,000, 500, 250, 125, 62.5 µg/mL of L. czeczottianus or L. nissolia methanolic extracts prepared in DMEM containing 0.1% dimethyl sulfoxide (DMSO). Cells were treated with 0.1% of DMSO as negative control and 5% DMSO as positive control in all experiments. Data analysis was performed using RTCA iCELLigence software and IC 50 values (µL/mL) were automatically calculated based on the final readings obtained at the 72 h time point. All experiments were repeated three times.

Morphological Evaluation of the Cells
Cells were plated at a density of 5 × 10 5 cells/well in 6well plate with 2 mL growth medium. Twenty-four hours after seeding, the medium was removed and replaced with a medium containing 500 µg/mL L. czeczottianus and 1,000 µg/mL L. nissolia methanolic extracts, which are close to IC 50 values. After 72 h, the cells were washed with 1X phosphate buffered saline (PBS) to remove cellular debris and photographed under an inverted microscope Leica DM IL LED using a DFC-290 camera (Leica, Wetzlar, Germany).

Statistical Analysis
For all the experiments, the assays were carried out in triplicate. The results are expressed as mean values and standard deviation (SD). The differences between the different extracts were analyzed using one-way analysis of variance (ANOVA) followed by Tukey's honestly significant difference post-hoc test with α = 0.05. This treatment was carried out using SPSS v. 14.0 program.

HPLC-ESI-MS n
The analysis of the phenolic composition in aerial parts of L. czeczottianus and L. nissolia was performed by HPLC-ESI-MS n using the negative ionization mode. Two independent assays were carried out for each sample, obtaining similar data regarding the nature and relative intensities of the detected fragments. The base peak chromatograms of the methanolic extracts are shown in Figure 1.
The initial step for the characterization of the phenolic compounds consisted in the determination of the molecular weight of each compound. In the negative ionization mode (ESI − ) MS 1 spectrum, the most intense peak usually corresponded to the deprotonated molecular ion [M-H] − . Compounds were numbered in both chromatograms by their order of elution. The tentative characterization of the detected compounds is shown in Table 1, and the discussion of the characterization is explained in the following sub-sections.

Phenolic Acids
Compound 3 exhibited the deprotonated molecular ion at m/z 315, and suffered the neutral loss of 162 Da (hexoside) to yield a fragment ion at m/z 153. It was characterized as dihydroxybenzoic acid hexoside (Engels et al., 2012).
Four caffeoylquinic acids were detected in the analyzed extracts. Compound 4, with [M-H] − at m/z 353, base peak ion at m/z 191, and low intensity at m/z 179, was characterized as 5-caffeoylquinic acid (Clifford et al., 2003). Compounds 42 and 49 exhibited the deprotonated molecular ions at m/z 515 and typical base peaks at m/z 353. Following the hierarchical key for the identification of di-caffeoylquinic acids (Clifford et al., 2003), 42 and 49 were identified as 3,5-dicaffeoylquinic acid and 4,5-dicaffeoylquinic acid, respectively. Finally, compound 55 was characterized as a tricaffeoylquinic acid due to its [M-H] − ion at m/z 677 and MS 2 and MS 3 base peak ions at m/z 515 and 353 (Han et al., 2008).
Compound 5, with [M-H] − at m/z 325, suffered the neutral loss of 162 Da, yielding a fragment ion at m/z 163 (coumaric acid, due to its typical fragment ion at m/z 119). Hence, it was identified as coumaric acid-O-hexoside.
Compound 10 was identified as epicatechin considering its mass spectrum and retention time.
Compound 13 exhibited the deprotonated molecular ion at m/z 339, and the 179→135 transition, typical from caffeic acid, was observed in the MS 3 mass spectrum. Without further information, it was tentatively characterized as a caffeic acid derivative.
Compound 17, with [M-H] − at m/z 367, displayed the base peak ion at m/z 191 and a secondary fragment ion at m/z 173. According to (Clifford et al., 2003), it was identified as 5feruloylquinic acid. Compound 25 was also characterized as a feruloylquinic acid. Compound 51 exhibited the deprotonated molecular ion at m/z 359 and presented fragment ions at m/z 197, 179, 161, and 133, which are typical of rosmarinic acid (Liu et al., 2007).

Flavonoids
Most of the observed compounds in the analyzed extracts were flavonoids. They will be discussed based on their aglycone. The mass spectra of the glycosides showed the aglycone ions as a result of the loss of moieties such as hexosyl (162 Da), pentosyl (132 Da), or rutinosyl (308 Da).
Seven kaempferol derivatives were detected, all of them showing the aglycone kaempferol at m/z 285 (kaempferol standard was analyzed to compare its mass spectrum). Compound 9, with [M-H] − at m/z 755, suffered consecutive losses of 162 and 308 Da, and was tentatively characterized as kaempferol-O-hexoside-O-rutinoside. Compounds 21 and 29 exhibited the deprotonated molecular ion at m/z 609 and suffered two losses of 162 Da, so they were characterized as kaempferol-O-hexoside-O-hexoside. Compounds 28 and 37 were merely characterized as kaempferol derivatives as the only fragment ion observed in MS 2 corresponded to kaempferol, and the identity of the attached moieties could not be elucidated. Compounds 41 and 47 exhibited neutral losses of hexoside and pentoside moieties, respectively, and were characterized as kaempferol-O-hexoside and kaempferol-O-pentoside, respectively.
Compound 54, with [M-H] − at m/z 301 and fragment ions at m/z 179 and 151, was identified as quercetin after comparison with an analytical standard. Six quercetin derivatives were detected in the analyzed extracts. The characterization of compounds 18, 20, 32, 35, and 40 (Table 1) was carried out considering the observed neutral losses of 162 Da (hexoside), 132 Da (pentoside), 146 Da (deoxyhexoside) and 308 Da (rutinoside), to yield the aglycone quercetin at m/z 301. Compound 31 was merely characterized as a quercetin derivative.
Several methylated flavonoids, as aglycones or as glycosides, were tentatively characterized. Compound 56, with [M-H] − at m/z 329, suffered two consecutive losses of methyl groups (15 Da), and was characterized as a dimethyl-flavonoid. Compound 50 suffered the neutral loss of 162 Da to yield the aglycone, which exhibited the same fragmentation pattern than 56, so it was characterized as its O-hexoside. Compounds 30 and 52 suffered neutral losses of 162 Da, yielding aglycones that presented the pattern of methylated flavonoids, so they were characterized as methyl-flavonoid-O-hexoside. In a similar way, 45 was assigned to a methyl-flavonoid-O-rutinoside. Finally, compounds 61 and 64 suffered three losses of methyl groups and were characterized as trimethyl-flavonoids.

Other Compounds
Compound 1, with the deprotonated molecular ion at m/z 191 and MS 2 base peak at m/z 111, was identified as citric acid (Flores et al., 2012).
By comparison of their fragmentation patterns and bibliographic information (Van Hoyweghen et al., 2014), 57 and 58 were characterized as oxo-dihydroxy-octadecenoic acids, and 60 as trihydroxy-octadecenoic acid.
Finally, two saponins were characterized in the analyzed extracts. Compounds 62 and 66, with deprotonated molecular ions at m/z 971 and 941, respectively, were identified after comparison of their fragmentation patterns ( Table 1) with bibliographic data (Pollier et al., 2011).

Antioxidant Properties
Phenolic compounds are known to serve as multifunctional bioactive components (anti-oxidant, anti-microbial, antiinflammatory and anti-cancer agents). Several studies also concluded that these components are main contributors on the antioxidant effects of plant extracts (Shahidi and Ambigaipalan, 2015). In this context, they were selected as target compounds in the present study. Total phenolic and flavonoid contents were detected by colorimetric methods, namely Folin-Ciocβlteu and AlCl 3 . The results are shown in Table 2. The total phenolic content in L. czeczottianus (63.16 mgGAE/g extract) was higher than in L. nissolia (26.47 mgRE/g extract). However, L. nissolia (20.97 mgRE/g extract) had higher concentrations of flavonoids than L. czeczottianus (14.16 mgRE/g extract). Apparently, the total flavonoids accounted for more than 50% of the total phenolics in L. nissolia extract. From these results, flavonoids were thus the major components in L. nissolia extract. This finding was confirmed also by HPLC-ESI-MS n analysis.
Free radicals trigger a wide variety of chronic and degenerative diseases, such as Alzheimer's disease, diabetes mellitus and cancer. At this point, DPPH and ABTS radicals are widely used to assess free radical scavenging ability of plant extracts or synthetic compounds. In these assays, the initial absorbances (blue for ABTS and purple for DPPH) are decreased by antioxidants, and these changes are spectrophotometrically measured (517 nm for DPPH and 734 nm for ABTS). As can be seen in Table 3, L. czeczottianus was found to exhibit potent free radical scavenging activities in both DPPH (IC 50 : 1.42 mg/mL) and ABTS (IC 50 : 1.80 mg/mL) assays. The observed remarkable free radical scavenging activity may be explained by the high level of phenolics in the extracts. This is in accordance with the data published by other authors, indicating that the higher the total phenolic contents, the better the free radical scavenging activity (Bannour et al., 2016;Zhang et al., 2016). Also, some phenolics were reported as effective hydrogen donors by means of their hydroxyl groups (Sevgi et al., 2015). Reducing power is considered as an important antioxidant feature and it is reflected by the electron-donating ability of antioxidants. For this purpose, CUPRAC (Cu 2+ →Cu + ) and FRAP (Fe 3+ →Fe 2+ ) assays were performed. L. czeczottianus exhibited stronger reducing abilities in both assays as similar to the free radical scavenging abilities (Table 3). From these results, the high contents of phenolics in the extracts indicated that these compounds contribute to the higher reducing abilities. Similar results were obtained by several researchers, who reported a strong relationship between total phenolic content and reducing abilities (Beghlal et al., 2016;Smeriglio et al., 2016). Phosphomolybdenum assay is based on the reduction of Mo (VI) to Mo (V) by antioxidants at acidic pH; then, green phosphate-Mo (V) complex is formed, which has its maximum absorbance at 695 nm. The phosphomolybdenum activity was observed as lower for L. czeczottianus (EC 50 : 1.29 mg/mL) than L. nissolia (EC 50 : 0.86 mg/mL). (Table 3). Because phosphomolybdenum assay is considered as a total antioxidant capacity assay, our findings for this assay might be attributed also to non-phenolic antioxidants such as tocopherol or vitamin C. These results were consistent with other reports (Albayrak et al., 2010;Zengin et al., 2015). Also, some researchers observed a weak correlation between total phenolic content and phosphomolybdenum assay (Nićiforović et al., 2010;Sarikurkcu et al., 2015).
Transition metals (iron and copper) play an important role in Fenton and Haber-Weiss reactions, which are main reactions involved in the production of free radicals. In this sense, the  chelating ability of these metals is considered to have a significant role among the antioxidant mechanisms. The metal chelating ability of Lathyrus extracts was determined by ferrozine method and tabulated in Table 3. In contrast to the free radical and reducing power assays, L. nissolia (IC 50 : 1.14 mg/mL) exhibited stronger metal chelating activity compared to L. czeczottianus (IC 50 : >5 mg/mL). Thus, the observed activity for L. nissolia may be caused from non-phenolic chelators, including peptides, polysaccharides or citric acid. In accordance with our results, some authors reported that the metal chelating activity of phenolics plays a minor role among antioxidant mechanisms, and is suggested by several researchers as a secondary antioxidant mechanism (Rice-Evans et al., 1996;Gursoy et al., 2009;Wang et al., 2009). To the best of our knowledge, this study is the first report on the antioxidant properties of these Lathyrus species.
In this direction, the obtained results could be useful to provide new perspectives on the Lathyrus genus, and its species could be regarded as potential candidates for developing novel phytopharmaceuticals.

Enzyme Inhibitory Properties
Alzheimer's disease (AD) and diabetes mellitus (DM) are major health problems due to the severe increase of their prevalence. For example, recent reports indicated DM prevalence to be 9% among adults, and it is supposed to be well above 15% until 2025 (Waltenberger et al., 2016). Today, more than 50 million people are also affected by AD, and this number is estimated to increase to 131.5 million people by 2,050 (Prince, 2015). Therefore, new effective strategies have vital importance in the treatment of these diseases (Rescigno et al., 2002;Mao et al., 2015;Zengin et al., 2016). Moreover, key enzyme inhibition theory is accepted as one of the most efficient strategies to counteract these diseases (Etxeberria et al., 2012;Murray et al., 2013;Mocan et al., 2016a,b). At this point, enzyme inhibitory assays are considered a very useful tool for investigating the biological properties of herbal extracts or isolated compounds Dinparast et al., 2016). Considering the mentioned aspects, the enzyme inhibitory properties of L. czeczottianus and L. nissolia toward cholinesterases, α-amylase, α-glucosidase, and tyrosinase were investigated, and the results are gathered in Table 4.
The inhibitory effects on cholinesterases of both species are modest, a higher effect being observed for L. nissolia (IC 50 : 1.22 mg/mL for AChE). Additionally, L. czeczottianus extract exhibited a high anti-tyrosinase effect, whereas L. nissolia was not active toward this enzyme. High anti-tyrosinase effects and modest inhibitory effects on cholinesterases were already reported in our previous studies concerning L. aureus and L. pratensis (Llorent-Martínez et al., 2016). Moreover, higher inhibitory effects toward α-glucosidase in comparison with αamylase were observed in this study, the results being in line with our previous findings (Llorent-Martínez et al., 2016;Mocan et al., 2016b). Nonetheless, previous reports indicated that Frontiers in Pharmacology | www.frontiersin.org phenolic compounds have lower α-amylase inhibitory activity and a stronger inhibition activity against yeast α-glucosidase (Apostolidis et al., 2011). According to our results, the Lathyrus species have similar α-amylase inhibition and but L. nissolia exerted stronger α-glucosidase inhibition as compared to L. czeczottianus. Also, additional computational studies (molecular docking) were performed to understand possible interactions between the target compounds [5-O-CQA (4), isoschaftoside (22) and vitexin (33)] and tested enzymes as well as to get further insights into the modulation of biological response.

Molecular Modeling
Enzymatic assays performed on the extract of L. czeczottianus evidenced a modest inhibition activity to AchE, and strong tyrosinase, whereas the L. nissolia extracts have shown relevant inhibition activity for AchE and α-glucosidase. The best and most representative enzyme-ligand complexes are reported in Figures 2-8. Among the identified bioactive compounds, vitexin (33) and 5-O-CQA (4), which were found to be dominant in both extracts (especially vitexin), and isoschaftoside (22) which is prevalent in the extract of L. nissolia, were selected for docking experiments. These compounds have been docked to AchE, BChE, α-glucosidase, α-amylase and tyrosinase.
The best pose found for 5-O-CQA (4) docked to AchE (docking score −9.01 Kcal/mol) was stabilized by several interactions with the residues surrounding the enzymatic pocket involving Trp84 with one π-π stack and His440 with two hydrogen bonds. The best ranked pose of vitexin (33) to the same enzyme (score = −8.30 Kcal/mol) interacts with the binding pocket of the AchE by two hydrogen bonds with Arg289, and one hydrogen bond with Ser286, Tyr334, Asp72, and Tyr70. The best pose found for iscoschaftoside (22) docked to AchE was stabilized by two hydrogen bonds to Ser286, two hydrogen bonds to Asp285, one hydrogen bond to Tyr121 and a π-π stack to Trp279 (score = −11.05 Kcal/mol). The best pose found for 5-O-CQA (4) in the enzymatic pocket of BchE was stabilized by the formation of one hydrogen bond with His438, Gln197 and Leu286 with a docking score = −9.39 Kcal/mol. The best pose for isoschaftoside (22) was stabilized by the formation of hydrogen bond with Glu197, Ser198, Leu286, Ser 287 and a π-π stack with Tyr332 (score = −10.08). The best pose found for vitexin (33) in the interaction to the enzymatic pocket of BchE was stabilized by the formation of three hydrogen bonds, respectively, with Gln197, Asp70, Thr120, and a much lower docking score: −11.24 Kcal/mol. These data are in agreement with the literature data in which isoschaftoside (22) has shown a high inhibitory activity to AchE and BChE with IC 50 of 0.23 µM and 3.11 µM respectively (Hung et al., 2015).
5-O-CQA (4) binds to the enzymatic pocket of the αamylase by forming a series of hydrogen bonds with Gln63, Val163, Asp197 residues (docking score = −8.50 Kcal/mol), whereas vitexin (33) interacted with the enzymatic pocket by   establishing hydrogen bonds with Tyr62, His201, Glu233, and Asp300 (docking score = −8.27 Kcal/mol). Isoschaftoside (22) best pose established several hydrogen bonds to Asp356, Asp330, Glu233, Asp197, and Trp59, and also two π-π stacks to Trp58 and Trp59. Among the three tested substances, it presented the most favorable docking score (score = −9.69 Kcal/mol) to this enzyme. Several papers were found in scientific literature reporting the antidiabetic properties of some medicinal plants containing isoschaftoside (Song et al., 2012;Khazneh et al., 2016), so supporting the hypothesis that isoschaftoside may be a good inhibitor of α-amylase. In a previous study, vitexin (33) proved to be a potent α-glucosidase inhibitor with a reported EC50 = 0.4 mg/mL, close to that of Acarbose (Yao et al., 2011) and it was found that its best pose interacted with the enzymatic pocket of α-glucosidase by forming three hydrogen bonds with Arg315, Asp215, Arg213, and a π-π stack with the aromatic ring of Phe178 (docking score = −7.29 Kcal/mol). 5-O-CQA (4) interacted only by forming hydrogen bonds with Asp69, Trp158, Glu343, and Asp307 (Docking score = −8.83 Kcal/mol); isoschaftoside (22) interact by forming hydrogen bonds to Glu411, Asp307, Thr310, Asp242 (score = −7.73 Kcal/mol). For tyrosinase, the docking of vitexin (33), isoschaftoside (22) and 5-O-CQA (4) were significantly different. Vitexin (33) and isoschaftoside (22) were both unable to bind directly to the copper atoms present in the enzymatic pocket of tyrosinase also the interactions to the enzymatic pocket are marginal (see Figure  1S in the supporting material); on the contrary, the best pose found for 5-O-CQA (4) evidenced the ability of the molecule to coordinate one Cu atom present at the binding pocket of the enzyme, as depicted in Figure 8. This finding is in full agreement with our enzymatic assays, and also with the reported inhibitory activity found for some caffeoylquininic acid derivatives toward the tyrosinase (Iwai et al., 2004). In conclusion, the docking experiments performed in this work can explain the literature data regarding the inhibitory activity of vitexin (33), isoshaftoside (22) and 5-O-QCA (4), being capable to interact with the residues present in the binding pocket of the considered enzymes. In particular, vitexin (33), which has been reported to be efficacious to control blood glucose in diabetic rats (Choo et al., 2012), is present in the extracts of both plants and can be addressed as one of the compounds responsible for the α-glucosidase inhibitory activity found for both extracts together with isoschaftoside (22) which has shown a remarkable docking score on alphaamylase, whereas 5-O-CQA (4), as demonstrated in the docking experiments, may efficaciously inhibit the tyrosinase activity by binding to one Cu atom present in the enzymatic pocket.

Cytotoxicity Assay
Toxicity is a major issue to consider in drug screening or drug development from natural product derivatives. Some of the active compounds obtained from these natural products have biomedical uses such as anti-neoplastic and anti-inflammatory agents, and the toxicity of these compounds needs to be clarified by cytotoxicity assays in biological systems. In this study, we have tested the cytotoxicity of two different Lathyrus species (L. czeczottianus and L. nissolia) on human embryonic kidney cells, HEK-293, since these cells are widely used in drug toxicity studies as kidney may be one of the major target organ of several toxins (Mantle et al., 2015;Sang et al., 2016). Previous studies reported the use of HEK-293 cell line for neurobiological research as well, based on the findings that HEK-293 cells have many properties in common with neuronal precursors (Shaw et al., 2002;Thomas and Smart, 2005;Lin et al., 2014). For monitoring cytotoxicity, we used iCELLigence system, which enables to get real-time, label-free and non-invasive measurements from living cells (Bird and Kirstein, 2009;Martinez-Serra et al., 2014). In recent years, this system has been one of the gold standards in cell biology research since it measures the effect of the drugs without any additional treatment on the cells. In addition to this, RTCA systems have own data analysis software which calculates IC 50 values automatically. In this study, IC 50 values of L. czeczottianus and L. nissolia were determined as 413.19 ± 16.3 µg/mL and 1040.83 ± 48.8 µg/mL, respectively, from three independent experiments in a time and dose dependent manner after treating cells with six different concentrations of these two extracts (Figure 9). Furthermore, very high doses were tested, close to IC 50 concentrations of the extracts, on HEK-293 cell morphology and almost 50% cellular viability were observed as seen in middle and right panels compared to control cells on the left (Figure 10). From these results, the studied extracts had no cytotoxic activities against HEK-293 cell lines. However, L. czeczottianus had higher cytotoxic effect with lower value of IC 50 as compared to L. nissolia. In this direction, the observed remarkable activity for L. nissolia may be explained by the higher concentration of flavonoids in the extract. In fact, several authors reported that some flavonoids exhibited no cytotoxic effects on various cell lines (Kumar and Pandey, 2013;Sak, 2014). To the extent of our knowledge, only one previous study has been conducted on the cytotoxic evaluation of the Lathyrus species (L. laxiflorus) (Spanou et al., 2010).

CONCLUSIONS
Lathyrus species are still an underestimated source of bioactive compounds, especially flavonoids, being poorly investigated upto-date. The present study is the first report on the biological properties (antioxidant, enzyme inhibitory and cytotoxic effects) and phytochemical profiles of selected Lathyrus species. On one hand, L. czeczottianus exhibited remarkable antioxidant abilities, presenting the highest concentration of phenolics. On the other hand, L. nissolia had promising inhibitory effects on the studied enzymes, except for tyrosinase. Additionally, the extracts exhibited no cytotoxic effects on HEK-293 cell lines. Molecular docking experiments provided new insights on some major components (vitexin, isoschaftoside, and 5-Ocaffeoylquininic acid) by studying their interactions with the enzyme pool, which gave the molecular basis of the biological response found for these species. According to our results, further studies on Lathyrus species could open new perspectives for developing novel health-promoting agents by pharmaceutical and nutraceutical industries.

AUTHOR CONTRIBUTIONS
EJLM, GZ, MLFC, OB, ArA, RC, AnM, AdM, and SU set up and carried out experiments. GOG and AbA executed data analysis.