Effects of Coumarinyl Schiff Bases against Phytopathogenic Fungi, the Soil-Beneficial Bacteria and Entomopathogenic Nematodes: Deeper Insight into the Mechanism of Action

Coumarin derivatives have been reported as strong antifungal agents against various phytopathogenic fungi. In this study, inhibitory effects of nine coumarinyl Schiff bases were evaluated against the plant pathogenic fungi (Fusarium oxysporum f. sp. lycopersici, Fusarium culmorum, Macrophomina phaseolina and Sclerotinia sclerotiourum). The compounds were demonstrated to be efficient antifungal agents against Macrophomina phaseolina. The results of molecular docking on the six enzymes related to the antifungal activity suggested that the tested compounds act against plant pathogenic fungi, inhibiting plant cell-wall-degrading enzymes such as endoglucanase I and pectinase. Neither compound exhibited inhibitory effects against two beneficial bacteria (Bacillus mycoides and Bradyrhizobium japonicum) and two entomopathogenic nematodes. However, compound 9 was lethal (46.25%) for nematode Heterorhabditis bacteriophora and showed an inhibitory effect against acetylcholinesterase (AChE) (31.45%), confirming the relationship between these two activities. Calculated toxicity and the pesticide-likeness study showed that compound 9 was the least lipophilic compound with the highest aquatic toxicity. A molecular docking study showed that compounds 9 and 8 bind directly to the active site of AChE. Coumarinyl Schiff bases are promising active components of plant protection products, safe for the environment, human health, and nontarget organisms.


Introduction
Fungal plant pathogens are dominant causal agents of plant diseases. Although in modern agriculture most fungal pathogens are controlled by modern crop management, fungal epidemics still occur, resulting in huge yield losses that have significant negative effects on the economy and society [1]. Nine Fusarium species and Microdohium nivale were identified from isolates of wheat plants in Croatia [2], while in eastern Croatia F. graminearum and F. culmorum are dominant (Wm. G. Sm.) [3].
Although the current use of pesticides in agriculture has led to a significant improvement in crop yields, future crop management will require the development of novel active agrochemicals that have ultra-high efficiency and are environmentally safe. Therefore, research on the multiple biological effects of newly synthesized compounds, including a test for the environmental properties, toxicity, and broad-spectrum mode of action, is necessary [4,5].
The development of new pesticides from the synthesis of active compounds to their commercialization is a costly and time-consuming process. Computer-aided molecular and reused [30]. The aim of this study was to evaluate the inhibitory effects of these compounds against plant pathogenic fungi.
However, the application of these compounds as active substances of plant protection products and their residues could threaten the environment and human health. To design efficient and environmentally friendly antifungal agents, their environmental impact was assessed against the soil-beneficial bacteria and nematodes. Their toxicity was estimated using a computer program based on QSAR methodology. Three enzymes responsible for the fungal growth and the three plant cell-wall-degrading enzymes were chosen as targets for molecular docking. To suggest the possible inhibitory mechanism of the tested compounds against plant pathogenic fungi, experimentally obtained results were compared with the results of molecular docking on these six enzymes. In vitro enzyme inhibition study should confirm the relationship of nematicidal activity with the inhibitory capability of AChE. Molecular docking was used to determine the binding affinity and interactions of the most active compounds with AChE and reveal their mode of action.

Results and Discussion
The structures of the analyzed compounds are presented in Table 1. Details of synthesis and characterization of compounds were described previously [30]. The invention of new active ingredients in plant protection products must be highly specific and environmentally and toxicologically acceptable. Crop protection by new chemicals demands extensive environmental impact assessment. Recently, we have synthesized novel coumarinyl Schiff bases using deep eutectic solvents (DESs) as environmentally friendly media. DESs are low-toxic, biodegradable solvents, which can be easily recycled and reused [30]. The aim of this study was to evaluate the inhibitory effects of these compounds against plant pathogenic fungi.
However, the application of these compounds as active substances of plant protection products and their residues could threaten the environment and human health. To design efficient and environmentally friendly antifungal agents, their environmental impact was assessed against the soil-beneficial bacteria and nematodes. Their toxicity was estimated using a computer program based on QSAR methodology. Three enzymes responsible for the fungal growth and the three plant cell-wall-degrading enzymes were chosen as targets for molecular docking. To suggest the possible inhibitory mechanism of the tested compounds against plant pathogenic fungi, experimentally obtained results were compared with the results of molecular docking on these six enzymes. In vitro enzyme inhibition study should confirm the relationship of nematicidal activity with the inhibitory capability of AChE. Molecular docking was used to determine the binding affinity and interactions of the most active compounds with AChE and reveal their mode of action.

