Molecular Dereplication and In Vitro and In Silico Pharmacological Evaluation of Coriandrum sativum against Neuroblastoma Cells

The aim of this study was to investigate the cytotoxic activity of the Coriandrum sativum (C. sativum) ethanolic extract (CSEE) in neuroblastoma cells, chemically characterize the compounds present in the CSEE, and predict the molecular interactions and properties of ADME. Thus, after obtaining the CSEE and performing its chemical characterization through dereplication methods using UPLC/DAD-ESI/HRMS/MS, PM6 methods and the SwissADME drug design platform were used in order to predict molecular interactions and ADME properties. The CSEE was tested for 24 h in neuroblastoma cells to the establishment of the IC50 dose. Then, the cell death was evaluated, using annexin-PI, as well as the activity of the effector caspase 3, and the protein and mRNA levels of Bax and Bcl-2 were analyzed by ELISA and RT-PCR, respectively. By UHPLC/DAD/HRMS-MS/MS analysis, the CSEE showed a high content of isocoumarins-dihydrocoriandrin, coriandrin, and coriandrones A and B, as well as nitrogenated compounds (adenine, adenosine, and tryptophan). Flavonoids (apigenin, hyperoside, and rutin), phospholipids (PAF C-16 and LysoPC (16:0)), and acylglicerol were also identified in lower amount as important compounds with antioxidant activity. The in silico approach results showed that the compounds 1 to 6, which are found mostly in the C. sativum extract, obey the “Five Rules” of Lipinski, suggesting a good pharmacokinetic activity of these compounds when administered orally. The IC50 dose of CSEE (20 µg/mL) inhibited cell proliferation and promoted cell death by the accumulation of cleaved caspase-3 and the externalization of phosphatidylserine. Furthermore, CSEE decreased Bcl-2 and increased Bax, both protein and mRNA levels, suggesting an apoptotic mechanism. CSEE presents cytotoxic effects, promoting cell death. In addition to the promising results predicted through the in silico approach for all compounds, the compound 6 showed the best results in relation to stability due to its GAP value.


