Abstract

This paper reports on the synthesis of new thiazole derivatives that could be profitably exploited in medical treatment of tumors. Molecular electronic structures have been modeled within density function theory (DFT) framework. Reactivity indices obtained from the frontier orbital energies as well as electrostatic potential energy maps are discussed and correlated with the molecular structure. X-ray crystallographic data of one of the new compounds is measured and used to support and verify the theoretical results.

1. Introduction

The thiazole has an important component effect of the pharmacophores of a large number of medicinal significance molecules and the evaluation of their biological activity, such as antibacterial [1], antiprotozoal [2], antitubercular [3], antifungal [4, 5], and anthelmintic [6], with emphasis on their potential medicinal applications, is desirable. Here we are interested to study newly synthetized aminothiazoles, especially 2-aminothiazole derivatives which represent a class of heterocyclic ring system possessing antiviral [7], antimicrobial [8], anticancer [9], and anti-inflammatory activities [10]. Previously, in vitro anticancer evaluation studies of different 2-aminothiazole analogs exhibited their potent and selective nanomolar inhibitory activity against a wide range of human cancerous cell lines such as breast, leukemia, lung, colon, CNS, melanoma, ovarian, renal, and prostate cell lines [1114]. Substitutions at 2-position benzothiazole have emerged in its usage as a core structure in the diversified therapeutic applications [1521]. The studies of structure-activity relationship interestingly reveals that change of the structure of substituent group at C-2 position commonly results in the change of its bioactivity. Though literature survey reports many therapeutic applications of 2-substituted benzothiazoles, their investigation for anti-inflammatory activity is limited [16, 2225]. Furthermore, thiazole derivatives have attracted a great deal of interest due to their wide applications in the field of pharmaceuticals. Thiazole derivatives display a wide range of biological activities such as cardiotonic, fungicidal, sedative, anesthetic, bactericidal, and anti-inflammatory [26, 27]. In addition, thiazole derivatives are reported to show a variety of biological activities. Depending on the substituents, this heterocycle possesses anthelmintic, antibiotic, and immunosuppressant activity [28]. Recent research indicates that some of 2-aminothiazoles derivatives are inhibitors of enzymes such as kynurenine-3-hydroxylase 29 or possess inhibitors activity against enzyme cyclin-dependent kinase [29].

Additionally, monoazo disperse dyes with thiazole-diazo components have been intensively investigated to produce bright and strong color shades ranging from red to greenish blue on synthetic fiber. Color Index described various basic, direct, vat, and disperse dyes wherein thiazole nucleus occurs [30]. Derivative of 2-aminothiazole has a long history of use as heterocyclic diazo components for disperse dyes [31].

In the present study, quantum chemical computations will be performed within DFT using WB97XD/6-31G(d) model to investigate the molecular structure, IR, and NMR of the newly synthetized molecules [3238]. X-ray crystallographic data of 3B will be obtained and used to support and verify the theoretical results. The energies of HOMO-LUMO frontier orbitals will be used to estimate molecular reactivation towards nucleophilic/electrophilic reagents. Electrostatic potential energy maps (ESP maps) will be graphically presented to locate binding sites of these new derivatives.

2. Experimental Section

2.1. Synthesis

All melting points were measured on an electric melting point apparatus and were uncorrected. The infrared spectra were recorded using potassium bromide disks on a Pye Unicam SP-3-300 infrared spectrophotometer; the established values of the gas phase frequencies are given between brackets. 1HNMR spectra were run at 300 MHz, on a Varian Mercury VX-300 NMR spectrometer and Brukeravance III 400 MHZ, using TMS as an internal standard in deuterated dimethylsulphoxide. Chemical shifts are quoted in ppm. The mass spectra were recorded on Shimadzu GCMS-QP-1000EX mass spectrometers at 70 eV. All the spectral measurements were carried out at the NMR Laboratory of Cairo University, Egypt, and the NMR Laboratory of Faculty of Pharmacy, Ain shams University, Egypt; the microanalytical data were measured in Central Lab of Cairo University, Egypt; the Ministry of Defense Chemical Laboratories, Egypt; and the Microanalytical Center of Ain Shams University, Egypt. All the chemical reactions were monitored by TLC. The bold values corresponded to values calculated from DFT.

2.1.1. General Procedure for the Preparation of Compounds 2ag

