Next Article in Journal
Design of κ-Opioid Receptor Agonists for the Development of Potential Treatments of Pain with Reduced Side Effects
Previous Article in Journal
One-Pot Hydrothermal Preparation of Hydroxyapatite/Zinc Oxide Nanorod Nanocomposites and Their Cytotoxicity Evaluation against MG-63 Osteoblast-like Cells
Previous Article in Special Issue
Impact of Dispersion Force Schemes on Liquid Systems: Comparing Efficiency and Drawbacks for Well-Targeted Test Cases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prototropy, Intramolecular Interactions, Electron Delocalization, and Physicochemical Properties of 1,8-dihydroxy-9-anthrone—DFT-D3 Study of Substituent Effects

by
Małgorzata Szymańska
* and
Irena Majerz
*
Faculty of Pharmacy, Wroclaw Medical University, Borowska 211a, 50-556 Wroclaw, Poland
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(1), 344; https://doi.org/10.3390/molecules28010344
Submission received: 17 August 2022 / Revised: 27 December 2022 / Accepted: 28 December 2022 / Published: 1 January 2023
(This article belongs to the Special Issue Molecular Sensitivity and Weak Interactions)

Abstract

:
1,8-dihydroxy-9-anthrone are tricyclic compounds with a ketone group in the middle ring and two hydroxyl groups substituted in the side-aromatic rings what results in formation of two intramolecular hydrogen bonds in which the oxygen atom from the ketone group is the proton acceptor. 1,8-dihydroxy-9-anthrones in which intramolecular proton transfer between C10 and CO in the middle ring occurs, can exist in a tautomeric keto-enol equilibrium. For anthralin, the most important representative of this group, this equilibrium has been studied previously, but it has not been studied for its derivatives. Substituents in the middle ring change the geometry of 1,8-dihydroxy-9-anthrones so they are also expected to affect the keto-enol equilibrium. It is also important to study the effect of intramolecular hydrogen bonds on the structure of both tautomeric forms. It was found that the nature of the substituent in the middle ring could affect the antioxidant properties of the investigated compound.

1. Introduction

In the previous work [1], we studied the influence of substituent on the structure of anthrones and anthraquinones-tricyclic compounds with a wide importance in biological processes. The change of the angle between the anthrone aromatic rings is associated with the change in electron density at the ring critical point of the central ring.
An important group of anthrone derivatives are 1,8-dihydroxy-9-anthrones with substituents in the 1 and 8 position (Scheme 1). Anthralin, the popular drug from the 1,8-dihydroxy-9-anthrones group, has two hydroxyl groups at position 1 and 8 on either side of the ketone group located at position 9 (Scheme 1) and is used primarily in the treatment of psoriasis [2]. In the ketone form, it is stabilized by hydrogen bonds which are formed between the OH groups and the ketone oxygen [3]. For this reason, according to Hellier et al. [4], anthraline exists in its entirety in the ketone form [4]. The release of the hydrogen atom from the C10 position initiates the formation of free radicals, and although it may have an effect on skin irritation, it can be important in the mechanism of action of the drug [3,5]. Anthralin reduction power can be directly modified by changes in the structure [6]. For this reason, scientists are looking for new analogues substituted at the C10 position, which will prove to be more effective and, additionally, will not have side effects [5,7].
Intramolecular interactions of polycyclic compounds affect their structure [8,9,10]. The structure of a compound is related to the pharmacological activity; therefore, it is important to understand the effects of substituents, intramolecular hydrogen bonds, and keto-enol equilibrium on the structures of 1,8-dihydroxy-9-anthrone derivatives.
Two side rings of the 1,8-dihydroxy-9-anthrones derivatives are aromatic. The central ring in the lowest-energy ketone structure contains one CH2 group, which influences the aliphatic character of the ring. However, the remaining fragment of the central ring further delocalized as a result of π-π cross conjugation [11] and affects the cross delocalization of the other rings. The sensitivity of the middle ring to substitution can shift the character of the middle ring toward aromaticity, which can be affected by its non-planarity [1]. The presence of a ketone group in the middle ring and two adjacent OH groups in the side aromatic rings at the 1 and 8 positions causes the hydrogen atoms to form intramolecular hydrogen bonds. Migration of the proton between these ketone and hydroxyl group results in formation of tautomers [12,13,14,15]. Compounds, in which the intramolecular proton transfer takes place, are important in medicine and pharmacy and designing new drugs [16,17,18,19].
Keto-enol equilibrium can exist in 1,8-dihydroxy-9-anthrones through to a single proton transfer from C10 to CO in the middle ring (Scheme 2). According to Marrero-Carballo et al. [15], for proton transfer to occur the central ring must be twisted into a boat-like conformation. This equilibrium is influenced by the solvents in which the substance is dissolved [13,20,21], the substituents [21,22,23,24], and the temperature [21]. Baba and Takemura [13] studied the keto-enol equilibrium between anthrone and anthranol-1 in isooctane using the spectrophotometric method. After dropping the anthrone in isooctane, no changes in the spectrum were observed even after two days. The reaction was much faster in the presence of a small amount of triethylamine, which is basic in nature [13]. Using computational methods, Marrero-Carballo et al. [15] compared the activation energy of intramolecular proton transfer in the chrysophanol molecule with the proton transfer to the pyridine molecule. They found that, in the case of intramolecular proton transfer in the anthrone molecule, the activation energy is much greater than in the case of proton transfer to the pyridine molecule. The higher activation energy is related to the deformation of the central ring in the anthrone molecule [15].
Laurella et al. [21] studied the effects of substituents on the keto-enol balance of β-ketoamides. They found that electron-withdrawing properties of chlorine atom causes an equilibrium shift towards the enol form, while methoxy groups that donate the electrons increases the content of the keto form. Additionally, they noticed that the equilibrium was also influenced by the position of the substituent, which is related to intramolecular hydrogen bonds in the molecule [21].
It is important to study the antioxidant properties of the analyzed compounds. Antioxidants have the ability to neutralize the harmful effects of free radicals. One of the methods of predicting the antioxidant properties is the determination of the enthalpy of the OH bond dissociation (BDE) and so the nature of the substituents in the ring affects the BDE value [25]. In the enol form, dissociation of the O-H bond at C9 is probable, which may have a positive effect on antioxidant properties. According to Lucarini et al. [26], electron-positive substituents decrease the BDE of O-H, while electron withdrawing substituents increase the BDE of O-H value relatively to unsubstituted phenol [26].
This work is a continuation of our previous [1] work in which we studied the influence of substituents on the structure of anthrones and anthraquinones and an analysis of 50 optimized compounds was performed. The substituents used in that work were differed in size, electron donating and electron accepting properties: NO2, CHO, COOH, CH3, CH2CH3, NH2, OH, Cl, C(CH3)3. In this work, each structure with a specific substituent has been optimized in the ketone form and in the enol form and the energy difference between them has been calculated. The lowest energy has been obtained for the keto structure with two hydrogen bonds and has been used as the reference energy. The electron density of the middle ring is sensitive to substitution as mentioned in the previous work [1,27] so also in this work it is used as a measure of substituent influence and aromaticity of the central ring. The process of intramolecular proton transfer is responsible for acid-base regulation in the cell or for enzymatic reactions, so it is important to study the effect of substituents on this reaction [12,15]. Moreover, the intramolecular hydrogen bonds have an influence on the geometry investigated compounds. If the substituents and hydrogen bonds affect the structure, then it is important to study their effect on the keto-enol equilibrium, which is the aim of the research undertaken.

Computational Details

