Next Article in Journal
Elucidation of Vasodilation Response and Structure Activity Relationships of N2,N4-Disubstituted Quinazoline 2,4-Diamines in a Rat Pulmonary Artery Model
Next Article in Special Issue
Cooperative Effects in Weak Interactions: Enhancement of Tetrel Bonds by Intramolecular Hydrogen Bonds
Previous Article in Journal
Synthesis of Novel Pyrazole Derivatives and Their Tumor Cell Growth Inhibitory Activity
Previous Article in Special Issue
Exotic Spectra and Lattice Vibrations of Ice X Using the DFT Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Parameter Suitability for the Strength Evaluation of Intramolecular Resonance Assisted Hydrogen Bonding in o-Carbonyl Hydroquinones

by
Maximiliano Martínez-Cifuentes
1,*,
Matías Monroy-Cárdenas
2,
Juan Pablo Millas-Vargas
2,
Boris E. Weiss-López
3 and
Ramiro Araya-Maturana
2,4,*
1
Programa Institucional de Fomento a la Investigación, Desarrollo e Innovación, Universidad Tecnológica Metropolitana, Ignacio Valdivieso 2409, Casilla 9845, Santiago 8940577, Chile
2
Instituto de Química de Recursos Naturales, Universidad de Talca, Av. Lircay s/n, Casilla 747, Talca 3460000, Chile
3
Departamento de Química, Facultad de Ciencias, Universidad de Chile, Las Palmeras 3425, Casilla 653, Santiago 7800003, Chile
4
Programa de Investigación Asociativa en Cáncer Gástrico (PIA-CG), Universidad de Talca, Talca 3460000, Chile
*
Authors to whom correspondence should be addressed.
Molecules 2019, 24(2), 280; https://doi.org/10.3390/molecules24020280
Submission received: 22 November 2018 / Revised: 27 December 2018 / Accepted: 8 January 2019 / Published: 14 January 2019
(This article belongs to the Special Issue Intramolecular Hydrogen Bonding 2018)

Abstract

:
Intramolecular hydrogen bond (IMHB) interactions have attracted considerable attention due to their central role in molecular structure, chemical reactivity, and interactions of biologically active molecules. Precise correlations of the strength of IMHB’s with experimental parameters are a key goal in order to model compounds for drug discovery. In this work, we carry out an experimental (NMR) and theoretical (DFT) study of the IMHB in a series of structurally similar o-carbonyl hydroquinones. Geometrical parameters, as well as Natural Bond Orbital (NBO) and Quantum Theory of Atoms in Molecules (QTAIM) parameters for IMHB were compared with experimental NMR data. Three DFT functionals were employed to calculated theoretical parameters: B3LYP, M06-2X, and ωB97XD. OH distance is the most suitable geometrical parameter to distinguish among similar IMHBs. Second order stabilization energies ΔEij(2) from NBO analysis and hydrogen bond energy (EHB) obtained from QTAIM analysis also properly distinguishes the order in strength of the studied IMHB. ΔEij(2) from NBO give values for the IMHB below 30 kcal/mol, while EHB from QTAIM analysis give values above 30 kcal/mol. In all cases, the calculated parameters using ωB97XD give the best correlations with experimental 1H-NMR chemical shifts for the IMHB, with R2 values around 0.89. Although the results show that these parameters correctly reflect the strength of the IMHB, when the weakest one is removed from the analysis, arguing experimental considerations, correlations improve significantly to values around 0.95 for R2.

1. Introduction

Hydrogen bonding (HB) is generally described as a stabilizing electrostatic interaction with a partly covalent character. Some of the most direct experimental evidence of the HB formation is the deshielding of the hydrogen atom involved in the interaction, which is easily observed by NMR spectroscopy [1,2]. It is well known that intramolecular hydrogen bonding (IMBH) has a key role in determining the structure and properties of biologically-active molecules, and their study in model compounds is relevant in drug design [3,4,5]. For example, IMBH formation is determinant for passive membrane permeability of cyclic peptides [6] and in lipophilicity [7].
Any intramolecular hydrogen bonding IMHB is governed by the equilibrium between closed (IMHB) and open (no IMHB) forms. Consequently, the implementation of IMHB considerations in drug discovery programs requires experimental methodologies to assess the position of these equilibria, which are dependent on the strength of each IMBH [8]. The 1H-NMR spectroscopy is the most common tool to measure the propensity of compounds to form IMHBs, nevertheless, up to now the throughput of NMR, and other traditional methods, is limited in obtaining structure–activity relationships (SAR) and driving rapid cycles of compound optimization [8]. However, current computational chemistry allows one to assess, theoretically, the strength of IMBH through parameters such as the distance hydrogen-acceptor, the distance donor-acceptor, and the hydrogen bond angle [9]. Among the theoretical estimated parameters for the HB interaction, those from Natural Bond Orbital (NBO) and Quantum Theory of Atoms in Molecules (QTAIM) analysis [9,10,11,12,13] are the most popular. In this field, DFT methodologies have proven to be valuable tools to calculate molecular properties with an appreciably low computational time demand compared with ab initio methodologies.
Previously we have studied the properties of some ortho-carbonyl hydroquinones with capability to cross membranes and to reach the inner mitochondrial membrane. These compounds act either as oxidative phosphorylation (OXPHOS) uncouplers, complex I inhibitors, or exert both activities; as a consequence they exhibit antitumor [14,15,16] or antiplatelet effects [17], or both. Interestingly, small structural changes on the hydroquinone scaffold determine their biological activities [15]. All these molecules exhibit strong resonance assisted intramolecular hydrogen bonding, assessed by 1H-NMR spectrometry [18]. These compounds act either as weak acids into the mitochondrial membrane, or by losing a hydrogen atom (or alternatively losing an electron followed by deprotonation) to afford the corresponding semiquinone radical, which can interfere in cell signaling. The formation of IMHB plays a key role [19,20] on both properties. IMHBs have shown appreciable effects on the antioxidant and electrochemical properties of hydroquinones and related phenols [21,22].
On the other hand, the o-carbonyl hydroquinone motif is present on natural products with different biological activities, such as the anticancer drugs doxorubicin and daunorubicin [23], acylnaphtohydroquinones [24], shikonin [25], and in peyssonol A [26].
The precise assessment of the strength of the IMBH interaction has many difficulties. Determination of the H-bond energy is possible only for the intermolecular hydrogen bond: the energy of the intermolecular interaction is usually equal to the energy required to separate the two interacting molecules [1]. However, there is no clear definition of the hydrogen bond energy for the IMHB because the interacting entities cannot be fully separated [27].
Moreover, it is relatively easy to roughly calculate the relative strength of a series of IMBHs when they take place in molecules with large structural differences; however, to assess the relative strength of IMBHs in very similar molecules, establish correlations with experimental parameters, and use these data to make correlations with biological activities, requires an adequate choice of theoretical methodologies.
The aim of this work is to study the IMHBs experimentally and theoretically in a series of structurally similar o-carbonyl hydroquinones. Calculated NBO and QTAIM parameters of the IMHBs obtained by different DFT functionals (B3LYP, M06-2X, and ωB97XD) will be correlated with experimental values obtained by NMR.

