The major cysteine protease of Trypanosoma cruzi, cruzain (CRZ), has been described as a therapeutic target for Chagas' disease, which affects millions of people worldwide. Thus, a series of CRZ inhibitors has been studied, including a new competitive inhibitor, Nequimed176 (NEQ176). Nevertheless, the structural and dynamic basis for CRZ inhibition remains unclear. Hoping to contribute to this ever-growing understanding of timescale dynamics in the CRZ inhibition mechanism, we have performed the first study using 100 ns of molecular dynamics (MD) simulations of two CRZ systems in an aqueous solvent under pH 5.5: CRZ in the apo form (ligand free) and CRZ complexed to NEQ176. According to the MD simulations, the enzyme adopts an open conformation in the apo form and a closed conformation in the NEQ176-CRZ complex. We also suggest that this closed conformation is related to the hydrogen-bonding interactions between NEQ176 and CRZ, which occurs through key residues, mainly Gly66, Met68, Asn69, and Leu160. In addition, the cross-correlation analysis shows evidence of the correlated motions among Ala110-Asp140, Leu160-Gly189, and Glu190-Gly215 subdomains, as well as, the movements related to Ala1-Thr59 and Asp60-Pro90 regions seem to be crucial for CRZ activity.
Keywords: Chagas’ disease; Nequimed176; cruzain; cysteine protease inhibitors; molecular dynamics simulation.