Ethylene Induction of Non-Enzymatic Metabolic Antioxidants in Matricaria chamomilla

Phytochemical investigations of Matricaria chamomilla L. (Asteraceae) stated the presence of several compounds with an established therapeutic and antioxidant potential. The chamomile non-enzymatic antioxidant system includes low molecular mass compounds, mainly polyphenols such as cinnamic, hydroxybenzoic and chlorogenic acids, flavonoids and coumarins. The objective of this work was to evaluate the role of the non-enzymatic antioxidant system after stimulation by ethylene in tetraploid chamomile plants. Seven days of ethylene treatment significantly increased the activity of phenylalanine ammonia-lyase, which influenced the biosynthesis of protective polyphenols in the first step of their biosynthetic pathway. Subsequently, considerable enhanced levels of phenolic metabolites with a substantial antioxidant effect (syringic, vanillic and caffeic acid, 1,5-dicaffeoylquinic acid, quercetin, luteolin, daphnin, and herniarin) were determined by HPLC-DAD-MS. The minimal information on the chlorogenic acids function in chamomile led to the isolation and identification of 5-O-feruloylquinic acid. It is accumulated during normal conditions, but after the excessive effect of abiotic stress, its level significantly decreases and levels of other caffeoylquinic acids enhance. Our results suggest that ethephon may act as a stimulant of the production of pharmaceutically important non-enzymatic antioxidants in chamomile leaves and thus, lead to an overall change in phytochemical content and therapeutic effects of chamomile plants, as well.


Introduction
In plants, the non-physiologic over-accumulation of reactive oxygen species (ROS) is related to the effect of biotic or abiotic stimuli. Simultaneously, other physiological defence mechanisms are triggered, which help the normal functioning of plants during stressful conditions. Apart from the enzymatic antioxidants (e.g., superoxide dismutase, peroxidase, or catalase) the non-enzymatic antioxidant system, which includes low molecular mass compounds such as polyphenols (phenolic acids, flavonoids, anthocyanins, lignans, and stilbenes), carotenoids and vitamins (vitamins C and E) is activated, as well [1,2].
From the group of chlorogenic acids, feruloylquinic and mono-, di-, and tri-caffeoylquinic derivatives have been postulated [10]. Information dealing with their specific role in chamomile plants is low, due to the missing detailed identification in plant extracts and subsequent monitoring of changes in their content during induced physiological reactions. The connection of caffeoylquinic acid derivatives in plant response to abiotic stimuli was described, e.g., the effect of short-term UV-B radiation [11] or nutrient imbalance [12]. Thanks to their polyphenol character, chlorogenic acid and its derivatives have the antioxidative potential [13], so they can participate in the non-enzymatic detoxification of oxygen radicals and lipid peroxidation products.
Ethylene is a plant hormone regulating the growth, development, and aging of plants. ROS-generated reagents are active during stress conditions, e.g., cold, salinity, or UV-B radiation. They induce the transcription of genes participating in the ethylene signaling pathway (such as ERF1 as well as LeEILs). Moreover, the ethylene-inducible production of transcription factor ERF1 modulates plant tolerance to abiotic stress, which enhances the transcription of antioxidant enzymes genes coding catalase or glutathione peroxidase [14,15].
The aim of this paper was to evaluate the effect of ethylene on the protective non-enzymatic system of tetraploid chamomile plants, including the accumulation of essential antioxidants, such as flavonoids, cinnamic acids derivatives, and coumarin-like compounds. Based on the ethephon induced changes in the accumulation of chlorogenic acid derivatives and due to the detailed identification of 5-O-feruloylquinic acid (5-FQA) isolated from chamomile, offered in this paper, we could determine the role of individual chlorogenic acids in chamomile plant physiology, which had been evaluated at minimum up to now. A comparison of their antioxidative activity and evaluation of changes in their ethylene-induced accumulation can support a hypothesis that 5-FQA represents a storage form of chlorogenic acid derivatives in chamomile plants. In addition, an evaluation of the major chemical constituents by quantum chemical studies (DFT B3LYP/3-31G**), such as maps of molecular electrostatic potential (MEP), frontier orbital's (HOMO and LUMO) of chlorogenic acid, 1,5-DCQA, and 5-FQA was conducted. The general scheme summarizing the methodological steps is presented in Figure 1.
Molecules 2020, 25, 5720 2 of 16 a lower content of polyphenols than the alcohol extracts (methanol or ethanol), so it proves a lower antioxidant activity [3,4]. The main constituents present in alcohol extracts are flavonoids (apigenin, luteolin, quercetin and related derivatives), dicycloethers [5], phenolic acids and their derivatives, e.g., 2-β-D-glucopyranosyloxy-4-methoxycinnamic acids (GMCA) [6], chlorogenic acid derivatives [7], and coumarins [8,9]. From the group of chlorogenic acids, feruloylquinic and mono-, di-, and tri-caffeoylquinic derivatives have been postulated [10]. Information dealing with their specific role in chamomile plants is low, due to the missing detailed identification in plant extracts and subsequent monitoring of changes in their content during induced physiological reactions. The connection of caffeoylquinic acid derivatives in plant response to abiotic stimuli was described, e.g., the effect of short-term UV-B radiation [11] or nutrient imbalance [12]. Thanks to their polyphenol character, chlorogenic acid and its derivatives have the antioxidative potential [13], so they can participate in the non-enzymatic detoxification of oxygen radicals and lipid peroxidation products.
Ethylene is a plant hormone regulating the growth, development, and aging of plants. ROSgenerated reagents are active during stress conditions, e.g., cold, salinity, or UV-B radiation. They induce the transcription of genes participating in the ethylene signaling pathway (such as ERF1 as well as LeEILs). Moreover, the ethylene-inducible production of transcription factor ERF1 modulates plant tolerance to abiotic stress, which enhances the transcription of antioxidant enzymes genes coding catalase or glutathione peroxidase [14,15].
The aim of this paper was to evaluate the effect of ethylene on the protective non-enzymatic system of tetraploid chamomile plants, including the accumulation of essential antioxidants, such as flavonoids, cinnamic acids derivatives, and coumarin-like compounds. Based on the ethephon induced changes in the accumulation of chlorogenic acid derivatives and due to the detailed identification of 5-O-feruloylquinic acid (5-FQA) isolated from chamomile, offered in this paper, we could determine the role of individual chlorogenic acids in chamomile plant physiology, which had been evaluated at minimum up to now. A comparison of their antioxidative activity and evaluation of changes in their ethylene-induced accumulation can support a hypothesis that 5-FQA represents a storage form of chlorogenic acid derivatives in chamomile plants. In addition, an evaluation of the major chemical constituents by quantum chemical studies (DFT B3LYP/3-31G**), such as maps of molecular electrostatic potential (MEP), frontier orbital's (HOMO and LUMO) of chlorogenic acid, 1,5-DCQA, and 5-FQA was conducted. The general scheme summarizing the methodological steps is presented in Figure 1.

