A Study about Regioisomeric Hydroquinones with Multiple Intramolecular Hydrogen Bonding

A theoretical exploration about hydrogen bonding in a series of synthetic regioisomeric antitumor tricyclic hydroquinones is presented. The stabilization energy for the intramolecular hydrogen bond (IHB) formation in four structurally different situations were evaluated: (a) IHB between the proton of a phenolic hydroxyl group and an ortho-carbonyl group (forming a six-membered ring); (b) between the oxygen atom of a phenolic hydroxyl group and the proton of an hydroxyalkyl group (seven membered ring); (c) between the proton of a phenolic hydroxyl group with the oxygen atom of the hydroxyl group of a hydroxyalkyl moiety (seven-membered ring); and (d) between the proton of a phenolic hydroxyl group and an oxygen atom directly bonded to the aromatic ring in ortho position (five-membered ring). A conformational analysis for the rotation around the hydroxyalkyl substituent is also performed. It is observed that there is a correspondence between the conformational energies and the IHB. The strongest intramolecular hydrogen bonds are those involving a phenolic proton and a carbonyl oxygen atom, forming a six-membered ring, and the weakest are those involving a phenolic proton with the oxygen atom of the chromenone, forming five-membered rings. Additionally, the synthesis and structural assignment of two pairs of regioisomeric hydroquinones, by 2D-NMR experiments, are reported. These results can be useful in the design of biologically-active molecules.


Introduction
Hydrogen bonding (HB) plays a central role determining the structure, properties, and functions of biologically-important molecules [1][2][3]. HB interaction is the driving factor in many chemical and biochemical processes. Thus, studying HB in model compounds, which are structurally related to molecules fundamental for life, can be useful for the design of new biologically-active molecules. An intramolecular hydrogen bond (IHB) is formed when the proton donor and acceptor are in the same molecule [4]. An intermolecular hydrogen bond (IEHB) is formed when the donor and acceptor are in different molecules and, unlike IHB, IEHB energy is obtained easily from the energy difference between the two molecules together and separated. The non-bonded situation for IHB is not as easy to assess [5]. The identification of the close contact of an attractive IHB can sometimes be ambiguous, for example, those bonds found in a crystal. Experimental and theoretical tools, such as NMR and quantum chemical calculations, have been employed to evaluate the existence and strength of this kind of interaction in a series of tricyclic antitumor hydroquinones [5][6][7][8][9].
p-Hydroquinones (p-HQ) and their oxidized form, p-quinones (p-Q) are naturally occurring compounds, commonly isolated from different living organisms [10][11][12][13][14][15]. The biological activity of p-HQs has been related to their capability to form the semiquinone radical, through a one-electron oxidation. These intermediates have been associated to biological properties, such as pro-oxidant activity, by interacting with several intracellular molecules, such as DNA and proteins [16][17][18]. For example, in the inhibitors of mitochondrial complex I (NADH:ubiquinone oxidoreductase), the active compounds usually exhibit aromatic rings with a hydroquinone/quinone motif, planarity being one of the geometrical requirements. Three classes of inhibitors share a common binding domain in mitochondrial complex I (NADH:ubiquinone oxidoreductase), [19] and its interaction with complex I is highly sensitive to structural modifications [20]. In effect, small structural changes on hydroquinone scaffolds can determine the complex I inhibition or uncoupling of tumoral oxidative phosphorylation [21].
Modulation of the electron-transfer capability in p-HQs and p-Qs is an important variable to explore in order to obtain new biologically active compounds. Among the interactions that play a central role in this issue, the formation of inter-or intramolecular hydrogen bonds is of central importance in these molecules [22,23]. A recent study about the electrochemistry of hydroxy quinones possessing intramolecular hydrogen bonds (IHBs) shows that this interaction stabilizes the anion radical structure, leading to a shift in reduction potentials toward less negative values, making the process more spontaneous, when compared with hydroxy-quinones without IHBs [24]. IHBs have shown appreciable effects on the electrochemistry [25] and the antioxidant properties of hydroquinone and related phenols [26]. Antioxidant properties and free radical scavenging reactivity of a family of hydroxynaphthalenones and dihydroxyanthracenones were previously reported [27].
The o-carbonyl hydroquinone moiety is an important structural feature of several natural products with different biological activities, such as doxorubicin, daunorubicin, [28] acylnaphthohydroquinone [29,30], shikonin [30], and peyssonol A [31]. In previous works, we have shown that antioxidant molecules, containing an o-carbonyl hydroquinone moiety, inhibit mitochondrial tumor cell respiration [32][33][34]. Theoretical and biochemical insights about the activity of these hydroquinones have been recently described highlighting the role of intramolecular hydrogen bonding [35]. In this work, we carry out a conformational study of previously-reported hydroquinones 1,2 [36] and some new hydroquinone analogues, 3-6, (Figure 1) analyzing the different types of intramolecular hydrogen bonding (IHB) displayed and their influences on the conformational preferences of these molecules. An experimental (NMR) and a theoretical approach (MO) are employed in this report. The IHBs are studied by natural-bond-orbital (NBO) and atoms-in-molecules (AIM) analyses. Hydroquinone 3 was theoretically examined using diastereomers (5S, 8R, 1 S, and 1 R) and 4 was studied using diastereomers (5S, 8R, 1 S, and 1 R). The synthetic route used for obtaining compounds 1 and 2 was tested to know if compounds 3-6 are accessible experimentally. In effect, compounds 3-6 are obtained using the same methodology, though the respective isolated compounds were not stereochemically assigned.

