Design, Synthesis and Evaluation of Novel Trichloromethyl Dichlorophenyl Triazole Derivatives as Potential Safener

The dominance of safener can unite with herbicides acquiring the efficient protection of crop and qualifying control of weeds in agricultural fields. In order to solve the crop toxicity problem and exploit the novel potential safener for fenoxaprop-P-ethyl herbicide, a series of trichloromethyl dichlorobenzene triazole derivatives were designed and synthesized by the principle of active subunit combination. A total of 21 novel substituted trichloromethyl dichlorobenzene triazole compounds were synthesized by substituted aminophenol and amino alcohol derivatives as the starting materials, using cyclization and acylation. All the compounds were unambiguously characterized by IR, 1H-NMR, 13C-NMR, and HRMS. A greenhouse bioassay indicated that most of the title compounds could protect wheat from injury caused by fenoxaprop-P-ethyl at varying degrees, in which compound 5o exhibited excellent safener activity at a concentration of 10 μmol/L and was superior to the commercialized compound fenchlorazole. A structure–activity relationship for the novel compounds was determined, which demonstrated that those compounds containing benzoxazine groups showed better activity than that of oxazole-substituted compounds. Introducing a benzoxazine fragment and electron-donating group to specific positions could improve or maintain the safener activity for wheat against attack by the herbicide fenoxaprop-P-ethyl. A molecular docking model suggested that a potential mechanism between 5o and fenoxaprop-P-ethyl is associated with the detoxication of the herbicide. Results from the present work revealed that compound 5o exhibited good crop safener activities toward wheat and could be a promising candidate structure for further research on wheat protection.


Introduction
Herbicides are a major tool for controlling weeds to protect crop yields. From the mid-1980s, more than 140 new herbicide active ingredients were commercialized and played an indispensable role in the rapid development of agriculture and rapid economic growth [1]. However, many were used for years and then eliminated from the market as unfavorable toxicology became apparent or safer and more effective replacements were found [2]. New challenges, such as herbicide toxicity, weed resistance, low crop selectivity, and high cost of discovery, production, and registration of herbicides, have been raised over the years [3]. Therefore, the design of more potent, selective, environmentally friendly and cost-effective herbicides is one of the main the focus on pesticide development [4]. In recent years, herbicides have entered the stage of efficient development and constantly updated varieties [5]. gramineous weeds in wheat fields and as an aryloxyphenoxypropionate postemergence herbicide inhibiting fatty acid synthesis in grasses through inhibition of acetyl-CoA carboxylase (ACCase, EC. 6.4.1.2) [6]. ACCase, the main target of FE, exerts a particularly important role in catalyzing fatty acid biosynthesis. Once initiated, the reaction commits malonyl-CoA to produce fatty acids. In plants, this reaction occurs mainly in plastids and, when inhibited by FE, eventually leads to plant leaf chlorosis death [7]. However, extensive application with the use of herbicides has increased phytotoxicity to the crop plants [8]. The pressure of multiple herbicides can lead to major practical problems, particularly because herbicide toxicity results in dramatic crop yield losses. Thus, it would be urgently necessary to investigate and develop potent safeners for obtaining crop safety and yield.
Herbicide safeners are a group of synthetic compounds that could protect selected plants against herbicide injury without decreasing herbicidal activity to target weeds. Alternatively, herbicide safeners may effectively improve herbicide selectivity [9][10][11]. Nearly 20 safeners have subsequently been developed and commercialized by famous agrochemical companies [12]. Fenchlorazole, one of the triazole herbicide safeners developed by Hoechst AG Company, has received much attention on account of its promising protection of wheat against FE [13]. It may enhance the metabolism of FE, not only by stimulating the activity of the key enzyme P450-dependent monooxygenase but also, by activating isozymes in the detoxification of FE [14]. The presence of the safener also enhances the activity of antiretroviral, glutathione-S-transferase (GST) and promotes the metabolism of herbicides [15]. Benoxacor, a chloroacetamide herbicide safener, can also protect crops against herbicides [16]. It protects maize from chloroacetanilide herbicides (e.g., alachlor) by increasing the expression of GST, which efficiently catalyzes the conjugation of herbicides with glutathione (GSH) [17].
The substructure combination method is very important in the search for bioactive materials, and it has been continuously used in the development of novel pesticides [18][19][20]. Recently, many successful cases have been reported. For example, isoxadifen-ethyl is developed as a herbicide safener with strong protection of rice by combining with the active structure of 5-phenyl-4,5dihydroisoxazole-3-ethyl ester and the safener diphenyl acid (Scheme 1) [21]. Diclofop-methyl, developed by Farbwerke Hoechst Company, is designed by pharmacophore combination of phenoxycarboxylic acid and diphenyl ether herbicides [22]. As we noted above and continuous of our research on the design of nitrogen-containing heterocyclic safeners [23][24][25][26], a series of the new safener skeleton structure compounds are designed based on active substructure combinations retaining the trichloromethyl dichlorobenzene triazole of fenchlorazole and the nitrogen-containing heterocyclic ring of benoxacor or R-28725 as the parent skeleton structure (Scheme 2). Herein, we report the synthesis of a series of novel substituted trichloromethyl dichlorobenzene triazole compounds via cyclization and acylation without any expensive reagent or catalyst. They are assessed as safeners to protect wheat against the herbicide fenoxaprop-P-ethyl. Furthermore, molecular docking analyses are extensively performed on the target compound to identify the possible detoxification mechanism for their safener potency. Biomolecules 2019, 9,