2. Results

The studied compounds were synthesized according to Scheme 1 [28,29]. Starting from hydroquinone (1) or 2,3-dimethylhydroquinone (2), the corresponding acylhydroquinones (35) were obtained by Fries rearrangement [30]. These compounds were oxidized with silver (I) oxide to obtain the corresponding quinones (68), which were not isolated. These quinones yield benzofurans (1217) by reaction with the enamines (911). The last step involves an acid rearrangement of these benzofurans, yielding the corresponding bicyclic hydroquinones (18, 20, 2224). Monomethyl ethers of hydroquinone 19 and 21 were obtained by methylation of hydroquinones 18 and 20, respectively; they were synthesized to assess the effect of blocking the hydroxyl group not involved in IMHB.
Three hybrid functional were employed in this work: B3LYP, composed of Becke three parameters hybrid functional B3 and the non-local correlation functional of Lee, Yang, and Parr (LYP) [34,35,36]. The second corresponds to a hybrid meta-generalized gradient-approximation (hybrid meta-GGAs) functional, M06-2x, which have been described as one the of the best Truhlar group functionals for non-covalent interactions [37]. The last functional corresponds to the long-range corrected functional ωB97XD, which is designed to avoid rapid die-off of the non-coulomb part of the exchange functionals, which affects the modelling of long distance processes [38].
In this first part, we study main geometrical parameters from the IMHBs and their correlation with 1H-NMR chemical shifts. The results are presented in Table 1. Calculations were carried out at DFT level using the B3LYP, M06-2X, and ωB97XD functionals with the 6-311++g(d,p) basis set. Among the geometrical parameters, only the O15H16 distances displayed enough differentiation to be used as comparison parameters for correlation with the chemical shifts of chelated protons. Compound 21 displayed the largest chemical shift and a more de-shielded proton, indicating the strongest IMHB, followed in strength by 20. The later indicates that methyl substituents on C1 and C2 increase the IMHB strength. Additionally, the methylation at O12 and also a methyl group at C9 increase the strength of the IHB in 21 and 22, respectively. The three functionals give the shortest O15H16 distance for 21, indicating the strongest IHB, which is in agreement with the largest downfield chemical shift. The smallest chemical shift for 18 indicates the weakest IHB in the series. Nevertheless none of the three functional gives the largest O15H16 distance for 18. These discrepancies can reflect the presence of intermolecular interactions competing with the IMHB, contributing in a non-negligible degree to the 1H-NMR spectra at the probe temperature, especially between pairs of molecules where one phenolic hydroxyl is replaced by a methoxy group. The same considerations can be applied to correlations with the next parameters. Global correlation among the chemical shift and O15H16 distance for all molecules were carried out (Figure 1). Results indicate that the ωB97XD functional gives the best correlation with an R2 = 0.89, followed by M06-2X with a R2 = 0.85 and B3LYP, and the worst correlation with R2 = 0.80.
In order to estimate the energy of IMHB, the open-close method was employed, which defines the H-bond energy as the difference between the open and close conformers [39,40]. Calculations were carried out with the functional ωB97XD, which presents the best correlation between experimental chemical shift and O16H15 distance. Table 2 shows that energy for IMHB in the studied hydroquinones are around 15 kcal/mol. Hydroquinone 20 presents the highest value, 16.37 kcal/mol, while 19 presents the lowest IMHB energy, 14.22 kcal/mol. These energy differences are large enough to discard the contribution of the open form to the experimental NMR shieldings. Figure 2 shows the correlation between Hydrogen bond energy (EHB) and 1H-NMR chemical shifts (δ), which presents a R2 = 0.71, a value lower than the correlation with the O15H16 distance (R2 = 0.89), calculated with the same functional.
The results from NBO analysis are presented in Table 3. Second order stabilization energies ΔEij(2) between donor (φi) and acceptor orbital (φj) allow one to study orbital interactions responsible for the IMHB. For all studied compounds, the HB donor atom is always an oxygen (two lone-pairs) and the acceptor is the antibonding sigma orbital of an O-H groups (σ*OH). As we found in the geometrical analysis, strong IMHB presence in 21 (the largest chemical shift) is in agreement with the largest ΔEij(2) obtained from the calculations with the three functionals. Nevertheless, ΔEij(2) did not predict well the weakest IMHB, corresponding to 18, and the lowest stabilization energy was not predicted by any of the three calculations using the studied functionals. Global correlation among the chemical shift and ΔEij(2) for all HQs were carried out (Figure 3). As in the geometrical analysis, correlation results indicate that the best correlation is obtained with the ωB97XD functional, with R2 = 0.89. However, in this case it was followed by B3LYP with R2 = 0.87 and then M06-2X, the one that gave the worst correlation with R2 = 0.85.
The QTAIM analysis is presented in Table 4. It shows negative values for ∇2ρ and H for all studied compounds, indicating the covalent nature of their IHBs. |V|/G ratio is also a sensitive parameter for evaluating the covalency [41,42]. V can be interpreted as the pressure exerted on the electrons by the HB system, while G can be interpreted as the pressure exerted (as a reaction) by the electrons on the HB system, at the critical point. For negative values of ∇2ρ, as occurred in all our cases, values of |V|/G > 2 were indicative of covalent interaction [2]. Hydrogen bond energy (EHB) was obtained from the potential energy density (V) [13]. All values were slightly above 30 kcal/mol, while NBO analysis gave all energies slightly below 30 kcal/mol for the same interaction. These values were approximately twice those calculated by the open-close method, which indicates that NBO and AIM overestimate the interaction energy for the IMHB. EHB values calculated by different methods have shown to be generally inconsistent [43]. As we found in both geometrical and NBO analysis, the strongest IHB is obtained for compound 21 and is well represented by QTAIM analysis for the three functionals, giving the highest EHB for 21. On the other hand, similarly to NBO analysis, here QTAIM EHB also did not adequately reflect the weakest IMHB present in 18. Global correlations between QTAIM EHB and chemical shifts are presented in Figure 4. A better performance is shown by ωB97XD and B3LYP with a R2 = 0.88. As in the above cases, M06-2X gave the worst correlation with R2 = 0.85.
The results shown in Figure 3 corroborates that correlations with the weakest IMHB of this series are not well represented by the computational calculations employed. In order to test what is the impact of the weakest IMHB on the correlations studied previously, compound 18 was eliminated from the data and correlations were obtained again. The results showed a significant improvement, without change in the relative position among different functionals and parameters studied.
Figure 5 shows correlations between 1H-NMR chemical shifts and O11-H16O15 distance for the IHBs. R2 for B3LYP changed from 0.80 to 0.88, from 0.85 to 0.91 for M06-2x, and from 0.89 to 0.95 for ωB97XD.
Similarly, correlations between 1H-NMR chemical shifts and NBO energy for the IHBs are presented in Figure 6. R2 changed from 0.88 to 0.95 for B3LYP, from 0.85 to 0.92 for M06-2X, and from 0.89 to 0.96 for ωB97XD.
Finally, correlations between 1H-NMR chemical shifts and QTAIM energy for the IHBs, presented in Figure 7, showed a change in R2 from 0.88 to 0.94 for B3LYP, from 0.85 to 0.91 to M06-2X, and from 0.88 to 0.94 for ωB97XD.

