Identification of Potential Allosteric Site Binders of Indoleamine 2,3-Dioxygenase 1 from Plants: A Virtual and Molecular Dynamics Investigation

Ligand and structure-based computational screenings were carried out to identify flavonoids with potential anticancer activity. Kushenol E, a flavonoid with proven anticancer activity and, at the same time, an allosteric site binder of the enzyme indoleamine 2,3-dioxygenase-1 (IDO1), was used as the reference compound. Molecular docking and molecular dynamics simulations were performed for the screened flavonoids with known anticancer activity. The following two of these flavonoids were identified as potential inhibitors of IDO1: dichamanetin and isochamanetin. Molecular dynamics simulations were used to assess the conformational profile of IDO1-flavonoids complexes, as well as for calculating the bind-free energies.


Tryptophan Metabolism and the Kynurenine Pathway
Studies show that tumor tissues within a chronic inflammatory context (presence of pro-inflammatory cytokines for a long time) catabolize in large proportions the amino acid L-tryptophan in the cell matrix and synthesize metabolites from one of the main degradation pathways of tryptophan, the pathway of kynurenine [1,2].
Almost all metabolites of the kynurenine pathway affect immune activity through various mechanisms. Figure 1 illustrates the effects of kynurenine pathway metabolites on the immune system. The determining step of the kynurenine pathway is the formation of N-formylkynurenine, regulated by the family of dioxygenase enzymes known as tryptophan 2,3-dioxygenase (TDO) and indolamine 2,3-dioxygenase (IDO), which may exist in the following two isoforms: IDO1 and IDO2 [3].
TDO, IDO1, and IDO2 are heme-bearing enzymes and catalyze the rate-limiting step of tryptophan catabolization. TDO is mainly expressed in the liver and oxidizes excess tryptophan in the body, generating ATP and especially NAD+. The expression of TDO is stimulated by the concentration of blood tryptophan as well as by the synthesis of the heme cofactor [4]. IDO1 can be expressed by many different cells, including antigen-presenting cells (APCs), such as monocyte-derived macrophages, dendritic cells (DCs), and fibroblasts. Its expression is mainly induced by inflammatory cytokines, such as IFN-γ, TNF-α, IL-1, and IL-2, secreted by Th1-type lymphocytes, as well as TGF-β, IL-10, and adenosine secreted by regulatory T lymphocytes [5]. IDO1 expression is further stimulated by its own product, kynurenine, via the aryl hydrocarbon receptor (AhR) [6,7].   IDO2 is significantly less active than IDO1 [8]. Similar to IDO1, its expression is stimulated by the activation of AhR [9]. Although IDO2 is expressed by cancer cells, when compared to IDO1, it contributes little to the accumulation of kynurenine pathway metabolites [10,11].
In healthy patients, IDO1 expression is restricted to endothelial cells in the placenta, lung, mature dendritic cells in secondary lymphoid organs, and in epithelial cells scattered in the female genital tract [12]. Whereas, under inflammatory conditions, IDO1 expression is strongly induced by INF-γ [13]. Furthermore, it is highly expressed at inflammatory sites, where it may contribute to negative feedback against local immune response activity [12]. Inflamed tumors may also express IDO1 as an adaptive resistance mechanism, in which it is induced by IFN-γ produced by dendritic cells and tumor-infiltrated macrophages [14].
Given the physiological importance of the IDO1 enzyme within the tumor microenvironment, the need for structural studies and the development of increasingly efficient inhibitors has become an attractive approach to treating cancer patients.

IDO1: Structure
The human IDO1 is a monomeric α-helical dioxygenase enzyme containing an N-terminal minor domain (residues 1-154) and a C-terminal major domain (residues 155-403), where the catalytic site and a cofactor heme are located ( Figure 2). IDO2 is significantly less active than IDO1 [8]. Similar to IDO1, its expression is stimulated by the activation of AhR [9]. Although IDO2 is expressed by cancer cells, when compared to IDO1, it contributes little to the accumulation of kynurenine pathway metabolites [10,11].
In healthy patients, IDO1 expression is restricted to endothelial cells in the placenta, lung, mature dendritic cells in secondary lymphoid organs, and in epithelial cells scattered in the female genital tract [12]. Whereas, under inflammatory conditions, IDO1 expression is strongly induced by INF-γ [13]. Furthermore, it is highly expressed at inflammatory sites, where it may contribute to negative feedback against local immune response activity [12]. Inflamed tumors may also express IDO1 as an adaptive resistance mechanism, in which it is induced by IFN-γ produced by dendritic cells and tumorinfiltrated macrophages [14].
Given the physiological importance of the IDO1 enzyme within the tumor microenvironment, the need for structural studies and the development of increasingly efficient inhibitors has become an attractive approach to treating cancer patients.