A mixture of 2-aminothiazole 1 (1 g, 10 mmol) and different electrophilic reagents, namely, 2-chloro-N-(4-sulfamoylphenyl) acetamide, ethyl chloroacetate, phenyl cyanate, chloroacetyl chloride, 2-chloro-N-(4-chlorophenyl) acetamide, phenyl cyanate, phenyl thiocyanate, chloroacetyl chloride and 2-chloro-N-(4-chlorophenyl) acetamide (10 mmol) in dimethylformamide (20 ml), and anhydrous potassium carbonate, was refluxed for 5–8 h. The reaction mixture was poured in ice water (200 ml); the formed ppt was filtered off, dried, and crystallized from ethanol afforded compounds 2ag. It should be noticed that quantum chemically calculated spectroscopic parameters in gaseous phase are given in for comparison with the experimentally obtained parameters.

N-(4-Sulfamoylphenyl)-2-(thiazol-2-ylamino) Acetamide (2a). Yield 65%; m.p. 118–120°C; orange crystals; (EtOH); IR (KBr) broad band at 3378, 3325, 3273 cm−1 (), 3022 cm−1 (), 2958 cm−1 (), 1691 cm−1 (). 1HNMR (300 MHz, DMSO-d6) δ ppm: 5.1 (s, 2H, CH2), 7.15 (d, 1H, thiazole H, J = 8.6 Hz), 7.4 (d, 1H, thiazole H, J = 8.6 Hz), 7.2 (s, 2H, NH2, D2O exchangeable), 7.7 (d, 2H, Ar-H), 7.8 (d, 2H, Ar-H), 10.8 (s, 1H, D2O exchangeable NH), 7.1 and 12 (br s, 1H, D2O exchangeable NH-OH tautomerism). Anal. Calculated for C11H12N4O3S2 (312.36): C, 42.30; H, 3.87; N, 17.94 Found: C, 42.01; H, 4.07; N, 17.84.

Ethyl thiazole-2-ylglycinate (2b). Yield 80%; m.p. 142–144°C; orange crystals; (EtOH); IR (KBr) 3346 cm−1 (), 3045 cm−1 (), 1730 cm−1 ester (). 1HNMR (300 MHz, DMSO-d6) δ ppm: 1.4 (trip, 3H, CH3), 3.2 (quar, 2H, CH2), 5.1 (S, 2H, CH2), 7.0 (d, 1H, thiazole H, J = 8.4 Hz), 7.4 (d, 1H, thiazole H, J = 8.4 Hz), 10.2 (s, 1H, NH, D2O exchangeable), Anal. Calculated for C7H10N2O2S (186.23): C, 45.15; H, 5.41; N, 15.04, Found: C, 45.01; H, 5.37; N, 14.84.

1-Phenyl-3-(thiazole-2-yl) Urea (2c). Yield 78%; m.p. 125–127°C; yellow crystals; (EtOH/benzene); IR (KBr) 3326, 3284 cm−1 (), 3063 cm−1 (), 1648 cm−1 (). 1HNMR (300 MHz, DMSO-d6) δ ppm: 7.00 (d, 2H, J = 7.2 HZ) 7.4 (m, 3H, ArH), 7.1 (d, 1H, thiazole H, J = 8.0 Hz), 7.8 (d, 1H, thiazole H, J = 8.0 Hz), 8.9 (s, 1H NH, D2O exchangeable,), 10.6 (s, 1H, NH, D2O exchangeable), 13CNMR (300 MHz, DMSO-d6) δ ppm: 90.3, 100.2, 105.7, 120.8, 138.5, 140.9, 159.3, 165.4, Anal. Calculated for C10H9N3OS (219.26): C, 54.78; H, 4.14; N, 19.16, Found: C, 54.52; H, 4.17; N, 18.94.

1-Phenyl-3-(thiazole-2-yl) Thiourea (2d). Yield 76%; m.p. 178–180°C; brownish red crystals; (butanol); IR (KBr) 3169, 3081 cm−1 (), 3009 cm−1 (). 1HNMR (300 MHz, DMSO-d6) δ ppm: 7.00 (d, 2H, J = 7.2 HZ), 7.4 (m, 3H, ArH), 7.4 (d, 1H, thiazole H, J = 7.8 Hz), 7.6 (d, 1H, thiazole H, J = 7.8 Hz), 10.4 (s, 1H, NH, D2O exchangeable), 12.4 (s, 1H, NH, D2O exchangeable), 13CNMR (300 MHz, DMSO-d6) δ ppm: 125.3, 126.4, 128.7, 129.4, 111.2, 137.4, 162.4, 176.6, Anal. Calculated for C10H9N3S2 (235.32): C, 51.04; H, 3.86; N, 17.86, Found: C, 51.01; H, 4.07; N, 17.74.

