Abstract
This work describes the first determination of the absolute configuration (AC) of (−)-cubebin by means of electronic circular dichroism (ECD), supported by quantum chemical calculations. The comparison of experimental ECD with the corresponding quantum chemical prediction for the proper diastereoisomer resulted in the definitive assignment of the AC of the naturally occurring (−)-cubebin as (8R,8aR,9S). The challenging determination of the relative configuration (RC) of cubebin based only on experimental nuclear magnetic resonance (NMR) methods is stressed. Computation of the 13C and 1H NMR chemical shifts for all the possible diastereoisomers leads to values of mean absolute error and root mean square deviation that do not allow distinguishing among them. Thus, errors in stereochemical determination can easily occur even when using two-dimensional methods, which clearly demonstrate the complexity of this special case. To determine the RC of this bioactive natural compound with high level of confidence, it was necessary to combine the DP4+ method with X-ray crystallography. Therefore, employing the commonly used empirical methods to determine the AC of (−)-cubebin can easily lead to misassignment of its stereochemistry.
Keywords:
cubebin; absolute configuration; NMR simulation; X-ray crystallography; ECD
Introduction
(−)-Cubebin is a bioactive natural dibenzylbutyrolactolic lignan found in several species of plants, especially in the genus Piper.11 Pan, J. Y.; Chen, S. L.; Yang, M. H.; Wu, J.; Sinkkonen, J.; Zou, K.; Nat. Prod. Rep. 2009, 26, 1251.,22 Pissurno, A. P. R.; Laurentiz, R. S.; Rev. Virtual Quim. 2017, 9, 2. Although this compound has been known for 150 years by now, its numerous biological properties, such as analgesic, anti-inflammatory, antimicrobial, antihistaminic, vasodilator and cytotoxicity on tumor cells, have attracted the attention of several research groups all over the world.22 Pissurno, A. P. R.; Laurentiz, R. S.; Rev. Virtual Quim. 2017, 9, 2. Recent research22 Pissurno, A. P. R.; Laurentiz, R. S.; Rev. Virtual Quim. 2017, 9, 2. with cubebin has focused on the discovery of new biological properties and in obtaining semisynthetic derivatives to potentiate their pharmacological properties.
Due to its promising therapeutic effects on schistosomiasis and Chagas disease, cubebin became extremely important as a potential new drug, justifying the great chemical and pharmacological interest.22 Pissurno, A. P. R.; Laurentiz, R. S.; Rev. Virtual Quim. 2017, 9, 2. We recently reported33 Macedo, A. L.; da Silva, D. P. D.; Moreira, D. L.; de Queiroz, L. N.; Vasconcelos, T. R. A.; Araujo, G. F.; Kaplan, M. A. C.; Pereira, S. S. C.; de Almeida, E. C. P.; Valverde, A. L.; Robbs, B. K.; Biomed. Pharmacother. 2019, 110, 342. the promising cytotoxic activity of a dichloromethane fraction of the methanolic extract from leaves of Piper cernuum Vell. In that previous work, Macedo et al.33 Macedo, A. L.; da Silva, D. P. D.; Moreira, D. L.; de Queiroz, L. N.; Vasconcelos, T. R. A.; Araujo, G. F.; Kaplan, M. A. C.; Pereira, S. S. C.; de Almeida, E. C. P.; Valverde, A. L.; Robbs, B. K.; Biomed. Pharmacother. 2019, 110, 342. found one of the stereoisomers of cubebin (Figure 1) as the major component in the bioactive fraction.
Despite of its relatively simple structure, cubebin presents 3 stereogenic centers. Several authors44 Umezawa, T.; Okunishi, T.; Shimada, M.; Wood Res., Bull. Wood Res. Inst., Kyoto Univ. 1997, 84, 62.
5 Umezawa, T.; Phytochemistry Rev. 2003, 2, 371.
6 da Silva, R.; de Souza, G. H.; da Silva, A. A.; de Souza, V. A.; Pereira, A. C.; Royo, V.; Silva, M. L.; Donate, P. M.; de Matos, A. L.; Carvalho, J. C.; Bastos, J. K.; Bioorg. Med. Chem. Lett. 2005, 15, 1033.
7 Zhang, J.; Chen, J.; Liang, Z.; Zhao, C.; Chem. Biodivers. 2014, 11, 1.-88 Marcotullio, M. C.; Pelosi, A.; Curini, M.; Molecules 2014, 19, 14862. have shown that its configuration are intrinsically related to the wide range of biological activities. In 2006, de Pascoliet al.99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735. proposed an absolute configuration (AC) for cubebin based on chemical transformations. The hydroxyl group in cubebin was oxidized, transforming it into hinokinin, with the AC well established.99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735. However, even using the chemical transformations, the configuration of the anomeric center in cubebin remains unclear, as seen in other butyrolactolic lignans.99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735.
10 Li, N.; Wu, J. L.; Sakai, J.; Ando, M.; J. Nat. Prod. 2003, 66, 1421.
11 Barrero, A. F.; Haïdour, A.; Muñoz Dorado, M.; Grávalos, D. G.; Quesada, T. G.; J. Nat. Prod. 1994, 57, 713.-1212 Gozler, B.; Rentsch, D.; Gozler, T.; Unver, N.; Hesse, M.; Phytochemistry 1996, 42, 695. Additionally, the use of correlations or spectral comparison with compounds of known AC, usually does not yield reliable results regarding the absolute configuration.1313 Polavarapu, P. L.; Chem. Rec. 2007, 7, 125. In this particular case, empirical methodologies such as chemical transformations can easily lead to misassignment, since the electronic circular dichroism (ECD) spectrum of different diastereoisomers of cubebin could show similar bands.99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735. Thus, revisiting the configuration of cubebin using a more reliable stereochemical method is highly recommended.
Nowadays, the combination of nuclear magnetic resonance (NMR), one of the most powerful techniques for relative configuration determination, with ECD, which explores the electronic properties of the molecules, both associated with high level quantum chemical calculations, is considered the best approach for reliable stereochemical assignments of natural products.1313 Polavarapu, P. L.; Chem. Rec. 2007, 7, 125.
14 Petrovic, A. G.; Navarro-Vázquez, A.; Alonso-Gómez, J. L.; Curr. Org. Chem. 2010, 14, 1612.
15 Batista, A. N. L.; dos Santos, F. M. J.; Batista, J. M. J.; Cass, Q. B.; Molecules 2018, 23, 492.
16 Batista, J. M. J.; Blanch, E. W.; Bolsani, S.; Nat. Prod. Rep. 2015, 32, 1280.-1717 Pescitelli, G.; Bruhn, T.; Chirality 2016, 28, 466. Thus, aiming the full identification of the compound present in the P. cernuum cytotoxic fraction,33 Macedo, A. L.; da Silva, D. P. D.; Moreira, D. L.; de Queiroz, L. N.; Vasconcelos, T. R. A.; Araujo, G. F.; Kaplan, M. A. C.; Pereira, S. S. C.; de Almeida, E. C. P.; Valverde, A. L.; Robbs, B. K.; Biomed. Pharmacother. 2019, 110, 342. we report the isolation and the first analysis of the stereostructure of naturally occurring cubebin by means of NMR and ECD, both supported by their corresponding quantum chemical prediction.
Results and Discussion
The isolated compound was obtained as a white amorphous solid, whose molecular formula was suggested as C20H20O6 based on [M]+ m/z 356 by gas chromatography-mass spectrometry (GC-MS). The 1H and 13C NMR were similar to those reported for cubebin (Table 1).99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735. The 13C NMR spectrum indicated 20 carbons and the 1H NMR showed signals that were attributed to the aromatic hydrogens (δ 6.66, 6.55, 6.49, 6.65, 6.51 and 6.44), two methylenedioxy groups (δ 5.83), two methynic hydrogens (δ 1.94 and 1.83), two methylenic hydrogens (δ 2.46, 2.30 and 2.39), and one anomeric hydrogen (δ 4.90). The 2D spectral data of HSQC (heteronuclear single quantum coherence) and HMBC (heteronuclear multiple bond correlation spectroscopy) allowed the full assignment of the cubebin basic structure.
Mass spectrometry and NMR spectroscopy data suggested that the compound has a structure compatible with the lignan cubebin. Thus, we focused to determine the relative configuration for the three stereogenic centers present in the structure. The eight possible stereoisomers comprise four pairs of enantiomers. As each pair of enantiomers displays exactly the same NMR chemical shifts in isotropic media the first approach was to compute the 13C and 1H NMR chemical shifts for the four possible diastereoisomers of cubebin (see Figure 1) and compare them with those obtained experimentally for the isolated compound (Tables 2 and 3, respectively).
Comparison of computed 13C NMR chemical shifts (δscal) of the four possible diastereoisomers of cubebin with experimental data (δexp, 125 MHz, DMSO-d6) for the corresponding isolated natural product
Comparison of theoretical 1H NMR chemical shifts (δscal) of the four diastereoisomers of cubebin with experimental data (δexp, 500 MHz, DMSO-d6) for the isolated natural product
Table 2 shows the comparison of the computed 13C NMR chemical shifts shifts for the four diastereoisomers of cubebin with experimental data. Mean absolute errors (MAE) and root mean square deviations (RMSD) are also given. Although diastereoisomers 2 and 4 show the lowest 13C NMR MAE and RMSD values, the differences for diastereoisomers 1 and 3 are small, not allowing for definitive determination of the relative configuration of the isolated compound, based only on these 13C NMR statistical parameters.
Table 3 shows the comparison of the simulated 1H chemical shifts obtained for the four diastereoisomers of cubebin with experimental data. Once again, 1H NMR MAE and RMSD values are similar and cannot discriminate between the possible diastereoisomers, although 1 and 2 show the lowest values. These results clearly indicate that using only data based on the 13C and 1H chemical shifts it is difficult to determine the relative configuration of cubebin. This could easily lead to misassignment if empirical methodologies are exclusively used to determine the stereochemistry of this potential biological lignan.
After the 13C and 1H NMR chemical shifts calculations, even with diastereoisomers 2 showing the lowest 13C and 1H NMR MAE and RMSD values, it is not possible to obtain a reliable discrimination between the possible diastereoisomers 1-4 (Figure 1). Trying to increase the reliability of the discrimination, we employed the DP4+ procedure, proposed by Sarotti and co-workers.1818 Grimblat, N.; Zanardi, M. M.; Sarotti, A. M.; J. Org. Chem. 2015, 80, 12526; sarotti-NMR.weebly.com, accessed in May 2020.
sarotti-NMR.weebly.com...
The basic idea of this procedure is to compare the experimental NMR data with all possible simulated structures. Taking into account the differences in chemical shifts obtained between the simulated and experimental values, each possible structure is ranked after the application of a Student’s t-test, based on statistical treatment.1818 Grimblat, N.; Zanardi, M. M.; Sarotti, A. M.; J. Org. Chem. 2015, 80, 12526; sarotti-NMR.weebly.com, accessed in May 2020.
sarotti-NMR.weebly.com...
This method has been established to provide a more reliable assignment procedure, when comparing with the statistical parameters,1818 Grimblat, N.; Zanardi, M. M.; Sarotti, A. M.; J. Org. Chem. 2015, 80, 12526; sarotti-NMR.weebly.com, accessed in May 2020.
sarotti-NMR.weebly.com...
19 Zanardi, M. M.; Suárez, A. G.; Sarotti, A. M.; J. Org. Chem. 2017, 82, 1873.-2020 Sarotti, A. M.; Org. Biomol. Chem. 2018, 16, 944. such as the MAE and RMSD values listed in Tables 2 and 3. The DP4+ value obtained using both the 13C and 1H NMR data in the scaled form (all data, sDP4+) enabled the assignment of diastereoisomer 2 as the most likely one, in overall confidence level of 97.9%, against 1.15, 0.87 and 0.08% for diastereoisomers 1, 3 and 4, respectively (see Figure 2).
sDP4+ values (all data in the scaled form) obtained by correlating the calculated 1H and 13C NMR chemical shifts of the four diastereoisomers with the experimental NMR data of the isolated natural product.
For further verification of the reliability of this procedure, we also carried out X-ray diffraction study of the isolated cubebin. White single crystals suitable for X-ray diffraction analysis were obtained by slow evaporation of a cubebin acetone solution. Based on the X-ray crystallographic analysis, the molecular structure was confirmed as the diastereoisomer 2 (Figure 1). Figure 3 shows the Oak Ridge thermal ellipsoid plot (ORTEP) diagram of the compound. Crystal data and refinements are provided in Table S1 of the Supplementary Information Supplementary Information Crystallographic data (excluding structure factors) for the structures in this work were deposited in the Cambridge Crystallographic Data Centre as supplementary publication number CCDC 1991277. Copies of the data can be obtained, free of charge, via https://www.ccdc.cam.ac.uk/structures/. Crystallographic data is also available for download in CIF format at http://jbcs.sbq.org.br. Supplementary information is available free of charge at http://jbcs.sbq.org.br as PDF file. (SI) section.
Asymmetric unit representation of dibenzylbutyrolactolic lignan. Displacement ellipsoids are drawn at the 50% probability level.
Therefore, the relative stereochemical configuration of the isolated compound was unambiguously established as (8R*,8aR*,9S*) by the robust combination of 1D and 2D NMR studies, comparison of 1H and 13C NMR chemical shifts simulation, application of the sDP4+ method and X-ray crystallography study. Next, using the proper relative configuration, we focused on the determination of the AC of the isolated compound by means of chiroptical spectroscopy.
As mentioned in the context of ECD spectra, de Pascoli et al.99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735. reported that the ECD spectrum of different diastereoisomers of cubebin could display similar bands, directly related to the stereochemistry of the stereogenic centers C-8 and C-8a.99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735. Thus, any uncertainty about the relative configuration of the anomeric carbons can easily lead to incorrect stereochemistry assignment if only the empirical methodology is used, e.g., by just comparing the ECD spectrum obtained experimentally with those reported in the literature.Cases like these have already been seen in the literature2121 Junior, F. M.; Velozo, L. S.; de Carvalho, E. M.; Marques, A. M.; Borges, R. M.; Trindade, A. P.; dos Santos, M. I.; de Albuquerque, A. C.; Costa, F. L.; Kaplan, M. A.; de Amorim, M. B.; Molecules 2013, 18, 13520.,2222 Junior, F. M. S.; Covington, C. L.; de Amorim, M. B.; Velozo, L. S. M.; Kaplan, M. A. C.; Polavarapu, P. L.; J. Nat. Prod. 2014, 77, 1881. and demand a more robust method to determine their relative configuration, preferentially by independent tools, including X-ray crystallography and NMR spectroscopy.1717 Pescitelli, G.; Bruhn, T.; Chirality 2016, 28, 466. Nevertheless, a combination of experimental and theoretical methods is highly recommended, to avoid errors in the stereochemical determination.1717 Pescitelli, G.; Bruhn, T.; Chirality 2016, 28, 466.,2121 Junior, F. M.; Velozo, L. S.; de Carvalho, E. M.; Marques, A. M.; Borges, R. M.; Trindade, A. P.; dos Santos, M. I.; de Albuquerque, A. C.; Costa, F. L.; Kaplan, M. A.; de Amorim, M. B.; Molecules 2013, 18, 13520.
22 Junior, F. M. S.; Covington, C. L.; de Amorim, M. B.; Velozo, L. S. M.; Kaplan, M. A. C.; Polavarapu, P. L.; J. Nat. Prod. 2014, 77, 1881.
23 Junior, F. M. S.; Covington, C. L.; de Albuquerque, A. C. F.; Lobo, J. F. R.; Borges, R. M.; de Amorim, M. B.; Polavarapu, P. L.; J. Nat. Prod. 2015, 78, 2617.
24 dos Santos Jr., F. M.; Bicalho, K. U.; Calisto, I. H.; Scatena, G. S.; Fernandes, J. B.; Cass, Q. B.; Batista Jr., J. M.; Org. Biomol. Chem. 2018, 16, 4509.-2525 Batista, N. L. A.; dos Santos Jr., F. M.; Valverde, A. L.; Batista Jr., J. M.; Org. Biomol. Chem. 2019, 17, 9772. Therefore, the use of ECD in this case is recommended to unambiguous and direct determination of the absolute stereochemistry, in combination with quantum chemical calculations and a high level of confidence in their relative configuration.1717 Pescitelli, G.; Bruhn, T.; Chirality 2016, 28, 466. It is noteworthy that this combined methodology has been deemed as one of the most powerful methodology for the stereochemical analysis of natural products and synthetic molecules.2121 Junior, F. M.; Velozo, L. S.; de Carvalho, E. M.; Marques, A. M.; Borges, R. M.; Trindade, A. P.; dos Santos, M. I.; de Albuquerque, A. C.; Costa, F. L.; Kaplan, M. A.; de Amorim, M. B.; Molecules 2013, 18, 13520.
22 Junior, F. M. S.; Covington, C. L.; de Amorim, M. B.; Velozo, L. S. M.; Kaplan, M. A. C.; Polavarapu, P. L.; J. Nat. Prod. 2014, 77, 1881.
23 Junior, F. M. S.; Covington, C. L.; de Albuquerque, A. C. F.; Lobo, J. F. R.; Borges, R. M.; de Amorim, M. B.; Polavarapu, P. L.; J. Nat. Prod. 2015, 78, 2617.
24 dos Santos Jr., F. M.; Bicalho, K. U.; Calisto, I. H.; Scatena, G. S.; Fernandes, J. B.; Cass, Q. B.; Batista Jr., J. M.; Org. Biomol. Chem. 2018, 16, 4509.
25 Batista, N. L. A.; dos Santos Jr., F. M.; Valverde, A. L.; Batista Jr., J. M.; Org. Biomol. Chem. 2019, 17, 9772.
26 Gimenes, L.; Batista Jr., J. M.; dos Santos Jr., F. M.; da Silva, S. M.; Luna-Dulcey, L.; Ellena, J.; Cominetti, M. R.; Silva, F. G. F.; Viera, P. C.; Fernandes, J. B.; Staerk, D.; Phytochemistry 2019, 163, 38.
27 Amorim, M. R.; Hilário, F.; Junior, F. M. S.; Batista, J. M.; Bauab, T. M.; Araújo, A. R.; Carlos, I. Z.; Vilegas, W.; dos Santos, L. C.; Planta Med. 2019, 85, 957.
28 Covington, C. L.; Junior, F. M. S.; Silva, J. H. S.; Kuster, R. M.; de Amorim, M. B.; Polavarapu, P. L.; J. Nat. Prod. 2016, 79, 530.-2929 Jimenez, D. E. Q.; Barreiro, J. C.; dos Santos Jr., F. M.; de Vasconcellos, S. P.; Porto, A. L. M.; Batista Jr., J. M.; Chirality 2019, 31, 534. Thus, to unequivocal determination of the absolute stereochemistry of (−)-cubebin, we employed a combination of the experimental ECD spectrum with the corresponding one computed by quantum chemical methods.
The experimental ECD spectrum in the 190-300 nm (Figure 4) region measured in acetonitrile (ACN) clearly shows a negative Cotton effect at approximately 200 nm, associated with the benzodioxole chromophore in cubebin.99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735.
Figure 5 shows the comparison of the experimental ECD spectrum of isolated cubebin in ACN with the simulated data for the Boltzmann average of the 14 lowest-energy conformers identified for diastereoisomer 2 (8R,8aR,9S) at the CAM-B3LYP/PCM(ACN)/TZVP (Coulomb-attenuating method Becke, 3-parameter, Lee-Yang-Parr/polarizable continuum model(acetonitrile)/valence triple-zeta polarization) level. Range-separated hybrid functional such as CAM-B3LYP, combined with the TZVP basis sets, were selected due to their good performance for ECD calculations.2929 Jimenez, D. E. Q.; Barreiro, J. C.; dos Santos Jr., F. M.; de Vasconcellos, S. P.; Porto, A. L. M.; Batista Jr., J. M.; Chirality 2019, 31, 534. The negative Cotton effect around 200 nm is perfectly reproduced in the calculations. The good agreement observed between the experimental and the time-dependent density functional theory (TDDFT) simulated spectrum allowed the unequivocal assignment of the absolute configuration of cubebin as (8R,8aR,9S).
Comparison of the observed ECD spectrum of cubebin in ACN (black trace) with the calculated [CAM-B3LYP/PCM(ACN)/TZVP] ECD spectrum of the Boltzmann average of the 14 lowest-energy conformers identified for diastereoisomer 2 (8R,8aR,9S) of cubebin (red trace).
Conclusions
In conclusion, with the proper assignment of the relative configuration of (−)-cubebin, we described the first determination of its absolute configuration by the use of ECD, supported by quantum chemical calculations. The comparison of experimental ECD with the corresponding quantum chemical prediction of the correct diastereoisomer, confirmed by X-ray analysis and NMR studies, resulted in the definitive assignment of the AC of cubebin. When taken together, the analyses made here confirmed the AC of the naturally occurring (−)-cubebin as (8R,8aR,9S).
Experimental
General
Solvents used for extraction and chromatography were acquired from Vetec® (Duque de Caxias, RJ, Brazil) and used without further purification. Acetonitrile used for ECD, dichloromethane for GC-MS and dimethyl sulfoxide (DMSO-d6) were acquired from Sigma-Aldrich (São Paulo, Brazil). Silica gel used for vacuum liquid chromatography, flash chromatography and preparative thin layer chromatography were bought from Merck (Darmstadt, Germany).
Plant material
Piper cernuum Vell. was collected in January 2015 in Parque Nacional da Serra dos Órgãos, Guapimirim, RJ, Brazil (22°29’36”S, 42°59’54”W) and identified by Dr Elsie Franklin Guimarães, a specialist in Piperaceae family. Voucher specimen was deposited in the Herbarium of the Botanical Garden Research Institute of Rio de Janeiro under the number RB457067.
Extraction and isolation
Dried and pulverized leaves of P. cernuum (600 g) were extracted by static maceration in methanol and concentrated under reduced pressure on rotary evaporator. The methanolic extract (80 g) was re-suspended in mixture of MeOH:H2O (7:3) and successively extracted with n-hexane, dichloromethane, ethyl acetate and n-butanol. Aiming at the isolation of cubebin, the dichloromethane fraction (15 g) was submitted to a vacuum liquid column chromatography (silica gel 60, 2-25 µm) using a polarity gradient with small changes in the mobile phase composition with different mixtures of hexane:ethyl acetate (99:1-0:1, 11 fractions), ethyl acetate:methanol (9:1-0:1, four fractions) and methanol:water (8:2 and 6:4, two fractions), yielding 17 fractions. Fraction 9 (1.219 g) was eluted in normal phase flash column chromatography (silica gel 60, 230-400 mesh) with a gradient of polarity that changed from hexane:ethyl acetate to methanol, yielding 26 fractions. Fraction 3 (255 mg) was subjected to reverse phase flash column chromatography (octadecyl-functionalized silica gel 60, 200-400 mesh) eluted in different proportions of methanol:water (6:4 and 8:2, five fractions each), methanol:dicloromethane (2:8 and 8:2, five fractions each) and dichloromethane (one fraction), affording 21 fractions. Fractions 8 and 9 were combined (89 mg) and subjected to preparative layer plate chromatography (silica gel 60 GF254, 2-25 µm), eluted with dichloromethane:methanol (95:5) and observed under 254 nm UV light to afford cubebin (27 mg). The basic structure of cubebin was confirmed by GC-MS and comparison of its 1H and 13C NMR spectra with the published data by de Pascoli et al.99 de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735. For the X-ray analysis, the amorphous powder was subjected to the crystallization: the powder was solubilized with acetone followed by addition of hexane until turbidity. The sample was dried at room temperature, furnishing the cubebin crystals.
GC-MS analysis
GC-MS analysis of cubebin was performed using a Hewlett-Packard 6890 GC coupled to a Hewlett-Packard 5972 MS equipped with a DB-5MS (5% diphenyl/95% dimethyl polysiloxane) fused silica capillary column (30 m × 0.25 µm ID (inner diameter) × 0.25 µm df (film thickness)). For GC-MS detection, an electron ionization system was operated in electron impact mode with an ionization energy of 70 eV. Helium gas (99.999%) was used as a carrier gas at a constant flow rate of 0.5 mL min−1) The solvent delay was 0 to 3.5 min, and the total GC-MS running time was 60 min. The injector temperature was maintained at 270 °C, the ion-source temperature was 250 °C and the oven temperature was programmed from 70 °C, with an increase of 4 °C min−1 to 315 °C. Mass spectra were taken at 70 eV; a scan interval of 0.5 s and fragments from 40 to 700 Da. Data acquisition and handling were done via ChemStation software.3030 GC/MSD ChemStation Software; Agilent, Santa Clara, CA, USA, 2005.
NMR analysis
The NMR spectra were obtained on a Varian VNMRS 500 MHz, operating at 500 MHz for 1H and 125 MHz for 13C frequencies, with tetramethylsilane (TMS) as the internal standard. Dimethyl sulfoxide-d6 was chosen as solvent to ensure good solubility of the sample.
Electronic circular dichroism
ECD spectrum was obtained for solutions (1 mg mL−1) in ACN on a Chirascan (Applied Photophysics Ltd., UK) in the 190-350 nm region in a 0.1 cm path length cell and a bandwidth of 1 nm at room temperature. The reported Δε (differential molar absorptivity) values are expressed in units of L mol−1 cm−1. Acetonitrile was chosen as solvent for the preparation of the sample due to its cutoff wavelength value of 190 nm, below which the solvent absorption becomes excessive.
Specific rotation
The optical rotations at 589 nm were measured for solutions (0.5 mg mL−1) in ACN on a Jasco P-200 polarimeter (Jasco Brasil, São Paulo, Brazil), at 25 °C, using a 2.5 mm cell. Specific rotation [α](D)25: −83°.
Calculations
Randomized conformational searches were performed for each of the four possible diastereoisomers of cubebin (Figure 1), using the Monte Carlo method and MMFF force field in the Spartan’10 software package.3131 Spartan; Wavefunction, Irvine, CA, USA, 2014. The more stable conformations, within a 10.0 kcal mol−1 energy range, were reoptimized using the mPW1PW91/6-31G* DFT method. The resulting conformers within the energy range of 2.5 kcal mol−1, accounting for more than 99.0% of the total Boltzmann population, were then used for NMR and ECD calculations. This theoretical approach computes NMR chemical shifts with high accuracy and low computational cost, being considered one of high accuracy methods available.3232 Costa, F. L. P.; de Albuquerque, A. C. F.; Borges, R. M.; Junior, F. M. S.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2014, 11, 219. For each selected structure, 13C and 1H NMR chemical shifts were calculated using the GIAO-mPW1PW91/6-31G(d)//mPW1PW91/6-31G(d) DFT approach. The Gaussian 09 software package3333 Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery Jr., J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, N. J.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J.; Gaussian 09, Revision A.01; Gaussian, Inc., Wallingford, CT, USA, 2009. was used. Calculated chemical shifts were scaled using a linear equation (δscal = 1.05δcalc − 1.22).3232 Costa, F. L. P.; de Albuquerque, A. C. F.; Borges, R. M.; Junior, F. M. S.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2014, 11, 219. This approach allows the cancellation of systematic errors as scaling factors are determined via linear regression analyses between the calculated and experimental chemical shifts of a pool of adequately chosen compounds. This is a good alternative to avoid more demanding calculations necessary to include solvent and rovibration effects.3232 Costa, F. L. P.; de Albuquerque, A. C. F.; Borges, R. M.; Junior, F. M. S.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2014, 11, 219. To apply the DP4+ parameter, it was used the Excel spreadsheet provided by the authors.1818 Grimblat, N.; Zanardi, M. M.; Sarotti, A. M.; J. Org. Chem. 2015, 80, 12526; sarotti-NMR.weebly.com, accessed in May 2020.
sarotti-NMR.weebly.com...
After assignment of the relative configuration, the correct diastereoisomer were reoptimized at B3LYP/6-31G(d) and used for simulation of the ECD spectrum with CAM-B3LYP/TZVP//B3LYP/6-31G(d) TD-DFT approach, using the PCM solvent (ACN) method, in order to take into account solvation effects (see SI section). The final ECD spectra were generated according to Boltzmann weighting of the lowest-energy conformers identified and plotted using Origin 8 software.3434 OriginPro 8, v. 80724; OriginLab Corporation, Northampton, MA, USA, 2020.
X-ray data collection and structure refinement
X-ray diffraction data for cubebin was carried out with radiation Mo Ka (λ = 0.71073 Å) in Bruker D8 Venture diffractometer.3535 Bruker APEX3, SAINT and SADABS; Bruker AXS Inc., Madison, Wisconsin, USA, 2016. Using Olex2,3636 Dolomanov, O. V.; Bourhis, L. J.; Gildea, R. J.; Howard, J. A. K.; Puschmann, H.; J. Appl. Crystallogr. 2009, 42, 339. the structure was solved with the olex2.solve3737 Bourhis, L. J.; Dolomanov, O. V.; Gildea, R. J.; Howard, J. A. K.; Puschmann, H.; Acta Crystallogr. 2015, A71, 59. structure solution program using Charge Flipping and refined with the ShelXL3838 Sheldrick, G. M.; Acta Crystallogr. 2015, C71, 3. refinement package using least squares minimization. All H atoms were refined with fixed individual displacement parameters [Uiso/(H) = 1.2 Ueq (Csp2 and Car) or 1.5 Ueq (Csp3)] using a riding model. All non-hydrogen atoms were refined anisotropically. Molecular graphics were drawn using ORTEP-3 for Mercury.3939 Macrae, C. F.; Edgington, P. R.; McCabe, P.; Pidcock, E.; Shields, G. P.; Taylor, R.; Towler, M.; van de Streek, J.; J. Appl. Crystallogr. 2006, 39, 453.
Legal information
The sample collection of P. cernuum was authorized by the Chico Mendes Institute of Biology (ICMBio) through the Biodiversity authorization and information system (SISBio) in document number 45566-1. Access was authorized by the Management Board of Genetic Heritage Council in process number 010771/ 2014-0.
Acknowledgments
The authors thank LDRX-UFF (www.ldrx.uff.br/) for the X-ray diffraction analysis and for the facilities. We wish to thank CNPq, and FAPERJ (grant 211.319-2019), for financial support. This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brazil (CAPES), finance code 001. We dedicate this paper to memory of Professor Mauro Barbosa de Amorim (1954-2019), a masterful teacher and mentor, and more than that, a father, a husband and a friend.
Supplementary Information
Crystallographic data (excluding structure factors) for the structures in this work were deposited in the Cambridge Crystallographic Data Centre as supplementary publication number CCDC 1991277. Copies of the data can be obtained, free of charge, via https://www.ccdc.cam.ac.uk/structures/. Crystallographic data is also available for download in CIF format at http://jbcs.sbq.org.br.
Supplementary information is available free of charge at http://jbcs.sbq.org.br as PDF file.
References
-
1Pan, J. Y.; Chen, S. L.; Yang, M. H.; Wu, J.; Sinkkonen, J.; Zou, K.; Nat. Prod. Rep. 2009, 26, 1251.
-
2Pissurno, A. P. R.; Laurentiz, R. S.; Rev. Virtual Quim. 2017, 9, 2.
-
3Macedo, A. L.; da Silva, D. P. D.; Moreira, D. L.; de Queiroz, L. N.; Vasconcelos, T. R. A.; Araujo, G. F.; Kaplan, M. A. C.; Pereira, S. S. C.; de Almeida, E. C. P.; Valverde, A. L.; Robbs, B. K.; Biomed. Pharmacother. 2019, 110, 342.
-
4Umezawa, T.; Okunishi, T.; Shimada, M.; Wood Res., Bull. Wood Res. Inst., Kyoto Univ. 1997, 84, 62.
-
5Umezawa, T.; Phytochemistry Rev. 2003, 2, 371.
-
6da Silva, R.; de Souza, G. H.; da Silva, A. A.; de Souza, V. A.; Pereira, A. C.; Royo, V.; Silva, M. L.; Donate, P. M.; de Matos, A. L.; Carvalho, J. C.; Bastos, J. K.; Bioorg. Med. Chem. Lett. 2005, 15, 1033.
-
7Zhang, J.; Chen, J.; Liang, Z.; Zhao, C.; Chem. Biodivers. 2014, 11, 1.
-
8Marcotullio, M. C.; Pelosi, A.; Curini, M.; Molecules 2014, 19, 14862.
-
9de Pascoli, I. C.; Nascimento, I. R.; Lopes, L. M. X.; Phytochemistry 2006, 67, 735.
-
10Li, N.; Wu, J. L.; Sakai, J.; Ando, M.; J. Nat. Prod 2003, 66, 1421.
-
11Barrero, A. F.; Haïdour, A.; Muñoz Dorado, M.; Grávalos, D. G.; Quesada, T. G.; J. Nat. Prod. 1994, 57, 713.
-
12Gozler, B.; Rentsch, D.; Gozler, T.; Unver, N.; Hesse, M.; Phytochemistry 1996, 42, 695.
-
13Polavarapu, P. L.; Chem. Rec. 2007, 7, 125.
-
14Petrovic, A. G.; Navarro-Vázquez, A.; Alonso-Gómez, J. L.; Curr. Org. Chem. 2010, 14, 1612.
-
15Batista, A. N. L.; dos Santos, F. M. J.; Batista, J. M. J.; Cass, Q. B.; Molecules 2018, 23, 492.
-
16Batista, J. M. J.; Blanch, E. W.; Bolsani, S.; Nat. Prod. Rep. 2015, 32, 1280.
-
17Pescitelli, G.; Bruhn, T.; Chirality 2016, 28, 466.
-
18Grimblat, N.; Zanardi, M. M.; Sarotti, A. M.; J. Org. Chem. 2015, 80, 12526; sarotti-NMR.weebly.com, accessed in May 2020.
» sarotti-NMR.weebly.com -
19Zanardi, M. M.; Suárez, A. G.; Sarotti, A. M.; J. Org. Chem. 2017, 82, 1873.
-
20Sarotti, A. M.; Org. Biomol. Chem. 2018, 16, 944.
-
21Junior, F. M.; Velozo, L. S.; de Carvalho, E. M.; Marques, A. M.; Borges, R. M.; Trindade, A. P.; dos Santos, M. I.; de Albuquerque, A. C.; Costa, F. L.; Kaplan, M. A.; de Amorim, M. B.; Molecules 2013, 18, 13520.
-
22Junior, F. M. S.; Covington, C. L.; de Amorim, M. B.; Velozo, L. S. M.; Kaplan, M. A. C.; Polavarapu, P. L.; J. Nat. Prod. 2014, 77, 1881.
-
23Junior, F. M. S.; Covington, C. L.; de Albuquerque, A. C. F.; Lobo, J. F. R.; Borges, R. M.; de Amorim, M. B.; Polavarapu, P. L.; J. Nat. Prod. 2015, 78, 2617.
-
24dos Santos Jr., F. M.; Bicalho, K. U.; Calisto, I. H.; Scatena, G. S.; Fernandes, J. B.; Cass, Q. B.; Batista Jr., J. M.; Org. Biomol. Chem. 2018, 16, 4509.
-
25Batista, N. L. A.; dos Santos Jr., F. M.; Valverde, A. L.; Batista Jr., J. M.; Org. Biomol. Chem. 2019, 17, 9772.
-
26Gimenes, L.; Batista Jr., J. M.; dos Santos Jr., F. M.; da Silva, S. M.; Luna-Dulcey, L.; Ellena, J.; Cominetti, M. R.; Silva, F. G. F.; Viera, P. C.; Fernandes, J. B.; Staerk, D.; Phytochemistry 2019, 163, 38.
-
27Amorim, M. R.; Hilário, F.; Junior, F. M. S.; Batista, J. M.; Bauab, T. M.; Araújo, A. R.; Carlos, I. Z.; Vilegas, W.; dos Santos, L. C.; Planta Med. 2019, 85, 957.
-
28Covington, C. L.; Junior, F. M. S.; Silva, J. H. S.; Kuster, R. M.; de Amorim, M. B.; Polavarapu, P. L.; J. Nat. Prod. 2016, 79, 530.
-
29Jimenez, D. E. Q.; Barreiro, J. C.; dos Santos Jr., F. M.; de Vasconcellos, S. P.; Porto, A. L. M.; Batista Jr., J. M.; Chirality 2019, 31, 534.
-
30GC/MSD ChemStation Software; Agilent, Santa Clara, CA, USA, 2005.
-
31Spartan; Wavefunction, Irvine, CA, USA, 2014.
-
32Costa, F. L. P.; de Albuquerque, A. C. F.; Borges, R. M.; Junior, F. M. S.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2014, 11, 219.
-
33Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery Jr., J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, N. J.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J.; Gaussian 09, Revision A.01; Gaussian, Inc., Wallingford, CT, USA, 2009.
-
34OriginPro 8, v. 80724; OriginLab Corporation, Northampton, MA, USA, 2020.
-
35Bruker APEX3, SAINT and SADABS; Bruker AXS Inc., Madison, Wisconsin, USA, 2016.
-
36Dolomanov, O. V.; Bourhis, L. J.; Gildea, R. J.; Howard, J. A. K.; Puschmann, H.; J. Appl. Crystallogr. 2009, 42, 339.
-
37Bourhis, L. J.; Dolomanov, O. V.; Gildea, R. J.; Howard, J. A. K.; Puschmann, H.; Acta Crystallogr. 2015, A71, 59.
-
38Sheldrick, G. M.; Acta Crystallogr. 2015, C71, 3.
-
39Macrae, C. F.; Edgington, P. R.; McCabe, P.; Pidcock, E.; Shields, G. P.; Taylor, R.; Towler, M.; van de Streek, J.; J. Appl. Crystallogr. 2006, 39, 453.
Publication Dates
-
Publication in this collection
09 Oct 2020 -
Date of issue
Oct 2020
History
-
Received
18 Mar 2020 -
Accepted
29 May 2020