Expected signs of contributions to ΔH and ΔS.
With more than 4000 compounds identified up to now, flavonoids are present in human diet since they can be found in fruits, vegetables, seeds, grains, and beverages, such as wine and tea. Over the years, medicinal properties of these polyphenolic compounds have been noticed. Consequently, the search for the biological targets and for the description of flavonoids action mechanism has been growing. Fluorescence spectroscopy and molecular docking are techniques based on physical theories which have been helping researchers to describe the interaction between flavonoids and biological targets. In this way, this chapter comes not only as an attempt to gather some works dedicated to explain flavonoid molecular mechanisms of action but also to introduce a brief theory of steady-state fluorescence spectroscopy and molecular docking.
- fluorescence spectroscopy
- molecular docking
- molecular mechanism
- physical pharmacy
More common than you might think, flavonoids are present in human diet since they can be found in fruits, vegetables, seeds, grains, and beverages, such as wine and tea . They are most famous to beautify fruits and vegetables with vivid colors, but flavonoids most powerful actions are still unknown for most of population, which is unaware the antioxidant [2, 3], anticarcinogenic [4, 5], anti-inflammatory [6, 7], antiviral [8, 9], and antimicrobial [10, 11] effects provoked by these potent compounds. But, what are those compounds, after all?
Flavonoids, with over 4000 compounds identified until now , comprehend a wide group of molecules synthesized by plants as secondary metabolites responsible for ensuring vascular plant colonization and surviving on earth’s environment . These molecules play crucial roles to plants’ life such as protection against insect attack and microbes invasion [14, 15] and by absorption of harmful ultraviolet radiation, attraction of insect pollinators by colorful anthocyanins synthesis , antioxidant action by inhibiting the generation of reactive oxygen species (ROS) , involvement in pollen germination , involvement in biological communication in the rhizosphere  and action as regulators, involved in auxin transport and catabolism .
In general, flavonoids are polyphenolic compounds with the flavan nucleus as the structure skeleton, which consists of 15 carbon atoms arranged in 3 rings (C6-C3-C6), labeled by A, B, and C (Figure 1) .
|< 0||< 0||van der Waals and hydrogen bonds|
|> 0||> 0||Hydrophobic interactions|
|≈ 0||> 0||Ionic process|
Flavonoids are classified into classes such as flavonol, flavone, flavanone, flavanol, isoflavone, and anthocyanidin, just to mention a few . The basic structure of each class is shown in Figure 2. The differences between flavonoid class structures are in the oxidation level and C-ring substitution pattern, while the difference among flavonoids belonging to the same class is in the pattern of A- and B-ring substitutions .
As reviewed by Aidyn Mouradov and German Spangenberg , flavonoid synthesis in plants occurs commonly from the aromatic amino acids phenylalanine and tyrosine. At the beginning of synthesis, the amino acid is converted into coumaroyl-CoA by a number of enzymatic reactions involving
Epidemiological studies reviewed by Romano and co-workers  have demonstrated an inverse relationship between dietary flavonoid intake and prevalence and risk of cardiovascular diseases and some types of cancer such as breast, colon, lung, prostate, and pancreas. Besides that, the authors emphasized the advances of the application of flavonoids in diseases related to the central nervous system, obesity, diabetes, inflammation, digestive system, and respiratory tract, as well as the effects of flavonoids in reproduction and the antimicrobial effects of these compounds.
Many efforts have been done to elucidate the mechanism of action of these compounds by using experimental approaches based on physical observable. Steady-state fluorescence spectroscopy is a technique widely used to follow the interaction of small molecules and biological targets. However, there are many others experimental techniques that have been used to describe biological systems, such as time-resolved fluorescence spectroscopy, nuclear magnetic resonance (NMR), isothermal titration calorimetry (ITC), circular dichroism (CD), and surface plasmon resonance (SPR).
Besides that, it has also used computational analysis in order to describe the interactions between flavonoids and molecular targets. Molecular docking and molecular dynamics are the most common computational techniques used to model biological systems.
In this way, this chapter comes not only as an attempt to gather some works dedicated to explain the molecular mechanisms of action of some flavonoids but also to introduce a brief theory of steady-state fluorescence spectroscopy and molecular docking.
2. Theory of steady-state fluorescence spectroscopy
The use of fluorescence spectroscopy experiments to understand the protein-ligand complex formation has become increasingly frequent in scientific community because many proteins have endogenous fluorescence probe such as tryptophan, tyrosine, and phenylalanine amino acids. Following the fluorescence signal of these probes, it is possible to characterize the molecular mechanism of complex formation .
As the name already suggests, this technique makes use of fluorescence as the physical observable. Fluorescence, by definition, is a photon emission mechanism which occurs by the decaying of molecule electrons from a higher energy singlet state to a lower energy singlet state with emission rate of the order of 108 s−1, resulting in a typically fluorescence lifetime close to 10 ns .
The fluorescence spectroscopy, when applied in protein-ligand system, consists in analyzing the quenching of fluorescence signal in the presence of different concentrations of ligand. Quenching may be static or collisional (dynamic) depending on the nature of interactions. In collisional quenching, the fluorophore, which is in the excited state, is deactivated returning to its ground state because of diffusive encounters with some other molecule of the solution, called quencher. On the other hand, the static quenching occurs when a fundamental and non-fluorescent complex occurs .
In a typical experiment of fluorescence spectroscopy, the way to distinguish the quenching mechanism is analyzing the Stern-Volmer plot (Eq. (1)) in different temperatures, where Fo is the fluorescence intensity in the absence of the quencher and F is the intensity in different concentrations of the quencher, [Lt] .
The increase in the constant K as a consequence of the increase in temperature is a strong indication that the quenching mechanism is collisional; in that case the constant is usually denoted by Kd. On the other hand, if the constant K decreases with the increasing in temperature, the characteristic mechanism is the static one, and the constant is usually denoted by Ksv. Figure 3 shows the different behaviors of the static and collisional quenching .
If the mechanism determined was the static one, there are couples of models that can be applied in the system in order to obtain the association constant Ka or also called binding constant. The most common and simplest model reported in literature to study protein-ligand is the binding equilibria for a first-order reaction, where the ligands (L) are entering one by one at the binding site of the protein (P) .
Assuming that the sites are equivalents and independents, the association constant is
where [P] is the free protein concentration, [L] is the free quencher concentration, and [LnP] is the concentration of protein-quencher complex. In this model, the protein is either free in solution or bound to quencher, then
where [Pt] is the concentration of total protein in the system . Considering that the fluorescence intensity is only emitted by free proteins
Eq. (2) can be rearranged to
Applying the logarithm function in Eq. (5) and rearranging it
where R is the universal gas constant (≈ 8.31 J/mol.K) and T is the temperature. In a graphic of ln Ka versus T, the enthalpy change (ΔH) and the entropy change (ΔS) are obtained by the slope of the linear function and the linear coefficient, respectively.
Besides that, Gibbs free energy variation says about the spontaneity of the complex formation process, where ∆G > 0 indicates a non-spontaneity process and ∆G < 0 indicates a spontaneity process. The enthalpy variation, on the other hand, indicates if the system is either an exothermic process (∆H < 0) or an endothermic process (∆H > 0) [24, 25].
Unlike the binding equilibrium method that is based upon a first-order reaction predicting a single binding site, there is another method that does not make use of a previous model, and in addition to determining the number of binding sites, it can still determine the cooperativity, if any. This method was developed by Scatchard , and it is not as popular as binding equilibrium method to characterize protein-ligand interaction in the scientific community yet.
This method makes use of another technique to obtain preliminary data known as binding density function (BDF) , where a physical observable is chosen to follow the interactions, such as absorbance and fluorescence intensity, among others. Supposing that the intensity of fluorescence emitted by the protein was chosen as physical observable, in which in an aqueous system containing protein is titrated by the ligand. Considering the system in equilibrium, the average number of quencher bound by protein Σνi is determined from a given free quencher concentration (L); in this context, if the free quencher concentration is the same for two or more solutions of different concentrations of total protein Pt, then the average distribution of the binding density of the quencher will also be the same, reaching in the expression for the mass conservation :
where [Lt] is the total quencher concentration, [L] is the free quencher concentration, and [Pt] is the total protein concentration. In order to obtain the average number of quenchers bound (Σνi), a graphic of total quencher concentration [Lt] versus the total protein [Pt] can be plotted, in which by the angular coefficient, the Σνi is obtained and by the linear coefficient of the linear function [L] is obtained. The values for [Pt] and [Lt] are obtained from the graphic of ΔF versus log [Lt] shown in Figure 4, where ΔF is the suppression percentage .
Once the values of Σνi and [L] were obtained from the graphic in the BDF method, the method developed by Scatchard  can be used to follow protein-ligand interaction. To determine all those information from the method, a graphic of Σνi/[L] versus Σνi is plotted, and from the function obtained in the graphic, it can be said that the sites have cooperativity or they are equivalents and independents.
Bordbar and co-workers  reported the mathematical functions that describe the behavior of the system in the Scatchard plots with the type of cooperativity. They determined that if the function is polynomial with the positively concavity (Figure 5(A)), the protein has negative cooperativity between the sites. If instead there is a polynomial function with negative concavity, the protein has a positive cooperativity among the sites (Figure 5(B)) .
For these cases where the Scatchard graph shows a cooperativity function shape, one can apply the equation developed by Hill  to determine the number of sites n of each set of equal sites and the Ka association constant of each set of interaction sites:
where [L] is the concentration of free quenchers and H is Hill’s index. If H = 1 the system is noncooperative, H > 1 the system has positive cooperativity, and H < 1 the system has negative cooperativity [24, 28].
In the case that the function, which describes the system, is linear, the protein has no cooperativity among the sites, and all sites are identical and independent. For this behavior, Scatchard  developed his own mathematical model to find the number of sites n and the association constant Ka:
To rearrange in the form of linear equation to model the Scatchard graphic, it is given as follows:
3. Theory of molecular docking
Molecular docking is a powerful technique that has been used in association with experimental data to determine ligand binding sites in targets with pharmacological interest. It can also predict the ligand conformation in the binding site and consequently the interactions that stabilize the complex. In addition, molecular docking has been also used as an efficient and a cheap technique for virtual screening large molecule databases and selects compounds which bind with specificity in pharmaceutical targets before carrying out in vitro and in vivo experiments [29, 30].
Molecular docking consists in determining the most probable conformations of the complex composed by receptor ligand based on an energy ranking of each conformation. To obtain this energy ranking, the ligand is put to interact in an environment under a protein force field including non-covalent potentials such as van der Waals force, hydrogen bonds, and electrostatic nature forces . There are several softwares to perform this kind of prediction such as AutoDock, GOLD, DOCK, and FlexX, just to mention a few.
Each of those mentioned software has its own potential mathematical function to elucidate the forces involved in the complex and calculate the energy score to be ranked. The following is an example of potential utilized by AutoDock04  software:
The weighting constants I, J, E(t), K, and W are those optimized to calibrate the empirical free energy based on a set of experimentally characterized complexes. The first term is the Lenard-Jones potential, in which parameters A and B are taken from the Amber force field . The second term refers to the hydrogen bond in which the parameters C and D are obtained to ensure a minimum energy of 5 kcal/mol in 1.9 Å for O-H and N-H and 1 kcal/mol in 2.5 Å for S-H. The function E(t) provides directionality based on the angle t of the geometry of an ideal hydrogen bond . The third term is a shielded Coulomb potential for electrostatic interaction. The last term is the desolvation potential based on the volume of the atoms surrounding a given atom and shelter it of the solvent .
To reduce the computational costs, the programs usually make use of a pre-calculation type. The software creates grid maps for each type of atom present in the ligand to be used in the docking. Grid maps consist of a three-dimensional array of regularly spaced points centered on the active site of the protein or macromolecule under study. Each point inside the grid maps records the energy interaction of a test atom with the protein . The complexity of finding the best conformation requires computational methods with the potential to effectively investigate a large number of possible solutions, aiming to find the best result. The search algorithms that are usually utilized by molecular docking softwares can be classified into three categories such as systematic, deterministic, and stochastic search methods .
In systematic search algorithms, each degree of freedom has a set of values, so that all degrees of freedom of the molecule are explored combinatorically during the search. Examples of systematic search algorithms are anchor-and-grow or incremental construction algorithm . The deterministic search methods are characterized by the fact that, given the same initial input state, they always produce the same output. It happens because the initial state determines the possible movement to generate the next state, which has to have the same or lower energy than the previous state. Examples of deterministic search are energy minimization methods and molecular dynamics (MD) simulations . Stochastic methods vary randomly in all degrees of freedom of the ligand (translational, rotational, and conformational) at each step, generating a great diversity of solutions. The solutions are evaluated according to a probabilistic criterion to decide if they will be accepted or not. In this way this method requires a large number of conformations to obtain a desired result. Monte Carlo, simulated annealing, and evolutionary methods are examples of stochastic search method, and they are the most common in molecular docking softwares [32, 33].
4. Molecular targets for flavonoids
Many efforts have been done to describe flavonoids’ main biological targets involved in cellular processes such as inflammation and cancer. Flavonoid mechanism of action in cellular responses includes the inhibition of proteins in cytoplasm medium such as IκB kinase (IKK) complex and mitogen-activated protein kinases (MAPKs), as well as in extracellular medium such as interleukins and other cytokines.
In the pathway that triggers the activation of nuclear factor κB (NF-κB) in tumor necrosis factor (TNF) stimulated cells (Figure 6), IKK complex is a target for many flavonoids such as apinegin , luteolin , acacetin , ochnaflavone , delphinidin , and morin .
The inhibition of IKK complex results in failure in phosphorylation of Iκβ. As a consequence, NF-κB is not released to translocate to the nucleus. In the case of the flavonoids apinegin and luteolin, they can also inhibit NF-κβ activation by inhibition of MAPKs .
In the pathway that triggers the activation of activator protein 1 (AP-1), c-Jun N-terminal kinase (JNK) is a target for flavonoids kaempferol and chrysin . Such flavonoids bind to the protein and inhibit the activation of AP-1. Sumaflavone also suppresses the AP-1 activation, but the main target for this flavonoid has not been identified yet .
In LPS-stimulated cells (Figure 7), epigallocatechin-3-gallate (EGCG)  and morin  inhibited the activation of NF-κB activity by IKK complex inhibition. Besides that, EGCG also inhibited p38 MAP kinase activity . Another molecular target in LPS-stimulated cells is JNK, which is inhibited for galangin .
The molecules genistein, kaempferol, quercetin, and daidzein inhibited the activation of STAT-1 and NF-kB in LPS-induced cells, while flavone, isorhamnetin, naringenin, and pelargonidin inhibited only the NF-kB activation. However the molecular targets for such flavonoids have not been identified yet [39, 40].
Cytokines can also be targets for flavonoids. As verified by Li and co-workers , baicalin is able to complex with a variety of chemokines such as SDF-1a, IL-8, MIP1-b, and MCP-2 and reduce the capacity of the cytokines to bind and activate their receptors.
Despite the works described above show which proteins are inhibited by flavonoids, the methodology used does not give information related to how strong is the affinity of the flavonoid for the protein, which amino acids of the protein participate in the interaction with flavonoids, which are the molecular driving forces involved in the interaction, the thermodynamic parameters of complex formation, how many binding sites there are in the protein for the flavonoids, and whether the protein presents cooperativity among their sites. In this way, the methodology of fluorescence spectroscopy and molecular docking described in Sections 2 and 3 come to help give this set of information that are of the utmost importance to drug discovery.
Among the target proteins mentioned above, NF-κβ (p50 and p65 subunits) has one tryptophan residue, while JNK have four tryptophan residues. Besides that, all the cytokines (SDF-1a, IL-8, and MCP-2) mentioned above have 1 tryptophan residue, except for MIP1-b which has 12. The presence of tryptophan makes the use of fluorescence spectroscopy possible, but the presence of more than one makes the interpretation of the spectral data more complicated. Different of the proteins mentioned above, AP-1 does not have tryptophan residue. On the other hand, there are tyrosine residues, which can also be followed in the experiments.
Besides the endogenous fluorescent probe, solubility is another important factor because fluorescence spectroscopy experiments are conducted with protein in solution. As the proteins mentioned above are found in cytoplasm or in extracellular medium, they are expected to be soluble in buffer.
5. Physical pharmacy in the description of molecular mechanism of flavonoids
As already discussed in the sections before, the stabilization of the complex is made by non-covalent interactions, such as hydrogen bonds, hydrophobic interactions, van der Waals and electrostatic forces, salt bridge, π stacking, and cation-π interactions. Consequently, the affinity and the thermodynamic parameters of each system are results of these interactions.
The complete characterization of a complex involves the description of these interactions, the affinity, and the thermodynamics. That is the reason why the use of techniques based on physical approaches in the pharmacy area has been growing.
Several authors characterized the complex formed by flavonoids and biological targets based on fluorescence spectroscopy data, as showed in Section 5.1. Some others used molecular docking to describe the complexes, as shown in Section 5.2. Nevertheless, the most accurate description of the complex is reached by the association of data from fluorescence and docking, as one can see in Section 5.3.
5.1 The use of steady-state fluorescence spectroscopy
Based on the theory of fluorescence quenching described in Section 2, several authors investigated the molecular mechanism of interaction between flavonoids and proteins. The most basic method used to estimate the affinity of the ligand toward the protein is the Stern-Volmer plot, where one can calculate the Stern-Volmer constant (Ksv) and associate the magnitude of this constant with the binding affinity.
Papadopoulou and co-workers  used Ksv to evaluate the interaction of two flavonoids, quercetin and rutin, and the bovine serum albumin (BSA). The authors verified that quercetin demonstrated stronger affinity toward BSA compared to rutin. Such difference was explained based on the structure of both flavonoids, once rutin is a glucoside of quercetin. The authors concluded that hydrophobicity was an important factor in the change of affinity, once the incorporated disaccharide rutinose made the flavonoid less hydrophobic. Besides that, another important factor was the steric hindrance caused by the incorporated disaccharide.
With a similar methodology, Cao and co-workers  investigated the influence of glycosylation of quercetin and baicalein in the affinity toward BSA. The authors used the Stern-Volmer plot to classify the quenching mechanism as static or dynamic. Then they used the double-logarithm plot to calculate the binding affinity of the ligands studied. They verified that the glycosylation made the interaction weaker and suggested that the decrease in binding constant was an effect of steric hindrance caused by glycoside groups.
The influence of glycosylation in the binding affinity observed by Papadopoulou and co-workers  and Cao and co-workers  was attributed to steric hindrance. The steric hindrance suggested by the authors could be verified if they had used molecular docking to model the system.
The influence of ions Cu2+, Al3+, Mg2+, and Zn2+ in the interaction of flavonoids and human serum albumin (HSA) was investigated by Bi and co-workers . The authors calculated the binding constant by double-logarithm plot (Eq. (6)), and they concluded that binding constants were 14.2–99.6% of the ones without these metal ions. Therefore, the ion concentration would shorten the storage time of the compounds in blood plasma and enhance the effectiveness of the flavonoids. The influence of ions in the binding constant was also observed by Hu and co-workers  during experiments with morin and BSA. In this case, the authors observed that the constant decreased in the presence of Ba2+ and Hg2+, but with the addition of the ions K+, Li+, Mn2+, Fe2+ and Sn2+, the constant increased.
Another experimental methodology has been used to characterize the interaction between protein and small ligands. As described in Section 2, thermodynamic parameters are calculated from fluorescence quenching data by using the van’t Hoff equation.
Li and co-workers  used fluorescence experiments to compare the inhibition ability of three flavonoids (quercetin, isoquercetin, and rutin) toward
Xi and co-workers  studied the interaction between hemoglobin and the flavonoids quercetin and rutin by using double-logarithm and Van’t Hoff plots. The authors observed that for both flavonoids, the values of ΔH and ΔS were negative and concluded that the acting forces are hydrogen bonding and van der Waals. As one can see, the results obtained by [45, 46] reveal that molecular forces that drive quercetin and rutin into
Another molecular target for flavonoids in biological systems is nucleic acid. The study of the interaction between flavonoids and DNA using fluorescence spectroscopy is made by following fluorescence signal from flavonoids once DNA exhibits a very weak intrinsic emission.
The article published by Jana and co-workers  described the affinity of 3-hydroxyflavone toward DNA, by following the flavonoid fluorescence signal. The same methodology was used by Sengupta and co-workers  to describe the affinity of fisetin toward DNA. The interaction of polyphenolic compounds with DNA has a protective role, once flavonoids exhibit antioxidant properties.
DNA is also a target for luteolin, as shown in the article published by Chowdhury and co-workers ; the interaction between luteolin and nucleic acid is one of the causes of the decrease in affinity between DNA and
5.2 The use of molecular docking
The influence of flavonoid structure, the addition of organic group as well as the main interaction between the small molecules and biological targets can be described through the molecular docking technique.
The potential of some flavonoids as inhibitors of
The inhibitory effect of some flavonoids was also observed for
Si and co-workers  showed that some flavonoids inhibiting CYP2C9 activity may increase the risk of toxicity from coadministered drugs that are CYP2C9 substrates. The authors used molecular docking to give details about the interaction of CYP2C9–6-hydroxyflavone complex, showing that the 6-hydroxyflavone is bound by a π-π stacking interaction with the phenyl group of Phe 100 and by two hydrogen bonds with Leu102 and Phe100.
The examples above used molecular docking to describe the interaction between flavonoids and molecular targets in order to explain the results obtained from experimental assays. However, one can also use molecular docking in order to select the best candidates and then submit the chosen molecules to the experimental assays.
In the work published by Salam and co-workers , the authors built a library of 200 natural compounds including some flavonoids. With virtual screening, they selected the 29 best candidates to bind to the peroxisome proliferator-activated receptors (PPARs) in the cell culture experiments. Besides that, docking studies showed that flavonoids are predicted to occupy the hydrophobic environment formed by residues Phe282, Phe360, and Phe360. They also showed that there is one specific OH group in all the flavonoids studied that is responsible to make a hydrogen bond with the receptor, except in the case of apigenin.
5.3 The use of fluorescence spectroscopy associated with molecular docking
Data from fluorescence spectroscopy offers a description of the complex in terms of binding constant, number of binding sites, cooperativity, and the thermodynamics of the complex formation. On the other hand, molecular docking offers details about the binding site environment, the ligand conformation, and the interactions that stabilize the complex. Some authors associated the fluorescence quenching data with the results from molecular docking in order to give a complete description of the biological system.
He and co-workers  verified the tryptophan fluorescence quenching of HSA caused by alpinetin. The authors used the van’t Hoff model to calculate the thermodynamic parameters and concluded that the hydrophobic interactions were predominant to stabilize the complex. With molecular docking technique, the authors selected the best binding site of alpinetin in the HSA and showed that the binding site was close to tryptophan, which explains the static quenching. The authors also showed that the binding site was a hydrophobic cavity, which explains the data obtained by van’t Hoff equation.
With a similar methodology, Kim and co-workers  showed the relationship between the number of OH groups in flavonoids and the affinity for mushroom
Caruso and co-workers  used fluorescence spectroscopy to describe the interaction between isovitexin and HSA. The authors calculated the contribution of enthalpy and entropy by van’t Hoff equation and concluded that the process was enthalpically driven. The complete description of the system was reached with molecular docking, as one can see in Figure 8. The oxygen atoms of Ala191 and Pro339 form hydrogen bonds with isovitexin, while neutral polar side chain residues like Gln221, Asn295, Tyr341, and Tyr452 are mainly taking part in electrostatic interactions (charge neutralization).
In this last article, the authors used the Scatchard method, which has been recently used to calculate the number of binding sites and to verify allosterism. With the Scatchard method, the authors concluded that there was one binding site for isovitexin in HSA, which was in agreement with the results obtained by the authors when they used the binding equilibrium model.
In another article, Caruso and co-workers  used a similar methodology in order to describe the interaction between 2-phenylchromone and HSA. In this case, the contributions of enthalpy and entropy calculated by van’t Hoff equation showed that the process was entropically driven. With the Scatchard method, the authors concluded that there were two cooperative binding sites for 2-phenylchromone in HSA with two different binding constants. Note that this type of information could not be obtained from Stern-Volmer or double-logarithm models; only the Scatchard method can offer information about the system without the use of any binding model a priori.
Also in this study, molecular docking showed that 2-phenylchromone binds at subdomain IIA (Site 1) and subdomain IIIA (Site 2). As shown in Figure 9, at Site 1 amino acids Lys199, Leu (219, 234, 238, and 260), Ala291, and Glu292 are residues involved in the interactions with the molecule. At Site 2 amino acids Leu387, Tyr411, Leu453, Arg485, and Phe488 compose the microenvironment of interaction.
This chapter presents the real possibilities of combinations, between physical experimental techniques with computational tools which can contribute significantly to the advancement of the proposal of new drugs within the point of view of molecular action. Moreover, it has been demonstrated how to carry out such work, combining the knowledge between fluorescence spectroscopy and molecular docking. These techniques may be added to pharmacology, pharmacokinetics, and pharmacognosy research fields, which represent the cradle of the search for new medicinal products and design.
The author A.P.R.P. gratefully acknowledges a CAPES scholarship. The author G.Z. gratefully acknowledges a CNPq scholarship (grant number: 141953/2017-9). The author M.L.C. gratefully acknowledges the financial support received from FAPESP (grant number: 2017/08834-9).
Conflict of interest
The authors declare that no conflict of interests exists.