key: cord-0848758-enr2q0xw
authors: Barakat, Assem; Haukka, Matti; Soliman, Saied M.; Ali, M.; Al-Majid, Abdullah Mohammed; El-Faham, Ayman; Domingo, Luis R.
title: Straightforward Regio- and Diastereoselective Synthesis, Molecular Structure, Intermolecular Interactions and Mechanistic Study of Spirooxindole-Engrafted Rhodanine Analogs
date: 2021-11-30
journal: Molecules
DOI: 10.3390/molecules26237276
sha: 6fd9344b43622f40d0aefd2a8c80508a5294490d
doc_id: 848758
cord_uid: enr2q0xw

Straightforward regio- and diastereoselective synthesis of bi-spirooxindole-engrafted rhodanine analogs 5a–d were achieved by one-pot multicomponent [3 + 2] cycloaddition (32CA) reaction of stabilized azomethine ylide (AYs 3a–d) generated in situ by condensation of L-thioproline and 6-chloro-isatin with (E)-2-(5-(4-chlorobenzylidene)-2,4-dioxothiazolidin-3-yl)-N-(2-morpholinoethyl)acetamide. The bi-spirooxindole-engrafted rhodanine analogs were constructed with excellent diastereo- and regioselectivity along with high chemical yield. X-ray crystallographic investigations for hybrid 5a revealed the presence of four contiguous stereocenters related to C11, C12, C19 and C22 of the spiro structure. Hirshfeld calculations indicated the presence of many short intermolecular contacts such as Cl...C, S...S, S...H, O...H, N...H, H...C, C...C and H...H interactions. These contacts played a very important role in the crystal stability. The polar nature of the 32CA reaction was studied by analysis of the conceptual DFT reactivity indices. Theoretical study of this 32CA reaction indicated that it takes place through a non-concerted two-stage one-step mechanism associated with the nucleophilic attack of AY 3a to the electrophilic ethylene derivative.

Rhodanines (2-thioxothiazolidin-4-ones) privileged structure has attracted much attention in pharmaceutical and medicinal chemistry [1, 2] . These substituted rhodanines have been discovered and reported for several pharmaceutical applications, including topoisomerase II inhibition potency [3] ; potential AChE inhibitors [4, 5] ; HIV-1 integrase inhibitors [6] ; anti-cancer agent towards MCF-7 breast cancer cells [7] and antidiabetic (T2DM) targeting the inhibition of enzymes such as aldose reductase (ALR) [8] , α-glucosidase [9] , α-amylase [10] and PTP1B enzyme [11] , and other therapeutic targets have also been reported, such as antibacterial ones [12] . Chemical modification, structural elucidation and molecular interactions investigations of these compounds are still of high interest.

On the other hand, one-pot multicomponent [3 + 2] cycloaddition (32CA) reaction has been established as a powerful and efficient approach for stereoselective synthesis of the spirooxindole pyrrolidine core structure in a wide range of natural products and pharmaceutical compounds [13] [14] [15] [16] [17] [18] [19] . The key step of this strategy is the generation of azomethine

The straightforward regio-and diastereoselective synthesis of the desired bi-spirooxindole-engrafted rhodanine analogs 5a-d is shown in Scheme 1. The requisite starting material named Scheme 1. Synthesis of compounds 5a-d via 32CA reaction of AYs 3a-d with ethylene derivative 4. 

