Antifungal and Antiparasitic Activities of Metallocene-Containing Fluconazole Derivatives.

The search for new anti-infectives based on metal complexes is gaining momentum. Among the different options taken by researchers, the one involving the use of organometallic complexes is probably the most successful one with a compound, namely, ferroquine, already in clinical trials against malaria. In this study, we describe the preparation and in-depth characterization of 10 new (organometallic) derivatives of the approved antifungal drug fluconazole. Our rationale is that the sterol 14α-demethylase is an enzyme part of the ergosterol biosynthesis route in Trypanosoma and is similar to the one in pathogenic fungi. To demonstrate our postulate, docking experiments to assess the binding of our compounds with the enzyme were also performed. Our compounds were then tested on a range of fungal strains and parasitic organisms, including the protozoan parasite Trypanosoma cruzi (T. cruzi) responsible for Chagas disease, an endemic disease in Latin America that ranks among some of the most prevalent parasitic diseases worldwide. Of high interest, the two most potent compounds of the study on T. cruzi that contain a ferrocene or cobaltocenium were found to be harmless for an invertebrate animal model, namely, Caenorhabditis elegans (C. elegans), without affecting motility, viability, or development.

Platinum and Palladium Organometallic Compounds: Disrupting the Ergosterol Pathway in Trypanosoma cruzi.

Current treatment for Chagas' disease is based on two drugs, Nifurtimox and Benznidazol, which have limitations that reduce the effectiveness and continuity of treatment. Thus, there is an urgent need to develop new, safe and effective drugs. In previous work, two new metal-based compounds with trypanocidal activity, Pd-dppf-mpo and Pt-dppf-mpo, were fully characterized. To unravel the mechanism of action of these two analogous metal-based drugs, high-throughput omics studies were performed. A multimodal mechanism of action was postulated with several candidates as molecular targets. In this work, we validated the ergosterol biosynthesis pathway as a target for these compounds through the determination of sterol levels by HPLC in treated parasites. To understand the molecular level at which these compounds participate, two enzymes that met eligibility criteria at different levels were selected for further studies: phosphomevalonate kinase (PMK) and lanosterol 14-α demethylase (CYP51). Molecular docking processes were carried out to search for potential sites of interaction for both enzymes. To validate these candidates, a gain-of-function strategy was used through the generation of overexpressing PMK and CYP51 parasites. Results here presented confirm that the mechanism of action of Pd-dppf-mpo and Pt-dppf-mpo compounds involves the inhibition of both enzymes.

Antioxidant Capacity and NF-kB-Mediated Anti-Inflammatory Activity of Six Red Uruguayan Grape Pomaces.

Grape pomaces have a wide and diverse antioxidant phenolics composition. Six Uruguayan red grape pomaces were evaluated in their phenolics composition, antioxidant capacity, and anti-inflammatory properties. Not only radical scavenging methods as DPPH· and ABTS·+ were employed but also ORAC and FRAP analyses were applied to assess the antioxidant potency of the extracts. The antioxidant reactivity of all extracts against hydroxyl radicals was assessed with ESR. The phenol profile of the most bioactive extract was analyzed by HPLC-MS, and a set of 57 structures were determined. To investigate the potential anti-inflammatory activity of the extracts, Nuclear Factor kappa-B (NF-κB) modulation was evaluated in the human colon cancer reporter cell line (HT-29-NF-κB-hrGFP). Our results suggest that Tannat grapes pomaces have higher phenolic content and antioxidant capacity compared to Cabernet Franc. These extracts inhibited TNF-alpha mediated NF-κB activation and IL-8 production when added to reporter cells. A molecular docking study was carried out to rationalize the experimental results allowing us to propose the proactive interaction between the NF-κB, the grape extracts phenols, and their putative anti-inflammatory bioactivity. The present findings show that red grape pomace constitutes a sustainable source of phenolic compounds, which may be valuable for pharmaceutical, cosmetic, and food industry applications.

Identification of potential biological target for trypanocidal sesquiterpene lactones derivatives.