2-Chloro-N-(thiazole-2-yl) Acetamide (2e). Yield 60%; m.p. 220–222°C; brownish red crystals; m.p. 162–164°C (butanol); IR (KBr) 3187 cm−1 (), 3041 cm−1 (), 1703 cm−1 (). 1HNMR (300 MHz, DMSO-d6): δ ppm: 4.5 (s, 2H, CH2), 7.2 (d, 1H, thiazole H, J = 8.2 Hz), 7.6 (d, 1H, thiazole H, J = 8.2 Hz), 12.4 (s, 1H, NH, D2O exchangeable), 13CNMR (300 MHz, DMSO-d6) δ ppm: 39.7, 114.5, 137.7, 157.8, 157.6, 164.4, Anal. Calculated for C5H5ClN2OS (176.62): C, 34.00; H, 2.85; N, 15.86, Found: C, 34.12; H, 3.07; N, 15.74.

N-(4-Chlorophenyl)-2-(thiazole-2-ylamino) Acetamide (2f). Yield 68%; m.p. 198–200°C; red crystals; (butanol); IR (KBr) 3292, 3184 cm−1 (), 3047 cm−1 (), 1694 cm−1 (). 1HNMR (300 MHz, DMSO-d6) δ ppm: 2.4, 2.6 (s, 2H, 2NH, D2O exchangeable), 7.1 (s, 2H, CH2), 7.2 (d, 2H, J = 7.6 Hz), 7.6 (m, 3H, ArH), 8.0 (d, 1H, thiazole H, J = 8.2 Hz), 8.2 (d, 1H, thiazole H, J = 8.2 Hz) Anal. Calculated for C11H10ClN3OS (267.73): C, 34.35; H, 3.76; N, 15.70, Found: C, 34.12; H, 3.67; N, 15.74.

2.1.2. General Procedure for the Preparation of Compounds 3ae

A mixture of 2e (1.76 g, 10 mmol) and different nucleophilic reagents, namely, p-toluidine, thiourea, anthranilic, aminothiophenol, and quinoxaline-2,3-diol (10 mmol) in dimethylformamide (20 ml) and anhydrous potassium carbonate was refluxed for 4-5 h. The reaction mixture was poured onto ice water (200 ml). The formed precipitate was filtered off and recrystallized from the suitable solvent to afford compounds 3a–e.

3-(Thiazole-2-ylamino) Benzo[e][1,4]oxazepin-5(1H)-one (3a). Recrystallized from ethanol to produce brownish red crystals, Yield 84%; m.p. 252–254°C; (butanol); IR (KBr) 3672, 3377 cm−1   asym, sym (), 3198 cm−1 asym  str (), 3020 cm−1 (), 2967.42 cm−1 sym  str () 1698 cm−1 (). 1HNMR (300 MHz, DMSO-d6) δ ppm: 7.2 (d, 1H, thiazole H, J = 8.2 Hz), 7.4 (d, 1H, thiazole H, J = 8.2 Hz), 7.6 (d, 2H, ArH, J = 8.1 Hz), 7.8 (m, 2H, ArH), 8.2 (s, 1H, NH, D2O exchangeable), 8.8 (s, 1H, 1-ethelyene), 11.4 (s, 1H, NH, D2O exchangeable), Anal. Calculated for C12H9N3O2S (259.28): C, 55.59; H, 3.50; N, 16.21, Found: C, 55.45; H, 3.67; N, 16.14.

N-(Thiazole-2-yl)-2-(p-tolylamino) Acetamide (3b). Oily product, solidified with diethyl ether, Yield 44%; m.p. 235–237°C; brownish red crystals; (butanol); IR (KBr) 3287, 3150 cm−1 (), 3045 cm−1 (), 1701 cm−1 (). 1HNMR (300 MHz, DMSO-d6) δ ppm: 2.11 (s, 3H, CH3), 3.5 (s, 2H CH2), 6.25 (6.32) (d, 2H, Ar benz, J = 6.8 Hz), 6.41 (6.52) (d, 1H, thiazole H, J = 8.2 Hz), 7.0 (7.1) (d, 2H, ArH, J = 8.1 Hz), 7.3 (7.4) (d, 1H, thiazole H, J = 8.2 Hz), 8.01 (7.8) (d, 2H, ArH, J = 6.8 Hz), 11.4 (s, 1H, NH, D2O exchangeable), Anal. Calculated for C12H13N3OS (247.32): C, 58.28; H, 5.30; N, 16.99, Found: C, 58.45; H, 5.67; N, 16.80.