Chemicals
The melting points of target compounds were measured using a Beijing Taike melting point apparatus (X-4) (Beijing Tech Instrument Co. Ltd., Beijing, China) and were uncorrected. Infrared (IR) spectra were recorded on a Bruker ALPHA-T spectrometer (BRUKER Inc., Beijing, China) using KBr disks. 1 H and 13 C-NMR nuclear magnetic resonance (NMR) spectra were obtained using a Bruker AV-400 spectrometer (BRUKER Inc., Beijing, China) in CDCl3 solution with tetramethylsilane as the internal standard. The high-resolution mass spectrometry (HRMS) was recorded on micrOTOF-Q II 10410 spectrometer of Bruker (BRUKER Inc., Beijing, China) and FT-ICR MS spectrometer of Bruker (BRUKER Inc., Beijing, China). The single-crystal diffraction experiment was carried out on a Rigaku R-AXIS RAPID area-detector diffractometer (Rigaku Company, Tokyo, Japan). All solvents and reagents were of analytical reagent grade. Column chromatography purification was carried out on silica gel. (3) 1-(2,4-Dichlorobenzyl)-5-trichloromethyl-1,2,4-triazole-3-carboxylate 1 (10 mmol) was dissolved in EtOH (20 mL) and NaOH aq. (1 mol/L, 10 mL), and the mixture was stirred at 25 °C for 5 h. Next, distilled water (80 mL) was added into the mixture. The reaction mixture was then acidified by concentrated HCl solution to pH = 1 to induce precipitation. The precipitate was collected through vacuum filtration, washed with water and dried to afford product 2 as white solid without purification. A mixture of compound 2 (37 mmol), sulfoxide chloride (137 mmol, 10 mL) and CH2Cl2 (20 mL) was stirred for 1 h at 40 °C. Then, five drops of N,N-dimethylformamide (DMF) was added. After reaction completed, the resulting mixture was concentrated in vacuum to give compound 3 as pale yellow solid (yield 87%), which was used in the next step without further purification.

General Procedure for the Synthesis of Intermediates Compounds (4)
The preparation of compounds 4a-4u were conducted as the method reported by references [27][28][29][30][31]. Substituted o-aminophenol (27 mmol) was added to a stirred solution of DMSO (40 mL) and K2CO3 (7.5 g, 54 mmol) of substituted dibromoethane (27 mmol) or (CH2Cl2, 41 mmol for 4p-4q) at reflux. The mixture was heated under reflux for 6-8 h until the reaction was completed (TLC monitored). The resulting mixture was then filtered, and the filtrate was concentrated under reduced pressure to obtain compounds 4a-4q.
Ketone (26 mmol) and substituted o-amino alcohol (26 mmol) were dissolved in toluene (30 mL). After that, the reaction mixture was stirred for another 8 min with microwave irradiation (85 °C, 800 W) to obtain compounds 4r-4u. All the compounds 4a-4u were directly subjected to the next reaction step.