Results and Discussion
Ethephon (ET) acts as a plant growth regulator converted into ethylene in plant metabolism. Ethylene is involved in most plant defensive responses induced by abiotic stress [16,17] or wounding [18]. Its elevated production affects many physiological processes in plants.
We mainly focused on the evaluation of the non-enzymatic antioxidant metabolite system, which occurred in leaves treated by 0.01% ET solution in triple spraying. Previously published results [19] using the same ET concentration point to the significant changes in the accumulation of several coumarins, total soluble phenols, and flavonoids. However, the work does not specifically address phenolic acids or flavonoids, which are more effective antioxidants. For this reason, we chose the given concentration and focused in detail on the selected metabolites. The HPLC-DAD-MS analysis of chamomile leaf methanol extract determined an induced accumulation of flavonoids aglycons, simple cinnamic acids, chlorogenic acid derivatives, and some coumarins with their related compounds (Figures 2 and 3).

Results and Discussion
Ethephon (ET) acts as a plant growth regulator converted into ethylene in plant metabolism. Ethylene is involved in most plant defensive responses induced by abiotic stress [16,17] or wounding [18]. Its elevated production affects many physiological processes in plants.
We mainly focused on the evaluation of the non-enzymatic antioxidant metabolite system, which occurred in leaves treated by 0.01% ET solution in triple spraying. Previously published results [19] using the same ET concentration point to the significant changes in the accumulation of several coumarins, total soluble phenols, and flavonoids. However, the work does not specifically address phenolic acids or flavonoids, which are more effective antioxidants. For this reason, we chose the given concentration and focused in detail on the selected metabolites. The HPLC-DAD-MS analysis of chamomile leaf methanol extract determined an induced accumulation of flavonoids aglycons, simple cinnamic acids, chlorogenic acid derivatives, and some coumarins with their related compounds (Figures 2 and 3).

