Insights on the mechanism of action of immunostimulants in relation to their pharmacological potency. The effects of imidazoquinolines on TLR8

Imidazoquinolines are powerful immunostimulants (IMMS) that function through Toll-like receptors, particularly TLR7 and TLR8. In addition to enhancing the immune response, IMMS also function as antineoplastic drugs and vaccine adjuvants. These small compounds display almost the same molecular structure, except in some cases in which atom in position 1 varies and changes the imidazole characteristics. A variable acyclic side chain is also always attached at atom in position 2, while another chain may be attached at atom in position 1. These structural differences alter immune responses, such as the production of interferon regulatory factor and nuclear factor-κB (IRF-NFκB). In this work, quantum mechanics theory and computational chemistry methods were applied to study the physicochemical properties of the crystal binding site of TLR8 complexed with the following six IMMS molecules: Hybrid-2, XG1-236, DS802, CL075, CL097 and R848 (resiquimod). The PDB IDs of the crystals were: 4R6A, 4QC0, 4QBZ, 3W3K, 3W3J, and 3W3N respectively. Thus, were calculated, the total energy, solvation energy, interaction energy (instead of free energy) of the system and interaction energy of the polar region of the IMMS. Additionally, the dipole moment, electrostatic potential, polar surface, atomic charges, hydrogen bonds, and polar and hydrophobic interactions, among others, were assessed. Together, these properties revealed important differences among the six TLR8-immunostimulant complexes, reflected as different interaction energies and therefore different electrostatic environments and binding energies. Remarkably, the interaction energy of a defined polar region composed of the highly polarized N3, N5 atoms and the N11 amino group, acted as a polar pharmacophore that correlates directly with the reported immunopharmacological potency of the six complexed molecules. Based on these results, it was concluded that accurate physicochemical analysis of the crystal binding site could reveal the binding energy (measured as interaction energy) and associated molecular mechanism of action between IMMS and TLR8. These findings may facilitate the development and design of improved small molecules with IMMS properties that are targeted to the TLR system and have enhanced pharmacological effectiveness and reduced toxicity.


Introduction
The discovery of interferon by Isaacs and Linderman [1] was a milestone in immunology. This cytokine is produced and released in response to the harmful effect of viruses, bacteria, and parasites, as well as other strange or invasive cells. However, interferon was not purified and thus available for therapy and research until 20 years later [2]. Consequently, those findings prompted the search for modulators of the immune system that are capable of inducing interferon synthesis. Interestingly, nucleosides appeared to be good candidates to activate interferon production and, accordingly, the immune system. Notably, an early work utilizing pyrimidine derivatives indicated a clear antibacterial activity of, particularly, thymine and purine compounds [3].
Several antiviral compounds derived from nucleosides have been developed to activate the immune response. Among the original modified nucleosides that initiated modern immunopharmacology were the following: the nucleoside iododeoxyuridine [4], the tri-fluorinated uracil nucleoside Trifluridine [5] and the iodinated compound 5-iodo-2-deoxyuridine named Vidarabine [6]. The guanosine derivative hydroxyethoxymethyl-guanine, (acyclovir) synthetized by Elion et al. in 1977 [7] may be the best example of a successful antiviral agent based on a purine.
Based on the adenine molecule, imidazoquinoline (IMZQ) compounds were developed also with the aim of fighting viral infections. In addition to the adenine molecule, the structure of IMZQ includes a third ring to yield the quinoline moiety joined to the five-membered aromatic heterocyclic imidazole ring. Among the first synthetized IMZQ was the compound S-25059, followed by imiquimod (S-26308, R837) in 1983 [8]. The analog resiquimod (S-28463, R848) was later synthetized and released, as well as antiviral compound [9]. The induction of interferon by imiquimod was demonstrated in early studies [10], and consequently both imiquimod and resiquimod have been typified as immunomodulators [11]. It is currently accepted that IMZQ is capable of modifying the immune response to induce antiviral and antitumor activity, basically via cytokine production, e.g., interferon-α induction, mediated partially by NF-κB activation, which in turn stimulates the innate immune response and acquired immunity [8,11].
However, an important issue is the activation of the immune response by Toll like receptors (TLRs) induced by IMZQ, particularly TLR7 and TLR8, which function in immune responses against viral infections [12]. Thus, the structural basis of the mechanism of action of IMZQ and, therefore, IMMS has been unveiled in a larger capacity by valuable studies applying X-ray crystallography. In this regard, several IMMS analogs have been crystallized, showing binding particularly to TLR8 and revealing that IMMS form a bridge between two TLR molecules. Based on these findings, we selected six IMMS complexed with TLR8: Hybrid-2 [13], XG1-236, DS802 [14], CL075, CL097, and R848 [15]. The aim was to explore the binding energy calculated as interaction energy, along with binding characteristics associated with their pharmacological activity. To achieve this goal, we studied the binding core of the six IMMS-TLR8 complexes from the reported X-ray structures. These cores were studied comparatively based on their electronic structure and molecular properties. For comparison, two unbound IMMS antagonists (compounds 15 and 16) reported by Shukla et al. [16] were also examined to characterize their molecular and electronic properties.

