Acessibilidade / Reportar erro

Conformational Analysis, Experimental and GIAO-DFT 13C NMR Chemical Shift Calculation on 2’-Hydroxy-3,4,5-trimethoxy-chalcone

Abstract

In this paper we investigated the ability of the GIAO-mPW1PW91/6-31G(d)//mPW1PW91/6-31G(d) level of theory to predict the 13C nuclear magnetic resonance (NMR) chemical shifts of the 2’-hydroxy-3,4,5-trimethoxy-chalcone molecule. Two different approaches were used. First: the absolute shieldings σ for all carbon atoms in each geometrically optimized conformers of the 2’-hydroxy-3,4,5-trimethoxy-chalcone molecule were calculated at the GIAO-mPW1PW91/6-31G(d)//mPW1PW91/6-31G(d) level of theory. This approach is further used to generate weighted average values for each atom considering the previously obtained conformational distribution. Second: only the σ for the lowest energetic conformer will be taken to account. The robustness of the method was evaluated for two other chalcones: (E)-1-(4-hydroxy-3-methoxyphenyl)-3-(3,5-di-tert-butyl-4-hydroxyphenyl)prop-2-en-1-one and (E)-1-(4-aminophenyl)-3-(3,4-dimethoxyphenyl)prop-2-en-1-one, corroborating the ability of the method in chemical shift prevision. Although, both approaches were able to reproduce the chemical shifts of the 2’-hydroxy-3,4,5-trimethoxy-chalcone, significant differences in the calculated values for C-4 and methoxy carbons were observed. The best results were obtained using the second approach (II).

Keywords:
chalcone; 2’-hydroxy-3,4,5-trimethoxy-chalcone; GIAO-NMR; scaling factor; NMR spectroscopy


Introduction

Chalcones (1,3-diaryl-2-propen-1-ones) are important intermediates for the synthesis of compounds such as flavonoids, isoflavonoids and their derivatives.11 Daskiewicz, J.-B.; Depeint, F.; Viornery, L.; Bayet, C.; Comte-Sarrazin, G.; Comte, G.; Gee, J. M.; Johnson, I. T.; Ndjoko, K.; Hostettmann, K.; Barron, D.; J. Med. Chem. 2005, 48, 2790. Besides having a physiological role in plants, flavonoids have also been reported to have a wide variety of biological activities, including antioxidant,22 Stevens, J. F.; Miranda, C. L.; Frei, B.; Buhler, D. R.; Chem. Res. Toxicol. 2003, 16, 1277.

3 Nishida, J.; Kawabata J.; J. Biosci. Biotechnol. Biochem. 2006, 70, 193.

4 Gacche, R. N.; Dhole, N. A.; Kamble, S. G.; Bandgar, B. P. J.; Enzyme Inhib. Med. Chem. 2007, 23, 28.

5 Vogel, S.; Ohmayer, S.; Brunner, G.; Heilmann, J.; Bioorg. Med. Chem. 2008, 16, 4286.
-66 Jung, J.-C.; Jang, S.; Lee, Y.; Min, Y.; Lim, E.; Jung, H.; Oh, M.; Oh, S.; Jung, M.; J. Med. Chem. 2008, 51, 4054. antibacterial,77 Sugamoto, K.; Kurogi, C.; Matsushita, Y.; Matsui, T.; Tetrahedron Lett. 2008, 49, 6639.,88 Avila, H. P.; Smania, E.; Monache, F. D.; Smania, A.; Bioorg. Med. Chem. 2008, 16, 9790. anticancer,99 Lawrence, N. J.; Patterson, R. P.; Ooi, L. L.; Cook, D.; Ducki, S.; Bioorg. Med. Chem. Lett. 2006, 16, 5844.,1010 Cabrera, M.; Simoens, M.; Falchi, G.; Lavaggi, M. L.; Piro, O. E.; Castellano, E. E.; Vidal, A.; Azqueta, A.; Monge, A.; de Cerain, A. L.; Sagrera, G.; Seoane, G.; Cerecettoand, H.; Gonzalez, M.; Bioorg. Med. Chem. 2007, 15, 3356. antiangiogenic ones.1111 Boumendjel, A.; Boccard, J.; Carrupt, P.-A.; Nicolle, E.; Blanc, M.; Geze, A.; Choisnard, L.; Wouessidjewe, D.; Matera, E.-L.; Dumontet, C.; J. Med. Chem. 2008, 51, 2307. The growing interest in these compounds and their potential use in medicinal applications is indicated by the increasing number of publications concerning the synthesis and biological evaluation of chalcone analogues.11 Daskiewicz, J.-B.; Depeint, F.; Viornery, L.; Bayet, C.; Comte-Sarrazin, G.; Comte, G.; Gee, J. M.; Johnson, I. T.; Ndjoko, K.; Hostettmann, K.; Barron, D.; J. Med. Chem. 2005, 48, 2790.