Sesquiterpene lactones are natural products of the Asteraceae family that have shown trypanocidal activity against Trypanosoma cruzi, even exceeding the effectiveness of drugs used in the treatment of American trypanosomiasis. However, there is no agreement on their mechanism of action and their specificity to interact with parasite proteins. For this reason, we aimed to find biological targets that can interact with these compounds by reverse virtual screening with ligand pharmacophores and putative binding sites and the use of bioinformatic databases. Therefore, 41 possible biological targets were found, and four of them (with crystallized proteins), interfering directly or indirectly in the trypanosomatid redox system, were studied in detail. As a first approach, we focused on the study of trypanothione reductase, and protein-ligand interaction fingerprint analyses were performed to find binding site determinants that promote a possible inhibition of the enzyme. This study contributes to the understanding of one of the putative mechanisms of action of sesquiterpene lactones on one of the numerous suggested targets.Communicated by Ramaswamy H. Sarma.

Modeling, molecular dynamics and docking studies of a full-length Echinococcus granulosus 2DBD nuclear receptor.

Nuclear receptors are ligand-activated transcription factors capable of regulating the expression of complex gene networks. The family includes seven subfamilies of protein with a wide phylogenetic distribution. A novel subfamily with two DNA-binding domains (2DBDs) has been first reported in Schistosoma mansoni (Platyhelminth, Trematoda). Employing an ab initio protocol and homology modeling methods, the full-length 3D structure of the Eg2DBDα.1 nuclear receptor from Echinococcus granulosus (Platyhelminth, Cestoda) was generated. The model analysis reveals the presence of the conserved three-layered alpha-helical sandwich structure in the ligand binding domain, and a particularly long and flexible hinge region. Molecular dynamics simulations were performed previous to dock a conformational library of fatty acids and retinoic acids. Our results indicate that oleic and linoleic acids are suitable ligands to this receptor. The ligand-protein complex is stabilized mainly by hydrogen bonds and hydrophobic interactions. The fact that 2DBD nuclear receptors have not been identified in vertebrates confers particular interest to these nuclear receptors, not only concerning their structure and function but as targets of new anthelmintic drugs.Communicated by Ramaswamy H. Sarma.

Mode of action of p-quinone derivatives with trypanocidal activity studied by experimental and in silico models.

Quinones are attractive pharmacological scaffolds for developing new agents for the treatment of different transmissible and non-transmissible human diseases due to their capacity to alter the cell redox homeostasis. The bioactivity and potential mode of action of 19 p-quinone derivatives fused to different aromatic rings (carbo or heterocycles) and harboring distinct substituents were investigated in infective Trypanosoma brucei brucei. All the compounds, except for a furanequinone (EC50=38 µM), proved to be similarly or even more potent (EC50 = 0.5-5.5 µM) than the clinical drug nifurtimox (EC50 = 5.3 µM). Three furanequinones and one thiazolequinone displayed a higher selectivity than nifurtimox. Two of these selective hits resulted potent inhibitors of T. cruzi proliferation (EC50=0.8-1.1 µM) but proved inactive against Leishmania infantum amastigotes. Most of the p-quinones induced a rapid and marked intracellular oxidation in T. b. brucei. DFT calculations on the oxidized quinone (Q), semiquinone (Q•-) and hydroquinone (QH2) suggest that all quinones have negative ΔG for the formation of Q•-. Qualitative and quantitative structure-activity relationship analyses in two or three dimensions of different electronic and biophysical descriptors of quinones and their corresponding bioactivities (killing potency and oxidative capacity) were performed. Charge distribution over the quinone ring carbons of Q and Q.- and the frontier orbitals energies of SUMO (Q.-) and LUMO (Q) correlate with their oxidative and trypanocidal activity. QSAR analysis also highlighted that both bromine substitution in the p-quinone ring and a bulky phenyl group attached to the furane and thiazole rings (which generates a negative charge due to the π electron system polarized by the nearby heteroatoms) are favorable for activity. By combining experimental and in silico procedures, this study disclosed important information about p-quinones that may help to rationally tune their electronic properties and biological activities.

A Computer Simulation Insight into the Formation of Apocarotenoids: Study of the Carotenoid Oxygenases BCO1 and BCO2 and Their Interaction with Putative Substrates.