Molecular Interactions and ADME Properties Predictions
The Molecular Electrostatic Potential Map (MEP) is a classic tool in the analysis of chemical activity, especially in the drug designing. In most cases, MEP is used as a qualitative approach. The MPE is based on the calculated properties of the charge density directly from the function of the molecular wave and measures the interaction of a positively charged point with the nuclei and electrons of a molecule. The interaction between the molecules occurs between regions of opposite electrostatic potentials. The improved accuracy of PM6 is particularly valuable for generating electronic descriptors for qualitative structure activity relationships (QSAR). When comparing the maps of electrostatic potentials of compounds 1, 2, and 3 ( Figure 3) found mostly in coriander, we can see the presence of ribose in compound 2 in terms of volume, as adenosine creates an important increase in the total volume when we compare it with the compounds 1 and 3. The negative charge concentration of compounds 1 and 2 (red region of the map) is very similar when we observe the pyrimidine and imidazole rings; however, the ribose found in compound 2 also has a high electronic density. Regarding compound 3, tryptophan, the negative charge concentration in the indole ring, and the higher electron density is present on the carbonyl of the carboxylic acid group. In addition, compound 3 has large differences in shape and volume. Regarding the electrostatic potential maps obtained for compounds 4, 5, and 6, it is possible to observe a small variation in relation to the shape and volume. The negative charge distribution is also relatively similar, mainly between compounds 4 and 5, where it is predominant over the carbonyl of isochron-1-one.
The Molecular Electrostatic Potential Map (MEP) is a classic tool in the analysis of chemical activity, especially in the drug designing. In most cases, MEP is used as a qualitative approach. The MPE is based on the calculated properties of the charge density directly from the function of the molecular wave and measures the interaction of a positively charged point with the nuclei and electrons of a molecule. The interaction between the molecules occurs between regions of opposite electrostatic potentials. The improved accuracy of PM6 is particularly valuable for generating electronic descriptors for qualitative structure activity relationships (QSAR). When comparing the maps of electrostatic potentials of compounds 1, 2, and 3 ( Figure 3) found mostly in coriander, we can see the presence of ribose in compound 2 in terms of volume, as adenosine creates an important increase in the total volume when we compare it with the compounds 1 and 3. The negative charge concentration of compounds 1 and 2 (red region of the map) is very similar when we observe the pyrimidine and imidazole rings; however, the ribose found in compound 2 also has a high electronic density. Regarding compound 3, tryptophan, the negative charge concentration in the indole ring, and the higher electron density is present on the carbonyl of the carboxylic acid group. In addition, compound 3 has large differences in shape and volume. Regarding the electrostatic potential maps obtained for compounds 4, 5, and 6, it is possible to observe a small variation in relation to the shape and volume. The negative charge distribution is also relatively similar, mainly between compounds 4 and 5, where it is predominant over the carbonyl of isochron-1-one.  The energies of the frontier orbitals, that is, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), were also calculated. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. The energies of the frontier orbitals, that is, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), were also calculated. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. The energies of the frontier orbitals, that is, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), were also calculated. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. The energies of the frontier orbitals, that is, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), were also calculated. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. The energies of the frontier orbitals, that is, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), were also calculated. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. (HOMO) and the lowest unoccupied molecular orbital (LUMO), were also calculated. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. (HOMO) and the lowest unoccupied molecular orbital (LUMO), were also calculated. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. (HOMO) and the lowest unoccupied molecular orbital (LUMO), were also calculated. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals. These results are quantitative descriptive and are widely used, which play an important role in the chemistry reaction and in the formation of several charge complexes. The energy of HOMO is directly related to the potential ionization of the compound and the ability of the molecule to create nucleophiles. The energy of LUMO is directly related to electronic affinity, characterized by the susceptibility of the compound to nucleophiles in relation to them. The difference between the energies of the HOMO-LUMO orbitals, called GAP, is an important indicator of molecular stability. Molecules with a low GAP value are generally reactive, while molecules with a high GAP value indicate high stability and low possibility to react with another compounds [17]. Table 3 presents the values calculated for the HOMO and LUMO orbitals of compounds 1-6, as well as the GAP value. The electronic density of frontier orbitals is a useful way to the detailed characterization of donor-acceptor interactions and the most of chemistry reactions occur at the location of the highest electronic density in the frontier orbitals.  Considering that the higher the energy of HOMO, the greater the electron-donor capacity, and that the lower the energy of LUMO, the lower the resistance to accept electrons, we can define that among the proposed new compounds, compound 6 has the highest GAP value and may have high molecular stability. Compound 2, which has the lowest GAP value, can be considered the least stable compound. Thus, it is possible to predict that compound 6 may demonstrate greater biological activity; however, this is not an indicator of low reactivity for other compounds. The pharmacokinetic properties of the compounds were evaluated, and the SwissADME toolbar was employed. The Swiss ADMET Predictor is a designed software to estimating pharmacokinetic parameters/properties of drug-like compounds from their molecular structures [18]. Lipinski proposed four ADMET properties called the "Rule of Five". This rule of five was the authentic and "most well-known rule-based filter" of drug-likeness, which is used to examine if the compound can be well absorbed orally or not. The rule of five includes: molecular weight (MW) ≤ 500; octanol/water partition coefficient (iLOGP) ≤ 5; number of hydrogen bond donors (HBDs) ≤ 5 (accounted in function of NH or OH groups in the molecule); and number of hydrogen bond acceptors (HBAs) ≤ 10 (accounted in function of N or O atoms in the molecule). Under the Rule of Five, a molecule can only be orally Considering that the higher the energy of HOMO, the greater the electron-donor capacity, and that the lower the energy of LUMO, the lower the resistance to accept electrons, we can define that among the proposed new compounds, compound 6 has the highest GAP value and may have high molecular stability. Compound 2, which has the lowest GAP value, can be considered the least stable compound. Thus, it is possible to predict that compound 6 may demonstrate greater biological activity; however, this is not an indicator of low reactivity for other compounds. The pharmacokinetic properties of the compounds were evaluated, and the SwissADME toolbar was employed. The Swiss ADMET Predictor is a designed software to estimating pharmacokinetic parameters/properties of drug-like compounds from their molecular structures [18]. Lipinski proposed four ADMET properties called the "Rule of Five". This rule of five was the authentic and "most well-known rule-based filter" of drug-likeness, which is used to examine if the compound can be well absorbed orally or not. The rule of five includes: molecular weight (MW) ≤ 500; octanol/water partition coefficient (iLOGP) ≤ 5; number of hydrogen bond donors (HBDs) ≤ 5 (accounted in function of NH or OH groups in the molecule); and number of hydrogen bond acceptors (HBAs) ≤ 10 (accounted in function of N or O atoms in the molecule). Under the Rule of Five, a molecule can only be orally Considering that the higher the energy of HOMO, the greater the electron-donor capacity, and that the lower the energy of LUMO, the lower the resistance to accept electrons, we can define that among the proposed new compounds, compound 6 has the highest GAP value and may have high molecular stability. Compound 2, which has the lowest GAP value, can be considered the least stable compound. Thus, it is possible to predict that compound 6 may demonstrate greater biological activity; however, this is not an indicator of low reactivity for other compounds. The pharmacokinetic properties of the compounds were evaluated, and the SwissADME toolbar was employed. The Swiss ADMET Predictor is a designed software to estimating pharmacokinetic parameters/properties of drug-like compounds from their molecular structures [18]. Lipinski proposed four ADMET properties called the "Rule of Five". This rule of five was the authentic and "most well-known rule-based filter" of drug-likeness, which is used to examine if the compound can be well absorbed orally or not. The rule of five includes: molecular weight (MW) ≤ 500; octanol/water partition coefficient (iLOGP) ≤ 5; number of hydrogen bond donors (HBDs) ≤ 5 (accounted in function of NH or OH groups in the molecule); and number of hydrogen bond acceptors (HBAs) ≤ 10 (accounted in function of N or O atoms in the molecule). Under the Rule of Five, a molecule can only be orally 6 Considering that the higher the energy of HOMO, the greater the electron-donor capacity, and that the lower the energy of LUMO, the lower the resistance to accept electrons, we can define that among the proposed new compounds, compound 6 has the highest GAP value and may have high molecular stability. Compound 2, which has the lowest GAP value, can be considered the least stable compound. Thus, it is possible to predict that compound 6 may demonstrate greater biological activity; however, this is not an indicator of low reactivity for other compounds. The pharmacokinetic properties of the compounds were evaluated, and the SwissADME toolbar was employed. The Swiss ADMET Predictor is a designed software to estimating pharmacokinetic parameters/properties of drug-like compounds from their molecular structures [18]. Lipinski proposed four ADMET properties called the "Rule of Five". This rule of five was the authentic and "most well-known rule-based filter" of drug-likeness, which is used to examine if the compound can be well absorbed orally or not. The rule of five includes: molecular weight (MW) ≤ 500; octanol/water partition coefficient (iLOGP) ≤ 5; number of hydrogen bond donors (HBDs) ≤ 5 (accounted in function of NH or OH groups in the molecule); and number of hydrogen bond acceptors (HBAs) ≤ 10 (accounted in function of N or O atoms in the molecule). Under the Rule of Five, a molecule can only be orally Considering that the higher the energy of HOMO, the greater the electron-donor capacity, and that the lower the energy of LUMO, the lower the resistance to accept electrons, we can define that among the proposed new compounds, compound 6 has the highest GAP value and may have high molecular stability. Compound 2, which has the lowest GAP value, can be considered the least stable compound. Thus, it is possible to predict that compound 6 may demonstrate greater biological activity; however, this is not an indicator of low reactivity for other compounds. The pharmacokinetic properties of the compounds were evaluated, and the SwissADME toolbar was employed. The Swiss ADMET Predictor is a designed software to estimating pharmacokinetic parameters/properties of drug-like compounds from their molecular structures [18]. Lipinski proposed four ADMET properties called the "Rule of Five". This rule of five was the authentic and "most well-known rule-based filter" of drug-likeness, which is used to examine if the compound can be well absorbed orally or not. The rule of five includes: molecular weight (MW) ≤ 500; octanol/water partition coefficient (iLOGP) ≤ 5; number of hydrogen bond donors (HBDs) ≤ 5 (accounted in function of NH or OH groups in the molecule); and number of hydrogen bond acceptors (HBAs) ≤ 10 (accounted in function of N or O atoms in the molecule). Under the Rule of Five, a molecule can only be orally Considering that the higher the energy of HOMO, the greater the electron-donor capacity, and that the lower the energy of LUMO, the lower the resistance to accept electrons, we can define that among the proposed new compounds, compound 6 has the highest GAP value and may have high molecular stability. Compound 2, which has the lowest GAP value, can be considered the least stable compound. Thus, it is possible to predict that compound 6 may demonstrate greater biological activity; however, this is not an indicator of low reactivity for other compounds. The pharmacokinetic properties of the compounds were evaluated, and the SwissADME toolbar was employed. The Swiss ADMET Predictor is a designed software to estimating pharmacokinetic parameters/properties of drug-like compounds from their molecular structures [18]. Lipinski proposed four ADMET properties called the "Rule of Five". This rule of five was the authentic and "most well-known rule-based filter" of drug-likeness, which is used to examine if the compound can be well absorbed orally or not. The rule of five includes: molecular weight (MW) ≤ 500; octanol/water partition coefficient (iLOGP) Considering that the higher the energy of HOMO, the greater the electron-donor capacity, and that the lower the energy of LUMO, the lower the resistance to accept electrons, we can define that among the proposed new compounds, compound 6 has the highest GAP value and may have high molecular stability. Compound 2, which has the lowest GAP value, can be considered the least stable compound. Thus, it is possible to predict that compound 6 may demonstrate greater biological activity; however, this is not an indicator of low reactivity for other compounds. The pharmacokinetic properties of the compounds were evaluated, and the SwissADME toolbar was employed. The Swiss ADMET Predictor is a designed software to estimating pharmacokinetic parameters/properties of drug-like compounds from their molecular structures [18]. Lipinski proposed four ADMET properties called the "Rule of Five". This rule of five was the authentic and "most well-known rule-based filter" of drug-likeness, which is used to examine if the compound can be well absorbed orally or not. The rule of five includes: molecular weight (MW) ≤ 500; octanol/water partition coefficient (iLOGP) ≤ 5; number of hydrogen bond donors (HBDs) ≤ 5 (accounted in function of NH or OH groups in the molecule); and number of hydrogen bond acceptors (HBAs) ≤ 10 (accounted in function of N or O atoms in the molecule). Under the Rule of Five, a molecule can only be orally active/absorb if it does not violate any two or more rules [19]. Table 4 represented some of the ADMET properties/parameters for the compounds 1-6. Physicochemical properties calculated on SwissADME: MW: molecular weight; HBDs: hydrogen bonding donor; HBAs: hydrogen bonding acceptor; iLogP: octanol/water partition coefficient; nV: number of violations; TPSA: total polar surface area; LogS: coefficient of solubility determined by the ESOL method; Class: insoluble < -10 < poor < -6 < moderately < -4 < soluble < -2 < very < 0 < highly; %ABS was expressed by the equation %ABS = 109 -(0.345 × TPSA).
According to the results obtained through the theoretical studies, all compounds obey the standards of the Lipinski Rule, thus indicating a good oral bioavailability during the administration of this substance when in an oral pharmaceutical form. The percentage of absorption showed results for all compounds between 60.85 and 92.20%, indicating that these compounds have a good permeability in the cellular plasmatic membrane. The LogS of many substances already approved for use in therapeutics have a value greater than −4.00; in this study the compounds showed a value between −0.68 and 3.44, and the compounds 1, 2, 3, and 4 were considered soluble, while compounds 5 and 6 were considered moderately soluble. The aqueous solubility of a compound significantly affects its absorption and distribution characteristics. Typically, a low solubility goes along with a bad absorption and, therefore, the general aim is to avoid poorly soluble compounds. The score of all compounds in relation to the drug score was determined by combining records of similarity with already approved drugs, as lipophilicity, solubility, molecular mass, and toxicity risks, and a single numeric value was given, which ranges from 0.0 to 1.0 and can be used to predict the global potential of a compound as a new drug candidate.