2 Stevens, J. F.; Miranda, C. L.; Frei, B.; Buhler, D. R.; Chem. Res. Toxicol. 2003, 16, 1277.

3 Nishida, J.; Kawabata J.; J. Biosci. Biotechnol. Biochem. 2006, 70, 193.

4 Gacche, R. N.; Dhole, N. A.; Kamble, S. G.; Bandgar, B. P. J.; Enzyme Inhib. Med. Chem. 2007, 23, 28.

5 Vogel, S.; Ohmayer, S.; Brunner, G.; Heilmann, J.; Bioorg. Med. Chem. 2008, 16, 4286.

6 Jung, J.-C.; Jang, S.; Lee, Y.; Min, Y.; Lim, E.; Jung, H.; Oh, M.; Oh, S.; Jung, M.; J. Med. Chem. 2008, 51, 4054.

7 Sugamoto, K.; Kurogi, C.; Matsushita, Y.; Matsui, T.; Tetrahedron Lett. 2008, 49, 6639.

8 Avila, H. P.; Smania, E.; Monache, F. D.; Smania, A.; Bioorg. Med. Chem. 2008, 16, 9790.

9 Lawrence, N. J.; Patterson, R. P.; Ooi, L. L.; Cook, D.; Ducki, S.; Bioorg. Med. Chem. Lett. 2006, 16, 5844.

10 Cabrera, M.; Simoens, M.; Falchi, G.; Lavaggi, M. L.; Piro, O. E.; Castellano, E. E.; Vidal, A.; Azqueta, A.; Monge, A.; de Cerain, A. L.; Sagrera, G.; Seoane, G.; Cerecettoand, H.; Gonzalez, M.; Bioorg. Med. Chem. 2007, 15, 3356.
-1111 Boumendjel, A.; Boccard, J.; Carrupt, P.-A.; Nicolle, E.; Blanc, M.; Geze, A.; Choisnard, L.; Wouessidjewe, D.; Matera, E.-L.; Dumontet, C.; J. Med. Chem. 2008, 51, 2307. Their properties are related, among other factors, to its great conformational freedom as well as to the several patterns of substitution of A and B rings.1212 Larsen, M.; Kromann, H.; Kharazmi, A.; Nielsen, S. F.; Bioorg. Med. Chem. Lett. 2006, 15, 4858. Thus, the correct structural determination and the knowledge about three-dimensional (3D) atomic structure of chalcones are crucial to understand their properties. In this context, the complete assignments of 13C nuclear magnetic resonance (NMR) chemical shifts for chalcones, even using computation methodology, can help in identifying the chemical structures of new chalcone derivatives. The 2’-hydroxy-3,4,5-trimethoxy-chalcone (Figure 1) has shown in vitro inhibitory effects on PGE2 (prostaglandin E2) production from RAW 264.7 cells induced by LPS (lipopolysaccharide).1313 Tran, T.-D.; Park, H.; Kim, H. P.; Ecker, G. F.; Thai, K.-M.; Bioorg. Med. Chem. 2009, 19, 1650. At 10 μM 2’-hydroxychalcones proved strong effect to inhibit the PGE2 production (102.3%). However, these compounds also indicated the effect on cell viability with potential for anti-inflammation and anticancer activity.1313 Tran, T.-D.; Park, H.; Kim, H. P.; Ecker, G. F.; Thai, K.-M.; Bioorg. Med. Chem. 2009, 19, 1650.