The straightforward regio-and diastereoselective synthesis of the desired bi-spirooxind ole-engrafted rhodanine analogs 5a-d is shown in Scheme 1. The requisite starting material named (E)-2-(5-(4-chlorobenzylidene)-2,4-dioxothiazolidin-3-yl)-N-(2-morpholinoethyl)ac etamide 4 was synthesized following the method reported in the literature [29] . Employing the one-pot multicomponent-based 32CA reaction approach of the arylidene rhodanine analog 4 with di-carbonyl compounds (6-chloroisatin 2a; 5-nitroisatin 2b; 5-fluoroisatin 2c; 5-chloroisatin 2d) and thioproline 1 under reflux conditions in MeOH for 2 h produced the target stereoselective bi-spirooxindole-engrafted rhodanine analogs 5a-d. The reaction proceeded initially with isatins 2a-d reacting with the secondary amino acid (L-thioproline) 1 to generate the azomethine ylides in situ (AYs 3a-d). In the second step, the generated azomethine ylides (AYs 3a-d) reacted with the arylidene rhodanine analog 4 with four possible approach modes: two regioisomeric approach modes of the aryl-substituted carbon of 4 to the two carbons of the azomethine ylides (AYs) framework of AYs 3a-d, named ortho and meta, and two stereoisomeric approach modes of the C-C double bond of 4 to AYs 3a-d, named endo and exo. The four possible alternative reaction paths are shown in Scheme 2. The formation of the bi-spirooxindole-engrafted rhodanine analogs 5a-d indicates that this 32CA reaction is completely ortho regioselective and exo stereoselective, taking place via the transition state structure TS-ox. The spectral data of the bi-spirooxindole-engrafted rhodanine analogs 5a-d were assigned and elucidated and were found to be in good agreement with the proposed structure. The chemical structure of 5a was further confirmed by single-crystal X-ray diffraction analysis. Scheme 2. Synthesis and plausible mechanism of the desired spirooxindole derivatives 5a-d.

The X-ray structure of 5a showing atom numbering and thermal ellipsoids drawn at the 30% probability level is shown in Figure 1 . The structure agrees very well with the spectral analyses. The crystal data and structure refinement details are depicted in Table  1 , while selected bond distances and angles are listed in Table 2 . The compound crystal-Scheme 2. Synthesis and plausible mechanism of the desired spirooxindole derivatives 5a-d. 

The X-ray structure of 5a showing atom numbering and thermal ellipsoids drawn at the 30% probability level is shown in Figure 1 . The structure agrees very well with the spectral analyses. The crystal data and structure refinement details are depicted in Table 1 , while selected bond distances and angles are listed in Table 2 . The compound crystallized in the monoclinic system and P2 1 space group with unit cell parameters of a = 14.0397(2) Å, b = 15.9918(2) Å, c = 14.1488(2) Å and β = 109.439(2) • . The unit cell volume was 2995.61(8) Å 3 and the number of molecules per unit cell was 4. The asymmetric unit comprised two crystallographically independent molecules of 5a and two water molecules. The geometrical parameters of the two units are marginally different. The reported Xray structure revealed very well the presence of four asymmetric carbon centers in the spiro-system, namely C11, C12, C19 and C22 ( Figure 1 ). Table 3 ). The hydrogen bond network is shown in Figure 2 . 

In light of the great importance of intermolecular interactions for crystal stability, Hirshfeld calculations were employed in order to analyze the molecular packing in the studied system. The Hirshfeld surfaces of 5a are shown in Figure S1 (Supplementary Materials). The crystal structure of 5 indicated the presence of two molecules as an asymmetric unit; hence, the Hirshfeld surfaces of both units are presented in this figure.

In Figure 3 , the d norm maps of both molecules are also presented. The red regions indicate the atomic sites included in strong intermolecular interactions with the neighboring molecules. It is clear that the molecular packing in the studied system is controlled by many short contacts such as Cl.. 

In light of the great importance of intermolecular interactions for crystal stability, Hirshfeld calculations were employed in order to analyze the molecular packing in the studied system. The Hirshfeld surfaces of 5a are shown in Figure S1 (Supplementary Materials). The crystal structure of 5 indicated the presence of two molecules as an asymmetric unit; hence, the Hirshfeld surfaces of both units are presented in this figure.

In Figure 3 , the dnorm maps of both molecules are also presented. The red regions indicate the atomic sites included in strong intermolecular interactions with the neighboring molecules. It is clear that the molecular packing in the studied system is controlled by many short contacts such as Cl.. In addition, fingerprint plots gave a quantitative summary for all intermolecular interactions. Decomposed fingerprint plots enabled us to calculate the percentage of each 

In light of the great importance of intermolecular interactions for crystal stability Hirshfeld calculations were employed in order to analyze the molecular packing in the studied system. The Hirshfeld surfaces of 5a are shown in Figure S1 (Supplementary Materials). The crystal structure of 5 indicated the presence of two molecules as an asymmetric unit; hence, the Hirshfeld surfaces of both units are presented in this figure.