3. Materials and Methods

3.1. Synthetic Procedures

Melting points are uncorrected. All NMR spectra were acquired using a Bruker AVANCE DRX 300 spectrometer (Bruker BioSpin GmbH, Rheinstetten, Germany) operating at 300.13 MHz (1H) or 75.47 MHz (13C). Measurements were carried out at a probe temperature of 300 K.

3.1.1. 1-(5-Hydroxy-3,3,6,7-tetramethyl-2-morpholino-2,3-dihydrobenzofuran-4-yl)ethan-1-one (14)

To a solution of 1-(2,5-dihydroxy-3,4-dimethylphenyl)ethan-1-one (400 mg, 2.2 mmol) in 20 mL of dichloromethane, Ag2O (1.28 g) was added and stirred with magnetic stirring for one hour. Then, the filtered solution was added over a solution of 4-(2-methylprop-1-en-1-yl)morpholine (314 mg, 2.2 mmol), in CH2Cl2 (10 mL), yielding 637 mg of the furan (14) (2.0 mmol, 90% yield). 1H-NMR δ(CDCl3): 1.37 (s, 3H, 3-CH3); 1.45 (s, 3H, 3-CH3); 2.14 (s, 3H, Ar-CH3); 2.19 (s, 3H, Ar-CH3); 2.33–2.60 (m, 4H, 2× CH2-O); 2.61 (s, 3H, CH3CO); 3.70 (m, 4H, 2× CH2-N); 4.68 (s, 1H, O-CH-N); 7.09 (s, 1H, 5-OH). 13C-NMR δ(CDCl3): 11.82; 12.51; 19.80; 31.17; 32.81; 45.07; 49.28; 66.75; 108.78; 120.94; 122.53; 123.74; 128.69; 147.01; 150.86; 205.43. M.P.: 140.1–142.2 °C. IR (KBr) ν/cm−1: 842.77; 1684.89; 2836.66; 2909.00; 2949.52; 2981.35; 3357.56. HRMS calculated: 319.1784; found: 319.1787.

3.1.2. 5,8-Dihydroxy-4,4,6,7-tetramethylnaphthalen-1(4H)-one (21)

A solution of furan 14 (451 mg, 1.4 mmol) in ethanol (15 mL) was refluxed for 2 h, then the solution was poured in a mixture ice and water, and the solid precipitate was filtered and dried yielding hydroquinone 21 (256 mg, 78% yield). 1H-NMR δ(CDCl3): 1.61 (s, 6H, 2 × 4-CH3); 2.23 (s, 3H, Ar-CH3); 2.26 (s, 3H, Ar-CH3); 4.55 (s, 1H, 5-OH); 6.23 (d, 1H, J = 10 Hz, 2-H), 6.82 (d, 1H, J = 10 Hz, 3-H), 13.20 (s, 1H, 8-OH). 13C-NMR δ(CDCl3): 11.56; 13.04; 25.17; 38.03; 112.68; 123.43; 124.04; 131.51; 131.64; 143.52; 155.12; 160.73; 191.11. P.F.: 200.5–201.6 °C. IR (KBr) ν/cm−1: 1658.84; 2963; 2998.71; 3406.75. HRMS calculated: 232.1099; found: 232.1113.