Interaction analysis of the IMMS-TLR8 complex
The interactions between the crystal structures of Hybrid-2, XG1-236, DS802, CL075, CL097 and R848 complexed to TLR8 were analyzed. Only the first core of residues belonging to the TLR8 binding pocket was included, and the cut off distance between the IMMS and the receptor was set at 4 Å. The included amino acids ranged from 9 to 12, among all TLR8-IMMS complexes. Residues showing hydrophilic (polar) interactions achieved a distance of approximately 3.3 Å, while those with hydrophobic interactions reached 4 Å with respect to the IMMS atoms. The Ligand Explorer program was used for this assessment. The geometrical parameters (Å) were set up as follows: hydrogen bonds: 3.3, water bridged hydrogen bonds: 3.3, and hydrophobic interactions: 4.0 [20]. Thus, important binding interactions were studied, namely: hydrogen bonds, water bridged hydrogen bonds, hydrophobic interactions (van der Waals) and electrostatic charges of atoms of the agonists interacting with those atoms in the receptor amino acids. The number and type of molecular interactions of the TLR8-agonist complex and the calculation of the interaction energy of each complex permitted the evaluation of the energy of the receptor-ligand interaction, which was considered an index of the extent of binding.
Two interaction energy values were calculated; the total interaction energies included the IMMS molecule, all adjacent binding pocket residues and molecules of water (hydrophilic and hydrophobic interactions) and the selected polar region interaction energy, that included N3, N5 and the N11 amino group as well as the C4 atom of the IMMS compound (this region was defined as: the polar pharmacophore), whereas from the TLR side, were included two residues Asp543 and Thr574, as well as, a fixed molecule of water. A simplified procedure to calculate the interaction energy is provided by the following formula: where IE is the interaction energy of binding, RLW is the energy of the complex, formed by the receptor amino acids, the ligand and water energies, R is the energy of the amino acids of TLR8, L is the energy of the ligand agonist, and W is the energy of molecules of water.
The binding interaction energies were determined by single point calculations using the DFT with the ωB97X-D function and the 6-31G Ã basis set level. Likewise, the volume of the binding pocket cavities of each TLR8 dimer was calculated using the Swiss-Pdb Viewer program [21]. All ab initio calculations were conducted using Spartan'14 software [22]. The computer consisted of an eight-core AMD processor at 4.0 GHz. A line regression and correlation statistical analysis, comparing the natural logarithm values of TLR8 agonist activity and pharmacophore interaction energy values was calculated with Origin 8 Pro (OriginLab, Northampton, MA, USA).
The atom at position 1 may be either C (XG1-236), O (DS802), or S (CL075). Likewise, an amino group (CL097) or a 2-methylpropan-2-ol side chain may be found at position 1 (Hybrid-2 and R848). At position 2, diverse side chains are found, namely butyl (I and III), N2-butyl (II), propyl (IV) and ethoxy methyl (V and VI). The antagonist molecules 15 and 16 have a nitrogen atom at position 1, an ethylaminomethyl side chain at position 2, and a 2-methylpropan-2-ol side chain at N3. Clearly, the structure of the antagonist modifies importantly the polar pharmacophore region because of the steric bulk of the N3 side chain in both antagonist molecules and the absence of the N11 amino group in compound 16. The potency of the IMMS agonists as well as the antagonistic IMMS seem to be related to the chemical structure and associated physicochemical properties.