In Figure 3 , the dnorm maps of both molecules are also presented. The red regions indicate the atomic sites included in strong intermolecular interactions with the neighboring molecules. It is clear that the molecular packing in the studied system is controlled by many short contacts such as Cl.. In addition, fingerprint plots gave a quantitative summary for all intermolecular interactions. Decomposed fingerprint plots enabled us to calculate the percentage of each In addition, fingerprint plots gave a quantitative summary for all intermolecular interactions. Decomposed fingerprint plots enabled us to calculate the percentage of each intermolecular interaction ( Figure 4 ). The distribution of all possible intermolecular interactions is shown in Figure 5 H2 (1.995 Å) hydrogen bonds are significantly short and hence have great importance in the molecular packing. The two spikes in the decomposed fingerprint plot of the O...H interactions indicate not only strong hydrogen bonds but also the molecule acting as both the hydrogen bond donor and acceptor. In the same table, all contacts with a distance shorter than the vdWs radii sum of the interacting atoms are also presented. The only exception is the N2 . . . H6BA (2.683Å) contact, which has a slightly longer contact distance than the vdWs radii sum of the N and H atoms (2.64 Å). 

The molecular geometry of 5 was calculated and the resulting optimized structure is shown in Figure 6 (right part). As can be seen from Table 5 , there are good agreements between the calculated and experimental bond distances. In addition, there are straight-line correlations between the calculated and experimental bond distances (R 2 = 0.9932) and angles (R 2 = 0.8759), as shown in Figure 7 . In addition, the structure matching between the calculated and experimental structures is shown in Figure 6 (left part). Generally, the structural parameters indicate good agreement between the calculated results and X-ray geometry, with some little variations which could be attributed to the packing effects. Figure 6 . The charges at the different atoms of 5 are collated in Table 6 . The results indicate the presence of a number of negatively charged atomic sites, including N, O and the majority of C sites. There are five oxygen atoms and five nitrogen sites in 5a where the exo carbonyl oxygen and the neighboring amidic nitrogen atoms are the most negative, respectively. All carbon atoms bonded to N or O atoms have highly positive natural charges. In The charges at the different atoms of 5 are collated in Table 6 . The results indicate the presence of a number of negatively charged atomic sites, including N, O and the majority of C sites. There are five oxygen atoms and five nitrogen sites in 5a where the exo carbonyl oxygen and the neighboring amidic nitrogen atoms are the most negative, respectively. All carbon atoms bonded to N or O atoms have highly positive natural charges. In addition, the NH protons are the most positive hydrogen atoms (0.4475 and 0.4196 e). These results are in good agreement with the molecular electrostatic potential map (MEP) shown in Figure 8 . In the MEP, the red colored areas are related to regions of high electron density while the blue colored areas are related to highly electron-deficient regions. The red regions are related to carbonyl oxygen while the blue regions are close to the NH protons. The studied compound has a polar nature with a net dipole moment of 2.5735 Debye, and the direction of the dipole moment vector is shown in Figure 8 . Figure 6 . 

The reactivity indices defined within the conceptual DFT (CDFT) [36, 37] have been shown to be powerful tools to understand the reactivity in polar reactions. The global reactivity indices, namely the electronic chemical potential µ, chemical hardness η, global electrophilicity ω and global nucleophilicity N, for AYs 3a-d and ethylene 4 involved in this 32CA reaction and the reduced models AY 6 and ethylene 7 are gathered in Table 7 . The electronic chemical potentials [36] µ of AYs 3a-d, ranking from µ = −3.26 (3c) to −3.71 (3b) eV, are higher than that of the ethylene derivative 4, µ = −4.30 eV, indicating that along a polar 32CA reaction, the global electron density transfer (GEDT) [38] will take place from AYs 3a-d to the ethylene derivative 4, with the 32CA reaction being classified as the forward electron density flux (FEDF) [39] .

The electrophilicity index ω [40] of AYs 3a-d ranges from 1.62 (3c) to 2.14 (2b) eV, thus being classified as strong electrophiles within the electrophilicity scale [37] . On the other hand, AYs 3a-d present a nucleophilicity index N [41] between 3.79 (3d) and 4.23 (3c) eV, thus being classified as strong nucleophiles within the nucleophilicity scale [37] . The strong nucleophilic character of AYs 3a, 3c and 3d, higher than 4.0 eV, allows for their classification as supernucleophiles [42] . Note that the presence of the strong electron-withdrawing NO 2 group in AY 3b markedly increases its electrophilicity index ω to 2.14 eV and decreases its nucleophilicity index N to 3.79 eV. In spite of these behaviors, AY 3b is a very strong nucleophile participating in polar 32CA reactions of FEDF.