Synthesis
In order to test if the synthetic route used in the synthesis of hydroquinones 1 and 2 is also adequate when a secondary dienol is used, a Diels-Alder reaction between quinone 7 and racemic (3E,5E)-hepta-3,5-dien-2-ol (8) (Figure 1, Scheme 1) was carried out. Though diasteromeric mixtures of regioisomers 3 and 4 must be obtained, only two major compounds were isolated. However, the stereochemical assignments of the isolated diasteromers were not carried out. Compounds 3 and 4 were obtained at 26% yield and 18% yield, respectively, in addition to hydroquinone 9 at 47% yield. The percentage of 9, a hydroquinone reduction product of 7, is very high when compared with the yield obtained with the primary dienol (15% yield) [37]. Racemic diene 8 was obtained by the reaction of hexadienal and methylmagnesium bromide.

Synthesis
In order to test if the synthetic route used in the synthesis of hydroquinones 1 and 2 is also adequate when a secondary dienol is used, a Diels-Alder reaction between quinone 7 and racemic (3E,5E)-hepta-3,5-dien-2-ol (8) (Figure 1, Scheme 1) was carried out. Though diasteromeric mixtures of regioisomers 3 and 4 must be obtained, only two major compounds were isolated. However, the stereochemical assignments of the isolated diasteromers were not carried out. Compounds 3 and 4 were obtained at 26% yield and 18% yield, respectively, in addition to hydroquinone 9 at 47% yield. The percentage of 9, a hydroquinone reduction product of 7, is very high when compared with the yield obtained with the primary dienol (15% yield) [37]. Racemic diene 8 was obtained by the reaction of hexadienal and methylmagnesium bromide. Dienophile 12 was synthesized following a described route starting from compound 10 (Scheme 2) [38].

Synthesis
In order to test if the synthetic route used in the synthesis of hydroquinones 1 and 2 is also adequate when a secondary dienol is used, a Diels-Alder reaction between quinone 7 and racemic (3E,5E)-hepta-3,5-dien-2-ol (8) (Figure 1, Scheme 1) was carried out. Though diasteromeric mixtures of regioisomers 3 and 4 must be obtained, only two major compounds were isolated. However, the stereochemical assignments of the isolated diasteromers were not carried out. Compounds 3 and 4 were obtained at 26% yield and 18% yield, respectively, in addition to hydroquinone 9 at 47% yield. The percentage of 9, a hydroquinone reduction product of 7, is very high when compared with the yield obtained with the primary dienol (15% yield) [37]. Racemic diene 8 was obtained by the reaction of hexadienal and methylmagnesium bromide. Dienophile 12 was synthesized following a described route starting from compound 10 (Scheme 2) [38]. Hydroquinones 5 and 6 were obtained by the same synthetic route used for compounds 3 and 4, a Diels-Alder reaction was performed between 2,2-dimethyl(2H)-chromen-4,5,8(3H)-trione (12) and Scheme 1. Synthesis of compounds 3 and 4.