3.1.3. 8-Hydroxy-5-methoxy-4,4,6,7-tetramethylnaphthalen-1(4H)-one (22)

To a solution of hydroquinone 21 (100 mg, 0.43 mmol) in acetone (20 mL), powdered K2CO3 (85 mg, 0.6 mmol) and dimethyl sulfate (0.1 mL, 0.6 mmol) were added and stirred with magnetic stirring at reflux for 3 h. To the filtered mixture, a solution of KOH (5% in methanol) was added and kept at room temperature overnight. Then, it was neutralized with HCl 10% and extracted with CH2Cl2 5 × 20 mL. The organic phase was washed with water and dried with anhydrous sodium sulfate, then was filtered and the solvent was eliminated at vacuum, by purification in flash column chromatography with hexane/ethyl acetate (5:1) as eluent, the product 22 was obtained (79 mg, 74% yield) 1H-NMR δ(CDCl3): 1.57 (s, 6H, 2 × 4-CH3); 2.20 (s, 3H, Ar-CH3); 2.31 (s, 3H, Ar-CH3); 3.72 (s, 3H, 8-OCH3); 6.24 (d, 1H, J = 10 Hz, 2-H); 6.78 (d, 1H, J = 10 Hz, 3-H); 13.44 (s, 1H, 5-OH). 13C-NMR δ(CDCl3): 11.39; 14.77; 27.09; 38.15; 61.49; 112.36; 124.13; 124.84; 137.88; 139.43; 149.29; 157.10; 160.76; 190.92. P.F.: 82.1–84.3 °C. IR (KBr) ν/cm−1: 1655.95; 2848.23; 2920.58; 2972.67; 3001.61; 3042.12; 3441.48. HRMS Calculated: 246.1256; found: 246.1265.

3.2. Computational Calculations

The calculations were carried out using the Gaussian 09 [44] program package (Revision a.01; Gaussian, Inc.: Wallingford, CT, USA). No symmetry constraints were imposed on the optimizations, which were performed at the DFT B3LYP/6-311++G (d,p) level. No imaginary vibrational frequencies were found at the optimized geometries, indicating that they are true minima of the potential energy surface. NBO analysis was performed using the NBOPro 6.0 [45] program package (NBO 6.0; University of Wisconsin: Madison, WI, USA). QTAIM analysis was performed using the AIM2000 [46] program package (Büro für innovative software, Bielefeld, Germany).

4. Conclusions

In this work, we found that OH distance is the most suitable geometrical parameter for distinguishing among IMHBs of similar strengths. Second order stabilization energies ΔEij(2) from NBO analysis and hydrogen bond energy (EHB) obtained from QTAIM analysis were also able to properly distinguish among IMHBs in the studied molecule series. In all cases, ΔEij(2) from NBO give values for the IMHB lower than EHB obtained from QTAIM analysis. ΔEij(2) are slightly below 30 kcal/mol, while EHB from QTAIM gave values somewhat above 30 kcal/mol, which are approximately twice of those calculated by the open-close method. In all cases, parameters calculated using ωB97XD gave the best correlations with experimental 1H-NMR chemical shifts for the IMHBs, with values of R2 around 0.89. Despite this, the results showed that the employed parameters correctly described the strength of the IMHB, when the weakest one was removed from the analysis, correlations improved significantly to values around 0.95 for R2. These results can help to select the most suitable parameters for modeling molecules with IMHBs of similar strengths and achieving an adequate distinction between them through computational calculations.

Author Contributions

Conceptualization, M.M.-C. (Maximiliano Martínez-Cifuentes) and R.A.-M.; methodology, M.M.-C. (Maximiliano Martínez-Cifuentes) and R.A.-M.; software M.M.-C. (Maximiliano Martínez-Cifuentes) and M.M.-C. (Matías Monroy-Cárdenas); validation, M.M.-C. (Maximiliano Martínez-Cifuentes), R.A.-M., and B.E.W.-L.; formal analysis, M.M.-C. (Maximiliano Martínez-Cifuentes), R.A.-M., and B.E.W.-L.; investigation, M.M.-C. (Matías Monroy-Cárdenas), J.P.M-V., and R.A.-M.; writing—original draft preparation, M.M.-C. (Maximiliano Martínez-Cifuentes) and R.A.-M.; writing—review and editing, M.M.-C. (Maximiliano Martínez-Cifuentes), R.A-M., and B.E.W.-L.; project administration, M.M.-C. (Maximiliano Martínez-Cifuentes) and R.A.-M.; funding acquisition, M.M.-C. (Maximiliano Martínez-Cifuentes), R.A.-M., and B.E.W.-L.

Funding

This research was funded by Fondo Nacional de Desarrollo Científico y Tecnológico (FONDECYT), grant number 11170142 (M.M.-C. (Maximiliano Martínez-Cifuentes)) and 1180069 (R.A.-M.).

Acknowledgments