Our research group has worked with quantum mechanics using computational tools to better understand the interpretation of the 13C NMR chemical shift experimental data. Recently, we published the applicability of the empirical equations for conversion of quantum mechanically calculated chemical shielding tensors into chemical shift values.1414 Costa, F. L. P.; de Albuquerque, A. C. F.; dos Santos Jr., F. M.; de Amorim, M. B.; J. Phys. Org. Chem. 2010, 23, 972.

15 Costa, F. L. P.; de Albuquerque, A. C. F.; Borges, R. M.; dos Santos Jr., F. M.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2014, 11, 219.
-1616 Costa, F. L. P.; de Albuquerque, A. C. F.; dos Santos Jr., F. M.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2015, 12, 2195. In these papers, it was pointed out the importance of a linear conversion formula in order to achieve great 13C NMR chemical shift experimental data reproduction and prediction. In this context, the goal of this work was to investigate the ability of the scaling factor protocol at the GIAO-mPW1PW91/6-31G(d)//mPW1PW91/6-31G(d) level of theory to predict the 13C NMR chemical shifts (d) of the 2’-hydroxy-3,4,5-trimethoxy-chalcone molecule (Figure 1). Moreover, two different approaches for determining the d of the 2’-hydroxy-3,4,5-trimethoxy-chalcone were compared.

Figure 1
Structure of the 2’-hydroxy-3,4,5-trimethoxy-chalcone.

Methodology

Chemicals

The 2’-hydroxy-3,4,5-trimethoxy-chalcone was obtained from Sigma-Aldrich Chemical Co.

NMR system and operating conditions

NMR analyses were acquired on a Bruker Avance III 11.75 Tesla spectrometer at 298 K using a 5 mm triple resonance broadband inverse (TBI) probehead. The spectra were obtained at 125.77 MHz for 13C using CDCl3 as the solvent. The 13C spectra (Figure S1, Supplementary Information (SI)) were acquired with spectral window of 37,878 Hz, 32,768 digitalized points and accumulation of 3,926 FIDs. The heteronuclear single-quantum correlation (HSQC, Figure S2, SI) and heteronuclear multiple-bond correlation (HMBC, Figure S3, SI) experiments were acquired with a spectral window of 10,000 and 37,731 Hz for 1H and 13C, respectively. The phase and baseline were corrected with the TopSpin software (version 3.2 Bruker BioSpin). NMR assignments are based on 1H,13C and 1H-13C HSQC/HMBC experiments. 1H NMR (500 MHz, CDCl3) d 12.85 (s, OH), 7.93 (dd, J 8.1, 1.6 Hz, H-6’), 7.84 (d, J 15.4 Hz, H-β), 7.53 (d, J 15.4 Hz, H-α), 7.50 (ddd, J 8.4, 7.2, 1.6 Hz, H-4’), 7.03 (dd, J 8.4, 1.2 Hz, H-3’), 6.95 (ddd, J 8.1, 7.2, 1.2 Hz, H-5’), 6.88 (s, H-2/H-6), 3.93 (s, OMe-3/OMe-5), 3.92 (s, OMe-4); HMBC/HSQC crosspeaks: OMe-3/OMe-5 (C-2, C-3, C-5, C-6), OMe-4 (C-4), H-2/H-6 (C-1, C-3, C-4, C-5, C-β), H-β (C-1, C-2, C-6, C-α, C=O), H-α (C-1, C=O), H-3’ (C-1’, C-2’, C-5’, C=O), H-4’ (C-2’, C-3’, C-5’, C-6’), H-5’ (C-2’, C-3’, C4’, C-6’), H-6’ (C-1’, C-2’, C-4’, C-5’, C=O).