4-(Thiazole-2-ylamino)-1,5-dihydro-2H-imidazole-2-thione (3c). Dried and crystallized from ethanol the product has brownish red crystals; m.p. 206–208°C; Yield 54%; (butanol); IR (KBr) 3374, 3275 cm−1 (), 3045 cm−1 (), 2600 cm−1 (), 1260 cm−1 (), 1HNMR (300 MHz, DMSO-d6) δ ppm: 4.22 (s, 2H, imidazole H), 6.9 (d, 1H, thiazole H, J = 7.4 Hz), 77.2 (d, 1H, thiazole H, J = 7.4 Hz), 9.5 (s, H, NH, D2O exchangeable), 13.2 (s, H, imidazole NH, D2O exchangeable), Anal. Calculated for C6H6N4S2 (198.00): C, 36.35; H, 3.05; N, 28.26 Found: C, 36.45; H, 3.25; N, 28.32.

N-(Thiazole-2-yl)-4H-benzo[b][1,4]thiazine-2-amine (3d). Dried and crystallized from butanol, Yield 62%; brownish red crystals; m.p. 262–266°C; (butanol); IR (KBr) 3445.47, 3310.31 cm−1 str,  hetero  ring,    str,   link  NH (), 3198.45 cm−1 (), 3063 cm−1 (), 1663.89 cm−1 (), 1HNMR (300 MHz, DMSO-d6) δ ppm: 5.1 (s, H, Thiazine H), 7.0 (d, 2H, ArH, J = 7.9 Hz), 7.4 (m, 2H, ArH), 7.5 (d, 1H, thiazole H, J = 8.2 Hz), 7.6 (d, 1H, thiazole H, J = 8.2 Hz), 10.8, 12.2 (s, 2H, 2NH, D2O exchangeable), 13CNMR (300 MHz, DMSO-d6) δ ppm: 114.6, 122.2, 123.7, 124.6, 126.8, 141.0, 153.6, 155.7, 165.4, Anal. Calculated for C11H9N3S2 (247.33): C, 53.42; H, 3.67; N, 16.99, Found: C, 36.45; H, 3.25; N, 28.32.

2-((3-Hydroxyquinoxalin-2-yl)oxy)-N-(thiazole-2-yl) Acetamide (3e). Crystallized from ethanol, Yield 58%; m.p. 298–300°C; brownish red crystals; (butanol); IR (KBr) broad band at 3450 cm−1 of  hetero  ring (), 3342 cm−1 of  hetero  ring (), 3190 cm−1 (), 3120 cm−1 (), 3048 cm−1 (), 1725 cm−1 (), 1658 cm−1 (). 1HNMR (300 MHz, DMSO-d6) δ ppm: 5.4 (s, 2H, CH2), 7.4 (s, 1H, NH, D2O exchangeable), 7.5 (d, 2H, ArH, J = 8.2 Hz), 7.7 (m, 2H, ArH), 8.0 (d, 1H, thiazole H, J = 8.56 Hz), 8.2 (d, 1H, thiazole H, J = 8.56 Hz), 13.2 (s, 1H, OH, D2O exchangeable), Anal. Calculated for C13H10N4O3S (302.31): C, 51.65; H, 3.33; N, 18.53, Found: C, 51.56; H, 3.21; N, 18.56.

2.2. X-Ray Crystallography

X-ray structure analysis offers perfect addition to our synthetic work. X-ray structures of the compound 3b were performed in the Central Service and X-Ray Laboratories, National Research Centre, Cairo, Egypt. Crystal and molecular structures were prepared by Maxus Computer Program for the Solution and Refinement of Crystal Structures. All diagrams and calculations were performed using maXus (Bruker Nonius, Delft; MacScience, Japan). There was no extinction correction. Atomic scattering factors were from Waasmaier and Kirfel, 1995. Data collection parameters are as follows: KappaCCD; cell refinement: HKL Scalepack; data reduction: Denzo Program(s) used to solve structure; SIR92 and Scalepak Program(s) used to refine structure; maXus: ORTEP Software which was used for molecular graphics. Crystal data, fractional atomic coordinates, and equivalent isotropic thermal parameters, anisotropic displacement parameters and geometric parameters of compounds 3b are given in Table 2. The additional data for the molecule 3b are alternatively available from the Cambridge Crystallographic Data Centre as CCDC1402910.