Results and Discussion
Ethephon (ET) acts as a plant growth regulator converted into ethylene in plant metabolism. Ethylene is involved in most plant defensive responses induced by abiotic stress [16,17] or wounding [18]. Its elevated production affects many physiological processes in plants.
We mainly focused on the evaluation of the non-enzymatic antioxidant metabolite system, which occurred in leaves treated by 0.01% ET solution in triple spraying. Previously published results [19] using the same ET concentration point to the significant changes in the accumulation of several coumarins, total soluble phenols, and flavonoids. However, the work does not specifically address phenolic acids or flavonoids, which are more effective antioxidants. For this reason, we chose the given concentration and focused in detail on the selected metabolites. The HPLC-DAD-MS analysis of chamomile leaf methanol extract determined an induced accumulation of flavonoids aglycons, simple cinnamic acids, chlorogenic acid derivatives, and some coumarins with their related compounds (Figures 2 and 3).

Effect of Ethephon on PAL Activity and Phenylpropanoids Accumulation
All the evaluated components of the non-enzymatic antioxidant system are related to the phenylpropanoid biosynthetic pathway, in which phenylalanine ammonia-lyase (PAL) has a crucial position [20]. In our work, the used ET treatment stimulated the PAL activity in leaves ( Figure 4). A correlation of PAL activity and ethylene has been already described in the work of Heredia et al. [21], in which the applied ethylene stimulated an increase in both the enzyme activity and antioxidants content, specifically 1,5-dicaffeoylquinic acid (1,5-DCQA), ferulic acid, and isocoumarin.

Effect of Ethephon on PAL Activity and Phenylpropanoids Accumulation
All the evaluated components of the non-enzymatic antioxidant system are related to the phenylpropanoid biosynthetic pathway, in which phenylalanine ammonia-lyase (PAL) has a crucial position [20]. In our work, the used ET treatment stimulated the PAL activity in leaves ( Figure 4). A correlation of PAL activity and ethylene has been already described in the work of Heredia et al. [21], in which the applied ethylene stimulated an increase in both the enzyme activity and antioxidants content, specifically 1,5-dicaffeoylquinic acid (1,5-DCQA), ferulic acid, and isocoumarin. One of the first products of the phenylpropanoid pathway after the formation of cinnamic acid are the active forms of its derivatives, such as p-coumaric, caffeic, sinapic, or ferulic acid [22]. Our results in Table 1 show an increased production of the mentioned phenolic acids, already 24 h after the last treatment, except for the p-coumaric acid. The most significant increment was recorded for syringic acid (about 170%). Other substantial elevations were monitored for the vanillic (50%) and caffeic (34%) acid. Ferulic acid acts as a precursor of vanillic and syringic acid [23], so its slightly increased pool may correlate to the significantly elevated levels of both related acids. Both of them expressed significant antioxidative properties, and they are useful as free radical scavengers [24]. The reason for p-coumaric acid's decreased content correlates with its precursor role in the biosynthesis of caffeic acid and its related derivatives (feulic, vanillic and syringic acid), in which the content was enhanced as mentioned above.  One of the first products of the phenylpropanoid pathway after the formation of cinnamic acid are the active forms of its derivatives, such as p-coumaric, caffeic, sinapic, or ferulic acid [22]. Our results in Table 1 show an increased production of the mentioned phenolic acids, already 24 h after the last treatment, except for the p-coumaric acid. The most significant increment was recorded for syringic acid (about 170%). Other substantial elevations were monitored for the vanillic (50%) and caffeic (34%) acid. Ferulic acid acts as a precursor of vanillic and syringic acid [23], so its slightly increased pool may correlate to the significantly elevated levels of both related acids. Both of them expressed significant antioxidative properties, and they are useful as free radical scavengers [24]. The reason for p-coumaric acid's decreased content correlates with its precursor role in the biosynthesis of caffeic acid and its related derivatives (feulic, vanillic and syringic acid), in which the content was enhanced as mentioned above.
In the methanolic extract, mainly inactive glycosides are presented, and related aglycones are the minority. The mutual ratio of glycosides and aglycons amounts indicates the actual physiological state of the plant, since inactive glycosides are considered as a storage form with minimal physiologic effects [6]. Flavonoids in chamomile leaves are usually determined as a sum of flavonols [11,19], while their detection in flowers or foods is more specified. Kováčik and Klejdus [12] offered a more detailed determination of flavonoids in fresh chamomile leaves. The authors established three flavonols (hyperoside, rutin, astragalin) and three flavones glycosides (luteolin-7-O-glucoside, luteolin-3 ,7-O-diglucoside, and apigenin-7-O-glucoside) in hydroponically cultivated plants affected by nutrient stress. From our results, the qualitative and quantitative determination of free quercetin, luteolin, and kaempferol aglycons in ET-treated leaves indicates the mobilization of active forms in conditions simulating stress stimuli. During normal conditions, their levels seemed to be in the following order: Quercetin < kaempferol < luteolin ( Figure 5). The accumulation of quercetin and luteolin arose substantially 24 h from the last ET-treatment and simultaneously, the kaempferol level decreased. This state was maintained during the experiment. In general, quercetin is the most abundant flavonoid in plants, acting as an effective protectant against the harmful effects of ROS [25]. Thanks to the dihydroxylated B-ring of luteolin and quercetin, both could be considered as better antioxidant agents than monohydroxylated flavonoid derivatives (such as kaempferol) [26]. Their ET-stimulated production and antioxidant potential destined them for participation in the active non-enzymatic antioxidant system of chamomile. The used ET treatment had an impact on the accumulation of other main metabolites which occurred in the methanol leaf extract (GMCAs, herniarin and daphnetine glucoside-daphnin), as well (Table 1). Decreased levels of glycosylated cinnamic acid derivatives were in correlation with higher levels of their active forms, such as cinnamic acid and related derivatives with a substantial antioxidant potential. From coumarin-related compounds, daphnin and herniarin were detectable. The production of both after the treatment increased, however, the amount of daphnetin was minimal. In our experiment, umbelliferon was not detected, but its increment in ET-treated chamomile leaves was shown in [19]. The production of free coumarin forms, mainly umbelliferon,  The used ET treatment had an impact on the accumulation of other main metabolites which occurred in the methanol leaf extract (GMCAs, herniarin and daphnetine glucoside-daphnin), as well (Table 1). Decreased levels of glycosylated cinnamic acid derivatives were in correlation with higher levels of their active forms, such as cinnamic acid and related derivatives with a substantial antioxidant potential. From coumarin-related compounds, daphnin and herniarin were detectable. The production of both after the treatment increased, however, the amount of daphnetin was minimal. In our experiment, umbelliferon was not detected, but its increment in ET-treated chamomile leaves was shown in [19]. The production of free coumarin forms, mainly umbelliferon, is connected with the reaction induced by wounding or severe abiotic stress [6,27]. The biosynthesis of coumarins in plants is still under investigation due to their different biological effects. Plants are natural sources of coumarins, but their low levels and instability during isolation lead to the development of other ways of acquisition. Bourgaud et al. [28] suggested the hydroxylation of umbelliferon to dihydroxycoumarins (e.g., daphnetin). However, it has been just confirmed in in vitro experiments. Similarly, Chu et al. [29] successfully produced three different molecules herniarin, esculetin, and skimmin from umbelliferon in in vitro conditions. Based on the results of the mentioned experiments, an increased level of umbelliferon in ET-treated chamomile leaves [19] could be correlated with increased levels of herniarin and daphnin in our results. Chamomile coumarins showed an antioxidant activity in the DPPH assay [30]. The position and type of substituent attached to the aromatic ring of coumarin molecules have a great influence on the antioxidant activity of these metabolites. The higher number of hydroxyl group on the ring structure of coumarins is correlated with higher ROS suppression effects in comparison with methoxylated or non-substituted coumarins.

