A quantum mechanics/molecular mechanics study of the reaction mechanism of the hepatitis C virus NS3 protease with the NS5A/5B substrate.
Combined quantum mechanics and molecular mechanics (QM/MM) calculations were carried out to characterize the reaction mechanism of the NS3 protease with its preferred substrate (NS5A/5B). The main purpose of this study was to locate the barrier states and intermediates along the distinguished coordinate path (DCP) involved in this process. These structures, and in particular the one corresponding to the first barrier state and intermediate (B1 and I1), could be a starting point for the synthesis of inhibitors of this protease, which could be used to treat hepatitis C. The two first steps of the reaction mechanism were studied, i.e., the acylation step and the breaking of the peptide bond. The first step takes place through a tetracoordinated intermediate, as suggested from previous works on other Serine proteases. The importance of the different amino acid residues was also considered (perturbation study where the MM charges of each residue were set to zero independently). The residues of the oxyanion hole were confirmed as the most important for the electrostatic stabilization of the tetracoordinate intermediate. Moreover, the role of other residues, e.g., Arg-155 and Asp-79, was also explained.