Cysteine proteases play a major role in many life processes and are the target of key drugs. The reaction mechanism of these enzymes is a complex process, which involves several steps that are divided into two main groups: acylation and deacylation. In this work, we studied the energy profile for the acylation and a part of the deacylation reaction of three different enzymes, cruzain, papain, and the Q19A-mutated papain with the benzyloxycarbonyl-phenylalanylarginine-4-methylcoumaryl-7-amide (CBZ-FR-AMC) substrate. The calculations were performed using the EVB and PDLD/S-LRA methods. The overall agreement between the calculated and observed results is encouraging and indicates that we captured the correct reaction mechanism. Finally, our finding indicates that the minimum of the reaction profile, between the acylation and deacylation steps, should provide an excellent state for the binding of covalent inhibitors.