Docking simulation between HIV peptidase inhibitors and Trypanosoma cruzi aspartyl peptidase
BMC Research Notes volume 11, Article number: 825 (2018)
The low investment in research, diagnosis and treatment are factors that contribute to the continuity of Chagas’ disease as a neglected tropical diseases (NTDs). In this context, the repositioning of drugs represents a useful strategy, in the search for new chemotherapeutic approaches for NTDs. HIV aspartic peptidase inhibitors (HIV IPs) are good candidates for drug repurposing. Here, we modeled the three dimensional structure of an aspartyl peptidase of Trypanosoma cruzi, the causative agent of Chagas’ disease, aligned it to the HIV aspartyl peptidase and performed docking binding assays with the HIV PIs.
The 3D structure confirmed the presence of acid aspartic residues, which are critical to enzyme activity. The docking experiment revealed that HIV IPs bind to the active site of the enzyme, being ritonavir and lopinavir the ones with greater affinity. Benznidazole presented the worst binding affinity, this drug is currently used in Chagas’ disease treatment and was included as negative control. These results together with previous data on the trypanocidal effect of the HIV PIs support the hypothesis that a T. cruzi aspartyl peptidase can be the intracellular target of these inhibitors. However, the direct demonstration of the inhibition of T. cruzi aspartyl peptidase activity by HIV PIs is still a goal to be persuaded.
Chagas disease, caused by the kinetoplastid parasite Trypanosoma cruzi, affects 8 million people worldwide and about 10,000 people die of complications linked to the disease . The classical contamination route involves a triatomine vector bite, usually associated with poor habitation conditions in Latin America. However, human immigration associated with blood transfusions contributed to the spread of the disease to North America, Europe and some Western Pacific countries. This change in epidemiological scenario led to an increased report of co-infected patients with Chagas disease and acquired immune deficiency syndrome (AIDS) [2,3,4]. Even tough more than a century has passed since the discovery of Chagas disease, the drug repertoire is still based on benznidazole and nifurtimox, which are compounds with severe side effects, questionable specificity and efficacy [5, 6]. Although some promising drug candidates are under development, such as ravuconazole, posoconazole and fexinidazole [7,8,9], the questionable economic return of investing in drug development for a disease mainly associated with poverty discourages pharmaceutical companies to search for new chemotherapeutics against this parasitic illness .
In this context, drug-repurposing strategy has the potential to facilitate an effective chemotherapeutic approach [10, 11]. The marked improvement in the life expectancy of AIDS sufferers after the incorporation of HIV aspartic peptidase inhibitors (HIV IPs) into the drug cocktail, the so-called highly active anti-retroviral therapy (HAART), was due to reduction in viremia, recovery of the immune response and a direct action on opportunistic pathogens [reviewed in 3]. The last effect has been extensively demonstrated in Trypanosoma cruzi by our research group [3, 12,13,14] and others , although the mode of action and intracellular target of the compounds are still unknown and have been a matter of extensive research [3, 14, 16, 17]. Here, we identified in the genome of T. cruzi, CL-Brener strain, a homologue of the HIV retroviral peptidase and selected a crystalized protein for generating a three dimensional model, which was aligned with HIV aspartyl peptidase. Then, we performed docking experiments between the modeled 3D protein and HIV PIs. We have also included pepstatin A and benznidazole as a reference for positive and negative controls, respectively, since the former is a classical aspartyl peptidase inhibitor  and the latter a classical drug used in Chagas’ disease treatment with unrelated mechanism of action and binding site .
The crystal structures of the retroviral aspartic peptidase domain of DNA-damage inducible proteins (DDI-1 like) from T. cruzi were not available in the Protein Data Bank (PDB) and a comparative modeling was performed to predict their 3D structures. The amino acid sequences of these target proteins of T. cruzi (strain CL brenner; entries Tc00.1047053510155.40 and Tc00.1047053511585.40 from TriTrypdb) were submitted to BLASTp search through the PDB to retrieve similar proteins with available 3D structure. The crystal structures selected as templates were the retroviral peptidase-like domains of the DDI-1 like proteins from Homo sapiens (PDB ID 4RGH, sequence identity 43% for Tc00.1047053510155.40 and 44% for Tc00.1047053511585.40) and Saccharomyces cerevisiae (PDB ID 4Z2Z, sequence identity 41% for both T. cruzi sequences). The comparative modeling was carried out using the selected templates and the auto-model module from MODELLER 9v16 . The two final homodimeric models of aspartyl peptidase domains (presenting the residues 227-353/354) of the target proteins of T. cruzi were selected based on the DOPE scoring function. These models were called TcRP-A and TcRP-B for the retroviral aspartic peptidase domain of Tc00.1047053510155.40 and Tc00.1047053511585.40, respectively. An energy minimization using the Amber ff14SB force field through UCSF CHIMERA interface  was performed to improve the overall structure geometry of the both TcRP-A and TcRP-B models. The stereochemical quality of the two refined models was evaluated using the ERRAT and VERIFY 3D from SAVES server (http://services.mbi.ucla.edu/SAVES/) and MOLPROBITY (http://molprobity.biochem.duke.edu/).
A multiple sequence alignment was performed using the ClustalW through Multiple Sequence Viewer of Maestro (Schrödinger Suite) in order to evaluate the sequence identity and similarity between the aspartyl peptidase domain of DDI-1 like proteins from T. cruzi and the HIV-1 peptidase. The superposition of the 3D structures of the refined models and the HIV-1 peptidase (PDB ID 3OXC) and the root mean square deviation (RMSD) between them were obtained with the “super” command in Pymol v126.96.36.199 (http://www.pymol.org/). Further alignments and comparisons with these structures were carried out with the TM-align algorithm (http://zhanglab.ccmb.med.umich.edu/TM-align/) .
In order to predict the inhibitory potential and binding modes of HIV-PIs towards the aspartyl peptidase domain of DDI-1 like protein of T. cruzi, molecular docking assays were performed with the constructed TcRP-A and TcRP-B models. Besides that, pepstatin A, a classical aspartyl peptidase inhibitor , was included as a positive control, while benznidazole was used as a negative control, since it has an unrelated mode of action. These ligands were retrieved in the sdf format from PubChem and prepared by the ligand preparation module (LigPrep) of Maestro. As an outcome, ionization states at a pH of 3.5 and 5.0 and tautomers were generated, the specified chiral centers were retained, and the resulted molecules were energetically minimized using the OPLS-2005 force field.
Both T. cruzi protein models were prepared using the protein preparation wizard of Maestro, in which hydrogen atoms were added, and protonation states were determined at pH 3.5 with PROPKA. Grids were generated around the potential active site of the two prepared protein models using the Receptor Grid Generation module of Glide, also from Maestro. The grid box was set to have 35Å of edge with coordinates x = − 220.53, y = 32.3, and z = 16.96 for TcRP-A and x = 66.35, y = − 8.35, and z = 13.47 for TcRP-B, both coordinates were determined using the Asp248A, Asp248B and Arg279B (one of the potential substrate binding residues) as centroids. Following ligand and protein preparation, as well as grid generation, docking simulations were performed with the extra precision (XP) protocol from Glide (GlideXP) (version 6.7) . The potential inhibitors were docked in the potential active site of TcRP-A and TcRP-B domains. The inhibitors were ranked based on their docking scores.
In a similar manner, the HIV protein (PDB ID 3OXC) was prepared with protein preparation wizard at pH 5.0, once the optimal pH range described for HIV-1 peptidase activity is 4–6 in contrast to eukaryotic aspartyl peptidases, which have, in general, optimal pH range of 2–4 . The grid box was centered on the crystal structure of inhibitor saquinavir complexed with the studied HIV protein (x = 5.04, y = − 2.74, and z = 14.8) and had an edge of 35 Å. Besides the redocking of saquinavir inhibitor with XP protocol from Glide, all other well know HIV inhibitors were also docked, and benznidazole compound was used as negative control.
Results and discussion
To perform structural analysis and molecular docking studies, two 3D structures of the aspartyl peptidase domain of DDI-1 like proteins from T. cruzi were constructed through comparative modeling. For this purpose, the homodimeric retroviral aspartyl peptidase domains of DDI-1 like proteins from humans (PDB ID 4RGH) and yeast (PDB ID 4Z2Z) were used as templates. The models TcRP-A and TcRP-B were selected by considering their DOPE scores of − 29,695.480 and − 29,901.357, the lowest among the set of models generated for Tc00.1047053510155.40 and Tc00.1047053511585.40, respectively. The Ramachandran plot of the two final refined models indicated that 96.4% and 96% of the amino acids for TcRP-A and TcRP-B, respectively, are in favored regions, and 99.2% (in both cases) are in allowed regions (data not shown). The Errat overall quality factor for TcRP-A was 90.783 and 87.037 for TcRP-B. The Verify 3D server estimated that 95.31 and 97.66% of the residues of TcRP-A and TcRP-B, respectively, had an averaged 3D-1D score ≥ 0.2. These results indicate that the two refined models have good quality and are reliable for carrying out further computational analysis.
The presence of a retroviral aspartyl peptidase domain in the DDI-1 like proteins of T. cruzi suggests the structural similarity between the same domains in proteins of different organisms. In order to visualize these similarities, a comparison between the primary, secondary and 3D structures of the HIV-1 peptidase domain (PDB ID: 3OXC) and the generated models were performed (Fig. 1). The sequence alignment of TcRP-A/B and HIV-1 peptidase showed an overall identity of 12% and a similarity of 28% (Fig. 1a). Despite the low identity, the superimposition of the 3D structures confirms the similarity between these domains with an RMSD of 2.67 Å for TcRP-A and 2.98 Å for TcRP-B in relation to HIV-1 peptidase (Fig. 1b). The sequence-order independent structure comparisons using the TM-align emphasize these structural similarities with a TM-score value of ~ 0.71 for both models towards HIV-1 peptidase, indicating that these domains have the same fold. The potential active site of the retroviral aspartyl peptidase domain of T. cruzi (TcRP-B model) is shown in Fig. 1c. The models for the retroviral aspartyl peptidase domain were generated as homodimers with the active sites containing residues from both subunits A and B such as the catalytic aspartic peptidases (Asp238A and B). Besides other active site residues, the potential substrate-binding loop (residues Ala274-Gly280) described for retroviral aspartyl peptidase domain of yeast  are also indicated in Fig. 1c.
After the design and validation of the 3D models, we performed molecular docking simulations aiming to predict inhibitory potential (docking scores) and the intermolecular interactions of the HIV aspartic peptidase inhibitors towards the active sites of the two modeled retroviral aspartyl peptidase domains (TcRP-A and TcRP-B) of the DDI-1 like proteins from T. cruzi. The ligands were classified based on their docking scores, as shown in Table 1, and only the lowest values for each compound were considered in this evaluation. As expected, pepstatin A, a well-described and potent aspartyl peptidase inhibitor , had the lowest docking scores in comparison with the other ligands (except ritonavir for TcRP-B). In both modeled domains, pepstatin A interacts through hydrogen bonding with Gln252B and the catalytic residues Asp238A and B (Fig. 2A, D). Besides that, in the case of TcRP-A, pepstatin makes a hydrogen bond with Arg279A, Arg253A and B, and Gly250A, while with TcRP-B, only one hydrogen bond with Asp312A is observed. Considering the HIV-PIs, ritonavir and lopinavir had lower docking scores to TcRPA and B (Table 1), it is conceivable that they have the potential to inhibit the activity of both domains. As observed for the positive control, these inhibitors form hydrogen bonds with the catalytic aspartates (Fig. 2B, C, E, F). On the other hand, atazanavir and benznidazole were the worst ranked, presenting higher docking scores due mainly to the poor interactions observed (Fig. 2G–J). These two compounds probably have no inhibitory activity or binding on the aspartyl peptidase domains of DDI-1 like protein from T. cruzi.
The structures of HIV and T. cruzi aspartyl peptidases present remarkable 3D similarities, to verify if the molecular docking outcomes are also paralleled, docking experiments were carried out using HIV peptidase. Although the docking scores of ligands in different proteins cannot be directly compared, the use of positive and negative controls allows the comparison of binding modes and intermolecular interactions. Therefore, we used benznidazole again as a negative control and all the know HIV aspartic peptidase inhibitors as positive controls, including pepstatin A, which is also described in the literature as an inhibitor of the HIV enzyme . As expected, benznidazole had the worst docking score, followed by pepstatin A. The top scored compounds were saquinavir (originally complexed in the studied HIV protein structure, presenting an RMSD of 1.03 Å in relation to the docking pose), atazanavir and nelfinavir. This result is in contrast to the T. cruzi protein that had pepstatin A, ritonavir and lopinavir as the better scored compounds (Table 1 and Additional file 1: Fig. S1). This contrast is expected since there was no rational design towards T. cruzi molecule. When comparing the predicted binding modes of these top scored inhibitors, similar key interactions can be observed in HIV and T. cruzi proteins, mainly with backbone and side chains of aspartates (including the catalytic ones) and glycines present in both active sites: Asp130B, Asp129B, Asp29A, Asp25A, Gly27A and Gly148B from HIV and Asp248A, Asp248B, Asp312A and Gly250A from T. cruzi (compare Fig. 2 and Additional file 1: Fig. S1). Therefore, the docking results predicted similar binding mechanisms, suggesting the binding of HIV PIs to T. cruzi DDI-1 protein in vitro.
DDI-1 gene product of T. cruzi, which shares homology to HIV aspartyl peptidase, is not biochemically characterized and has no crystallography structure resolved yet. In HIV and several fungi, HIV PIs target is an acidic aspartyl peptidase . Curiously, although the intracellular target in trypanosomatids is still unknown, HIV PIs have a strong trypanocidal effect [3, 13,14,15,16,17], and have also a strong leishmanicidal activity [reviewed in 26]. The HIV PIs are capable of inhibiting the enzymatic cleavage of cathepsin D fluorogenic peptide substrates by crude extracts of either T. cruzi or Leishmania spp. [12, 14, 27]. To evaluate the structure–function relationship and the possible binding of HIV PIs to T. cruzi aspartyl peptidase, we realized homology modeling and molecular docking studies of the catalytic domain of the DDI-1 like protein from T. cruzi. The 3D model generated exhibited structural similarities with the proteins from the A2 aspartyl peptidase family, suggesting that the retroviral domain of the DDI-1 like must be functional. Molecular docking revealed that eight out of nine HIV IPs tested presented binding affinity similar to the positive control. The intracellular target of HIV IPs in T. cruzi is still unknown and multifactorial effects exerted by these inhibitors may be responsible for disturbance in parasite homeostasis that culminates in cell death. Recently, we have demonstrated on Leishmania amazonensis that one of the possible mechanisms of action of HIV IPs is related to lipids metabolism disturbance [27, 28]. Therefore, studies aiming to demonstrate the direct binding of HIV IPs in vitro to the DDI-1 like protein from T. cruzi can provide further information on the intracellular target of these inhibitors.
An important step towards the determination of the intracellular target of the HIV PIs in T. cruzi would be the direct demonstration of the binding or the aspartyl peptidase activity inhibition by these compounds towards the purified enzyme. Up to now, our attempts to generate active recombinant T. cruzi DDI-1-like aspartyl peptidase in Escherichia coli BL-21 or to purify the enzyme from T. cruzi crude extracts through pepstatin-A affinity chromatography were unsuccessful.
acquired immune deficiency syndrome
- DDI-1 like:
DNA-damage inducible proteins
highly active anti-retroviral therapy
- HIV IPs:
HIV aspartic peptidase inhibitors
ligand preparation module
neglected tropical diseases
Protein Data Bank
root mean square deviation
Rassi A Jr, Rassi A, Marcondes de Rezende J. American trypanosomiasis (Chagas disease). Infect Dis Clin N Am. 2012;26:275–91.
Almeida EA, Ramos Júnior AN, Correia D, Shikanai-Yasuda MA. Co-infection Trypanosoma cruzi/HIV: systematic review (1980–2010). Rev Soc Bras Med Trop. 2011;44:762–70.
Santos AL, d’Avila-Levy CM, Kneipp LF, Sodré CL, Sangenito LS, Branquinha MH. The widespread anti-protozoal action of HIV aspartic peptidase inhibitors: focus on Plasmodium spp., Leishmania spp. and Trypanosoma cruzi. Curr Top Med Chem. 2017;17:1303–17.
Bern C. Chagas disease in the immunosuppressed host. Curr Opin Infect Dis. 2012;25:450–7.
Coura JR. Chagas disease: control, elimination and eradication. Is it possible? Mem Inst Oswaldo Cruz. 2013;108:962–7.
Salomon CJ. First century of Chagas’ disease: an overview on novel approaches to nifurtimox and benzonidazole delivery systems. J Pharm Sci. 2012;101:888–94.
Torreele E, Bourdin Trunz B, Tweats D, Kaiser M, Brun R, Mazué G, Bray MA, Pecoul B. Fexinidazole—a new oral nitroimidazole drug candidate entering clinical development for the treatment of sleeping sickness. PLoS Negl Trop Dis. 2010;4:e923.
Urbina JA. Ergosterol biosynthesis and drug development for Chagas disease. Mem Inst Oswaldo Cruz. 2009;104(Suppl 1):311–8.
Bahia MT, De Andrade IM, Martins TA, Do Nascimento AF, Diniz LF, Caldas IS, Talvani A, Trunz BB, Torreele E, Ribeiro I. Fexinidazole: a potential new drug candidate for Chagas disease. PLoS Negl Trop Dis. 2012;6:e1870.
Ekins S, Williams AJ, Krasowski MD, Freundlich JS. In silico repositioning of approved drugs for rare and neglected diseases. Drug Discov Today. 2011;16:298–310.
Ennes-Vidal V, Menna-Barreto RF, Branquinha MH, Dos Santos AL, d’Avila-Levy CM. Why calpain inhibitors are interesting leading compounds to search for new therapeutic options to treat leishmaniasis? Parasitology. 2017;144:117–23.
Sangenito LS, d’Avila-Levy CM, Branquinha MH, Santos ALS. Nelfinavir and lopinavir impair Trypanosoma cruzi trypomastigote infection in mammalian host cells and show anti-amastigote activity. Int J Antimicrob Agents. 2016;48:703–11.
Sangenito LS, Gonçalves DS, Seabra SH, d’Avila-Levy CM, Santos AL, Branquinha MH. HIV aspartic peptidase inhibitors are effective drugs against the trypomastigote form of the human pathogen Trypanosoma cruzi. Int J Antimicrob Agents. 2016;48:440–4.
Sangenito LS, Menna-Barreto RF, Avila-Levy CM D, Santos AL, Branquinha MH. Decoding the anti-Trypanosoma cruzi action of HIV peptidase inhibitors using epimastigotes as a model. PLoS ONE. 2014;9:e113957.
Bellera CL, Balcazar DE, Vanrell MC, Casassa AF, Palestro PH, Gavernet L, Labriola CA, Gálvez J, Bruno-Blanch LE, Romano PS, Carrillo C, Talevi A. Computer-guided drug repurposing: identification of trypanocidal activity of clofazimine, benidipine and saquinavir. Eur J Med Chem. 2015;93:338–48.
Sangenito LS, de Guedes AA, Gonçalves DS, Seabra SH, d’Avila-Levy CM, Santos ALS, Branquinha MH. Deciphering the effects of nelfinavir and lopinavir on epimastigote forms of Trypanosoma cruzi. Parasitol Int. 2017;66:529–36.
Sangenito LS, Menna-Barreto RFS, Oliveira AC, d’Avila-Levy CM, Branquinha MH, Santos ALS. Primary evidence of the mechanisms of action of HIV aspartyl peptidase inhibitors on Trypanosoma cruzi trypomastigote forms. Int J Antimicrob Agents. 2018;52:185–94.
Umezawa H, Aoyagi T, Morishima H, Matsuzaki M, Hamada M. Pepstatin, a new pepsin inhibitor produced by actinomycetes. J Antibiot (Tokyo). 1970;23:259–62.
Sali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234:779–815.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25:1605–12.
Zhang Y, Skolnick J. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005;33:2302–9.
Friesner RA, Murphy RB, Repasky MP, Frye LL, Greenwood JR, Halgren TA, Sanschagrin PC, Mainz DT. Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J Med Chem. 2006;49:6177–96.
Ido E, Han HP, Kezdy FJ, Tang J. Kinetic studies of human immunodeficiency virus type 1 protease and its active-site hydrogen bond mutant A28S. J Biol Chem. 1991;266:24359–66.
Feng X, Shenker S, Svoboda M, Kožíšek M, Konvalinka J, Gehring K. Structural studies of the yeast DNA damage-inducible protein Ddi1 reveal domain architecture of this eukaryotic protein family. Sci Rep. 2016;6:33671.
Tyagi SC. Inhibitors of human immunodeficiency virus-1 protease. Biochem Cell Biol. 1992;70:309–15.
Santos ALS. HIV aspartyl protease inhibitors as promising compounds against Candida albicans. World J Biol Chem. 2010;1:21–30.
Santos LO, Vitório BS, Branquinha MH, Pedroso e Silva CM, Santos AL, Avila-Levy CM. Nelfinavir is effective in inhibiting the multiplication and aspartic peptidase activity of Leishmania species, including strains obtained from HIV-positive patients. J Antimicrob Chemother. 2013;68:348–53.
Rebello KM, Andrade-Neto VV, Zuma AA, Motta MCM, Gomes CRB, de Souza MVN, Atella GC, Branquinha MH, Santos ALS, Torres-Santos EC, d’Avila-Levy CM. Lopinavir, an HIV-1 peptidase inhibitor, induces alteration on the lipid metabolism of Leishmania amazonensis promastigotes. Parasitology. 2018;145:1304–10.
Conceived the experiments: MHB, ALS, RFSB, ACG, CMD; Performed the experiments: KCSG, VVSC, KMR, LPRB, LSS; Analyzed the data: KCSG, VVSC, LPRB, HLCS, MHB, ALS, RFSB, ACG, CMD; Wrote the manuscript: KCSG, VVSC, LPRB, HLCS, MHB, ALS, RFSB, ACG, CMD. All authors read and approved the final manuscript.
The authors are grateful to the research team leaded by Dr. Jaime M Santana (Universidade de Brasilia) for the extensive efforts on the cloning and expression of T. cruzi aspartyl peptidase, namely: Flavia N. Motta, Isabela M. Dourado and Carla N. de Araújo. The research team leaded by Dr. Philip Noel Suffys (Instituto Oswaldo Cruz, Fiocruz) for the efforts on peptidase expression conditions variation and peptidase purification, namely Dr. Harrison Magdinier Gomes, Raquel L. F. Teixeira and Dr. Márcia Quinhones Pires Lopes. We are also grateful to Dr. Laurent Dardenne (Laboratório Nacional de Computação Científica) for kindly allowing access to the software.
The authors declare that they have no competing interests.
Availability of data and materials
Data are presented in the main manuscript.
Consent for publication
Ethics approval and consent to participate
This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)—Finance Code 001, Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Fundação de Amparo a Pesquisa do Rio de Janeiro (FAPERJ) and Fundação Oswaldo Cruz (Fiocruz). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Intermolecular interactions obtained from the molecular docking pose of the assayed compounds and HIV aspartyl peptidase. We selected three compounds with the highest hits from the molecular docking simulation: saquinavir, atazanavir, and nefilnavir (A–C). In A, the picture on the square represents the superposition of co-crystal (green) and redocked saquinavir structures (orange).
About this article
Cite this article
Castilho, V.V.S., Gonçalves, K.C.S., Rebello, K.M. et al. Docking simulation between HIV peptidase inhibitors and Trypanosoma cruzi aspartyl peptidase. BMC Res Notes 11, 825 (2018). https://doi.org/10.1186/s13104-018-3927-z
- Aspartic peptidase
- Chagas’ disease
- Neglected tropical diseases