Role of Chlorogenic Acid Derivatives in Chamomile Plant Physiology and Their Antioxidant Potential
In normal physiological conditions, an average 10-fold higher accumulation of leaf 5-FQA than chlorogenic acid and 1,5-DCQA was shown ( Figure 6). However, 24 h after the ET treatment, the content of 5-FQA significantly decreased with a simultaneous increment in the chlorogenic acid and 1,5-DCQA content. This state was maintained without changes through the next 24 h after the last spraying. Mutual changes of quinic acid derivatives in chamomile have not been described in detail yet. A similar trend in changes of chlorogenic acids accumulation was monitored after the short-term effect of UV-B radiation on chamomile leaves (unpublished data for 5-FQA; Figures A1 and A2 in Appendix A) [11]. Moreover, a significant enhancement of the chlorogenic acid and 1,5-DCQA effect was recorded here. The derivatives of caffeic acid, such as chlorogenic (3-caffeoylquinic acid) and 1,5-DCQA, contain several free hydroxyl groups and express strong antioxidative properties, which were demonstrated in [13,31].
One possible explanation for the enhanced accumulation of caffeoylquinic derivatives in conditions of higher ROS concentration induced by ethephon or UV-B radiation could be their better antioxidant potential in comparison with feruloyl derivatives. Kono et al. [32] declared substantial antioxidant properties of the chlorogenic acid (3-caffeoylquinic acid), and Chan et al. [33] proved the higher effectivity of caffeoylquinic acids to scavenge free radicals than any other phenols. Testing 16 different phenolic compounds, including caffeic and ferulic acids, on the DPPH scavenging activity showed higher effectiveness of caffeic acid to scavenge radicals than ferulic acid [34]. It could be supposed that chlorogenic acids with a higher number of caffeoyl groups could prove a more significant antioxidative capacity and reactivity, and physiological efficacy than feruloyl derivatives. This assumption was demonstrated in the DPPH testing of the antioxidative capacity of standards of chlorogenic acid derivatives confirmed in chamomile plants and isolated 5-FQA (Figure 7). The DPPH assay is widely used to evaluate the free radical scavenging effectiveness, which is the antioxidant activity of various substances. It is expressed as a percentage of reduction of DPPH radicals by testing substances or as an EC50 value, which is the concentration of the extract determining a 50% decrease of absorbance of the DPPH solution [35]. From our results (Figure 7), the highest EC50 value showed 5-FQA, but it reflected the lowest antioxidative capacity (AC) in comparison with the AC of chlorogenic acid and 1,5-DCQA. Yang et al. [36] obtained similar results, which showed a higher AC of metabolites with the caffeoyl group presented in the molecule than these containing only feruloyl. Moreover, they concluded that the antioxidant activity is determined not only by the number of phenolic OH-groups, but also by their position in the aromatic rings. Compounds exhibiting vicinal OH-groups on the aromatic rings, such as dihydrocaffeic acid, are particularly potent antioxidants. radicals by testing substances or as an EC50 value, which is the concentration of the extract determining a 50% decrease of absorbance of the DPPH solution [35]. From our results (Figure 7), the highest EC50 value showed 5-FQA, but it reflected the lowest antioxidative capacity (AC) in comparison with the AC of chlorogenic acid and 1,5-DCQA. Yang et al. [36] obtained similar results, which showed a higher AC of metabolites with the caffeoyl group presented in the molecule than these containing only feruloyl. Moreover, they concluded that the antioxidant activity is determined not only by the number of phenolic OH-groups, but also by their position in the aromatic rings. Compounds exhibiting vicinal OH-groups on the aromatic rings, such as dihydrocaffeic acid, are particularly potent antioxidants.  A detailed identification of 5-FQA contributed to the elucidation of the role of chlorogenic acid derivatives in chamomile plant physiology. In chamomile plants, several quinic acids (such as caffeoyl and feruloyl derivatives) were recorded mostly by LC-MS methods using the related chemical standards [7,10]. Up to now, 5-FQA was determined similarly. Here, its LC isolation and following NMR identification was conducted ( Figure 8). Its high content in normal physiological conditions and decreased amount after the stress indicate its storage role of chamomile quinic acid derivatives. However, for an absolute confirmation of the hypothesis dealing with other changes of chlorogenic acid derivatives, more detailed analyses of participating enzymatic and molecular components of corresponding chamomile biosynthetic pathways are necessary.

