Mohamed Elagawanyab,
Lamees Hegazya,
Feng Caoc,
Maureen J. Donlind,
Nigam Rathe,
John Tavisf and
Bahaa Elgendy*ag
aDepartments of Pharmacology and Physiology, Saint Louis University School of Medicine, St. Louis, MO, USA
bDepartment of Pharmaceutical Chemistry, Faculty of Pharmacy, Damanhour University, Damanhour, Egypt
cJohn Cochran Division, Department of Veteran's Affairs Medical Center, 915 North Grand Blvd., St. Louis, MO 63106, USA
dEdward A. Doisy Department of Biochemistry and Molecular Biology, Saint Louis University School of Medicine, St. Louis, Missouri, USA
eDepartment of Chemistry and Biochemistry, Center for Nanoscience, University of Missouri – St. Louis, One University Boulevard, St. Louis, USA
fDepartment of Molecular Microbiology and Immunology, The Saint Louis University Liver Center, Saint Louis University School of Medicine, St. Louis, Missouri, USA
gChemistry Department, Faculty of Science, Benha University, Benha 13518, Egypt. E-mail: belgendy@fsc.bu.edu.eg
First published on 23rd August 2018
We have synthesized and separated tosylated thujaplicin isomers for the first time, and elucidated their structures using 1D, 2D-NMR techniques and X-ray crystallography. The tosylated isomers were used to synthesize 4-isopropyl–thiotropolone and 6-isopropyl–thiotropolone in a regioselective manner. 1H and 13C Chemical shifts of synthesized isomers were fully assigned using several NMR experiments, and their isotropic magnetic shielding was calculated using the GIAO (Gauge Including Atomic Orbitals) method and the B3LYP def2-TZVPP level of theory. The calculated chemical shift values were in a good agreement with the experimental results. The biological activity of all synthesized compounds was evaluated against the fungal pathogen Cryptococcus neoformans and four different bacterial strains (Staphylococcus aureus (ATCC 29213), E. coli (ATCC 35218), Acinetobacter baumannii and Pseudomonas aeruginosa (ATCC 27853)). 4-Isopropyl–thiotropolone was found to inhibit Staphylococcus aureus in a low micro molar range and exhibit good therapeutic index and ADME properties. This compound can be used for future lead optimization to design inhibitors against Staphylococcus aureus (ATCC 29213).
Thiotropolones have not been widely explored compared to structurally related hydroxytropolones. Some examples exist in literature that shows anticancer activity such as thiocolchicin (3), which can act as tubulin12 or topoisomerase II inhibitor.13 Thiotropocin (4) is another structurally related antibiotic that was isolated from Pseudomonas sp. CB-104 in 1984.14,15 This antibiotic was found to exhibit antibacterial, antifungal, and antiprotozoal activities in vitro. Recently, Donlin and co-workers4 reported that thiotropolone 5 can inhibit the growth of the fungus Cryptococcus neoformans with a minimal 80% inhibitory concentration (MIC80) of 0.25 μM. In our effort to design and synthesize new thiotropolones to treat pathogenic fungi and multidrug resistant bacteria, we used the natural product β-thujaplicin (2) as starting point (Fig. 1).
Here, we report for the first time the isolation of the intractable mixture of tosylated thujaplicines, 8 and 9, their full NMR assignments and the X-ray structure of 8. Moreover, we report the synthesis, full NMR chemical shift assignments of 4-isopropyl and 6-isopropyl-2-mercaptocyclohepta-2,4,6-trien-1-one, 10 and 11, and their anti-Cryptococcus neoformans and antibacterial activity. Moreover, we have evaluated the ADME (adsorption, distribution, metabolism, excretion and toxicity) properties of tested compounds.
Tosylated thujaplicines, 8 and 9, are important synthetic intermediates in the synthesis of many biologically active tropolones. For example, they were used as an intermediate in the preparation of copper aminotropones, which were found to be more effective than a commercial toothpaste formulation in inhibiting plaque formation.18 The inability to separate 8 and 9 and their use as a mixture of isomers complicates the purification of the desired ligands. Obtaining such mixture of isomers is undesirable, and selective tosylation or facile separation of isomers is crucial from the synthetic point of view.18
Hinokitiol (β-thujaplicin) (2) was treated with p-toluenesulfonyl chloride in presence of pyridine in dichloromethane (DCM) at room temperature to give a mixture of 5-isopropyl-7-oxocyclohepta-1,3,5-trien-1-yl-4-methylbenzenesulfonate (8) and 3-isopropyl-7-oxocyclohepta-1,3,5-trien-1-yl-4-methylbenzenesulfonate (9) in ≈1:1 ratio. The isolation of this mixture was reported to be intractable.19 To obtain 9 (the isomer of the natural product) exclusively, it has to be through multistep synthesis. The first step involves iodination of 2 in position 7 followed by tosylation of hydroxyl group at position 2. The final step is to remove the iodine by hydrogenation using Pd/C (10%) in presence of anhydrous sodium acetate.19 Nevertheless, we were able to separate the two isomers using a Reveleris X2 flash chromatography system. The structures of 8 and 9 (Scheme 1) were elucidated using 1D and 2D NMR spectroscopy and single crystal X-ray diffraction.
Compound 8 was obtained as colorless sharp needles. Its molecular formula C17H18O4S was determined by HRESIMS data (m/z 341.0816 [M + Na]+, calc. 341.0818). The analysis of 1H, 13C, COSY, HSQC, and HMBC spectra revealed the presence of a singlet proton H-7 (δC/δH 137.7/7.05, C-7). This proton shows two bond correlations with C-1 (δC 179.1) and C-6 (δC 157.3) and strong three bond correlations with C-2 (δC 154.5) and C-8 (δC 38.3). Furthermore, the one methine septet H-8 (δC/δH 38.3/2.73, C-8) shows two bond correlations to C-6 and three bond correlations to C-7. The two-methyl doublet H3-9 and H3-10 (δC/δH 22.6/1.17, C-9 and C-10) shows a strong three bond correlation to C-6. H-4 (δC/δH 129.1/6.88, C-4) shows strong three bond correlations to C-2 and C-6, which further confirms the location of tosyl and isopropyl groups (Table 1).
Atom | δExp 8 | δCalc 8 | δExp 9 | δCalc 9 |
---|---|---|---|---|
a 13C assignments were based on HSQC and HMBC spectra. 1H assignments were based on 1D-1H NMR and COSY spectra. | ||||
C-1 | 179.1 | 188.9 | 178.9 | 189.1 |
C-2 | 154.5 | 173.1 | 154.6 | 171.4 |
C-3 | 129.2 | 140.9 | 131.7 | 145.6 |
C-4 | 129.6 | 142.1 | 152.4 | 167.8 |
C-5 | 136.6 | 143.4 | 130.5 | 138.2 |
C-6 | 157.3 | 167.1 | 137.0 | 144.5 |
C-7 | 137.7 | 155.5 | 138.9 | 152.4 |
C-8 | 38.3 | 51.2 | 38.0 | 48.1 |
C-9 and C-10 | 22.6 | 28.9 | 22.7 | 28.3 |
C-11 | 133.4 | 151.9 | 133.4 | 151.7 |
C-12 and C-16 | 128.4 | 141.4 | 128.5 | 141.4 |
C-13 and C-15 | 129.2 | 138.5 | 129.6 | 138.7 |
C-14 | 145.3 | 157.3 | 145.6 | 157.2 |
C-17 | 21.7 | 27.3 | 21.8 | 27.4 |
H-3 | 7.29 | 8.12 | 7.36 | 8.31 |
H-4 | 6.88 | 7.69 | — | — |
H-5 | 6.98 | 7.73 | 6.93 | 7.68 |
H-6 | — | — | 7.14 | 7.94 |
H-7 | 7.05 | 8.23 | 7.02 | 8.044 |
H-8 | 2.73 | 3.47 | 2.80 | 3.64 |
H-9 and H-10 | 1.17 | 2.01 | 1.21 | 2.05 |
H-12 and H-16 | 7.89 | 10.15 | 7.90 | 10.07 |
H-13 and H-15 | 7.31 | 8.27 | 7.33 | 8.28 |
H-17 | 2.42 | 3.35 | 2.43 | 3.3445 |
We have calculated both 1H and 13C isotropic chemical shielding of 8 at the B3LYP/def2-TZVPP level of theory. Correlation plots between experimental and calculated chemical shifts shows very good correlations with correlation coefficients (R2) of 0.9946 and 0.9755 for 13C and 1H NMR chemical shifts, respectively (Fig. S50†).
Compound 9 was obtained as a white powder. Its molecular formula C17H18O4S was determined by HRESIMS data (m/z 341.0818 [M + Na]+, calc. 341.0815). The analysis of 1H, 13C, COSY, HSQC, and HMBC spectra revealed the presence of a singlet proton (δC/δH 131.7/7.36, C-3). This proton shows two bond correlations with C-2 (δC 154.5) and C-4 (δC 152.4) and strong three bond correlations with C-1 (δC 178.9) and C-8 (δC 38.0). Furthermore, the one methine septet H-8 (δC/δH 38.0/2.80, C-8) shows two bond correlations to C-4 (δC 152.4) and three bond correlations to C-3. The two-methyl doublet H3-9 and H3-10 (δC/δH 22.7/1.21, C-9 and C-10) shows a strong three bond correlation to C-6. H-6 (δC/δH 137.0/7.14, C-6) shows strong three bond correlations to C-1 and C-4, which further confirms the location of carbonyl and isopropyl groups. Besides, H-7 (δC/δH 138.9/7.02, C-7) shows a strong three bond correlation to C-2, which confirms the position of tosyl group (Table 1).
We have calculated both 1H and 13C isotropic chemical shielding of 9 at the B3LYP/def2-TZVPP level of theory. Correlation plots between experimental and calculated chemical shifts shows very good correlations with R2 of 0.9962 and 0.9789 for 13C and 1H NMR chemical shifts, respectively (Fig. S51†).
Our NMR assignments were further confirmed by obtaining an X-ray crystal structure of compound 8 (Fig. 3). The X-ray shows that the tosyl group is in position 2 while isopropyl group is in position 6 of the tropolone ring. Compound 8 was crystallized from ethyl acetate (0.2 mL), ether (0.5 mL), and hexanes (10 mL) by solvent evaporation. Crystallographic data is given in Table S1.† The crystals of compound 8 crystallize in the monoclinic space group P21/c with four molecules in the unit cell (Z = 4) of dimensions a = 5.62, b = 28.31, c = 9.51 Å and β = 92.4°. V = 1512.31(11) Å3.
The 7-membered ring is disordered over 2 orientations (80:20%). The disorder was modelled with partial occupancy atoms and it indicates that the 7-membered ring was “breathing” by moving from less hindered to more hindered confirmation which facilitate the migration of tosyl group in substitution reactions. The tosyl group in 8 and 9 were shown to be fluxional and undergo thermally-induced shift between the two thujaplicin oxygen atoms. Molecular mechanics calculations by the MMFF94 force field suggests that this process takes place most likely through bipyramidal intermediate.20
6-Isopropyl–thiotropolone (10) and 4-isopropyl–thiotropolone (11) were readily accessible in good yields by reacting 8 and 9 with sodium hydrosulfide in ethanol at room temperature (Scheme 2). The structures of 10 and 11 were fully assigned using 1D and 2D NMR spectroscopy.
2-Hydroxy-2,4,6-cycloheptatrien-1-one (tropolone 12) possesses hydroxyl group alpha to the carbonyl group and is expected to exhibit keto–enol tautomerism. Interestingly, this compound exists almost exclusively in the enol form (12A), which possess resonance stabilization because of its aromaticity. Density functional theory calculations of the keto–enol equilibrium constants of 2-hydroxy-2,4,6-cycloheptatrien-1-one (12A) and 3,5- and 3,6-cycloheptadiene-1,2-dione (12B and 12C) showed that the equilibrium is tilted in favour of the enol (12A) (Fig. 4).21
The corresponding thiotropolones have seven possible tautomers. The tautomeric equilibrium of all tautomers of thiotropolone was studied using Hartree Fock (HF) and density functional theory (B3LYP) at different levels of theory.22 Calculating electronic energy and Gibbs free energy showed that the enol forms were more stable than the keto forms because of its aromaticity. This observation was supported by calculating the nucleus independent chemical shifts of NMR chemical shifts of the studied thiotropolones.22 Compound 13A (2-hydroxy-2,4,6-cycloheptatriene-1-thione) was the most stable followed by compound 13B (2-mercapto-2,4,6-cycloheptatrien-1-one). The low energy of these two isomers is attributed to the favourable interactions between H5-S1 in 13A and H5-O1 in 13B.22
Nozoe et al.23 was the first to synthesize 2-mercaptotropone and it was suggested to exist in two tautomeric forms. Sulfur-substituted derivatives were obtained upon alkylation and acylation of this compound, which suggests that 2-mercaptotropone (13B) is the prevailing tautomer. However, studying the physicochemical properties of 2-mercaptotropone and its structural analogues 6-isopropyl-2-mercaptotropone (10) and 4-isopropyl-2-mercaptotropone (11) in depth suggested that the major tautomer is the 2-hydroxyl tropothione (13A).24
Compound 10 was obtained as a dark red oil. Its molecular formula (C10H11OS)2Na+was determined by HRESIMS data (m/z 381.0951 [M2 + Na]+, calc. 381.0953). The 1H NMR showed that this compound exist in CDCl3 as two tautomers 10A (≈87%) and 10B (≈13%). We were able to fully assign the chemical shifts of the major tautomer 10A (Table 2). The analysis of 1H, 13C, COSY, HSQC, and HMBC spectra revealed the presence of a singlet proton H-3 (δC/δH 119.6/7.45, C-3). This proton shows two bond correlations with C-2 (δC 172.5) and C-4 (δC 161.2) and strong three bond correlations with C-1 (δC 183.4), C-5 (δC 132.1) and C-8 (δC 38.8). Furthermore, the one methine septet H-8 (δC/δH 38.8/2.91, C-8) shows two bond correlations to C-4, C-9 and C-10 (δC 23.4) and three bond correlations to C-3 and C-5. The two-methyl doublet H3-9 and H3-10 (δC/δH 23.4/1.27, C-9 and C-10) shows two bond correlations to C-8 and a strong three bond correlation to C-4. H-5 (δC/δH 132.1/7.21, C-5) shows two bond correlations to C-6, and strong three bond correlations to C-3, C-7 and C-8. H-6 (δC/δH 133.8/7.14, C-3) shows strong three bond correlations to C-1 and C-4. H-7 (δC/δH 143.0/8.45, C-7) shows strong three bond correlations to C-2 and C-5 (Fig. 5).
Atom | δExp 10A | δCalc 10A | δCalc 10B | Atom | δExp 10A | δCalc 10A | δCalc 10B |
---|---|---|---|---|---|---|---|
C-1 | 183.4 | 210.3 | 180.0 | H-3 | 7.45 | 6.71 | 7.48 |
C-2 | 172.5 | 179.0 | 185.9 | H-4 | — | — | — |
C-3 | 119.6 | 118.2 | 141.2 | H-5 | 7.21 | 6.98 | 7.07 |
C-4 | 161.2 | 164.9 | 164.1 | H-6 | 7.14 | 6.88 | 7.18 |
C-5 | 132.1 | 132.6 | 133.9 | H-7 | 8.45 | 8.98 | 8.08 |
C-6 | 133.8 | 130.9 | 136.8 | H-8 | 2.91 | 2.62 | 2.91 |
C-7 | 143.0 | 159.8 | 130.8 | H3-9 | 1.27 | 1.46 | 1.42 |
C-8 | 38.8 | 47.8 | 47.9 | H3-10 | 1.27 | 1.45 | 1.47 |
C-9 | 23.4 | 26.9 | 25.2 | OH or SH | 9.50 | 8.98 | 4.04 |
C-10 | 23.4 | 24.3 | 24.0 |
We have calculated both 1H and 13C isotropic chemical shielding of 10A and 10B at the B3LYP/def2-TZVPP level of theory. Correlation plots between experimental and calculated chemical shifts shows very good correlations in case of 10A with R2 of 0.9828 and 0.99 for 13C and 1H NMR chemical shifts, respectively (Fig. 6). The calculated chemical shift values of isomer 10B correlated to less extent with the experimental values (see Fig. S52†).
Fig. 6 (A) Experimental 13C chemical shifts (δ) of 10A (vs.) calculated. (B) Experimental 1H chemical shifts (δ) of 10A (vs.) calculated. |
Compound 11 was obtained as orange microcrystals. Its molecular formula (C10H11OS)2Na+ was determined by HRESIMS data (m/z 381.0950 [M2 + Na]+, calc. 381.0953). This compound exists in CDCl3 as a mixture of two tautomers 11A (≈96%) and 11B (≈4%). We were able to fully assign the chemical shifts of the major tautomer 11A (Table 3). The analysis of 1H, 13C, COSY, HSQC, and HMBC spectra revealed the presence of a singlet proton H-7 (δC/δH 144.0/8.62, C-7). This proton shows two bond correlations with C-1 (δC 184.1) and C-6 (δC 156.5) and strong three bond correlations with C-2 (δC 172.2) and C-5 (δC 133.0). Furthermore, the one methine septet H-8 (δC/δH 39.0/2.88, C-8) shows two bond correlations to C-6 (δC 156.5) and three bond correlations to C-7. The two-methyl doublet H3-9 and H3-10 (δC/δH 23.1/1.27, C-9 and C-10) shows two bond correlations to C-8 and a strong three bond correlation to C-6. H-5 (δC/δH 133.0/7.19, C-5) shows two bond correlations to C-6, strong three bond correlations to C-3 and C-8, and four bond correlations to C-2. H-3 (δC/δH 118.2/7.34, C-3) shows two bond correlations to C-2 and C-4 and strong three bond correlations to C-1 and C-5. H-4 (δC/δH 137.9/7.35, C-4) shows two bond correlations to C-5 and strong three bond correlations to C-2 and C-6 which further confirms the location of the isopropyl groups (Fig. 5).
Atom | δExp 11A | δCalc 11A | δCalc 11B | Atom | δExp 11A | δCalc 11A | δCalc 11B |
---|---|---|---|---|---|---|---|
C-1 | 184.1 | 202.0 | 181.5 | H-3 | 7.34 | 7.50 | 7.29 |
C-2 | 172.2 | 179.8 | 185.5 | H-4 | 7.35 | 7.54 | 7.51 |
C-3 | 118.2 | 117.1 | 138.3 | H-5 | 7.19 | 7.34 | 7.07 |
C-4 | 137.9 | 142.0 | 140.4 | H-6 | — | — | — |
C-5 | 133.0 | 134.3 | 131.1 | H-7 | 8.62 | 9.32 | 8.27 |
C-6 | 156.5 | 161.1 | 163.3 | H-8 | 2.88 | 2.94 | 2.91 |
C-7 | 144.0 | 152.8 | 135.2 | H3-9 | 1.27 | 1.47 | 1.51 |
C-8 | 39.0 | 47.8 | 47.7 | H3-10 | 1.27 | 1.51 | 1.44 |
C-9 | 23.1 | 26.4 | 24.0 | OH or SH | 9.47 | 11.46 | 4.10 |
C-10 | 23.1 | 25.4 | 28.1 |
1H and 13C isotropic chemical shielding of 11A and 11B were calculated at the B3LYP/def2-TZVPP level of theory (Table 3). Correlation plots between experimental and calculated chemical shifts shows very good correlations in the case of 11A, with R2 of 0.994 and 0.9861 for 13C and 1H NMR chemical shifts, respectively (Fig. 7). The calculated chemical shift values of isomer 11B correlated to less extent with the experimental values (see Fig. S53†).
Fig. 7 (A) Experimental 13C chemical shifts (δ) of 11A (vs.) calculated. (B) Experimental 1H chemical shifts (δ) of 11A (vs.) calculated. |
Comp. | Bacteria MIC80 (μM) | C. neoformans MIC80 | CC50 | ||||
---|---|---|---|---|---|---|---|
S. aureus (29213) | E. coli (35218) | A. baumannii (from a patient) | P. aeruginosa (27853) | CC50 (neutral red) | CC50 (MTS) | ||
24 | 66.7 | 44.4 | 51.2 | >100 | 21 (ref. 4) | na | 66.4 |
8 | >100 | >100 | >100 | >100 | 24 | 31.8 | 33.6 |
9 | 8.8 | >100 | >100 | >100 | 24 | 19.5 | 34.2 |
10 | >100 | >100 | >100 | >100 | 50 | 39.7 | 95 |
11 | 16 | >100 | >100 | >100 | 24 | 30.4 | 64.2 |
The mercapto derivative of 8, compound 10, has no antibacterial activity against all four tested bacteria. The activity of 10 against C. neoformans (MIC80 = 50 μM) is >2-fold less than the corresponding hydroxy tropolone 2. The mercapto derivative of 9, compound 11, which is the direct isomer of 2, showed promising antibacterial activity against S. aureus (MIC80 = 16 μM), but no activity against the other tested bacteria. Compound 11 showed similar activity as 2, 8, and 9 against C. neoformans (MIC80 = 24 μM) (Table 4).
We tested the cytotoxicity of the synthesized compounds in hepatoblastoma cells using an MTS cytotoxicity assay that measures mitochondrial function and a neutral red assay that measures lysosomal function and calculated the CC50, which is the concentration of inhibitor required to reduce cell viability 50% relative to untreated cells (Table 4). Compound 10 have the lowest toxicity among all synthesized compounds in both assays. The toxicity of compound 11 is relatively low, and comparable to the natural product β-thujaplicin (2).
Comp. | amol_MW | bQPlogPo/w | cHBD | dHBA | e% Human oral absorption | fPSA | gQplogS | hQPPCaco | i#Metab | jQplogBB | kQPPMDCK |
---|---|---|---|---|---|---|---|---|---|---|---|
a Molar weight in Daltons (130–725).b Logarithm of partitioning coefficient between n-octanol and water phases (range for 95% of drugs: −2 to 6).c Number of hydrogen bonds donors (0–6).d Number of hydrogen bond acceptors (2–20).e Predicted human oral absorption on a 0–100% scale, based on a multiple linear regression model (<25% low, >80% high).f Polar surface area (7–200).g Predicted aqueous solubility, logS. S in mol dm−3 (−6.5 to 0.5).h Predicted apparent Caco-2 cell permeability in nm s−1 as a model for the gut-blood barrier (<25 poor, >500 excellent).i Number of possible metabolic reactions (2–8).j Predicted brain/blood partition coefficient (−3 to 1.2).k Predicted apparent MDCK cell permeability in nm s−1 as a mimic for blood/brain barrier (<25 poor, >500 excellent). Qikrop predictions are for non-active transport. | |||||||||||
2 | 164.204 | 1.443 | 1 | 2.75 | 89.394 | 50.125 | −2.193 | 1040.346 | 2 | −0.396 | 516.318 |
8 | 318.387 | 2.1 | 0 | 6.5 | 92.748 | 68.778 | −2.742 | 975.826 | 2 | −0.615 | 490.424 |
9 | 318.387 | 2.216 | 0 | 6.5 | 94.321 | 68.77 | −2.952 | 1095.291 | 2 | −0.591 | 555.56 |
10A | 180.264 | 2.351 | 0.8 | 2.5 | 100 | 28.102 | −2.794 | 3286.638 | 2 | 0.207 | 4437.612 |
10B | 180.264 | 2.4 | 1 | 2.25 | 100 | 29.939 | −2.881 | 3010.296 | 2 | 0.168 | 3962.203 |
11A | 180.264 | 2.36 | 0.8 | 2.5 | 100 | 28.166 | −2.843 | 3222.618 | 2 | 0.195 | 4334.35 |
11B | 180.264 | 2.398 | 1 | 2.25 | 100 | 30.01 | −2.884 | 2998.094 | 2 | 0.166 | 3935.7 |
Crystal data and intensity data collection parameters are listed in Table 1S.† Structure solution and refinement were carried out using the SHELXTL-PLUS software package.40 The structure was solved and refined successfully in the monoclinic space group P21/c. Full matrix least-squares refinements were carried out by minimizing w(Fo2 − Fc2)2. The non-hydrogen atoms were refined anisotropically to convergence. The 7-membered ring is disordered, and the two orientations were refined with geometrical restraints to 82:18%. All hydrogen atoms were treated using appropriate riding model (AFIX m3). The final residual values and structure refinement parameters are listed in Table 1S.†
Complete listings of positional and isotropic displacement coefficients for hydrogen atoms, anisotropic displacement coefficients for the non-hydrogen atoms and other geometrical parameters are listed as ESI (Tables 2S to 7S†). Table of calculated and observed structure factors are available in electronic format.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 1856450. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c8ra06297h |
This journal is © The Royal Society of Chemistry 2018 |