Structural Assignments
Tricyclic hydroquinones 3-6 are closely related to the structure of a series of previously-studied tricyclic regioisomeric antitumor hydroquinones [32,34]. An ortho-carbonyl substituted hydroquinone derivative is an anticancer agent that acts by inhibiting mitochondrial bioenergetics and by inducing G2/M-phase arrest in mouse mammary adenocarcinoma TA3/Ha [21]. The complete assignments of the 13 C resonances by concerted use of Heteronuclear Multiple-Quantum Correlation (HMQC), Heteronuclear Multiple Bond Correlation (HMBC), and theoretical calculations have been previously reported [36,39]. Intramolecular hydrogen bonds in some of these o-carbonyl hydroquinones have been recently studied [40]. The key feature, in which NMR assignments are based, is the different nature of the phenolic hydroxyl groups. The hydroxyl group at C-9 is engaged in a strong IHB with the carbonyl oxygen atom, and this is evidenced by its downfield resonance at δ ranging from 11.34 to 13.17. The non-chelated phenolic hydroxyl protons, from compounds 1-4, resonate between 4.35 and 8.02 ppm. The non-chelated phenolic protons exhibit a chemical shift concentration-dependent behavior, because of the intermolecular interactions. The 10-OH resonances for compounds 5 and 6 are 5.05 and 5.70 ppm, respectively. Figure 2 shows the main HMBC correlations used for the assignments of hydroquinones 3-6.

Structural Assignments
Tricyclic hydroquinones 3-6 are closely related to the structure of a series of previously-studied tricyclic regioisomeric antitumor hydroquinones [32,34]. An ortho-carbonyl substituted hydroquinone derivative is an anticancer agent that acts by inhibiting mitochondrial bioenergetics and by inducing G2/M-phase arrest in mouse mammary adenocarcinoma TA3/Ha [21]. The complete assignments of the 13 C resonances by concerted use of Heteronuclear Multiple-Quantum Correlation (HMQC), Heteronuclear Multiple Bond Correlation (HMBC), and theoretical calculations have been previously reported [36,39]. Intramolecular hydrogen bonds in some of these o-carbonyl hydroquinones have been recently studied [40]. The key feature, in which NMR assignments are based, is the different nature of the phenolic hydroxyl groups. The hydroxyl group at C-9 is engaged in a strong IHB with the carbonyl oxygen atom, and this is evidenced by its downfield resonance at δ ranging from 11.34 to 13.17. The non-chelated phenolic hydroxyl protons, from compounds 1-4, resonate between 4.35 and 8.02 ppm. The non-chelated phenolic protons exhibit a chemical shift concentration-dependent behavior, because of the intermolecular interactions. The 10-OH resonances for compounds 5 and 6 are 5.05 and 5.70 ppm, respectively. Figure 2 shows the main HMBC correlations used for the assignments of hydroquinones 3-6. (2E,4E)-hexa-2,4-dien-1-ol (13), followed by enolization with silica gel (Scheme 3). This method allowed obtaining compounds 5 and 6 in 39% yield and 20% yield, respectively. Additionally, hydroquinone 11 was obtained with 36% yield. Scheme 3. Synthesis of compounds 5 and 6.