Molecular Modelling
The map of the molecular electrostatic potential (MEP) shows the regions which are specifically coloured depending on the potential/charges. The blue/green colour indicates the positive regions, the parts of the molecule with a higher probability to suffer a nucleophilic attack. The negative regions visualized by a green/red colour indicate the regions which are probably performing a nucleophilic attack [37]. MEPs of the studied quinic derivatives are presented in Figure S1  A detailed identification of 5-FQA contributed to the elucidation of the role of chlorogenic acid derivatives in chamomile plant physiology. In chamomile plants, several quinic acids (such as caffeoyl and feruloyl derivatives) were recorded mostly by LC-MS methods using the related chemical standards [7,10]. Up to now, 5-FQA was determined similarly. Here, its LC isolation and following NMR identification was conducted (Figure 8). Its high content in normal physiological conditions and decreased amount after the stress indicate its storage role of chamomile quinic acid derivatives. However, for an absolute confirmation of the hypothesis dealing with other changes of chlorogenic acid derivatives, more detailed analyses of participating enzymatic and molecular components of corresponding chamomile biosynthetic pathways are necessary. A detailed identification of 5-FQA contributed to the elucidation of the role of chlorogenic acid derivatives in chamomile plant physiology. In chamomile plants, several quinic acids (such as caffeoyl and feruloyl derivatives) were recorded mostly by LC-MS methods using the related chemical standards [7,10]. Up to now, 5-FQA was determined similarly. Here, its LC isolation and following NMR identification was conducted (Figure 8). Its high content in normal physiological conditions and decreased amount after the stress indicate its storage role of chamomile quinic acid derivatives. However, for an absolute confirmation of the hypothesis dealing with other changes of chlorogenic acid derivatives, more detailed analyses of participating enzymatic and molecular components of corresponding chamomile biosynthetic pathways are necessary.

