New Substituted Benzoylthiourea Derivatives: From Design to Antimicrobial Applications

The increasing threat of antimicrobial resistance to all currently available therapeutic agents has urged the development of novel antimicrobials. In this context, a series of new benzoylthiourea derivatives substituted with one or more fluorine atoms and with the trifluoromethyl group have been tested, synthesized, and characterized by IR, NMR, CHNS and crystal X-ray diffraction. The molecular docking has provided information regarding the binding affinity and the orientation of the new compounds to Escherichia coli DNA gyrase B. The docking score predicted the antimicrobial activity of the studied compounds, especially against E. coli, which was further demonstrated experimentally against planktonic and biofilm embedded bacterial and fungal cells. The compounds bearing one fluorine atom on the phenyl ring have shown the best antibacterial effect, while those with three fluorine atoms exhibited the most intensive antifungal activity. All tested compounds exhibited antibiofilm activity, correlated with the trifluoromethyl substituent, most favorable in para position.

In the last 15 years the fluoro-derivatives class has flourished. In 2018 FDA approved 17 new molecules, among which fostamatinib, baloxavir, marboxil, doravirine, ivosidenib, and apalutamide [28]. The isosteric modification, consisting of replacing hydrogen with fluorine, modifies the steric and electronic properties of the resulting molecules, affecting their metabolic transformations because of the relative stability given by the C-F energy bond (485.7 kJ/mol) [29]. The presence of the fluorine atom modulates the pharmacodynamic and pharmacokinetic properties [25]. Drug fluorination is commonly used to modulate pKa values, to enhance membrane permeation, and to increase liposolubility and, thus, their in vivo absorption and transport [30,31]. The thiourea skeleton favors other chemical modulations (e.g., incorporation of fluorine), these derivatives being thus promising candidates, both as active molecules and ligands. For instance, the fluoroquinolone derivatives having thioureido moiety and piperazin-1-yl groups have shown a significant improvement of the activity against Gram-positive (S. aureus, B. subtilis) and Gram-negative (E. coli, P. aeruginosa) bacterial strains [32]. The 2-(6'-fluorobenzothiazol-2'-ylamino)-4,6-(disubstituted thioureido)-1,3-pyrimidine derivatives and the 1-(isomeric fluorobenzoyl)-3-(isomeric fluorophenyl)thiourea compounds have shown inhibitory activity against the above-mentioned four species, correlated with the presence of nitro groups and with the presence and position of fluorine substituent on aroyl and aryl rings [33,34]. The last ones exhibited a better antifungal activity against Thiophanate and thiophanate methyl are thiourea-based fungicides largely used on different fruits and vegetable crops. The 1,3-disubstituted aliphatic and aromatic symmetrical and unsymmetrical thioureas proved to be effective against Pyricularia oryzae and Drechslera oryzae phytopathogens [26], while N-(o-fluorophenoxyacetyl)thiourea derivatives exhibited herbicidal activity against Amaranthus retroflexus L. [27].
In the last 15 years the fluoro-derivatives class has flourished. In 2018 FDA approved 17 new molecules, among which fostamatinib, baloxavir, marboxil, doravirine, ivosidenib, and apalutamide [28]. The isosteric modification, consisting of replacing hydrogen with fluorine, modifies the steric and electronic properties of the resulting molecules, affecting their metabolic transformations because of the relative stability given by the C-F energy bond (485.7 kJ/mol) [29]. The presence of the fluorine atom modulates the pharmacodynamic and pharmacokinetic properties [25]. Drug fluorination is commonly used to modulate pKa values, to enhance membrane permeation, and to increase liposolubility and, thus, their in vivo absorption and transport [30,31]. The thiourea skeleton favors other chemical modulations (e.g., incorporation of fluorine), these derivatives being thus promising candidates, both as active molecules and ligands. For instance, the fluoroquinolone derivatives having thioureido moiety and piperazin-1-yl groups have shown a significant improvement of the activity against Gram-positive (S. aureus, B. subtilis) and Gram-negative (E. coli, P. aeruginosa) bacterial strains [32]. The 2-(6'-fluorobenzothiazol-2'-ylamino)-4,6-(disubstituted thioureido)-1,3-pyrimidine derivatives and the 1-(isomeric fluorobenzoyl)-3-(isomeric fluorophenyl)thiourea compounds have shown inhibitory activity against the above-mentioned four species, correlated with the presence of nitro groups and with the presence and position of fluorine substituent on aroyl and aryl rings [33,34]. The last ones exhibited a better antifungal activity against Rhizopus oryzae, Aspergillus terreus, Fusarium oxysporum, A. niger, and A. fumigatus, probably due to the presence of fluorine, which increased their lipophilicity, penetration rate, and intracellular transport to the microbial target. An improved receptor binding could be achieved due to the higher polarizability assured by the C-F bond [34].
On the other hand, the trifluoromethyl group is an important substituent in many drugs or agrochemical compounds because it increases the lipophilicity and metabolic stability [29]. Therefore, a number of synthesis methods for such compounds have been developed [35] so that it can be expected that in the future a lot of new molecules with biological activity from the class of trifluoromethyl derivatives will be developed [36].
The common isotope of fluorine ( 19 F) is the favorite heteroatom for incorporation into drugs [26], while the isotope ( 18 F) is used for radiolabelling of some ligands used in HD PET imaging for Huntington disease diagnosis [38]. Recently, fluorodopa (FDOPATM), the fluorine-18 isotopologue of L-DOPA, was approved as a radiotracer in PET for the visualization of nerve cells in patients with symptoms of Parkinson disease [39].
Presently, the emergence of antimicrobial resistance in human, animal, and plant pathogens, correlated with the unavailability of newer drugs challenges the efficiency of current antimicrobial agents [40]. Therefore, comprehensive and multidisciplinary efforts across human medicine, veterinary, agricultural, and environmental sectors are needed to reduce this global threat by investigating the mechanisms of resistance, improving antimicrobial stewardship, and developing antimicrobial agents [41][42][43].
Our aim was to exploit the multiple advantages of 2-((4-ethylphenoxy)methyl)-N-(fluoro/trifluorosubstituted phenylcarbamothioyl)benzamides to obtain potent antibacterial and antifungal agents, potentially useful for both clinical and agricultural sectors.