The ethylene derivative 4 presents an electrophilicity index ω of 2.51 eV, being classified as a strong electrophile within the electrophilicity scale. On the other hand, it presents a nucleophilicity index N of 3.09 eV, being classified also as a strong nucleophile within the nucleophilicity scale.

The similar electrophilicity ω and nucleophilicity N indices of the reduced models AY 6 and ethylene 7 to those of the experimental AYs 3a-d and ethylene 4 (see Table 7 ) indicate that the corresponding 32CA reactions will be very similar.

The highly nucleophilic character of AYs 3a-d and 6 together with the strong electrophilic character of the ethylene derivatives 4 and 7 indicate that the corresponding 32CA reactions will have a highly polar character, being classified as FEDF [39] .

In a polar reaction involving non-symmetric species, the most favorable reaction path involves the two-center interaction between the most electrophilic and the most nucleophilic centers [43] . Many studies have shown that analysis of the electrophilic P k + and nucleophilic P k − Parr functions [44] , resulting from the excess of spin electron density gathered via the GEDT [38] , is one of the most accurate and insightful tools for the analysis of the local reactivity in polar and ionic processes. Hence, according to the characteristics of the reagents, the P k − nucleophilic Parr functions of AY 3a and the P k + electrophilic Parr functions of the ethylene derivative 4 were analyzed (see Figure 9 ).

The two carbons of AY 3a are nucleophilically activated by P k − = 0.35 (C22) and 0.30 (C19), the exocyclic C22 carbon being slightly more activated. Note that the N5 nitrogen is deactivated by P k − = −0.11. On the other hand, the conjugated C12 carbon of the ethylene derivative 4 is electrophilically activated, P k + = 0.30, while the exocyclic C11 carbon in marginally electrophilically activated, P k + = 0.09. Analysis of the Parr functions suggest that the meta reaction paths will be slightly more favorable than the ortho ones. However, steric hindrances present along then meta reaction paths will be responsible of the meta regioselectivity experimentally observed. These steric hindrances can also explain that the ortho/exo reaction path will be favored with respect to the ortho/endo one. The two carbons of AY 3a are nucleophilically activated by Pk − = 0.35 (C22) and 0.30 (C19), the exocyclic C22 carbon being slightly more activated. Note that the N5 nitrogen is deactivated by Pk − = −0.11. On the other hand, the conjugated C12 carbon of the ethylene derivative 4 is electrophilically activated, Pk + = 0.30, while the exocyclic C11 carbon in marginally electrophilically activated, Pk + = 0.09.

Analysis of the Parr functions suggest that the meta reaction paths will be slightly more favorable than the ortho ones. However, steric hindrances present along then meta reaction paths will be responsible of the meta regioselectivity experimentally observed. These steric hindrances can also explain that the ortho/exo reaction path will be favored with respect to the ortho/endo one.

For the study of the mechanism of the 32CA reactions of AYs 3a-d with ethylene 6, a reduced model was selected (see Scheme 3) . In this reaction model, the substituents present at the aromatic ring of AYs 3a-d were removed, while the chain present at the amide N3 nitrogen of the ethylene derivative 4 was replaced by a methyl group in ethylene 7. Analysis of the CDFT indices of AY 6 and ethylene 7 showed that they have a similar nucleophilic and electrophilic character compared to the experimental AYs 3a-d and ethylene 4, indicating that the former will present a similar reactivity (see Table 7 ). Analysis of the reaction path indicated that this 32CA reaction takes place though a one-step mechanism. The B3LYP/6-31G(d) relative energies are given in Scheme 3. The activation energy associated to this 32CA reaction is 8.6 kcal·mol −1 , with the reaction being exothermic by −21.7 kcal·mol −1 . 