Chemicals
The melting points of target compounds were measured using a Beijing Taike melting point apparatus (X-4) (Beijing Tech Instrument Co. Ltd., Beijing, China) and were uncorrected. Infrared (IR) spectra were recorded on a Bruker ALPHA-T spectrometer (BRUKER Inc., Beijing, China) using KBr disks. 1 H and 13 C-NMR nuclear magnetic resonance (NMR) spectra were obtained using a Bruker AV-400 spectrometer (BRUKER Inc., Beijing, China) in CDCl 3 solution with tetramethylsilane as the internal standard. The high-resolution mass spectrometry (HRMS) was recorded on micrOTOF-Q II 10410 spectrometer of Bruker (BRUKER Inc., Beijing, China) and FT-ICR MS spectrometer of Bruker (BRUKER Inc., Beijing, China). The single-crystal diffraction experiment was carried out on a Rigaku R-AXIS RAPID area-detector diffractometer (Rigaku Company, Tokyo, Japan). All solvents and reagents were of analytical reagent grade. Column chromatography purification was carried out on silica gel. (3) 1-(2,4-Dichlorobenzyl)-5-trichloromethyl-1,2,4-triazole-3-carboxylate 1 (10 mmol) was dissolved in EtOH (20 mL) and NaOH aq. (1 mol/L, 10 mL), and the mixture was stirred at 25 • C for 5 h. Next, distilled water (80 mL) was added into the mixture. The reaction mixture was then acidified by concentrated HCl solution to pH = 1 to induce precipitation. The precipitate was collected through vacuum filtration, washed with water and dried to afford product 2 as white solid without purification. A mixture of compound 2 (37 mmol), sulfoxide chloride (137 mmol, 10 mL) and CH 2 Cl 2 (20 mL) was stirred for 1 h at 40 • C. Then, five drops of N,N-dimethylformamide (DMF) was added. After reaction completed, the resulting mixture was concentrated in vacuum to give compound 3 as pale yellow solid (yield 87%), which was used in the next step without further purification.

General Procedure for the Synthesis of Intermediates Compounds (4)
The preparation of compounds 4a-4u were conducted as the method reported by references [27][28][29][30][31]. Substituted o-aminophenol (27 mmol) was added to a stirred solution of DMSO (40 mL) and K 2 CO 3 (7.5 g, 54 mmol) of substituted dibromoethane (27 mmol) or (CH 2 Cl 2 , 41 mmol for 4p-4q) at reflux. The mixture was heated under reflux for 6-8 h until the reaction was completed (TLC monitored). The resulting mixture was then filtered, and the filtrate was concentrated under reduced pressure to obtain compounds 4a-4q.
Ketone (26 mmol) and substituted o-amino alcohol (26 mmol) were dissolved in toluene (30 mL). After that, the reaction mixture was stirred for another 8 min with microwave irradiation (85 • C, 800 W) to obtain compounds 4r-4u. All the compounds 4a-4u were directly subjected to the next reaction step.

General Procedure for the Synthesis of Title Compounds (5)
Compound 3 (15 mmol) and appropriate compound 4 (12.5 mmol) were added to DMAC (30 mL) in succession and the mixture was stirred at room temperature for 2 h. NaOH aq. (2 mol/L, 3 mL) Biomolecules 2019, 9, 438 4 of 14 was added to the mixture and the reaction was monitored by TLC. The solvent was evaporated, and the residue was dissolved with ethyl acetate, washed with water, dried and concentrated. The crude product was purified on silica gel with petroleum ether/ethyl acetate (6:1, v/v) to obtain the target compounds.

X-ray Diffraction
A suitable single crystal of compound 5j (CCDC 1909489) was obtained by dissolving the compound in a solution of ethyl acetate, and the solvent was placed at room temperature to slowly evaporate the solvent in air for several days. Bright white crystals of 5j (0.13 mm × 0.12 mm × 0.10 mm) were selected and mounted on a Rigaku R-AXIS RAPID area detector diffractometer equipped with graphite-monochromated Mo-Kα radiation (λ = 0.71073Å); θmax = 27.58. The crystal was kept at 293 K during data collection. Direct methods were used for structure solving by the SHELXS-97 program.