Results and Discussion
The structures of the analyzed compounds are presented in Table 1. Details of synthesis and characterization of compounds were described previously [30]. Heterorhabditidae may be attributed to the differences in nematode's AChE concentration [29]. The invention of new active ingredients in plant protection products must be highly specific and environmentally and toxicologically acceptable. Crop protection by new chemicals demands extensive environmental impact assessment. Recently, we have synthesized novel coumarinyl Schiff bases using deep eutectic solvents (DESs) as environmentally friendly media. DESs are low-toxic, biodegradable solvents, which can be easily recycled and reused [30]. The aim of this study was to evaluate the inhibitory effects of these compounds against plant pathogenic fungi.
However, the application of these compounds as active substances of plant protection products and their residues could threaten the environment and human health. To design efficient and environmentally friendly antifungal agents, their environmental impact was assessed against the soil-beneficial bacteria and nematodes. Their toxicity was estimated using a computer program based on QSAR methodology. Three enzymes responsible for the fungal growth and the three plant cell-wall-degrading enzymes were chosen as targets for molecular docking. To suggest the possible inhibitory mechanism of the tested compounds against plant pathogenic fungi, experimentally obtained results were compared with the results of molecular docking on these six enzymes. In vitro enzyme inhibition study should confirm the relationship of nematicidal activity with the inhibitory capability of AChE. Molecular docking was used to determine the binding affinity and interactions of the most active compounds with AChE and reveal their mode of action.

Results and Discussion
The structures of the analyzed compounds are presented in Table 1. Details of synthesis and characterization of compounds were described previously [30].
The results of antifungal, antibacterial, and nematicidal activity evaluation of nine coumarinyl Schiff bases are presented in Table 2. Heterorhabditidae may be attributed to the differences in nematode's AChE concentration [29]. The invention of new active ingredients in plant protection products must be highly specific and environmentally and toxicologically acceptable. Crop protection by new chemicals demands extensive environmental impact assessment. Recently, we have synthesized novel coumarinyl Schiff bases using deep eutectic solvents (DESs) as environmentally friendly media. DESs are low-toxic, biodegradable solvents, which can be easily recycled and reused [30]. The aim of this study was to evaluate the inhibitory effects of these compounds against plant pathogenic fungi.
However, the application of these compounds as active substances of plant protection products and their residues could threaten the environment and human health. To design efficient and environmentally friendly antifungal agents, their environmental impact was assessed against the soil-beneficial bacteria and nematodes. Their toxicity was estimated using a computer program based on QSAR methodology. Three enzymes responsible for the fungal growth and the three plant cell-wall-degrading enzymes were chosen as targets for molecular docking. To suggest the possible inhibitory mechanism of the tested compounds against plant pathogenic fungi, experimentally obtained results were compared with the results of molecular docking on these six enzymes. In vitro enzyme inhibition study should confirm the relationship of nematicidal activity with the inhibitory capability of AChE. Molecular docking was used to determine the binding affinity and interactions of the most active compounds with AChE and reveal their mode of action.

Results and Discussion
The structures of the analyzed compounds are presented in Table 1. Details of synthesis and characterization of compounds were described previously [30]. The results of antifungal, antibacterial, and nematicidal activity evaluation of nine coumarinyl Schiff bases are presented in Table 2.
Heterorhabditidae may be attributed to the differences in nematode's AChE concentration [29]. The invention of new active ingredients in plant protection products must be highly specific and environmentally and toxicologically acceptable. Crop protection by new chemicals demands extensive environmental impact assessment. Recently, we have synthesized novel coumarinyl Schiff bases using deep eutectic solvents (DESs) as environmentally friendly media. DESs are low-toxic, biodegradable solvents, which can be easily recycled and reused [30]. The aim of this study was to evaluate the inhibitory effects of these compounds against plant pathogenic fungi.
However, the application of these compounds as active substances of plant protection products and their residues could threaten the environment and human health. To design efficient and environmentally friendly antifungal agents, their environmental impact was assessed against the soil-beneficial bacteria and nematodes. Their toxicity was estimated using a computer program based on QSAR methodology. Three enzymes responsible for the fungal growth and the three plant cell-wall-degrading enzymes were chosen as targets for molecular docking. To suggest the possible inhibitory mechanism of the tested compounds against plant pathogenic fungi, experimentally obtained results were compared with the results of molecular docking on these six enzymes. In vitro enzyme inhibition study should confirm the relationship of nematicidal activity with the inhibitory capability of AChE. Molecular docking was used to determine the binding affinity and interactions of the most active compounds with AChE and reveal their mode of action.