For the study of the mechanism of the 32CA reactions of AYs 3a-d with ethylene 6, a reduced model was selected (see Scheme 3) . In this reaction model, the substituents present at the aromatic ring of AYs 3a-d were removed, while the chain present at the amide N3 nitrogen of the ethylene derivative 4 was replaced by a methyl group in ethylene 7. Analysis of the CDFT indices of AY 6 and ethylene 7 showed that they have a similar nucleophilic and electrophilic character compared to the experimental AYs 3a-d and ethylene 4, indicating that the former will present a similar reactivity (see Table 7 ). Analysis of the reaction path indicated that this 32CA reaction takes place though a one-step mechanism. The B3LYP/6-31G(d) relative energies are given in Scheme 3. The activation energy associated to this 32CA reaction is 8.6 kcal·mol −1 , with the reaction being exothermic by −21.7 kcal·mol −1 . The two carbons of AY 3a are nucleophilically activated by Pk − = 0.35 (C22) and 0.30 (C19), the exocyclic C22 carbon being slightly more activated. Note that the N5 nitrogen is deactivated by Pk − = −0.11. On the other hand, the conjugated C12 carbon of the ethylene derivative 4 is electrophilically activated, Pk + = 0.30, while the exocyclic C11 carbon in marginally electrophilically activated, Pk + = 0.09.

Analysis of the Parr functions suggest that the meta reaction paths will be slightly more favorable than the ortho ones. However, steric hindrances present along then meta reaction paths will be responsible of the meta regioselectivity experimentally observed. These steric hindrances can also explain that the ortho/exo reaction path will be favored with respect to the ortho/endo one.

For the study of the mechanism of the 32CA reactions of AYs 3a-d with ethylene 6, a reduced model was selected (see Scheme 3) . In this reaction model, the substituents present at the aromatic ring of AYs 3a-d were removed, while the chain present at the amide N3 nitrogen of the ethylene derivative 4 was replaced by a methyl group in ethylene 7. Analysis of the CDFT indices of AY 6 and ethylene 7 showed that they have a similar nucleophilic and electrophilic character compared to the experimental AYs 3a-d and ethylene 4, indicating that the former will present a similar reactivity (see Table 7 ). Analysis of the reaction path indicated that this 32CA reaction takes place though a one-step mechanism. The B3LYP/6-31G(d) relative energies are given in Scheme 3. The activation energy associated to this 32CA reaction is 8.6 kcal·mol −1 , with the reaction being exothermic by −21.7 kcal·mol −1 . The geometry of the transition state structure TS is given in Figure 10 . The C-C distances between the two pairs of interacting carbons are 2.116 (C12-C19) and 2.875 (C11-C22) Å. These distances indicate that this TS is associated with a highly asynchronous C-C single-bond formation process, in which the shorter C12-C19 distance corresponds to that involving the most electrophilic carbon C12 of the ethylene derivative 7. These C-C distances, which are higher than 2.1 Å, indicate that the formation of the two C-C single bonds had not begun at the TS [38] . Analysis of the intrinsic reaction coordinates [45] associated to the highly asynchronous TS indicated that this 32CA reaction takes place through a non-concerted two-stage one-step mechanism [46] in which formation of the second C11-C22 single bond begins when the first C12-C19 single bond is completely formed.

These C-C distances, which are higher than 2.1 Å, indicate that the formation of the two C-C single bonds had not begun at the TS [38] . Analysis of the intrinsic reaction coordinates [45] associated to the highly asynchronous TS indicated that this 32CA reaction takes place through a non-concerted two-stage one-step mechanism [46] in which formation of the second C11-C22 single bond begins when the first C12-C19 single bond is completely formed. Finally, analysis of the GEDT [38] at the TS allowed for assessment of the polar character of this 32CA reaction. GEDT values lower than 0.05 e correspond with non-polar processes, while values higher than 0.20 e correspond with polar processes. The GEDT value at the TS was 0.26 e. This high GEDT value is a consequence of the supernucleophilic character of AY 6 and the strong nucleophilic character of ethylene 7. The flux of the electron density, which goes from AY 6 to ethylene 7, classifies this 32CA reaction as FEDF, in clear agreement with the analysis of the CDFT indices.