2.3. Computations

Computations were carried out using Gaussian 16 revision A.03 package [32] and/or Spartan’16 parallel QC program [Wavefunction, Inc., USA]. Optimized structures and spectroscopic data derived from quantum chemical calculations have been used within the WB97DX/6-31G(d) model. A Broadberry (UK) 40-core workstation and/or MAC Pro 12-core computers were used.

3. Result and Discussion

3.1. Synthesis and Spectroscopic Properties

In our study, 2-aminothiazole 1 was used as a key starting material. Reaction of 1 with chloro-N-(4-sulfamoylphenyl) acetamide afforded the amide derivative 2a (Scheme 1). The structure of 2a is substantiated from its spectral data. The IR spectrum shows appearance of absorption band of C=O group for the amide at 1691 cm−1, as well as the presence of OH-NH tautomerization at δ 7.1 ppm and 12 ppm. On the other hand, when 1 was refluxed in dimethylformamide with ethyl chloroacetate, the ester 2b was abstained and its structure was confirmed with different spectral data: the presence of the ester C=O at 1730 cm−1 in IR, for example, and the presence of CH2CH3 in H-NMR as quartet and triplet at δ 3.2, 1.4 ppm, respectively.

In addition, urea and thiourea derivatives 2c, 2d were obtained from the reaction of the aminothiazole with phenyl isocyanate and phenyl isothiocyanate. The IR spectrum revealed the absence of doublet bands of NH2 in both compounds, the appearance of band that is attributed to C=O for 2c at 1648 cm−1, and the appearance of four peaks that is attributed to phenyl ring in C13-NMR. The most important compound in this work is compound 2e that resulted from interaction of 1 with chloroacetyl chloride; the structure was proved by appearance of C=O at 1703 cm−1 as well as absence of NH2 doublets. The amide derivative 2f is obtained from reaction of thiazole derivative with 2-chloro-N-(4-chlorophenyl) acetamide; the IR spectrum shows the presence of C=O band and the presence of double doublets of para-substituted-benzene ring of chlorophenyl in H-NMR.

The compound 2e was the key start for many other reactions; refluxing 2e with anthranilic acid afforded the oxazipin-one 3a (Scheme 2). The cyclic structure was proved from IR spectrum which showed the absence of broad OH band that is attributed to open structure and the appearance of C=O band at 1698 cm−1.

On the other hand, the amide derivatives 3b were obtained from reaction of 2e with p-toluidine for five hours; the open structure was confirmed with many tools as IR which show two NH bands at 3287, 3150 cm−1, as well as the X-Ray crystallography; furthermore, refluxing of chloroderivative 2e with thiourea afforded cyclic structure 3c, which had been proved with absence of C=O band in IR and appearance of 3374 and 3275 cm−1 for 2 NH. Furthermore, the appearance of weak band as 2600 cm−1 is attributed to thione-thiol SH tautomerization. The thiazine 3d is another cyclic compound resulting from refluxing 2e with 2-aminophenol; the structure was proved by disappearance of C=O, as well as the appearance of peak in H-NMR for thiazine H at δ 5.1 ppm and appearance for extra peak at δ 114 ppm for thiazine ring in C13-NMR.

At last, refluxing 2e with quinoxaline-2,3-diol in DMF/anhydrous carbonate produced the ether 3e, whose structure was evaluated from IR by peaks at 1725 cm−1 for C=O, appearance of broad band at 3350 cm−1 that is attributed to OH, and, in addition, aromatic peaks in H-NMR at δ 7.5 ppm and 7.7 ppm.

3.2. X-Ray Crystallography and Optimized Molecular Structure

X-ray results are depicted in Table 1. X-ray structure analysis offers perfect addition to our synthetic work. X-ray structures of the compound 3b (Figure 1) showed that the molecule is planar. Table 2 shows the agreement between the optimized parameters and the experimentally obtained geometry of 3B molecule.