Computational details

The 13C NMR chemical shifts of the three chalcones were calculated with two different approaches. In the first one, the absolute shieldings (σ) of all carbon atoms in each geometrically optimized conformer of molecule were calculated using the GIAO (gauge-independent atomic orbital) approximation at the mPW1PW91/6-31G(d)//mPW1PW91/6-31G(d) (NMR//optimization) level of theory, and further used to generate weighted average values for each atom considering the previously obtained conformational distribution, σaver. In the second one, only the σ for the lowest energetic conformer was taken to account, σlowe. Molecular mechanics calculations were performed using the Spartan’08 modeling software,1717 Spartan'08, Wavefunction Inc., Irvine, California, USA, 2010. whereas DFT (density functional theory) calculations were performed using the Gaussian 09 W software package.1818 Frisch, A. E.; Frisch, M. J.; Trucks, G.; Gaussian 09 User's Reference, Gaussian Inc., Pittsburgh, USA, 2009. Solvent effects were not taken into account in any calculation.

Statistical validation

In order to perform a statistical validation of our results the mean deviation (MD) and the root mean square deviation (RMSD) errors (in ppm) were calculated.

Results and Discussion

A randomized conformational search of the 2’-hydroxy-3,4,5-trimethoxy-chalcone molecule was performed (Figure 2) using the Monte Carlo (MC) method with a search limit of 200 structures. Merck molecular force field (MMFF) as implemented in the Spartan’08 software package, considering an initial energy cutoff of 10 kcal mol-1, was employed.

Figure 2
Superposition of the three lowest energy conformations of the 2’-hydroxy-3,4,5-trimethoxy-chalcone molecule.

The 28 more significant conformations of 2’-hydroxy-3,4,5-trimethoxy-chalcone molecule were saved, which are responsible for more than 99.99% of the total Boltzmann population in the first 10 kcal mol-1. This was followed by single-point energy calculations at the B3LYP/6-31G(d) and level of theory. The 21 more significant conformations within the range of 0.0-5.0 kcal mol-1 were selected by energy minimization calculations carried out at the mPW1PW91/6-31G(d) level of theory. The relevant results are given in Table 1. Frequency calculations carried out at the mPW1PW91/6-31G(d) level of theory confirmed the optimized geometries to be local minima and delivered values of free energy at 298 K and 1 atm. In this step, the three most significant conformations within the range of 0.0-3.0 kcal mol-1 were selected (Figure 3).

Table 1
DFT, thermo chemical parameter and DFT population for the 28 more significant conformations of 2'-hydroxy-3,4,5-trimethoxy-chalcone molecule

Figure 3
DFT-optimized structures of the three more stable conformers, 1, 2 and 3, respectively, of 2’-hydroxy-3,4,5-trimethoxy-chalcone. For relative energies and abundances of the conformers, see Table 1.

In the approach (I), for each optimized conformer geometry, the 13C atomic chemical shielding tensors (σ) were computed at the mPW1PW91/6-31G(d)//mPW1PW91/6-31G(d) level of theory. Isotropic atomic chemical shifts d in units of ppm were computed as differences between the atomic isotropic shielding of the solutes and corresponding reference atoms in tetramethylsilane (TMS). Thus, the population-averaged chemical shifts for the selected conformers were computed assuming Boltzmann statistics, according to equation 1, based on mPW1PW91/6-31G(d) free energies. Finally, the 13C NMR chemical shifts were scaled according to Costa et al. protocols.1515 Costa, F. L. P.; de Albuquerque, A. C. F.; Borges, R. M.; dos Santos Jr., F. M.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2014, 11, 219.

(1) δ = i δ i e Δ E i / k T i e Δ E i / k T

DEi is the relative energy of the ith conformer to the lowest energy, k is the Boltzmann constant, and the temperature T is set to 298 K. In the approach (II), only the lowest-energetic conformer was used to obtain the scaled chemical shifts.