Structural Assignments
Tricyclic hydroquinones 3-6 are closely related to the structure of a series of previously-studied tricyclic regioisomeric antitumor hydroquinones [32,34]. An ortho-carbonyl substituted hydroquinone derivative is an anticancer agent that acts by inhibiting mitochondrial bioenergetics and by inducing G2/M-phase arrest in mouse mammary adenocarcinoma TA3/Ha [21]. The complete assignments of the 13 C resonances by concerted use of Heteronuclear Multiple-Quantum Correlation (HMQC), Heteronuclear Multiple Bond Correlation (HMBC), and theoretical calculations have been previously reported [36,39]. Intramolecular hydrogen bonds in some of these o-carbonyl hydroquinones have been recently studied [40]. The key feature, in which NMR assignments are based, is the different nature of the phenolic hydroxyl groups. The hydroxyl group at C-9 is engaged in a strong IHB with the carbonyl oxygen atom, and this is evidenced by its downfield resonance at δ ranging from 11.34 to 13.17. The non-chelated phenolic hydroxyl protons, from compounds 1-4, resonate between 4.35 and 8.02 ppm. The non-chelated phenolic protons exhibit a chemical shift concentration-dependent behavior, because of the intermolecular interactions. The 10-OH resonances for compounds 5 and 6 are 5.05 and 5.70 ppm, respectively. Figure 2 shows the main HMBC correlations used for the assignments of hydroquinones 3-6.
The hydroxymethyl group in compounds 1 ( Figure 4, Table 1) and 2 ( Figure 5, Table 2) exhibits a markedly different conformational behavior. When it is bonded to C-5 (compound 2) it exists as a single rotamer around the C5-CH 2 OH single bond, as evidenced by the very different values of the vicinal coupling constants of the methylene protons with HC-5, about 10 and 3 Hz, assigned to the anti and gauche protons, respectively. Their preferred conformation may be a consequence, in turn, of the preferred conformation of the hydroxyl group bonded to C-10, whose proton is oriented towards the gem-dimethyl group, thus avoiding the repulsion between the non-bonded electrons of oxygen atoms and the alkyl groups. Theoretical calculations support these considerations showing that this is the only conformation populated at room temperature. Compound 2 mainly shows two intramolecular hydrogen bonds, in which one of the phenolic protons are bonded to the carbonyl forming a six-membered ring, and the other is hydrogen bonded to the hydroxymethyl oxygen atom forming a seven-membered ring [36]. The hydroxymethyl group in compounds 1 ( Figure 4, Table 1) and 2 ( Figure 5, Table 2) exhibits a markedly different conformational behavior. When it is bonded to C-5 (compound 2) it exists as a single rotamer around the C5-CH2OH single bond, as evidenced by the very different values of the vicinal coupling constants of the methylene protons with HC-5, about 10 and 3 Hz, assigned to the anti and gauche protons, respectively. Their preferred conformation may be a consequence, in turn, of the preferred conformation of the hydroxyl group bonded to C-10, whose proton is oriented towards the gem-dimethyl group, thus avoiding the repulsion between the non-bonded electrons of oxygen atoms and the alkyl groups. Theoretical calculations support these considerations showing that this is the only conformation populated at room temperature. Compound 2 mainly shows two intramolecular hydrogen bonds, in which one of the phenolic protons are bonded to the carbonyl forming a six-membered ring, and the other is hydrogen bonded to the hydroxymethyl oxygen atom forming a seven-membered ring [36].     The hydroxymethyl group in compounds 1 ( Figure 4, Table 1) and 2 ( Figure 5, Table 2) exhibits a markedly different conformational behavior. When it is bonded to C-5 (compound 2) it exists as a single rotamer around the C5-CH2OH single bond, as evidenced by the very different values of the vicinal coupling constants of the methylene protons with HC-5, about 10 and 3 Hz, assigned to the anti and gauche protons, respectively. Their preferred conformation may be a consequence, in turn, of the preferred conformation of the hydroxyl group bonded to C-10, whose proton is oriented towards the gem-dimethyl group, thus avoiding the repulsion between the non-bonded electrons of oxygen atoms and the alkyl groups. Theoretical calculations support these considerations showing that this is the only conformation populated at room temperature. Compound 2 mainly shows two intramolecular hydrogen bonds, in which one of the phenolic protons are bonded to the carbonyl forming a six-membered ring, and the other is hydrogen bonded to the hydroxymethyl oxygen atom forming a seven-membered ring [36].      The hydroxymethyl group in compounds 1 ( Figure 4, Table 1) and 2 ( Figure 5, Table 2) exhibits a markedly different conformational behavior. When it is bonded to C-5 (compound 2) it exists as a single rotamer around the C5-CH2OH single bond, as evidenced by the very different values of the vicinal coupling constants of the methylene protons with HC-5, about 10 and 3 Hz, assigned to the anti and gauche protons, respectively. Their preferred conformation may be a consequence, in turn, of the preferred conformation of the hydroxyl group bonded to C-10, whose proton is oriented towards the gem-dimethyl group, thus avoiding the repulsion between the non-bonded electrons of oxygen atoms and the alkyl groups. Theoretical calculations support these considerations showing that this is the only conformation populated at room temperature. Compound 2 mainly shows two intramolecular hydrogen bonds, in which one of the phenolic protons are bonded to the carbonyl forming a six-membered ring, and the other is hydrogen bonded to the hydroxymethyl oxygen atom forming a seven-membered ring [36].     In compound 1, with the hydroxymethyl group at C-8, instead of C-5, the vicinal coupling constant exhibits values in the range 4-6 Hz, evidencing the presence of a mixture of rotamers. This fact shows that the hydrogen bond between a phenolic proton and the oxygen atom from the carbinol group, as observed for compounds 2, 4(R), and 4(S), is stronger than those bonded through a carbinol proton and a phenolic oxygen atom (compounds 1, 3(R), and 3(S)).
However, calculations on compound 3 show a markedly different behavior between epimers S ( Figure 6, Table 3) and R (Figure 7, Table 4) at carbinolic carbon. The analysis of the 3S rotamers show that 3Sa exhibits a hydrogen bond between the phenolic oxygen atom and the proton of the secondary hydroxyl group, and it has a population higher than 95% at 298 K. On the other hand, 3(R) conformers show very different results, where the 3Rc rotamer is the most stable one, with a 77% population, but without hydrogen bonding between the hydroxyalkyl substituent and the phenolic hydroxyl group, because that interaction imposes a steric repulsion between the methyl group, located at C5, and the methyl group of the hydroxyethyl substituent (1,4-diaxial repulsion in the boat like conformation of cyclohexane). The lower energy of the 3Rc conformer relies on the non-conventional stabilizing interaction between the hydroxyl oxygen atom and the methyl protons. Secondary alcohol analogs of 4 (epimers R and S at carbinolic center) show a similar conformational behavior, despite their absolute configuration (Figures 8 and 9, Tables 5 and 6). In compound 1, with the hydroxymethyl group at C-8, instead of C-5, the vicinal coupling constant exhibits values in the range 4-6 Hz, evidencing the presence of a mixture of rotamers. This fact shows that the hydrogen bond between a phenolic proton and the oxygen atom from the carbinol group, as observed for compounds 2, 4(R), and 4(S), is stronger than those bonded through a carbinol proton and a phenolic oxygen atom (compounds 1, 3(R), and 3(S)).
However, calculations on compound 3 show a markedly different behavior between epimers S ( Figure 6, Table 3) and R (Figure 7, Table 4) at carbinolic carbon. The analysis of the 3S rotamers show that 3Sa exhibits a hydrogen bond between the phenolic oxygen atom and the proton of the secondary hydroxyl group, and it has a population higher than 95% at 298 K. On the other hand, 3(R) conformers show very different results, where the 3Rc rotamer is the most stable one, with a 77% population, but without hydrogen bonding between the hydroxyalkyl substituent and the phenolic hydroxyl group, because that interaction imposes a steric repulsion between the methyl group, located at C5, and the methyl group of the hydroxyethyl substituent (1,4-diaxial repulsion in the boat like conformation of cyclohexane). The lower energy of the 3Rc conformer relies on the non-conventional stabilizing interaction between the hydroxyl oxygen atom and the methyl protons. Secondary alcohol analogs of 4 (epimers R and S at carbinolic center) show a similar conformational behavior, despite their absolute configuration (Figures 8 and 9, Tables 5 and 6).      In compound 1, with the hydroxymethyl group at C-8, instead of C-5, the vicinal coupling constant exhibits values in the range 4-6 Hz, evidencing the presence of a mixture of rotamers. This fact shows that the hydrogen bond between a phenolic proton and the oxygen atom from the carbinol group, as observed for compounds 2, 4(R), and 4(S), is stronger than those bonded through a carbinol proton and a phenolic oxygen atom (compounds 1, 3(R), and 3(S)).
However, calculations on compound 3 show a markedly different behavior between epimers S ( Figure 6, Table 3) and R (Figure 7, Table 4) at carbinolic carbon. The analysis of the 3S rotamers show that 3Sa exhibits a hydrogen bond between the phenolic oxygen atom and the proton of the secondary hydroxyl group, and it has a population higher than 95% at 298 K. On the other hand, 3(R) conformers show very different results, where the 3Rc rotamer is the most stable one, with a 77% population, but without hydrogen bonding between the hydroxyalkyl substituent and the phenolic hydroxyl group, because that interaction imposes a steric repulsion between the methyl group, located at C5, and the methyl group of the hydroxyethyl substituent (1,4-diaxial repulsion in the boat like conformation of cyclohexane). The lower energy of the 3Rc conformer relies on the non-conventional stabilizing interaction between the hydroxyl oxygen atom and the methyl protons. Secondary alcohol analogs of 4 (epimers R and S at carbinolic center) show a similar conformational behavior, despite their absolute configuration (Figures 8 and 9, Tables 5 and 6).          A different behavior takes place in chromenone derivatives 5 ( Figure 10, Table 7) and 6 ( Figure 11, Table 8). Compound 5 exhibits two hydrogen bonds in all calculated structures, the same already discussed before between one phenolic proton with the carbonyl oxygen atom forming a six membered ring, and a different hydrogen bond between the other phenolic proton and the oxygen atom of the chromenone ring, forming a five membered ring. The 1 H-NMR spectrum evidences these interactions by displaying chemical shifts of phenolic protons at δ = 5.05 and 11.34 for OH-5 and OH-10, respectively. The conformation of the hydroxymethyl group here has a similar behavior to that described for compound 1, displaying a mixture of three conformers. The preferred rotamer       A different behavior takes place in chromenone derivatives 5 ( Figure 10, Table 7) and 6 ( Figure 11, Table 8). Compound 5 exhibits two hydrogen bonds in all calculated structures, the same already discussed before between one phenolic proton with the carbonyl oxygen atom forming a six membered ring, and a different hydrogen bond between the other phenolic proton and the oxygen atom of the chromenone ring, forming a five membered ring. The 1 H-NMR spectrum evidences these interactions by displaying chemical shifts of phenolic protons at δ = 5.05 and 11.34 for OH-5 and OH-10, respectively. The conformation of the hydroxymethyl group here has a similar behavior to that described for compound 1, displaying a mixture of three conformers. The preferred rotamer  A different behavior takes place in chromenone derivatives 5 ( Figure 10, Table 7) and 6 ( Figure 11, Table 8). Compound 5 exhibits two hydrogen bonds in all calculated structures, the same already discussed before between one phenolic proton with the carbonyl oxygen atom forming a six membered ring, and a different hydrogen bond between the other phenolic proton and the oxygen atom of the chromenone ring, forming a five membered ring. The 1 H-NMR spectrum evidences these interactions by displaying chemical shifts of phenolic protons at δ = 5.05 and 11.34 for OH-5 and OH-10, respectively. The conformation of the hydroxymethyl group here has a similar behavior to that described for compound 1, displaying a mixture of three conformers. The preferred rotamer exhibits a hydrogen bond between the alkyl hydroxyl proton and the oxygen atom of the phenolic hydroxyl group. Compound 6 also exhibits the two hydrogen bonds of the phenolic protons described for regioisomer 5. In this case, practically two rotamers of the hydroxymethyl group are present, though with the hydrogen-bonded structure of the hydroxymethyl group as the most populated one. exhibits a hydrogen bond between the alkyl hydroxyl proton and the oxygen atom of the phenolic hydroxyl group. Compound 6 also exhibits the two hydrogen bonds of the phenolic protons described for regioisomer 5. In this case, practically two rotamers of the hydroxymethyl group are present, though with the hydrogen-bonded structure of the hydroxymethyl group as the most populated one.

