Nourhène Chouchène1, Amani Toumi1, Sarra Boudriga1, Hayet Edziri2, Mansour Sobeh3, Mohamed A O Abdelfattah4, Moheddine Askri1, Michael Knorr5, Carsten Strohmann6, Lukas Brieger6, Armand Soldera7. 1. Laboratory of Heterocyclic Chemistry Natural Product and Reactivity (LR11ES39), Department of Chemistry, Faculty of Science of Monastir, University of Monastir, Monastir 5019, Tunisia. 2. Laboratoire des Maladies Transmissibles et des Substances Biologiquement Actives, Faculté de Pharmacie, Monastir 5000, Tunisia. 3. AgroBioSciences Research, Mohammed VI Polytechnic University, Lot 660-Hay MoulayRachid, Ben Guerir 43150, Morocco. 4. College of Engineering and Technology, American University of the Middle East, Kuwait. 5. Institut UTINAM-UMR CNRS 6213, Université Bourgogne Franche-Comté, 16 Route de Gray, 25030 Besançon, France. 6. Faculty of Chemistry, Inorganic Chemistry, Technical University Dortmund, Otto-Hahn-Strasse 6, 44227 Dortmund, Germany. 7. Laboratory of Physical Chemistry of Matter, Department of Chemistry, Université de Sherbrooke, Sherbrooke, QC J1K 2R1, Canada.
Abstract
A novel series of 14 spiropyrrolidines bearing thiochroman-4-one/chroman-4-one, and oxindole/acenaphthylene-1,2-dione moieties were synthesized and characterized by spectroscopic techniques, as well as by three X-ray diffraction studies, corroborating the stereochemistry. Quantum chemical calculations studies, using the DFT approach, were performed to rationalize the stereochemical outcome. These N-heterocycles were evaluated for their antibacterial and antifungal activities against some pathogenic organisms. Several compounds displayed moderate to excellent activity towards the screened microbe strains in the study compared to Amoxicillin (AMX), Ampicillin (AMP), and Amphotericin B. Furthermore, a structural activity relationship (SAR) was established considering the synthesized compounds. Pharmacokinetic studies reveal that these derivatives exhibit an acceptable predictive ADMET profile (Absorption, Distribution, Metabolism, Excretion and Toxicity) and good drug-likeness.
A novel series of 14 spiropyrrolidines bearing thiochroman-4-one/chroman-4-one, and oxindole/acenaphthylene-1,2-dione moieties were synthesized and characterized by spectroscopic techniques, as well as by three X-ray diffraction studies, corroborating the stereochemistry. Quantum chemical calculations studies, using the DFT approach, were performed to rationalize the stereochemical outcome. These N-heterocycles were evaluated for their antibacterial and antifungal activities against some pathogenic organisms. Several compounds displayed moderate to excellent activity towards the screened microbe strains in the study compared to Amoxicillin (AMX), Ampicillin (AMP), and Amphotericin B. Furthermore, a structural activity relationship (SAR) was established considering the synthesized compounds. Pharmacokinetic studies reveal that these derivatives exhibit an acceptable predictive ADMET profile (Absorption, Distribution, Metabolism, Excretion and Toxicity) and good drug-likeness.
The emergence of drug-resistant pathogens has threatened human ability to treat common infections, and, thus, antimicrobial resistance (AMR) became a health threat all over the globe. The rapid spread of Staphylococcus aureus, for instance, is especially alarming, since it is considered to be among the most common causes of infections in health-care facilities [1,2]. Furthermore, patients suffering with infections caused by methicillin-resistant Staphylococcus aureus (MRSA) are more likely to die, compared to others infected with drug-sensitive infections, by about 64% [3,4]. Consequently, antimicrobial medications including antibiotics are getting more and more ineffective and, thus, infections become difficult or nearly impossible to treat. This lowered efficacy urgently increases the demand for new therapeutic agents.In recent years, spiropyrrolidine scaffolds, especially those fused to oxindole or acenaphthylene-1(2H)-one moieties, and occurring in a huge number of pharmacological entities and natural compounds with significant bioactive properties, have attracted the attention of synthetic chemists [5,6,7,8,9,10,11,12,13] (Figure 1).
Figure 1
(A) Biologically relevant thiochroman, thiochromanone, and chromanone derivatives, (C) Biologically relevant spiropyrrolidines bearing oxindole or acenaphthylene-1(2H)-one moieties, and (B) Our target compounds.
On the other hand, 4-chromanone derivatives represent privileged scaffolds in heterocyclic chemistry and for drug discovery. They are used as versatile intermediates for the synthesis of many natural products [14,15,16,17,18,19] and constitute interesting building blocks in drug design and organic synthesis [20,21,22,23]. They also exhibit significant biological activities allowing their use as anticancer [24,25,26,27], antifungal and antibacterial [28,29,30,31], anti-inflammatory [32,33], antidiabetic [34,35], anti-leishmanial (caused by protozoan parasites) [36,37], and insecticidal agents [38] (Figure 1).Thiochromanones, the thio-analogs of chromanones, feature interesting scaffolds, which have been reported to exert a plethora of pronounced biological and pharmacological properties including antimicrobial [39,40,41,42], anticancer [43,44,45,46], human steroid sulfatase inhibition [47,48], α-adrenergic antagonist [49], antiviral [50], and insecticidal activities [51] (Figure 1).In light of the significance of spiropyrrolidine and thiochromanone/chromanone heterocycles in medicinal chemistry, we felt a great potential to combine these two moieties into one single scaffold to design a novel class of hybrid spiroheterocycles, which may exhibit interesting biological profiles. The latter were obtained via multicomponent reactions (MCRs) of arylidene thiochroman-4-ones/chroman-4-ones, glycine ethyl ester, and cyclic 1,2-diketones (isatin derivatives or acenaphthenequinone) (Scheme 1). We are aware that related spiropyrrolidine(s) fused with thiochroman-4-one/chroman-4-one moieties have been prepared by Subramaniyan [52] and Rani [53], however their pharmacological profile has, up to now, never been explored.
Scheme 1
Synthesis scheme for the target compounds 4 and 9.
In continuation of our research interests in the synthesis of novel bioactive N-spiroheterocycles [5,54,55,56,57], we herein report the synthesis of diverse thiochromanone/chromanone substituted spiropyrrolidines (Figure 1). The latter were obtained via multicomponent reactions (MCRs) of arylidene thiochroman-4-ones/chroman-4-ones, glycine ethyl ester, and cyclic 1,2-diketones (isatin derivatives or acenaphthenequinone) (Scheme 1). These new scaffolds were tested in vitro for their antimicrobial potential. DFT calculations were carried out to obtain the optimized molecular structures of the resulting spiranic thiochromanone/chromanone-linked spiropyrrolidines. Furthermore, frontiers molecular orbitals (FMO) and molecular electrostatic potential (MEP) calculations were performed for the structural exploration of these compounds. A drug-likeness analysis was also performed.
2. Results and Discussion
2.1. Synthetic Chemistry
The synthetic strategy for our target compounds 4 and 9 is shown in Scheme 1. It is based on a multicomponent process involving the 1,3-dipolar cycloaddition reaction between exocyclic arylated alkenes derived from thiochromanone 1 or chromanone 6 and azomethine ylides. The latter were generated in situ from condensation of glycine ethyl ester 2 and a 1,2-diketone, namely isatin 3 or acenaphthenequinone 3′ (Scheme 1).At the onset of our investigation, we attempted to access the hybrid thiochromanone-grafted spiropyrrolidines according to our recently established reaction conditions for MCR of arylidene rhodanines, glycine ethyl ester, and cyclic 1,2-diketones [5]. For that purpose, an equivalent amount of 3-benzylidenethiochroman-4-one 1a, glycine ethyl ester 2, and isatin 3 was refluxed in acetonitrile in the presence of triethylamine for 2 h. The reaction smoothly afforded the desired spiroheterocycle 4a in 80% yield (Table 1). Subsequently, a variety of dipolarophiles 1 were subjected to these reaction conditions. As shown in Table 1, this process is applicable to a wider series of arylidene thiochromanones 1 bearing electronically different substituents, or groups, on the phenyl ring, affording in good yields the targeted spiropyrrolidine products 4b–d in a high regio- and diastereoselective manner as racemic mixtures. To expand the scope of this protocol, we further examined the multicomponent reaction with acenaphthenequinone 3′ as an alternative 1,2-diketone, instead of isatin 3 (Table 1), to synthetize an original spiropyrrolidines scaffold bearing an acenaphthylen-1(2H)-one core. The azomethine ylide derived from acenaphthenequinone seemed to be well appropriate for this process, affording the corresponding exo-spirropyrrolidine-chromanones 4e and 4f with good yields. Since this straightforward multicomponent [3+2] cycloaddition worked efficiently, we then studied the scope of the reaction using an oxygen surrogate, namely 5-arylidene chroman-4-ones 6, under the same reaction conditions (Table 2).
Table 1
Synthesis of spiropyrrolidine-thiochroman-4-ones hybrids 4.
All reactions were carried out using 1 (1 mmol), 2 (1 mmol), 3, or 3′ (1 mmol), NEt3 (1 mmol) in CH3CN (5.0 mL) at reflux for 2 h. Yields of the isolated cycloadducts are given.
Table 2
Synthesis of acenaphthylen-1(2H)-one-spiropyrrolidine-chromanone hybrids 9.
All reactions were carried out using 6 (1 mmol), 2 (1 mmol), 3, or 3′ (1 mmol), NEt3 (1 mmol) in CH3CN (5.0 mL) at reflux for 2 h. Yields of the isolated cycloadducts are given.
The MCR proceeded readily providing the spiropyrrolidines exo-9a–h as the only detectable diastereoisomers (TLC monitoring) with good to excellent yield (70–94%), regardless of the electronic nature of the p-substituent at the aryl group (H, OCH3, CH3, and Cl) of the dipolarophile 6. On the other hand, the steric property of the acenaphthenequinone 3′ had almost a negligible effect on the efficiency of the 1,3-dipolar cycloaddition reaction, which proceeded smoothly producing the acenaphthylen-1(2H)-one spiropyrrolidine-chromanones 9e–h in excellent yields. All compounds are stable and were obtained as colorless or yellow solids.
2.1.1. Spectroscopic and Crystallographic Characterization of Cycloadducts 4 and 9
The composition and stereochemistry of the chromanone/thiochromanone-grafted spiropyrrolidines was unambiguously elucidated by their spectroscopic data and elemental analyses. The assignments of the chemical shifts from the 1H and 13C NMR spectra of the studied compounds were done based on the literature data [5,53,55] and calculations at the B3LYP/6-311++G(2d,p) level of compounds 4a and 9c (Supplementary Materials; Table S1). As representative examples, relevant 1H and 13C chemical shifts of spiro-compounds 4a and 9c are illustrated in Figure 2 and Figure 3, respectively.
Figure 2
1H (black) and 13C {1H} NMR (green) chemical shifts (ppm) of the selected compound 4a (see also Figure S3 in the Supporting Material).
Figure 3
1H (black) and 13C {1H} NMR (green) chemical shifts (ppm) of the selected compound 9c (see also Figure S20 in the Supporting Material).
In the 1H NMR spectra of 4a and 9c, the pyrrolidinyl protons H-4 and H-5 appear as two mutually coupled doublets at δ 5.14 and 5.02 ppm, and δ 5.02 and 4.79 ppm, respectively. This multiplicity unambiguously confirms the regiochemistry proposed in Figure 2 and Figure 3. If the hypothetical alternative regioisomers 4′,9′ (Scheme 1) would have been formed, the pyrrolidinyl protons should give rise to two singlets in the 1H NMR spectrum. Furthermore, the value of the 3J coupling constants of 10.5 and 10.8 Hz indicate that H-4 and H-5 protons are in trans relationship, in accordance with our earlier reports on related compounds [5,55]. A broadened singlet at δ 8.18 and 7.89 ppm can be assigned to the oxindole N-H proton of 4a and 9c, and a multiplet in the region between 6.59–8.17 ppm and 6.06–7.94 ppm indicates the presence of aromatic protons. Additionally, the 1H NMR spectra show two mutually coupled doublets at δ 2.73 and 3.42 ppm for 4a, and 3.56 and 4.84 ppm for 9c, corresponding to the diastereotopic 2′–CH2 group. The 3J values of 12.3 Hz for 4a and 14.1 Hz for 9c are in the same range with those previously reported values of similar spiropyrrolidine(s) fused with thiochroman-4-one/chroman-4-one moieties [52,53].The proton decoupled 13C NMR spectrum of 9c reveals the presence of two signals at δ 191.8 and δ 178.8 ppm, attributed to the carbonyl groups of chromanone and oxindole moieties, respectively. The resonances corresponding to the two spirocarbons C-2 and C-3 are observed at δ 71.2 and δ 61.4 ppm, respectively. For the assignment of the methoxy, methylene, methine, spiro, and quaternatry carbons, the DEPT-135 spectra shown in Figures S3 and S20 were recorded.The relative configuration of the stereogenic carbons in spiropyrrolidines 4 and 9 was corroborated after determination of the single-crystal X-ray structures of cycloadducts 4a, 4e, and 9c, whose molecular structures are shown in Figure 4, Figure 5 and Figure 6, respectively. The structural analysis reveals that the carbonyl carbon of acenaphthenequinone moiety and (i) the carbonyl carbon group of thiochromanone/chromanone part are in trans-relationship, (ii) the pyrrolidinyl proton attached at C-5 (C-11, C-4, and C-18 in the crystal structures of 4a, 4e, and 9c, respectively) are in cis-relationship (note that the crystallographic atom numbering scheme differs from that used in Figure 2 and Figure 3).
Figure 4
Ball and sticks presentation of the molecular structure of 4a recorded at 100 K. Only one of the two independent molecules is shown. Apart from stereochemically relevant H atoms, all other hydrogen atoms are omitted for clarity. Furthermore, an EtOH molecule of crystallization, present in the asymmetric unit, is not depicted.
Figure 5
Ball and sticks presentation of the molecular structure of 4e recorded at 100 K. Only stereochemically relevant H atoms are shown.
Figure 6
Ball and sticks presentation of the molecular structure of 9c recorded at 100 K. Only one of the two independent molecules is shown. Apart from stereochemically relevant H atoms, all other hydrogen atoms are omitted for clarity.
It is noteworthy that the packing of 4a contains two independent molecules with slightly different metric parameters, and the individual molecules are associated through intermolecular N2-H2····O6 and N4-H4…O2 hydrogen bonding with H····O bond lengths of 1.94 and 2.00 Å, respectively (see Figure S31). A similar supramolecular pairwise N-H····O association occurs also in the crystal structure of 9c (see Figure S32).Accordingly, we propose that the cycloadducts 4 and 9 are formed through an exo-approach between the arylidene thiochromanone/chromanone and the Z,E-dipole, as outlined above in Table 1 and Table 2.
2.1.2. DFT Mechanistic Studies
To better grasp the experimentally observed high regio- and diastereoselectivity in the MCR, DFT calculations were performed using the 6–31G(d,p) basis set and the ωB97xd functional [58]. The effect of solvent (acetonitrile) was taken into consideration using the polarizable continuum model (PCM) approach [59,60]. We have chosen the reaction of dipolarophile 1a, glycine ethyl ester 2, and isatin 3 as model case study, as illustrated in Scheme 2.
Scheme 2
Proposed mechanism for the regio- and stereoisomeric 1,3-dipolar cycloaddition reaction of dipolarophile 1d with azomethine ylide d
In line with our previous study [55], condensation of a glycine ethyl ester and isatin, or acenaphthenequinone, followed by a [1,2]-prototropic sequence furnishes (E,E) and (Z,E) azomethine ylides are the two most stable isomers. Note that the subsequent 1,3-dipolar cycloaddition of these intermediates can a priori occur through both endo-TS and exo-TSs (Scheme 2). Kinetic and thermodynamic calculations of the Gibbs free energy (ΔG), enthalpy (ΔH), and entropy (ΔS), corresponding to all possible isomeric pathways of the reaction and relative energies for the TSs, are summarized in Table 3.
Table 3
Difference in the Gibbs free energies (ΔG, in kcal mol−1), enthalpies (ΔH, in kcal mol−1), entropies (ΔS, in cal K−1 mol−1), and electronic energies (ΔU, in kcal mol−1) for TSs and products calculated at the ωB97xd/6-31G(d,p) level of theory.
ΔG
ΔH
ΔS
ΔU
TS-exo-4d
0.4
−17.6
−50.9
−16.9
TS-endo-4d
2.0
15.4
−48.9
−14.7
TS-exo-5d
15.1
−3.6
−52.7
−2.9
TS-endo-5d
14.9
−0.3
−42.7
0.4
exo-4d
−21.8
−42.4
−58.1
−41.7
endo-4d
−14.6
−34.8
57.0
−34.2
exo-5d
−16.8
−37.9
−59.5
−37.2
endo-5d
−16.3
−36.2
−56.3
−35.5
The analysis of the four located transition states reveals that TS-exo-4d, corresponding to the exo-approach between the (E)-p-bromobenzylidenethiochromanone 1d and the (Z,E)-dipole, presents the lowest barriers with a ΔG value of 0.4 kcal mol−1, favoring the formation of the observed exo-4d regioisomer. Moreover, ΔG values indicate that the product exo-4d, identified as the sole product is more stable (ΔG, −21.8 kcal mol−1) than the other theoretical stereoisomers endo-4d, exo-5d, and endo-5d. Thus, the spiroadduct 4d is kinetically and thermodynamically preferred, which agrees well with our experimental observations.
2.2. Biological Evaluation of the Synthesized Compounds
2.2.1. Antimicrobial Activity
The dispiropyrrolidine derivatives 4 and 9 were evaluated for their antibacterial activity against an assortment of five Gram-positive bacteria (Bacillus subtilis ATCC 6633, Staphylococcus epidermidis CI1232, Staphylococcus aureus ATCC 29213, Staphylococcus aureus ATCC 25923, and Enterococcus faecalis ATCC 29212) and four Gram-negative bacteria (Escherichia coli ATCC 25922, Klebsiella pneumoniae ATCC 4352, Salmonella enterica 800390, and Pseudomonas aeruginosa ATCC 9023). The antifungal activity of the synthesized compounds was also investigated against Candida albicans ATCC 90028, Candida glabrata ATCC 90030, and Candida krusei ATCC 6258. As references, to evaluate and to compare the potency of the tested compounds under the same conditions, the antibiotics Amoxicillin (AMX) and Ampicillin (AMP) were chosen as antibacterial agents. Macrocyclic Amphotericin B was used as antifungal reference.The in vitro antimicrobial activity of the novel dispiropyrrolidines has been assessed by the determination of the minimal inhibitory concentration (MIC) using the microdilution method, as described by Rattan [61]. The results are summarized in Table 4 and Table 5.
Table 4
MIC values (in µg/mL) of the target compounds 4 and 9 against Gram-positive and Gram-negative bacteria activities.
Gram Positive Bacteria
Gram Negative Bacteria
Comp.
B. subtilisATCC 6633
S. epidermidisCI1232
S. aureusATCC 29213
S. aureusATCC 25923
E. faecalisATCC 29212
E. coliATCC 25922
K. pneumoniaeATCC 4352
S. enterica800390
P. aeruginosaATCC 9023
4a
32
32
64
32
32
125
64
64
64
4b
32
32
64
32
32
125
64
64
125
4c
32
32
64
64
32
125
64
64
64
4d
32
32
64
64
32
125
64
64
64
4e
32
32
64
32
64
250
64
125
64
9a
125
250
125
125
125
250
250
550
250
9b
125
250
250
250
125
125
125
550
125
9c
125
125
125
250
64
250
64
550
125
9d
64
64
500
125
64
125
64
550
125
9e
125
250
125
125
64
250
64
550
125
9f
64
64
125
125
125
125
250
550
250
9g
64
64
125
125
125
125
500
500
500
9h
64
64
125
125
250
250
500
500
500
AMXa
64
-
-
-
1
256
256
-
256
AMPa
-
78
50
50
-
-
-
100
-
a Positive control.
Table 5
Antifungal activity of synthesized compounds 4 and 9 (MIC was determined in µg/mL).
Compound
Candida albicans ATCC 90028
Candida glabrata ATCC 90030
Candida krusei ATCC 6258
4a
64
64
64
4b
64
64
64
4c
64
64
64
4d
64
32
32
4e
64
64
64
9a
250
250
250
9b
250
125
125
9c
125
125
125
9d
125
125
250
9e
500
500
500
9f
125
250
125
9g
250
125
250
9h
125
250
250
Amphotericin Ba
500
500
500
a Positive control.
Most of the screened compounds show variable antibacterial activities ranging from poor to excellent against Gram-positive and Gram-negative bacteria with MIC values in the range of 32–250 µg/mL. As shown in Table 4, compounds 4a–e exhibit, in most cases, interesting antibacterial potential. Compounds 4a–d displayed the highest activity (MIC = 32 µg/mL) towards B. subtilus and S. epidermis compared to the reference antibiotics Amoxicillin (MIC = 64 µg/mL) and Ampicillin (MIC = 78 µg/mL).Moreover, all these compounds were found to be more active than the standard drugs Amoxicillin (MIC = 64 µg/mL) and Ampicillin (MIC = 78 µg/mL) against the screened Gram-negative bacteria with MIC values ranging from 64 to 125 µg/mL. For the series 9a–h, results presented in Table 4 indicate that the most performant derivatives exhibit an inferior activity compared to the spirocompounds of series 4. They show only a poor to moderate activity with MIC values ranging from 64 to 250 µg/mL. It is noteworthy that compounds 4a, 4b, 4d, and 4e, with MIC value of 64 µg/mL, perform with an excellent activity, up to four times higher towards P. aeruginosa than Amoxicillin (MIC = 256 µg/mL).The compounds 9d and 9f–h (MIC = 64 µg/mL) are more potent than Amoxicillin (MIC = 78 µg/mL) against S. epidermis and exhibit an activity against B. subtilus similar to that of the reference antibiotic Ampicillin (MIC = 64 µg/mL). The results of antifungal screening (Table 5) reveal that compound 9e was the only one exhibiting an activity comparable to that of reference antifungal Amphotericin B (MIC = 500 µg/mL) against the tested fungi.The other synthesized compounds of series 4 were found to be more active than the standard drug with MICs ranging between 32–500 µg/mL. Notably, compound 4d was found to be the most potent antifungal agent (MIC = 32 µg/mL) towards C. krusei and C. glabrata compared to the reference antifungal Amphotericin B (MIC =500 µg/mL).In contrast, compounds 9a–h were revealed to be less active than their analogues, 4, against the three screened yeasts with MIC values ranging from of 125–500 µg/mL.
2.2.2. Structure-Activity Relationship (SAR)
Another aim of these studies was to analyze qualitatively and quantitatively the structure-activity relationships (SAR) of our dispiropyrrolidines. As described in the literature, spiropyrrolidine derivatives are known to be excellent antimicrobial agents [8,62,63]. For the first screening of the in vivo antimicrobial results, it is clearly outlined that the presence of a thiochromanone moiety in spiropyrrolidines 4a–e, independently of the nature of the para-substituent attached on the phenyl ring, greatly influences the inhibitory activity against all of the tested bacterial and fungal strains in comparison to their analogues containing a chromanone moiety 9a–h. These results conclude the importance of the introduction of the thiochromanone moiety to improve the activity of spiropyrrolidines. Compounds 4a–e especially feature a broad-spectrum bacterial and antifungal action and exerted the more potent antimicrobial effect as compared to the rest of the synthesized compounds.Concerning the series 9a–d, bearing an oxindole core, it is interesting to note that the introduction of an electron-withdrawing bromine substituent (9d) at the para position of the aromatic cycle of the chromanone, caused an increase of the antibacterial activity towards B. subtilus, S. epidermis, E.coli, K. pneumoniae, and P. aeruginosa. In contrast, the introduction of an electron-donating group, such as a methyl (9b) or methoxy group (9c), decreased the activity potency or retained it against the same bacteria.Among the spiropyrrolidines incorporating an acenaphthylene-1(2H)-one core 9e–h, compound 9e, with an unsubstituted phenyl ring, is the most active one against K. pneumoniae and P. aeruginosa. Compounds 9f (4-MePh), 9g (4-MeOPh), and 9h (4-BrPh) are found to exert a superior antibacterial effect against B. subtilus and S. epidermis. The inductive or mesomeric electronic effects exerted by these groups allow, probably, each bacterium to interact in its own specific way with the corresponding compounds.According to the antifungal activity results and based on the study of the structure-activity relationships (SAR), we noticed that compounds containing a thiochromanone ring, 4a–e, display significant antifungal activity against all the tested fungal strains, superior to those of the chromanone-grafted spiropyrrolidines, 9a–h, and the reference antifungal Amphotericin B. In addition, the compound 4d is the most potent one against C. krusei and C. glabrata. This finding shows that the presence of a bromine atom in para position of the phenyl group of chromanone seems to be responsible for the enhanced activity.
2.3. DFT Computational Studies
2.3.1. Optimized Molecular Structures and HOMO-LUMO Energies
The ground state geometric optimization of spiropyrrolidines 4a–f and 9a–h was carried out in the gas phase with the ωB97xd functional [58] and the standard 6–31g(d,p) basis set in Gaussian 16 package [64]. The simulated electron densities distribution of the FMOs, including both HOMO and LUMO orbitals of all compounds, are indicated in Table S2 (see Supporting Information). Frontier molecular orbital energies and chemical reactivity descriptor values of the optimized geometries are summarized in Table 6 and Table 7.
Table 6
The HOMO, LUMO energies and molecular descriptors of compounds 4a–f.
Compound
EHOMO
ELUMO
ΔEgap
ղ
μ
ω
(eV)
(eV)
4a
−7.463
−5.502
1.961
0.981
−6.483
21.422
4b
−7.462
−5.501
1.961
0.981
−6.482
21.415
4c
−7.459
−5.494
1.965
0.983
−6.477
21.338
4d
−7.462
−5.507
1.955
0.978
−6.485
21.501
4e
−7.497
−5.860
1.637
0.819
−6.679
27.234
4f
−7.496
−5.860
1.636
0.818
−6.678
27.259
Table 7
The HOMO, LUMO energies and molecular descriptors of compounds 9a–h.
Compound
EHOMO
ELUMO
ΔEgap
ղ
μ
ω
(eV)
(eV)
9a
−8.386
−5.499
2.887
1.444
−6.943
16.692
9b
−8.387
−5.497
2.890
1.445
−6.942
16.675
9c
−8.380
−5.492
2.888
1.444
−6.936
16.658
9d
−8.387
−5.499
2.888
1.444
−6.943
16.692
9e
−8.449
−5.875
2.574
1.287
−7.162
19.928
9f
−8.447
−5.874
2.573
1.287
−7.161
19.922
9g
−8.426
−5.874
2.552
1.276
−7.150
20.032
9h
−8.448
−5.878
2.570
1.285
−7.163
19.964
The molecular frontier orbitals provide important clues to the chemical reactivity of molecules. The FMOs energy gap describes the charge transfer character from HOMO to LUMO within the studied compounds. The HOMO energy represents the tendency of a molecule to donate electrons, while the LUMO energy indicates the ability to accept electrons [65].Thus, increasing the HOMO energy and decreasing the LUMO energy of the ligand lead to large intermolecular interactions with the LUMO and the HOMO of the receptor. The extent of these stabilizing interactions between ligand-receptor interacting orbitals correlates with the HOMO-LUMO energy gap. The softest molecules are determined with small energy gap values (ΔEgap = ELUMO − EHOMO) and are also identified as more stable and, hence, more reactive, according to Pearson’s “Hard and Soft Acids and Bases” principle [66].As shown from Table 4, the HOMO and LUMO electron-density of spiropyrrolidines 4a–f is delocalized on the thiochromanone ring except in 4e and 4f, where the electron-density of the LUMO, which indicates the electrophilic attack site, is mainly distributed over the acenaphtene nucleus. For the designed compounds 9a–f, the HOMO orbital was found to be distributed all over the pyrrolidine ring and isatin or acenaphtenequinone nucleus. However, the LUMO orbitals are dispersed over the chromanone ring in the series of compounds 9a–d and on the acenaphtenequinone ring in compounds 9e–h. Interestingly, the spiropyrrolidines 4a–f bearing a thiochroman-4-one moiety show smaller energy gap values, ranging between 1.636 and 1.965 eV, compared to spiropyrrolidine derivatives 9a–h featuring a chroman-4-one scaffold (ΔEgap in the range of 2.552 and 2.890 eV). Consequently, compounds 4a–f are the most reactive, compared to compounds 9a–h, which show the lowest reactivity (the most chemically stable). These results may explain the higher biological activity of compounds 4a–f, providing a good match with the experimental antifungal and antibacterial data.The reactivity descriptors, based on the analysis of the electronic chemical potential (μ), the chemical hardness (η), and the electrophilicity index (ω), provide useful insight into the chemical reactivity and stability of the molecules. The chemical potential describes the electron transfer capacity that occurs in the molecule in the fundamental state and the propensity of electrons to escape from an equilibrium system, whereas the chemical hardness indicates the resistance to charge transfer [67,68]. As shown in Table 4, compounds 4a–f showed higher chemical potential (ranging from −6.477 to −6.679 eV) and lower hardness values (in the range of 0.818 to 0.983 eV). Consequently, the spiropyrrolidines 4a–f are softer and more reactive than the chromanone-grafted spiropyrrolidines 9a–h. The biological activity of compounds 4a–f can, thus, be affected by the strong charge transfer interaction. The electrophilicity (ω), which defines the tendency of molecules to attract electrons [69], is also superior in the series of spiropyrrolidines 4a–d. Hence, the bioactivity of these compounds may be explained by the ability of the biological target to receive electrons from neighboring molecules, which may be important for stabilization of the active site.
2.3.2. Molecular Electrostatic Potential (MEP)
The MEP has emerged as a powerful approach in drug design and molecular biology to understand the intermolecular interactions between molecules and their biological receptors (proteins, enzymes). The MEP is applied to explore the chemical affinity and most likely binding modes with a molecular receptor, thus providing insight into the structure activity relationship of bioactive compounds. MEP-mapped surfaces inform on the chemical reactivity of the molecules and provide a reliable description of the charge distributions in a pictorial form. They also indicate the most likely electrophilic and nucleophilic reactive sites, polarization of molecules, as well as hydrogen bonding interactions [70,71,72,73].The sites related to the polar and nonpolar regions of the molecule are visualized by color variations. Electrophilic regions are indicated by blue coloration (electron-deficient regions), while nucleophilic regions are shown in red color (electron-rich regions). The electrostatic potential maps of the selected compounds 4a, 4e, 9c, 9d, 9e, and 9f are shown in Figure 7.
Figure 7
Molecular electrostatic potential (MEP) formed by mapping of the total density over the electrostatic potential for compounds 4a, 4e, 4c, 9d, 9e, and 9f.
2.4. Drug-Likeness Analysis
Analyzing the physicochemical properties of the developed drug hits is a crucial step to analyze and determine their drug-likeness potential. In this regard, quantitative structure-activity relationships (QSAR) descriptors, which can be calculated by means of a variety of software packages, may be a very useful tool. In this work, we evaluated the drug-likeness potential of the six selected compounds (4a, 4d, 4e, 9c, 9d, and 9e) through computing various QSAR descriptors, according to the Lipinski rule of five [74] and the Veber’s parameters (Table 8) [75].
Table 8
Results of Lipinski and Veber rules for the tested compounds.
Compound
4d
9d
9e
9c
4a
4e
lip_acc
6
7
6
8
6
5
lip_don
3
3
2
3
3
2
lip_druglike
1
1
1
1
1
0
lip_violation
1
1
1
0
0
2
logP(o/w)
4.46
3.83
4.81
2.99
3.67
5.75
Weight
564.48
548.41
504.56
499.54
485.58
534.66
TPSA
89.08
98.31
86.28
107.54
89.08
77.05
logS
−8.03
−7.06
−7.95
−6.02
−6.94
−9.4
vol
478.5
464.88
476
464.88
451.25
507.5
b_rotN
4
4
4
5
4
4
mr
14.37
13.8
14.32
13.68
13.61
15.34
a_hyd
27
26
29
26
26
31
Molecular weights of the compounds range from 485.58 to 564.48 Da. Derivatives 4a and 9c, having molecular weights less than 500 Da, obey the first Lipinski rule for effective and safe drug delivery. Lipinski’s second rule stipulates that drug-like compounds should not possess more than five hydrogen bond donating groups. All of the six compounds comply with this rule. Hydrogen bond accepting groups in all compounds are in the range of 5–8, which is less than ten, thus also meeting Lipinski’s third rule. The determined logP value was 2.99–5.75. All compounds, therefore, satisfy Lipinski’s rule regarding logP with a value inferior to five. An exception is compound 4e, whose logP value exceeds five. Compounds 4a and 9c showed zero Lipinski violations, compounds 4d, 9d, and 9e violated one rule, while compound 4e violated two rules and, thus, failed to show a drug-likeness potential (Table 8). Veber’s parameters, such as molar refractivity (mr) and total polar surface area (TPSA), demonstrate the polarizability of the compounds and suggest their oral bioavailability [75]. The mr values of the six compounds did not fall in the recommended range of 40–130 cm3/mol, however, this was offset by their TPSA value, which was not superior to 140 for any of the compounds. This suggests, for derivatives, an appreciable oral bioavailability. Besides, the number of rotatable bonds is four or five, which is inferior to 10. Except for compound 4e, which failed in drug-likeness parameters, some additional parameters, such as solubility, van der Waals volume, and number of hydrophobic atoms (a-hyd) were found to be in the acceptable range of −6 to −8, 451–478, and 26–31, respectively.
3. Materials and Methods
3.1. Apparatus and General Information
The 1H NMR spectra were recorded at 300 and 400 MHz using a Bruker Avance 300 or Bruker Avance III-400 machine (Rheinstetten, Germany). The chemical shifts were recorded in ppm relative to TMS and with the solvent resonance as the internal standard. Data were reported as follows: chemical shift, multiplicity (bs = broad singlet, s = singlet, d = doublet, t = triplet, m = multiplet), coupling constants (Hz), integration. 13C {1H} NMR data were collected at 75 or 100 MHz with complete proton decoupling with the solvent resonance as the internal standard. Elemental analyses were performed on a Perkin Elmer 2400 Series II Elemental CHNS analyzer (Waltham, MA, USA). Materials: thin-layer chromatography (TLC): TLC plates (Merck, silica gel 60 F254 0.2 mm, 200 × 200 nm) (Darmstadt, Germany); substances were detected using UV light at 254 nm.
3.2. General Procedure for Preparation of Cycloadducts 4 and 9
A mixture of 3-arylidenethiochroman-4-ones/3-arylidenecroman-4-ones (1 mmol), ethyl glycinate hydrochloride (1 mmol) 2, isatin 3 or acenaphthenequinone 3′ (1 mmol), and Et3N (1 mmol) in acetonitrile (5 mL) was heated under reflux for 2h. After completion of the reaction (TLC monitoring), the solvent was removed under vacuum. The residue was chromatographed on silica gel employing ethylacetate-cyclohexane (3:7 v/v) as eluent to obtain the pure products 4 and 9.
X-ray suitable crystals were grown in EtOH. The crystal structure determination was accomplished on a Bruker D8 Venture four-circle diffractometer using a PHOTON 100 CMOS (4a, 9c), or PHOTON II CPAD (4e), detector by Bruker AXS GmbH. X-ray radiation was generated by microfocus source IμS-Mo (λ = 0.71073 Å) by Incoatec GmbH with HELIOS mirror optics and a single-hole collimator by Bruker AXS GmbH. The crystals of 4a, 4e, and 9c were covered with an inert oil (perfluoropolyalkylether) and mounted on a MicroMount from MiTeGen. For the data collection, the programs APEX 3 Suite (v.2018.7-2), with the integrated programs SAINT (integration) and SADABS (adsorption correction) by BrukerAXS GmbH, were used. The processing and finalization of the crystal structures was done with the program Olex2 [76]. The crystal structure was solved with the ShelXT [77] structure solution program using Intrinsic Phasing and refined with the ShelXL [78] refinement package using Least Squares minimization. The non-hydrogen atoms were refined anisotropically. The C-bound and H atoms were placed in geometrically calculated positions and each was assigned a fixed isotropic displacement parameter based on a riding-model: C–H =0.95–0.99 Å with Uiso(H) = 1.5Ueq(CH3) and 1.2Ueq(CH2, CH) for other hydrogen atoms. The N-bound hydrogen atoms were located in the difference-Fourier-map and refined independently. The crystal and structure refinement data of 4a, 4e, and 9c are gathered in Table S3. Crystallographic data for the structures of 4a, 4e, and 9c have been deposited with the Cambridge Crystallographic Data Centre as supplementary publication number 2116000 (4a), 2116001 (4e), and 2115999 (9c). Copies of these data can be obtained, free of charge, on application to CCDC, 12 Union Road, Cambridge, CB2 IEZ, UK; fax: 144-(0)1223-336033 or e-mail: deposit@ccdc.cam.ac.uk.
3.4. Antimicrobial Activity
3.4.1. Determination of the Antibacterial and Antifungal Activity
Various concentrations of the screened compounds have been used to determine their antimicrobial activity. Both cultures have been set to yield approximately 1 × 106 CFUmL−1 of bacteria or yeast overnight. The minimum inhibitory concentrations (MICs) have been determined on brain heart infusion (BHI) agar plates (Bio-Rad, Marne la Coquette, France) by a standard method (Edziri et al., 2012). One milliliter of each of the extracts, previously dissolved in 10% dimethyl sulfoxide (DMSO), have been mixed for each concentration with 19 mL of BHI agar at 40 °C and poured over Petri dishes. The resulting DMSO concentration was approximately 0.5%. Plates containing only medium, or medium with 0.5% DMSO, were used as controls to ensure that DMSO did not affect the growth. Standard antibiotics Amoxicillin (AMX) and Ampicillin (AMP), and the antifungal Amphotericin B were used to control the sensitivity of the tested microorganism. After 18 h of incubation at 37 °C, the MIC was defined as the lowest concentration of the extract inhibiting the visible growth of each microorganism. Each test was carried out in triplicate.
3.4.2. Microorganisms
The antibacterial potential of the dispiropyrrolidine derivatives has been investigated against nine microorganisms, namely the Gram-positive Staphylococcus epidermidis (CI1232), Bacillus subtilis (ATCC 6633), Staphylococcus aureus (ATCC 29213), Staphylococcus aureus (ATCC 25923), and Enterococcus faecalis (ATCC 29212) and the Gram-negative rods Escherichia coli (ATCC 25922), Klebsiella pneumoniae (ATCC 4352), Salmonella enterica 800390, and Pseudomonas aeruginosa (ATCC 9023). For the evaluation of the antifungal activity, Candida glabrata (ATCC 90030), Candida albicans (ATCC 90028), and Candida krusei (ATCC 6258) have been used. As references, to evaluate and to compare the potency of the tested compounds under the same conditions, the antibiotics Amoxicillin (AMX) and Ampicillin (AMP) were chosen as antibacterial agents. Macrocyclic Amphotericin B was used as antifungal reference (both of the previously identified bacterial and fungal strains were obtained from the laboratory of Transmissible Diseases and Biologically Active Substances, Faculty of Pharmacy, Tunisia).
3.5. Computational Methods
DFT calculations were performed using the Gaussian 16 suite of programs [64]. All geometry optimizations and frequency analyses were carried out using the ωB97xd functional with the standard 6-31G(d,p) basis set [58]. Vibrational analyses were performed to identify the stationary points as either local minima, or transition states (with one imaginary vibrational mode) on the potential energy surface. Solvent effects were considered using the polarizable continuum model (PCM) with acetonitrile as the solvent medium at a temperature of 355 K [59,60]. The frontier molecular orbitals analysis (FMOs) and molecular electrostatic potential (MEP) surfaces for compounds 4a–f and 9a–h were calculated at the same level of theory and visualized with the graphical interface Gauss View 6 [79] (Gaussian Inc., Wallingford, CT, USA). The chemical reactivity descriptors [80,81,82] (electronic chemical potential (μ), electronegativity (χ), chemical hardness (η), and electrophilicity index ω) are calculated as follows: μ = (EH + EL)/2; η = (EL− EH)/2; ω = (μ2/2η).
4. Conclusions
A one-pot three-component 1,3-dipolar cycloaddition reaction of a stabilized azomethine ylide, generated in situ by thermal [1,2]-prototropy of the corresponding iminoesters with (E)-3-arylidene-thiochroman-4-one/chroman-4-one as dipolarophiles afforded the dispiro-oxindole-pyrrolidine-thiochroman-4-ones, 4a–f and dispiro-oxindole-pyrrolidine-chroman-4-ones, 9a–h in satisfactory yields along with high regio- and stereoselectivity. DFT calculations correctly predicted the observed experimental results, evidencing that the reaction is favorable according to kinetic and thermodynamic parameters. The antimicrobial activity of the novel N-heterocyclic compounds was examined against different pathogenic bacteria and fungi. According to the antimicrobial activity results, and on the basis of a structure-activity relationships (SAR) study, it follows that compounds 4a–d containing a thiochromanone ring display a significant antimicrobial activity against most of the tested bacterial and fungal strains, superior than those of chromanone-grafted spiropyrrolidines 9a–h. Even the activity of reference antibiotics is outperformed in some cases. Theoretical chemical reactivity, taking into account the FMOs and MEP maps, predicted a relatively higher chemical reactivity and softer character of spiropyrrolidines 4a–f, which seems to have a prominent effect on their higher antimicrobial activity compared to the 9a–h derivatives. The predicted absorption, distribution, metabolism, and excretion (ADME) profiles of some derivatives were in line with Lipinski rules. These promising results encourage further study of the antimicrobial activity of these derivatives by measuring their toxicity. Advanced docking studies are in progress and will be published elsewhere.
Authors: Leyla Yurttaş; Halide Edip Temel; Mehmet Onur Aksoy; Emre Fatih Bülbül; Gülşen Akalin Çiftçi Journal: Drug Dev Res Date: 2021-09-17 Impact factor: 4.360
Authors: Dipakkumar D Chudasama; Manan S Patel; Jaydeepkumar N Parekh; Harsh C Patel; Chetan V Rajput; Navin P Chikhaliya; Kesur R Ram Journal: Mol Divers Date: 2022-08-01 Impact factor: 3.364