- 1Department of Chemistry, College of Science, King Faisal University, Al Hofuf, Saudi Arabia
- 2Department of Chemistry, Organic Chemistry Division, College of Science, Mansoura University, Mansoura, Egypt
- 3Department of Chemistry, College of Science, IMSIU (Imam Mohammad Ibn Saud Islamic University), Riyadh, Saudi Arabia
- 4Chemistry Department, Faculty of Science, Sohag University, Sohag, Egypt
- 5Institute of Chemistry, Martin Luther University Halle-Wittenberg, Halle (Saale), Germany
- 6Department of Chemistry, Faculty of Science, Cairo University, Giza, Egypt
- 7Toxic and Narcotic Drug, Forensic Medicine Department, Mansoura Laboratory, Medicolegal Organization, Ministry of Justice, Cairo, Egypt
Reduction of 4,4′-diselanediyldianiline (1) followed by the reaction with bromo-4-(bromomethyl)benzene afforded the corresponding 4-((4-bromobenzyl)selanyl)aniline (2) in 85% yield. N-Maleanilic acid 3 was obtained in 94% yield via the reaction of selenoamine 2 with toxilic anhydride. Subsequent dehydration of N-maleanilic acid 3 using acetic anhydride furnished the unexpected isomaleimide 5-((4-((4-bromophenyl)selanyl)phenyl)imino)furan-2(5H)-one (4) instead of the maleimide 5. The molecular structure of compound 4 was confirmed by mass spectrometry, 1H- and 13C-NMR spectroscopy, and X-ray diffraction analysis. Their cytotoxicity was assessed against two oligodendrocytes, and their respective redox properties were evaluated using 2′,7′-dichlorodihydrofluorescein diacetate (H2-DCFDA) assay. Furthermore, their antiapoptotic potential was also evaluated by flow cytometry. The compound crystallizes in triclinic P-1 space group with unit cell parameters a = 5.7880 (4) Å, b = 9.8913 (6) Å, c = 14.5951 (9) Å, V = 1731.0 (3) Å3 and Z = 2. The crystal packing is stabilized by intermolecular hydrogen bonding, π···π, C-Br···π stacking interactions, and other non-covalent interactions. The mapping of different Hirshfeld surfaces and 2D-fingerprint were used to investigate intermolecular interactions. The interaction energies that stabilize the crystal packing were calculated and graphically represented as framework energy diagrams. We present a computational investigation of compound 4’s molecular structure at the Density Functional Theory level using the B3LYP method and the 6-31G ++ basis set in this paper. The optimized structure matches the experimental outcome. The global reactivity descriptors and molecular electrostatic potential (M.E.P.) map emphasize the molecule’s reactive locations, allowing reactivity prediction. The charge transfer properties of molecules can be estimated by examining Frontier molecular orbitals.
1 Introduction
Organoselenium compounds (OSe) have recently acquired significant interest as an exciting family of organic molecules with diverse applications in medicinal and organic chemistry (Phadnis, 2021; Chuai et al., 2021). They also possess enormous applications in advanced materials (Liao and Zhao, 2021). These special activities are due to the extraordinary properties of the selenium (Se) center (Handy et al., 2021). The non-metal bio-trace element Se is essential for the immune system’s normal function and protects cells from oxidative damage (Liao and Zhao, 2021; Makhal et al., 2021; Nogueira et al., 2021). It presents in most the living organisms as part of the selenoproteins and the antioxidant enzymes [e.g., thioredoxin reductases (TrxR) and glutathione peroxidase (GPx)] (Xu et al., 2020; Makhal et al., 2021; Radomska et al., 2021). Furthermore, the Se lower electronegativity (2.55), larger size (1.17 Å), and higher polarizability (3.8 Å) compared to sulfur (2.58, 1.02 Å, and 2.9 Å, respectively) made OSe compounds generally better nucleophiles (Sarma and Mugesh, 2005). Therefore, OSe compounds can react with O2-free radicals and thus attenuate oxidative stress-related disease progression (Rathore et al., 2019). Recently, OSe compounds have become promising candidates in cancer therapeutics (He et al., 2020). They have also manifested good histone deacetylase inhibitor activity (Adimulam et al., 2021). On the other hand, OSe compounds were extensively used as semiconductors in advanced materials, including photovoltaic cells and sodium-ion batteries and catalysts for the H2 evolution (Arora et al., 2021). Additionally, the pharmacological properties of OSe naturally occurring and drug molecules are attributed to the presence of the Se atom as a part of their scaffolds. Within this context, the selenocysteine (I), selenomethionine (II), and selenocystine (III) amino acids are present in the structure of several selenoproteins and selenoenzymes essential for the maintenance of metabolic rate, and immune responses, and oxidative homeostasis (Scheme 1) (Barbosa et al., 2017; Bartolini et al., 2017; Nogueira et al., 2021).
Furthermore, ebselen (IV) and ethaselen (VI) are among the most investigated OSe compounds with exciting GPx- and TrxR-like activities, respectively (Scheme 1) (Wang et al., 2012; Benelli et al., 2021). Indeed, they have recently entered clinical trial II as possible drugs for Meniere’s disease and non-small lung cancer, respectively (Wang et al., 2012; Shaaban et al., 2016a). We have recently directed our research toward synthesizing OSe-base maleimides (Shaaban et al., 2015a; Shaaban et al., 2015b; Shaaban et al., 2015c; Shaaban et al., 2018; Cherkaoui-Malki et al., 2019). The latter are privileged scaffolds found in many pharmacologically active drug molecules such as oxaleimide A (VI) and farinomalein (VII) (Viveki et al., 2021). Moreover, they are widely used in material science with industrial relevance, such as optoelectronic devices, resins, adhesives, rubber, and aerospace applications (Tonglairoum et al., 2016; Ravasco et al., 2019).
Toxilic anhydride reaction with amine is the common method used to prepare N-maleamic acids (Chen et al., 2007; Alam et al., 2022). Subsequent dehydration of the N-maleamic acids and ring closure affords the corresponding maleimides and isomaleimides, depending on the reaction conditions (Haval et al., 2006). Maleimides are the typical thermodynamically controlled product, whereas isomaleimides are usually kinetically controlled (Haval et al., 2006; Alam et al., 2022). Despite the numerous application and attention paid to maleimides, the synthesis of isomaleimides has gained little concern. Since the report of the 1st naphthyl-isomaleimide (Tsou et al., 1955) and new protocols are emerging describing the synthesis of such potential scaffolds; however, most of these evolving methods are limited to the dehydration of the respective maleamic acid with dehydrating agents under specified conditions. The dehydrating agents used include cyanuric chloride, oxalyl chloride, ethyl chloroformate, N,N′-dicyclohexylcarbodiimide, and 2-chloro-1,3-dimethylimidazolinium chloride, as well as triflouroacetic anhydride and propanephosphonic acid anhydride (Tsou et al., 1955; Burkitt and Gilbert, 1989; Haval et al., 2006; Guevara-Salazar et al., 2011; Haratake et al., 2011; Tonglairoum et al., 2016; Gao et al., 2019; Ravasco et al., 2019; Alam et al., 2022). Unfortunately, these protocols are limited by the lack of dehydrating agent with wide substrate scope and harsh conditions and long reaction time as well as the poor selectivity often associated with the formation of the undesired maleimide side product. Therefore, the development of an alternative, efficient, and general synthetic methods of isomaleimides remains a challenge.
We herein report the accidental synthesis and crystal structure of isomaleimide 4. Its cytoprotective activity was evaluated against oligodendrocytes. Furthermore, its antioxidant and antiapoptotic actives were also assessed. The Hirshfeld surface analysis description was used to characterize the nature of intermolecular interactions in crystal packing. Moreover, DFT simulations were used to optimize compound 4’s structure in its isolated condition. Furthermore, global reactivity descriptors and complementary interaction sites in isomaleimide 4 were identified using Frontier molecular orbitals (F.M.O.) and M.E.P. mapping studies.
2 Materials and methods
2.1 Synthesis of isomaleimide 4
Isomaleimide 4 was synthesized in good yield (77%) over three steps by the reduction of the diselenide 1 using NaBH4 and subsequent reaction with bromo-4-(bromomethyl)benzene followed by reaction with toxilic anhydride and subsequent ring closure using acetic anhydride and sodium acetate (Scheme 2) (Shaaban et al., 2018). Diselenide 1 was synthesized starting from aniline according to our reported method (Shaaban et al., 2018).
SCHEME 2. Reagents and conditions: (i) 4,4′-diselanediyldianiline (1) (2.5 mmol), bromo-4-(bromomethyl)benzene (3 mmol), EtOH (30 ml), NaBH4 (12.5 mmol); (ii) 4-((4-bromobenzyl)selanyl)aniline (2) (2.5 mmol), toluene (15 ml), toxilic anhydride (2.5 mmol), 4 h, r.t.; (iii) N-maleanilic acid 3 (2.5 mmol), Ac2O (8 ml), 250 mg NaOAc, 3 h, 40–50°C.
2.1.1 Synthesis of 4-((4-bromobenzyl)selanyl)aniline (2)
Compound 2 was synthesized from the reaction of diselenide 1 (344 mg, 1 mmol), 4-bromobenzyl bromide (550 mg, 2.2 mmol), aliquat 336 (45 mg, 5% mol) and sodium tetrahydridoborate (189.15 mg, 5 mmol) in ethylactetae and water mixed solvent (40 ml, 1:1) under reflux for 3 h. The formation of compound 2 was followed by TLC (petroleum ether: EtOAc = 6:1). The solvent was evaporated and the product was purified by column chromatography (petroleum ether: EtOAc = 6:1.5) as white solid; Yield: 310.31 mg (91%); mp 156–158°C (see supporting information for the analytical details).
2.1.2 Synthesis of 4-((4-((4-bromobenzyl)selanyl)phenyl)amino)-4-oxobut-2-enoic acid (3)
Compound 3 was synthesized from the reaction of compound 2 (344 mg, 1 mmol) and maleic anhydride (98 mg, 1 mmol) in dry toluene (5 ml) with stirring at r.t. for 3 h. The formation of compound 3 was monitored by TLC (chloroform: methanol = 8:1). The product was purified by column chromatography (chloroform: methanol = 6:1) to give yellow solid; Yield: 412.66 mg (94%); mp 215–217°C (see supporting information for the analytical details).
2.1.3 Synthesis of 4-((4-((4-Bromobenzyl)selanyl)phenyl)imino)-1H-pyrrol-2(5H)-one (4)
Compound 4 was synthesized from the gentle heating of compound 3 (439 mg, 1 mmol), acetic anhydride (3 ml), and sodium acetate (100 mg) for 2 h at 50–60°C. Water was added and the resulting mixture was extracted with CH2Cl2 (200 ml), dried with Mg2SO4, CH2Cl2 was evaporated and the residue was purified by silica gel chromatography. The progress of the product formation was followed by TLC petroleum ether: EtOAc = 8:1, Rf = 0.45, purified by column silica gel chromatography with petroleum ether: EtOAc = 6:1. Yellow solid; Yield: 324.17 mg (77%), mp 185–187°C. 1H NMR (300 MHz, CDCl3) δ 7.47–7.29 (m, 4H, Ar-H), 7.25–7.14 (m, 2H, Ar-H), 7.04–6.96 (m, 2H, Ar-H), 6.76 (d, J = 5.5 Hz, 1H, CH=), 6.62 (d, J = 5.5 Hz, 1H, CH=), 4.01 (s, 2H, SeCH2); 13C NMR (75 MHz, CDCl3) δ 166.94, 150.34, 143.16, 142.77, 137.50, 134.29, 134.02, 131.59, 130.49, 127.93, 125.94, 120.83, 31.52; MS (ESI): m/z = found 478.92 [M+ + CH3COO−]; calcd. 420.92 [M+]; HRMS calcd. for C17H12BrNO2Se [M+ + 1]: 475.93688, found 477.93289 [M+ + CH3COO−].
2.2 Biological evaluation
The antiproliferative activities of OSe compounds 2, 3, and 4 were assessed against two oligodendrocyte cell lines, namely, 158 N and 158 JP, by the M.T.T. assay using 7-ketocholesterol as reference control (Supporting info, Supplementary Table S1) (Vejux et al., 2005; Nury et al., 2014). The redox profile of the compounds was evaluated using the H2-DCFDA assay, and their antiapoptotic potential was also assessed by flow cytometry (Wilhelm et al., 2009; Karlsson et al., 2010; Shaaban et al., 2018).
2.3 Crystal structure measurement
A suitable crystal of compound 4 was mounted on a glass fiber loop. At 100 K, X-ray data were collected using a Bruker D8 VENTURE diffractometer and graphite-monochromated Mo(Kα) radiation (λ = 0.71073 Å). The full data set was used to calculate and refine unit cell parameters. S.A.D.A.B.S. Krause et al. (2015) was used to scale reflections and apply absorption corrections. The structure was solved using the ShelXT (Wilhelm et al., 2009; Karlsson et al., 2010) structure solution program using Direct Methods and refined using the ShelXL (Wilhelm et al., 2009) refinement package using Least Squares minimization using Olex2 (Krause et al., 2015). Hydrogen atoms were placed in calculated positions (C-H = 0.95–0.99 Å) and were constrained to ride on their parent atoms, with Uiso(H) = 1.2Ueq(C). The crystal packing and molecular packing were drawn using the program Diamond 3 (Sheldrick, 2008) and MERCURY (Macrae et al., 2020). The details of structural refinement and crystal data are reported in Table 1. Crystallographic data from the structural analysis have been deposited with the Cambridge Crystallographic Data Center, Nos CCDC-1472956 Copies of this information may be obtained free of charge from: The Director, C.C.D.C., 12 Union Road, Cambridge CB2 1EZ, United Kingdom. Fax: +44(1223)336-033, e-mail: deposit@ccdc.cam.ac.uk, or http://www.ccdc.cam.ac.uk.
2.4 Hirshfeld surface, 2D-fingerprint plot, and interaction energy calculations
Hirshfeld surfaces (Spackman and Jayatilaka, 2009) were mapped with property dnorm, and 2D-fingerprint plots were created with CrystalExplorer 17 (Spackman and Jayatilaka, 2009). Hirshfeld surfaces and 2D-fingerprint plots are excellent visualization tools for comparing intermolecular interactions in constructing various supramolecular motifs in the crystal structure. The 2D-fingerprint plot decomposes Hirshfeld surfaces into the contributions of different intermolecular interactions found in the crystal structure. The red, white, and blue colors on the Hirshfeld surface indicate shorter, equal, and longer contacts than the sum of the van der Waals radii, respectively. Furthermore, the crystal structure was analyzed using TONTO (Spackman and Jayatilaka, 2009; Spackman et al., 2021), with the CE-HF···HF/3-21G energy model (Sloot et al., 2003), starting with the .cif files derived from single-crystal X-ray diffraction data. An energy framework is a one-of-a-kind tool for visualizing crystal structures’ supramolecular architecture. To allow for comparison, total energy (Etot) is divided into electrostatic (Eele), polarization (Epol), dispersion (Edis), and repulsion (Erep) contribution energies, with cylinders representing the relative strength of the molecular packing and fixed at a scale factor of 200 and cutoff energy of 5 kJ/mol.
2.5 Theoretical calculations
The Gaussian 09 software (Shaaban et al., 2016b) was used for all calculations. The B3LYP functional (Becke’s three-parameter nonlocal exchange function with the Lee-Yang-Parr correlation function) (Hughes and Appel, 2019; Shaaban et al., 2021) was employed with the Density Functional Theory (DFT) approach. The geometrical optimizations were carried out in the gas phase and at the minima, and frequency calculations were used to confirm them; they accorded well with the experimental structure data. As a result, it was possible to compare energy and other physicochemical parameters with confidence. Chemical potential (μ), electronegativity (χ), electrophilicity index (ω) and chemical hardness (η), and softness (S) were all calculated using the corresponding HOMO and LUMO energies (Allen, 2002; Shaaban et al., 2019).
3 Results and discussion
3.1 Design and synthesis of isomaleimide 4
4,4′-Diselanediyldianiline (1) was used as the starting building block to prepare the target isomaleimide 4 (Shaaban et al., 2016b). In this regard, reduction of 1 with NaBH4 furnished the respective sodium phenylselenolate, which in turn trapped in ethanol and reacted with bromo-4-(bromomethyl)benzene to give Ose amine 2 in 85% yield (Scheme 2). The reaction of Ose amine 2 with toxilic anhydride in toluene at room temperature furnished the respective N-maleanilic acid 3 in excellent yield (94%). Surprisingly, gentle heating of compound 3 with Ac2O in the presence of NaOAc afforded the unexpected kinetically controlled isomaleimide 4 instead of thermodynamically controlled 1-(4-((4-bromobenzyl)selanyl)phenyl)-1H-pyrrole-2,5-dione (5) (Scheme 2) (Shaaban et al., 2018).
3.2 Charachterization of compound isomaleimide 4
The structure of isomaleimide 4 was established by 1HNMR and 13CNMR spectroscopy. 1HNMR of isomaleimide 4 showed the characteristic upfield singlet signal of the methylene fragment (CH2) at δ 4.01 ppm with a coupling constant (J) at 5.5 Hz. The characteristic isomaleimide vinylic protons appeared at δ 4.01 ppm δ 6.76 and 6.62 ppm. The rest eight aromatic protons appeared as multiplet signals at δ 7.47, 7.25, and 7.04 ppm. Furthermore, the isomaleimide 4 showed 13 carbon signals in the 13CNMR spectroscopy. The aliphatic methylene carbon appeared upfield at δ 31.52 ppm. On the other hand, one carbonyl signal appeared downfield at δ 166.94 ppm and the azomethine (N=C) at δ 150.34 ppm. The vinylic carbons appeared at δ 137.50 and δ 127.93 ppm. The rest eight aromatic carbons appeared at δ 143.16–120.83 ppm.
3.3 Biology
Oligodendrocytes are susceptible to deterioration by oxygen and nitrogen reactive species (R.O.S. and R.N.S.), affecting the transmission of the neuronal signal and the proper axon function (Shaaban et al., 2018; Hughes and Appel, 2019; Shaaban et al., 2021). OSe compounds have recently manifested potential chemoprotective and antioxidant properties (Shaaban et al., 2015a; Shaaban et al., 2016b; Shaaban et al., 2019). The cytoprotective properties of OSe compounds 2, 3, and 4 were assessed in 158N and 158 JP cells. Interestingly, OSe compounds 2 and 4 did not show any apparent toxicity (IC50 ≥ 100 µM), whereas the N-maleanilic acid 3 showed moderate-low cytotoxicity (Supporting info, Supplementary Table S1). These results encourage the further assessment of their redox properties. Considering that OSe compounds are good nucleophilic reductants, they can react with R.O.S. and R.N.S. to protect tissues from oxidative damage (OD) (Sarma and Mugesh, 2005; Chuai et al., 2021; Nogueira et al., 2021). Oligodendrocytes are usually vulnerable to OD; thus, our main objective is further to estimate the antioxidant potential of compounds 2, 3, and 4 using the H2-DCFDA (Supporting info, Supplementary Table S1). Interestingly, all compounds significantly diminished the R.O.S. levels in 158N and were even more significant than vitamin E.
Furthermore, the antiapoptotic activities of compounds 2, 3, and 4 were also evaluated by flow cytometry in 158N cells stained with P.I. using 7Kc as apoptosis stimulator and the positive reference (Supporting info, Supplementary Figure S1). Among the tested compounds, the N-maleanilic acid 3 lowered, in a concentration-dependent manner, the formation of SubG1 peak. To this point, N-maleanilic acid 3 showed antioxidant and antiapoptotic activities probably via diminishing the R.O.S. levels. The N-maleanilic acid 3 amphiphilic characters favor its crossing through the cell membrane and, therefore, its better activity.
3.4 Analysis of the molecular packing
The isomaleimide 4 crystallized in the triclinic with space group P-1. Its molecular structure is shown in Figure 1. The main body of the (A) structure consists of a Furan ring (Cg1) which has an α,β-unsaturated carbonyl group, and two benzene rings (Cg2 and Cg3). The planarity of the Cg1 (C14-C17-O1/O2), Cg2 (C2-C7/Br1/C1), and Cg3 (C8-C13/Se1) is evident by the root-mean-square (rms) deviation of 0.0096 of Cg1, 0.0065 of Cg2, and 0.0084 Å of Cg3. The dihedral angle between Cg2 and Cg3 is 72.32 (9)°. The central Se atom shows a bent geometry [C1-Se1-C8 = 98.29 (10)°]. All bond distances and angles in isomaleimide 4 are within the acceptable ranges (Table 2) (Allen, 2002; Gouda et al., 2022; Shaaban et al., 2022).
The molecules in isomaleimide 4 are linked together by weak non-classical Hydrogen-bonding chains (C16-H16···O1iii and C15-H15···Nii [symmetry codes: (ii) −x−1, −y+1, −z; (iii) −x, −y+2, −z]), (Figure 2; Table 3). It is worth noting that the proximity of donor atoms C12 and acceptor atom O2 allows for the formation of one intramolecular hydrogen bond C12-H12···O2. Bromine atoms play an essential role in crystal packing stabilization because they are involved in C12-H12···Bri hydrogen bonding [symmetry code: (i) −x+2, −y+1, −z+1], (Figure 2). The main characteristics of the hydrogen bonds data are listed in Table 3. Slipped π-stacking interactions strengthen the chains of molecules along a axis. The centroid-to-centroid separation of Cg1 and Cg3 is 3.564(1) Å (Figure 3; Table 3). The crystal packing of 4 is also dominated by C-Br···Cg1 (dBr···Cg1 = 3.7321(11) Å), Se···Br and Se···O interactions (Figure 4; Table 3). The closest intermolecular Se···Br distance is 3.8251(5) Å, which is greater than the sum of the van der Waals radii for these two elements (3.75 Å) (Bondi, 1964).
FIGURE 2. Isomaleimide 4 packing diagram, viewed along a axis, showing the hydrogen-bonded chains structure (dashed line). C-H···O (Red), C-H···N (Cyan), and C-H···Br (Green).
FIGURE 4. (a) C-Br···Cg (Green dashed lines), Se···O (Purple dashed lines), and Se···Br (Pink dashed lines) interactions in isomaleimide 4.
3.5 Hirshfeld surfaces, fingerprint plots, and energy framework
The Hirshfeld surface of isomaleimide 4 is mapped over dnorm (Figure 5A) in the ranges −0.2503 to 1.4928 Å, shape index (Figure 5B) in the range (−0.9976 to 0.9980 Å), curvedness (Figure 5C) in the range −4.1687 to 0.2406 Å and fragment patches (Figure 5D) in the range 0.000–14.000 Å.
FIGURE 5. (A) dnorm mapped on Hirshfeld surfaces for visualizing the intermolecular interactions, (B) shape-index, (C) curvedness, and (D) fragment patches of isomaleimide 4.
Around the C.H. and the carbonyl group in the Furan ring and the N-atom of isomaleimide 4, the H.S. has intense red spots, which indicates that these specific atoms are promising in the H-bonding interactions (C15-H15···N and C16-H16···O1). By plotting H.S., we can not only visualize H-bonding interactions but also investigate other non-covalent interaction types, such as stacking interactions. To visualize this interaction, the H.S. is plotted on a shape index and curvedness. The presence of triangular regions of red and blue around the furan and aromatic rings on the shape index (Figure 5B) and a flat region around the phenyl rings and furan rings in curvedness (Figure 5C) indicates the presence of π-stacking interactions. The nearest molecules’ neighbor environment is determined by the color patches (Figure 5D) on the Hirshfeld surface and their proximity to adjacent molecules.
Figure 6 shows the two-dimensional fingerprint of the significant contacts contributing to the Hirshfeld surface of isomaleimide 4 is showed in expanded mode. Interatomic contact percentage contributions in isomaleimide 4 are calculated using two-dimensional fingerprint plots. The C∙∙∙H/H∙∙∙C contacts contribute the most (26.2%), corresponding to C-H···N and C-H···Br interactions, as shown in the 2D fingerprint plot by a pair of pointed spikes representative of strong hydrogen-bonding interactions (Figure 6A).
FIGURE 6. Fingerprint plots with fragment patches of (a) C∙∙∙H, (b) H∙∙∙H, (c) O∙∙∙H and (d) H∙∙∙Br contacts in isomaleimide 4.
The H∙∙∙H contact appears in the middle of the scattered points in the two-dimensional fingerprint plots (Figure 6B) with a contribution to the overall Hirshfeld surface of 25.4%. Because of the abundance of hydrogen on the molecular surface, they are the second most common interactions (55%). O∙∙∙H/H∙∙∙O contacts make up the third-largest contribution to the Hirshfeld surface (Figure 6C), accounting for 13%. This contact indicates the presence of intermolecular C-H···O hydrogen bonds. The Br∙∙∙H/H∙∙∙Br contacts (Figure 6D), which raise to N-H∙∙∙Br interactions, are the third most crucial interaction on the surface, accounting for approximately 9.2% of the Hirshfeld surfaces. The relative percentage contributions to the overall Hirshfeld surface are shown in Figure 7. Finally, the Hirshfeld surface analysis yields the same results as the X-ray crystal structure analysis and provides a new visual explanation for intermolecular interactions.
FIGURE 7. Different intermolecular contacts’ relative contributions to the Hirshfeld surface area in isomaleimide 4.
The intermolecular interaction energies are calculated using the HF/3-21G energy model with scale factors to determine Etot: kele = 1.019, kpol = 0.651, kdis = 0.901, krep = 0.811, where a cluster of molecules is generated by applying crystallographic symmetry operations with respect to a chosen central molecule within a radius of 3.8 by default. The total interaction energies (Etot = −263 kJ/mol) can be divided into total electrostatic (Eele = −54.7 kJ/mol), total polarization (Epol = −33.6 kJ/mol), total dispersion (Edis = −327 kJ/mol), and total repulsion (Erep = 137 kJ/mol) (Table 4). Obviously, the interactions between neighboring molecules contribute the most to the stability of this structure, and in most cases, the dispersion energy component dominates these interactions. Additionally, even for the hydrogen-bonded molecule pairs, the dispersion energy component is dominant or comparable to the electrostatic components which due to the presence of π-π and C-Br···π interactions (the highest importance of electrostatic components is observed in a pair having Eele = −22.6 kJ mol−1, Epol = −4.5 kJ mol−1 and Edisp = −54.8 kJ mol−1). Figure 8 depicts the Coulomb interaction energy (red), dispersion energy (green), and total interaction energy (blue) between molecular pairs along the a axis with respect to the selected molecule.
FIGURE 8. Energy-framework diagrams of isomaleimide 4 along the a axis: Electrostatic interaction energies (red); dispersion interaction energies (green); total interaction energies (blue).
3.6 Density functional theory calculations
3.6.1 Geometric structures
The optimized structural parameters of the isomaleimide 4 are depicted in Figure 9. In this theoretical study, the hybrid functional B3LYP with the basis 6-311 ++ G(d,p) has been used in all the calculations made by the Gaussian 09 program. The optimization of the structure of the studied compound has been carried out, starting from the geometry of X-rays. A comparison of experimental results with theoretical ones (Table 5) reveals that most of the calculated values of bond lengths and angles are very close to those experimental data, which shows that the choice of the base 6-311 ++ G(d,p) is suitable for this theoretical study. However, the slight difference observed can be attributed to the environment of the molecule studied, being isolated in phase gaseous for the theoretical study and subjected to interactions with solid-state intermolecular molecules in the experimental study. Figure 10 depicts the visual HOMO and LUMO of compound 4.
Furthermore, we believe it is worthwhile, to begin with the results gained from the geometrical parameters of the investigated molecules. In Table 4, the results of this calculation are grouped by the following numbering scheme. Compare the experimental results from a crystallographic investigation with these theoretically computed geometric characteristics (Klocker et al., 2003). The relation:
Where, the theoretical value of the quantity X is Xtheo, while the experimental value is Xexp. The average variance of the distances and angles produced by the DFT approach is on the order of less than 2%, indicating that the diverse theoretical conclusions obtained are in excellent agreement with those found experimentally by crystallographic analysis.
3.6.2 Study of overall global reactivities
The amount of energy required to remove an electron from a molecule is known as ionization potential. Furthermore, high ionization energy denotes great stability and thus chemical inertness, whereas low ionization energy shows a molecule’s reactivity. The energy generated when an electron is introduced to a neutral molecule is characterized as electron affinity isomaleimide 4. A significant value isomaleimide 4 shows the molecule’s tendency to keep its electrons. A negative chemical potential (μ) reflects the molecule’s molecular stability or difficulty breaking it down into its constituent parts (Table 6). The resistance of the cloud of molecular electrons to deformation during tiny perturbations is measured in hardness (η). A big HOMO-LUMO energy gap implies a complicated molecule with low polarizability and chemical and biological activities but high kinetic sensitivity (Table 6).
In contrast, a small HOMO-LUMO energy gap indicates a soft molecule with high polarizability and activities but low kinetic sensitivity. Chemical and biological sensitivities are low, whereas kinetic sensitivity is heightened. The overall electrophilicity index (ω) of a molecule measures its stabilization energy or resistance to exchange electrons with the system after the addition of an external electronic charge (Parr et al., 1999).
The gap energy, which is the energy difference of the two preceding molecular orbitals (Eg = E.L.U.M.O.–E.H.O.M.O.) is 2.90 eV. Chemical potential, electronegativity, hardness, softness, Global Softness, and electrophilicity for the isomaleimide 4 were determined as −4.52, 4.52, 1.45, 0.73, 0.69, and 7.04 eV, respectively, using Eqs 2–7 (Rajan and Muraleedharan, 2017):
3.6.3 Surfaces with molecular electrostatic potential
Isomaleimide 4’s M.E.P. was computed using the DFT-B3LYP/6-31G ++ optimized geometry, and its surface map is shown in Figure 11. This diagram uses a color scheme to indicate the electrostatic potential values. The highest negative value is represented by the red color, indicating the most likely areas for electrophilic assault. The most positively charged areas appear in dark blue color, indicating potential nucleophilic attack sites. The determined limits are −4.64e-2 (deepest red) and +4.64e-2 (deepest blue), with the intermediate color scale flowing from red through orange, yellow, green, and blue in that order, as illustrated in Figure 11. The most substantial negative potential is centered near the oxygen atoms, while positive potentials are scattered throughout the TTF unit, particularly around the outer H atoms. Finally, the M.E.P. primarily suggests an electrophilic attack on oxygen atoms, with the possibility of a nitrogen atom attack. On the other hand, a strong base may be able to destroy one of the H atoms as a proton.
FIGURE 11. Molecular electrostatic potential (M.E.P.) map of isomaleimide 4 calculated at the 6-311++G(d,p) level.
4 Conclusion
Isomaleimide 4 was accidentally obtained in 77% yield instead of the maleimide 5 via dehydration and cyclization of the respective N-maleanilic acid 3 upon heating with acetic anhydride. The molecular structure of isomaleimide 4 was confirmed by X-ray diffraction analysis. The cytotoxicity was assessed against two oligodendrocytes and the antioxidant properties were evaluated using H2-DCFDA assay. The intermolecular interactions in the crystal packing are quantified and visualized using Hirshfeld surfaces, 2D fingerprint plots, and 3D energy frameworks. The two-dimensional fingerprint revealed that the most significant contributions to these surfaces come from C∙∙∙H/H∙∙∙C (26.2%), H···H (25.4%), O∙∙∙H/H∙∙∙O (13%) and Br∙∙∙H/H∙∙∙Br (9.2%) interactions. The energy-framework analysis reveals that the dispersive energies are the most important forces in the crystal. In parallel with the experimental study, we carried out a theoretical study detail using quantum chemical methods to determine the properties structural of compound 4. We carried out a comparison between the theoretical geometrical parameters and those obtained by X-ray diffraction. It appears significantly that the calculations obtained are in good agreement with the experimental data.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
Conceptualization, SS, TY, MA, and HF; methodology, SS, TY, and HF; software, TY and HF; validation, MG, HA, and MK; formal analysis, SS, TY, and HF; investigation, SS, TY, and HF; resources, MG, HA, and MK; data curation, SS, TY, and HF; writing—original draft preparation, SS, TY, and HF; writing—review and editing, SS, TY, MA, and HF; visualization, MG, HA, and MK; supervision, SS, TY, MA, and HF; project administration, SS; funding acquisition, MG, HA, and MK. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research at King Faisal University, Saudi Arabia, for financial support under the annual funding track (Grant No. 810). Also, the authors extend their appreciation to the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University for funding this work through Research Group No. RG-21-09-72.
Acknowledgments
The authors acknowledge the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research at King Faisal University, Saudi Arabia, for financial support under annual funding track (GRANT810). Also, the authors extend their appreciation to the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University for funding this work through Research Group No. RG-21-09-72.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2022.961787/full#supplementary-material
References
Adimulam, T., Arumugam, T., Foolchand, A., Ghazi, T., and Chuturgoon, A. A. (2021). The effect of organoselenium compounds on histone deacetylase inhibition and their potential for cancer therapy. Int. J. Mol. Sci. 22, 12952. doi:10.3390/ijms222312952
Alam, K., McFee, E. C., and Croatt, M. P. (2022). Rapid and facile synthesis of isomaleimides: Dehydration of maleamic acids using methanesulfonyl chloride. Synthesis 54 (13), 3085–3092. doi:10.1055/s-0041-1737414
Allen, F. H. (2002). The Cambridge structural database: a quarter of a million crystal structures and rising. Acta Crystallogr. B 58, 380–388. doi:10.1107/s0108768102003890
Arora, A., Oswal, P., Kumar Rao, G., Kumar, S., and Kumar, A. (2021). Organoselenium ligands for heterogeneous and nanocatalytic systems: development and applications. Dalton Trans. 50, 8628–8656. doi:10.1039/d1dt00082a
Barbosa, N. V., Nogueira, C. W., Nogara, P. A., de Bem, A. F., Aschner, M., Rocha, J. B. T., et al. (2017). Organoselenium compounds as mimics of selenoproteins and thiol modifier agents. Metallomics 9, 1703–1734. doi:10.1039/c7mt00083a
Bartolini, D., Sancineto, L., Fabro de Bem, A., Tew, K. D., Santi, C., Radi, R., et al. (2017). Selenocompounds in cancer therapy: An overview. Adv. Cancer Res. 136, 259–302. doi:10.1016/bs.acr.2017.07.007
Benelli, J. L., Poester, V. R., Munhoz, L. S., Melo, A. M., Trapaga, M. R., Stevens, D. A., et al. (2021). Ebselen and diphenyl diselenide against fungal pathogens: A systematic review. Med. Mycol. Open. Access. 59, 409–421. doi:10.1093/mmy/myaa115
Bondi, A. (1964). van der Waals volumes and radii. J. Phys. Chem. 68, 441–451. doi:10.1021/j100785a001
Burkitt, M. J., and Gilbert, B. C. (1989). The control of iron-induced oxidative damage in isolated rat-liver mitochondria by respiration state and ascorbate. Free Radic. Res. Commun. 5, 333–344. doi:10.3109/10715768909073416
Chen, X., Zheng, Y., and Shen, Y. (2007). Natural products with maleic anhydride structure: nonadrides, tautomycin, chaetomellic anhydride, and other compounds. Chem. Rev. 107, 1777–1830. doi:10.1021/cr050029r
Cherkaoui-Malki, M., Shaaban, S., Tahri-Joutey, M., Elshobaky, A., Saih, F.-E., Vervandier-Fasseur, D., et al. (2019). Cytoprotective and antioxidants in peroxisomal neurodegenerative diseases. Multidiscip. Digit. Publ. Inst. Proc. 11, 33.
Chuai, H., Zhang, S. Q., Bai, H., Li, J., Wang, Y., Sun, J., et al. (2021). Small molecule selenium-containing compounds: Recent development and therapeutic applications. Eur. J. Med. Chem. 223, 113621. doi:10.1016/j.ejmech.2021.113621
Gao, X., Tang, L., Huang, L., Huang, Z. S., Ma, Y., Wu, G., et al. (2019). Oxidative aminoarylselenation of maleimides via copper-catalyzed four-component cross-coupling. Org. Lett. 21, 745–748. doi:10.1021/acs.orglett.8b03980
Gouda, M., Ferjani, H., Abd El-Lateef, H. M., Khalaf, M. M., Shaaban, S., Yousef, T. A., et al. (2022). A competition between hydrogen, stacking, and halogen bonding in N-(4-((3-Methyl-1, 4-dioxo-1, 4-dihydronaphthalen-2-yl) selanyl) phenyl) acetamide: Structure, Hirshfeld surface analysis, 3D energy framework approach, and DFT calculation. Int. J. Mol. Sci. 23, 2716. doi:10.3390/ijms23052716
Guevara-Salazar, J. A., Quintana-Zavala, D., Jiménez-Vázquez, H. A., and Trujillo-Ferrara, J. (2011). Synthesis of Diels–Alder adducts of N-arylmaleimides by a multicomponent reaction between maleic anhydride, dienes, and anilines. Monatsh. Chem. 142, 827–836. doi:10.1007/s00706-011-0515-5
Handy, D. E., Joseph, J., and Loscalzo, J. (2021). Selenium, a micronutrient that modulates cardiovascular health via redox enzymology. Nutrients 13, 3238. doi:10.3390/nu13093238
Haratake, M., Sakano, T., Fuchigami, T., and Nakayama, M. (2011). Thiol-targeted introduction of selenocysteine to polypeptides for synthesis of glutathione peroxidase mimics. Metallomics 3, 702. doi:10.1039/c1mt00001b
Haval, K. P., Mhaske, S. B., and Argade, N. P. (2006). Cyanuric chloride: decent dehydrating agent for an exclusive and efficient synthesis of kinetically controlled isomaleimides. Tetrahedron 62, 937–942. doi:10.1016/j.tet.2005.10.027
He, X., Zhong, M., Li, S., Li, X., Li, Y., Li, Z., et al. (2020). Synthesis and biological evaluation of organoselenium (NSAIDs-SeCN and SeCF3) derivatives as potential anticancer agents. Eur. J. Med. Chem. 208, 112864. doi:10.1016/j.ejmech.2020.112864
Hughes, A. N., and Appel, B. (2019). Oligodendrocytes express synaptic proteins that modulate myelin sheath formation. Nat. Commun. 10, 4125. doi:10.1038/s41467-019-12059-y
Karlsson, M., Kurz, T., Brunk, U. T., Nilsson, S. E., and Frennesson, C. I. (2010). What does the commonly used DCF test for oxidative stress really show? Biochem. J. 428, 183–190. doi:10.1042/bj20100208
Klocker, J., Karpfen, A., and Wolschann, P. (2003). Surprisingly regular structure–property relationships between C–O bond distances and methoxy group torsional potentials: An ab initio and density functional study. J. Mol. Struct: THEOCHEM. 635, 141–150.
Krause, L., Herbst-Irmer, R., Sheldrick, G. M., and Stalke, D. (2015). Comparison of silver and molybdenum microfocus X-ray sources for single-crystal structure determination. J. Appl. Crystallogr. 48, 3–10. doi:10.1107/s1600576714022985
Liao, L., and Zhao, X. (2021). Modern organoselenium catalysis: Opportunities and challenges. Synlett 32 (13), 1262–1268. doi:10.1055/a-1506-5532
Macrae, C. F., Sovago, I., Cottrell, S. J., Galek, P. T., McCabe, P., Pidcock, E., et al. (2020). Mercury 4.0: From visualization to analysis, design and prediction. J. Appl. Crystallogr. 53, 226–235. doi:10.1107/s1600576719014092
Makhal, P. N., Nandi, A., and Kaki, V. R. (2021). Insights into the recent synthetic advances of organoselenium compounds. ChemistrySelect 6, 663–679. doi:10.1002/slct.202004029
Nogueira, C. W., Barbosa, N. V., and Rocha, J. B. T. (2021). Toxicology and pharmacology of synthetic organoselenium compounds: an update. Arch. Toxicol. 95, 1179–1226. doi:10.1007/s00204-021-03003-5
Nury, T., Zarrouk, A., Vejux, A., Doria, M., Riedinger, J. M., Delage-Mourroux, R., et al. (2014). Induction of oxiapoptophagy, a mixed mode of cell death associated with oxidative stress, apoptosis and autophagy, on 7-ketocholesterol-treated 158N murine oligodendrocytes: impairment by alpha-tocopherol. Biochem. Biophys. Res. Commun. 446, 714–719. doi:10.1016/j.bbrc.2013.11.081
Parr, R. G., Szentpály, L. V., and Liu, S. (1999). Electrophilicity index. J. Am. Chem. Soc. 121, 1922–1924.
Phadnis, P. P. (2021). “Synthesis strategies for organoselenium compounds and their potential applications in human life,” in Handbook on synthesis strategies for advanced materials (Germany: Springer), 537–641.
Radomska, D., Czarnomysy, R., Radomski, D., and Bielawski, K. (2021). Selenium compounds as novel potential anticancer agents. Int. J. Mol. Sci. 22, 1009. doi:10.3390/ijms22031009
Rajan, V. K., and Muraleedharan, K. (2017). A computational investigation on the structure, global parameters and antioxidant capacity of a polyphenol, gallic acid. Food Chem. 220, 93–99.
Rathore, V., Jose, C., and Kumar, S. (2019). Organoselenium small molecules as catalysts for the oxidative functionalization of organic molecules. New J. Chem. 43, 8852–8864. doi:10.1039/c9nj00964g
Ravasco, J. M., Faustino, H., Trindade, A., and Gois, P. M. (2019). Bioconjugation with maleimides: A useful tool for chemical biology. Chem. Eur. J. 25, 43–59. doi:10.1002/chem.201803174
Sarma, B. K., and Mugesh, G. (2005). Glutathione peroxidase (GPx)-like antioxidant activity of the organoselenium drug ebselen: unexpected complications with thiol exchange reactions. J. Am. Chem. Soc. 127, 11477–11485. doi:10.1021/ja052794t
Shaaban, S., Gaffer, E. H., Alshahd, M., and Elmorsy, S. S. (2015a). Cytotoxic symmetrical thiazolediselenides with increased selectivity against MCF-7 breast cancer cells. Int. J. Res. Dev. Pharm. Life Sci. 4, 1654–1668.
Shaaban, S., Gaffer, E. H., Jabar, Y., and Elmorsy, S. S. (2015b). Cytotoxic naphthalene based-symmetrical diselenides with increased selectivity against MCF-7 breast cancer cells. Int. J. Pharm. 5, 738–746.
Shaaban, S., Negm, A., Sobh, M. A., and Wessjohann, L. A. (2015c). Organoselenocyanates and symmetrical diselenides redox modulators: Design, synthesis and biological evaluation. Eur. J. Med. Chem. 97, 190–201. doi:10.1016/j.ejmech.2015.05.002
Shaaban, S., Negm, A., Ashmawy, A. M., Ahmed, D. M., and Wessjohann, L. A. (2016a). Combinatorial synthesis, in silico, molecular and biochemical studies of tetrazole-derived organic selenides with increased selectivity against hepatocellular carcinoma. Eur. J. Med. Chem. 122, 55–71. doi:10.1016/j.ejmech.2016.06.005
Shaaban, S., Negm, A., Sobh, M. A., and Wessjohann, L. A. (2016b). Expeditious entry to functionalized pseudo-peptidic organoselenide redox modulators via sequential ugi/SN methodology. Anticancer. Agents Med. Chem. 16, 621–632. doi:10.2174/1871520615666150916092035
Shaaban, S., Vervandier-Fasseur, D., Andreoletti, P., Zarrouk, A., Richard, P., Negm, A., et al. (2018). Cytoprotective and antioxidant properties of organic selenides for the myelin-forming cells, oligodendrocytes. Bioorg. Chem. 80, 43–56. doi:10.1016/j.bioorg.2018.05.019
Shaaban, S., Ashmawy, A. M., Negm, A., and Wessjohann, L. A. (2019). Synthesis and biochemical studies of novel organic selenides with increased selectivity for hepatocellular carcinoma and breast adenocarcinoma. Eur. J. Med. Chem. 179, 515–526. doi:10.1016/j.ejmech.2019.06.075
Shaaban, S., Zarrouk, A., Vervandier-Fasseur, D., Al-Faiyz, Y. S., El-Sawy, H., Althagafi, I., et al. (2021). Cytoprotective organoselenium compounds for oligodendrocytes. Arabian J. Chem. 14, 103051. doi:10.1016/j.arabjc.2021.103051
Shaaban, S., Ferjani, H., Yousef, T., and Abdel-Motaal, M. (2022). Supramolecular self-assembly built by hydrogen, stacking and Br··· Br interactions in 4-((4-bromobenzyl) selanyl) aniline: Structure, Hirshfeld surface analysis, 3D energy framework approach and global reactivity descriptors. J. Inorg. Organomet. Polym. Mat. 32, 1878–1890. doi:10.1007/s10904-022-02284-2
Sheldrick, G. M. (2008). A short history of SHELX. Acta Crystallogr. A 64, 112–122. doi:10.1107/S0108767307043930
Sheldrick, G. M. (2015). SHELXT - integrated space-group and crystal-structure determination. Acta Crystallogr. A Found. Adv. 71, 3–8. doi:10.1107/S2053273314026370
Sloot, P. M., Abramson, D., Bogdanov, A. V., Dongarra, J. J., Zomaya, A. Y., and Gorbachev, Y. E. (2003). Computational science—ICCS 2003: International conference melbourne, Australia and st. Petersburg, Russia june 2–4, 2003 proceedings, Part I. Germany: Springer, 2657.
Spackman, M. A., and Jayatilaka, D. (2009). Hirshfeld surface analysis. CrystEngComm 11, 19–32. doi:10.1039/b818330a
Spackman, P. R., Wolff, S. K., Grimwood, D. J., Jayatilaka, D., Spackman, M. A., Turner, M. J., et al. (2021). CrystalExplorer: a program for Hirshfeld surface analysis, visualization and quantitative analysis of molecular crystals. J. Appl. Crystallogr. 54, 1006–1011. doi:10.1107/s1600576721002910
Tonglairoum, P., Brannigan, R. P., Opanasopit, P., and Khutoryanskiy, V. V. (2016). Maleimide-bearing nanogels as novel mucoadhesive materials for drug delivery. J. Mat. Chem. B 4, 6581–6587. doi:10.1039/c6tb02124g
Tsou, K.-C., Barrnett, R. J., and Seligman, A. M. (1955). Preparation of some N-(1-Naphthyl)-maleimides as sulfhydryl group Reagents1a, 1b. J. Am. Chem. Soc. 77, 4613–4616. doi:10.1021/ja01622a051
Vejux, A., Kahn, E., Dumas, D., Bessede, G., Menetrier, F., Athias, A., et al. (2005). 7-Ketocholesterol favors lipid accumulation and colocalizes with Nile Red positive cytoplasmic structures formed during 7-ketocholesterol-induced apoptosis: Analysis by flow cytometry, FRET biphoton spectral imaging microscopy, and subcellular fractionati. Cytom. A 64, 87–100. doi:10.1002/cyto.a.20124
Viveki, A. B., Pol, M. D., Halder, P., Sonavane, S. R., and Mhaske, S. B. (2021). Annulation of enals with carbamoylpropiolates via NHC-catalyzed enolate pathway: Access to functionalized maleimides/iso-maleimides and synthesis of Aspergillus FH-X-213. J. Org. Chem. 86, 9466–9477. doi:10.1021/acs.joc.1c00782
Wang, L., Yang, Z., Fu, J., Yin, H., Xiong, K., Tan, Q., et al. (2012). Ethaselen: a potent mammalian thioredoxin reductase 1 inhibitor and novel organoselenium anticancer agent. Free Radic. Biol. Med. 52, 898–908. doi:10.1016/j.freeradbiomed.2011.11.034
Wilhelm, J., Vytasek, R., Ostadalova, I., and Vajner, L. (2009). Evaluation of different methods detecting intracellular generation of free radicals. Mol. Cell. Biochem. 328, 167–176. doi:10.1007/s11010-009-0086-5
Keywords: isomaleimide, organoselenium, crystal structure, antioxidant, hirshfeld surface analysis, DFT calculations, 3D energy framework
Citation: Shaaban S, Ferjani H, Abd El-Lateef HM, Khalaf MM, Gouda M, Alaasar M and Yousef TA (2022) Unexpected kinetically controlled organoselenium-based isomaleimide: X-ray structure, hirshfeld surface analysis, 3D energy framework approach, and density functional theory calculation. Front. Chem. 10:961787. doi: 10.3389/fchem.2022.961787
Received: 05 June 2022; Accepted: 30 June 2022;
Published: 05 August 2022.
Edited by:
Prasat Kittakoop, Chulabhorn Graduate Institute, ThailandCopyright © 2022 Shaaban, Ferjani, Abd El-Lateef, Khalaf, Gouda, Alaasar and Yousef. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Saad Shaaban, sibrahim@kfu.edu.sa, dr_saad_chem@mans.edu.eg; Mohamed Alaasar, mohamed.alaasar@chemie.uni-halle.de; Tarek A. Yousef, tayousef@imamu.edu.sa