Powered@NLHPC: this research was partially supported by the supercomputing infrastructure of the NLHPC (ECM-02). R.A-M thanks Programa de Investigación Asociativa en Cáncer Gástrico (PIA-CG), Universidad de Talca, Talca.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arunan, E.; Desiraju, G.R.; Klein, R.A.; Sadlej, J.; Scheiner, S.; Alkorta, I.; Clary, D.C.; Crabtree, R.H.; Dannenberg, J.J.; Hobza, P.; et al. Defining the hydrogen bond: An account (IUPAC Technical Report). Pure Appl. Chem. 2011, 83, 1619–1636. [Google Scholar] [CrossRef] [Green Version]
  2. Grabowski, S.J. What Is the Covalency of Hydrogen Bonding? Chem. Rev. 2011, 111, 2597–2625. [Google Scholar] [CrossRef]
  3. Pairas, G.N.; Tsoungas, P.G. H-Bond: The Chemistry-Biology H-Bridge. Chemistryselect 2016, 1, 4520–4532. [Google Scholar] [CrossRef]
  4. Giordanetto, F.; Tyrchan, C.; Ulander, J. Intramolecular Hydrogen Bond Expectations in Medicinal Chemistry. ACS Med. Chem. Lett. 2017, 8, 139–142. [Google Scholar] [CrossRef] [PubMed]
  5. Kuhn, B.; Mohr, P.; Stahl, M. Intramolecular Hydrogen Bonding in Medicinal Chemistry. J. Med. Chem. 2010, 53, 2601–2611. [Google Scholar] [CrossRef]
  6. Rezai, T.; Bock, J.E.; Zhou, M.V.; Kalyanaraman, C.; Lokey, R.S.; Jacobson, M.P. Conformational flexibility, internal hydrogen bonding, and passive membrane permeability: Successful in silico prediction of the relative permeabilities of cyclic peptides. J. Am. Chem. Soc. 2006, 128, 14073–14080. [Google Scholar] [CrossRef] [PubMed]
  7. Alex, A.; Millan, D.S.; Perez, M.; Wakenhut, F.; Whitlock, G.A. Intramolecular hydrogen bonding to improve membrane permeability and absorption in beyond rule of five chemical space. Medchemcomm 2011, 2, 669–674. [Google Scholar] [CrossRef]
  8. Caron, G.; Vallaro, M.; Ermondi, G. High throughput methods to measure the propensity of compounds to form intramolecular hydrogen bonding. Medchemcomm 2017, 8, 1143–1151. [Google Scholar] [CrossRef]
  9. Afonin, A.V.; Sterkhova, I.V.; Vashchenko, A.V.; Sigalov, M.V. Estimating the energy of intramolecular bifurcated (three-centered) hydrogen bond by X-ray, IR and H-1 NMR spectroscopy, and QTAIM calculations. J. Mol. Struct. 2018, 1163, 185–196. [Google Scholar] [CrossRef]
  10. Shahi, A.; Arunan, E. Hydrogen bonding, halogen bonding and lithium bonding: An atoms in molecules and natural bond orbital perspective towards conservation of total bond order, inter- and intra-molecular bonding. Phys. Chem. Chem. Phys. 2014, 16, 22935–22952. [Google Scholar] [CrossRef]
  11. Tang, S.S.; Tsona, N.T.; Du, L. Ring-Size Effects on the Stability and Spectral Shifts of Hydrogen Bonded Cyclic Ethers Complexes. Sci. Rep.-UK 2018, 8, 11. [Google Scholar] [CrossRef] [PubMed]
  12. Feng, L.; Bai, F.Q.; Wu, Y.; Zhang, H.X. Does the molecular structure of CaH2 affect the dihydrogen bonding in CaH2 center dot center dot center dot HY (Y = CH3, C2H3, C2H, CN, and NC) complexes? A quantum chemistry study using MP2 and B3LYP methods. Sci. China Chem. 2012, 55, 262–269. [Google Scholar] [CrossRef]
  13. Martinez-Cifuentes, M.; Cardona, W.; Saitz, C.; Weiss-Lopez, B.; Araya-Maturana, R. A Study about Regioisomeric Hydroquinones with Multiple Intramolecular Hydrogen Bonding. Molecules 2017, 22, 593. [Google Scholar] [CrossRef] [PubMed]
  14. Urra, F.A.; Munoz, F.; Cordova-Delgado, M.; Ramirez, M.P.; Pena-Ahumada, B.; Rios, M.; Cruz, P.; Ahumada-Castro, U.; Bustos, G.; Silva-Pavez, E.; et al. FR58P1a; a new uncoupler of OXPHOS that inhibits migration in triple-negative breast cancer cells via Sirt1/AMPK/beta 1-integrin pathway. Sci. Rep.-UK 2018, 8, 16. [Google Scholar] [CrossRef]
  15. Urra, F.A.; Córdova-Delgado, M.; Lapier, M.; Orellana-Manzano, A.; Acevedo-Arévalo, L.; Pessoa-Mahana, H.; González-Vivanco, J.M.; Martínez-Cifuentes, M.; Ramírez-Rodróguez, O.; Millas-Vargas, J.P.; et al. Small structural changes on a hydroquinone scaffold determine the complex I inhibition or uncoupling of tumoral oxidative phosphorylation. Toxicol. Appl. Pharmacol. 2016, 291, 46–57. [Google Scholar] [CrossRef] [PubMed]
  16. Urra, F.A.; Martínez-Cifuentes, M.; Pavani, M.; Lapier, M.; Jaña-Prado, F.; Parra, E.; Maya, J.D.; Pessoa-Mahana, H.; Ferreira, J.; Araya-Maturana, R. An ortho-carbonyl substituted hydroquinone derivative is an anticancer agent that acts by inhibiting mitochondrial bioenergetics and by inducing G2/M-phase arrest in mammary adenocarcinoma TA3. Toxicol. Appl. Pharmacol. 2013, 267, 218–227. [Google Scholar] [CrossRef] [PubMed]
  17. Fuentes, M.; Araya-Maturana, R.; Palomo, I.; Fuentes, E. Platelet mitochondrial dysfunction and mitochondria-targeted quinone-and hydroquinone-derivatives: Review on new strategy of antiplatelet activity. Biochem. Pharmacol. 2018, 156, 215–222. [Google Scholar] [CrossRef]
  18. Dobado, J.A.; Gomez-Tamayo, J.C.; Calvo-Flores, F.G.; Martinez-Garcia, H.; Cardona, W.; Weiss-Lopez, B.; Ramirez-Rodriguez, O.; Pessoa-Mahana, H.; Araya-Maturana, R. NMR assignment in regioisomeric hydroquinones. Magn. Reson. Chem. 2011, 49, 358–365. [Google Scholar] [CrossRef]
  19. Salazar, R.; Vidal, J.; Martínez-Cifuentes, M.; Araya-Maturana, R.; Ramírez-Rodríguez, O. Electrochemical characterization of hydroquinone derivatives with different substituents in acetonitrile. New J. Chem. 2015, 39, 1237–1246. [Google Scholar] [CrossRef]
  20. Martinez-Cifuentes, M.; Salazar, R.; Escobar, C.A.; Weiss-Lopez, B.E.; Santos, L.S.; Araya-Maturana, R. Correlating experimental electrochemistry and theoretical calculations in 2’-hydroxy chalcones: The role of the intramolecular hydrogen bond. RSC Adv. 2015, 5, 50929–50937. [Google Scholar] [CrossRef]
  21. Rodriguez, J.; Olea-Azar, C.; Cavieres, C.; Norambuena, E.; Delgado-Castro, T.; Soto-Delgado, J.; Araya-Maturana, R. Antioxidant properties and free radical-scavenging reactivity of a family of hydroxynaphthalenones and dihydroxyanthracenones. Bioorgan. Med. Chem. 2007, 15, 7058–7065. [Google Scholar] [CrossRef] [PubMed]
  22. Li, X.C.; Chen, B.; Xie, H.; He, Y.H.; Zhong, D.W.; Chen, D.F. Antioxidant Structure-Activity Relationship Analysis of Five Dihydrochalcones. Molecules 2018, 23, 1162. [Google Scholar] [CrossRef] [PubMed]
  23. Lown, J.W. Molecular mechanisms of action of anticancer agents involving free radical intermediates. Adv. Free Radic. Biol. Med. 1985, 1, 225–264. [Google Scholar] [CrossRef]
  24. Pedroza, D.A.; De Leon, F.; Varela-Ramirez, A.; Lema, C.; Aguilera, R.J.; Mito, S. The cytotoxic effect of 2-acylated-1,4-naphthohydroquinones on leukemia/lymphoma cells. Bioorg. Med. Chem. 2014, 22, 842–847. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Yang, J.T.; Li, Z.L.; Wu, J.Y.; Lu, F.J.; Chen, C.H. An Oxidative Stress Mechanism of Shikonin in Human Glioma Cells. PLoS ONE 2014, 9, 12. [Google Scholar] [CrossRef] [PubMed]
  26. Loya, S.; Bakhanashvili, M.; Kashman, Y.; Hizi, A. Peyssonol-A and Peyssonal-B, 2 novel inhibitors of the reverse transcriptases of human-immunodeficiency-virus type-1 and type-2. Arch. Biochem. Biophys. 1995, 316, 789–796. [Google Scholar] [CrossRef] [PubMed]
  27. Arunan, E.; Desiraju, G.R.; Klein, R.A.; Sadlej, J.; Scheiner, S.; Alkorta, I.; Clary, D.C.; Crabtree, R.H.; Dannenberg, J.J.; Hobza, P.; et al. Definition of the hydrogen bond (IUPAC Recommendations 2011). Pure Appl. Chem. 2011, 83, 1637–1641. [Google Scholar] [CrossRef] [Green Version]
  28. Castro, C.G.; Santos, J.G.; Valcarcel, J.C.; Valderrama, J.A. Kinetic study of the acid-catalyzed rearrangement of 4-acetyl-3,3-dimethyl-5-hydroxy-2-morpholino-2,3-dihydrobenzo[b]furan. J. Org. Chem. 1983, 48, 3026–3029. [Google Scholar] [CrossRef]
  29. Araya-Maturana, R.; Cassels, B.K.; Delgado-Castro, T.; Valderrama, J.A.; Weiss-Lopez, B.E. Regioselectivity in the Diels-Alder reaction of 8,8-dimethylnaphthalene-1,4,5(8H)-trione with 2,4-hexadien-1-ol. Tetrahedron 1999, 55, 637–648. [Google Scholar] [CrossRef]
  30. Viglianisi, C.; Bartolozzi, M.G.; Pedulli, G.F.; Amorati, R.; Menichetti, S. Optimization of the Antioxidant Activity of Hydroxy-Substituted 4-Thiaflavanes: A Proof-of-Concept Study. Chem. Eur. J. 2011, 17, 12396–12404. [Google Scholar] [CrossRef]
  31. Araya-Maturana, R.; Delgado-Castro, T.; Garate, M.; Ferreira, J.; Pavani, M.; Pessoa-Mahana, H.; Cassels, B.K. Effects of 4,4-dimethyl-5,8-dihydroxynaphtalene-1-one and 4,4-dimethyl-5,8-dihydroxytetralone derivatives on tumor cell respiration. Bioorgan. Med. Chem. 2002, 10, 3057–3060. [Google Scholar] [CrossRef]
  32. Cassis, R.; Scholz, M.; Tapia, R.; Valderrama, J.A. Dienone-phenol rearrangement of naphthalenetriones-a route to 10-acetoxy-5,6,7,8-tetrahydrophenanthrene-1,4-diones. J. Chem. Soc. Perk. T 1987, 2855–2859. [Google Scholar] [CrossRef]
  33. Almodovar, I.; Ramirez-Rodriguez, O.; Barriga, A.; Rezende, M.C.; Araya-Maturana, R. Electrospray ionization mass spectrometric fragmentation of hydroquinone derivatives. Rapid. Commun. Mass. Spectrom. 2011, 25, 370–378. [Google Scholar] [CrossRef]
  34. Becke, A.D. Density-functional thermochemistry. 3. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  35. Lee, C.T.; Yang, W.T.; Parr, R.G. Development of the colle-salvetti correlation-energy formula into a functional of the electron-density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef]
  36. Miehlich, B.; Savin, A.; Stoll, H.; Preuss, H. Results obtained with the correlation energy density functionals of becke and Lee, Yang and Parr. Chem. Phys. Lett. 1989, 157, 200–206. [Google Scholar] [CrossRef]
  37. Zhao, Y.; Truhlar, D.G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar] [CrossRef]
  38. Chai, J.D.; Head-Gordon, M. Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615–6620. [Google Scholar] [CrossRef]
  39. Schuster, P.; Zundel, G.; Sandorfy, C. The Hydrogen Bond, Recent Development in Theory and Experiment; North Holland Publ. Co.: Amsterdam, The Netherland, 1976. [Google Scholar]
  40. Jablonski, M. Full vs. constrain geometry optimization in the open-closed method in estimating the energy of intramolecular charge-inverted hydrogen bonds. Chem. Phys. 2010, 376, 76–83. [Google Scholar] [CrossRef]
  41. Espinosa, E.; Molins, E.; Lecomte, C. Hydrogen bond strengths revealed by topological analyses of experimentally observed electron densities. Chem. Phys. Lett. 1998, 285, 170–173. [Google Scholar] [CrossRef]
  42. Espinosa, E.; Alkorta, I.; Elguero, J.; Molins, E. From weak to strong interactions: A comprehensive analysis of the topological and energetic properties of the electron density distribution involving X-H center dot center dot center dot F-Y systems. J. Chem. Phys. 2002, 117, 5529–5542. [Google Scholar] [CrossRef]
  43. Rusinska-Roszak, D. Energy of Intramolecular Hydrogen Bonding in ortho-Hydroxybenzaldehydes, Phenones and Quinones. Transfer of Aromaticity from ipso-Benzene Ring to the Enol System(s). Molecules 2017, 22, 481. [Google Scholar] [CrossRef] [PubMed]
  44. 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.; et al. Gaussian 09; revision A.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  45. Glendening, E.D.; Badenhoop, J.K.; Reed, A.E.; Carpenter, J.E.; Bohmann, J.A.; Morales, C.M.; Landis, C.R.; Weinhold, F. NBO 6.0; University of Wisconsin: Madison, WI, USA, 2013. [Google Scholar]
  46. Biegler-Konig, F.; Schonbohm, J. Update of the AIM2000-program for atoms in molecules. J. Comp. Chem. 2002, 23, 1489–1494. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Samples of the compounds are available from the authors.