Biological Assay
The concentration of the safener and compounds applied in the bioassay was determined after a preliminary screening. Wheat seeds (College of Agriculture, Northeast Agricultural University, Harbin, China) were moistened with warm water for 30 min, and then soaked with 0.6% carbendazim for about 30 min. Afterwards, the seeds were washed with distilled water and soaked by title compounds (10 µmol/L) for 12 h, and then germinated for 24 h at 26.5 ± 1 • C. The commercial safener fenchlorazole was chosen as a positive control. The seeds were planted 1.0 cm deep in paper cups, and then sprayed by FE until the two-leaf stage was reached. After 7 days, the effects of the title compounds on the detoxification of FE in soil were evaluated by testing the growth index [32,33].

Computational Methods
The steric structure of compound 5o and commercial safener fenchlorazole were generated with the sketch module of SYBYL-X 2.0 (Tripos Inc., Saint Louis, MO, USA). Subsequently, the Gasteiger-Huckel charges were calculated, and the tested molecules were optimized. The crystallographic structure of ACCase was obtained from the Protein Data Bank (PDB: 1UYR). Subsequently, docking calculations were performed using the CDOCKER method in Accelrys Discovery Studio 3.0 (BIOVIA Inc., San Diego, CA, USA) [34,35]. Before the docking, the CHARMM force field was added to the protein structure, and other co-crystallized small molecules and water were isolated. After that, the active site was defined, with a subset region of 13.0 Å by the center of the known ligand. In the docking process, the top 10 conformations were reserved for every ligand based on -CDOCKER_ENERGY and the remaining parameters were set to the default values.

Synthesis
The preparation of compounds 3 involved the hydrolysis reaction and acylating chlorination by fenchlorazole-ethyl (Scheme 3). The intermediates 4 were synthesized with substituted aminophenol or amino alcohol derivatives as the starting materials (Scheme 4). The synthetic procedures for the title compounds are depicted in Scheme 5. Compounds 5 were obtained by the compound 3 with compounds 4 (yields in range of 50-82%) ( Table 1).

Synthesis
The preparation of compounds 3 involved the hydrolysis reaction and acylating chlorination by fenchlorazole-ethyl (Scheme 3). The intermediates 4 were synthesized with substituted aminophenol or amino alcohol derivatives as the starting materials (Scheme 4). The synthetic procedures for the title compounds are depicted in Scheme 5. Compounds 5 were obtained by the compound 3 with compounds 4 (yields in range of 50-82%) ( Table 1).  All the title products were characterized by IR, 1 H-NMR, 13 C-NMR, and HRMS spectroscopies. Taking compound 5j as an example. IR analysis showed the existence of a carbonyl group at approximately 1665 cm −1 and 1253 cm −1 for compound 5j. In the 1 H-NMR spectrum analysis, the signals in the region of δ 7.62-6.90 ppm related to benzene ring-bound protons. The signals at δ 1.57-1.23 ppm were characteristic of two methyl groups of the oxazine ring. In the 13 C-NMR spectra, the signals observed at δ 158.25 ppm were attributed to the carbon of the carbonyl group. The treated HRMS information also corroborated the accuracy of the proposed structures. Meanwhile, the molecular structure of 5j was additionally determined by X-ray analysis, as seen from Figure 1. All the title products were characterized by IR, 1 H-NMR, 13 C-NMR, and HRMS spectroscopies. Taking compound 5j as an example. IR analysis showed the existence of a carbonyl group at approximately 1665 cm −1 and 1253 cm −1 for compound 5j. In the 1 H-NMR spectrum analysis, the signals in the region of δ 7.62-6.90 ppm related to benzene ring-bound protons. The signals at δ 1.57-1.23 ppm were characteristic of two methyl groups of the oxazine ring. In the 13 C-NMR spectra, the signals observed at δ 158.25 ppm were attributed to the carbon of the carbonyl group. The treated HRMS information also corroborated the accuracy of the proposed structures. Meanwhile, the molecular structure of 5j was additionally determined by X-ray analysis, as seen from Figure 1.