Chemistry
The new compounds, 5a-g, were obtained with good yields using the synthesis route presented in Scheme 1. The reaction between compound 3 with ammonium thiocyanate in dry acetone resulted in in situ 2-(4-ethylphenoxymethyl)benzoyl isothiocyanate (4). The appropriate primary aromatic amines were directly added to the reaction mixture, in dry acetone, resulting in new derivatives 5a-g. CH  The obtained compounds were solid, white, or light-yellow crystals, exhibiting solubility in chloroform and acetone at room temperature, soluble in lower alcohols, benzene, toluene, and xylene after heating and insoluble in water.
The melting point, elemental analysis, infrared (IR), and nuclear magnetic resonance (NMR) spectral studies were used for the characterization and structure confirmation of the new compounds.
In the IR spectra, the stretching bands due to νN-H of the amide group can be found to the highest values of the wave numbers. These were sharp peaks with a medium intensity occurring in the region 3401 cm −1 . The thioamide group showed a less intense stretching band at 3155-3183 cm −1 and, with a high probability, the band located at 1374-1392 cm −1 can be attributed to the thioamide group. For the observed antisymmetric stretching vibrations, methyl and methylene groups gave saturated νC-H stretching bands at about 2959-2970 cm −1 and 2912-2933 cm −1 , respectively; These bands are typical for aromatic compounds containing saturated carbons. A medium or intense sharp stretching band, shown in the IR spectrum of these compounds in the region 1676-1695 cm −1 is due to the νC=O vibrations. The νC=O band was lower than that observed with the ordinary carbonyl absorption (1730 cm −1 ). This low frequency value can be explained as a result of conjugated resonance with the phenyl ring and the formation of intramolecular hydrogen bonding with N-H. Near to this peak lies the intense band of δN-H, with a maximum at 1506-1514 cm −1 , which overlaps the aromatic core vibrations. These compounds also show a typical alkyl-aryl ether at 1237-1256 Scheme 1. Synthesis of the new derivatives (5a-g).
The reaction between compound 3 with ammonium thiocyanate in dry acetone resulted in in situ 2-(4-ethylphenoxymethyl)benzoyl isothiocyanate (4). The appropriate primary aromatic amines were directly added to the reaction mixture, in dry acetone, resulting in new derivatives 5a-g.
The obtained compounds were solid, white, or light-yellow crystals, exhibiting solubility in chloroform and acetone at room temperature, soluble in lower alcohols, benzene, toluene, and xylene after heating and insoluble in water.
The melting point, elemental analysis, infrared (IR), and nuclear magnetic resonance (NMR) spectral studies were used for the characterization and structure confirmation of the new compounds.
In the IR spectra, the stretching bands due to νN-H of the amide group can be found to the highest values of the wave numbers. These were sharp peaks with a medium intensity occurring in the region 3401 cm −1 . The thioamide group showed a less intense stretching band at 3155-3183 cm −1 and, with a high probability, the band located at 1374-1392 cm −1 can be attributed to the thioamide group. For the observed antisymmetric stretching vibrations, methyl and methylene groups gave saturated νC-H stretching bands at about 2959-2970 cm −1 and 2912-2933 cm −1 , respectively; These bands are typical for aromatic compounds containing saturated carbons. A medium or intense sharp stretching band, shown in the IR spectrum of these compounds in the region 1676-1695 cm −1 is due Molecules 2020, 25, 1478 5 of 20 to the νC=O vibrations. The νC=O band was lower than that observed with the ordinary carbonyl absorption (1730 cm −1 ). This low frequency value can be explained as a result of conjugated resonance with the phenyl ring and the formation of intramolecular hydrogen bonding with N-H. Near to this peak lies the intense band of δN-H, with a maximum at 1506-1514 cm −1 , which overlaps the aromatic core vibrations. These compounds also show a typical alkyl-aryl ether at 1237-1256 cm −1 for the antisymmetric vibration, and 1013-1045 cm −1 for the symmetric one. The presence of fluorine in the new compounds is proven by the stretching bands situated at 1013-1079 cm −1 .
The structure of the new compounds was also confirmed by the NMR spectra. The new compounds were dissolved in DMSO-d6 (hexadeuterio-dimethyl sulfoxide) and the chemical shifts' values, expressed in parts per million (ppm), were referenced downfield to tetramethylsilane for 1 H-NMR and 13 C-NMR and upfield to trichloro-fluoro-methane for 19 F-NMR and the constants (J) values in Hertz. The chemical shifts for hydrogen and carbon atoms were established also by GCOSY, GHMBC, GHSQC experiments. The 1 H-NMR data were reported in the following order: Chemical shifts, multiplicity, coupling constants, number of protons, and signal/atom attribution. Spin multiplets were given as s (singlet), d (doublet), t (triplet), q (quartet), m (multiplet), dd (double doublet), dt (double triplet), dq (double quartet), ddd (doublet of a double doublet), td (triple doublet), tdd (triplet of a double doublet), and br (broad) signal. The 13 C-NMR data were reported in the following order: chemical shifts, coupling constants in some cases, and signal/atom attribution ((Cq) quaternary carbon). The chemical shift of N-H protons was found in the range of δ 12.14-12.59 ppm for hydrogen bonds, and other protons of the N-H group were found in the range of δ 11.56-12.14 ppm. The hydrogen atom of the N-H group was observed downfield because of the presence of resonance and hydrogen bonds. The 13 C-NMR signals of the carbonyl group, appearing at δ 169.75-170.61 ppm were due to the existence of the intramolecular hydrogen bond related to the carbonyl oxygen atom. The carbon atom of the thiocarbonyl group at δ 179.06-182.27 ppm showed the highest values.
The molecular structure of 5d was analyzed by single crystal X-ray diffraction method, and crystallographic data are: C 23 (Table S1 from Supplementary Materials). The molecular structure ( Figure S1 from Supplementary Materials) shows that 2-(4-ethylphenoxymethyl)benzoyl fragment is bonded to the 2,4,6-trifluorophenyl thiourea moiety. The C=O and C=S bonds are oriented in a trans configuration in a plane stabilized by an intramolecular hydrogen bond N2-H···O2 ( Figure 2). This plane forms a dihedral angle of 52.58 • with the central phenyl ring. This characteristic appears in similar structures involving the benzoyl fragment connected to a thiourea fragment, as was seen by searching the Cambridge Crystallographic Data Centre (CCDC) (e.g., ethyl 4-[3-(2-methylbenzoyl)thioureido]benzoate (FUGZOX code in CCDC) [46], or N-(2-methylbenzoyl)-N -(4-nitrophenyl)thiourea (GOJNAV code in CCDC) [47].

Computational Studies
For the structure-activity relationship (SAR) studies, the following electronic properties were determined: Highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values, HOMO and LUMO orbital coefficients distribution, molecular dipole moment (for 5d compound 5.15 debye), polar surface area (for 5d compound PSA = 34.930 Å 2 ) (a descriptor that has been shown to correlate well with passive molecular transport through membranes and, therefore, allows thar predictions for the membrane transport properties), the ovality (for 5d compound 1.62), polarizability (for 5d compound 74.42 10 -30 m 3 ), and the octanol water partition coefficient (for 5d compound logP = 6.29) [48].
Properties relating to 3D structure of 5d compound, bond lengths and angles, are presented in correlation with the experimental data (Tables S4 and S5). NMR spectra of 5d compound have been View of the asymmetric unit for 5d with the atom numbering scheme (50% probability ellipsoid).

Computational Studies
For the structure-activity relationship (SAR) studies, the following electronic properties were determined: Highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values, HOMO and LUMO orbital coefficients distribution, molecular dipole moment (for 5d compound 5.15 debye), polar surface area (for 5d compound PSA = 34.930 Å 2 ) (a descriptor that has been shown to correlate well with passive molecular transport through membranes and, therefore, allows thar predictions for the membrane transport properties), the ovality (for 5d compound 1.62), polarizability (for 5d compound 74.42 10 -30 m 3 ), and the octanol water partition coefficient (for 5d compound logP = 6.29) [48].
Properties relating to 3D structure of 5d compound, bond lengths and angles, are presented in correlation with the experimental data (Tables S4 and S5). NMR spectra of 5d compound have been

Computational Studies
For the structure-activity relationship (SAR) studies, the following electronic properties were determined: Highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values, HOMO and LUMO orbital coefficients distribution, molecular dipole moment (for 5d compound 5.15 debye), polar surface area (for 5d compound PSA = 34.930 Å 2 ) (a descriptor that has been shown to correlate well with passive molecular transport through membranes and, therefore, allows thar predictions for the membrane transport properties), the ovality (for 5d compound 1.62), polarizability (for 5d compound 74.42 10 −30 m 3 ), and the octanol water partition coefficient (for 5d compound logP = 6.29) [48].
Properties relating to 3D structure of 5d compound, bond lengths and angles, are presented in correlation with the experimental data (Tables S4 and S5). NMR spectra of 5d compound have been calculated with Spartan 14 software. The experimental and calculated spectra proved to be correlated (Table S6).

Frontier Molecular Orbitals Analysis
Frontier molecular orbitals (FMOs) (the most important being HOMO and LUMO) analysis are predicting for the chemical stability and interactions between atoms and dictating the optical properties and biological activities of a certain molecule. HOMO indicates one molecule's capacity to donate, while LUMO indicates one molecule's capacity to accept an electron [49,50]. The HOMO and LUMO, calculated at the B3LYP/6-31G* level can be seen in Figure 4 for compound 5d. The "blue and red" regions of the graphic correspond to the orbital's positive and negative values. The frontier orbital gap (∆E) characterizes the chemical reactivity. The high value of the HOMO-LUMO gap (∆E) calculated for the 5d compound refers to a chemically stable molecule [51].
Molecules 2020, 25, x FOR PEER REVIEW 7 of 20 calculated with Spartan 14 software. The experimental and calculated spectra proved to be correlated (Table S6).

Frontier Molecular Orbitals Analysis
Frontier molecular orbitals (FMOs) (the most important being HOMO and LUMO) analysis are predicting for the chemical stability and interactions between atoms and dictating the optical properties and biological activities of a certain molecule. HOMO indicates one molecule's capacity to donate, while LUMO indicates one molecule's capacity to accept an electron [49,50]. The HOMO and LUMO, calculated at the B3LYP/6-31G* level can be seen in Figure 4 for compound 5d. The "blue and red" regions of the graphic correspond to the orbital's positive and negative values. The frontier orbital gap (E) characterizes the chemical reactivity. The high value of the HOMO-LUMO gap (E) calculated for the 5d compound refers to a chemically stable molecule [51].

Molecular Electrostatic Potential (MEP)
MEP has been evaluated using the B3LYP method with the basic set 6-31G* to investigate the chemical reactivity of a molecule, indicating the reactive sites of nucleophilic or electrophilic attack in hydrogen-bonding interactions involved in the biological recognition process [52,53]. The putative electrostatic map of the compound 5d shows hydrophilic regions in red (negative potential) and blue (positive potential) and hydrophobic regions in green (Figure 5a). The local ionization map is the result of the overlay of the energy of electron removal (ionization) on the electron density and can also indicate the electrophilic addition (Figure 5b). The LUMO map, representing a superposition of the absolute value of LUMO on the electron density indicates the nucleophilic addition (Figure 5c).

Molecular Electrostatic Potential (MEP)
MEP has been evaluated using the B3LYP method with the basic set 6-31G* to investigate the chemical reactivity of a molecule, indicating the reactive sites of nucleophilic or electrophilic attack in hydrogen-bonding interactions involved in the biological recognition process [52,53]. The putative electrostatic map of the compound 5d shows hydrophilic regions in red (negative potential) and blue (positive potential) and hydrophobic regions in green (Figure 5a). The local ionization map is the result of the overlay of the energy of electron removal (ionization) on the electron density and can also indicate the electrophilic addition (Figure 5b). The LUMO map, representing a superposition of the absolute value of LUMO on the electron density indicates the nucleophilic addition (Figure 5c).

Docking Studies
Molecular docking allowed us to establish an accurate prediction of the optimized conformation of new compounds (as ligands) and their target receptor protein in order to achieve a stable complex.
The score and hydrogen bonds formed with the amino acids from group interaction atoms are used to predict the binding modes, the binding affinities and the orientation of the docked substances (Figure 6a-i and Figure S4a-z) in the active site of the protein receptor (Table S7). The docking score used in the Drug Discovery Workbench is the PLANTSPLP score [52].

Docking Studies
Molecular docking allowed us to establish an accurate prediction of the optimized conformation of new compounds (as ligands) and their target receptor protein in order to achieve a stable complex.
The score and hydrogen bonds formed with the amino acids from group interaction atoms are used to predict the binding modes, the binding affinities and the orientation of the docked substances (Figure 6a-i and Figure S4a-z) in the active site of the protein receptor (Table S7). The docking score used in the Drug Discovery Workbench is the PLANTS PLP score [52].
The ligand orientation and position suggested by the molecular docking studies and the binding modes of ligands have been validated by redocking (Figure 6c,e,h, and Figure S4). According to the docking score evaluation, the best docking score was obtained for 5e: −63.91 (RMSD 0.10). The compound 5e shows three hydrogen bonds with ASN 45 (2.870 Å), ILE 78 (3.244 Å) and GLY 77 (2.924 Å) and its orientation is different from the co-crystallized one. The same orientation of the co-crystallized is shown by the the compounds 5b, 5c, 5d, and 5f (Figure 7a). The compounds 5a and 5g show also the same orientation as the 5e compound (Figure 7b).
The number of hydrogen bond donors, the number of hydrogen bond acceptors, and log P were calculated and are presented in Table 1. The calculated parameters can predict if a molecule possesses properties that might turn it into an active drug, according to the Lipinski's rule of five [53]: The number of hydrogen donors <5, the number of hydrogen acceptors <10, molecular weight <500 Da, and the octanol-water partition coefficient (log P) < 5. The log P calculation was performed using the XLOGP3-AA method [54].
Molecular docking allowed us to establish an accurate prediction of the optimized conformation of new compounds (as ligands) and their target receptor protein in order to achieve a stable complex.
The score and hydrogen bonds formed with the amino acids from group interaction atoms are used to predict the binding modes, the binding affinities and the orientation of the docked substances (Figure 6a-i and Figure S4a-z) in the active site of the protein receptor (Table S7). The docking score used in the Drug Discovery Workbench is the PLANTSPLP score [52]. The ligand orientation and position suggested by the molecular docking studies and the binding modes of ligands have been validated by redocking (Figure 6c, 6e, 6h, and S4). According to the docking score evaluation, the best docking score was obtained for 5e: −63.91 (RMSD 0.10). The compound 5e shows three hydrogen bonds with ASN 45 (2.870 Å), ILE 78 (3.244 Å) and GLY 77 (2.924 Å) and its orientation is different from the co-crystallized one. The same orientation of the co-crystallized is shown by the the compounds 5b, 5c, 5d, and 5f (Figure 7a). The compounds 5a and 5g show also the same orientation as the 5e compound (Figure 7b)  The ligand orientation and position suggested by the molecular docking studies and the binding modes of ligands have been validated by redocking (Figure 6c, 6e, 6h, and S4). According to the docking score evaluation, the best docking score was obtained for 5e: −63.91 (RMSD 0.10). The compound 5e shows three hydrogen bonds with ASN 45 (2.870 Å), ILE 78 (3.244 Å) and GLY 77 (2.924 Å) and its orientation is different from the co-crystallized one. The same orientation of the co-crystallized is shown by the the compounds 5b, 5c, 5d, and 5f (Figure 7a). The compounds 5a and 5g show also the same orientation as the 5e compound (Figure 7b) The results for 5d compound revealed a docking score -56.80 (RMSD 0.20) and two hydrogen bonds with GLY 101 (3.144 Å and 3.381 Å) (Figure 6d).

Antimicrobial Activity Assay
The synthesized compounds were tested by a quantitative method, allowing us to establish the MIC (minimum inhibitory concentration) value ( Table 2). The compounds 5c and 5f proved to have a low inhibitory effect against the tested strains. The compounds encoded 5a and 5b proved to be active against E. coli and P. aeruginosa strains, exhibiting the same MIC value. The compound 5d exhibited both antibacterial and antifungal activity, being active on E. coli and C. albicans strains. The compound 5a exhibited the largest spectrum of antibacterial activity, being active both against Gram-positive and −negative strains. Interestingly, the largest spectrum of activity of the compound 5a was correlated with the presence of one fluorine atom as substituent on the phenyl ring. Table 2. The inhibitory activity of the new F-compounds 5a-g on the planktonic growth of the tested microbial strains. The minimum inhibitory concentration (MIC) expressed in µg/mL of the new fluoro compounds (5a-g) are presented as the average of three different experiments. The most active compounds are highlighted. Compared to the Gram-positive strains, the Gram-negative bacterial strains proved to be more susceptible to the tested compounds, with five of the seven tested compounds being active against the E. coli strain. S. aureus and B. subtilis proved to be resistant to all tested compounds, while E. faecalis and C. albicans were susceptible to only one compound, respectively, 5a and 5d ( Table 2).
The presence of multiple fluorine atoms (e.g., three atoms in the compounds 5b, 5c, 5d, and 5e) located at different positions in the compounds did not correlate with the intensity of the antibacterial effect, with three of them exhibiting the same anti-E. coli activity similar to 5a, with the other one displaying a lower activity against the same strain. All of them were inactive on E. faecalis and all but 5b against P. aeruginosa when compared with 5a. The compounds substituted with the trifluoromethyl group (5e and 5g) were active on E. coli. Further, the presence of three fluorine atoms in the positions 2, 4, and 6 conferred antifungal properties for the compound 5d.
The antibiofilm activity spectrum was different from that obtained on planktonic strains, with all tested compounds being able to inhibit the formation of biofilms formed by at least two of the tested strains. The E. coli, S. aureus, and B. subtilis biofilms were the most recalcitrant to the tested compounds, with a maximum of three compounds exhibiting antibiofilm activity ( Table 3). The E. faecalis and C. albicans biofilms were susceptible to the tested compounds, excepting 5c and 5g. The most active compounds on the development of microbial biofilms proved to be those substituted by the trifluoromethyl group, the most favorable substitution being in para position (Table 3).

Discussion
The emergence and spread of antimicrobial resistance have become some of the major worldwide health problems today, raising the acute need for developing novel antimicrobial agents. Considering the valuable pharmacological properties of thiourea derivatives, our aim was to design, synthesize and evaluate new compounds with fluorinated and trifluoromethylated aryl thiourea scaffolds. Fluorine substitution can improve the efficacy of drugs by influencing the absorption, tissue distribution, the rate and the route of biotransformation, the pharmacodynamic and the toxicological properties. Our interest was to optimize the substitution atoms or the functional groups present on the phenyl bound to nitrogen atom of thiourea in order to increase the antimicrobial activity of the obtained derivatives.
The best antibacterial activity (against E. coli, P. aeruginosa, and E. faecalis) was recorded for the compounds bearing one fluorine atom, while the antifungal effect was favored by the isomeric substitution with three fluorine atoms. The 5d exhibited the lowest recorded MIC value against the C. albicans fungal strain. This effect could be explained by the enhancement of lipid solubility and permeability of the fungal cell wall, with a particular structure owing to the eukaryotic nature of this microorganism, for the antimicrobial compound.
The spectrum of the antibiofilm activity was different from that recorded on planktonic cells, with all tested compounds inhibiting the ability of some microbial strains to form biofilms.The intensity of this activity was correlated with the presence of the trifluoromethyl substituent, most favorably in para position and the compound 5g being the most active.
Despite the great susceptibility of planktonic E. coli cells to the majority of the tested compounds, the adherent E. coli cells were more resistant, remaining susceptible to only two compounds (5e, 5g). This could account for the increased tolerance of microbial biofilms to the usual doses of antibiotics.
It is worth noting that P. aeruginosa biofilms were highly susceptible to the all tested compounds, except 5a.
Some electronic and molecular properties were discussed in order to assess the flexibility and the binding ability of the studied conformer to bind to the receptor protein. The docking studies using E. coli DNA gyrase B revealed that all compounds showed a good docking score, indicating the inhibition of DNA replication as a potential mechanism of their antimicrobial activity that is worth to be studied further.

Chemistry
Highest quality phthalide, p-ethylphenol, ammonium thiocyanate, primary aromatic amines, and thionyl chloride used for the synthesis were purchased from Sigma-Aldrich Chemical Co. Analytical grade solvents (1,2-dichloroethane, xylene, isopropanol, acetone) supplied by Sigma-Aldrich were used and were dried when necessary. Potassium hydroxide and hydrochloric acid were obtained from Merck, and potassium carbonate from Fluka Chemical.
Acetone was dried over K 2 CO 3 and then distilled. Ammonium thiocyanate was heated at 100 • C before use. The liquid amines were dried with potassium hydroxide and then distilled.
The reaction progress and the purity were investigated by thin layer chromatography on silica gel 60 F254 (0.5 mm thick) plates (Merck, Germany) and unidimensional migration using ethyl acetate/cyclohexane (4:6 v/v) as eluent with visualization by ultraviolet light (λ = 254 nm) and exposed to iodine vapors.
Melting points were measured by heating a small amount of sample, by the capillary tube method, using an Electrothermal 9100 apparatus (Bibby Scientific Ltd., Stone, UK), without thermometer correction. The tube was observed with a magnifying glass, with a temperature range of 0.5-1.0 • C indicating a relatively high level of purity.
C, H, N, and S analyses were performed on a Perkin Elmer 2400 Series II CHNS/O Elemental Analyzer (Waltham, MA, USA), using the classical Pregl-Dumas method, and considering an accepted deviation of elemental analysis results from the calculated one of 0.4%.
The FTIR spectra were recorded with a Bruker Vertex-70 spectrophotometer (Bruker Corporation, Billerica, MA, USA) with an optical system with diamond. The IR bands were given as w (weak), m (medium), s (intense), and vs (very intense) and were obtained using attenuated total reflection Fourier-transformed infrared (ATR-FTIR) spectra at room temperature.
NMR spectra were recorded on a Varian Unity Inova 400 instrument operating at 400 MHz for 1 H, 100 MHz for 13  As solvent, we used perdeuterated dimethyl sulfoxide with a min. 98% isotopic purity. The spectra were recorded at a temperature of 20 ± 1 • C. In order to perform the right attribution of spectral signals, we used different techniques to simplify the spectra, like the homonuclear decoupling and the specific deuteration.

General Synthesis Procedure of the New Benzoylthiourea Derivatives
To a solution of ammonium thiocyanate (0.01 mol) in acetone (5 mL) we added a solution of 2-(4-ethylphenoxymethyl)benzoyl chloride (3) (0.01 mol) in acetone (15 mL) to obtain in situ arylisothiocyanate (4). The reaction mixture was refluxed for 1 h, and then cooled at room temperature. A primary amine solution (0.01 mol) in acetone (2 mL) was added to the mixture and refluxed for 1 h. After cooling, the compound was precipitated and poured into 500 mL cold water, the solid was filtered, dried, and then recrystallized from isopropanol with active carbon.
The steps to go through to explore protein-ligand interaction are: Set up the binding site in a molecule project, dock the ligands imported to a molecule table, and inspect the docking results.

In Vitro Assessment of the Antimicrobial Activity of the Newly Synthesized F-Benzoylthiourea Derivatives on Planktonic Microbes (Grown in Suspension)
The inhibitory activity upon microbial growth was tested on Gram-positive (S. aureus ATCC 25923, E. faecalis ATCC 29212, B. subtilis ATCC 6633) and Gram-negative (E. coli ATCC 25922, P. aeruginosa ATCC 27853) bacteria and one fungal (C. albicans ATCC 10231) strain. The compounds were solubilized in DMSO to a final concentration of 10 mg/mL and further sterilized by filtration (using 0.22 µm filter membranes). The quantitative assay of the minimal inhibitory concentration (MIC, µg/mL) was performed by the microdilution method in 96 multi-well plates [59]. In this purpose, serial binary dilutions (ranging between 1024 and 8 µg/mL) were achieved in nutrient broth/YPG (yeast peptone glucose) and each well was seeded with microbial inoculum of 0.5 McFarland density. The plates were incubated for 24 hours at 37 • C in the case of bacterial strains, and 48 hours at 28 • C in the case of the fungal strain. The respective MIC values were reported as being the lowest concentration of the tested compound which inhibited the visible microbial growth.

Antimicrobial Activity of Newly Synthesized F-Benzoylthiourea Derivatives on Microbial Adherence and Biofilms' Development on the Inert Substratum
The antibiofilm activity was assessed by the microtiter method. For this purpose, after reading the MIC, the microplates were emptied and the microbial biofilms formed on the plastic walls were fixed with cold methanol, stained with 1% violet crystal solution, and resuspended in 30% acetic acid. The absorbance of the colored suspension was measured at 490 nm with an ELISA reader (Apollo LB 911). The inhibition of the microbial biofilm development was interpreted as the decrease of the absorbance below the value registered for the positive control [60][61][62][63][64].

Conclusions
We report here the synthesis, characterization, and bioevaluation of new fluorinated compounds as potential antimicrobials. The antibacterial and antifungal activity was correlated with the position of fluorine or trifluoromethyl substituents on the phenyl ring. The molecular docking simulation was performed to predict the binding modes of the obtained benzoylthiourea derivatives to E. coli DNA gyrase B (chosen because E. coli strain was the most susceptible to the highest number of the tested compounds), the binding affinities, and the orientation. The docking score, predicting the antimicrobial activity of the studied compounds, was correlated with the experimental data obtained from the evaluation of the antimicrobial activity against the E. coli strain.

Conflicts of Interest:
The authors declare no conflict of interest