Carotenoids are isoprenoid pigments, and sources of vitamin A in humans. The first metabolic pathway for their synthesis is mediated by the enzymes ß,ß-carotene-15,15'-dioxygenase (BCO1) and ß,ß-carotene-9',10'-dioxygenase (BCO2), which cleave carotenoids into smaller compounds, called apocarotenoids. The objective of this study is to gain insight into the interaction of BCO1 and BCO2 with carotenoids, adding structural diversity and importance in the agro-food and/or health sectors. Homology modeling of BCO1 and BCO2, and the molecular dynamics of complexes with all carotenoids were performed. Interaction energy and structures were analyzed. For both enzymes, the general structure is conserved with a seven beta-sheet structure, and the ß-carotene is positioned at an optimal distance from the catalytic center. Fe2+ forms in an octahedral coordination sphere with four perfectly conserved histidine residues. BCO1 finds stability in a structure in which the ß-carotene is positioned ready for enzymatic catalysis at the 15-15' bond, and BCO2 in positioning the bond to be cleaved (C9-C10) close to the active site. In BCO1 the carotenoids interact with only seven residues with aromatic rings, while the interaction of BCO2 is much more varied in terms of the type of interaction, with more residues of different chemical natures.

Characterization of the B-Cell Epitopes of Echinococcus granulosus Histones H4 and H2A Recognized by Sera From Patients With Liver Cysts.

Cystic echinococcosis (CE) is a zoonotic disease worldwide distributed, caused by the cestode Echinococcus granulosus sensu lato (E. granulosus), with an incidence rate of 50/100,000 person/year and a high prevalence in humans of 5-10%. Serology has variable sensitivity and specificity and low predictive values. Antigens used are from the hydatid fluid and recombinant antigens have not demonstrated superiority over hydatid fluid. A cell line called EGPE was obtained from E. granulosus sensu lato G1 strain from bovine liver. Serum from CE patients recognizes protein extracts from EGPE cells with higher sensitivity than protein extracts from hydatid fluid. In the present study, EGPE cell protein extracts and supernatants from cell colonies were eluted from a protein G affinity column performed with sera from 11 CE patients. LC-MS/MS proteomic analysis of the eluted proteins identified four E. granulosus histones: one histone H4 in the cell extract and supernatant, one histone H2A only in the cell extract, and two histones H2A only in the supernatant. This differential distribution of histones could reflect different parasite viability stages regarding their role in gene transcription and silencing and could interact with host cells. Bioinformatics tools characterized the linear and conformational epitopes involved in antibody recognition. The three-dimensional structure of each histone was obtained by molecular modeling and validated by molecular dynamics simulation and PCR confirmed the presence of the epitopes in the parasite genome. The three histones H2A were very different and had a less conserved sequence than the histone H4. Comparison of the histones of E. granulosus with those of other organisms showed exclusive regions for E. granulosus. Since histones play a role in the host-parasite relationship they could be good candidates to improve the predictive value of serology in CE.

Structure-Based Bioisosterism Design, Synthesis, Biological Evaluation and In Silico Studies of Benzamide Analogs as Potential Anthelmintics.

A recent screen of 67,012 compounds identified a new family of compounds with excellent nematicidal activity: the ortho-substituted benzamide families Wact-11 and Wact-12. These compounds are active against Caenorhabditis elegans and parasitic nematodes by selectively inhibiting nematode complex II, and they display low toxicity in mammalian cells and vertebrate organisms. Although a big number of benzamides were tested against C. elegans in high-throughput screens, bioisosteres of the amide moiety were not represented in the chemical space examined. We thus identified an opportunity for the design, synthesis and evaluation of novel compounds, using bioisosteric replacements of the amide group present in benzamides. The compound Wact-11 was used as the reference scaffold to prepare a set of bioisosteres to be evaluated against C. elegans. Eight types of amide replacement were selected, including ester, thioamide, selenoamide, sulfonamide, alkyl thio- and oxo-amides, urea and triazole. The results allowed us to perform a structure-activity relationship, highlighting the relevance of the amide group for nematicide activity. Experimental evidence was complemented with in silico structural studies over a C. elegans complex II model as a molecular target of benzamides. Importantly, compound Wact-11 was active against the flatworm Echinococcus granulosus, suggesting a previously unreported pan-anthelmintic potential for benzamides.

Modelling of the Citrus CCD4 Family Members: In Silico Analysis of Membrane Binding and Substrate Preference.