The Structure-Activity Relationships
All compounds 5 (10 μmol/L) were evaluated for their protection of wheat in vivo against FE ( Table 2). The growth index was measured after 7 days of treatment including recovery of root length, recovery of plant weight and chlorophyll content.
The herbicide FE was absorbed by the plant and converted into the corresponding acid, which reached the growth point of the plant meristem and the root growth point, affecting plant growth [36,37]. FE provoked a significant chlorosis and inhibition of growth for maize, but recovered when the safener fenchlorazole was included [38]. The chlorophyll content of plant leaves directly affects photosynthesis associated with plant growth [39]. A positive correlation between growth level and recovery of root length, recovery of plant weight and chlorophyll content has been observed in previous research [40]. As seen from Table 2, it was observed that compounds 5 showed different recoveries of wheat growth. Compound 5o possessed good protection against FE, as indicated by the 70.76% recovery of root length, 143.51% recovery of plant weight, and 99.77% recovery of chlorophyll, respectively, which was superior to the commercial safener fenchlorazole. Compounds 5a-5o showed a more significant improvement of growth index than compounds 5p-5u due to the introduction of substituted benzoxazine. For the preliminary assessment of the substituents on the benzoxazine derivatives, the electron-donating group in the benzoxazine ring was crucial for good safener activity. It was noticeable that compounds 5c and 5i exhibited superior detoxification activity compared to compound 5e due to the replacement of the Br atom with the methyl or ethyl of the substitution at 4-position of benzene ring. Compound 5c exerted a strong protective effect against FE with recovery rates of root length 71.2% and recovery of plant weight 109.87%, respectively. The result indicated that compounds with a benzoxazine fragment and electron-donating groups were greatly associated with the protective effect of wheat.

The Structure-Activity Relationships
All compounds 5 (10 µmol/L) were evaluated for their protection of wheat in vivo against FE ( Table 2). The growth index was measured after 7 days of treatment including recovery of root length, recovery of plant weight and chlorophyll content. The herbicide FE was absorbed by the plant and converted into the corresponding acid, which reached the growth point of the plant meristem and the root growth point, affecting plant growth [36,37]. FE provoked a significant chlorosis and inhibition of growth for maize, but recovered when the safener fenchlorazole was included [38]. The chlorophyll content of plant leaves directly affects photosynthesis associated with plant growth [39]. A positive correlation between growth level and recovery of root length, recovery of plant weight and chlorophyll content has been observed in previous research [40]. As seen from Table 2, it was observed that compounds 5 showed different recoveries of wheat growth. Compound 5o possessed good protection against FE, as indicated by the 70.76% recovery of root length, 143.51% recovery of plant weight, and 99.77% recovery of chlorophyll, respectively, which was superior to the commercial safener fenchlorazole. Compounds 5a-5o showed a more significant improvement of growth index than compounds 5p-5u due to the introduction of substituted benzoxazine. For the preliminary assessment of the substituents on the benzoxazine derivatives, the electron-donating group in the benzoxazine ring was crucial for good safener activity. It was noticeable that compounds 5c and 5i exhibited superior detoxification activity compared to compound 5e due to the replacement of the Br atom with the methyl or ethyl of the substitution at 4-position of benzene ring. Compound 5c exerted a strong protective effect against FE with recovery rates of root length 71.2% and recovery of plant weight 109.87%, respectively. The result indicated that compounds with a benzoxazine fragment and electron-donating groups were greatly associated with the protective effect of wheat.
An improvement safener activity of 5j-5m against FE was observed relative to the 4-position of the benzene. Additionally, compounds in which substitution at 4-position of benzene ring is by two methyl groups were always favorable to the safener potency than that of one methyl. Most interesting, comparing the protective effects of compounds 5j-5m showed that 4-position substituted with at least one methyl group led to a remarkable improvement in safener activity.
With the above encouraging results, it can be concluded that the activity of compounds with the same parent skeleton was even more related to the various substituents ( Figure 2). Moreover, it was found that substitution at the 4-position of benzene was essential for maintaining safener activity. Particularly, the protection of compound 5o was superior to that of fenchlorazole, which is a commercial safener and effectively reduced injury from FE herbicides.