Molecular Modelling
The map of the molecular electrostatic potential (MEP) shows the regions which are specifically coloured depending on the potential/charges. The blue/green colour indicates the positive regions, the parts of the molecule with a higher probability to suffer a nucleophilic attack. The negative regions visualized by a green/red colour indicate the regions which are probably performing a The frontier orbitals provide information on the highest occupied molecular orbital (HOMO) energies and the lowest unoccupied molecular orbital (LUMO) energies, in addition to acting as electron-donors and electron-acceptors, respectively. An important electronic parameter for describing the antioxidant ability of the studied molecules is the HOMO energy. Therefore, the higher the HOMO energy, the higher the capacity to realize nucleophilic attacks, i.e., donating electrons and scavenging the radical [37]. 1,5-DCQA is presented by the highest HOMO value, −5.484 eV, followed by 5-FQA (−5.5648 eV) and chlorogenic acid (−5.888 eV). HOMO regions are located mainly over aromatic rings and double bounds of caffeic and ferrulic parts of chlorogenic acid and 5-FQA, respectively. In 1,5-DCQA, HOMO is located only over one aromatic ring of caffeic parts.
On the other hand, the analyses of LUMO orbitals bring information on the region susceptible to suffer nucleophilic attacks. Therefore, the lower the LUMO energy, the lower the resistance to accept electrons. The minimum value was represented by 1,5-DCQA (−1.937 eV). The distribution of LUMO regions is similar to HOMO, over the aromatic rings and double bonds. In 1,5-DCQA, it is mainly located on the aromatic ring of the other caffeic parts.
For each molecule, the gap energy value (∆E = E LUMO − E HOMO ) indicating the chemical reactivity and molecular stability was calculated. Molecules with high HOMO values and a low gap value show higher reactivity and consequently less stability, which predetermines them to have an antioxidant activity by scavenging ROS. Here, the highest HOMO and the lowest gap value (3.92 eV) was observed for 5-FQA, although Urbaniak et al. [38] showed better properties of the caffeic acid derivatives than ferulic acid and its derivatives.
Gallic acid represents a molecule with significant antioxidative properties and therefore, it is widely used as a standard compound in many antioxidative assays. A comparison of our results with the data on gallic acid available in [37] can provide information on the degree of antioxidant effects of the investigated quinic derivatives. The HOMO, LUMO, and gap values for gallic acid were 1.2327, 0.2414, and 1.474 eV, respectively. All these values meet the criteria significantly better, which suggests better antioxidant properties compared to the studied substances.
For in silico analyses of the studied quinic acid derivatives, the Molinspiration Cheminformatics and Swiss Target Prediction servers were used. The molecules considered as a potential orally active compound had to satisfy the Lipinski's rule of five of the drug-like properties. These physicochemical parameters are associated with an acceptable aqueous solubility and intestinal permeability and comprise the first steps in oral bioavailability. For example, the lipophilicity coefficient (LogP) is a property that quantitatively measures the lipophilia of the compounds, one of the most important molecular properties for drug absorption [39]. It was observed that only 5-FQA did not violate the rule of five (Table S1), i.e., ≤ 5 H-bond donors, ≤ 10 H-bond acceptors, < 500 g mol −1 of molecular weight, < 5 logP coefficient, < 10 rotatable bonds, and topological polar surface area (TPSA) < 140. 1,5-DCQA violates three rules due to the higher molecular weight and number of H-bond donors and acceptors. Detailed in silico studies were performed only for the chlorogenic acid (no available data for 5-FQA and 1,5-DCQA), which on the basis of the results, such as identification of its thermodynamically stable conformation or nucleophilic site, appears to be a promising drug [40]. Although the chlorogenic acid violates one rule (number of H-bond donors), its pharmaceutical properties are extensively studied due to the fact that it is the most abundant isomer among caffeoylquinic acid isomers and one of the most available acids among phenolic acid compounds, which can be naturally found in green coffee extracts and tea [41].
The Swiss target prediction estimates the most probable macromolecular targets of the tested molecules, which were assumed as bioactive. In this server, the targets prediction is obtained by the combination of molecular 2D and 3D similarity with a library of 370,000 known actives on more than 3000 proteins from three different species. The Swiss prediction for all three molecules was similar, but the score differed. The highest score was found for two enzymes: Aldose reductase and aldo-keto reductase family 1 member B10, and protease matrix metalloproteinase 2 (see Table S1).