Antioxidant Activity
The free radical scavenging activity of the extracts were measured by the radical scavenging ability, using 1.1-diphenyl-2-picrylhydrazil (DPPH). The quantification of the sample was performed, with an EC 50 of 27.27 ± 1.21 µg/mL.

CSEE Decreased the Viability of SH-SY5Y Cells, and Induces Apoptosis via the Intrinsic Pathway
In order to verify the cytotoxic activity of CSEE, the IC50 of the extract was initially determined (Figure 4). Thus, 20 µg/mL of CSEE was defined as the study concentration.  Aiming to identify cell death modalities induced by CSEE in SH-SY5Y cells were stained with FITC-labeled AnnexinV, which binds to PS, and PI, a small fluorescent dye that is excluded from live cells but stains the nuclei of cells with a ruptured or otherwise permeabilized plasma membrane. Hence, the combination of AnnexinV/PI cell staining was used to distinguish apoptotic and necrotic cell death ( Figure 5A), CSEE (20 µg/mL, 24 h), decreased SH-SY5Y cell viability ( Figure 5B), and increased the percentage of cells in early or late apoptosis after 24 h of incubation ( Figure 5C,D, respectively). However, there was no significant increase in cell death by necrosis when compared with untreated cells ( Figure 5E). Aiming to identify cell death modalities induced by CSEE in SH-SY5Y cells were stained with FITC-labeled AnnexinV, which binds to PS, and PI, a small fluorescent dye that is excluded from live cells but stains the nuclei of cells with a ruptured or otherwise permeabilized plasma membrane. Hence, the combination of AnnexinV/PI cell staining was used to distinguish apoptotic and necrotic cell death ( Figure 5A), CSEE (20 µg/mL, 24 h), decreased SH-SY5Y cell viability ( Figure 5B), and increased the percentage of cells in early or late apoptosis after 24 h of incubation ( Figure 5C,D, respectively). However, there was no significant increase in cell death by necrosis when compared with untreated cells ( Figure 5E).
Aiming to identify cell death modalities induced by CSEE in SH-SY5Y cells were stained with FITC-labeled AnnexinV, which binds to PS, and PI, a small fluorescent dye that is excluded from live cells but stains the nuclei of cells with a ruptured or otherwise permeabilized plasma membrane. Hence, the combination of AnnexinV/PI cell staining was used to distinguish apoptotic and necrotic cell death ( Figure 5A), CSEE (20 µg/mL, 24 h), decreased SH-SY5Y cell viability ( Figure 5B), and increased the percentage of cells in early or late apoptosis after 24 h of incubation ( Figure 5C,D, respectively). However, there was no significant increase in cell death by necrosis when compared with untreated cells ( Figure 5E).  An evaluation of the effects of CSEE treatment on the caspase-3 activity was determined and it was possible to observe that CSEE (20 µg/mL) induced a significant regulation of cleaved caspase-3 in treated samples in almost half of the treated cells ( Figure 6), which strongly suggests apoptosis induction properties for CSEE. An evaluation of the effects of CSEE treatment on the caspase-3 activity was determined and it was possible to observe that CSEE (20 µg/mL) induced a significant regulation of cleaved caspase-3 in treated samples in almost half of the treated cells ( Figure 6), which strongly suggests apoptosis induction properties for CSEE.