IDO1: Structure
The human IDO1 is a monomeric α-helical dioxygenase enzyme containing an Nterminal minor domain (residues 1-154) and a C-terminal major domain (residues 155-403), where the catalytic site and a cofactor heme are located ( Figure 2). The minor domain or N-terminal domain, shown in green in Figure 2, contains residues Gly11 to Asp158, consisting of six α-helices (involving amino acids Tyr36 to His45, His45 to Ser52, Gln54 to Lys61, Asp72 to Gly93, Pro104 to Glu119 and Val125 to Val130 or helix A, contributing to the active site), four 310-helices (Ser12 to His16, Pro33 to Phe35, Ser66 to Leu70 and Thr144 to Glu146) and a small antiparallel β-sheet formed by residues Lys135 to Lys136 and Met148 to Asp149. Such a domain makes up part of the upper region of the catalytic site and acts as structural support for the enzyme [15].
The major domain or C-terminal domain, shown in grey in Figure 2, contains residues Cys159 to Gly403, being constituted by 9 α-helices [15]. This domain performs functions related to conversion, positioning, and displacement of substrates/products in the enzyme and is where the catalytic site is located. There are also regions with specific functionalities that are important for the enzymatic reaction to occur. Among them, there is the "access" tunnel for the entrance of small molecules from the external environment The minor domain or N-terminal domain, shown in green in Figure 2, contains residues Gly11 to Asp158, consisting of six α-helices (involving amino acids Tyr36 to His45, His45 to Ser52, Gln54 to Lys61, Asp72 to Gly93, Pro104 to Glu119 and Val125 to Val130 or helix A, contributing to the active site), four 3 10 -helices (Ser12 to His16, Pro33 to Phe35, Ser66 to Leu70 and Thr144 to Glu146) and a small antiparallel β-sheet formed by residues Lys135 to Lys136 and Met148 to Asp149. Such a domain makes up part of the upper region of the catalytic site and acts as structural support for the enzyme [15].
The major domain or C-terminal domain, shown in grey in Figure 2, contains residues Cys159 to Gly403, being constituted by 9 α-helices [15]. This domain performs functions related to conversion, positioning, and displacement of substrates/products in the enzyme and is where the catalytic site is located. There are also regions with specific functionalities that are important for the enzymatic reaction to occur. Among them, there is the "access" tunnel for the entrance of small molecules from the external environment (ligand delivery Pharmaceuticals 2022, 15, 1099 4 of 18 tunnel) and the JK loop (red), with the function of conducting and maintaining the substrate (tryptophan) inside the catalytic site [16].
IDO1 has the following two non-competitive natural substrates: tryptophan and molecular oxygen (O 2 ). Moreover, the enzyme has the following four states: • apo-IDO1: without the heme cofactor; • holo-IDO1: containing the iron atom of the heme group in the ferric state (Fe 3+ ); • holo-IDO1 with the iron atom of the heme group in the ferrous state (Fe 2+ ); Despite not being the active form of the enzyme, apo-IDO1 is a promising macromolecular target since cellular studies have shown that this form is the most abundant in the cell medium, about 85% [18]. Therefore, it is likely that there is a binding site where potential inhibitors of the apo form can interact with this form of the enzyme. The main objective of this project was to apply molecular modeling methods in order to investigate this possibility.

IDO1: Inhibitors
Over the past decade, the scientific academy has been making efforts for the development of IDO1 inhibitors [2]. Röhrig and collaborators classified existing test-phase inhibitors into the following four categories: type I, II, III, and IV inhibitors [19].
Type I inhibitors are classified as competitive inhibitors of tryptophan. In general, type I inhibitors target the molecular oxygen-bound holo-IDO1 form and do not form a direct bond with heme iron [17]. The 1-LMT and 1-DMT are classified within this category of inhibitors [20].
Type II inhibitors are classified as competitive inhibitors of molecular oxygen by binding to the iron atom in its ferrous state of the heme group, acting in the ferrous holo-IDO1 form. Inhibitors such as β-carboline [21], and Epacadostat [22] fall into this category.
Type III inhibitors are noncompetitive inhibitors of tryptophan. New inhibitors in this class are being studied, such as 4PI, Navoximod, and MMG-0358 [19].
Type IV inhibitors are the most recently described class of inhibitors in the literature, including inhibitors such as the clinical compound BMS-986205, which targets the apo-IDO1 form [23].
The structures of the mentioned inhibitors are shown in Figure 3. (ligand delivery tunnel) and the JK loop (red), with the function of conducting and maintaining the substrate (tryptophan) inside the catalytic site [16]. Finally, there is the DE fragment, consisting of the DE-hairpin (Ser235 to Tyr249) and the DE loop (Glu250 to Gly262) (magenta), the EF loop (Gly278 to His287) (blue), and the JK loop (Gln360 to Gly380) [17].
IDO1 has the following two non-competitive natural substrates: tryptophan and molecular oxygen (O2). Moreover, the enzyme has the following four states: • apo-IDO1: without the heme cofactor; • holo-IDO1: containing the iron atom of the heme group in the ferric state (Fe 3+ ); • holo-IDO1 with the iron atom of the heme group in the ferrous state (Fe 2+ ); • O2-linked holo-IDO1 complex.
Despite not being the active form of the enzyme, apo-IDO1 is a promising macromolecular target since cellular studies have shown that this form is the most abundant in the cell medium, about 85% [18]. Therefore, it is likely that there is a binding site where potential inhibitors of the apo form can interact with this form of the enzyme. The main objective of this project was to apply molecular modeling methods in order to investigate this possibility.

