Attila Mándia,
Jun Wu*b and
Tibor Kurtán
*a
aDepartment of Organic Chemistry, University of Debrecen, P. O. Box 400, 4002 Debrecen, Hungary. E-mail: kurtan.tibor@science.unideb.hu; Fax: +36-52-512-744
bSchool of Pharmaceutical Sciences, Southern Medical University, 1838 Guangzhou Avenue North, Guangzhou 510515, P. R. China. E-mail: wwujun2003@yahoo.com
First published on 1st September 2020
TDDFT-ECD calculations were utilized to explain the mirror image or different ECD spectra of previously reported homochiral natural products thaigranatins A–E and granatumin L, the simple comparison of which would result in a wrong stereochemical conclusion. The configurational assignment was confirmed independently and geometrical parameters of the chromophores governing the ECD spectra were identified in the structurally related natural products by analyzing the ECD spectra and geometries of the low-energy computed conformers obtained by different methods. Different conformations of the furan-2-yl-δ-lactone subunit were found responsible for the mirror image ECD spectra of the homochiral thaigranatins C–E. Two DFT 13C NMR chemical shift calculation methods and DP4+ analysis were performed on the C-6 epimers of thaigranatin D, which together with the ECD calculation, could determine the absolute configuration of C-6 as (R).
For instance, the pyridine alkaloids penipyridones C (D) and F (E) had the same chromophoric system and absolute configuration but E had an additional N-2-hydroxyethyl group, which resulted in near mirror image ECD spectra and oppositely signed specific rotation (Fig. 1).10 TDDFT-ECD calculations and analysis of the low-energy conformers revealed that in the presence of the achiral N-2-hydroxyethyl substituent, the additional 15-OH group formed an intramolecular hydrogen bond with the 7-OH, which changed the relative orientation of the phenyl and γ-pyridone chromophores.10 The replacement of oxygen heteroatom with sulfur in the benzene-condensed spiro derivatives F and G containing benzo[1,4]oxazin-3-one and benzo[1,4]thiazin-3-one chromophores resulted in significantly different ECD spectra with mirror image transitions in the high-wavelength region, although the computed conformers had very similar geometries (Fig. 1).9
Aflaquinolone I (H) and aniduquinolone A (I) shared a very similar substituted 4-aryl-3,4-dihydroquinolin-2(1H)-one chromophore but in the presence of the 5-OH group in the peri position of I, the C-4 phenyl group adopted an axial orientation, while in H, having a hydrogen atom in the peri position, the C-4 aryl group preferred the equatorial orientation in all the low-energy computed conformers (Fig. 1).6,7 Due to the different conformations of the condensed heterocycles, the ECD spectra of H and I were fairly different and comparison of ECDs could not be used to determine the AC of H.7 Psammaplysins A (J) and B (K), first representatives of the psammaplysin alkaloid family, differed in the additional C-19 benzylic hydroxyl group of K, which implied a C-19 chirality center in the side-chain of K (Fig. 1).14 Regardless this benzylic chirality center, the experimental ECD spectra of J and K were found near identical, which confirmed that the side-chain with an aryl ether chromophore had negligible contribution to the ECD spectra. Thus a truncated model compound could be utilized for the ECD calculations of J, which reduced the number of conformers and enabled the configurational assignment of J.14 The above examples may justify to explore the structural and chiroptical background for the different ECD profiles of structurally related derivatives by ECD calculations, which would allow us avoiding wrong configurational assignments for analogous derivatives.
Thaigranatins A–E (1–5), limonoid derivatives containing a pentacyclic skeleton with a furan-3-yl substituent at the δ-lactone ring E, have been recently isolated together with the previously described granatumin L15 (6) from the seeds of the Thai Xylocarpus granatum and reported by Ren et al. (Fig. 2).16
![]() | ||
Fig. 2 Structures of thaigranatins A–E (1–5) and granatumin L (6) isolated from Xylocarpus granatum.16 |
Related limonoids with the same skeleton have been isolated recently from Xylocarpus granatum,17–20 Xylocarpus moluccensis,21,22 Chisocheton ceramicus23 and Chisocheton erythrocarpus.24 Limonoids were suggested to have even potent inhibitory activity against SARS-CoV-2 by in silico docking studies.25 Since limonoids exhibit a wide range of bioactivities such as antiviral, anticancer activity and they are potential endogenous regulators in the nervous system, they have also attracted considerable synthetic interest.26
The AC of thaigranatin A (1) was elucidated by single crystal X-ray analysis, while those of 2 and 6 were determined on the basis of their similar ECD spectra to that of thaigranatin A (1).16 The experimental solution ECD spectra of 1 and 6 were quite similar, while that of 2 had a broad positive ECD band with a maximum at 217 nm, which were missing from those of 1 and 6. The ECD spectra of 3–5 were significantly different from the others, which did not enable to utilize the simple comparison of the ECD spectra for the configurational assignment. Moreover, the ECD spectrum of 5 was found near mirror image of those of 3 and 4, which is quite puzzling considering their analogous structure with minor structural differences (Fig. 3). The absolute configurations of 3–5 were proposed on the basis of the common biosynthetic origin with thaigranatin A (1).16 In this work, we utilized TDDFT-ECD calculations of 1–6 with different combinations of methods to explain the ECD spectral differences of 1–6 and the mirror image ECD spectra of 3–5. This approach was tested whether it is suitable to prove the proposed homochirality of 1–6 and thus it can provide an independent configurational assignment for 3–5. Furthermore, the AC of C-6 in the side-chain of 4 remained undetermined.16 In this work we applied DFT-NMR 13C chemical shift calculations of the C-6 epimers to elucidate the AC of this chirality center.27–29
Semi-synthetic derivatives obtained by the derivatization of granatumin L (6) showed potent human immunodeficiency virus 1 (HIV-1) and influenza A virus (IAV) inhibitory activities, which may also justify to explore relationship between stereochemistry and characteristic ECD features and to determine solution geometries and confirm the absolute configuration.16
In all the CAM-B3LYP/TZVP PCM/MeCN conformers of 1, the C-17 furan-3-yl group adopted an equatorial arrangement with syn coplanar orientation of the 17-H and C-21 atoms (ω17-H,C-17,C-20,C-21 = −10.2°) and the 14-H had axial orientation. The δ-lactone ring had a distorted half-chair conformation with the C-13 being out of plane, while the other five atoms located in nearly one plane. The geometry of the low-energy conformers differed in the relative arrangements of the C-5 and C-3 substituents by rotation along the C-6 to C-7 and C-32 to C-33 sigma bonds (Fig. 5). In the lowest-energy conformer (conf A), the C-32 carbonyl oxygen was syn coplanar with the C-33 methyl group (ωO,C-32,C-33,C-36 = −0.1°, s-trans conformation), while the C-5 had a gauche orientation with the C-7 carbonyl oxygen (ωC-5,C-6,C-7,
O = 42.2°). Conformers A and D differed only in the arrangements of the C-1 hydroxyl proton and they had identical computed ECD spectra (Fig. 6). Positive rotatory strength values were computed above 200 nm as shown in Fig. 4 but these were overcompensated by the negative ones. Conformers C and F had also positive computed transitions above 200 nm with larger intensity than those of conformers A and D manifesting in a positive CE (Fig. 6), and their geometries also differed only in the orientation of the 1-OH proton. Similarly to conformers A and D, conformers C and F had syn coplanar arrangement of the C-32 carbonyl oxygen (s-trans enone) and the C-33 methyl group but the C-5 had gauche orientation with the C-7 methoxy oxygen (ωC-5,C-6,C-7,
O = −137.4°). The positive CEs of conformers A, C, D, F were cancelled by the strong negative CEs of conformers B and E, in which the C-32 carbonyl oxygen was anti periplanar with the C-33 methyl group (ωO
,C-32,C-33,C-36 = −178.4°, s-cis enone). Thus the rotation along the C-32 to C-33 bond of the α,β-unsaturated ester moiety and around the C-6 to C-7 bond was mostly responsible for the differences in the ECD spectra of the low-energy conformers of 1.
![]() | ||
Fig. 5 Overlapped geometries and selected torsional angles of the six lowest-energy conformers of (1R,2S,3R,4S,5S,6R,9S,10R,13R,14S,17R)-1. Level of optimization: CAM-B3LYP/TZVP PCM/MeCN. |
Compared to compound 1, 2 did not contain the conjugating Δ33,34 double bond and the C-6 chirality center in the other side chain and the experimental ECD spectrum had a strong negative CE below 210 nm and a weaker positive CE at 217 nm with a positive plateau up to 260 nm. Compound 2 had 43 initial MMFF conformers, the gas-phase and solvent model re-optimizations of which yielded 6 and 10 low-energy conformers over 1% Boltzmann-population. The first five low-energy CAM-B3LYP/TZVP PCM/MeCN conformers had differences in the orientation of the C-3 ester group and the furan ring while the lactone ring showed only marginal differences (Fig. 7). The individual computed ECDs of these conformers were similar but that of conformer E showed a significant red shift due to the opposite orientation of the furan ring (Fig. 8 and S2†). This red shift of the minor conformers obviously contributes to the high-wavelength positive shoulder of the experimental ECD spectrum. Similarly to the conformers of 1, the δ-lactone ring had a distorted half-chair conformation in the low-energy conformers. The Boltzmann-weighted ECD spectra computed at all the applied combinations reproduced well the experimental spectrum including the positive CE and shoulder (Fig. 9), which confirmed the (1R,2S,3R,4S,5S,9S,10R,13R,14S,17R) absolute configuration of 2 and hence the homochirality with 1.16 The difference in the experimental ECD spectra of 1 and 2 can be attributed to the absence of the Δ33,34 double bond and the C-6 chirality center.
![]() | ||
Fig. 7 Overlapped geometries of the five lowest-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (1R,2S,3R,4S,5S,9S,10R,13R,14S,17R)-2. |
![]() | ||
Fig. 8 Experimental ECD spectrum of 2 in MeCN compared with the B3LYP/TZVP PCM/MeCN ECD spectra of the individual five lowest-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (1R,2S,3R,4S,5S,9S,10R,13R,14S,17R)-2. Experimental ECD spectrum of 2 was reproduced with permission from ref. 16. |
Compound 3 contained the same C-3 α,β-unsaturated ester moiety as 1 but it missed the C-6 chirality center due to the exchange of the 6-OH group with a hydrogen atom. Moreover, the Δ8,30 double bond of 1 was shifted to Δ8,14 position in 3 and C-30 became a chirality center due to the presence of the sec-hydroxyl group. Besides the intense negative CE below 211 nm, the experimental ECD spectrum of 3 had a positive CE at 220 nm accompanied with a flat positive shoulder having maximum at 286 nm.
Re-optimization of the initial 51 MMFF conformers of (1R,2S,3R,4S,5S,9S,10R,13R,17R,30R)-3 resulted in 15 and 24 low-energy conformers over 1% Boltzmann-population in gas-phase and solvent model calculations, respectively. In the first five low-energy CAM-B3LYP/TZVP PCM/MeCN conformers, there were only minor difference in orientations of the two side-chains and the lactone ring (Fig. 10). Similarly to the lowest-energy conformer of 1, the C-32 carbonyl oxygen of conformers A–E was syn coplanar with the C-33 methyl group (ωO,C-32,C-33,C-36 = 34.8–37.8° range for confs A–E, s-trans enone) and the δ-lactone ring had a distorted half-chair conformation even with the sp2-hybridized C-14. For the first five low-energy CAM-B3LYP/TZVP PCM/MeCN conformers, the computed ECD spectra had similar profile for the main transitions (Fig. 11). The different ECD profiles of compounds 1 and 3 derived from the different distributions of conformers differing in the C-3 side-chain conformation.
![]() | ||
Fig. 10 Overlapped geometries of the five lowest-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (1R,2S,3R,4S,5S,9S,10R,13R,17R,30R)-3 (left: conformers A and D; right: conformers B, C and E). |
Boltzmann-weighted ECD spectra obtained at all combinations of the applied levels reproduced well the two major transitions at 220 and 199 nm while underestimated the positive shoulder at 286 nm (Fig. 12). Reproduction of the major transitions allowed verification of the absolute configuration of 3 as (1R,2S,3R,4S,5S,9S,10R,13R,17R,30R).
Similarly to compound 1, 4 had a chirality center at C-6 and the only difference from 1 was the position of the Δ8,14 double bond. The absolute configuration of C-6 has not been determined independently in the original paper16 thus we used TDDFT-ECD and 13C NMR DFT calculations of the C-6 epimers to assign this chirality center.
Although chiroptical methods are usually applied for the elucidation of absolute configuration and conformation by exploiting the mirror image ECD profile of the enantiomers, there are reports where they were utilized successfully to distinguish diastereomers or more than two stereoisomers.34–38 MMFF conformational search of the epimeric (1R,2S,3R,4S,5S,6R,9S,10R,13R,17R)-4 and (1R,2S,3R,4S,5S,6S,9S,10R,13R,17R)-4 produced 39 and 120 conformers, respectively, in a 21 kJ mol−1 energy window indicating considerable flexibility difference between the two diastereomers. After DFT re-optimization, ECD spectra were computed for 21/23 and 28/28 low-energy conformers in the gas-phase and the solvent-model calculations, respectively, and both epimers produced practically the same agreement with the experimental ECD (Fig. 13, S5 and S7†). This results allow verification of the absolute configuration of the core part as (1R,2S,3R,4S,5S,9S,10R,13R,17R) but it did not allow determining the AC of C-6. The C-3 substituent had s-trans conformation in most of the low-energy conformers (Fig. S4 and S6†). The configuration of the C-6 chirality center had apparently only minor contribution to the overall experimental ECD spectrum, and thus DFT 13C NMR calculation was utilized to distinguish the epimers.
The DFT-NMR methods were validated first on compounds 1–3 and 6. The widely applied mPW1PW91/6-311+G(2d,p)//B3LYP/6-31+G(d,p) level27,39 was tested on 1–3 and 6 giving good agreement for most carbons (see ESI Tables S1–S4†). Then 13C chemical shifts were computed for the epimeric (1R,2S,3R,4S,5S,6R,9S,10R,13R,17R)-4 and (1R,2S,3R,4S,5S,6S,9S,10R,13R,17R)-4. The corrected mean absolute error (CMAE)40 values were found rather similar (1.80 vs. 1.89) and the (6R) epimer appeared to be only slightly better in average (see Table S5†). A strong indication could be observed, however, for the C-5 carbon adjacent to the C-6 chirality center with more than 3 ppm difference suggesting (6R) configuration. Application of the DP4+ statistical method41,42 gave a 95.23% confidence for the (6R) epimer which is above the empirical limit suggested by the developers42 of the method for a solid determination of relative configuration.
In order to confirm further this result, the mPW1PW91/6-311+G(2d,p) SMD/CHCl3//mPW1PW91/6-311+G(2d,p) SMD/CHCl3 level40 was also tested on 1. This method including solvent model both at the DFT optimization and the NMR calculation step was not recommended by the CHESHIRE database,43 since it was developed from a very limited test molecule set for a special case but we could apply it successfully lately for a different type of natural product together with the VCD and coupling constant results.44 The solvent model NMR calculation of 1 gave similar results to those of the well-tested mPW1PW91/6-311+G(2d,p)//B3LYP/6-31+G(d,p) combination (see Table S6†). Thus it was also applied to the epimeric (1R,2S,3R,4S,5S,6R,9S,10R,13R,17R)-4 and (1R,2S,3R,4S,5S,6S,9S,10R,13R,17R)-4. CMAE values obtained at this mPW1PW91/6-311+G(2d,p) SMD/CHCl3//mPW1PW91/6-311+G(2d,p) SMD/CHCl3 combination for the two epimers were found 1.69 for the (6R) and 1.96 for the (6S) diastereomer (see Table S7†). The DP4+ analysis showed a 99.81% confidence for the (6R) epimer. The two 13C chemical shift DFT-NMR calculations produced the same conclusion for the AC of C-6 and together with the TDDFT-ECD calculations allowed elucidating the absolute configuration of 4 as (1R,2S,3R,4S,5S,6R,9S,10R,13R,17R).
The structure of compound 5 contained a Δ8,9 double bond instead of the Δ8,14 double bond of 3, while their rings A–C had identical substitution pattern. Despite their closely related structures, compound 5 had oppositely signed CEs to the corresponding ones of 3 and 4, which resulted in near mirror image experimental ECD spectra. The geometrical background for the mirror image ECD spectra of these homochiral derivatives was explored by TDDFT-ECD calculations.
Re-optimization of the 120 MMFF conformers of (1R,2S,3R,4S,5S,10S,13R,14R,17R,30R)-5 yielded 19 and 25 low-energy conformers over 1% Boltzmann-population in the gas-phase and the solvent model calculations, respectively. In all the low-energy CAM-B3LYP/TZVP PCM/MeCN conformers, the δ-lactone ring adopted a distorted boat conformation with C-15 and C-17 being out of plane (Fig. 14). This conformation of the δ-lactone ring was different from the preferred half-chair conformation of 1–4, and it was the reason of the mirror image ECD curves. The ECD spectra computed for the low-energy CAM-B3LYP/TZVP PCM/MeCN conformers reproduced well the experimental ECD curve of 5 (Fig. 15 and S10†), which confirmed that the geometries of the conformers were predicted properly. Interestingly, the B3LYP/6-31+G(d,p) gas-phase re-optimization provided different low-energy conformers, in which the δ-lactone ring had the half-chair conformation (Fig. 14) and the computed ECD gave mirror image agreement with the two high-energy ECD transitions (Fig. 15 and S11†). This result confirmed that the shift of the position of the double bond in 3 and 5 (from Δ8,14 to Δ8,9) changed the preferred conformation of the δ-lactone ring from half-chair to boat, which is manifested in near mirror ECD spectra, while the corresponding chirality centers had identical absolute configuration.
![]() | ||
Fig. 14 Comparison of the lowest-energy CAM-B3LYP/TZVP PCM/MeCN (12.7%) and B3LYP/6-31+G(d,p) gas-phase (10.8%) conformers of (1R,2S,3R,4S,5S,10S,13R,14R,17R,30R)-5. |
This is an instructive example how chiroptical properties of structurally related compounds can differ and why it is important to apply more than one functional in both the DFT optimization and the ECD calculation steps.1,2,45,46 The solvent model results reproduced the experimental ECD spectrum allowing verification of the homochiral nature of 5 with the AC of (1R,2S,3R,4S,5S,10S,13R,14R,17R,30R) despite its nearly mirror image experimental ECD.
Compared to the structure of 1, compound 6 does not contain the C-6 chirality center and since this center does not have significant contribution to the ECD spectrum, the experimental ECD spectra of 1 and 6 were very similar. The initial 34 MMFF conformers of (1R,2S,3R,4S,5S,9S,10R,13R,14S,17R)-6 were re-optimized at the B3LYP/6-31+G(d,p) and the CAM-B3LYP/TZVP PCM/MeCN levels resulting in 4 and 12 low-energy conformers, respectively. The δ-lactone ring adopted half-chair conformation in all the low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of 6, while the α,β-unsaturated ester moiety had s-trans conformation in conformers A and C and s-cis in conformers B and D (Fig. 16). The computed ECD spectra of conformers A and C had a weak positive CE at 220 nm, which did not match the experimental curve (Fig. 17). However, the intense negative computed ECD bands of conformers B and D below 225 nm overrode the positive CE of conformers A and C.14 The negative CEs of the experimental ECD spectra below 220 nm are mainly determined by the contribution of the higher-energy conformers B and D containing an s-cis enone chromophore.
![]() | ||
Fig. 16 Overlapped geometries and selected torsional angle of the four lowest-energy conformers of (1R,2S,3R,4S,5S,9S,10R,13R,14S,17R)-6. Level of optimization: CAM-B3LYP/TZVP PCM/MeCN. |
![]() | ||
Fig. 17 Experimental ECD spectrum of 6 in MeCN compared with the BH&HLYP/TZVP PCM/MeCN ECD spectra of the individual four lowest-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (1R,2S,3R,4S,5S,9S,10R,13R,14S,17R)-6. Experimental ECD spectrum of 6 was reproduced with permission from ref. 16. |
All combinations of the applied levels reproduced nicely the experimental ECD spectrum (Fig. 18) allowing verification of the homochiral (1R,2S,3R,4S,5S,9S,10R,13R,14S,17R) AC proposed on the similarity of the experimental ECD spectra of 1 and 6.16
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ra03725g |
This journal is © The Royal Society of Chemistry 2020 |