NBO and AIM Analysis
The NBO and AIM analysis of results for 1 and 2 are given in Table 9 and Table 10, respectively. The main rotamers for each molecule were studied (>5% population). The results of NBO analysis for molecule 1 shows that rotamer 1a presents the highest values for the IHB O1 … HO9 stabilization energy, compared with 1b and 1c. In addition, 1a presents stabilization energy of 5.38 kcal/mol for the cooperative IHB O9 … HO1′, not present in the other rotamers. Regioisomer 2 displays only one major populated rotamer, 2a, which possess two IHB, O1 … HO9 with 30.59 kcal/mol and O1′ … HO10 with 16.26 kcal/mol. The strength of the last IHB is remarkably higher than the weakest IHB O9 … H-O1′ present in 1a, and it explains the difference in rotamer populations for compound 2 when compared with regioisomer 1. The higher acidity of a phenolic proton, and the concomitant lower basicity of a phenolic oxygen compared with the corresponding properties of the proton and oxygen in a hydroxyl group bonded to an alkyl chain can explain these results.
AIM analysis shows negative values for both, the Laplacian of electron density at the critical point (∇ 2 ρ) and the total electron energy density H (a.u.), which indicates the covalent nature of all IHB studied for 1 and 2. Hydrogen bond energy (EHB) obtained from the potential energy density (V)  exhibits a hydrogen bond between the alkyl hydroxyl proton and the oxygen atom of the phenolic hydroxyl group. Compound 6 also exhibits the two hydrogen bonds of the phenolic protons described for regioisomer 5. In this case, practically two rotamers of the hydroxymethyl group are present, though with the hydrogen-bonded structure of the hydroxymethyl group as the most populated one.