Results and Discussion
The structures of the analyzed compounds are presented in Table 1. Details of synthesis and characterization of compounds were described previously [30]. The results of antifungal, antibacterial, and nematicidal activity evaluation of nine coumarinyl Schiff bases are presented in Table 2.
The results of antifungal, antibacterial, and nematicidal activity evaluation of nine coumarinyl Schiff bases are presented in Table 2. Table 2. Results of biological activity assays of nine coumarinyl Schiff bases: a inhibition rate %, 48 h after inoculation at the concentration 0.08 µmol/mL; b minimum inhibitory concentration (MIC/µg mL −1 ; c percentage corrected mortality, %, 48 h after inoculation at the concentration 500 µg/mL), and acetylcholinesterase (AChE) inhibition assay. Tested compounds were most active against Macrophomina phaseolina. The inhibition rates 48 h after inoculation varied from 65.17% (8) to 71.51% (3) compared to control ( Table 2). The most active compounds against M. phaseolina, 3 and 4, are pyrazole-based Schiff bases with the aromatic rings as substituents. A study of Guo et al. [19] also showed that Schiff bases with an aromatic nucleus have enhanced antifungal activity. The substitution of bromine (compound 3, 71.51% inhibition) or methoxy group (compound 4, 70.36% inhibition) from position R 1 (Table 1) with the stronger electron-withdrawing group, chlorine atom (compound 2), lowered the inhibition against M. phaseolina to 67.47%, as well as two strong electron-withdrawing groups at the position R 1 and R 2 (compound 1, 67.47%). The most effective compounds against Sclerotinia sclerotiorum were 9 (56.69%) and 4 (55.33%). Compound 9 is Schiff base with 4-methoxybenzylidene. Substitution of 4-methoxy group on benzylidene with the stronger electron-withdrawing group, dimethylamine, lowers the inhibition of compound 8 to 20.49 %. Coumarins with an electron-withdrawing group (NO 2 or acetate) that favor the inhibition have shown an opposite effect on the human pathogen strain Aspergilus [31]. Fungal growth inhibitions against the Fusarium oxysporum f. sp. lycopersici, and F. culmorum were lower than 37.26%. None of the compounds stimulated growth of mycelia of pathogenic fungi. The highest inhibitory effect (30.34%) on F. oxysporum f. sp. lycopersici had compound 5 with chlorine and iodine atoms on two phenyl rings.