IDO1: Inhibitors
Over the past decade, the scientific academy has been making efforts for the development of IDO1 inhibitors [2]. Röhrig and collaborators classified existing test-phase inhibitors into the following four categories: type I, II, III, and IV inhibitors [19].
Type I inhibitors are classified as competitive inhibitors of tryptophan. In general, type I inhibitors target the molecular oxygen-bound holo-IDO1 form and do not form a direct bond with heme iron [17]. The 1-LMT and 1-DMT are classified within this category of inhibitors [20].
Type II inhibitors are classified as competitive inhibitors of molecular oxygen by binding to the iron atom in its ferrous state of the heme group, acting in the ferrous holo-IDO1 form. Inhibitors such as β-carboline [21], and Epacadostat [22] fall into this category.
Type III inhibitors are noncompetitive inhibitors of tryptophan. New inhibitors in this class are being studied, such as 4PI, Navoximod, and MMG-0358 [19].
Type IV inhibitors are the most recently described class of inhibitors in the literature, including inhibitors such as the clinical compound BMS-986205, which targets the apo-IDO1 form [23].
The structures of the mentioned inhibitors are shown in Figure 3.  Kwon and collaborators [24] determined the amino acid residues that define an allosteric site for the action of apo-IDO1 inhibitors (Ser12, Tyr15, Ile17, Ile178, Lys179, Ile181, Pro182, Phe185, Lys186, Phe306, Ser309, Leu310, Ser312, Asn313, and Pro314), and found that the flavonoid obtained from plants of the species Sophora flavescens, known as Kushenol E, acts as an inhibitor of this site ( Figure 4). The authors described the amino acids Pro182 and Phe185, located in the larger domain of the enzyme, as the main components of the allosteric site for inhibitory activity. Site-directed mutagenesis experiments were performed on these residues, and in vitro enzyme inhibition assays were performed. The results proved the contribution of these amino acids to enzymatic function and indicate this allosteric site as a target for new inhibitors [24]. Kwon and collaborators [24] determined the amino acid residues that define an allosteric site for the action of apo-IDO1 inhibitors (Ser12, Tyr15, Ile17, Ile178, Lys179, Ile181, Pro182, Phe185, Lys186, Phe306, Ser309, Leu310, Ser312, Asn313, and Pro314), and found that the flavonoid obtained from plants of the species Sophora flavescens, known as Kushenol E, acts as an inhibitor of this site ( Figure 4). The authors described the amino acids Pro182 and Phe185, located in the larger domain of the enzyme, as the main components of the allosteric site for inhibitory activity. Site-directed mutagenesis experiments were performed on these residues, and in vitro enzyme inhibition assays were performed. The results proved the contribution of these amino acids to enzymatic function and indicate this allosteric site as a target for new inhibitors [24]. Most IDO1 inhibitors are obtained from a natural source; however, as shown by Kwon, flavonoids may inhibit IDO1 [24]. The diverse variety of natural products can make the screening of these compounds a time-consuming and expensive task. Due to the development of more powerful computers, virtual screening of chemical compounds and pharmacodynamic validation of ligand-receptor interactions can accelerate the process of discovering potential new drugs. In this context, we focused our work on the identification of potential inhibitors of the naturally occurring apo-IDO1 allosteric site.