NBO and AIM Analysis
The NBO and AIM analysis of results for 1 and 2 are given in Table 9 and Table 10, respectively. The main rotamers for each molecule were studied (>5% population). The results of NBO analysis for molecule 1 shows that rotamer 1a presents the highest values for the IHB O1 … HO9 stabilization energy, compared with 1b and 1c. In addition, 1a presents stabilization energy of 5.38 kcal/mol for the cooperative IHB O9 … HO1′, not present in the other rotamers. Regioisomer 2 displays only one major populated rotamer, 2a, which possess two IHB, O1 … HO9 with 30.59 kcal/mol and O1′ … HO10 with 16.26 kcal/mol. The strength of the last IHB is remarkably higher than the weakest IHB O9 … H-O1′ present in 1a, and it explains the difference in rotamer populations for compound 2 when compared with regioisomer 1. The higher acidity of a phenolic proton, and the concomitant lower basicity of a phenolic oxygen compared with the corresponding properties of the proton and oxygen in a hydroxyl group bonded to an alkyl chain can explain these results.
AIM analysis shows negative values for both, the Laplacian of electron density at the critical point (∇ 2 ρ) and the total electron energy density H (a.u.), which indicates the covalent nature of all IHB studied for 1 and 2. Hydrogen bond energy (EHB) obtained from the potential energy density (V) Figure 11. Rotamers for 6.