The structure produced (Figure 1 and Table 2) is in excellent match with the optimized structure obtained by quantum chemical calculations within the density functional theory (DFT) [33, 34] using WB97XD/6-31G(d) model.

3.3. Molecular Reactivities

Chemical reactivity theory quantifies the reactive propensity of isolated species through the introduction of a set of reactivity indices or descriptors. Its roots go deep into the history of chemistry, as far back as the introduction of such fundamental concepts as acid, base, Lewis acid, and Lewis base. It pervades almost all of chemistry.

The most relevant indices defined within the conceptual DFT [33] for the study of the organic reactivity are discussed elsewhere [3539]. Molecular reactivity indices [3539] such as chemical potential (), hardness (), and electrophilicity () were computed from the energies of frontier orbitals and defined as follows:(1)Chemical potential is given byor simply = (LUMO + HOMO).(2)Hardness is given byor simply (LUMO − HOMO). The chemical hardness can be thought as a resistance of a molecule to exchange electron density with the environment.(3)Electrophilicity: in 1999, Parr defined the electrophilicity index [40] , which measures the total ability to attract electrons. The electrophilicity index gives a measure of the energy stabilization of a molecule when it acquires an additional amount of electron density from the environment. The electrophilicity index comprises the tendency of an electrophile to acquire an extra amount of electron density, given by and the resistance of a molecule to exchange electron density with the environment, given by . Therefore, a good electrophile is a species characterized by a high absolute value and a low value. The electrophilicity index has become a powerful tool for the study of the reactivity of organic molecules [36].(4)Nucleophilicity (N): while the electrophilicity of the molecules accounts for the reactivity towards nucleophiles, it has been shown by Domingo and his coworkers [3639] that a simple index chosen for the nucleophilicity, , based on the HOMO energy, within DFT, is useful to explain the reactivity of these new compounds towards electrophiles.Nucleophilicity index is defined as , where -9.12 is the energy of the HOMO of tetracyanoethylene (TCE). It is noteworthy to mention that this nucleophilicity scale is referred to tetracyanoethylene (TCE) taken as a reference, because it presents the lowest HOMO energy in a large series of molecules investigated [36].

The numerical parameters reflect the tendency of transferring electronic charges during chemical interactions between molecules (Table 3). Electrophilicity is an important reactivity descriptor that is considered as a measure of a compound’s willingness to participate as an electron acceptor during a chemical reaction, or, in other words, the electron deficiency of a compound. 3E is the molecule with the largest electrophilicity , whereas 3D is the one showing smallest electrophilicity indicating lower susceptibility towards nucleophilic reaction. Table 3 shows that chemical potential value, which is the negative of molecular electronegativity reflecting the escaping tendency of electrons, decreases in the following order: 3D > 2B > 3B > 2C > 3A > 2F > 2D > 2E > 2A > 3C > 3E.

By examining the nucleophilicity descriptor (Table 3) for these molecules, we found that 2E ( = 0.78 eV) is one of the poorest nucleophiles of this series, while 3D ( = 2.97 eV) represents the best nucleophile. Generally, nucleophilicity increases in the following order: 2E < 2A < 3C < 3E < 2C < 2B < 2F < 2D < 3B < 3A < 3D

These results are consistent with the expected reactivity pattern.

Investigation of a molecule’s surface is probably a good start for considerations of the molecule’s reactivity since this is where two approaching molecules would first interact. ESP maps are depicted in Figure 2. The results should improve our knowledge about the binding sites, which are of importance in chemical reactivities and medical applications. Color codes point to the binding sites when interacting with other reagents [4143] (see caption of Figure 2).

4. Conclusions

One-step syntheses of 12 thiazole derivatives of medicinal importance are performed. Optimized structures, reactivity indices, electrostatic potential energy maps, and spectroscopic properties such as IR and NMR of the newly reported molecules are computed within DFT using WB97XD/6-31G(d) model. Satisfactory agreement between experiment and theory is observed. Trends in chemical reactivities are investigated. Molecules 3E and 3D have the largest and smallest electrophilicity, respectively. Generally, based on relative nucleophilicity index , nucleophilicity increases in the following order: 2E < 2A < 3C < 3E < 2C < 2B < 2F < 2D < 3B < 3A < 3D. These results are consistent with the expected reactivity pattern.

The graphically visualized ESP maps enable locating the binding sites of these molecules.

Conflicts of Interest

The authors declare that they have no conflicts of interest.