Ligand-Based Virtual Screening
In order to identify potential IDO1 inhibitors capable of interacting with its allosteric site, the following two virtual screening tools were used: the Mcule platform [25] and the ZINC 15 database [26]. At this stage of the work, the flavonoid Kushenol E ( Figure 5) was used as the reference compound, due to its proven effectiveness as a non-competitive inhibitor that binds to the allosteric site of the enzyme [24]. Tanimoto's coefficient was used as the sampling criterium [27], and 172 compounds were sampled.  Most IDO1 inhibitors are obtained from a natural source; however, as shown by Kwon, flavonoids may inhibit IDO1 [24]. The diverse variety of natural products can make the screening of these compounds a time-consuming and expensive task. Due to the development of more powerful computers, virtual screening of chemical compounds and pharmacodynamic validation of ligand-receptor interactions can accelerate the process of discovering potential new drugs. In this context, we focused our work on the identification of potential inhibitors of the naturally occurring apo-IDO1 allosteric site.

Ligand-Based Virtual Screening
In order to identify potential IDO1 inhibitors capable of interacting with its allosteric site, the following two virtual screening tools were used: the Mcule platform [25] and the ZINC 15 database [26]. At this stage of the work, the flavonoid Kushenol E ( Figure 5) was used as the reference compound, due to its proven effectiveness as a non-competitive inhibitor that binds to the allosteric site of the enzyme [24]. Tanimoto's coefficient was used as the sampling criterium [27], and 172 compounds were sampled. Kwon and collaborators [24] determined the amino acid residues that define an allosteric site for the action of apo-IDO1 inhibitors (Ser12, Tyr15, Ile17, Ile178, Lys179, Ile181, Pro182, Phe185, Lys186, Phe306, Ser309, Leu310, Ser312, Asn313, and Pro314), and found that the flavonoid obtained from plants of the species Sophora flavescens, known as Kushenol E, acts as an inhibitor of this site ( Figure 4). The authors described the amino acids Pro182 and Phe185, located in the larger domain of the enzyme, as the main components of the allosteric site for inhibitory activity. Site-directed mutagenesis experiments were performed on these residues, and in vitro enzyme inhibition assays were performed. The results proved the contribution of these amino acids to enzymatic function and indicate this allosteric site as a target for new inhibitors [24]. Most IDO1 inhibitors are obtained from a natural source; however, as shown by Kwon, flavonoids may inhibit IDO1 [24]. The diverse variety of natural products can make the screening of these compounds a time-consuming and expensive task. Due to the development of more powerful computers, virtual screening of chemical compounds and pharmacodynamic validation of ligand-receptor interactions can accelerate the process of discovering potential new drugs. In this context, we focused our work on the identification of potential inhibitors of the naturally occurring apo-IDO1 allosteric site.

Ligand-Based Virtual Screening
In order to identify potential IDO1 inhibitors capable of interacting with its allosteric site, the following two virtual screening tools were used: the Mcule platform [25] and the ZINC 15 database [26]. At this stage of the work, the flavonoid Kushenol E ( Figure 5) was used as the reference compound, due to its proven effectiveness as a non-competitive inhibitor that binds to the allosteric site of the enzyme [24]. Tanimoto's coefficient was used as the sampling criterium [27], and 172 compounds were sampled.

Structure-Based Virtual Screening and Molecular Docking
Besides the ligand-based virtual screening mentioned above, a structure-based virtual screening was also carried out ( Figure 6). In this approach, the previously sampled compounds, including Kushenol E, were docked into the allosteric site of the IDO1.

