Theoretical Studies on the Structure and Intramolecular Interactions of Fagopyrins—Natural Photosensitizers of Fagopyrum

Simple Summary The study determines the spatial structure and intramolecular interactions of fagopyrins—natural photosensitizers of Fagopyrum species. In silico calculations show many fagopyrin conformers characterized by the formation of strong intramolecular interactions. Abstract Compounds characterized by a double-anthrone moiety are found in many plant species. One of them are fagopyrins—naturally occurring photosensitizers of Fagopyrum. The photosensitizing properties of fagopyrins are related to the selective absorption of light, which is a direct result of their spatial and electronic structure and many intramolecular interactions. The nature of the interactions varies in different parts of the molecule. The aim of this study is to determine the structure and intramolecular interactions of fagopyrin molecules. For this purpose, in silico calculations were used to perform geometry optimization in the gas phase. QTAIM and NCI analysis suggest the formation of the possible conformers in the fagopyrin molecules. The presence of a strong OHO hydrogen bond was shown in the anthrone moiety of fagopyrin. The minimum energy difference for selected conformers of fagopyrins was 1.1 kcal∙mol−1, which suggested that the fagopyrin structure may exist in a different conformation in plant material. Similar interactions were observed in previously studied structures of hypericin and sennidin; however, only fagopyrin showed the possibility of brake the strong OHO hydrogen bond in favor of forming a new OHN hydrogen bond.


Introduction
Fagopyrins are a group of compounds of natural origin found in plants of the genus Fagopyrum. Parts of these plants are commonly consumed by humans and animals throughout the world [1]. There are many species of Fagopyrum; however, the most consumed and studied are Fagopyrum esculentum, F. tataricum, and F. cymosum [2,3]. Parts of these plants provide a low-calorie, gluten-free food and a source of many elements and organic compounds of biological interest, such as rutin, quercetin, and fagopyrins [4][5][6][7][8].
Fagopyrins are anthraquinone derivatives characterized by a polycyclic system, which is interesting from the chemical point of view. The structure and intramolecular interactions of polycyclic compounds affect the physical and chemical properties and can influence potential applications of the compounds in pharmacy and medicine [9][10][11][12]. As shown in earlier work on hypericin [13] and sennidines [14], a highly substituted polycyclic system can exhibit a non-planar structure due to a variety of intramolecular interactions. At the time of writing this paper, data unambiguously define the spatial and electronic structures of fagopyrins. An explanation of this may be found in difficult and multi-stage processes needed to obtain the pure substance from plant material [15]. Additionally, the existence of unstable protofagopyrins [16,17] and the possibility of the existence of many derivatives [17,18] can be a problem.
Conformers representing different arrangements and interactions of substituent hydroxyl groups, carbonyl oxygen, piperidine, and pyrrolidine rings were compiled based on a structure with minimum energy. The close position of the substituents allowed for the formation of intramolecular interactions at the "peri"

Conformational Analysis of Fagopyrins
The structure of fagopyrin consists of a polycyclic system of eight rings. As shown in Scheme 2, each ring, named A-H , consists of six carbon atoms. Rings A-H are characterized by the presence of at least one substituent. In the A + B + C and F + G + H regions of the molecule, there are two hydroxyl groups and carbonyl oxygen. Such a close position of the substituents allows the formation of a hydrogen bond system in which hydroxyl groups are directed to the centrally located carbonyl oxygen. As shown for the hypericin molecule [13], such formation of strong hydrogen bonds is energetically preferred, and these hydrogen bonds are difficult to break. What is new in the structure of fagopyrin is Conformers representing different arrangements and interactions of substituent hydroxyl groups, carbonyl oxygen, piperidine, and pyrrolidine rings were compiled based on a structure with minimum energy. The close position of the substituents allowed for the formation of intramolecular interactions at the "peri" (O-H···O···H-O) and "bay" region of the molecules (O-H···O-H). Additionally, the presence of piperidine and pyrrolidine substituents allowed for the formation of new interactions with the nitrogen atom.