Plant Material and Treatments
Seeds of tetraploid plants of Matricaria chamomilla L. (cv. 'Lutea') were germinated in the sand under laboratory conditions. Three-week-old seedlings were transferred into the pot filled with 130 g of garden substrate with a 60% water-holding capacity. They were cultivated at 25 • C with a 12-h photoperiod under fluorescent white light (TLD 36 W/33 Philips fluorescent tubes with a photon flux density 210 µmol m −2 s −1 ). After the next 5 weeks, plants were divided into two groups: Controls and treated by ethephon (2-chloroethyl phosphonic acid) (Sigma, Canada). Control plants were chemically untreated, cultivated under the same conditions such as the treated ethephon, and collected in time harvesting of the first treated group. During the 7 days, the 0.01% aq. solution of ethephon was triple sprayed on chamomile leaves and plant material was analyzed 24 and 48 h after the last spraying. Ethephon is usually used as an ethylene releaser due to its gradual degradation to ethylene after its contact with the plant surface.

Isolation and Identification of 5-O-feruloylquinic Acid
An ethanolic extract of leaves (50.2 g) was fractionated in column chromatography. Silica gel (70-140 µm) was deactivated by 11% of water and eluted with mobile phase CHCl 3 contained from 5 to 20% of methanol. The desired fraction (1.8 g) was re-chromatographed and eluted with 1% of methanol in CHCl 3 . After that, 40 mg of brown compound was obtained and used for NMR analysis. For the separation and purification of 5-FQA, the preparative HPLC with DAD detection with a column Kromasil100 C18 (250 × 4.6 mm, 5 µm) was used (Agilent Technologies 1260 Infinity device, Santa Clara, CA, USA).
The structure of 5-FQA was determined by 1D ( 1 H, 13 C, selective NOESY, TOCSY) and 2D (gCOSY, NOESY, gHSQCAD, gHMBCAD) NMR spectroscopy (Figure 8). The residue was dissolved in deuterated methanol (CD 3 OD) and chemical shifts were referenced to the residual solvent peak ( 1 H NMR 3.31 ppm and 13 C NMR 49.1 ppm). The NMR spectra were recorded on a 600 MHz Varian VNMRS NMR spectrometer (Agilent, Palo Alto, USA, CA) equipped with a 5 mm One NMR probeat at 25 • C. The melting point of the substance was set at 117-121 • C.
The 7-point calibration curve was constructed by analyzing various amounts of 5-FQA (21.57 mg/100 mL diluted in 75% (v/v) methanol, from which 1, 2, 3, 4, 5, 6, 7 µL were taken and analyzed). Every point of the curve is the result of two-injection analyses. The whole range of tested concentrations showed good linearity with a regression coefficient (R2) better than 0.999. Limits of detection (LOD) and quantification (LOQ) were not determined due to the high metabolite accumulation in plants.
The activity of phenylalanine ammonia-lyase (PAL, EC 4.3.1.5) was determined as the production of trans-cinnamic acid (t-CA) from phenylalanine using the HPLC method [43] in homogenates prepared using a sodium borate buffer (pH 8.8). The HPLC conditions: Kromasil SGX C18 7 µm (150 × 4.6 mm) column, mobile phase was 55% acetonitrile in the flow rate of 0.8 ml min −1 . The detection of t-CA was carried out by the retention time and performed at 275 nm. The activity of PAL was expressed as nmol t-CA min −1 g −1 fresh wt. The t-CA (Sigma-Aldrich, Steinheim, Germany) standard compound was used for the quantification.

Antioxidant Capacity Assay
The antioxidant capacity (AC) of quinic acid derivatives, 5-FQA, chlorogenic acid, and 1,5-dicaffeoylquinic acid (1,5-DCQA), was assayed by measuring the decrease in absorbance at 515 nm of the stable free radical DPPH [34,45]. The purple-coloured free radical reacts with the scavenger to yield the colourless product 1,1-diphenyl-2-picrylhydrazine. From the quinic acid derivatives stock solutions (1 mg of chlorogenic and 1,5-DCQA in 2 mL of 75% methanol; 1.43 mg of 5-FQA in 2 mL of 75% methanol) 100, 50, 25, 12.5, 6.25, 3.125, 1.075 µL were taken and mixed with 100 µL 0.2 mM DPPH. The absorption was recorded after 30 min on the fluorescent reader (Synergy HtBiotek, Bad Friedrichshall, Germany). AC was determined as the percentage of radical scavenging calculated according to the following equation: [%] = 100 × (Ac − As)/Ac Where As is the absorbance of the sample and Ac is the absorbance of the control. The percentage of DPPH scavenging versus the concentration of the sample was plotted. The concentration of the sample necessary to decrease the DPPH concentration by 50% was obtained by interpolation from the linear regression analysis and denoted a EC50 value (µg mL −1 ). All determinations were carried out in triplicate.