NBO and AIM Analysis
The NBO and AIM analysis of results for 1 and 2 are given in Tables 9 and 10, respectively. The main rotamers for each molecule were studied (>5% population). The results of NBO analysis for molecule 1 shows that rotamer 1a presents the highest values for the IHB O1 . . . HO9 stabilization energy, compared with 1b and 1c. In addition, 1a presents stabilization energy of 5.38 kcal/mol for the cooperative IHB O9 . . . HO1 , not present in the other rotamers. Regioisomer 2 displays only one major populated rotamer, 2a, which possess two IHB, O1 . . . HO9 with 30.59 kcal/mol and O1 . . . HO10 with 16.26 kcal/mol. The strength of the last IHB is remarkably higher than the weakest IHB O9 . . . H-O1 present in 1a, and it explains the difference in rotamer populations for compound 2 when compared with regioisomer 1. The higher acidity of a phenolic proton, and the concomitant lower basicity of a phenolic oxygen compared with the corresponding properties of the proton and oxygen in a hydroxyl group bonded to an alkyl chain can explain these results.  NBO and AIM results for 3 and 4 are given in Tables 11 and 12, respectively. Epimer 3S presents a major rotamer 3Sa with a population higher than 95%. 3Sa has a high stabilization energy for the IHB O1 . .   AIM analysis reveals that rotamers of 5 and 6 present IHB of a covalent nature, similarly to the other studied molecules (Table 14). However, no bond critical points were found for IHB O1 . . . H-O10, which was identified in the NBO analysis with stabilization energy around 1 kcal/mol (weak IHB). The IHB O4 . . . H-O5, present in all rotamers studied for 5 and 6, present a lower AIM E HB value compared with their NBO stabilization energies. Table 14. Atoms-in-molecules parameters for 5 and 6. 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).