Immunostimulant potency of IRF-NFκB transcription factors
The studied IMMS show different extents of IRF-NFκB transcription factor stimulation of TLR7 and TLR8. Thus, the EC 50 data show different pharmacological potency for TLR7 and TLR8 in the same IMMS. However, the pharmacological potency of the agonists Hybrid-2 [13] and XG1-236 [14] acts first and second in both TLR7 and TLR8, respectively.

IMMS crystal structures
The six IMMS agonists isolated from TLR8 crystals have two monomers, namely, A and B. The pairs of monomers for each agonist were compared by measuring the root mean square deviation (RMSD) of the atomic positions. The largest RMSD difference was found in the CL097-A and B pair due to the 2-ethoxymethyl side chain, which adopts two opposite conformations that interlace when superimposed. Fig 2 shows the IMMS superimposition and RMSD values.

Assessment of the physicochemical properties of IMMS
The molecular and physicochemical properties of both A and B crystal monomers were calculated by single point ab initio calculations. The monomers with the highest energy were discarded, leaving for further analysis only the monomers with the lowest energy. The physicochemical values of both agonists and antagonists show some differences due to the varied substitution of atoms and inherent side chains. The properties hardness (η) and electronegativity (χ) were calculated as follows: η = (εLUMO − εHOMO) / 2 and χ = − (εLUMO + εHOMO) / 2 respectively. Both antagonists were studied; however, the crystal of antagonist compound 15 was not available. Therefore, compound 15 was modeled using the crystal of compound 16 as a model and adding an amino group at position 11. Table 2 shows the values for the electronic and physicochemical properties of the lowest energy IMMS conformers from single point ab initio calculations. The molecular properties of the antagonists are also included.
The three agonists with the lower total energy were compounds I, IV and VI because of the sulfur atom in compound IV and the large number of atoms in I, IV and VI. In addition to the lowest energy, CL075 (IV) also has the lowest area, volume, polarizability, polar surface, hardness, HOMO and LUMO, and the highest solvation energy and electronegativity. The

Fig 2. RMSD values and superimposed views of monomers A and B of the IMMS agonist. Superimposition of compound CL097
showing the interlaced 2-ethoxymethyl side chain that markedly increases the RMSD value. The main differences are due to the side chains because the five-sided ring system alone shows significantly lower values.
https://doi.org/10.1371/journal.pone.0178846.g002 molecular properties followed roughly three trends; one trend includes total energy, HOMO and LUMO excluding compound I. A second trend includes most of the following properties: area, volume, polar surface area, DM, polarizability and, to some extent, hardness. Finally, the third trend is consistent with the solvation energy and electronegativity.
Interestingly, the compound 16 antagonist lacking the N11 amino group is more potent than its related compound 15. Coincidentally, compound 16 has 7 of the 11 physical and electronic properties with lower values: area, volume, polar surface, HOMO, LUMO, hardness and polarizability. Because of its large number of atoms, compound 15 has higher total and solvation energies, resulting in a larger DM and increased electronegativity.  show an intense negative EP and a marked positive EP in the N11 amino group. In contrast, the EP of the antagonists exhibited a different pattern, particularly compound 16, which lacks the N11 amino group. Indeed, the five-membered EP ring is the most versatile because of its atomic variety at position 1 and marked differences in the attached side chains. Molecules I and VI differ only in one carbon atom, in which an oxygen is interchanged in the side chain attached to atom at position 2. This feature produces an extra negative EP in VI and somewhat decreases the magnitude of the dipole moment and direction ranking for the second value as compared with compound I, which has the largest DM among the six IMMS. Additionally, molecules I and VI present the hydroxyl group of the side chain attached to atom at position 1 are in opposite position. These differences, however, are sufficient to situate these compounds as the most and less potent, respectively. Molecules II to V have a very similar EP pattern with a variable DM value but a similar direction. The EP of the antagonists are larger than those of the agonists and are also quite different in the polar pharmacophore region, also showing a yellow region in the benzen ring of the quinoline moiety, which indicates a tendency toward a negative charge. Notably, the antagonists have larger DM values and point opposite to agonist VI.