Structure-Based Virtual Screening and Molecular Docking
Besides the ligand-based virtual screening mentioned above, a structure-based virtual screening was also carried out ( Figure 6). In this approach, the previously sampled compounds, including Kushenol E, were docked into the allosteric site of the IDO1.  (Table S1). The range of calculated energies varied from −8.1 kcal/mol to −5.3 kcal/mol. The docking energy for Kushenol E is equal to −6.6 kcal/mol. We realized that most of the sampled compounds were not from natural sources. Since the focus of this research was the identification of substances obtained from natural products with the potential ability to inhibit the biological response of the enzyme IDO1 through interaction at the allosteric site, a new virtual screening was carried out. The strategy was to use each compound that presented interaction energies equal to or lower than −7.1 kcal/mol (19 compounds; represented in red in Supplementary Materials, Table S1) as reference compounds in new ligand-based virtual screenings. The ZINC 15 database was the chosen tool. In this context, 28 natural products were selected, which were later docked to the allosteric site of the IDO1 enzyme. The result of this second-round structure-based virtual screening is shown in Table 1. While, in this investigation, Kushenol E (shown in bold red) works as the positive control, it is also important to consider a compound that would work as a negative control. In this context, as in the paper by Kwon [24], steppogenin was chosen as such a compound. As shown in the table, among all compounds, steppogenin shows the least favorable docking energy (shown in light red).
In this second molecular docking calculation, two additional criteria were used for the final screening as follows: (a) compounds that presented docking energies equal to or lower than −7.1 kcal/mol; (b) compounds that showed some anticancer activity (as described in the literature). Consequently, the following three flavonoids were selected: dichamanetin, isochamanetin, and chamaejasmin B, shown in bold black.
Dichamanetin is a C-benzylated flavanone that can be obtained from several South- The calculated docking energies of each compound are shown in Supplementary Material (Table S1). The range of calculated energies varied from −8.1 kcal/mol to −5.3 kcal/mol. The docking energy for Kushenol E is equal to −6.6 kcal/mol. We realized that most of the sampled compounds were not from natural sources. Since the focus of this research was the identification of substances obtained from natural products with the potential ability to inhibit the biological response of the enzyme IDO1 through interaction at the allosteric site, a new virtual screening was carried out. The strategy was to use each compound that presented interaction energies equal to or lower than −7.1 kcal/mol (19 compounds; represented in red in Supplementary Materials, Table S1) as reference compounds in new ligand-based virtual screenings. The ZINC 15 database was the chosen tool. In this context, 28 natural products were selected, which were later docked to the allosteric site of the IDO1 enzyme. The result of this second-round structure-based virtual screening is shown in Table 1. While, in this investigation, Kushenol E (shown in bold red) works as the positive control, it is also important to consider a compound that would work as a negative control. In this context, as in the paper by Kwon [24], steppogenin was chosen as such a compound. As shown in the table, among all compounds, steppogenin shows the least favorable docking energy (shown in light red).
In this second molecular docking calculation, two additional criteria were used for the final screening as follows: (a) compounds that presented docking energies equal to or lower than −7.1 kcal/mol; (b) compounds that showed some anticancer activity (as described in the literature). Consequently, the following three flavonoids were selected: dichamanetin, isochamanetin, and chamaejasmin B, shown in bold black. Dichamanetin is a C-benzylated flavanone that can be obtained from several Southeast Asian plant species, such as the following: Uvaria alba, Uvaria chamae, Piper sarmentosum, and Xylopia pierrei [28][29][30]. Isochamanetin, another C-benzylated flavanone, was obtained from several plant species such as Sphaeranthus amaranthoides, Xylopia pierrei, and Ulvaria chamae [28,31,32]. Chamaejasmin B is a dimeric biflavonoid formed by two isosacuranetin subunits linked by C3-C3 carbons, and obtained from Southwest Asian plant species, such as Stellera chamaejasme L. [33]. Table 2 and Figure 7 show the intermolecular interactions formed with specific residues of the allosteric site of IDO1.
As shown in Table 2 and Figure 7, the flavonoids interact with amino acid residues from the allosteric site, as defined by Kwon [24] (Ser12, Tyr15, Ile17, Ile178, Lys179, Ile181, Pro182, Phe185, Lys186, Phe306, Ser309, Leu310, Ser312, Asn313, Pro314). Since flavonoids are compounds that have several aromatic rings, the main interactions formed with the enzyme are of the pi and van der Waals types. The screened flavonoids (Table 2) have a significant number of oxygens; consequently, they can work as hydrogen bond donors and acceptors. Furthermore, according to Kwon [24], residues Pro182 and Phe185 are essential for the structural and functional integrity of the IDO1 enzyme system.
Dichamanetin and isochamanetin A rings interact via pi-sigma interaction with Pro182. According to Kwon [24], this residue would be important for the structural and functional integrity of the IDO1 system. In chamaejasmin B, which is a biflavonoid, one of the A rings interacts with Pro182 (pi-alkyl) and Phe185 (pi-pi) with the assistance of van der Waals interactions with neighboring residues; the second ring A interacts with Phe185 (hydrogen bond) and with Pro314 (pi-alkyl), in addition to van der Waals interactions.
Focusing on substituents attached to the A ring of flavonoids, some structural features must be considered. The hydroxyl group attached to carbon 7 of the A rings of isochamanetin, dichamanetin, and Kushenol E forms strong hydrogen bond interactions with Ser12. Phenyl groups linked to the A ring of dichamanetin and isochamanetin probably help in the stability of the interaction of these compounds with the allosteric site. As shown in Table 2, dichamanetin interacts more strongly with IDO1 than with isochamanetin (−7.8 kcal/mol and −7.5 kcal/mol, respectively), with the first compound forming two pi-alkyl interactions with Pro182. It is also noted that the participation of a weak hydrogen bond interaction with Gly11. Isochamanetin, on the other hand, forms only one interaction. As will be shown in the next section, referring to molecular dynamics simulations, both Ser12 and Pro182 residues contribute significantly to the energetic stabilization of the formed complex.        Rings B of both dichamanetin and isochamanetin interact via pi-pi stacking with Phe185. Ring B of chamaejasmin B interacts with Pro182 and Lys186 (pi-alkyl interaction). The methoxyl group attached to the C4' carbon of the ring B of chamaejasmin B forms an alkyl interaction with Pro182. Of all the flavonoids, the only one that interacts with IDO1 from the C ring is chamaejasmin B. This occurs through a hydrogen bond with a carbonyl group on carbon 4 of the C ring. It is worthy of note the fact that steppogenin is the only compound that shows unfavorable bump and acceptor-acceptor interactions.
Pro182 interacts through pi interactions or with the A ring of flavanones or with phenolic groups attached to the same ring. It is known that the proline ring faces are partially positively charged [34]. As shown in Figure 8, hydrogens adjacent to the carbonyl and nitrogen of the amide (Hα and Hδ, respectively) are the most partially positive. The proline side chain is also conformationally restricted, allowing interaction with aromatic residues with minimal entropic or steric penalty [34].
alkyl interaction with Pro182. Of all the flavonoids, from the C ring is chamaejasmin B. This occurs throu group on carbon 4 of the C ring. It is worthy of note compound that shows unfavorable bump and accep Pro182 interacts through pi interactions or with nolic groups attached to the same ring. It is known t positively charged [34]. As shown in Figure 8, hyd nitrogen of the amide (Hα and Hδ, respectively) are line side chain is also conformationally restricted, al idues with minimal entropic or steric penalty [34]. Phe185, an aromatic residue (Figure 9), interacts with flavonoids through pi-pi and pi-pi/T stacking, and van der Waals interactions [35], which are very important for the stabilization of the structural energy of protein complexes with flavonoids.
compound that shows unfavorable bump and accep Pro182 interacts through pi interactions or with nolic groups attached to the same ring. It is known t positively charged [34]. As shown in Figure 8, hyd nitrogen of the amide (Hα and Hδ, respectively) ar line side chain is also conformationally restricted, a idues with minimal entropic or steric penalty [34]. Comparing our docking modeling of Kushenol Comparing our docking modeling of Kushenol E (Figure 7) with the docking modeling proposed by Kwon [24] (Figure 10), we can see that our simulation shows this compound in a different orientation in the cavity of the allosteric site, with a deviation of the longitudinal axis towards the B direction for the prenyl group attached to the C6 carbon of the A ring. According to Kwon, the B ring of the flavonoid is positioned inside the cavity formed by the amino acids Leu178, Ile181, and Phe306 ( Figure 10). Differently, we found the same ring located next to residues Leu178, Ile181, and Phe306. Furthermore, according to Kwon, A and C rings interact with Ser309 and Asn313. Instead, in our model, this interaction takes place with the Pro182, and there is assistance from a few van der Waals interactions. The Kwon model shows the prenyl group attached to carbon 6 of the ring interacting with Pro314, whereas our model shows interaction with Pro182. Since the structure of the IDO1-Kushenol E complex has not been experimentally determined to date, a definitive answer about the correct orientation of Kushenol E will only be inferred after the structural elucidation study by X-ray crystallography, nuclear magnetic resonance, or cryo-electron microscopy. after the structural elucidation study by X-ray crystallography, nuclear magnetic resonance, or cryo-electron microscopy.