Coloring is one of the most important characteristics in commercial flowers and fruits, generally due to the accumulation of carotenoid pigments. Enzymes of the CCD4 family in citrus intervene in the generation of ß-citraurin, an apocarotenoid responsible for the reddish-orange color of mandarins. Citrus CCD4s enzymes could be capable of interacting with the thylakoid membrane inside chloroplasts. However, to date, this interaction has not been studied in detail. In this work, we present three new complete models of the CCD4 family members (CCD4a, CCD4b, and CCD4c), modeled with a lipid membrane. To identify the preference for substrates, typical carotenoids were inserted in the active site of the receptors and the protein-ligand interaction energy was evaluated. The results show a clear preference of CCD4s for xanthophylls over aliphatic carotenes. Our findings indicate the ability to penetrate the membrane and maintain a stable interaction through the N-terminal α-helical domain, spanning a contact surface of 2250 to 3250 Å2. The orientation and depth of penetration at the membrane surface suggest that CCD4s have the ability to extract carotenoids directly from the membrane through a tunnel consisting mainly of hydrophobic residues that extends up to the catalytic center of the enzyme.

Peptides derived from the gastrointestinal digestion of amaranth 11S globulin: Structure and antioxidant functionality.

The relationship between structural and physicochemical properties and antioxidant activity of peptides from amaranth 11S-globulin was studied. Peptides AWEEREQGSR, TEVWDSNEQ, IYIEQGNGITGM and YLAGKPQQEH had the greatest in vitro activity (ORAC, HORAC). GDRFQDQHQ, HVIKPPSRA and KFNRPETT were the most active ones against Cu+2/H2O2-induced-LDL oxidation. In a cellular system (H2O2-induced-Caco2-TC7), TEVWDSNEQ, IYIEQGNGITGM, GDRFQDQHQ, LAGKPQQEHSGEHQ and KFNRPETT were the most effective in decreasing ROS, while the effects on SOD, GPx, and GSH were variable. To understand the structure-antioxidant activity relationships, the content of aromatic and acidic amino acids, the degree of hydrophobicity and the charge distribution on the accessible surface of peptides structures obtained by molecular dynamics were analysed. The low correlation between in vitro, ex vivo and cellular activities could be explained by the influence of physicochemical and structural properties on the interaction with complex systems (LDL/cells), peptide modifications and/or mechanisms other than direct ROS inhibition in the cells.

New aryloxy-quinone derivatives with promising activity on Trypanosoma cruzi.

Continuing with a program to develop new quinone derivatives as biologically active compounds, we designed and synthesized a new series of aryloxy-quinones, which were evaluated in vitro against Trypanosoma cruzi in epimastigote form. Chemical modifications in three specific moieties on the aryloxy-quinone core were considered for developing new anti-T. cruzi agents. The majority of our new quinones showed higher potency (IC50 values of <0.70 µM) than nifurtimox, a known pharmaceutical used as a baseline drug (IC50 values of 7.00 µM); however, only two of them elicited higher selectivity than nifurtimox against Vero cells. A structure-activity relationship analysis provided information about the stereoelectronic features of these compounds, which are responsible for an increase in trypanosomicidal activity. Using a pharmacophore model, we mapped the substitution patterns of the five pharmacophoric features of trypanosomicidal activity. We chose the Epc1 compounds and found no relationship with the trypanosomicidal effects. These results provided useful information about the structural characteristics for developing new aryloxy-quinones with higher potency against the protozoan parasite T. cruzi.

Bioinformatic analysis of a novel Echinococcus granulosus nuclear receptor with two DNA binding domains.

Nuclear receptors are ligand-activated transcription factors capable of regulating the expression of complex gene networks. The family includes seven subfamilies of proteins with a wide phylogenetic distribution. A novel subfamily with two DNA-binding domains (2DBDs) has been reported in Schistosoma mansoni (Platyhelminth, Trematoda). This work describes the cDNA cloning and bioinformatics analysis of Eg2DBDα, a 2DBD nuclear receptor isoform from the parasite Echinococcus granulosus (Platyhelminth, Cestoda). The Eg2DBDα gene coding domain structure was analysed. Although two additional 2DBD nuclear receptors are reported in the parasite database GeneDB, they are unlikely to be expressed in the larval stage. Phylogenetic relationships between these atypical proteins from different cestodes are also analysed including S. mansoni 2DBD nuclear receptors. The presence of two DNA binding domains confers particular interest to these nuclear receptors, not only concerning their function but to the development of new antihelminthic drugs.

Combined molecular modelling and 3D-QSAR study for understanding the inhibition of NQO1 by heterocyclic quinone derivatives.