Figure 3 shows the three most significant conformations conformers of 2’-hydroxy-3,4,5-trimethoxy-chalcone according to the geometry optimization calculations carried out at the mPW1PW91/6-31G(d) level of theory. The two more significant factors that determine the stability of the 2’-hydroxy-3,4,5-trimethoxy-chalcone conformers are apparently the intramolecular hydrogen bond between the OH and C=O and steric effects of methoxyl groups in the B-ring.

Table 2 shows the calculated (dcalc), scaled (dscal) and experimental (dexp) 13C NMR chemical shifts of the 2’-hydroxy-3,4,5-trimethoxy-chalcone as well as the differences (≠) among them. The data comparison demonstrated a great agreement between experimental and calculated NMR chemical shifts.

Table 2
Calculated (scaled) and differences between experimental 13C NMR chemical shifts (in ppm) in the 2'-hydroxy-3,4,5-trimethoxy-chalcone relative to TMS, for approaches (I) and (II)

The main differences between experimental and scaled chemical shifts, for both approaches, were observed in carbons C-4, OMe-3, OMe-4 and OMe-5. These differences are consistent with the three selected conformers (Figure 3) to obtain the chemical shifts at the approach (I). The reason is that they have different orientations of the OCH3 group (atoms C-3, C-4 and C-5) with respect to the aromatic ring and their chemical shifts are systematically dependent on the orientation of the methoxyl group.

For evaluating the robustness of the method, the chemical shift for two other chalcones: (E)-1-(4-hydroxy-3-methoxyphenyl)-3-(3,5-di-tert-butyl-4-hydroxyphenyl)prop-2-en-1-one and (E)-1-(4-aminophenyl)-3-(3,4-dimethoxyphenyl)prop-2-en-1-one (Figure S4, SI), were calculated. The results presented in Table 3 demonstrate an excellent predictive ability of the method. The experimental and theoretical calculation of the NMR data for these two chalcones are presented in the Supplementary Information.

Table 3
Statistical data illustrating the performance of two approaches for the GIAO calculation of 13C NMR chemical shifts (δ, in ppm) for the 2'-hydroxy-3,4,5-trimethoxy-chalcone (a), (E)-1-(4-hydroxy-3-methoxyphenyl)-3-(3,5-di-tert-butyl-4 hydroxyphenyl)prop-2-en-1-one (b) and (E)-1-(4-aminophenyl)- 3-(3,4-dimethoxyphenyl)prop-2-en-1-one (c)

Conclusions

In this work, we investigated the ability of the scaling factor protocol at the GIAO-mPW1PW91/6-31G(d)//mPW1PW91/6-31G(d) level of theory to predict the 13C NMR chemical shifts of the 2’-hydroxy-3,4,5-trimethoxy-chalcone. The robustness of the method was evaluated for two other chalcones: (E)-1-(4-hydroxy-3-methoxyphenyl)-3-(3,5-di-tert-butyl-4-hydroxyphenyl)prop-2-en-1-one and (E)-1-(4-aminophenyl)-3-(3,4-dimethoxyphenyl)prop-2-en-1-one, corroborating the ability of the method in chemical shift prevision. After using the scale factor, the two approaches were able to correctly reproduce the chemical shifts, despite the differences in the conditions of the experimental measurements and computational predictions. Both approaches had similar performance. These findings suggest that, in this case, either the approach (I) or the approach (II) could be chosen and that further analysis of the methoxychalcones is therefore justified.

Supplementary Information

Supplementary information is available free of charge at http://jbcs.sbq.org.br as PDF file.

Acknowledgments

The authors are grateful to Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Fundação de Amparo à Pesquisa do Estado de Goiás (FAPEG) and Financiadora de Estudos e Projetos (FINEP) for the financial support. The authors also thanks Prof Dr Guilherme Roberto de Oliveira for the chalcone samples.