Analysis of the crystal IMMS-TLR8 Interaction
The crystals of the six IMMS-TLR8 complexes show similar binding patterns that include both A and B monomers of the receptor, forming a dimer bridged by two molecules of the ligand. Notably, one TLR8 monomer binds the ligand with polar residues, and the second monomer makes mainly hydrophobic contacts with the quinoline ring system and the side chains. For the second ligand molecule, the inverse occurs; therefore, the atomic arrangement and the residues involved are similar in both ligands. Fig 5 shows a typical example of the TLR8 dimer, including the region between the two molecules of CL075 (PDB ID: 3W3K).

Binding pocket area and volume
Each of the two IMMS bound to the TLR8 dimer is enclosed in a cavity that forms the binding pocket. Thus, the volume of each binding pocket of the studied TLR8 agonists was calculated. Some cavities showed important differences among the six IMMS, particularly compounds IV and VI. However, the volume of the IMMS was quite similar among them. However, as expected, compounds I and VI had the highest volume because of its additional side chain. Two examples of the shape and position of the binding pocket for each TLR8 dimers are shown in Fig 6, and values of the binding pocket volume and IMMS physical properties are shown in Table 3.

The inner shell arrangement of the IMMS-TLR8 complex
The shell of the interacting residues in the binding pocket shared by the two monomers of the TLR8-IMMS complex at a distance of 4 Å, shows seven residues in common in the six studied agonists. Thus, in compounds I, II and IV four of them are in monomer A (Tyr348, Tyr353, Val378 and Phe405) while the remaining three (Asp543, Asp545 and Thr574) are in monomer B. In compounds III, V and VI the opposite occurs. Additionally, seven residues are found, to various extents, among the binding pockets bound to the IMMS: Phe346; Glycines 351, 376, and 572; Ser352; Arg429 and Val573. Thus, compounds II, IV and VI are bound to nine, compounds I and III to 11 and compound V, to 12 residues, respectively. Table 4 shows the distribution of residues by monomer for the studied TLR8 dimers. The two TLR8 monomers are colored differently to show their relationship with the two CL075 ligands in between. The insets show a magnified view that includes the two ligand molecules bound to residues Asp543 and Thr574. The (blue) atoms N3 and N5 and the N11 amino group that make up the polar pharmacophore form in all cases hydrogen bonds with the oxygen atoms (red) of the two residues, (for atom numbering see Fig 1). The ligands clearly have an inverted position with respect to the two TLR8 monomers. Marked angles show the rotation differences between the ligands.   In all the IMMS-TLR8 complexes, the two polar amino acids Asp543 and Thr574 interact with N3, N5 and N11 atoms of the polar pharmacophore. It is important to recall that one TLR8 monomer interacts with its polar residues binding the ligand polar pharmacophore, and the second TLR8 monomer interacts with the same ligand, mostly with hydrophobic residues through van der Waals forces. The second ligand interacts in the opposite manner, similarly forming a bridge between both TLR8 monomers. Fig 7 shows the interaction at a distance of 4 Å of three typical IMMS-TLR8 complexes. The residues are colored to be distinguished by their respective monomers.