Scheme 1. Synthesis of the bicyclic hydroquinones studied in this work.
Scheme 1. Synthesis of the bicyclic hydroquinones studied in this work.
Molecules 24 00280 sch001

HydroquinoneR1R2R3R4
18 [31]HHMeH
19a [31]HHMeMe
20MeHMeH
21bMeHMeMe
22 [32]HMeMeH
23 [32]HH-CH2-(CH2)-CH2-H
24 [15,33]HHEtH
Note: a Obtained by methylation of 18, b obtained by methylation of 20.
Figure 1. Correlations between H16O15 distance and 1H-NMR chemical shifts (δ) for the hydrogen bond.
Figure 1. Correlations between H16O15 distance and 1H-NMR chemical shifts (δ) for the hydrogen bond.
Molecules 24 00280 g001
Figure 2. Correlation between Hydrogen bond energy (EHB) and 1H-NMR chemical shifts (δ).
Figure 2. Correlation between Hydrogen bond energy (EHB) and 1H-NMR chemical shifts (δ).
Molecules 24 00280 g002
Figure 3. Correlations between second order stabilization energies ΔEij(2) LP O15 → σ* O11-H16 and 1H-NMR chemical shifts (δ) for IMHB.
Figure 3. Correlations between second order stabilization energies ΔEij(2) LP O15 → σ* O11-H16 and 1H-NMR chemical shifts (δ) for IMHB.
Molecules 24 00280 g003
Figure 4. Correlations between hydrogen bond energy EHB (kcal/mol) and experimental chemical shifts δ (ppm) of proton involved in IMHB.
Figure 4. Correlations between hydrogen bond energy EHB (kcal/mol) and experimental chemical shifts δ (ppm) of proton involved in IMHB.
Molecules 24 00280 g004
Figure 5. Correlations between O-HO distance and 1H-NMR chemical shifts (δ) for the IMHB, excluding the weakest value.
Figure 5. Correlations between O-HO distance and 1H-NMR chemical shifts (δ) for the IMHB, excluding the weakest value.
Molecules 24 00280 g005
Figure 6. Correlations between second order stabilization energies ΔEij(2) LP O15 → σ* O11-H16 and experimental chemical shifts (δ) of protons involved in IMHB, excluding the weakest value.
Figure 6. Correlations between second order stabilization energies ΔEij(2) LP O15 → σ* O11-H16 and experimental chemical shifts (δ) of protons involved in IMHB, excluding the weakest value.
Molecules 24 00280 g006
Figure 7. Correlations between hydrogen bond energy EHB (kcal/mol) and experimental chemical shifts δ (ppm) of protons involved in IMHB, excluding the weakest value.
Figure 7. Correlations between hydrogen bond energy EHB (kcal/mol) and experimental chemical shifts δ (ppm) of protons involved in IMHB, excluding the weakest value.
Molecules 24 00280 g007
Table 1. Experimental chemical shifts (δ) in ppm of protons involved in hydrogen bonding, and main geometric calculated parameters. Distances in angström, (Å), and angles in degree (°).
Table 1. Experimental chemical shifts (δ) in ppm of protons involved in hydrogen bonding, and main geometric calculated parameters. Distances in angström, (Å), and angles in degree (°).
Moleculeδ ppmDistance (O15H16)Distance (O11-H16)Angle (O15H16-O11)
B3LYPM062xωB97XDB3LYPM062xωB97XDB3LYPM062xωB97XD
1812.5401.6571.6941.6740.990.980.982148146147
1912.6901.6571.6951.6750.990.980.981148146147
2013.2001.6391.6621.6450.9920.9830.984149147148
2113.4401.6201.6581.6400.9940.9840.985149148149
2212.9301.6601.6911.6660.990.980.982148146147
2312.7401.6611.7011.680.990.980.982148146147
2412.8001.6541.691.6720.990.9810.982148146147
Table 2. Hydrogen bond energy (EHB), calculated as the difference between open and close conformers. Energies in kcal/mol.
Table 2. Hydrogen bond energy (EHB), calculated as the difference between open and close conformers. Energies in kcal/mol.
MoleculeEHB
1814.34
1914.22
2016.37
2115.71
2214.44
2314.36
2414.52
Table 3. Second order stabilization energies (ΔEij(2)) between donor LP O15 and acceptor σ* O11-H16 orbitals involved in the IMHB. Energies in kcal/mol.
Table 3. Second order stabilization energies (ΔEij(2)) between donor LP O15 and acceptor σ* O11-H16 orbitals involved in the IMHB. Energies in kcal/mol.
MoleculeΔEij(2) LP O15 → σ* O11-H16
B3LYPM06-2XωB97XD
1827.1423.1225.11
1927.0222.9524.87
2028.9926.1327.99
2131.1526.628.62
2227.6523.2325.68
2326.6722.5424.59
2427.4723.4425.27
Table 4. Atoms-in-molecule parameters for IHB of 1824. Electron density at the critical point ρBCP (a.u), its Laplacian ∇2ρ (a.u.), electron kinetic energy density G (a.u.), potential energy density V (a.u.), total electron energy density H (a.u.), and hydrogen bond energy EHB (kcal/mol).
Table 4. Atoms-in-molecule parameters for IHB of 1824. Electron density at the critical point ρBCP (a.u), its Laplacian ∇2ρ (a.u.), electron kinetic energy density G (a.u.), potential energy density V (a.u.), total electron energy density H (a.u.), and hydrogen bond energy EHB (kcal/mol).
FunctionalMoleculeρBCP2ρGV|V|/GHEHB
B3LYP180.0537−0.03770.0459−0.10132.2070−0.147231.91
190.0537−0.03780.0460−0.10152.2065−0.147531.96
200.0560−0.03830.0478−0.10512.1987−0.152933.12
210.0587−0.03920.0501−0.11002.1956−0.160134.65
220.0546−0.03830.0469−0.10332.2025−0.150232.54
230.0533−0.03740.0455−0.10032.2044−0.145731.59
240.0541−0.03780.0463−0.10202.2030−0.148332.14
M06-2X180.0454−0.04080.0444−0.09902.2297−0.143431.19
190.0454−0.04090.0444−0.09902.2297−0.143431.19
200.0490−0.04270.0479−0.10662.2255−0.154533.57
210.0495−0.04280.0483−0.10742.2236−0.155733.83
220.0458−0.04120.0449−0.10012.2294−0.145031.53
230.0446−0.04020.0435−0.09702.2299−0.140530.57
240.0459−0.04100.0448−0.09992.2299−0.144831.48
ωB97XD180.0510−0.03800.0447−0.09892.2125−0.143531.14
190.0508−0.03810.0446−0.09872.2130−0.143431.10
200.0546−0.03940.0479−0.10572.2067−0.153633.30
210.0553−0.03960.0485−0.10692.2041−0.155333.66
220.0518−0.03860.0457−0.10102.2100−0.146631.81
230.0502−0.03760.0439−0.09732.2164−0.141230.64
240.0512−0.03810.0449−0.09932.2116−0.144231.29