A combination of three-dimensional quantitative structure-activity relationship (3D-QSAR), and molecular modelling methods were used to understand the potent inhibitory NAD(P)H:quinone oxidoreductase 1 (NQO1) activity of a set of 52 heterocyclic quinones. Molecular docking results indicated that some favourable interactions of key amino acid residues at the binding site of NQO1 with these quinones would be responsible for an improvement of the NQO1 activity of these compounds. The main interactions involved are hydrogen bond of the amino group of residue Tyr128, π-stacking interactions with Phe106 and Phe178, and electrostatic interactions with flavin adenine dinucleotide (FADH) cofactor. Three models were prepared by 3D-QSAR analysis. The models derived from Model I and Model III, shown leave-one-out cross-validation correlation coefficients (q2LOO ) of .75 and .73 as well as conventional correlation coefficients (R2 ) of .93 and .95, respectively. In addition, the external predictive abilities of these models were evaluated using a test set, producing the predicted correlation coefficients (r2pred ) of .76 and .74, respectively. The good concordance between the docking results and 3D-QSAR contour maps provides helpful information about a rational modification of new molecules based in quinone scaffold, in order to design more potent NQO1 inhibitors, which would exhibit highly potent antitumor activity.

Structural evidence of quercetin multi-target bioactivity: A reverse virtual screening strategy.

The ubiquitous flavonoid quercetin is broadly recognized for showing diverse biological and health-promoting effects, such as anti-cancer, anti-inflammatory and cytoprotective activities. The therapeutic potential of quercetin and similar compounds for preventing such diverse oxidative stress-related pathologies has been generally attributed to their direct antioxidant properties. Nevertheless, accumulated evidence indicates that quercetin is also able to interact with multiple cellular targets influencing the activity of diverse signaling pathways. Even though there are a number of well-established protein targets such as phosphatidylinositol 3 kinase and xanthine oxidase, there remains a lack of a comprehensive knowledge of the potential mechanisms of action of quercetin and its target space. In the present work we adopted a reverse screening strategy based on ligand similarity (SHAFTS) and target structure (idTarget, LIBRA) resulting in a set of predicted protein target candidates. Furthermore, using this method we corroborated a broad array of previously experimentally tested candidates among the predicted targets, supporting the suitability of this screening approach. Notably, all of the predicted target candidates belonged to two main protein families, protein kinases and poly [ADP-ribose] polymerases. They also included key proteins involved at different points within the same signaling pathways or within interconnected signaling pathways, supporting a pleiotropic, multilevel and potentially synergistic mechanism of action of quercetin. In this context we highlight the value of quercetin's broad target profile for its therapeutic potential in diseases like inflammation, neurodegeneration and cancer.

Trypanothione Reductase: A Target for the Development of Anti- Trypanosoma cruzi Drugs.

BACKGROUND & OBJECTIVE: Chagas disease or American trypanosomiasis is a major parasitic disease in Latin America with restricted available treatment: nifurtimox and benznidazole. These two drugs are ineffective in the chronic phase of the disease; therefore, there is a need for the development of new, efficient and safe drugs for the treatment of this pathology. With this goal, one of the promising targets is trypanothione reductase (TR), a key enzyme in the metabolism of Trypanosoma cruzi. CONCLUSION: In this review, we analyse the importance of TR as a drug target, as well as the well-known and new inhibitors reported in the last decade as potential therapeutic agents for Chagas disease.

Structural analysis and molecular docking of trypanocidal aryloxy-quinones in trypanothione and glutathione reductases: a comparison with biochemical data.

A set of aryloxy-quinones, previously synthesized and evaluated against Trypanosoma cruzi epimastigotes cultures, were found more potent and selective than nifurtimox. One of the possible mechanisms of the trypanocidal activity of these quinones could be inhibition of trypanothione reductase (TR). Considering that glutathione reductase (GR) is the equivalent of TR in humans, biochemical, kinetic, and molecular docking studies in TR and GR were envisaged and compared with the trypanocidal and cytotoxic data of a set of aryloxy-quinones. Biochemical assays indicated that three naphthoquinones (Nq-h, Nq-g, and Nq-d) selectively inhibit TR and the TR kinetic analyses indicated that Nq-h inhibit TR in a noncompetitive mechanism. Molecular dockings were performed in TR and GR in the following three putative binding sites: the catalytic site, the dimer interface, and the nicotinamide adenine dinucleotide phosphate-binding site. In TR and GR, the aryloxy-quinones were found to exhibit high affinity for a site near it cognate-binding site in a place in which the noncompetitive kinetics could be justified. Taking as examples the three compounds with TR specificity (TRS) (Nq-h, Nq-g, and Nq-d), the presence of a network of contacts with the quinonic ring sustained by the triad of Lys62, Met400', Ser464' residues, seems to contribute hardly to the TRS. Compound Nq-b, a naphthoquinone with nitrophenoxy substituent, proved to be the best scaffold for the design of trypanocidal compounds with low toxicity. However, the compound displayed only a poor and non-selective effect toward TR indicating that TR inhibition is not the main reason for the antiparasitic activity of the aryloxy-quinones.