Electrostatic potential of the IMMS-TLR8 complex
The EP maps of the IMMS-TLR8 complex at a distance of 4 Å show important differences in the polar region, mainly due to the different number of fixed water molecules forming the hydrogen bonds network. Two EP patterns were observed; one pattern occurs in IMMS I, III and V, with two fixed water molecules and 11, 11 and 12 residues, respectively. A second pattern appears in IMMS II, IV and VI, with one, four and two fixed waters and nine residues in each. The first pattern includes Val573. The variety of atoms at position 1 somewhat modifies the EP pattern in neighboring atoms of the imidazole ring.
Similarities are observed for the different side chains of IMMS. Additionally, the hydrophobic interactions show important variations because of differences in the number of interacting residues. However, the most hydrophobic region (the quinolone-benzene ring) showed low variation in all cases. Notably, a neat and well defined highly polar IMMS pharmacophore region was identified, which included N3, N5 and the N11 amino group and the C4 atom. In that region, Asp545 and Thr574 form in all the studied IMMS-TLR8 dimer complexes a remarkable hydrogen bonds pattern with one fixed water molecule.

Hydrophobic interactions of IMMS in the TLR8 binding pocket
The IMMS-TLR8 complex has seven to ten hydrophobic residues, but only five of them are common to all IMMS. However, the number of hydrophobic contacts of the same residue may vary depending on the interacting IMMS. Aromatic residues make abundant and different contacts, particularly with the imidazole benzene ring. Thus, the number of van der Waals interactions varies from 32 to 37.  IMMS with their respective interacting residues from the TLR8-dimer binding pocket. Additionally, Table 5 shows the interacting residues and number of hydrophobic interactions.

Atomic charges of the IMMS-TLR8 complex
Atomic charges were calculated under unbound and bound conditions for both IMMS and the corresponding atoms of the binding pocket residues. In the unbound condition, the atomic charges of the N atoms 3, 5 and the N11 amino group, which form the polar pharmacophore, are highly negative, especially the N11 atom. In contrast, the two H atoms of the N11 amino group and the adjacent C4 atom are highly positive. As expected, differences in atomic polarization are observed among the six IMMS. The atoms at position 1 exhibit diverse charge values because of the variety of atoms at this position; the same occurs with its contiguous number 2, 9a and 10 atoms. The remaining atoms of the quinoline moiety have either negative (C6 to C9)  or positive charges (C3a and C5a). The binding of the IMMS atoms with the TLR8 residues polarizes the quinoline benzene ring atoms as well as the polar pharmacophore atoms, excluding the N11-H atoms, which remain almost unchanged.
In contrast, the interacting carboxylic oxygen atoms (O-δ1 and O-δ2) of Asp543 show different extents of polarization that favors the atom interacting with the N5 atom of the quinoline group. However, the atomic charges of NH, O, and O-δ1 of Thr574 remain almost unchanged. The ring atoms of Phe405 were the most interactive with the quinoline-benzene ring, showing a crowded hydrophobic ring-to-ring interaction. Fig 10 shows the curve tendencies for the groups of atoms in IMMS, and the atoms of the TLR8 residues under both unbound and bound conditions.

Interaction energy assessment of the IMMS-TLR8 complex
The total interaction energy for each of the six IMMS-TLR8 complexes was calculated, as well as the pharmacophore (polar region) interaction energy. The interaction energy of the IMMS-pharmacophore-TLR8 polar residues was calculated considering the following as a system: C4, N3, N5 and the N11 amino group of IMMS as the pharmacophore and the concurrent Thr574 and Asp543 residues of TLR8, with one fixed water molecule (see Fig 8).
Notably, the tendency of the polar pharmacophore interaction energy curve was opposite to that of the total interaction energy. Interestingly, the interaction energy curve elicited by the six complexed polar pharmacophores largely resembles the curve formed by the IMMS experimental pharmacological potency. In fact, the line regression and correlation analysis for both curves showed a linear trend with statistically significant values of r = 0.838 and p = 0.037 (for the pharmacological potency values was used their natural logarithm). Fig 11 shows the two interaction energy curves calculated from the binding pocket of the IMMS-TLR8 complex and the polar pharmacophore-TLR8 -Thr574 and Asp543 interactions. In addition, the potency of the IMMS agonist activity stimulation curve toward IRF-NFκB transcription factors is included for correlation (for values, see Table 1). In panel Mechanism of action of immidazoquinolines on TLR8 D is shown the line regression of TLR8 agonist activity (ln values) and pharmacophore interaction energy. The associated data for both interaction energy calculations are shown in Table 6.