Molecular Dynamics Simulations and Free Energy Calculation
To further refine and investigate the ability of molecular docking simulation to predict the conformational profile of the IDO1-flavonoid complexes, and to calculate the respective binding free energies, molecular dynamics simulations were carried out.
Three-dimensional structures of the IDO1-flavonoid complexes at 20 ns, 40 ns, 60 ns, 80 ns, and 100 ns are shown in Figure 11. With the exception of chamaejasmin B, which detaches from the allosteric site cavity after 40 ns, all other flavonoids remained docked. Interestingly, after detaching from the allosteric site, chamaejasmin B stabilizes next to the JK loop. Furthermore, only one of the isosacuranetin subunits interacts with the allosteric site, while the other subunit is exposed to the solvent. Moreover, Figure 11 also shows that after 80 ns, steppogenin is almost completely detached from the allosteric site.

Molecular Dynamics Simulations and Free Energy Calculation
To further refine and investigate the ability of molecular docking simulation to predict the conformational profile of the IDO1-flavonoid complexes, and to calculate the respective binding free energies, molecular dynamics simulations were carried out.
Three-dimensional structures of the IDO1-flavonoid complexes at 20 ns, 40 ns, 60 ns, 80 ns, and 100 ns are shown in Figure 11. With the exception of chamaejasmin B, which detaches from the allosteric site cavity after 40 ns, all other flavonoids remained docked. Interestingly, after detaching from the allosteric site, chamaejasmin B stabilizes next to the JK loop. Furthermore, only one of the isosacuranetin subunits interacts with the allosteric site, while the other subunit is exposed to the solvent. Moreover, Figure 11 also shows that after 80 ns, steppogenin is almost completely detached from the allosteric site.
It is worthy of note that although the interactions of chamaejasmin B were not strong enough to keep it docked in place during 100 ns, this flavonoid showed the "best" docking energy. It is possible that its isosacuranetin subunit ( Figure 12) be a good candidate for future studies of docking and molecular dynamics. It is worthy of note that although the interactions of chamaejasmin B were not strong enough to keep it docked in place during 100 ns, this flavonoid showed the "best" docking energy. It is possible that its isosacuranetin subunit ( Figure 12) be a good candidate for future studies of docking and molecular dynamics. In order to verify how flavonoids could alter the conformational profile (packing/folding) of IDO1, radii of gyration (Rg) plots were calculated not only for IDO1 in its apo form but also for each IDO1-flavonoid complex. The corresponding plots are shown in Supplementary Material (Figure S1). In molecular dynamics, the radius of gyration is a measure of the stability of a system consisting of several particles and describes the   In order to verify how flavonoids could alter the ing/folding) of IDO1, radii of gyration (Rg) plots were cal apo form but also for each IDO1-flavonoid complex. The In order to verify how flavonoids could alter the conformational profile (packing/ folding) of IDO1, radii of gyration (Rg) plots were calculated not only for IDO1 in its apo form but also for each IDO1-flavonoid complex. The corresponding plots are shown in Supplementary Material ( Figure S1). In molecular dynamics, the radius of gyration is a measure of the stability of a system consisting of several particles and describes the variation in density as a function of the center of mass over time. According to Figure S1, the IDO1-steppogenin complex is the one with the highest degree of energetic-structural destabilization when compared to the other complexes. The IDO1-dichamanetin complex was the most stable. These stability data are supported by the values of relative molecular docking energies, shown in Table 1, and by the values of free energy of interaction between the flavonoids and IDO1, shown in Table 3. Based on Figure S2, one can assess how each flavonoid changes the conformation of IDO1 (its apo conformation). In this context, Kushenol E and isochamanetin distort the IDO1 structure more significantly. Whereas dichamanetin is the flavonoid that less intensely distorts IDO1.
In addition to the radius of gyration and RMSD plots, root means square fluctuation (RMSF) for Cα of each enzyme residue was calculated ( Figure S3). This measure is described by means of a two-dimensional graph where the abscissa shows each Cα associated with a given residue and the ordinate shows the average value of the fluctuations. It is verified that helices B (residues Cys159 to Met190) and G (residues Pro300 to Asn313), important for the formation of the allosteric site, do not undergo significant fluctuations due to the interaction of flavonoids. The EF loop (residues Gly278 to His287), significatively fluctuates due to the interaction of Kushenol E, steppogenin, and isochamanetin. Such an effect is less significant when dichamanetin is the interacting flavonoid. As mentioned before, residues between 251 and 300 and the EF loop are part of the molecular oxygen access tunnel structure. Figure S3 indicates that, of the three flavonoids, only dichamanetin does not "stiffen" the tunnel. The JK loop did not fluctuate significantly under the action of flavonoids.
RMSF plots were also constructed for each docked flavonoid into the allosteric site of IDO1 along with the molecular dynamic simulations ( Figure S4). As can be seen, dichamanetin and isochamanetin showed lower atomic fluctuations than Kushenol E and steppogenin do. These results reinforce, once again, the relative stability of the interactions of each flavonoid with the allosteric site.
The complexes were also analyzed for the frequency of hydrogen bonds formed during the molecular dynamics simulation ( Figure S5). Dichamanetin has a higher frequency of hydrogen bonds than isochamanetin and Kushenol E. These data demonstrate that hydrogen bonds are important for the stabilization of the IDO1-flavonoid complexes, mainly for dichamanetin and isochamanetin. However, steppogenin shows an "anomalous" hydrogen bond frequency. Perhaps that should be due to the highest degree of freedom ("pose fluctuations"), since, as shown in Figure 11, that compound seems to be slightly bound to the allosteric site.
The MM-PBSA method evaluates conformational fluctuations as a function of the simulation trajectory and, consequently, it is a useful approach for inferring the entropic profile of the formed complexes, as well as for calculating the free energy of binding of each flavonoid with the IDO1. The results are shown, respectively, in Tables 3 and 4. According to Table 3, the affinity of each flavonoid with IDO1 increases in the following order: steppogenin (−5.73 kcal/mol); Kushenol E (−21.65 kcal/mol); isochamanetin (−22.09 kcal/mol); dichamanetin (−25.93 kcal/mol). It is worth noting the calculated affin-ity value for steppogenin (the negative control), testifying once again, and in accordance with the experimental data [24], that steppogenin forms the least stable complex with IDO1. The free energy contributions of each constituent residue of the allosteric site are shown in Table 4.