Conformational Analysis of Fagopyrins
The structure of fagopyrin consists of a polycyclic system of eight rings. As shown in Scheme 2, each ring, named A-H, consists of six carbon atoms. Rings A-H are characterized by the presence of at least one substituent. In the A + B + C and F + G + H regions of the molecule, there are two hydroxyl groups and carbonyl oxygen. Such a close position of the substituents allows the formation of a hydrogen bond system in which hydroxyl groups are directed to the centrally located carbonyl oxygen. As shown for the hypericin molecule [13], such formation of strong hydrogen bonds is energetically preferred, and these hydrogen bonds are difficult to break. What is new in the structure of fagopyrin is the close position of piperidine and pyrrolidine substituents. These substituents can occur at positions named R1 and R2 (Scheme 2). The rings containing the substituents with nitrogen atoms give an additional possibility to form of intramolecular OHN hydrogen bond and possible breaking of the strong OHO hydrogen bonds. So far, no studies have been found on the arrangement of these substituents in fagopyrin molecules. Another interesting part of the fagopyrin molecule is the "bay" region consisting of the A + D + F ring system. In hypericin molecule, the preferred arrangement of the substituents in the "bay" region is to form an OHO hydrogen bond between two hydroxyl groups. In the fagopyrin molecule, the addition of piperidine and pyrrolidine rings allows the interactions to be directed to the nitrogen atom forming new OHN interactions. Another part of fagopyrin molecule that may affect the overall structure is the presence of R3 and R4 substituents. Depending on the type of fagopyrin A-F molecule, these parts can be substituted by protons or methyl groups. As shown for hypericin, the close distance of two methyl groups may cause strain in the entire molecule and can strongly affect the planarity of the polycyclic system. Such a variety of substituents and possible strain effects from methyl groups make the structure and intramolecular interactions in fagopyrin molecules worth describing. the close position of piperidine and pyrrolidine substituents. These substituents can occur at positions named R1 and R2 (Scheme 2). The rings containing the substituents with nitrogen atoms give an additional possibility to form of intramolecular OHN hydrogen bond and possible breaking of the strong OHO hydrogen bonds. So far, no studies have been found on the arrangement of these substituents in fagopyrin molecules. Another interesting part of the fagopyrin molecule is the "bay" region consisting of the A + D + F ring system. In hypericin molecule, the preferred arrangement of the substituents in the "bay" region is to form an OHO hydrogen bond between two hydroxyl groups. In the fagopyrin molecule, the addition of piperidine and pyrrolidine rings allows the interactions to be directed to the nitrogen atom forming new OHN interactions. Another part of fagopyrin molecule that may affect the overall structure is the presence of R3 and R4 substituents. Depending on the type of fagopyrin A -F molecule, these parts can be substituted by protons or methyl groups. As shown for hypericin, the close distance of two methyl groups may cause strain in the entire molecule and can strongly affect the planarity of the polycyclic system. Such a variety of substituents and possible strain effects from methyl groups make the structure and intramolecular interactions in fagopyrin molecules worth describing. Six low-energy conformers of fagopyrin A have been obtained, and structure 2 is the minimum energy conformer ( Figure 1). The hydrogen bonds in the "peri" region of the molecule show alignment to the carbonyl oxygen. Two OHN hydrogen bonds in the "bay" region are preferred; however, breaking one OHN hydrogen bond in the "bay" region and forming OHO hydrogen bond between hydroxyl groups results in a total energy change of only 3.8 kcal•mol −1 (structure 1). A similar change in the total energy of the system (ΔE ≈ 4.9 kcal•mol −1 ) is caused by breaking a strong OHO hydrogen bond in the "peri" region and the formation of an OHN bond with the pyrrolidine substituent (structure 3). Breaking the OHO hydrogen bond in the "peri" region without the formation of another interaction destabilizes the fagopyrin structure and raises its energy (structure 4, 6). Six low-energy conformers of fagopyrin A have been obtained, and structure 2 is the minimum energy conformer ( Figure 1). The hydrogen bonds in the "peri" region of the molecule show alignment to the carbonyl oxygen. Two OHN hydrogen bonds in the "bay" region are preferred; however, breaking one OHN hydrogen bond in the "bay" region and forming OHO hydrogen bond between hydroxyl groups results in a total energy change of only 3.8 kcal·mol −1 (structure 1). A similar change in the total energy of the system (∆E ≈ 4.9 kcal·mol −1 ) is caused by breaking a strong OHO hydrogen bond in the "peri" region and the formation of an OHN bond with the pyrrolidine substituent (structure 3). Breaking the OHO hydrogen bond in the "peri" region without the formation of another interaction destabilizes the fagopyrin structure and raises its energy (structure 4, 6).
For fagopyrin B, six low-energy conformers ( Figure 2) have been obtained. Structure 8 showing the lowest energy is characterized by the OHO hydrogen bond arrangement in the "peri" region typical for anthrones. In the "bay" region, the OHN hydrogen bonds linking the hydroxyl group and the nitrogen atom are formed. The energy differences between the structures 7, 8, and 9 show the energy difference up to 10.0 kcal·mol −1 . The energy difference for these conformers is larger than the analogous difference for structures 1, 2, and 3 of fagopyrin A. Formation of OHN hydrogen bonds with the piperidine ring (fagopyrin A) shows larger energy differences than the formation of OHN interactions with the pyrrolidine ring (fagopyrin B). Additionally, it can be seen that the piperidine ring in the fagopyrin B prefers a "chair" conformation; however, interaction with the hydroxyl substituent in the "peri" and "bay" region can disrupt the chair conformation (structure [7][8][9][10][11][12]. The presence of a free hydroxyl group (structure 10, 12) results in a significant increase in the energy of fagopyrin B. In contrast, the lack of the methyl groups brings the double anthrone system closer to planarity.  For fagopyrin B, six low-energy conformers ( Figure 2) have been obtained. Structure 8 showing the lowest energy is characterized by the OHO hydrogen bond arrangement in the "peri" region typical for anthrones. In the "bay" region, the OHN hydrogen bonds linking the hydroxyl group and the nitrogen atom are formed. The energy differences between the structures 7, 8, and 9 show the energy difference up to 10.0 kcal•mol −1 . The energy difference for these conformers is larger than the analogous difference for structures 1, 2, and 3 of fagopyrin A. Formation of OHN hydrogen bonds with the piperidine ring (fagopyrin A) shows larger energy differences than the formation of OHN interactions with the pyrrolidine ring (fagopyrin B). Additionally, it can be seen that the piperidine ring in the fagopyrin B prefers a "chair" conformation; however, interaction with the hydroxyl substituent in the "peri" and "bay" region can disrupt the chair conformation (structure 7-12). The presence of a free hydroxyl group (structure 10, 12) results in a significant increase in the energy of fagopyrin B. In contrast, the lack of the methyl groups brings the double anthrone system closer to planarity.  Six conformers that were obtained for fagopyrin C (Figure 3) are characterized by low energy. The structure with the lowest energy (structure 14) favors the formation of OHN hydrogen bonds and the breaking of the OHO hydrogen bonds in the "bay" region of the molecule. The piperidine ring shows a "chair" conformation for all the obtained structures (structure 13-18). Breaking of the OHN hydrogen bond located in the "bay" region results in leaving the piperidine ring free and increasing the energy of the molecule by 7.0 kcal•mol −1 (structure 13). Breaking of the OHO hydrogen bond in the "bay" region OHN hydrogen bonds and the breaking of the OHO hydrogen bonds in the "bay" region of the molecule. The piperidine ring shows a "chair" conformation for all the obtained structures (structure [13][14][15][16][17][18]. Breaking of the OHN hydrogen bond located in the "bay" region results in leaving the piperidine ring free and increasing the energy of the molecule by 7.0 kcal·mol −1 (structure 13). Breaking of the OHO hydrogen bond in the "bay" region together with the formation of the OHN hydrogen bond with a hydroxyl group located in the "peri" region (structure 15) is associated with an energy increase of 9.6 kcal·mol −1 . Leaving the "free" hydroxyl group in the "peri" region results in a significant increase in the energy ∆E ≈ 27.7 kcal·mol −1 (structure 16) and ∆E ≈ 65.8 kcal·mol −1 (structure 18).  Six low-energy conformers ( Figure 4) were obtained for fagopyrin D. The lowest energy structure (structure 20) is characterized by the formation of an OHN hydrogen bond in the "bay" region. Structure 19 is characterized by a "hypericin-like" arrangement of the hydroxyl groups in the "bay", and the "peri" region differs in energy by 7.0 kcal•mol −1 from the lowest energy structure. The "chair" conformation is preferred for both piperidine rings in fagopyrin D. Structure 21 is characterized by the breaking of the strong OHO hydrogen bond in the "peri" region and the formation of an OHN hydrogen bond to the piperidine ring. Such transfer of the hydrogen interaction results in the energy difference of 9.6 kcal•mol −1 to the minimum energy structure (structure 20). As in the fagopyrin A-C structure, the "free" hydroxyl group (22,24) increases the energy of the conformer; however, in such a polycyclic system, this may not be a direct expression of breaking the OHN hydrogen bond but also due to possible structural changes of the multi-ring molecule. The formation of the OHN hydrogen bond in the "peri" region stabilizes the fagopyrin D molecule. Six low-energy conformers ( Figure 4) were obtained for fagopyrin D. The lowest energy structure (structure 20) is characterized by the formation of an OHN hydrogen bond in the "bay" region. Structure 19 is characterized by a "hypericin-like" arrangement of the hydroxyl groups in the "bay", and the "peri" region differs in energy by 7.0 kcal·mol −1 from the lowest energy structure. The "chair" conformation is preferred for both piperidine rings in fagopyrin D. Structure 21 is characterized by the breaking of the strong OHO hydrogen bond in the "peri" region and the formation of an OHN hydrogen bond to the piperidine ring. Such transfer of the hydrogen interaction results in the energy difference of 9.6 kcal·mol −1 to the minimum energy structure (structure 20). As in the fagopyrin A-C structure, the "free" hydroxyl group (22,24) increases the energy of the conformer; however, in such a polycyclic system, this may not be a direct expression of breaking the OHN hydrogen bond but also due to possible structural changes of the multi-ring molecule. The formation of the OHN hydrogen bond in the "peri" region stabilizes the fagopyrin D molecule.  Six conformers ( Figure 5) were obtained for fagopyrin E. The lowest-energy conformer (structure 26) shows hydrogen bonding in the "bay" region of the molecule. The OHN hydrogen bonds are formed by the hydroxyl groups to both nitrogen atoms in the piperidine and pyrrolidine substituent. In the minimum-energy conformer, the hydrogen bonds in the "peri" region are directed to the carbonyl oxygen. The "chair" conformation of the piperidine substituent is preferred. Conformer characterized by the "free" piperidine group (structure 25) differs in the energy of 7.5 kcal•mol −1 . Additionally, breaking of OHO hydrogen bond in the "peri" region and transferring it to the "free" piperidine substituent (structure 27) raises the energy relative to conformer 25 by 2.1 kcal•mol −1 . As in the fagopyrin structures described previously, breaking of a strong OHO hydrogen bond in the "peri" region and leaving the hydroxyl group unbound raises the total energy of the polycyclic system (structure 28 and structure 30). Six conformers ( Figure 5) were obtained for fagopyrin E. The lowest-energy conformer (structure 26) shows hydrogen bonding in the "bay" region of the molecule. The OHN hydrogen bonds are formed by the hydroxyl groups to both nitrogen atoms in the piperidine and pyrrolidine substituent. In the minimum-energy conformer, the hydrogen bonds in the "peri" region are directed to the carbonyl oxygen. The "chair" conformation of the piperidine substituent is preferred. Conformer characterized by the "free" piperidine group (structure 25) differs in the energy of 7.5 kcal·mol −1 . Additionally, breaking of OHO hydrogen bond in the "peri" region and transferring it to the "free" piperidine substituent (structure 27) raises the energy relative to conformer 25 by 2.1 kcal·mol −1 . As in the fagopyrin structures described previously, breaking of a strong OHO hydrogen bond in the "peri" region and leaving the hydroxyl group unbound raises the total energy of the polycyclic system (structure 28 and structure 30).
For fagopyrin F, six low-energy conformers were obtained. The lowest-energy structure again is characterized by the formation of the OHN hydrogen bonds in the "bay" region (structure 32). The chair conformation of the piperidine substituents is preferred. The energetically similar conformers 31 and 33 are characterized by an energy difference of 7.5 and 9.7 kcal·mol −1 , relatively to the minimum. As in conformers of fagopyrin E, it is possible to break the OHN hydrogen bond in the "bay" region and form an OHN hydrogen bond in the "peri" region. Breaking of the strong OHO hydrogen bond system in the "peri" region causes the deformation of the polycyclic system and deviates the molecule from planarity (36).
In general, the structure of fagopyrin tends to form OHN hydrogen bonds in the "bay" region. Energetically preferred formation of strong OHO hydrogen bonds to carbonyl oxygen in the "peri" region is evident in most conformers, and breaking of these interactions has the consequence of raising the energy of the system. Nevertheless, it is possible to break the strong OHO hydrogen bonds in the "peri" region in favor of the formation of an OHN hydrogen bond with the piperidine or pyrrolidine substituent. In summary, the introduction of piperidine and pyrrolidine substituents into the hypericin system provides an opportunity to form an OHN hydrogen bond instead of the strongest OHO.  For fagopyrin F, six low-energy conformers were obtained. The lowest-energy structure again is characterized by the formation of the OHN hydrogen bonds in the "bay" region (structure 32). The chair conformation of the piperidine substituents is preferred. The energetically similar conformers 31 and 33 are characterized by an energy difference of 7.5 and 9.7 kcal•mol −1 , relatively to the minimum. As in conformers of fagopyrin E, it is possible to break the OHN hydrogen bond in the "bay" region and form an OHN hydrogen bond in the "peri" region. Breaking of the strong OHO hydrogen bond system in the "peri" region causes the deformation of the polycyclic system and deviates the molecule from planarity (36).
In general, the structure of fagopyrin tends to form OHN hydrogen bonds in the "bay" region. Energetically preferred formation of strong OHO hydrogen bonds to carbonyl oxygen in the "peri" region is evident in most conformers, and breaking of these interactions has the consequence of raising the energy of the system. Nevertheless, it is possible to break the strong OHO hydrogen bonds in the "peri" region in favor of the formation of an OHN hydrogen bond with the piperidine or pyrrolidine substituent. In summary, the introduction of piperidine and pyrrolidine substituents into the hypericin system provides an opportunity to form an OHN hydrogen bond instead of the strongest OHO.

Analysis of Geometry of Fagopyrin Structures
A parameter that describes the geometry of the fagopyrin conformers is the angle between the planes formed by the peripheral rings A-C, F-H, A-F, and C-H (Scheme 2). For hypericin, (Table 1)

Analysis of Geometry of Fagopyrin Structures
A parameter that describes the geometry of the fagopyrin conformers is the angle between the planes formed by the peripheral rings A-C, F-H, A-F, and C-H (Scheme 2). For hypericin, (Table 1)  . These conformers are characterized by different substitutions at the R1-R4 position but the hydroxyl groups in the "bay" and "peri" regions are oriented as in the hypericin molecule and form the same type OHO hydrogen bonds. The difference in the angle between the ring plane is the greatest for fagopyrin C and D. These fagopyrins have an asymmetric substitution with a methyl group and a proton at the R3 and R4 position.
Transfer of the OHO hydrogen bond in the "peri" region to the piperidine or pyrrolidine results in the formation of an OHN hydrogen bond (structures 3, 9, 15, 21, 27, and 33) and causes little change in the angles between the A-C and A-F planes. Formation of another OHN hydrogen bond in molecules 5, 11, 17, 23, 29, and 35 causes more significant changes in the polycyclic system. The changes are visible in the angle between F-H and C-H planes, so the effect of the OHN hydrogen bonds in the "peri" region on the geometry of the fagopyrin molecule is apparent and may have a real impact on the electron structure. Additionally, breaking the OHO hydrogen bond and leaving the hydroxyl group in the "peri" region as free causes deformation of the polycyclic system (structures 4, 10, 16, 22, 28, and 34). Larger differences can be observed when two free hydroxyl groups in the molecule are present (structures 6, 12, 18, 24, 30, and 36). This arrangement of the hydroxyl groups causes strong deformation of the polycyclic system of fagopyrins, which is reflected in the high energy of these conformers. So far, two similarities of the fagopyrin molecules to the hypericin molecule can be given. These are a strong influence on the geometrical structure of substituents at the R3 and R4 positions (methyl groups) and the preferred formation of the OHO hydrogen bonds in the "peri" region formed by hydroxyl groups and carbonyl oxygen. The structures corresponding to the energy minima (2,8,14,20,26, and 32) differ from hypericin in the "bay" region. The hydroxyl groups in the "bay" region are directed to the nitrogen atom in the piperidine and pyrrolidine substituents. In the case of hypericin, the hydroxyl groups prefer the OHO hydrogen bonding. The introduction of piperidine or pyrrolidine rings at the R1 and R2 positions favors the formation of an OHN hydrogen bond and decreasing of the fagopyrin energy to a minimum.
Changes in the angles between the plane rings of the peripheral rings of fagopyrins result from a number of substituents and the intramolecular interactions. As fagopyrin F is the major form in the plant material [18], an analysis of hydrogen bond parameters has been performed for the structures shown in Figure 6. The length and angles of OHO and OHN hydrogen bonds are summarized in Table 2. The parameters of OHN hydrogen bonds directed to the pyrrolidine and piperidine rings are highlighted in bold. The results calculated for the fagopyrin F conformers have been compared with hypericin. The letters in parentheses in Table 2 identify the hydrogen bond location described according to Scheme 2. Changes in the angles between the plane rings of the peripheral rings of fagopyrins result from a number of substituents and the intramolecular interactions. As fagopyrin F is the major form in the plant material [18], an analysis of hydrogen bond parameters has been performed for the structures shown in Figure 6. The length and angles of OHO and OHN hydrogen bonds are summarized in Table 2. The parameters of OHN hydrogen bonds directed to the pyrrolidine and piperidine rings are highlighted in bold. The results calculated for the fagopyrin F conformers have been compared with hypericin. The letters in parentheses in Table 2 identify the hydrogen bond location described according to Scheme 2.     Structure 31 is characterized by an arrangement of substituents similar to hypericin. The length of hydrogen bonds in the "peri" region is similar to the length of analogous bonds in hypericin. In the "bay" region, the OHO bond length is shorter than in the hypericin molecule. The OHN hydrogen bond is characterized by a length of 1.5386 Å and an angle of 156.683 • . The lowest energy conformer (structure 32) is characterized by "peri"-OHO hydrogen bond lengths similar to hypericin. In structure 32, two OHN hydrogen bonds are presented in the "bay" region. The bond labeled as C(F)-O-H····N(R2) is characterized by length and angle similar to the OHO "peri" bonds. The bond labeled as C(A)-O-H···N(R1) is elongated up to 1.7233 Å. Additionally, the "peri" OHN hydrogen bond in structure 33 is longer (1.7164 Å) than the "peri" OHO hydrogen bonds.

Aromaticity of Fagopyrin
The aromaticity of polycyclic compounds is related to their structure and reactivity. There are many indices describing aromaticity; however, the classical HOMA index (Harmonic Oscillator Measure of Aromaticity) is convenient for the description of aromaticity in organic compounds and relates it directly to the structure [36]. In typical aromatic compounds, the values of the HOMA index are in the range from 0 to 1, where 0 corresponds to a non-aromatic ring, and the value of 1 corresponds to a fully delocalized benzene structure, and only in special cases, the HOMA value can exceed the 0-1 range. Figure 7  The HOMA values of rings A and F are lower relatively to the hypericin; thus, the presence of the piperidine ring decreases the aromaticity of the rings. However, the D ring in structure 31 gains aromaticity relative to the hypericin moiety. The formation of the strong OHO hydrogen bonds in the "peri" region stabilizes the polycyclic system and increases the aromaticity of the central B and G rings. In general, structure 31 shows aromaticity of the rings similar to hypericin, with the influence of piperidine substituents on the aromaticity of rings A, D, and F. The lowest energy structure 32 is characterized by the formation of two OHN hydrogen bonds to the piperidine substituent in the "bay" region. This arrangement increases the aromaticity of the A, D, and F ring. Structure 33 is characterized by the breaking of the strong OHO bond in the "peri" region and the formation of an OHN hydrogen bond to the piperidine ring. Such conformation causes an increase in the energy of the system, an increase in the HOMA value of the F ring, and a decrease in the HOMA value of the B ring up to 0.1644. Breaking of another OHO bond in the "peri" region (structure 34) deepens the loss of aromaticity of the B ring. The HOMA parameter below zero indicates a complete loss of aromaticity of the ring. In structures 35 and 36, two OHN bonds in the "peri" region are present and such conformation of the hydroxyl groups causes an increase in the HOMA value in the A and F ring with significant aromaticity decreasing in the B and G ring. These changes cause increasing the total energy of the molecular system ( Figure 6).

Analysis of Intramolecular Interactions in Fagopyrin Derivatives
Changes in aromaticity must be related to the changes in electron density of the polycyclic system. To describe the possible intermolecular interactions and arrangement of the electron density, the QTAIM [37] (Quantum Theory of Atoms in Molecules) analysis for the fagopyrin F conformers has been performed. In the frame of the QTAIM theory, a molecule consists of maximum, minimum, and saddle points of the electron density ρ(r). The saddle point indicates bond-critical points (BCPs) or ring-critical points (RCPs). The points representing the maximum electron density correspond to atoms. Figure 8 shows QTAIM graphs of fagopyrin F conformers. The structure of fagopyrin F is characterized by the presence of numerous substituents in a polycyclic system. Such structure allows for the occurrence of numerous intramolecular interactions of diverse nature [38]. QTAIM analysis confirms the presence of the hydrogen bond interactions in the "peri" region of fagopyrin F. Hydroxyl groups directed to carbonyl oxygen form a stable moiety as in the case of hypericin and sennidin [13,14]. QTAIM analysis indicates that OHN interaction can be formed in both the "peri" and "bay" region. The strong OHO hydrogen bonds in the anthrone moiety can be broken and replaced by weaker OHN hydrogen bonds. The electron density values at the bond critical points ρ(r) presented in Figure 8 reflect the strength of the OHO and OHN hydrogen bonds. The structures 33, 34, 35, and 36 show the formation of OHN hydrogen bonds characterized by lower values of ρ relative to the OHO hydrogen bonds in the "peri" region. However, conformer 32 (b) is characterized by the formation of OHN hydrogen bonds in the "bay" region. These bonds are characterized by a similar value of ρ(r) relative to the strong OHO hydrogen bonds in the "peri" region. This conformer shows the lowest energy; thus, the formation of strong OHN hydrogen bonds to the piperidine substituents stabilizes the anthrone system. The proximity of the piperidine ring affects the adjacent hydroxyl groups even if they do not form a

Analysis of Intramolecular Interactions in Fagopyrin Derivatives
Changes in aromaticity must be related to the changes in electron density of the polycyclic system. To describe the possible intermolecular interactions and arrangement of the electron density, the QTAIM [37] (Quantum Theory of Atoms in Molecules) analysis for the fagopyrin F conformers has been performed. In the frame of the QTAIM theory, a molecule consists of maximum, minimum, and saddle points of the electron density ρ(r). The saddle point indicates bond-critical points (BCPs) or ring-critical points (RCPs). The points representing the maximum electron density correspond to atoms. Figure 8 shows QTAIM graphs of fagopyrin F conformers. The structure of fagopyrin F is characterized by the presence of numerous substituents in a polycyclic system. Such structure allows for the occurrence of numerous intramolecular interactions of diverse nature [38]. QTAIM analysis confirms the presence of the hydrogen bond interactions in the "peri" region of fagopyrin F. Hydroxyl groups directed to carbonyl oxygen form a stable moiety as in the case of hypericin and sennidin [13,14]. QTAIM analysis indicates that OHN interaction can be formed in both the "peri" and "bay" region. The strong OHO hydrogen bonds in the anthrone moiety can be broken and replaced by weaker OHN hydrogen bonds. The electron density values at the bond critical points ρ(r) presented in Figure 8 reflect the strength of the OHO and OHN hydrogen bonds. The structures 33, 34, 35, and 36 show the formation of OHN hydrogen bonds characterized by lower values of ρ relative to the OHO hydrogen bonds in the "peri" region. However, conformer 32 (b) is characterized by the formation of OHN hydrogen bonds in the "bay" region. These bonds are characterized by a similar value of ρ(r) relative to the strong OHO hydrogen bonds in the "peri" region. This conformer shows the lowest energy; thus, the formation of strong OHN hydrogen bonds to the piperidine substituents stabilizes the anthrone system. The proximity of the piperidine ring affects the adjacent hydroxyl groups even if they do not form a direct bond. In addition, the close position of the methyl groups also introduces intermolecular interactions. To confirm the presence of the interactions, the non-covalent interactions [32] (NCI) analysis was performed. Figure 9 shows the conformers of fagopyrin F, showing multiple intramolecular interactions. The blue isosurfaces in the "peri" and "bay" moiety confirm the presence of strong hydrogen bonds in the fagopyrin F molecule. Conformer 31 (a) shows similarity to hypericin in the formation of strong OHO hydrogen bonds. The lowest energy conformer 32 (b) confirms the formation of OHN hydrogen bonds in the "bay" region. Structure 33 (c) confirms the possibility of breaking the strong OHO hydrogen bond in favor of OHN hydrogen bond formation with the nitrogen atom of the piperidine ring. The interactions between the methyl groups can be described as dispersive. To confirm the presence of the interactions, the non-covalent interactions [32] (NCI) analysis was performed. Figure 9 shows the conformers of fagopyrin F, showing multiple intramolecular interactions. The blue isosurfaces in the "peri" and "bay" moiety confirm the presence of strong hydrogen bonds in the fagopyrin F molecule. Conformer 31 (a) shows similarity to hypericin in the formation of strong OHO hydrogen bonds. The lowest energy conformer 32 (b) confirms the formation of OHN hydrogen bonds in the "bay" region. Structure 33 (c) confirms the possibility of breaking the strong OHO hydrogen bond in favor of OHN hydrogen bond formation with the nitrogen atom of the piperidine ring. The interactions between the methyl groups can be described as dispersive.

UV-VIS Spectra of Fagopyrin Conformers
Different conformation of the investigated fagopyrins is reflected in their electron structures. In Figure 10 are presented the HOMO and LUMO orbitals for the conformers of fagopyrin F-the most popular in the plant material. For other fagopyrins, the HOMO and LUMO orbitals are collected in Supplementary Materials. It is characteristic that for all conformers, the HOMO orbital is located mainly on the outer A, C, F, and H rings and on the oxygen atoms of the hydroxyl group. Only for conformers 34 and 36, the HOMO orbital is more concentrated on rings A and F than on C and H. The transfer of electrons to the LUMO orbital is connected with the shifting of electrons to the B and G rings, the oxygen of the carbonyl group, and outer bonds of the A, C, F, and H rings. Since the arrangement of the HOMO and LUMO orbitals is similar for all conformers, the HOMO-LUMO gap energy is also similar. This is persistent for all the analyzed fagopyrins. In Table 3 are collected the HOMO-LUMO gap energies for all the analyzed fagopyrins.

UV-VIS Spectra of Fagopyrin Conformers
Different conformation of the investigated fagopyrins is reflected in their electron structures. In Figure 10 are presented the HOMO and LUMO orbitals for the conformers of fagopyrin F -the most popular in the plant material. For other fagopyrins, the HOMO and LUMO orbitals are collected in Supplementary Materials. It is characteristic that for all conformers, the HOMO orbital is located mainly on the outer A, C, F, and H rings and on the oxygen atoms of the hydroxyl group. Only for conformers 34 and 36, the HOMO orbital is more concentrated on rings A and F than on C and H. The transfer of electrons to the LUMO orbital is connected with the shifting of electrons to the B and G rings, the oxygen of the carbonyl group, and outer bonds of the A, C, F, and H rings. Since the arrangement of the HOMO and LUMO orbitals is similar for all conformers, the HOMO-LUMO gap energy is also similar. This is persistent for all the analyzed fagopyrins. In Table 3 are collected the HOMO -LUMO gap energies for all the analyzed fagopyrins.     The UV spectra for fagopyrin F shown in Figure 11 are characterized by the presence of two intense bands. For the conformers 31, 32, and 33, the most intensive band shifts from 556 nm to 552 nm. The second intensive band is located at 457, 455, and 459 nm. For conformer 34, except for the most intensive band at 546 nm, two bands with similar intensity at 477 and 427 nm are present. The last band at 427 nm is visible in the UV-VIS spectra of 31, 32, 34, and 36 conformers; however, it is significantly lower compared to other bands. For conformers 31, 32, 33, and 34, the most intensive band is related to HOMO-LUMO transition. For conformer 35, this band is shifted to 539 nm, for conformer 36 to 503 nm, and the intensity of this band is lower than the bands at 437 and 421, respectively. For the conformers 35 and 36, transitions from lower energy orbitals are more intense than for the HOMO-LUMO transition. The electron transition participating in the bands for fagopyrin F conformers are collected in Table 4. The shape of the orbitals involved in the electron transitions in structure 32 (the lowest energy structure of fagopyrin F) is shown in Figure 12. The shape of the orbitals involved in UV-VIS transitions for fagopyrin F conformers except the presented in the text (Figure 12) is shown in Supplementary Materials.  A comparison of calculated fagopyrin F and experimental spectra [19] for the plant material suggests that in the plant material, many fagopyrin structures may be present. It is not clear which version of fagopyrin in the experimental spectra was registered; how-  A comparison of calculated fagopyrin F and experimental spectra [19] for the plant material suggests that in the plant material, many fagopyrin structures may be present. It is not clear which version of fagopyrin in the experimental spectra was registered; however, in the experimental spectra, the most intensive theoretically calculated bands are visible.

Conclusions
Theoretical calculations can provide information on the molecular structure when the structure is unknown, which is often the case with plant material. Fagopyrin compounds may exist as conformers characterized by a different energy. The presence of the piperidine and pyrrolidine ring in fagopyrin introduces novel intramolecular interactions compared to the double anthrone molecules. Fagopyrin A-F structures are characterized by the presence of a number of substituents and strong hydrogen bonds in the anthrone moiety. Although the OHO hydrogen bonds in the anthrone moiety are characterized as very strong, both the OHO and OHN hydrogen bonds may exist in the fagopyrin A-F structure. It is possible to break the strong OHO hydrogen bonds in the anthrone moiety in favor of interactions with the nitrogen atom in piperidine or pyrrolidine substituent. Changes in the molecular geometry are related to the changes in the orbital localization, which is reflected in the UV-VIS spectra of fagopyrin conformers.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.