Molecular Docking Studies
It is reasonable to expect that the physical properties of a molecule on its ability to express activity are important in the selection of potential molecular candidates. Then, the fenoxaprop-Pethyl, fenchlorazole and compound 5o were compared and analyzed. It was observed that the log p and HBAs of compound 5o were consistent with that of the safener fenchlorazole (Table 3). There was little difference in ARs and RBs between fenchlorazole and compound 5o. Molecular docking experiment showed that the fenchlorazole and compound 5o exited multiple interacts with active site such as hydrogen bond and Van der Waals forces etc. However, fenchlorazole and compound 5o were mainly dominated by hydrogen bonds, therefore the slight difference between ARs and RBs has negligible effect to the safety activity.

Molecular Docking Studies
It is reasonable to expect that the physical properties of a molecule on its ability to express activity are important in the selection of potential molecular candidates. Then, the fenoxaprop-P-ethyl, fenchlorazole and compound 5o were compared and analyzed. It was observed that the log p and HBAs of compound 5o were consistent with that of the safener fenchlorazole (Table 3). There was little difference in ARs and RBs between fenchlorazole and compound 5o. Molecular docking experiment showed that the fenchlorazole and compound 5o exited multiple interacts with active site such as hydrogen bond and Van der Waals forces etc. However, fenchlorazole and compound 5o were mainly dominated by hydrogen bonds, therefore the slight difference between ARs and RBs has negligible effect to the safety activity.

Molecular Docking Studies
It is reasonable to expect that the physical properties of a molecule on its ability to express activity are important in the selection of potential molecular candidates. Then, the fenoxaprop-Pethyl, fenchlorazole and compound 5o were compared and analyzed. It was observed that the log p and HBAs of compound 5o were consistent with that of the safener fenchlorazole (Table 3). There was little difference in ARs and RBs between fenchlorazole and compound 5o. Molecular docking experiment showed that the fenchlorazole and compound 5o exited multiple interacts with active site such as hydrogen bond and Van der Waals forces etc. However, fenchlorazole and compound 5o were mainly dominated by hydrogen bonds, therefore the slight difference between ARs and RBs has negligible effect to the safety activity. To further explore the possible and preliminary effect of the relation mode of compound 5o involved in FE herbicides, a molecular docking experiment was conducted ( Figure 3). In the docking modeling research, the whole molecular skeleton of fenoxaprop acid, in which the ethyl ester is a procide of fenoxaprop acid, could be effectively embedded in the target active site by stronger Hbond interactions with ALA-1627 and ILE-1735, and two π-π interactions with TYR-1738, which gave rise to herbicide activity [41]. Interestingly, both compound 5o and fenchlorazole could also enter the active channel. In contrast, both compound 5o and fenchlorazole were located below the active binding site, compared with that of fenoxaprop acid, which just partially obstructs the entrance and interacts with very weak forces such as Van der Waals forces etc. The small substrate had a greater

Molecular Docking Studies
It is reasonable to expect that the physical properties of a molecule on its ability to express activity are important in the selection of potential molecular candidates. Then, the fenoxaprop-Pethyl, fenchlorazole and compound 5o were compared and analyzed. It was observed that the log p and HBAs of compound 5o were consistent with that of the safener fenchlorazole (Table 3). There was little difference in ARs and RBs between fenchlorazole and compound 5o. Molecular docking experiment showed that the fenchlorazole and compound 5o exited multiple interacts with active site such as hydrogen bond and Van der Waals forces etc. However, fenchlorazole and compound 5o were mainly dominated by hydrogen bonds, therefore the slight difference between ARs and RBs has negligible effect to the safety activity. To further explore the possible and preliminary effect of the relation mode of compound 5o involved in FE herbicides, a molecular docking experiment was conducted ( Figure 3). In the docking modeling research, the whole molecular skeleton of fenoxaprop acid, in which the ethyl ester is a procide of fenoxaprop acid, could be effectively embedded in the target active site by stronger Hbond interactions with ALA-1627 and ILE-1735, and two π-π interactions with TYR-1738, which gave rise to herbicide activity [41]. Interestingly, both compound 5o and fenchlorazole could also enter the active channel. In contrast, both compound 5o and fenchlorazole were located below the active binding site, compared with that of fenoxaprop acid, which just partially obstructs the entrance and interacts with very weak forces such as Van der Waals forces etc. The small substrate had a greater 5o 6.02 3 431 3 4