Share and Cite

MDPI and ACS Style

Martínez-Cifuentes, M.; Monroy-Cárdenas, M.; Millas-Vargas, J.P.; Weiss-López, B.E.; Araya-Maturana, R. Assessing Parameter Suitability for the Strength Evaluation of Intramolecular Resonance Assisted Hydrogen Bonding in o-Carbonyl Hydroquinones. Molecules 2019, 24, 280. https://doi.org/10.3390/molecules24020280

AMA Style

Martínez-Cifuentes M, Monroy-Cárdenas M, Millas-Vargas JP, Weiss-López BE, Araya-Maturana R. Assessing Parameter Suitability for the Strength Evaluation of Intramolecular Resonance Assisted Hydrogen Bonding in o-Carbonyl Hydroquinones. Molecules. 2019; 24(2):280. https://doi.org/10.3390/molecules24020280

Chicago/Turabian Style

Martínez-Cifuentes, Maximiliano, Matías Monroy-Cárdenas, Juan Pablo Millas-Vargas, Boris E. Weiss-López, and Ramiro Araya-Maturana. 2019. "Assessing Parameter Suitability for the Strength Evaluation of Intramolecular Resonance Assisted Hydrogen Bonding in o-Carbonyl Hydroquinones" Molecules 24, no. 2: 280. https://doi.org/10.3390/molecules24020280

APA Style

Martínez-Cifuentes, M., Monroy-Cárdenas, M., Millas-Vargas, J. P., Weiss-López, B. E., & Araya-Maturana, R. (2019). Assessing Parameter Suitability for the Strength Evaluation of Intramolecular Resonance Assisted Hydrogen Bonding in o-Carbonyl Hydroquinones. Molecules, 24(2), 280. https://doi.org/10.3390/molecules24020280

Article Metrics

Back to TopTop