Molecule
Bond From these results, it is clear that the strongest intramolecular hydrogen bonds are those involving a phenolic proton and carbonyl oxygen atom forming a six-membered ring, and the weakest are those involving a phenolic proton with the oxygen atom of the chromenone ring, which forms five-membered rings. An intermediate situation is observed with hydrogen bonds formed by protons of alkyl hydroxyl groups and phenolic oxygen atoms, forming seven-membered rings. The size of the ring in the last two hydrogen bonds appears to be a main factor determining the energetics of these IHB.

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 ( 1 H) or 75.47 MHz ( 13 C). Measurements were carried out at a probe temperature of 300 K.
Cycloadition of quinone 7 with (3E,5E)-hepta-3,5-dien-2-ol (8). A solution of quinone 7 (132 mg, 0.65 mmol) and diene 8 (105 mg, 0.94 mmol) in C 6 H 6 (20 mL) was left at room temperature for ten days. A further addition of 2 g of silica gel and magnetic stirring overnight of the mixture allowed a solid mixture after filtration to be obtained, which was repeatedly washed with MeOH and the solvent removed. Column chromatography on silica gel with an 1:3 light petrol:EtOAc mixture allowed both hydroquinones to be separated. 9,10-Dihydroxy-8-

Computational Calculation
The calculations were carried out using the Gaussian 09 [42] program package. 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 [43] program package. AIM analysis were performed using the AIM2000 [44] program package.

Conclusions
A series of synthetic regioisomeric antitumor tricyclic hydroquinones were obtained, structurally characterized, and the multiple IHBs present theoretically studied. The synthesis of compounds 3, 4, 5, and 6 were carried out in a similar way than previously-reported compounds 1 and 2. Both yields and regioisomeric ratios for 4-5 and 5-6 were similar to those found in 1-2.
A ranking of the energy involved in the IHB formation between the proton of a phenolic hydroxyl group with an ortho-carbonyl group, forming a six-membered ring, in the IHB formation between the oxygen atom of a phenolic hydroxyl group and the proton of an hydroxyalkyl group, forming a seven-membered ring, in the IHB formation between the proton of a phenolic hydroxyl group with the oxygen atom of the hydroxyl group of a hydroxyalkyl group, forming a seven-membered ring and in the IHB between a proton of a phenolic hydroxyl group and an oxygen atom directly bonded to aromatic ring in ortho position, forming a five-membered ring, was constructed.
The strongest intramolecular hydrogen bonds are those involving a phenolic proton and carbonyl oxygen atom forming a six-membered ring, and the weakest are those involving a phenolic proton with the oxygen atom of the chromenone ring, which forms five-membered rings.
NBO and AIM analysis allowed us to rationalize, in most cases, the conformational preferences of these molecules in terms of the hyperconjugative donor-acceptor stabilization energies involved in the multiple IHBs present in each molecule. An exception is represented by 3R, where steric repulsions seem to predominate on the main rotamer, instead of IHB stabilization energies.
These results can be useful in the design of biologically-active molecules. Additionally, the two pairs of synthesized regioisomeric analogues of biologically-active hydroquinones 1 and 2, allowed the assessment of the possibility to experimentally access them.