Antifungal
Principal component analysis (PCA) was used for easier interpretation of the inhibitory effects of compounds on individual fungi. The results of PCA analysis are presented in Tables 3 and 4, and visualized by biplot ( Figure 1). A biplot represents the scores of the observations (compounds) on the principal components, and it uses vectors (fungi species) to represent the coefficients of the variables on the principal components [32]. PCA analysis shows that the first two components explain 87.20% of the total standard variation. PC 1 is associated with M. phaseolina, F. oxysporum f. sp. Lycopersici, and F. culmorum, while PC 2 is linked to S. sclerotiorum. The biplot of the first two principal components (PCs) reveals that compounds 2, 4, 5, and 7, which are positioned on the top right part of the PCA biplot, influence the two fungi F. oxysporum f. sp. lycopersici and M. phaseolina. F. culmorum is situated at the right bottom part and is correlated with compounds 1 and 3. S. sclerotiorum is an independent vector forming a large angle regarding a vector of F. culmorum, and angles close to 90 • considering vectors of M. phaseolina, Fusarium oxysporum f. sp. lycopersici. Compounds that are arranged on the left top part of biplot (9 and 6) were highly influenced on inhibition of S. sclerotiorum. Compound 8, which is located alone on the left bottom quadrant, did not show significant inhibition of any fungi.   To elucidate the possible mechanism of action of coumarinyl Schiff bases against pathogenic fungi, a molecular docking study was performed on two groups of enzymes: enzymes responsible for the fungal growth (demethylase [21], chitinase [35], transferase [36]) and plant cell-wall-degrading enzymes (endoglucanase I [37]; proteinase K [38]; pectinase [39]). The compounds were ranked by an energy-based scoring function. The results of molecular docking on enzymes related to the antifungal activity are presented in Table  5. Previously studies also have proved the antifungal activity of coumarins and Schiff bases. Chitosan-derived Schiff bases inhibited Colletotrichum lagenarium and Botrytis cinerea mycelial growth by 26-33% and 35-38%, respectively [33]. Chen et al. [34] showed that the inulin derivatives modified by Schiff bases had significantly better antifungal activity against important plant pathogens Botrytis cinerea, Phomopsis asparagi, and Fusarium oxysporum f. sp. cucumerinum compared to activity of pure inulin. In addition, the antifungal activity of Schiff bases of inulin that bears a pyridine ring exhibited good antifungal activities against Botrytis cinerea, but the different position of the N atom on the pyridine ring did not show significant influence on the antifungal activity [20]. The QSAR study of eight substituted coumarin derivatives of the antifungal activity against phytopathogenic fungi Valsa mali revealed that small electron-withdrawing substituents of coumarin's phenyl ring and hydrophilic electron-donating groups on the coumarin's pyrone ring could enhance the antifungal activity [8].
To elucidate the possible mechanism of action of coumarinyl Schiff bases against pathogenic fungi, a molecular docking study was performed on two groups of enzymes: enzymes responsible for the fungal growth (demethylase [21], chitinase [35], transferase [36]) and plant cell-wall-degrading enzymes (endoglucanase I [37]; proteinase K [38]; pectinase [39]). The compounds were ranked by an energy-based scoring function. The results of molecular docking on enzymes related to the antifungal activity are presented in Table 5.
A comparison of the results shown in Table 5 with the antifungal activity against M. phaseolina ( Table 2), reveals that the best docking scores were obtained for the walldegrading enzymes endoglucanase I. In addition, interactions of the best scored compound 3 with amino acid residuals of endoglucanase I generated the lowest total binding energy of all docked enzymes (−160.23 kcal/mol). Compound 3 exhibited the highest inhibition against M. phaseolina. In addition, compounds 7 and 4, which are among the first three compounds with the highest inhibition activity against M. phaseolina, have binding energies which are among the first four binding energies on endoglucanase I (−146.66 and −145.73 kcal/mol, respectively). Compound 3 also released the highest binding energy by forming complexes with the pectinase. The same results were obtained in our previous study of antifungal activity of coumarin derivatives, where the compounds that inhibited M. phaseolina most strongly had high binding energies on the plant cell-walldegrading enzymes, endoglucanase and pectinase [17]. Endoglucanase I is an enzyme involved in the degradation of cellulose by fungal species that catalyzes the hydrolysis of the α-1,4-glycosidic linkages of cellulose [37]. An ultrastructural microscopic observation showed that host cell penetration by M. phaseolina is very precise [40]. M. phaseolina excretes endoglucanase, which is rapidly transported through the xylem of an infected plant [41]. Endoglucanase activity from M. phaseolina was detected after 11 days of cultivation [24]. The amino acid sequence revealed a unique pathogen-specific endoglucanase (egl1) gene in addition to the endoglucanase gene (egl2) commonly found in M. phaseolina [42]. Endonogluconasa is also present in Fusarium oxysporum, especially in strain H57-1 [43]. Endo-1,4-glucanase (Cel12A), which belongs to the glycoside hydrolase family 12, was isolated from Trichoderma reesei Cel12A showed a strong hydrolysis against xyloglucan and (1 → 3,1 → 4)-β-glucan, the major polysaccharides of the cell wall [44]. It is important to note that a large amount of other hydrolytic enzymes for degrading cell wall polysaccharides, such as exocellobiohydrolases, β-glucosidases, as well as enzymes for lignin degradation (laccases, lignin peroxidases, galactose oxidases, chloroperoxidases, haloperoxidases, and heme peroxidases), which allows penetrating into the host tissue [45], was detected in M. phaseolina.    Table 6. Table 6. The energies of the main interactions between compound 3 and endoglucanase I (pdb ID: 2OVW).