Molecular Modelling
Molecular modelling started with the construction of chlorogenic acid, 1,5-DCQA, and 5-FQA in the program GaussView 6.0 [46], followed with computational calculations in the Gaussian 16 program [47]. The compounds were optimized in the density functional theory method (DFT) in theory level B3LYP/6-31G** and frequencies. The constructions of the MEP and the frontier orbitals (HOMO and LUMO) were visualized with the aid of the Avogadro program [48].
The Molinspiration Cheminformatics and Swiss Target Prediction (Lausanne, Swiss) servers were used for determination of in silico predictions.

Statistics
The obtained data were statistically analyzed using ANOVA and Tukey's test in Minitab Release 11 (Minitab Inc., State College, PA, USA). The number of replications (n) in the figures denotes individual plants measured for each parameter.

Conclusions
In summary, the application of ethephon (as a source of ethylene) on chamomile plants influenced their non-enzymatic antioxidant system involving compounds of the phenylpropanoid pathway with high antioxidant potentials, such as hydroxycinnamic, hydroxybenzoic and chlorogenic acids, flavonoids, and coumarins. Using ethylene could be understood as a way to change the phytochemical content and to ameliorate the therapeutic effects of chamomile plants. Moreover, this paper indicated the functions of chamomile chlorogenic acid derivatives. Each of the other changes in chlorogenic acid derivatives levels induced by the abiotic factor has shown a need for a more detailed study of the activity participating in enzymatic and non-enzymatic components of this branch of phenylpropanoid pathway. However, this should be preceded by a molecular analysis and identification. Our results indicated the assumed antioxidant properties of chlorogenic acid derivatives, which are based on molecular modelling including the analysis of MEP, frontiers orbitals (HOMO, LUMO), and in silico predictions using the Swiss Target Prediction and the Molinspiration servers. However, a more detailed determination of antioxidant properties in the sense of violating a rule for oral bioavailability, log P, mainly for 5-FQA and 1,5-DCQA, requires more precise ways of modelling.
Supplementary Materials: The following are available online, Figure S1: Potential electrostatic maps (MEP) of selected quinic acid derivatives; Figure S2: The frontier orbital HOMO and LUMO of selected quinic acid derivatives; Table S1: The predicted molecular properties confined to the drug-like properties (based on Lipinski's rule of five).  Acknowledgments: We thank Rafael Dolezal and Tomas Kucera for their valuable technical assistance in molecular modelling.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
Decreased accumulation of 5-O-feruloylquinic acid after the effect of short-term UV-B dosage in Matricaria chamomilla leaves (unpublished data). Short-term effect of UV-B radiation dosage (comparable to the highest possible eco-physiological daily UV level in Slovakia) on tetraploids chamomile plants increased levels of quinic acid derivatives, namely chlorogenic acid (3) and 1,5-dicaffeoylquinic acid (7) [11]. Metabolite (5) on Figure A1 was not identified during that time, but now, the presented spectral data offer the identification as 5-FQA.
Molecules 2020, 25, 5720 13 of 16 Acknowledgments: We thank Rafael Dolezal and Tomas Kucera for their valuable technical assistance in molecular modelling.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
Decreased accumulation of 5-O-feruloylquinic acid after the effect of short-term UV-B dosage in Matricaria chamomilla leaves (unpublished data). Short-term effect of UV-B radiation dosage (comparable to the highest possible eco-physiological daily UV level in Slovakia) on tetraploids chamomile plants increased levels of quinic acid derivatives, namely chlorogenic acid (3) and 1,5dicaffeoylquinic acid (7) [11]. Metabolite (5) on Figure A1 was not identified during that time, but now, the presented spectral data offer the identification as 5-FQA. Unpublished results obtained from the UV experiment showed a significant decrease in 5-FQA content through the 48 h ( Figure A2). Considering the quinic acid derivatives content changes and after the stress, with the statement of higher antioxidant capacity of caffeoyl derivatives, support the hypothesis of physiological function of 5-FQA as a native storage form of the quinic acid derivatives in chamomile plants. Unpublished results obtained from the UV experiment showed a significant decrease in 5-FQA content through the 48 h ( Figure A2). Considering the quinic acid derivatives content changes and after the stress, with the statement of higher antioxidant capacity of caffeoyl derivatives, support the hypothesis of physiological function of 5-FQA as a native storage form of the quinic acid derivatives in chamomile plants.