Optimization of the 1,8-dihydroxy-9-anthrone derivatives was performed with a Gaussian 16 package [28] at DFT-D3 B3LYP/6-311++G** level [29,30]. Grimme dispersion [31] was included to reproduce correctly the hydrogen bond in the investigated molecules. Vibrational frequencies were calculated to confirm that the optimized molecule reached the minimum of energy. QTAIM parameters were calculated with the AIMALL program [32] using the wave function for the optimized molecule.

2. Results

2.1. 1,8-Dihydroxy-9-anthrone Tautomers

To investigate the tautomeric preference, the Gibbs energy of different isomers in ketone form for four different substituents NH2, NO2, OH, and H has been compared (Table 1). There are six possible different tautomeric structures of the ketone form (Scheme 3). For NH2, NO2, H, and OH substituted structures, five isomers with different energies have been obtained. In all cases, the K3 isomer converges to the K6 isomer, which is confirmed by the C9-O (1.261 Å), C9-C9a (1.462 Å), 4a-9a (1.412 Å), and 8a-10a (1.412 Å) bond lengths equal and respectively identical as these bond length for the K6 tautomer. The lowest value of the Gibbs energy for the K6 isomer has been obtained. Isomers with relative Gibbs energies higher than 10 kcal/mol can be neglected in the isomeric mixture owing to their exceptionally low percentage contents (<5 ppm). The smallest differences in Gibbs energy of 10 kcal/mol is obtained between isomers of NH2 and OH substituted structures. In all cases, the α angle is the greatest for the lowest energy structures. The HOMA values in the side rings of the preferred isomeric forms are always close to 1.

2.2. Crystal Structures

It can be expected that the tautomerism of dihydroxyantrone will be reflected in the structures of compounds deposited in the CSD database [33]. Since the structure of the model compounds is known [34], the analysis of the bond lengths in the 1,8-dihydroxy-9-anthrone derivatives in the crystalline state should allow to indicate which tautomer is realized in the crystal of individual derivatives. Table S1 summarizes the C-C bond lengths for the 1,8-dihydroxy-9-anthrone derivatives present in the CSD database. The table also includes CO bond lengths, HOC angles, and torsion angles allowing for the indication of proton deviation from the plane defined by the system of intramolecular hydrogen bonds typical of 1,8-dihydroxy-9-anthrones. The geometrical parameters should indicate which tautomeric forms are realized in the crystalline state.
The analysis of the CO bond system allows for a preliminary indication of the tautomeric form present in the crystal structures. In the structures K1, K2, and K3, the double bond of CO occurs in the side ring. The comparison of CO bonds in Table S1 clearly shows that all solid state structures can be of the K4, K5, or K6 type, in which the CO double bond connects oxygen to the central ring. According to the model structures (Scheme 3), individual tautomeric forms should differ from each other in the arrangement of single and double bonds. The analysis of the CC bond lengths in the 1,8-dihydroxy-9-anthrone rings allows for the unequivocal elimination of the K5 tautomer, because all C3-C4 bonds are longer than those typical for benzene and are typical single bond. Since all the C2-C3 lengths in the analyzed compounds are shorter than the length typical for benzene (1.399) and the C1-C9a bond length is also shortened for a number of compounds, it can be considered that this structure is typical for solid-state 1,8-dihydroxy-9-anthrones. The analysis of C4a-C9a and C8a-C10a bond lengths indicates their shortening below the value typical for benzene in some dihydroxyantrone derivatives, which proves the possibility of the presence of the K6 tautomeric structure, characterized by the lowest energy. The linkage lengths of the dihydroxyanntron derivatives therefore indicate the possible presence of the K6 tautomer, with the K4 tautomer being more likely. However, the precision in determining the bond length should be taken into account. The comparison of the bond lengths in the central ring clearly shows that the K6 isomer dominates in the 1,8-dihydroxy-9-anthrone crystals. The shortest bonds, similar in length to benzene, are the bonds 4a-9a and 8a-10a.
It is a noteworthy fact that, in some structures, the protons of the hydroxyl groups remain unconnected with the carbonyl oxygen of the middle ring (BOLPEX [35], CARMYC11 [36], DHANQU03 [37], DHANQU04 [37], DHANQU08 [37], JUKREM [38], PIRFIH [39], QEGXUV [40], VURHEV [41]), as confirmed by the values of the torsion angles in Table S1. The deviation of the proton from the plane convenient for the formation of intramolecular hydrogen bonds is due to the participation of the oxygen atoms in other interactions resulting from intermolecular interactions in the crystal.

2.3. Geometry of the Investigated Compounds