Residue
Energy ( Figure 2b shows a 2D diagram of the main interaction with amino acid residues. The substrate-binding groove is located deeply at the bottom of the canyon, which is approximately 50 Å long and formed by the β-strands. The nucleophile, GLU-197, hydrogen bonds to an adjacent aspartate residue, ASP-199. Prior to a catalytic cycle, ASP-199 is in the protonated form, while GLU-197 is negatively charged [36]. Compound 3 is docked into the hydrophilic groove forming a conventional hydrogen bond TYR-171 through a hydrazide nitrogen atom (2.82 Å). One of the phenol rings creates strong π-σ interaction with TYR-177 (3.94 Å). Compound 3 forms the strongest van der Waals interaction with ASP-199, TYR-171, ASP-173, TYR-145, and ASP-199. Aromatic rings attached to the pyrazole ring, as well as benzopyran, generate several π-cation, π-σ interactions, π-π interactions, and π-alkyl interactions with surrounding amino acid residues, supporting the findings of the antifungal activity assay and literary evidence about the importance of aromatic rings in the structure for the enhanced activity [19].   Figure 2b shows a 2D diagram of the main interaction with amino acid residues. The substrate-binding groove is located deeply at the bottom of the canyon, which is approximately 50 Å long and formed by the β-strands. The nucleophile, GLU-197, hydrogen bonds to an adjacent aspartate residue, ASP-199. Prior to a catalytic cycle, ASP-199 is in the protonated form, while GLU-197 is negatively charged [36]. Compound 3 is docked into the hydrophilic groove forming a conventional hydrogen bond TYR-171 through a hydrazide nitrogen atom (2.82 Å). One of the phenol rings creates strong π-σ interaction with TYR-177 (3.94 Å). Compound 3 forms the strongest van der Waals interaction with ASP-199, TYR-171, ASP-173, TYR-145, and ASP-199. Aromatic rings attached to the pyrazole ring, as well as benzopyran, generate several π-cation, π-σ interactions, π-π interactions, and π-alkyl interactions with surrounding amino acid residues, supporting the findings of the antifungal activity assay and literary evidence about the importance of aromatic rings in the structure for the enhanced activity [19]. (green = conventional hydrogen bond; light green = van der Waals; brown = π-cation; dark purple = π-σ interactions; light purple = π-π interactions; pink = alkyl and π-alkyl interactions). light green = van der Waals; brown = π-cation; dark purple = π-σ interactions; light purple = π-π interactions; pink = alkyl and π-alkyl interactions).
Tested coumarin derivatives did not inhibit the growth of members of the beneficial bacterial soil population, Bacillus mycoides and Bradyrhizobium japonicum ( Table 2). This result is contrary to similar research by Dekić et al. [12] and Chen et al. [10], where the antibacterial efficacy of coumarin and coumarin derivatives with all tested bacteria was determined. This difference in antibacterial activity can be attributed to different coumarin structures.
Nematicidal activities of tested coumarin derivatives (expressed as % of inhibition) are presented in Table 2. Most of the compounds did not exhibit nematicidal activity against infective juveniles (IJs) of two beneficial nematode species, H. bacteriophora and S. feltiae. However, compound 9 was lethal for 46.25% IJs H. bacteriophora and 40.00% IJs S. feltiae after 48 h of incubation. Compound 2 was lethal for 8.75% IJs H. bacteriophora and 20.00% IJs S. feltiae. Compound 6 caused mortality for only 2.5% IJs S. feltiae. Previously, we identified several coumarin derivates thatcould potentially serve as new effective nematicides [17]. Specific phenolic compounds and coumarins have been reported by other authors as nematotoxic or affecting the development stages of the entomopathogenic nematodes [46][47][48][49][50]. The most active compounds potentially pose a threat to beneficial organisms in the rhizosphere, since entomopathogenic nematodes are con-sidered bioinsecticides. Further bioassays should include other nematode-trophic groups to identify selectivity of the most active coumarins. To better understand the mechanism of the nematicidal effect of analyzed compounds, AchE inhibition assay was performed. The AChE inhibition assay showed a low potency of the tested compounds as AChE inhibitors compared to standardized inhibitor donepezil (Table 2). However, compound 9 was an exception with 31.45% inhibition of the enzyme. Certain derivatives of Schiff bases showed pronounced activity as AChE inhibitors [51] as well as butyrylcholinesterase (BuChE) and monoamine oxidase (MAO) [52,53]. Obtained results are in correlation with the nematicidal activity assays where the same compound displayed the highest activity. As expected, obtained results confirmed a relation between inhibition of AChE and nematicidal activity, since this enzyme plays a critical role in terminating nerve impulses by hydrolyzing the neurotransmitter, acetylcholine (ACh) in the cholinergic nervous system of nematodes [54].
To understand the binding modes of tested coumarinyl Schiff bases and interactions with AChE at a molecular level, we performed blind-docking studies using AutoDock Vina software. The analysis of the lowest-energy poses of docked compounds showed their affinity to bind to the enzyme at two different sites. Compounds 9 and 8 bound to the active site with the binding energy of −7.8 kcal/mol and −7.6 kcal/mol, respectively, while compounds 7, 3, 1, 2, 6, 4, and 5 bound to the sites located away from (distal to) the enzyme active site with the binding energies of −9.6 kcal/mol, −9.4 kcal/mol, −9.1 kcal/mol, −9.1 kcal/mol, −9.0 kcal/mol, −8.9 kcal/mol, and −8.5 kcal/mol, respectively. Compound 9, with the lowest-energy pose at the enzyme active site, was used for further understanding of the binding mode with AChE. The binding orientation of compound 9 showed that it is placed near the bottom of the active site gorge (Figure 3a), interacting with amino acid residues of the anionic site (AS), acyl pocket (AP), and peripheral anionic site (PAS) of AChE. The AS is responsible for positioning acetylcholine in the enzyme cavity and as a site for possible inhibitor binding. The AP guarantees the specificity of the enzyme to the substrate by preventing the entry of larger molecules into the catalytic site, while PAS portion represents a region that is important for the binding of many inhibitors [55]. According to Figure 3b, PAS residues TYR-70, ASP-72, and TYR-334 formed a van der Waals interaction, an attractive charge interaction, and a π-π T-shaped interaction with compound 9, respectively. Amino acid residues PHE-330 and PHE-331 from the AChE AS interacted with compound 9 via van der Waals and carbon hydrogen interactions. Only one residue, PHE-288 from the AP, formed a van der Waals interaction with the phenyl ring of compound 9. Furthermore, interactions with compound 9 included hydrogen bonds with VAL-71, GLU-73, GLN-74, and GLN-272 which probably contributed to binding site stabilization. Additionally, van der Waals interactions were noted between compound 9 and residues ILE-275, ASP-276 and ILE-287. The involvement of the mentioned residues from the AS, AP, and PAS was confirmed to contribute stability to the complex between AChE and synthetic inhibitors such as coumarin-triazole-amino acid hybrids [55], tyrosol 1,2,3-triazole analogs [56], coumarin-3-carboxamide-N-morpholine hybrids [57], chromone derivatives [28], and chromenyl coumarate [58].
According to the rules for pesticide-likeness defined by Hao et al. [59], six molecular descriptors describe the distribution of pesticides. According to that rule, a pesticide-like compound should have: MW ≤ 435 Da; LOGP ≤ 6; HBA ≤ 6; HBD ≤ 2; RB ≤ 9; AB ≤ 17. Calculated descriptors for the nine coumarinyl Schiff bases are presented in the Table 7.
ing site stabilization. Additionally, van der Waals interactions were noted between compound 9 and residues ILE-275, ASP-276 and ILE-287. The involvement of the mentioned residues from the AS, AP, and PAS was confirmed to contribute stability to the complex between AChE and synthetic inhibitors such as coumarin-triazole-amino acid hybrids [55], tyrosol 1,2,3-triazole analogs [56], coumarin-3-carboxamide-N-morpholine hybrids [57], chromone derivatives [28], and chromenyl coumarate [58]. According to the rules for pesticide-likeness defined by Hao et al. [59], six molecular descriptors describe the distribution of pesticides. According to that rule, a pesticide-like compound should have: MW ≤ 435 Da; LOGP ≤ 6; HBA ≤ 6; HBD ≤ 2; RB ≤ 9; AB ≤ 17. Calculated descriptors for the nine coumarinyl Schiff bases are presented in the Table 7. None of the compounds completely satisfied pesticide-likeness criteria. Only compounds 8 and 9 had satisfactory molecular weight, but their number of hydrogen-bond acceptors was higher for one than the recommended six. The number of aromatic bonds of these compounds was also acceptable. Other compounds have high molecular weight and number of aromatics bond. The pesticide-likeness study revealed that a decrease of MW is associated with a toxic reduction of the pesticide. Aromatic double bonds in the  None of the compounds completely satisfied pesticide-likeness criteria. Only compounds 8 and 9 had satisfactory molecular weight, but their number of hydrogen-bond acceptors was higher for one than the recommended six. The number of aromatic bonds of these compounds was also acceptable. Other compounds have high molecular weight and number of aromatics bond. The pesticide-likeness study revealed that a decrease of MW is associated with a toxic reduction of the pesticide. Aromatic double bonds in the structure are usually associated with the photostability of pesticides, one of the most important properties of pesticide molecules [59,60]. From 1970 until the year 2000, the toxicity of pesticides was reduced by introducing the leading compounds with increased MW, which indirectly induced an increase in values for the other constitutive properties.
The toxicity of compounds estimated by program T.E.S.T. is presented in Table 8. According to the acute systemic toxicity classification based on oral LD 50 values for rats and recommended by the Organization for Economic Co-operation and Development (OECD), all compounds have estimated LD 50 values into the range of 300-2000 mg/kg, and characterized as "harmful if swallowed", except compound 5, which is the least harmful (LD 50 = 2206.90 mg/kg) and characterized as "may be harmful if swallowed" [61]. Only two, the most lipophilic compounds, 5 and 7 (Table 7), were evaluated as no mutagenic. Highest aquatic toxicity evaluated against the ciliate model organism Tetrahymena pyriformis was estimated for the least lipophilic compound 9 (Table 7). Compound 9 also proved to be the only lethal compound for nematode species H. bacteriophora and S. feltiae. Highly molecular weight and lipophilic compounds 1-5, and 7, proved to be highly toxic for the fathead minnow. Bioaccumulation is a process of absorption of compounds in an organism from the natural environment. All compounds had a logBAF of less than three, and their bioaccumulation was not significant [62].