References

  • 1
    Daskiewicz, J.-B.; Depeint, F.; Viornery, L.; Bayet, C.; Comte-Sarrazin, G.; Comte, G.; Gee, J. M.; Johnson, I. T.; Ndjoko, K.; Hostettmann, K.; Barron, D.; J. Med. Chem. 2005, 48, 2790.
  • 2
    Stevens, J. F.; Miranda, C. L.; Frei, B.; Buhler, D. R.; Chem. Res. Toxicol. 2003, 16, 1277.
  • 3
    Nishida, J.; Kawabata J.; J. Biosci. Biotechnol. Biochem. 2006, 70, 193.
  • 4
    Gacche, R. N.; Dhole, N. A.; Kamble, S. G.; Bandgar, B. P. J.; Enzyme Inhib. Med. Chem. 2007, 23, 28.
  • 5
    Vogel, S.; Ohmayer, S.; Brunner, G.; Heilmann, J.; Bioorg. Med. Chem. 2008, 16, 4286.
  • 6
    Jung, J.-C.; Jang, S.; Lee, Y.; Min, Y.; Lim, E.; Jung, H.; Oh, M.; Oh, S.; Jung, M.; J. Med. Chem. 2008, 51, 4054.
  • 7
    Sugamoto, K.; Kurogi, C.; Matsushita, Y.; Matsui, T.; Tetrahedron Lett. 2008, 49, 6639.
  • 8
    Avila, H. P.; Smania, E.; Monache, F. D.; Smania, A.; Bioorg. Med. Chem 2008, 16, 9790.
  • 9
    Lawrence, N. J.; Patterson, R. P.; Ooi, L. L.; Cook, D.; Ducki, S.; Bioorg. Med. Chem. Lett. 2006, 16, 5844.
  • 10
    Cabrera, M.; Simoens, M.; Falchi, G.; Lavaggi, M. L.; Piro, O. E.; Castellano, E. E.; Vidal, A.; Azqueta, A.; Monge, A.; de Cerain, A. L.; Sagrera, G.; Seoane, G.; Cerecettoand, H.; Gonzalez, M.; Bioorg. Med. Chem. 2007, 15, 3356.
  • 11
    Boumendjel, A.; Boccard, J.; Carrupt, P.-A.; Nicolle, E.; Blanc, M.; Geze, A.; Choisnard, L.; Wouessidjewe, D.; Matera, E.-L.; Dumontet, C.; J. Med. Chem. 2008, 51, 2307.
  • 12
    Larsen, M.; Kromann, H.; Kharazmi, A.; Nielsen, S. F.; Bioorg. Med. Chem. Lett. 2006, 15, 4858.
  • 13
    Tran, T.-D.; Park, H.; Kim, H. P.; Ecker, G. F.; Thai, K.-M.; Bioorg. Med. Chem. 2009, 19, 1650.
  • 14
    Costa, F. L. P.; de Albuquerque, A. C. F.; dos Santos Jr., F. M.; de Amorim, M. B.; J. Phys. Org. Chem. 2010, 23, 972.
  • 15
    Costa, F. L. P.; de Albuquerque, A. C. F.; Borges, R. M.; dos Santos Jr., F. M.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2014, 11, 219.
  • 16
    Costa, F. L. P.; de Albuquerque, A. C. F.; dos Santos Jr., F. M.; de Amorim, M. B.; J. Comput. Theor. Nanosci. 2015, 12, 2195.
  • 17
    Spartan'08, Wavefunction Inc., Irvine, California, USA, 2010.
  • 18
    Frisch, A. E.; Frisch, M. J.; Trucks, G.; Gaussian 09 User's Reference, Gaussian Inc., Pittsburgh, USA, 2009.

Publication Dates

  • Publication in this collection
    Nov 2017

History

  • Received
    26 Dec 2016
  • Accepted
    31 Mar 2017
Sociedade Brasileira de Química Instituto de Química - UNICAMP, Caixa Postal 6154, 13083-970 Campinas SP - Brazil, Tel./FAX.: +55 19 3521-3151 - São Paulo - SP - Brazil
E-mail: office@jbcs.sbq.org.br