Isatin derivatives 2a-d and thioproline 1 are commercially available (Sigma-Aldrich Chemie GmbH, Riedstr, Germany; and Alfa Aesar GmbH & Co KG, Karlsruhe, Germany). 1 H-NMR and 13 C-NMR spectra were recorded in DMSO-d6 (Jeol Spectrometer (400 MHz), (Jeol, Tokyo, Japan). X-ray diffraction data were collected on a Rigaku Oxford Diffraction Supernova diffractometer and processed with CrysAlisPro software v. 1.171.41.93a (Rigaku Oxford Diffraction, Yarnton, UK, 2020) using Cu Kα radiation.

The desired starting material 4 was synthesized according to the reported literature [29] .

A mixture of three components including isatin derivatives 2a-d (0.5 mmol), L-thioproline 1 (66.5 mg, 0.5 mmol) and the compound 4 (204.5 mg, 0.5 mmol) was re- Finally, analysis of the GEDT [38] at the TS allowed for assessment of the polar character of this 32CA reaction. GEDT values lower than 0.05 e correspond with non-polar processes, while values higher than 0.20 e correspond with polar processes. The GEDT value at the TS was 0.26 e. This high GEDT value is a consequence of the supernucleophilic character of AY 6 and the strong nucleophilic character of ethylene 7. The flux of the electron density, which goes from AY 6 to ethylene 7, classifies this 32CA reaction as FEDF, in clear agreement with the analysis of the CDFT indices.

Isatin derivatives 2a-d and thioproline 1 are commercially available (Sigma-Aldrich Chemie GmbH, Riedstr, Germany; and Alfa Aesar GmbH & Co KG, Karlsruhe, Germany). 1 H-NMR and 13 C-NMR spectra were recorded in DMSO-d 6 (Jeol Spectrometer (400 MHz), (Jeol, Tokyo, Japan). X-ray diffraction data were collected on a Rigaku Oxford Diffraction Supernova diffractometer and processed with CrysAlisPro software v. 1.171.41.93a (Rigaku Oxford Diffraction, Yarnton, UK, 2020) using Cu Kα radiation.

The desired starting material 4 was synthesized according to the reported literature [29] .

A mixture of three components including isatin derivatives 2a-d (0.5 mmol), Lthioproline 1 (66.5 mg, 0.5 mmol) and the compound 4 (204.5 mg, 0.5 mmol) was refluxed in an oil bath for 2 h. After completion of the reaction as evident from tlc (TLC Eluent: Ethyl acetate: n-Hexane 40%), the reaction mixture was kept at room temperature overnight for slow evaporation, and the solid crystalline materials were filtered off to obtain compounds 5a-d without any further purification as faint yellow solid compounds with 80-90% chemical yield. According to the general method, the 6-chloro-isatin 2a (90.5 mg, 0.5 mmol) was utilized and the final compound 5a was obtained as a faint yellow solid compound with 85% chemical yield. 1 

The crystal of 5a was immersed in cryo-oil, mounted in a loop and measured at a temperature of 120 K. The X-ray diffraction data were collected on a Rigaku Oxford Diffraction Supernova diffractometer using Cu Kα radiation. The CrysAlisPro [47] software package was used for cell refinement and data reduction. An analytical absorption correction (CrysAlisPro [47] ) was applied to the intensities before structure solution. The structure was solved by the intrinsic phasing (SHELXT [48] ) method. Structural refinement was carried out using SHELXL [49] software with the SHELXLE [50] graphical user interface. The crystal was solved in the space group P2 1 as an inversion twin. The BASF value was refined to 0.1478. The asymmetric unit contained two independent organic molecules and two water molecules. The ring C19B-C29B-S6B-C21B-N4B was slightly disordered, leading to small positive and negative electron densities near the atom S6B. No disorder model was used in the final refinement. The NH hydrogen atom was located from the difference Fourier map and refined isotropically. The H 2 O hydrogen atoms were also located from the difference Fourier map but constrained to ride on their parent oxygen with U iso = 1.5 U eq (parent atom). All other hydrogen atoms were positioned geometrically and constrained to ride on their parent atoms, with C-H = 0.95-1.00 Å and U iso = 1.2·U eq (parent atom).

The topology analyses were performed using the Crystal Explorer 17.5 program [51] .

All DFT calculations were performed using the Gaussian 09 software package [52] utilizing the B3LYP/6-31G(d,p) method. The geometries were visualized with the GaussView program [53] . Global and local CDFT indices were calculated using the equations given in [37] .

The straightforward synthesis of rhodanine-substituted spiro[pyrrolidine-oxindole]carboxamide/morpholine derivatives 5a-d using the 32CA reaction methodology was successfully achieved. The supramolecular structure of the studied compound 5a was analyzed using Hirshfeld calculations. Additionally, the DFT-calculated geometry was found to agree with the reported X-ray structure. The CDFT analysis of the reagents indicated that this 32CA reaction will have a highly polar character due to the supernucleophilic character of AYs 3a-d and the strong electrophilic character of the ethylene derivative 4, being classified as FEDF. The analysis of the Parr functions suggested that the meta reaction paths will be slightly more favorable than the ortho ones. However, steric hindrances present along the meta reaction paths are responsible for the meta regioselectivity experimentally observed. Finally, the MEDT analysis of these 32CA reactions indicated that they take place through a polar non-concerted two-stage one-step mechanism associated with the nucleophilic attack of the least substituted carbon of AYs 3a-d on the β-conjugated position of the ethylene derivative 4 as a consequence of the supernucleophilic character of AYs 3a-d and the strong electrophilic character of the ethylene 4.

Supplementary Materials: The following are available online, Figure S1 : Hirshfeld surfaces of 5a, Figure S2 -17: Copy of spectrum data (NMR, IR, MS) for the synthesized compounds 5a-d; Tables S1-S2: computational data of compound 5a for the intermolecular interactions and their percentages, also for the bond angles.

77 (d, J = 7.8 Hz, 1H), 2.37-2.30 (m, 6H); 13 C-NMR (100 Hz

-yl)-N-(2-morpholinoethyl)aceta mide 5b According to the general method, the 5-nitro

400 Hz, DMSO-d 6 ): δ 11.70 (s, 1H), 8.32-8.23 (m, 1H), 8.07 (t, J = 5.3 Hz, 2H), 7.97 (s, 1H), 7.68 (d, J = 8.6 Hz, 1H), 7.62 (d, J = 8.6 Hz, 1H), 7.52-7.41 (m, 4H), 7.08 (d, J = 8.7 Hz, 1H)

100 Hz, DMSO-d 6 ): δ 13 C-NMR (101 MHz

Chemical Formula: C 29 H 29 ClN 6 O 7 S 2

-chlorophenyl)-5-fluoro-2,2 ,4 -trioxo-7 ,7a -dihydro-1 H,3 H-disp iro