Identification of novel CAP superfamily protein members of Echinococcus granulosus protoscoleces.

Echinoccocus granulosus is the causative agent of Cyst Echinococcosis, a zoonotic infection affecting humans and livestock representing a public health and an economic burden for several countries. Despite decades of investigation an effective vaccine still remains to be found. Parasitic cysteine-rich secretory proteins, antigen 5 and pathogenesis-related 1 proteins (CAPs) have been proposed as vaccine candidates against helmith's infection. In this work we have identified two novel proteins of this superfamily expressed at the protoescoleces larval stage named EgVAL1 and EgVAL2. The open reading frame sequences were deduced. The aminoacidic sequence was analyzed and confronted against already known vertebrate' and helminth's proteins sequences in order to infer putative functions. Immunolocalization studies were also performed. The obtained data supported by immunolocalization studies and homology models suggest that these proteins could be involved in protease activity inhibition.

Toward the understanding of the molecular basis for the inhibition of COX-1 and COX-2 by phenolic compounds present in Uruguayan propolis and grape pomace.

Propolis and grape pomace have significant amounts of phenols which can take part in anti-inflammatory mechanisms. As the cyclooxygenases 1 and 2 (COX-1 and COX-2) are involved in said mechanisms, the possibility for a selective inhibition of COX-2 was analyzed in vitro and in silico. Propolis and grape pomace from Uruguayan species were collected, extracted in hydroalcoholic mixture and analyzed. Based on phenols previously identified, and taking as reference the crystallographic structures of COX-1 and COX-2 in complex with the commercial drug Celecoxib, a molecular docking procedure was devised to adjust 123 phenolic molecular models at the enzyme-binding sites. The most important results of this work are that the extracts have an overall inhibition activity very similar in COX-1 and COX-2, i.e. they do not possess selective inhibition activity for COX-2. Nevertheless, 10 compounds of the phenolic database turned out to be more selective and 94 phenols resulted with similar selectivity than Celecoxib, an outcome that accounts for the overall experimental inhibition measures. Binding site environment observations showed increased polarity in COX-2 as compared with COX-1, suggesting that polarity is the key for selectivity. Accordingly, the screening of molecular contacts pointed to the residues: Arg106, Gln178, Leu338, Ser339, Tyr341, Tyr371, Arg499, Ala502, Val509, and Ser516, which would explain, at the atomic level, the anti-inflammatory effect of the phenolic compounds. Among them, Gln178 and Arg499 appear to be essential for the selective inhibition of COX-2.

Analysis of cyclosporin A and a set of analogs as inhibitors of a T. cruzi cyclophilin by docking and molecular dynamics.

Cyclophilins (CyPs) are enzymes involved in protein folding. In Trypanosoma cruzi (T. cruzi), the most abundantly expressed CyP is the isoform TcCyP19. It has been shown that TcCyP19 is inhibited by the immunosuppressive drug cyclosporin A (CsA) and analogs, which also proved to have potent trypanosomicidal activity in vitro. In this work, we continue and expand a previous study on the molecular interactions of CsA, and a set of analogs modeled in complexes with TcCyP19. The modeled complexes were used to evaluate binding free energies by molecular dynamics (MD), applying the Linear Interaction Energy (LIE) method. In addition, putative binding sites were identified by molecular docking. In our analysis, the binding free energy calculations did not correlate with experimental data. The heterogeneity of the non-bonded energies and the variation in the pattern of hydrogen bonds suggest that the systems may not be suitable for the application of the LIE method. Further, the docking calculations identified two other putative binding sites with comparable scoring energies to the active site, a fact that may also explain the lack of correlation found. Kinetic experiments are needed to confirm or reject the multiple binding sites hypothesis. In the meantime, MD simulations at the alternative sites, employing other methods to compute binding free energies, might be successful at finding good correlations with the experimental data.