Bcl-2 and Bax Expression in SH-SY5Y Cells
We investigated the relationship between Bcl-2 and Bax expression in SH-SY5Y cells. Figure 7A,B shows the effects of CSEE on the protein and mRNA levels of Bcl-2 and Bax. The protein and mRNA levels of Bcl-2, which promote anti-apoptosis effects, showed significant reduction in the CSEE-treated cells compared with the control group (non-

Bcl-2 and Bax Expression in SH-SY5Y Cells
We investigated the relationship between Bcl-2 and Bax expression in SH-SY5Y cells. Figure 7A,B shows the effects of CSEE on the protein and mRNA levels of Bcl-2 and Bax. The protein and mRNA levels of Bcl-2, which promote anti-apoptosis effects, showed significant reduction in the CSEE-treated cells compared with the control group (non-treated cells); in contrast, the levels of Bax protein and mRNA, which promote pro-apoptosis effects, were increased in the CSEE-treated cells. Figure 6. Effects of CSEE on caspase-3 activity. The SH-SY5Y cells were treated with 20 µg/mL of CSEE for 24 h. The results showed significant increase in caspase-3 activity on CSEE-treated cells vs. control group (non-treated cells) * p < 0.05. Data shown are representative of three independent experiments. The values are expressed as mean ± SEM.

Bcl-2 and Bax Expression in SH-SY5Y Cells
We investigated the relationship between Bcl-2 and Bax expression in SH-SY5Y cells. Figure 7A,B shows the effects of CSEE on the protein and mRNA levels of Bcl-2 and Bax. The protein and mRNA levels of Bcl-2, which promote anti-apoptosis effects, showed significant reduction in the CSEE-treated cells compared with the control group (nontreated cells); in contrast, the levels of Bax protein and mRNA, which promote proapoptosis effects, were increased in the CSEE-treated cells.

Discussion
The benefits of medicinal plants have been reported by literature as an important support to the quality of life of patients with cancer when used in association with conventional antineoplastic therapies [20]. Thus, the need for further research in order to

Discussion
The benefits of medicinal plants have been reported by literature as an important support to the quality of life of patients with cancer when used in association with conventional antineoplastic therapies [20]. Thus, the need for further research in order to investigate potential candidates for new therapeutic applications may start with its popular use [21].
Drug therapy has long known of the importance of screening secondary metabolites in the development of new drugs, and screenings already carried out show that the phytochemical profile of C. sativum indicates the presence of polyphenols, terpenes, reducing sugars, alkaloids, fatty acids, and sterols [22].
In this study, we verified the presence of flavonoids and phenolic compounds in the extract of C. sativum, which was then confirmed by mass spectrometry.
Isocoumarins and nitrogen compounds were also found. Coriandrone A and B were isolated from C. sativum along with other isocoumarins (coriandrin and dihydrocoriandrin) and their structures, which were determined by spectroscopic means and X-ray crystallography [24]. Herein, the compounds identified at CSEE through UPLC/DAD-ESI/HRMS/MS analysis confirm these findings, evidencing the presence of Isofuranocoumarins-coriandrin and coriandrones A and B [25,26]. Besides that, studies have demonstrated the presence of these compounds in extracts of C. sativum leaves, especially the ethyl acetate fractions, in which the contents of dihydrocoriandrin (34.5%) and coriandrin (14.4%) were evaluated [27]. Products derived from C. sativum, as essential oil and extracts, exhibit antiproliferative [10,28,29] and antioxidant activities [30,31].
To verify the antioxidant activity of CSEE, an antioxidant activity test was performed using the DDPH method. The quantification of the sample was performed, with an EC 50 of 27.27 ± 1.21 µg/mL [4], using the Supercrit Fluid Extraction (SFE, 300 bar, 40 • C) to extract C. sativum compounds, and obtained an EC 50 28.71 ± 2.41 µg/mL, similar to our results. However, the antioxidant activity of a methanolic extract and ethanol: water (80:20) from C. sativum was evaluated, with EC 50 values of 2.2 ± 0.2 and 5.6 ± 0.6, respectively, which are much lower than ours. The authors reported that production procedures and climatic changes, such as average precipitation, harvesting time, altitude, and storage conditions, significantly influence the composition of phytochemicals in plants, directly reflecting on the antioxidant activity [31]. Some authors described that the antioxidant activity of C. sativum could be attributed to the action of (E)-2-decenal and linalool [32,33].
Confirming the CSEE antioxidant activity, assays were performed with SH-SY5Y cells in order to verify the ability of CSEE to cause cell death. Thus, after determining the IC 50 (20 µg/mL), we verified that CSEE-induced apoptosis in SH-SY5Y cells, since the results showed that CSEE promoted greater phosphatidylserine externalization in treated SH-SY5Y cells (20 µg/mL), when compared to untreated cells.
The results obtained suggest that the likely mechanism involved in CSEE apoptotic induction is associated with the intrinsic pathway, where caspase-3 can be activated by a mitochondrial apoptotic pathway, contributing to cell apoptosis [34]. Our results showed that caspase-3 activity was significantly higher in treated SH-SY5Y cells compared to untreated cells. Exploring markers of the intrinsic apoptotic pathway, we observed a reduction in the levels of protein and mRNA of Bcl-2, an anti-apoptotic marker, and an increase in the levels of protein and mRNA Bax, a pro-apoptotic marker in SH -SY5Y cells treated with CSE. Data obtained in the literature describe the cytotoxic activity of C. sativum on SH-SY5Y cells [8,35], although the same cell line has been used to investigate the neuroprotective role of C. sativum [12]. Our results showed in vitro antitumor activity of CSEE, which may be associated with a cytotoxic effect of CSEE, due, at least in part, to the presence of isofurnacoumarins, choriandrin, and choriandrones A and B in CSEE, as demonstrated by UPLC/DAD-ESI /HRMS/MS. Our findings are in agreement with studies that showed antitumor activity of C. sativum against a variety of cells such as MCF-7 [8], K562 and HL60 leukemic cells [36], and HeLa cells [37].
According to the results obtained through the theoretical studies, all compounds indicated a good oral bioavailability during the administration in an oral pharmaceutical form. The score of all compounds was determined by combining records of similarity with already approved drugs, as lipophilicity, solubility, molecular mass, and toxicity risks, and a single numeric value was determined and can be used to predict the global potential of a compound as a new drug candidate.
Nevertheless, even though the in silico methods used in the development of new drugs are promising, it is important to emphasize their limitations, which can sometimes reveal up to 40% differences in pharmacokinetic parameters, such as area under the curve or maximum serum concentration, between observed and simulated data [38,39].

Reagents
Dulbecco's Modified Eagle's Medium (DMEM) and all cell culture reagents were purchased from Life Technologies (CA, USA). Fetal bovine serum (FBS) and MTT were obtained from Himedia. FITC AnnexinV Apoptosis Detection Kit was bought from Becton Dickinson (NJ, USA). All other chemicals were of analytical grade. A colorimetric assay kit was used to measure caspase-3 activity (R&D Systems, Inc., Minneapolis, MN, USA). Enzyme-linked immunonosorbent assay (ELISA) kits were used to measure Bax and Bcl-2 protein levels (Elabscience, Houston, TX, USA). Ethanol solution (60%), for extraction procedures, was prepared using ethanol PA from Synth (Brazil). Sodium carbonate and aluminum chloride were obtained from Dinâmica (Brazil) and Synth, respectively, while Folin-Ciocalteu was obtained from Merck ® . Gallic acid and quercetin standards were purchased from Sigma (Annapolis, MD, USA). Methanol was used for UPLC/MS procedures, acetonitrile and formic acid were HPLC grade from Merck ® , and sodium formate was prepared from sodium hydroxide (Dinâmica, Brazil) and formic acid (Merck ® ).

Preparation of Plant Materials and Phytochemical Analysis
To obtain 500 g of totally dried plant, including roots, leaves, and stems, 3 kg of coriander were used, washed in running water, and dried in an oven between 30 • C and 37 • C for 6 days. After drying, the plant was processed in a knife mill and screened in mesh 80. In order to obtain the fluid C. sativum ethanolic extract-CSEE, 500 g of the dried and grinded plant was moistened with 100mL of 60% ethanol solution, i.e., in a proportion of 20% (w/v), for 48 h. After humidification the solution was transferred to a percolator, with 3 cm of hydrophilic cotton added to the bottom, together with more 100 mL of ethanol solution (60%). The macerate was added and lightly pressed. Filtration paper, porcelain chips, and 500 mL of 60% ethanol were added over the macerate, adjusting the tap to 5 drops/minute, as soon as the column of the macerate was impregnated with the extracting liquid. The percolation was maintained for 36 h, under controlled oven temperature. The fluid extract was subjected to evaporation at 45 • C, and at the end of 5 days of forced evaporation, the crude dry CSEE was collected and subjected to trituration in a mortar and sieved in an 80-mesh sieve, then kept in a dry environment with packages of silica gel [40,41].

Total Soluble Solids
Next, 5.0 mL of fluid CSEE was added to pre-calibrated beakers and was taken to oven to dry at 70 • C. After drying, it was placed in the desiccator for cooling and then being weighed. The amount of soluble solids in the extract was calculated as follows: where: b = weight of the beaker; w = final mass of the extract, after drying; Va = volume of 5 mL.

Estimation of Total Phenolic Content
In a 3 mL vessel we added 200 µL of the CSEE in 1:10 dilution and then added 1400 µL of water. Then, 160 µL of the sodium carbonate solution at 20% and then 240 µL of Folin-Ciocalteu was added, and the sample was shaken for few seconds and incubated in the absence of light for 2 h at room temperature [42]. The sample was shaken for a few seconds again, and the absorbance was measured at a wavelength of 760 nm (Cary-50 spectrophotometer, Varian-Inc, CA, USA). The total amount of phenols (in mg/mL) was calculated by a standard curve prepared with gallic acid in the same conditions of the assay. The procedure was performed in triplicate.

Flavonoids
The CSEE was diluted 1:5 in methanol, and then 100 µL was transferred to the appropriate vessel and diluted in 1500 µL of methanol. Then, 400 µL of aluminum chloride was added in a concentration of 5% and stirred for a few seconds. The solution was incubated for 30 min in the absence of light at room temperature. The sample was agitated again for a few seconds and the absorbance of the sample was measured at wavelength of 425 nm [43,44]. The amount of flavonoids (in mg/mL) was calculated by a standard curve prepared with quercetin in the same conditions of the assay. The procedure was performed in triplicate.

Chemical Composition-Molecular Dereplication (UPLC/DAD-ESI/HRMS/MS)
First, 1.0 mg of CSEE was dissolved in 1 mL of methanol (HPLC grade), filtered at a 0.45 µm modified Polytetrafluoroethylene (PTFE) filter-Millex (Merck-Millipore), and subjected to a UHPLC/DAD-ESI/HRM/SMS analysis, according to the previously described protocol by de Carvalho [45]. A Shimadzu chromatographic system was used with a Kinetex 2.6 µ, C 18 column, as the stationary phase, and kept in a column oven at 55 • C. Acetonitrile and ultrapure water, both with 0.1% of formic acid, were used as the mobile phase. The analytical conditions started from an exploratory gradient with a low concentration of strong elution solvent, 15% (Acetonitrile-ACN) to 95% of this same solvent for 12 min, and then the ACN concentration was kept constant for 4 min, returning the initial concentration after 1 min, and kept constant for 4 min. The UPLC effluent was subjected to electrospray ionization (ESI) and analyzed in a positive mode in high-resolution mass spectrometer, equipped with a QToF mass analyzer (Bruker Daltonics, Billerica, MA, EUA). The temperature of the drying gas was defined as 200 • C at a flow of 9 L min −1 , 2 bar for the pressure of the nebulizer, and 4500 V of capillary voltage (kV). For the characterization of the compounds detected by UHPLC/DAD-ESI/HRMS/MS, fragments from well-resolved chromatographic bands were selected, with a mass-charge ratio (m/z) between 50 and 1200 Da. Sodium formate was used as a calibrator. The data were treated using the DataAnalysis 4.4 software (Bruker), and Extracted Ion Chromatograms (EIC) were generated using the Target Analysis 1.3 software as well as an "in house" list of target candidates, within the compound name and molecular formula, according to C. sativum chemical composition literature and "Coriander Genomics Database" [46,47]. For MS/MS spectra confirmation, fragmentation profiles of identified compounds were compared with those available at FOODB database (https://foodb.ca/, accessed on 20 July 2022) and with predicted mass spectra [48]. Additionally, the Global Natural Products Social Molecular Network (GNPS-open access platform) was used. The data were filtered by removing all MS/MS fragment ions within ±17 Da of the precursor m/z. MS/MS spectra were window filtered by choosing only the top 6 fragment ions in the ±50 Da window throughout the spectrum. The precursor ion mass tolerance was set to 0.01 Da and a MS/MS fragment ion tolerance of 0.02 Da. A network was then created where edges were filtered to have a cosine score above 0.6 and more than 5 matched peaks. The spectra in the network were then searched against GNPS' spectral libraries. The library spectra were filtered in the same manner as the input data. All matches kept between network spectra and library spectra were required to have a score above 0.7 and at least 6 matched peaks [16,46].

Molecular Interactions and ADME Properties Prediction
Structural optimization was performed by the PM6 methods implemented in the semi-empirical quantum chemistry MOPAC 2016 software package. The files containing the three-dimensional information of each structure as well as the charge information of electrostatic potential were the input files viewed by the JMol software. The frontier orbital energy was calculated by the Chem3D Ultra 8.0. Data of orbital energies, dipole moment partial charges, and atomic volumes were used in the comparative analysis by which it was intended to carry out qualitative studies of chemical structure and biological activity (SAR). To evaluate the pharmacokinetic properties of the designed compounds, the 2D structure of the compounds were drawn on Chemdraw Ultra 12.0. Each structure was imported, and the structure smiley was entered at the interface of the website (http://swissadme.ch/, accessed on 20 July 2022). The SwissADME drug design study was run and the ADME properties/parameters were generated [18].

Antioxidant Activity
The activity of CSEE against free radicals was performed using DPPH (1,1-diphenyl-2picrylhydrazyl radical, Sigma, USA). The extract was diluted so that the concentration of soluble solids was equal to 0.1%, and free radical decay curve was performed, where the initial absorbance was approximately 0.7300. The curve was made with eleven points, with 4.9. ELISA (Enzyme Linked ImmunonoSorbent Assay) Bax and Bcl-2 protein levels in SH-SY5Y cells were measured using ELISA kit Bax (Elabscience, Houston, USA) and Bcl-2 (Elabscience, Houston, TX, USA), in accordance with the manufacturer's protocols. Basically, after treatment of SH-SY5Y cells with CSEE (20 µg/mL), samples were added to plates pre-coated with antibodies to Bcl-2 or Bax. After incubation, biotinylated antibodies to Bcl-2 or Bax and avidin horseradish peroxidase conjugate were added and incubated, followed by the addition of the substrate. After incubation, the enzyme-substrate reaction was stopped, followed by spectrophotometric determination (OD450nm), and the values were converted to protein level according to the standard curve.

Reverse Transcription-Quantitative PCR (RT-qPCR)
After treatment with 20 µg/mL of CSEE, total RNA extracted from SH-SY5Y cells samples was converted to cDNA using a SuperScript ® III RT kit (Invitrogen, Carlsbad, CA, USA), according to the manufacturer's protocol. The concentration of RNA was detected using a NanoDrop 2000 (Thermo Fisher Scientific, Inc., Waltham, MA, USA). GAPDH was used as the internal control. The thermocycling conditions were as follows: 95 • C for 10 min followed by 35 cycles of 95 • C for 15 s and 55 • C for 40 s. The 2-∆∆Cq method was used to quantify the relative gene expression levels of the target genes. Relative standard curves were generated by serial dilutions and all samples were run in triplicates. For PCR quantities analysis, the primers of Bcl-2, Bax, and GAPDH genes were designed using Allele ID software. Thus, the sense and anti-sense sequences of primers used in qRT-PCR analysis were: Bcl-2 Forward: '5-TTGTGGCCTTCTTTGAGTTCGGTG-3'; and Reverse: '5-GGTGCCGGTTCAGGTACTCAGTCA-3'; Bax Forward: '5-CTGTGCACCAAGGTGCCG GAACT-3'; and Reverse: '5-CACCCTGGTCTTGGATCCAGCCC-3'; GAPDH Forward: '5-ACCCAGAAGACTGTGGATGG-3'; and Reverse: '5-TCTAGACGGCAGGTCAGGTC-3'.

Statistical Analysis
For treatment of the statistical analysis of the data, the results were expressed as mean ± standard deviation, analyzed independently. The tests were analyzed by one-way ANOVA followed by application of the Tukey test using the program Graph Pad version 5.0. Significant differences were considered as p < 0.05.

Conclusions
In this study, we reported for the first time the potential of CSEE as an inducer of cell death in an in vitro cell model of neuroblastoma. The antiproliferative activity of CSEE may be associated with compounds such as coumarins present in the studied extract. Our results suggest a cytotoxic potential of CSEE, in addition to initial pharmacokinetic parameters, obtained by in silico method. Further studies are needed to elucidate the antiproliferative mechanism presented here in order to better understand the cytotoxic potential of C. sativum.  Data Availability Statement: All raw data will be available upon a reasonable justification request.