According to the general method, the 5-fluoro-isatin 2c (82.5 mg, 0.5 mmol) was

400 Hz, DMSO-d 6 ): δ 11.05 (s, 1H), 8.08 (t, J = 5.7 Hz, 1H), 7.45 (q, J = 8.5 Hz, 4H), 7.20 (td, J = 8.7, 2.8 Hz, 1H)

Chemical Formula: C 29 H 29 ClFN 5 O 5 S 2

-chlorophenyl)-2,2 ,4 -trioxo-7 ,7a -dihydro-1 H,3 H-disp iro

According to the general method, the 5-chloro-isatin 2d (90.5 mg, 0.5 mmol) was

400 Hz, DMSO-d 6 ): δ 11.16 (s, 1H), 8.10 (t, J = 5.6 Hz, 1H), 7.51-7.36 (m, 5H), 7.22 (s, 1H), 6.90 (d, J = 8.1 Hz, 1H)

Recent developments with rhodanine as a scaffold for drug discovery

A conceptual review of rhodanine: Current applications of antiviral drugs, anticancer and antimicrobial activities

Synthesis and biological evaluation of novel carbazole-rhodanine conjugates as topoisomerase II inhibitors

Synthesis and in vitro evaluation of novel rhodanine derivatives as potential cholinesterase inhibitors

Synthesis and biological evaluation of 5-benzylidenerhodanine-3-acetic acid derivatives as AChE and 15-LOX inhibitors

Design, synthesis and structure-activity studies of rhodanine derivatives as HIV-1 integrase inhibitors

Synthesis of new spirooxindole-pyrrolothiazoles derivatives: Anti-cancer activity and molecular docking

Clinical investigation of epalrestat, an aldose reductase inhibitor, on diabetic neuropathy in Japan: Multicenter study

Synthesis, α-glucosidase inhibition and molecular docking studies of novel thiazolidine-2,4-dione or rhodanine derivatives

Synthesis and biological evaluation of thiazolidine-2,4-dionepyrazole conjugates as antidiabetic, anti-inflammatory and antioxidant agents

Discovery of 1,3-diphenyl-1H-pyrazole derivatives containing rhodanine-3-alkanoic acid groups as potential PTP1B inhibitors

New spiro-oxindole constructed with pyrrolidine/thioxothiazolidin-4-one derivatives: Regioselective synthesis, x-ray crystal structures, Hirshfeld surface analysis, DFT, docking and antimicrobial studies