The α angles and relative energies for the optimized structures with different substituents are collected in Table 2. The lowest energy structure in the ketone form has been obtained for the K6 tautomer and remaining discussions in the study apply only to it. The Ea structure represents the enol form with one hydrogen bond and Eb with two hydrogen bonds. The Ka structure represents ketone form without hydrogen bonds, Kb with one hydrogen bond, and Kc with two hydrogen bonds. As in the previous works [1,27], the angle between the planes of the two side aromatic rings has been used as a measure of the substituent influence on the geometry of the ring system (Scheme 4. Depending on the character of the substituent and its size as well as the presence of hydrogen bonds, the α angle changes from 0 to 52°. The H···O distance in the hydrogen bonds is close to 1.7 Å. Molecules in the ketone form, without hydrogen bonds, are characterized by the highest α angle. The angle decreases with the increasing number of the hydrogen bonds in the molecule. Deviations for this rule were observed for the structures of 9K6c, 9K6b and 5K6b, 5K6a. This is due to the different orientation of the substituent in the molecule. The α angle in the enol form is smaller than in the ketone form. The substituent always causes an increase of the α angle in the keto and enol form relatively to the unsubstituted molecule. The highest values of the α angle in keto form has been obtained for structure 4K6a and in enol form for the structure 4Ea. Both structures have a tertbutyl substituent so the increase of the α angle is connected with steric interaction of a bulky substituent.
Table 2 shows the energy differences between the keto and enol form. For each substituent, the lowest energy has been obtained for the structure in the keto form with two hydrogen bonds.
In order to better understand the influence of substituents and hydrogen bonds on the aromaticity of the ring, the HOMA index was determined [42]. For the benzene aromatic ring, the HOMA index is equal to 1; HOMA=0 for the hypothetical structure of 1,3,5-cyclohexatriene with the reference C-C and C=C bonds of buta-1,3-diene [42]; and for the antiaromatic ring, it is negative.
H O M A = 1 α n R o R i 2
  • Ro—the optimized CC bond length of a perfectly aromatic system equal to 1.388 Å.
  • Ri—determined bond length.
  • α—standardization constant of 257.7.
  • n—number of bonds.
The determined HOMA parameters for every ring of the investigated molecules are presented in Table 2. Additionally, as it was done in the previous work [1], to show the dependence of the particular geometry on the electron density, the value of the HOMA parameter of the middle ring has been correlated with the α angle. In the previous work [1], we determined the HOMA parameter for anthrones and anthraquinones with the same substituents in the middle ring. All compounds, irrespective of the nature of the substituent in the middle ring, which have OH groups in the side rings which do not form hydrogen bonds with the adjacent carbonyl group (enones), are characterized by the lowest value of the HOMA parameter for the middle ring. For the compounds investigated in this work the HOMA values for the middle ring are lower, if the OH groups in the aromatic side rings are not present. The presence of one or two hydrogen bonds in ketone form shifts the HOMA value of the middle ring towards higher aromaticity. The highest values of HOMA have been observed for all structures in the enol form with two hydrogen bonds. The electron density of the middle ring is related to the α angle and both parameters can be correlated (Figure 1). In Figure 1, the E and K6 structures have been used according to Scheme 3. As the HOMA shifts towards higher aromaticity, the α angle decreases and hence the molecule flattens out. Similar conclusions were obtained in the previous work [1]. The highest values of the alpha angle for the enol form were obtained for structures with a bulky tert-butyl substituent.
According to Ośmiałowski et al. [34], aromaticity and high π-electron delocalization in the whole system in the enol form plays an important role and determines the tautomeric preference of monohydroxyarenes. Nevertheless, high electron delocalization and stability in individual rings of the condensed systems also affect tautomeric equilibrium in the gas phase. The high HOMA in the side rings of the compounds affects the shift of the tautomeric equilibrium toward the ketone form.

2.4. Electron Density and Ellipticity at Bond Critical Point

The ellipticity of the electron density at the bond critical point (BCP) gives information about the nature of the C-C bond. A correlation between ellipticity and bond length for C8a-C9 and C9a-C9 has been found (Figure 2). The ellipticity and the length of the C-C bond in the middle ring is influenced by the form of the molecule is in and the number of the hydrogen bonds. No effect of the substituent on ellipticity and bond length has been noted. The highest ellipticity is observed for bonds in enol form. Suitable bond length of the C8a-C9 bond is in the range of 1.406–1.419 Å. The length of the C9a-C9 bond is in a slightly larger range of 1.405–1.424 Å. The smallest ellipticity and the highest C8a-C9 and C9a-C9 bond lengths are characteristic for the structures in the ketone form without intramolecular hydrogen bonds. The bond lengths of C8a-C9 and C9a-C9 have also been correlated with electron density at the BCP and the obtained correlation equations are: y = −2.1207x + 2.0518, R² = 0.9986; y = −2.2251x + 2.0799, R2 = 0.9971.

2.5. Selected Bands in the Theoretical IR Spectra of Dihydroanthrones

To study the effect of the substituent on the IR spectra, vibrational frequencies have been calculated for the favored structure in ketone form. In Table 3, the bands characteristic for C=O and OH groups are collected. For the NO2 substituent, the strongest νasOH band is observed at 3399 cm−1. This band repeats at a similar wavenumber for other substituents. Similar wavenumber repetition is observed for the vibration of νCO and δOH and the bands are located close to 1636 and 1670 cm−1. For the structure with the NH2 substituent, the vibration of νCO; δOH; δCH (C3, C4, C5, C6, C10) are observed at 1333 cm−1, which is not observed for other substituents.

2.6. Dipole Moment, Average Local Ionization Energy, and Electrostatic Potential

The dipole moment of the molecule provides important information about its structure. Table 4 shows the dipole moments of enol and ketone form with two hydrogen bonds. The dipole moments are in the range of 0.78–2.84 D. The highest value is obtained for the 9Eb structure. The COOH group significantly increases the polarity of the molecules. However, with the increase of the non-polar carbon chain, the polarity of the molecules decreases which can be observed in the structure with a tertbutyl substituted. Only in some cases (-Cl, CHO, NO2 substituents), the effect of the form of the compound on the dipole moment is observed. For the compounds with other substituents the values of the dipole moment in keto and enol form are very similar.
Average local ionization energy (ALIE) is a function to reveal regions containing highly reactive electrons in chemical system. The surface used for presentation of the local ionization energy is plotted on ρ = 0.0005 a.u. isosurface. Light blue spheres correspond to ALIE minima on the isosurface, revealing favorable sites of electrophilic attack. Darker blue color reveals relatively low ALIE regions. In this region, the electrons have relatively high reactivity.
In the structure of the hydrogen-substituted enol form, electrophilic attack is possible on C10 and on carbon atoms in the side rings, which are not common with the middle ring. In the enol form, the C10 atom is not subject to electrophilic attack. The electron donor substituent in the 2Eb structure directs the electrophilic attack to the side ring atoms that are not in common with the middle ring. In the ketone form 2K6c, the electrophilic attack is directed at the NH2 substituent. The OH substituent in the 3Eb structure revealed a favorable site of electrophilic attack on the carbon atoms of the side rings, which are not shared with the middle ring. In the enol structure, the site of electrophilic attack is the carbon atoms in the side rings that are not shared with the central ring. The site of electrophilic attack becomes less sensitive in the keto structure.
Electrostatic potential (ESP) is popular to visualize the electrostatic nature of molecules. It gives information about the chemical reactivity of a molecule indicating positively and negatively charged fragments of the molecule. For the structures in Table 4, the blue surface corresponds to minimal and red to maximal value of ESP.
For the enol form, the minimum of ESP is located at the oxygen atom not accepting the proton what is connected with the oxygen lone pairs. In the ketone forms, the minimum of ESP is distributed on the oxygen atoms in the side rings OH groups at the and at CO in the middle ring. In addition, the minimum of ESP is observed on COOH, NO2, and Cl substituents and on the oxygen from CHO group. In the enol form, the maximum of ESP is accumulated on one of the hydrogen atoms in the OH group at the side ring.

2.7. Prototropy

Structures with OH and CHO substituents have labile protons, as they move around the molecule, can form a combination of different types of prototropic transformations [43]. Figure 3 shows all possible prototropic forms for OH-substituted 1,8-dihydroxy-9-anthrone. Five prototropic forms are obtained for the OH substituent and one for the CHO substituent (Figure 4). ΔE values, which are calculated relative to the isomer with the lowest energy, are placed under the structures. The highest energy for the OH substituent is obtained for the 3K5-C9a4aH structure, with a single bond between the C9a-C4a atoms. The most close in energy to the 3K6c structure is 3K6-C23H, which has a single bond between C2-C3 atoms. The only one prototropic transformation that has been obtained for the CHO substituent is very close energetically to the 8K6c structure (ΔE = 2.2 kcal/mol).

2.8. Antioxidant Activities

An important parameter for predicting antioxidant properties is the bond dissociation enthalpy (BDE) of the OH bond at C9 in enol form. The BDE of OH has been determined using the calculated total enthalpies (2). Additionally, BDE for C-H at C10 in keto form has been determined (3).
B D E O H = H ( e n o l O · ) H H · + H e n o l O H
H = B D E C H B D E O H

Investigation of the Substituents Effect on BDE (O-H) and BDE (C-H)

Table 5 shows the thermodynamic data for the tautomerization reaction. BDE values are given for C-H at C10 in keto form and O-H bond at C9 in the enol form. Structures with two and one intramolecular hydrogen bond have been considered. The low BDE(O-H) ensures easier detachment of hydrogen from the hydroxyl group. A low BDE(C-H) causes that the hydrogen attached to the C10 is easily removed by radicals [12]. A strongly electron-withdrawing substituent in the central ring like NO2 causes a significant increase of BDE(O-H) and BDE(C-H) in a structure with two intramolecular hydrogen bonds comparing to the unsubstituted structure. In structure with NO2 and one hydrogen bond slight increase of BDE(O-H) has been observed as well as a slight decrease of BDE(C-H) comparing to the unsubstituted structure. This means that the presence of the NO2 substituent and two intramolecular hydrogen bonds makes more difficult detaching the hydrogen atom from the OH and CH group. The electron donating substituent (NH2) causes a slight decrease in the BDE(O-H) value relatively to the unsubstituted molecule in structures with two and one intramolecular hydrogen bond. Thus, it is easier to detach the hydrogen from the OH group at C9 in enol form relatively to the unsubstituted structure. However, according to Korth and Mulder, a low BDE(O-H) alone does not make a compound a good antioxidant [12]. The opposite situation has been observed for BDE(C-H). The NH2 group caused a slight increase in the BDE(C-H) which means that it is more difficult to detach the hydrogen atom relatively to the unsubstituted structure. Similar BDE(O-H) results were obtained by Lucarini and Pedulli [26].

2.9. Hydrogen Bonds

Due to the presence of intramolecular hydrogen bonds, it is important to study their influence on the antioxidant potential and the tautomeric balance. In the conducted analysis, it is noticed that the BDE(O-H) and BDE(C-H) values are lower for the structures with two hydrogen bonds than with one hydrogen bond. Therefore, the presence of two hydrogen bonds increases the antioxidant properties of the compounds.
Hydrogen bonds have a significant impact on the enthalpy difference between ketone and enol structure. The free enthalpy difference between the ketone and enol form is greater for structures with two hydrogen bonds than for structures with one hydrogen bond. This means that two intramolecular hydrogen bonds have a greater effect on shifting the equilibrium towards the ketone form than the presence of one hydrogen bond.
To better illustrate the differences in energy associated with the presence of hydrogen bonds, Scheme 5 shows a cycle of changes in energy for 1,8-dihydroxy-9-anthrone. When, in ketone form, the hydrogens of the OH groups in the side aromatic rings are directed towards the center ring, hydrogen bonds are formed. The energy differences between 1K6a 1K6b, 1K6b 1K6c, and 1Ea 1Eb are the hydrogen bond energies. Higher intramolecular hydrogen bond energy is observed for structures in the ketone form with two intramolecular hydrogen bonds.
The energy of intramolecular hydrogen bonds is influenced by the substituents in the middle ring. For most compounds in the ketone form, the hydrogen bond energy is close to 12 kcal/mol. Form compounds in the enol form it is 8 kcal/mol. The exceptions are the compounds in the ketone form with COOH, CH3, Cl, and CH2CH3 substituents, in which the energies of the hydrogen bond are within the range of 8–19 kcal/mol.

2.10. Influence of the Substituent on the Transition State

To better understand the proton transfer between the ketone and enol form, the most probably transition state has been optimized. According to Rubén Marrero-Carballo et al. [15], proton transfer in the ketone form occurs from the methyl group in the middle ring to the oxygen in the carbonyl group to form the enol form. In the current work, the same proton pathway has been used but the effect of three different substituents in the middle ring on the transition state geometry and energy has additionally been investigated. Substituents that accept electrons from the middle ring or donate them to the ring have been chosen.
All substituents have subtle effects on the transition state energies (Table 6). The differences in the transition state energies of the hydrogen atom-substituted structure and the NO2 group-substituted structure are very similar. Despite the difference in the nature of the selected substituents, no direct interactions between the traveling proton and the substituents are observed.
For every transition state, the hydrogen bonds between the OH groups at the side rings and the CO group in the middle ring are broken. In the case of NH2 and OH substituents, the hydrogens from the side OH groups are directed toward each other (Figure 5). Only in the case of the NO2 substituent, both hydrogens from the OH groups face one direction. This is most likely due to the mesomeric effect of the NO2 group. In addition, the angle of the CO bond to the plane of the central ring constructed with four carbon atoms in common with the side aromatic rings has been investigated. The smallest angle is for the structure with the NO2 substituent, which pulls electrons out of the ring. On the other hand, the highest angle has been obtained when the substituent is the NH2 group, which is an electron donor.

2.11. Transition State for the Keto-Enol Reaction with Pyridine

Another proton transfer pathway has been carried out with the help of a pyridine molecule as a carrier [15] for the keto-enol reaction of 1,8-dihydroxy-9-anthrone substituted with NO2, OH, and H (Table 7). The lowest energy has been obtained for the initial state-pyridine with a substituted anthralin structure in the ketone form K + Pyr. In the first step, the proton from C10 is transferred to the nitrogen atom in the pyridine molecule. ΔG for the NO2 substituent is 17.60 kcal/mol and for the OH substituent ΔG is equal to 17.34 kcal/mol. The lowest transition state energy TS1 of 8.37 kcal/mol has been obtained for the hydrogen atom substituted structure. The next step on the reaction progress is the formation of the INT1 ion. The nitrogen in the pyridine moves toward the oxygen in the middle ring of the substituted anthralin to form a hydrogen bond. This results in the formation of an ion pair (INT2). In order for the protonated pyridine to release a proton into oxygen yielding the enol form E + Pyr, a low free energy is needed. In the case of the NO2 substituent, the ΔG is 13.93 kcal/mol, in the case of the OH group it is only 3.74 kcal/mol and for the H substituent-10.18 kcal/mol. For a better illustration of the reaction pathway involving pyridine, are provided figures showing each step of the reaction (Figure 6).

2.12. Theoretical Reaction Rate Constants

The rate constants in Table 8 illustrate how slow the conversion of the keto form to the enol form in 1,8-dihydroxy-9-anthrone molecules is. Comparison of the rate constants for the keto-enol reaction in substituted 1,8-dihydroxy-9-anthrones indicates that this reaction is not sensitive to substitution, except for the substitution of the NH2 group, which accelerates the investigated reaction.
The addition of pyridine changes the reaction mechanism, making it a two-stage process, with the second stage clearly faster. In the case of both stages of the reaction, the effect of the substituent on the reaction rate is visible (Table 9), which is particularly evident in the case of the second stage of the reaction of the OH-substituted 1,8-dihydroxy-9-anthrone.

3. Conclusions

The lowest energy structure is the K6 tautomer, which represents the ketone form.
The intramolecular hydrogen bond affects the geometry of 1,8-dihydroxy-9-anthrone. The changes of the α angle, which is the angle between the planes of the side rings, and so representing the general shape of 1,8-dihydroxy-9-anthrone, is sensitive to the OHO hydrogen bonds.
The electronodonor and electronoacceptor properties of the substituent in the middle ring as well as the presence of intramolecular hydrogen bonds affects the antioxidant properties of 1,8-dihydroxy-9-anthrone. The greatest decrease in BDE(O-H) has been obtained for the structure with the NH2 substituent in the middle ring. Therefore, this compound can be expected to have the highest antioxidant properties.
The energy of the transition state in the keto-enol reaction in 1,8-dihydroxy-9-anthrone is unsensitive to the properties of the substituents in the central ring.
Taking into account the ΔH and the HOMA of the side aromatic rings, the preferred structure is the ketone form regardless the nature of the substituent in the middle ring of 1,8-dihydroxy-9-anthrone.
The keto-enol transformation in 1,8-dihydroxy-9-anthrone is significantly faster in the presence of pyridine.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28010344/s1, Table S1: Bond lengths in crystal structures of 1,8-dihydroxy-9-anthrone derivatives.

Author Contributions

All authors (M.S. and I.M.) contributed to the conceptualization, methodology and writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Xyz or wfn files are available on request from the corresponding author.

Acknowledgments

The Wroclaw Center for Networking and Supercomputing is acknowledged for generous allocations of computer time. This work is supported by Wroclaw Medical University grant No SUBZ.D050.22.025.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Szymańska, M.; Majerz, I. Effect of Substitution of Hydrogen Atoms in the Molecules of Anthrone and Anthraquinone. Molecules 2021, 26, 502. [Google Scholar] [CrossRef]
  2. Marsden, J.; Coburn, P.; Marks, J.; Shuster, S. Measurement of the response of psoriasis to short-term application of anthralin. Br. J. Dermatol. 1983, 109, 209–218. [Google Scholar] [CrossRef]
  3. Ashton, R.E.; Andre, P.; Lowe, N.J.; Whitefield, M. Anthralin: Historical and current perspectives. J. Am. Acad. Dermatol. 1983, 9, 173–192. [Google Scholar] [CrossRef]
  4. Hellier, F.F.; Whitefield, M. The treatment of psoriasis with triacetoxyanthracene. Br. J. Dermatol. 1967, 79, 491–496. [Google Scholar] [CrossRef]
  5. Müller, K. Antipsoriatic and proinflammatory action of anthralin: Implications for the role of oxygen radicals. Biochem. Pharmacol. 1997, 53, 1215–1221. [Google Scholar] [CrossRef]
  6. Pečar, S.; Schara, M.; Müller, K.; Wiegrebe, W. Reduction of nitroxides by anthralin and some of its derivatives. Free Radic. Biol. Med. 1995, 18, 459–465. [Google Scholar] [CrossRef]
  7. Mueller, K.; Guerster, D.; Piwek, S.; Wiegrebe, W. Antipsoriatic anthrones with modulated redox properties. 1. Novel 10-substituted 1,8-dihydroxy-9(10H)-anthracenones as inhibitors of 5-lipoxygenase. J. Med. Chem. 1993, 36, 4099–4107. [Google Scholar] [CrossRef] [Green Version]
  8. Szymanski, S.; Majerz, I. Aromaticity and Electron Density of Hypericin. J. Nat. Prod. 2019, 82, 2106–2115. [Google Scholar] [CrossRef]
  9. Szymanski, S.; Majerz, I. In Silico Studies on Sennidines—Natural Dianthrones from Senna. Biology 2021, 10, 468. [Google Scholar] [CrossRef]
  10. Szymański, S.; Majerz, I. Theoretical Studies on the Structure and Intramolecular Interactions of Fagopyrins—Natural Photosensitizers of Fagopyrum. Molecules 2022, 27, 3689. [Google Scholar] [CrossRef]
  11. Limacher, P.A.; Lüthi, H.P. Cross-conjugation. WIREs Comput. Mol. Sci. 2011, 1, 477–486. [Google Scholar] [CrossRef]
  12. Korth, H.-G.; Mulder, P. Anthrone and Related Hydroxyarenes: Tautomerization and Hydrogen Bonding. J. Org. Chem. 2013, 78, 7674–7682. [Google Scholar] [CrossRef]
  13. Takemura, T.; Baba, H. Spectrophotometric investigations of the tautomeric reaction between anthrone and anthranol—II. Tetrahedron 1968, 24, 5311–5321. [Google Scholar] [CrossRef]
  14. Fain, V.Y.; Zaitsev, B.E.; Ryabov, M.A. Tautomerism of anthraquinones: IV. 1-Hydroxy-9,10-anthraquinone and its substituted derivatives. Russ. J. Org. Chem. 2006, 42, 1469–1472. [Google Scholar] [CrossRef]
  15. Marrero-Carballo, R.; Tun-Rosado, F.J.; Mena-Rejón, G.J.; Cáceres, D.; Barroso, J.; Murillo, F.; Merino, G.; Quijano-Quiñones, R.F. The base-catalyzed keto-enol tautomerism of chrysophanol anthrone. A DFT investigation of the base-catalyzed reaction. Mol. Simul. 2019, 45, 716–723. [Google Scholar] [CrossRef]
  16. Radi, S.; Tighadouini, S.; Feron, O.; Riant, O.; Bouakka, M.; Benabbes, R.; Mabkhot, Y.N. Synthesis of Novel β-Keto-Enol Derivatives Tethered Pyrazole, Pyridine and Furan as New Potential Antifungal and Anti-Breast Cancer Agents. Molecules 2015, 20, 20186–20194. [Google Scholar] [CrossRef] [Green Version]
  17. Gad, S.F.; El-Demerdash, S.H.; El-Mehasseb, I.M.; El-Nahas, A.M. Structure, stability and conversions of tautomers and rotamers of azulene-based uracil analogue. J. Mol. Struct. 2019, 1182, 271–282. [Google Scholar] [CrossRef]
  18. Brovarets’, O.O.; Hovorun, D.M. Prototropic tautomerism and basic molecular principles of hypoxanthine mutagenicity: An exhaustive quantum-chemical analysis. J. Biomol. Struct. Dyn. 2013, 31, 913–936. [Google Scholar] [CrossRef]
  19. Samijlenko, S.P.; Yurenko, Y.P.; Stepanyugin, A.V.; Hovorun, D.M. Tautomeric Equilibrium of Uracil and Thymine in Model Protein−Nucleic Acid Contacts. Spectroscopic and Quantum Chemical Approach. J. Phys. Chem. B 2010, 114, 1454–1461. [Google Scholar] [CrossRef]
  20. Ziółek, M.; Kubicki, J.; Maciejewski, A.; Naskrecki, R.; Grabowska, A. Enol-keto tautomerism of aromatic photochromic Schiff base N,N′-bis(salicylidene)-p-phenylenediamine: Ground state equilibrium and excited state deactivation studied by solvatochromic measurements on ultrafast time scale. J. Chem. Phys. 2006, 124, 124518. [Google Scholar] [CrossRef]
  21. Laurella, S.L.; Sierra, M.G.; Furlong, J.J.P.; Allegretti, P.E. Substituent, Temperature and Solvent Effects on the Keto-Enol EQUILIBRIUM in β-Ketoamides: A Nuclear Magnetic Resonance Study. Open J. Phys. Chem. 2013, 3, 138–149. [Google Scholar] [CrossRef] [Green Version]
  22. Wu, C.-C.; Lien, M.-H. Ab Initio Study on the Substituent Effect in the Transition State of Keto−Enol Tautomerism of Acetyl Derivatives. J. Phys. Chem. 1996, 100, 594–600. [Google Scholar] [CrossRef]
  23. Antonov, L.; Fabian, W.M.F.; Nedeltcheva, D.; Kamounah, F.S. Tautomerism of 2-hydroxynaphthaldehyde Schiff bases. J. Chem. Soc. Perkin Trans. 2000, 2, 1173–1179. [Google Scholar] [CrossRef]
  24. Taylor, P.J.; van der Zwan, G.; Antonov, L. Tautomerism: Introduction, History, and Recent Developments in Experimental and Theoretical Methods. Tautomerism 2013, 2013, 1–24. [Google Scholar] [CrossRef]
  25. Nazarparvar, E.; Zahedi, M.; Klein, E. Density Functional Theory (B3LYP) Study of Substituent Effects on O–H Bond Dissociation Enthalpies of trans-Resveratrol Derivatives and the Role of Intramolecular Hydrogen Bonds. J. Org. Chem. 2012, 77, 10093–10104. [Google Scholar] [CrossRef]
  26. Lucarini, M.; Pedulli, G.F. Free radical intermediates in the inhibition of the autoxidation reaction. Chem. Soc. Rev. 2010, 39, 2106–2119. [Google Scholar] [CrossRef]
  27. Szymańska, M.; Majerz, I. Geometry and electron density of phenothazines. J. Mol. Struct. 2020, 1200, 127095. [Google Scholar] [CrossRef]
  28. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian Inc. 16, Revision A.03; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  29. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef] [Green Version]
  30. Lee, C.; Yang, W.; 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] [Green Version]
  31. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104–154119. [Google Scholar] [CrossRef]
  32. Keith, T.A. AIMALL (Version 19.10.12); TK Gristmill Software: Overland Park, KS, USA, 2019. [Google Scholar]
  33. Allen, F.H. The Cambridge Structural Database: A quarter of a million crystal structures and rising. Acta Crystallogr. Sect. B Struct. Sci. 2002, 58, 380–388. [Google Scholar] [CrossRef]
  34. Osmialowski, B.; Raczyńska, E.D.; Krygowski, T.M. Tautomeric Equilibria and Pi Electron Delocalization for Some MonohydroxyarenesQuantum Chemical Studies. J. Org. Chem. 2006, 71, 3727–3736. [Google Scholar] [CrossRef]
  35. Brown, C.J.; Colclough, M.L. 1,8-Dinitro-4,5-dihydroxyanthraquinone, C14H6N2O8. Acta Crystallogr. Sect. C Cryst. Struct. Commun. 1983, 39, 300–302. [Google Scholar] [CrossRef] [Green Version]
  36. Von Dreele, R.B.; Einck, J.J. The crystal and molecular structure of carminomycin I hydrochloride monohydrate. Acta Crystallogr. Sect. B Struct. Crystallogr. Cryst. Chem. 1977, 33, 3283–3288. [Google Scholar] [CrossRef]
  37. Rohl, A.; Moret, M.; Kaminsky, W.; Claborn, K.; McKinnon, J.; Kahr, B. Hirshfeld Surfaces Identify Inadequacies in Computations of Intermolecular Interactions in Crystals: Pentamorphic 1,8-Dihydroxyanthraquinone. Cryst. Growth Des. 2008, 8, 4517–4525. [Google Scholar] [CrossRef]
  38. Schmidt-Bäse, K.; Noltemeyer, M.; Egert, E.; Eigelt, E.; Zeeck, A. Structure of the anthracycline antibiotic aranciamycin. Acta Crystallogr. Sect. C Cryst. Struct. Commun. 1993, 49, 250–253. [Google Scholar] [CrossRef]
  39. Armaghan, M.; Amini, M.M.; Ng, S.W.; Tiekink, E.R.T. Co-crystals of 1,8-dihydroxy-2,4,5,7-tetranitro-9,10-anthraquinone with dibenzothiophene and 4,6-dimethyldibenzothiophene. Z. Kristallographie-Crystalline Mater. 2013, 228, 598–606. [Google Scholar] [CrossRef]
  40. Hernandez-Medel, M.D.R.; O Ramirez-Corzas, C.; Rivera-Dominguez, M.; Ramirez-Mendez, J.; Santillan, R.; Rojas-Lima, S. Diastereomeric C-glycosyloxanthrones from picramnia antidesma. Phytochemistry 1999, 50, 1379–1383. [Google Scholar] [CrossRef]
  41. Aoyama, T.; Naganawa, H.; Muraoka, Y.; Nakamura, H.; Aoyagi, T.; Takeuchi, T.; Iitaka, Y. Benastatins A and B, new inhibitors of glutathione S-transferase, produced by Streptomyces sp. MI384-DF12. II. Structure determination of benastatins A and B. J. Antibiot. 1992, 45, 1391–1396. [Google Scholar] [CrossRef] [Green Version]
  42. Krygowski, T.M. Crystallographic studies of inter- and intramolecular interactions reflected in aromatic character of .pi.-electron systems. J. Chem. Inf. Comput. Sci. 1993, 33, 70–78. [Google Scholar] [CrossRef]
  43. Raczyńska, E.D.; Kosińska, W.; Ośmiałowski, B.; Gawinecki, R. Tautomeric Equilibria in Relation to Pi-Electron Delocalization. Chem. Rev. 2005, 105, 3561–3612. [Google Scholar] [CrossRef]