IDO1
The crystallographic structure of IDO1 in its holo form, with a resolution of 2.44 Å (PDB ID: 7A62) [15], was obtained from the Protein Data Bank [36]. The enzyme preparation steps were carried out as follows: (1) non-essential water molecules were removed; (2) polar hydrogens were added to the enzyme; (3) partial charges were calculated using both Kollman and Gasteiger's approaches [37,38].

Molecular Docking
Molecular docking simulations were performed with the AutoDock Vina 1.1.2 program [39]. This tool uses a hybrid scoring function, i.e., it is a combination of empirical and knowledge-based functions [39]. The binding energy is predicted as the sum of distancedependent atom pair interactions as follows: Here d is the surface distance calculated as follows, where r is the interatomic distance and R i and R j are the radii of the atoms in the pairs: Every atom pair interacts through a steric interaction given by the first three terms of the equation shown below. Moreover, depending on the atom type, there could be hydrophobic and non-directional H-bonding interactions, given by the last two terms. The combination of an attractive Gaussian function with a repulsive parabolic function reproduces the general shape of a typical Lennard-Jones interaction, provided the Gaussian term is negative and the parabolic positive.

Molecular Dynamics Simulations and Free Energy Calculation
All molecular dynamics simulations were performed with the CHARMM36 force field [40] implemented in the GROMACS 2018.1 program [41]. The parameterization of flavonoids was performed using the CHARMM General Force Field (CGenFF) program [42]. Subsequently, a dodecahedral simulation box, including solvent inside (TIP3P water model) was created. Sodium and chlorine ions were added to the system, aiming to electronically neutralize it. Periodic boundary conditions were used.
Each complex was submitted to energy minimization using the steepest descent algorithm, involving 50,000 steps and the convergence criterion of less than 2.39 kcal/mol. Subsequently, two 100 ps molecular dynamics simulations were performed, aiming to balance the IDO1-flavonoid complexes. In the first simulation, the NVT ensemble was used, and in the second one, the NPT ensemble was used. In both cases, the simulation temperature was kept constant at 300 K and, when performing the NPT ensemble, the pressure was kept constant at 1 bar.
After the equilibration of the complexes, 100 ns molecular dynamics simulations were performed, aiming to calculate the free energies of interaction of each flavonoid with IDO1 (production stage). The running conditions were the following: ensemble NPT, where temperature was maintained, using the V-rescale implementation of Berendsen's thermostat [43]. A molecular frame was sampled every 10 ps. To keep the pressure constant, the Parrinello-Rahman pressure coupling method was used [44]; for the longrange treatment, the PME method was used [45].
Free energy calculations were calculated with the g_mmpbsa tool [46]. The binding energy can be individually decomposed as a function of the constituent residues of the macromolecular target. Initially, the E MM , G polar, and G nonpolar energy components of the individual atoms of each residue are calculated in the bound, as well as the unbound form. Subsequently, the contribution to the interaction energy ∆R BE x of the residue x is calculated as follows: ∆R BE where, A bound i and A f ree i are the energy of atom i of residue x in bonded and unbonded states, respectively, and n is the total number of atoms in the residue. The energy contribution added to all residues is equal to the interaction energy, that is, ∆G binding = ∑ m x=0 ∆R BE x , where m is the total number of residues that make up the ligand-protein complex [46].

Conclusions
IDO1, one of the enzymes that participate in the immune modulation process, is an interesting macromolecular target for anticancer action. From a known drug (Kushenol E), with proven anticancer activity, extensive virtual screening studies were carried out in digital databases. Both ligand-and structure-based approaches were performed. In this context, three natural products, isolated from plants, were identified as potential allosteric site binders of IDO1. Two of which, dichamanetin and isochamanetin, were shown to be promising IDO1 inhibitors and, in our opinion, should be considered in future studies of biological activity, molecular optimization, organic synthesis, and enzymatic assays. From the point of view of structural biology (conformational characteristics of IDO1), it was found that helices constituting the allosteric site did not undergo significant modifications; the same happened for the JK loop (located at the entrance of the catalytic site). In order to investigate, if the presence of the identified allosteric binders could weaken the binding of orthosteric binders (O 2 , tryptophan), as well as weaken the binding of the heme cofactor, further simulations, involving multiscale quantum mechanics/molecular mechanics (QM/MM) are being performed in the laboratory.