Synthesis and characterization of a spiroindolone pyrothiazole analog via X-ray, biological, and computational studies

Synthesis of novel dispiropyrrolothiazoles by three-component 1,3-dipolar cycloaddition and evaluation of their antimycobacterial activity

Antiproliferative effects of mitraphylline, a pentacyclic oxindole alkaloid of Uncaria tomentosa on human glioma and neuroblastoma cell lines

Substituted spirooxindole derivatives as potent anticancer agents through inhibition of phosphodiesterase 1

A strategic approach to the synthesis of functionalized spirooxindole pyrrolidine derivatives: In vitro antibacterial, antifungal, antimalarial and antitubercular studies

Novel spirooxindole-pyrrolidine compounds: Synthesis, anticancer and molecular docking studies

A facile ionic liquid promoted synthesis, cholinesterase inhibitory activity and molecular modeling study of novel highly functionalized spiropyrrolidines

Molecular hybridization design and synthesis of novel spirooxindole-based MDM2 inhibitors endowed with BCL2 signaling attenuation; a step towards the next generation p53 activators

Design and synthesis of new substituted spirooxindoles as potential inhibitors of the MDM2-p53 interaction

Design, synthesis, chemical and biochemical insights into novel hybrid spirooxindole-based p53-MDM2 inhibitors with potential Bcl2 signaling attenuation

Regio-and stereoselective synthesis of a new series of spirooxindole pyrrolidine grafted thiochromene scaffolds as potential anticancer agents

Three-component access to functionalized spiropyrrolidine heterocyclic scaffolds and their cholinesterase inhibitory activity

Synthesis, X-ray single crystal, conformational analysis and cholinesterase inhibitory activity of a new spiropyrrolidine scaffold tethered benzo

Synthesis of a new class of spirooxindole-benzo[b]thiophene-based molecules as acetylcholinesterase inhibitors

Drug repurposing of lactoferrin combination in a nanodrug delivery system to combat severe acute respiratory syndrome coronavirus-2 infection

Synthesis, antidiabetic activity and molecular docking study of rhodanine-substitued spirooxindole pyrrolidine derivatives as novel α-amylase inhibitors

Synthesis and antimicrobial activity of a new series of thiazolidine-2,4-diones carboxamide and amino acid derivatives

Contribution of hydrogen bonds to protein stability

Hydrogen bonds as design elements in organic chemistry

Molecular electron density theory: A modern view of reactivity in organic chemistry

Unravelling the mysteries of the [3 + 2] cycloaddition reactions

Understanding the high reactivity of the azomethine ylides in [3 + 2] Cycloaddition Reactions

Unveiling the reactivity of cyclic azomethine ylides in [3 + 2] cycloaddition reactions within the molecular electron density theory

Density-Functional Theory of Atoms and Molecules

Applications of the conceptual density functional theory indices to organic chemistry reactivity

A new C-C bond formation model based on the quantum chemical topology of electron density

A molecular electron density theory study of the participation of tetrazines in aza-Diels-Alder reactions

Understanding the reactivity of captodative ethylenes in polar cycloaddition reactions. A theoretical study

A close look to the oxaphosphetane formation along the Wittig reaction: A [2 + 2] cycloaddition?

A theoretical study on the regioselectivity of 1,3-dipolar cycloadditions using DFT-based reactivity indexes

Understanding the local reactivity in polar organic reactions through electrophilic and nucleophilic Parr functions

Formulation of the reaction coordinate

Understanding the participation of quadricyclane as nucleophile in polar [2σ + 2σ + 2π] cycloadditions toward electrophilic π molecules

Rikagu Oxford Diffraction

SHELXT-Integrated space-group and crystal-structure determination

Crystal structure refinement with SHELXL

ShelXle: A Qt graphical user interface for SHELXL

Gaussian 09, Revision A02

The authors would like to extend their sincere appreciation to the Researchers Supporting Project (RSP-2021/64), King Saud University, Riyadh, Saudi Arabia, and the Ministerio de Ciencias, Innovación y Universidades of the Spanish Government, project PID2019-110776GB-I00 (AEI/FEDER, UE).

The authors declare no conflict of interest.Sample Availability: Samples of the compounds 5a-d are available from the authors.