Scheme 1. Molecular structure and atom numbering of 1,8-dihydroxy-9-anthrones.
Scheme 1. Molecular structure and atom numbering of 1,8-dihydroxy-9-anthrones.
Molecules 28 00344 sch001
Scheme 2. The keto-enol equilibrium in 1,8-dihydroxy-9-anthrone.
Scheme 2. The keto-enol equilibrium in 1,8-dihydroxy-9-anthrone.
Molecules 28 00344 sch002
Scheme 3. Tautomeric equilibrium in 1,8-dihydroxy-9-anthrone.
Scheme 3. Tautomeric equilibrium in 1,8-dihydroxy-9-anthrone.
Molecules 28 00344 sch003
Scheme 4. Definition of the α angle for 1,8-dihydroxy-9-anthrones as the angle between the side rings planes.
Scheme 4. Definition of the α angle for 1,8-dihydroxy-9-anthrones as the angle between the side rings planes.
Molecules 28 00344 sch004
Figure 1. Correlations between the HOMA parameter for the central ring and α angle: yellow—enols, blue—ketones.
Figure 1. Correlations between the HOMA parameter for the central ring and α angle: yellow—enols, blue—ketones.
Molecules 28 00344 g001
Figure 2. Correlations between ellipticity and bond length for (a)—C8a-C9 bond, (b)—C9a-C9 bond. Region A—enol, B—ketone form with two hydrogen bonds, C—ketone form with one hydrogen bond, D—ketone form without hydrogen bonds.
Figure 2. Correlations between ellipticity and bond length for (a)—C8a-C9 bond, (b)—C9a-C9 bond. Region A—enol, B—ketone form with two hydrogen bonds, C—ketone form with one hydrogen bond, D—ketone form without hydrogen bonds.
Molecules 28 00344 g002
Figure 3. Possible prototropic forms for OH-substituted anthrone.
Figure 3. Possible prototropic forms for OH-substituted anthrone.
Molecules 28 00344 g003
Figure 4. Prototropic tautomers and ΔE at 298 K: (a)—OH substituent, (b)—CHO substituent.
Figure 4. Prototropic tautomers and ΔE at 298 K: (a)—OH substituent, (b)—CHO substituent.
Molecules 28 00344 g004
Scheme 5. Energy differences [kcal/mol] between the structures of the unsubstituted 1,8-dihydroxy-9-anthrone.
Scheme 5. Energy differences [kcal/mol] between the structures of the unsubstituted 1,8-dihydroxy-9-anthrone.
Molecules 28 00344 sch005
Figure 5. Transition state for 1,8-dihydroxy-9-anthrone without substituent—(a), with OH substituent—(b), with NH2—(c), and with NO2—(d).
Figure 5. Transition state for 1,8-dihydroxy-9-anthrone without substituent—(a), with OH substituent—(b), with NH2—(c), and with NO2—(d).
Molecules 28 00344 g005
Figure 6. Steps of proton transfer reaction involving pyridine for 1,8-dihydroxy-9-anthrone with NO2 substituent, (a)—Reaction progress ΔG [kcal/mol] at 298 K, (b)—optimized structures.
Figure 6. Steps of proton transfer reaction involving pyridine for 1,8-dihydroxy-9-anthrone with NO2 substituent, (a)—Reaction progress ΔG [kcal/mol] at 298 K, (b)—optimized structures.
Molecules 28 00344 g006
Table 1. G at 298 K, α and HOMA parameter for optimized isomers.
Table 1. G at 298 K, α and HOMA parameter for optimized isomers.
XIsomerΔG [kcal/mol]αHOMA (A)HOMA (B)HOMA (C)
HE13.80.230.61880.66660.5957
K19.10.030.79270.6583−0.7815
K210.70.010.77690.7006−0.8107
K30.06.730.9053−0.71590.9053
K412.40.010.85500.4044−0.5983
K510.20.010.86020.3176−0.5043
K60.06.720.9053−0.71590.9053
NH2E10.33.510.61900.61780.5807
K14.93.040.79940.5773−0.7667
K25.62.790.77660.6521−0.7804
K30.012.010.9095−0.91300.9098
K47.33.040.84800.3583−0.5692
K56.14.290.86000.2443−0.5002
K60.012.020.9019−1.02180.9019
NO2E16.80.880.62230.65180.6053
K113.00.730.79670.6358−0.7255
K212.92.530.77600.6861−0.7420
K30.014.930.9176−0.62550.9176
K414.03.210.84460.4079−0.4888
K513.80.800.85560.3158−0.4320
K60.014.970.9176−0.62550.9176
OHE10.82.710.61590.68070.6067
K16.41.950.81100.6111−0.7856
K27.60.890.79450.6725−0.8222
K30.013.810.9074−0.96320.9074
K49.30.810.86600.3664−0.6312
K57.51.530.87230.2799−0.5133
K60.013.830.9074−0.96320.9074
Table 2. ΔG at 298 K, α and HOMA parameter for optimized structures of enol and ketone forms. Angle defined according to Scheme 1. For each substituent, the lowest energy has been obtained for the structure in the keto form with two hydrogen bonds.
Table 2. ΔG at 298 K, α and HOMA parameter for optimized structures of enol and ketone forms. Angle defined according to Scheme 1. For each substituent, the lowest energy has been obtained for the structure in the keto form with two hydrogen bonds.
XIsomerΔG [kcal/mol]αHOMA (A)HOMO (B)HOMA (C)
H1Ea21.80.450.61290.63070.6091
1Eb13.80.230.61880.66660.5957
1K6a24.628.080.9567−1.32190.9567
1K6b12.113.700.9300−0.93600.9181
1K6c0.06.720.9053−0.71590.9053
NH22Ea18.83.740.60880.58880.5984
2Eb10.33.510.61900.61780.5807
2K6a27.240.590.9629−1.65490.9629
2K6b16.231.660.9385−1.27960.9229
2K6c0.012.020.9019−1.02180.9019
OH3Ea19.74.280.60360.62650.6220
3Eb10.82.710.61590.66220.6252
3K6a23.628.700.9583−1.56920.9583
3K6b11.517.620.9329−1.20030.9196
3K6c0.013.830.9074−0.96320.9074
t−Bu4Ea24.522.570.57240.51000.5744
4Eb17.221.170.56490.53010.5467
4K6a21.052.000.9593−2.14100.9595
4K6b10.946.080.9314−1.76950.9175
4K6c0.041.590.8914−1.50240.8914
Et5Ea23.73.720.58490.56610.5549
5Eb15.43.380.58540.60030.5647
5K6a27.926.530.9549−1.69220.9549
5K6b18.836.690.9332−1.44960.9248
5K6c0.020.800.9083−0.80270.9098
Me6Ea23.20.960.60670.56820.5523
6Eb15.10.510.60660.59180.5384
6K6a27.538.320.9599−1.69600.9599
6K6b16.128.470.9305−1.34030.9194
6K6c0.015.000.9097−0.80680.9097
Cl7Ea22.80.680.59900.58540.5947
7Eb14.90.430.60770.62580.5854
7K6a27.030.390.9585−1.53860.9585
7K6b11.320.470.9409−0.77150.9275
7K6c0.017.120.9156−0.52780.9156
CHO8Ea21.96.800.63960.50810.6384
8Eb15.04.970.62240.53790.6198
8K6a23.427.590.9567−1.36380.9567
8K6b11.416.990.9324−1.00070.9207
8K6c0.012.140.9088−0.78170.9088
COOH9Ea19.01.870.62110.60450.6146
9Eb11.81.450.61610.63610.6585
9K6a25.830.290.9600−1.64220.9621
9K6b7.915.430.9315−1.07560.9169
9K6c0.017.850.9209−0.82430.9209
NO210Ea23.80.740.63040.61630.6264
10Eb16.80.880.62230.65180.6053
10K6a23.326.590.9620−1.21320.9699
10K6b11.318.680.9451−0.85500.9309
10K6c0.014.970.9176−0.62550.9176
Table 3. Selected bands without a scaling factor for the optimized dihydroanthrones. Only the favored keto form (K6) with different substituents has been included.
Table 3. Selected bands without a scaling factor for the optimized dihydroanthrones. Only the favored keto form (K6) with different substituents has been included.
SubstituentWavenumber [cm−1]Int. [km/mol]Description
NO2 1381109.2νCO; δOH; δCH; νCN; νNO; νCC
1399189.7νCO; δOH; δCH; νCN; νNO; νCC
1519230.8νCO; δOH; δCH (side rings); νCC
1636144.5νCO; δOH; δCH; νCC
1670261.2νCO; δOH; δCH; νCC
3399520.1νasOH
342750.2νsOH
t-Bu1502213.1νCO; δOH; δCH; νCC (rings)
1634286.5νCO; δOH; νCC (rings); δCH (rings)
1662193.3νCO; δOH; νCC (rings); δCH (rings)
3386436.3νasOH
341567.5νsOH
COOH1394144.3νCO; δOH; δCH; νCC
1518213.9νCO; δOH; δCH; νCC
1635186.8νCO; δOH; δCH; νCC
1669236.9νCO; δOH; δCH; νCC
3389561.5νasOH
341851.6νsOH
CH3138996.2νCO; δOH; δCH; νCC (rings)
1515236.5νCO; δOH; δCH (rings); νCC (rings)
1633254.5νCO; δOH; δCH (rings); νCC (rings)
1668192.3νCO; δOH; δCH (rings); νCC (rings)
3372601.1νasOH
340345.1νsOH
CHO1514209.6νCO; δOH; δCH (side rings); νCC
1629279.4νCO; δOH; δCH (side rings); νCC
1667224.9νCO; δOH; δCH (side rings); νCC
3381574.2νasOH
341146.4νsOH
Cl1395118.7νCO; δOH; δCH (C2, C7); νCC
1517237.3νCO; δOH; δCH (side rings); νCC
1634171.1νCO; δOH; δCH (C3, C4, C5, C6); νCC
1669231.5νCO; δOH; δCH (side rings); νCC
3381504.5νasOH
341749.1νsOH
Et1391109.6νCO; δOH; δCH (C2, C7); νCC (rings)
1514223.9νCO; δOH; δCH; νCO; νCC (rings)
1633264.5νCO; δOH; δCH (C3, C4, C5, C6); νCC (rings)
1666193.6νCO; δOH; δCH (C2, C4, C5, C7); νCC (rings)
3377559.1νasOH
340751.4νsOH
OH1391129.5νCO; δOH; δCH (C2, C7, C10); νCC (rings)
1511269.6νCO; δOH; δCH (rings); νCC
1633205.3νCO; δOH (side rings); δCH (C3, C4, C5, C6); νCC
1668200.7νCO; δOH (side rings); δCH (C2, C4, C5, C7); νCC
3377586.9νasOH
340745.7νsOH
NH21333101.7νCO; δOH; δCH (C3, C4, C5, C6, C10); δNH; νCC
1510244.9νCO; δOH; δCH; δNH; νCC; νCN
1635252.2νCO; δOH; δCH; δNH; νCC
1666207.2νCO; δOH; δCH; δNH; νCC
3372 νasOH
3402 νsOH
H1387106.5νCO; δOH; δCH (C2, C7, C10); νCC
1518244.9νCO; δOH; δCH; νCC
1633242.2νCO; δOH; δCH; νCC
1668188.4νCO; δOH; δCH (C2, C4, C5, C7); νCC
3690657.3νasOH
340140.7νsOH
Table 4. Dipol Moment, Average Local Ionization Energy, Electrostatic Potential of enol and ketone form with two hydrogen bonds.
Table 4. Dipol Moment, Average Local Ionization Energy, Electrostatic Potential of enol and ketone form with two hydrogen bonds.
XIsomerDipol Moment [D]ALIEESP
H1Eb1.79Molecules 28 00344 i001Molecules 28 00344 i002
1K6c1.76Molecules 28 00344 i003Molecules 28 00344 i004
NH22Eb1.80Molecules 28 00344 i005Molecules 28 00344 i006
2K6c1.74Molecules 28 00344 i007Molecules 28 00344 i008
OH3Eb1.33Molecules 28 00344 i009Molecules 28 00344 i010
3K6c1.13Molecules 28 00344 i011Molecules 28 00344 i012
t-Bu3Eb1.33Molecules 28 00344 i013Molecules 28 00344 i014
3K6c1.13Molecules 28 00344 i015Molecules 28 00344 i016
Et4Eb1.73Molecules 28 00344 i017Molecules 28 00344 i018
4K6c1.76Molecules 28 00344 i019Molecules 28 00344 i020
Me6Eb1.76Molecules 28 00344 i021Molecules 28 00344 i022
6K6c1.76Molecules 28 00344 i023Molecules 28 00344 i024
Cl7Eb1.95Molecules 28 00344 i025Molecules 28 00344 i026
7K6c1.41Molecules 28 00344 i027Molecules 28 00344 i028
CHO8Eb2.90Molecules 28 00344 i029Molecules 28 00344 i030
8K6c0.78Molecules 28 00344 i031Molecules 28 00344 i032
COOH9Eb2.84Molecules 28 00344 i033Molecules 28 00344 i034
9K6c2.41Molecules 28 00344 i035Molecules 28 00344 i036
NO210Eb2.66Molecules 28 00344 i037Molecules 28 00344 i038
10K6c1.54Molecules 28 00344 i039Molecules 28 00344 i040
Table 5. Enthalpy (ΔH) and Free Energy (ΔG) at 298 K for Tautomerization and the Bond Dissociation Enthalpy, BDE, in the enones (C-H) and in the enols (O-H). Structures without substituents are marked gray.
Table 5. Enthalpy (ΔH) and Free Energy (ΔG) at 298 K for Tautomerization and the Bond Dissociation Enthalpy, BDE, in the enones (C-H) and in the enols (O-H). Structures without substituents are marked gray.
ΔH [kcal/mol]ΔG [kcal/mol]TΔSBDE(C-H)BDE(O-H)
10 K 6 c 10Eb−17.3−16.8−7.82 × 10479.096.3
K 6 c 4Eb−18.2−17.2−1.61 × 10359.878.0
9 K 6 c 9Eb−12.0−11.8−2.94 × 10469.281.1
6 K 6 c 6Eb−17.1−15.1−3.18 × 10361.878.9
8 K 6 c 8Eb−15.4−15.0−6.81 × 10465.580.9
7 K 6 c 7Eb−15.9−14.9−1.58 × 10364.079.9
5 K 6 c 5Eb−16.7−15.4−2.15 × 10362.679.4
3 K 6 c 3Eb−11.9−10.8−1.73 × 10365.677.5
2 K 6 c 2Eb−7.0−5.6−2.11 × 10368.775.6
1 K 6 c 1Eb−14.6−13.8−1.27 × 10365.680.1
10 K 6 b 10Ea−12.9−12.5−5.57 × 10476.088.9
4 K 6 b 4Ea−14.8−13.6−1.89 × 10369.784.5
9 K 6 b 9Ea−11.5−11.2−5.28 × 10476.788.2
6 K 6 b 6Ea−9.3−7.1−3.51 × 10376.485.7
8 K 6 b 8Ea−10.6−10.4−2.90 × 10478.489.0
7 K 6 b 7Ea−12.1−11.4−1.04 × 10374.886.8
5 K 6 b 5Ea−6.1−4.9−1.83 × 10380.286.3
3 K 6 b 3Ea−8.7−8.1−1.01 × 10375.784.4
2 K 6 b 2Ea−4.1−2.7−2.25 × 10378.582.6
1 K 6 b 1Ea−10.1−9.8−5.64 × 10476.987.0
Table 6. Energies at 298 K for optimized structures of the transition state, the reactants, and products for 1,8-dihydroxy-9-anthrone with different substituents.
Table 6. Energies at 298 K for optimized structures of the transition state, the reactants, and products for 1,8-dihydroxy-9-anthrone with different substituents.
ΔE [kcal/mol]
SubstituentReactantTransition StateProduct
H0111.79314.828
OH0105.88112.503
NH2095.3537.501
NO20111.64817.800
Table 7. Relative enthalpies, entropies, and Gibbs energies (kcal/mol) at 298 K for the stationary points involved in the pirydyne catalyzed tautomeric reaction of antralin with two different substituent (in gas phase).
Table 7. Relative enthalpies, entropies, and Gibbs energies (kcal/mol) at 298 K for the stationary points involved in the pirydyne catalyzed tautomeric reaction of antralin with two different substituent (in gas phase).
SubstituentCompoundΔG [kcal/mol]ΔH [kcal/mol]TΔS
NO2K + Pyr0.000.000.00
TS117.6015.75−1.85
INT18.835.42−3.41
INT214.7712.91−1.86
TS223.5219.98−3.54
E + Pyr9.5911.081.49
OHK + Pyr0.000.000.00
TS117.3415.48−1.86
INT114.3811.73−2.65
INT25.595.02−0.57
TS26.936.07−0.86
E + Pyr3.195.131.94
HK + Pyr0.000.000.00
TS18.377.89−0.48
INT119.3417.04−2.30
INT210.8610.68−0.18
TS217.6017.600.00
E + Pyr7.4210.082.66
Table 8. Thermodynamic parameters for the keto-enol reaction of 1,8-dihydroxy-9-anthrone.
Table 8. Thermodynamic parameters for the keto-enol reaction of 1,8-dihydroxy-9-anthrone.
SubstituentΔrH° (298 K) [kcal/mol]ΔrG° (298 K) [kcal/mol]k (298 K)
H−14.58−13.783.13 × 10−66
NH2−6.96−5.644.85 × 10−54
NO2−17.32−16.838.49 × 10−66
OH−11.92−10.831.21 × 10−61
Table 9. Thermodynamic parameters for the keto-enol reaction of 1,8-dihydroxy-9-anthrone with participation of pyridine.
Table 9. Thermodynamic parameters for the keto-enol reaction of 1,8-dihydroxy-9-anthrone with participation of pyridine.
SubstituentΔrH1 (298 K)
[kcal/mol]
ΔrG1 (298 K)
[kcal/mol]
k1 (298)ΔrH2 (298 K)
[kcal/mol]
ΔrG2 (298 K)
[kcal/mol]
k2 (298 K)
H−17.04−19.341.42 × 10−20.593.447.09 × 107
NO2−5.43−8.837.80 × 10−11.835.192.38 × 106
OH−11.73−14.381.21 × 100−0.112.413.66 × 1011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Szymańska, M.; Majerz, I. Prototropy, Intramolecular Interactions, Electron Delocalization, and Physicochemical Properties of 1,8-dihydroxy-9-anthrone—DFT-D3 Study of Substituent Effects. Molecules 2023, 28, 344. https://doi.org/10.3390/molecules28010344

AMA Style

Szymańska M, Majerz I. Prototropy, Intramolecular Interactions, Electron Delocalization, and Physicochemical Properties of 1,8-dihydroxy-9-anthrone—DFT-D3 Study of Substituent Effects. Molecules. 2023; 28(1):344. https://doi.org/10.3390/molecules28010344

Chicago/Turabian Style

Szymańska, Małgorzata, and Irena Majerz. 2023. "Prototropy, Intramolecular Interactions, Electron Delocalization, and Physicochemical Properties of 1,8-dihydroxy-9-anthrone—DFT-D3 Study of Substituent Effects" Molecules 28, no. 1: 344. https://doi.org/10.3390/molecules28010344

APA Style

Szymańska, M., & Majerz, I. (2023). Prototropy, Intramolecular Interactions, Electron Delocalization, and Physicochemical Properties of 1,8-dihydroxy-9-anthrone—DFT-D3 Study of Substituent Effects. Molecules, 28(1), 344. https://doi.org/10.3390/molecules28010344

Article Metrics

Back to TopTop