Synthesis of Coumarinyl Schiff Bases
The synthesis of coumarin derivatives was performed in environmentally safe organic solvents [30].

Antifungal Assays
For the preparation of stock solutions of compound concentration 4 µmol/mL, a corresponding mass was dissolved in 2.5 mL of DMSO and 2.5 mL of distilled water. The volume of 1 mL of stock solution was added to the mixture to produce the final compound concentration (0.08 µmol/mL) and keep the amount of DMSO in the mixture at 1%. As a control, untreated potato dextrose agar (PDA) was used. The antifungal assay was performed on four cultures of phytopathogenic fungi from the fungal collection of the Faculty of Agrobiotehnical Sciences Osijek, (Fusarium oxysporum f. sp. lycopersici, Fusarium culmorum, Macrophomina phaseolina, and Sclerotinia sclerotiourum). The test was carried out according to the method of Siber et al. [63]. Petri dishes were kept in a growth chamber at 22 ± 1 • C, with a 12 h light/12 h dark regime. Each measurement consisted of four replicates. The radial growth of the fungal colonies was measured 48 h after inoculation. The in vitro inhibiting effects of the test compounds on the fungi were calculated by the antifungal index (% inhibition) [64].

Antibacterial Assays
Broth microdilution method was applied to test antibacterial activity and determine the minimum inhibitory concentration (MIC) of Bacillus mycoides (Gram-negative) and Bradyrhizobium japonicum (Gram-positive) (collection of the Faculty of Agrobiotehnical Sciences Osijek), which are included in the assay. Pure bacterial cultures were obtained by inoculation on nutrient (Liofilchem, Roseto degli Abruzzi, Italy) and Vincent agar. The stock solution of coumarinyl Schiff bases was prepared by dissolving 1.024 mg in 20 µL of DMSO with the addition of sterile water up to 200 µL. The solution was diluted from 512 to 1 µg/mL with sterile broth in 96-well plate. Inoculation was performed by adding 50 µL of pure bacterial cells with a density of 1.5 × 10 8 colony forming units (CFU/mL). Results were recorded after 48 h of incubation. The assay was set up in four repetitions.

