Cruzipains are the main papain-like cysteine proteases of Trypanosoma cruzi, the protozoan parasite that causes Chagas disease. Encoded by a multigenic family, previous studies have estimated the presence of dozens of copies spread over multiple chromosomes in different parasite strains. Here, we describe the complete gene repertoire of cruzipain in three parasite strains, their genomic organization, and expression pattern throughout the parasite life cycle. Furthermore, we have analyzed primary sequence variations among distinct family members as well as structural differences between the main groups of cruzipains. Based on phylogenetic inferences and residue positions crucial for enzyme function and specificity, we propose the classification of cruzipains into two families (I and II), whose genes are distributed in two or three separate clusters in the parasite genome, according with the strain. Family I comprises nearly identical copies to the previously characterized cruzipain 1/cruzain, whereas Family II encompasses three structurally distinct sub-types, named cruzipain 2, cruzipain 3, and cruzipain 4. RNA-seq data derived from the CL Brener strain indicates that Family I genes are mainly expressed by epimastigotes, whereas trypomastigotes mainly express Family II genes. Significant differences in the active sites among the enzyme sub-types were also identified, which may play a role in their substrate selectivity and impact their inhibition by small molecules.

Cruzipains are the main papain-like cysteine proteases of Trypanosoma cruzi, the protozoan parasite that causes Chagas disease. Encoded by a multigenic family, previous studies have estimated the presence of dozens of copies spread over multiple chromosomes in different parasite strains. Here, we describe the complete gene repertoire of cruzipain in three parasite strains, their genomic organization, and expression pattern throughout the parasite life cycle. Furthermore, we have analyzed primary sequence variations among distinct family members as well as structural differences between the main groups of cruzipains. Based on phylogenetic inferences and residue positions crucial for enzyme function and specificity, we propose the classification of cruzipains into two families (I and II), whose genes are distributed in two or three separate clusters in the parasite genome, according with the strain. Family I comprises nearly identical copies to the previously characterized cruzipain 1/cruzain, whereas Family II encompasses three structurally distinct sub-types, named cruzipain 2, cruzipain 3, and cruzipain 4. RNA-seq data derived from the CL Brener strain indicates that Family I genes are mainly expressed by epimastigotes, whereas trypomastigotes mainly express Family II genes. Significant differences in the active sites among the enzyme sub-types were also identified, which may play a role in their substrate selectivity and impact their inhibition by small molecules.

The cysteine protease cruzain is a Chagas disease target, exploited in computational studies. However, there is no consensus on the protonation states of the active site residues Cys25, His162, and Glu208 at the enzyme's active pH range. We evaluated the impact of different protonation states of these residues on docking calculations. Through a retrospective study with cruzain inhibitors and decoys, we compared the performance of virtual screening using four grids, varying protonation states of Cys25, His162, and Glu208. Based on enrichment factors and ROC plots, docking with the four grids affected compound ranking and the overall charge of top-ranking compounds. Different grids can be complementary and synergistic, increasing the odds of finding different ligands with diverse chemical properties.