Molecular Docking Studies
It is reasonable to expect that the physical properties of a molecule on its ability to express activity are important in the selection of potential molecular candidates. Then, the fenoxaprop-Pethyl, fenchlorazole and compound 5o were compared and analyzed. It was observed that the log p and HBAs of compound 5o were consistent with that of the safener fenchlorazole (Table 3). There was little difference in ARs and RBs between fenchlorazole and compound 5o. Molecular docking experiment showed that the fenchlorazole and compound 5o exited multiple interacts with active site such as hydrogen bond and Van der Waals forces etc. However, fenchlorazole and compound 5o were mainly dominated by hydrogen bonds, therefore the slight difference between ARs and RBs has negligible effect to the safety activity. To further explore the possible and preliminary effect of the relation mode of compound 5o involved in FE herbicides, a molecular docking experiment was conducted ( Figure 3). In the docking modeling research, the whole molecular skeleton of fenoxaprop acid, in which the ethyl ester is a procide of fenoxaprop acid, could be effectively embedded in the target active site by stronger Hbond interactions with ALA-1627 and ILE-1735, and two π-π interactions with TYR-1738, which gave rise to herbicide activity [41]. Interestingly, both compound 5o and fenchlorazole could also enter the active channel. In contrast, both compound 5o and fenchlorazole were located below the active binding site, compared with that of fenoxaprop acid, which just partially obstructs the entrance and interacts with very weak forces such as Van der Waals forces etc. The small substrate had a greater a The log p was were predicted by ChemBioOffice 2015. b The numbers of aromatic rings (ARs), the surface areas (SA), number of rotatable bonds (RBs) and hydrogen bond acceptors (HBAs) were predicted by Discovery Studio 2.5. c The electronegativity was predicted by Sybyl-X 2.0 (Tripos Inc., St. Louis, MO, USA).
To further explore the possible and preliminary effect of the relation mode of compound 5o involved in FE herbicides, a molecular docking experiment was conducted ( Figure 3). In the docking modeling research, the whole molecular skeleton of fenoxaprop acid, in which the ethyl ester is a procide of fenoxaprop acid, could be effectively embedded in the target active site by stronger H-bond interactions with ALA-1627 and ILE-1735, and two π-π interactions with TYR-1738, which gave rise to herbicide activity [41]. Interestingly, both compound 5o and fenchlorazole could also enter the active channel. In contrast, both compound 5o and fenchlorazole were located below the active binding site, compared with that of fenoxaprop acid, which just partially obstructs the entrance and interacts with very weak forces such as Van der Waals forces etc. The small substrate had a greater chance to propel itself into the channel leading to the active site and catalyze the active site, preventing the binding of fenoxaprop acid with ACCase [42].  This result showed that the effect of competing with the herbicide at the active site might be ruled out and demonstrated the validity of the previously active fragment we designed, which was closely related to the results of bioassay test. When compound 5o as safener was used as safener synchronously or before FE, antagonism potentially existed with FE in wheat, which protected the wheat from injury by FE herbicide. This result was similar to that of previous research [43,44], and might be the detoxification mechanism of the novel target compounds. This result showed that the effect of competing with the herbicide at the active site might be ruled out and demonstrated the validity of the previously active fragment we designed, which was closely related to the results of bioassay test. When compound 5o as safener was used as safener synchronously or before FE, antagonism potentially existed with FE in wheat, which protected the wheat from injury by FE herbicide. This result was similar to that of previous research [43,44], and might be the detoxification mechanism of the novel target compounds.

Conclusions
In summary, a series of new substituted trichloromethyl dichlorobenzene triazole compounds were rationally designed by the substructure combination method. All the title compounds were confirmed by 1 H-NMR, 13 C-NMR, and HRMS. The bioassay showed that most of them exhibited varying amounts of safening activity to wheat from FE. Notably, compound 5o exhibited better recovery rate of root length, plant height and chlorophyll with 70.76%, 143.51%, and 99.77%, respectively, which was even superior to the commercial safener fenchlorazole. The molecular docking experiment illustrated that the protective effect of compound 5o might be attributed to antagonism with FE at the active site. The present work indicates that compound 5o could be used as a candidate for further optimization of a potential safener.