Nematicidal Assays
For the preparation of 500 µg/mL stock solutions, 2 mg of each compound was dissolved in 20 µL of DMSO and 3980 µL of distilled water containing 0.1% Triton. Inhibition of nematode motility and mortality was tested for all compounds in maximum concentration 500 µg/ mL with four repetitions in a 24-well plate.
An aliquot of approximately 100 IJs entomopathogenic nematodes, Heterorhabditis bacteriophora same as Steinernema feltiae (collection of the Faculty of Agrobiotehnical Sciences Osijek) were placed in each well containing 250 µL of the working solution. Distilled water containing DMSO and Triton was used as a control. Well plates were incubated in the dark at 24 • C. The number of motile, dead, and live nematodes was recorded after 48 h of incubation. Nematodes were observed under a microscope (40×) and considered dead when they failed to respond to physical stimuli with a probe. The values were determined as percentage corrected mortality according to the Schneider-Orelli formula.
Acetylcholinesterase inhibition assay was conducted using the Ellman protocol [65]. Tested compounds were dissolved in DMSO and diluted with a phosphate buffer solution (pH 8.0), yielding the final concentration of DMSO in the assay less than 1%. The tested concentration of compounds in the final reaction mixture was 0.1 mM. Donepezil was used as a standard inhibitor. Enzyme stock solution was prepared by dissolving AChE in a phosphate buffer solution (pH 8.0) containing 0.1% (w/v) BSA. Reaction mixture was comprised of 10 µL of test compound solution, 130 µL of phosphate buffer, 20 µL of AChE solution (0.25 U/mL), 20 µL of DTNB solution, and 20 µL of ATChI solution. Prior to the addition of the substrate (ATChI), the mixture was preincubated for 15 min at 25 • C. The enzymatic activity measurements were conducted using a microplate reader (iMark, Bio-Rad, Hercules, CA, USA) in a span of 10 min at a wavelength of 412 nm. Inhibition activity was calculated using enzyme inhibition absorbance response using the formula: where A 1 is the absorbance of the tested compound, and A 0 is the absorbance of positive control.