Discussion
The chemical structure of the studied IMMS is quite similar, showing only few differences in the side chains attached to atoms at position 1 and 2 of the imidazol ring. Likewise, atom at position 1 also differs in some cases, showing the following trends for pharmacological effectivity as follows: C<O<S<N-H. The side chain at position 2 is always present and seems to be quite important in compounds I and VI, differing by one atom but exhibiting the highest and lowest pharmacological potency, respectively [13,24]. Notably, the highest pharmacological potency of compound I correlates quite well with its lower total and polar pharmacophore interaction energy values. A comparison of the binding interactions between compounds I and VI shows that two more residues participate in the binding pocket, two bridge the water bonds and two more atomic interactions favor compound I. These subtle differences account for better interaction energy and increase the pharmacological effectiveness of compound I. Similarly, compounds II and III have different chemical structure only in atom at position 1 with an analogous side chain, which changes their interaction energy and pharmacological potency. Thus, varying the number of fixed water molecules, bridged hydrogen bonds, and number of van der Waals interactions account for the pharmacological effectiveness. Likewise, the low effectiveness of compound V may be due to the ethoxy methyl side chain, which is similar to that of compound VI [24]. Interestingly some physicochemical properties correlated to a large extent; for example, the total energy pattern of the IMMS molecules resembles the total interaction energy curve. Both   curves exhibit lower values for compounds I, IV and VI. Additionally, the interaction energy curve emphasize that the binding pocket environment differs for each IMMS-TLR8 complex. Likewise, both the solvation energy and electronegativity curves mimic the polar pharmacophore interaction energy curve, suggesting that the properties are directly associated with that pharmacophoric activity. In contrast, the frontier orbital (HOMO and LUMO) energy showed an inverse tendency compared with the polar pharmacophore energy. The polar surface area correlates quite well with the polar pharmacophore area excluding compound CL075 (IV), which has a sulfur atom at position 1 [15]. In fact, because of its physicochemical properties, the S atom modifies the overall properties of compound IV. Notably, the atoms of the polar pharmacophore were the most polarized: N3, N5 and N11, followed by C5a and C4 to a lesser extent, reinforcing the importance of the polar pharmacophore. Regarding the residues in the binding pocket, only a few carbon atoms of Phe405 and one oxygen atom of Asp843 were slightly polarized. Nevertheless, those residues appear to be essential to activate NF-κB [15].
According to our analysis of the IMMS mechanism of action, the interaction energy (binding) of the polar pharmacophore is the main issue related to the pharmacological potency in the six studied IMMS-TLR8 crystals. The polar pharmacophore region is somewhat different from the four region model proposed by Musmuca et al. [26] for the IMMS molecule and related immunomodulators. In an interesting 3-D QSAR analysis of 156 interferon-inducing agents, those authors described a pharmacophoric model with four regions, including the adenine moiety with two regions: the HA region and a polarized area in which those authors suggested the presence of hydrogen bonds and polar interactions. Both regions are concurrent with the N3, N5 and N11 amino group of the IMMS, which we are denoting as the polar pharmacophore, that also includes the hydrogen bonds network with a fixed water molecule and the two polar residues of the TLR8 binding pocket. The other two pharmacophoric regions proposed by Musmuca et al. [26] are a fillable steric pocket and a hydrophobic area.
The two studied IMMS antagonists, compounds 15 and 16 [16], have important molecular modifications, particularly in the polar pharmacophoric region. Thus, both molecules contain the N3 atom covalently bound to a large 2-methyl-propan-2-ol side chain. Additionally, compound 16, the most effective antagonist, lacks the pharmacophoric N11 amino group. These molecular changes, which block or delete the polar region, clearly show that this region is essential to activate the immune response through the IMMS-TLR complex. Therefore, the molecular mechanism of the antagonist is likely due to the absence of the fixed water molecule, the absence of the hydrogen bonds network and null or a low contribution of residues Asp545 and Thr574. However, both antagonist compounds are able to occupy the binding pocket and interact, likely by increased van der Waals forces by their two side chains. Additionally, the two antagonists have the DM vector directed toward the hydrophobic region; the compound 15 vector is slightly shorter, and both possess a high solvation energy, implying low solvation. Those features apparently favor hydrophobicity and, likely because of a distinct electronic environment, results in the antagonist effect. In contrast, compounds I and VI have the DM vector toward the polar region, decreasing the solvation energy and, therefore, increasing the solvation of the molecule. In this regard, compound I has the lowest solvation energy and the largest DM pointing toward N5. Likewise, a comparison of compounds II and III shows that compound II has the lowest solvation energy and a higher DM, whereas compound III has the shortest DM. All these data correlate quite well with the pharmacological effectiveness of these compounds.
Analysis of the pharmacological potency of IMMS reveals almost the same trend for the activation of both TLR7 and TLR8. However, compounds I, V and VI are clearly more capable of stimulating TLR7 at lower doses. In contrast, TLR8 responds better to compounds II, III and IV. This difference occurs regardless of the close phylogenetic relationship of TLR7 and TLR8 to the same TLR subfamily [27], showing a high level of sequence homology [28,29]. In fact, a theoretical model predicts different residues for the binding site of TLR7 with respect to TLR8 [30]. Thus, the differences in IMMS potency may be related to the intrinsic structural disparities between the TLR7 and TLR8 binding pocket, which yield functional differences for both receptors, as previously proposed by Gorden et al. [23]. Additionally, because of their different hydrophobic side chains, IMMS may possess TLR selectivity similar to that reported for R837 (imiquimod), which specifically activates TLR7 [15].
The small ligands IMMS activate both TLR7 and TLR8 similarly to the natural pathogens to yield a cascade of biochemical events that initially stimulate several type of cells, such as monocytes and NK/NKT, T, B, mast and tumor cells. Most of these cells in turn undergo migration, proliferation and apoptosis, along with the production of considerable amounts of interferons, interleukins and tumor necrosis factors, among others [31]. The biochemical reactions of the immune system elicited by IMMS have been confirmed to be effective against viral infections [32]. In fact, TLR7 and TLR8 sense single-stranded viral RNA [33]. The pleiotropic characteristics of the IMMS that trigger the immune response also play an important role in improving cancer therapy, regarding this, see for instance [31] and the important review of Schön and Schön [34]. Besides, IMMS show an angiogenesis inhibitor effect [35] and likely act as vaccine adjuvants [36]. According to our data, all of these pharmacological effects may be correlated, to a large extent, to the regional interaction energy of IMMS. We further validate the usefulness of the crystal analysis as a tool to unveil new data that can be used to improve the design of the IMMS molecule to produce better immunostimulant responses by lowering toxicity and enhancing effectiveness.

Conclusions
The present results show that the application of quantum mechanics theory and computational chemistry methods to analyze the crystal structure of TLR8 interacting with different small synthetic IMMS, focusing particularly on the binding site of the receptor and the IMMS atomic interactions, enabled elucidation of the molecular properties and calculation of the interaction energy of the IMMS-TLR8 complex, which correlates quite well with the immunopharmacological activity.
The IMMS molecules clearly bears a defined polar region that acts as a regional pharmacophore, which seems to have variable selectivity for TLR7 and TLR8. This finding might facilitate the rational design and development of new agonist and antagonist IMMS that function specifically or selectively in the TLR system or toward other receptors capable of detecting pathogen-associated molecular patterns (PAMP), improving the therapeutic modulation of natural immunity and therapeutic efficiency toward infectious and neoplastic diseases.