Calculation of Toxicity
The toxicity of compounds was calculated entering their SMILES notation into the program T.E.S.T. v.4.1 using nearest neighbor method. The toxicity was estimated by taking an average of the three chemicals in the training set that were most similar to the test chemical [66]. A lethal dose for rats (oral rat LD 50 ) was expressed as the mg of the chemical per bodyweight of the rat (mg/kg bw). Aquatic toxicity of the compound was presented as the concentration of the test chemical in water in mol/L that caused 50% growth inhibition to Tetrahymena pyriformis after 48 h (48 h T. pyriformis IGC 50 ) [67] as well as by concentration in water, which killed half of the fathead minnows (Pimephales promelas) in 96 h (European Chemicals Agency, ECHA-11-R-004.2-EN, 2011) [68]. The Ames test estimates mutagenicity of a compound that induces revertant colony growth of Salmonella typhimurium. Its results represent the alert for the potential carcinogenicity and/or teratogenicity [69]. Bioaccumulation factor (BAF), is the ratio of the concentration of a chemical in the tissue of an aquatic organism (fish) to its concentration in water (in liters per kilogram of tissue), expressed as logarithmic value [70].
The molecular docking of compounds was performed using iGEMDOCK (BioXGEM, Taiwan). Applying the generic evolutionary method, each compound was docked into the binding site of a cocrystalized standard inhibitor using the following parameters: population size was 200, generations were 70, and the number of poses was 3. Program iGEMDOCK generates protein-compound interaction profiles based on electrostatic (E), hydrogenbonding (H), and van der Waals (V) interactions. Compounds were ranked by combining the pharmacological interactions and energy-based scoring function: Total Energy (kcal/mol) = vdW + Hbond + Elec. Blind docking was carried out to investigate the modes of AChE inhibition for tested compounds using AutoDock Vina 1.1.2 software [73]. For this purpose, the crystal structure of AChE, PDB ID: 2C5G [74] was extracted from the PDB, while MGL Tools 1.5.6 [75] was employed to prepare structures for molecular docking. The blind docking site for the ligands on AChE was defined by a grid box with the dimensions 81 × 82 × 84 Å3, and center set at x = 4.98, y = 64.57 and z = 55.94. Docking simulation was performed with the standard 0.375Å resolution and 20 conformations were generated. Receptor-ligand interactions were visualized with BIOVIA Discovery Studio Visualizer 4.5 (Dassault Systèmes, San Diego, CA, USA).

Statistical Analysis
Principal component analysis was performed by Statistica 14 (TIBCO Software Inc. 2020, Palo Alto, CA, USA).

Conclusions
Coumarinyl Schiff bases proved to be promising candidates in the inhibition of Macrophomina phaseolina. Pyrazole-based Schiff bases with the two aromatic rings with the weaker electron-withdrawing groups enhanced antifungal activity. The most effective antifungal compounds, 3 and 4, did not show negative effects against beneficial soil bacteria and nematodes. Their moderate lipophilicity and small molar weight do not pose a hazard for bioaccumulation in animals and humans. The possible inhibitory mechanism of the tested compounds against M. phaseolina is related to the inhibition of enzyme en-doglucanase, which participates in hydrolysis of cellulose. The most active compound 3, released the highest energy binding in the substrate-binding groove of endonuclease. In addition, the relation between nematicidal activity and inhibition of acetylcholinesterase was demonstrated. Compound nine was identified as nematotoxic and an AChE inhibitor in nematodes. This compound, as the lowest lipophilic, had the highest toxicity for the aquatic model organism. However, further toxicological studies are needed to support the development of new effective nematicides or anthelminthic drugs safe for mammals. Future bioassays should involve diverse species of target organisms to identify selectivity of the most active compounds.
Based on the knowledge gained by this study, the further development of coumarinyl Schiff bases as future plant protection products will continue. The future compounds should be more effective against pathogenic fungus and other pathogenic microorganisms without negative effects on beneficial organisms, environment pollution and human health.