Molecular mechanics calculations on Rous sarcoma virus protease with peptide substrates

Molecular models of Rous sarcoma virus (RSV) protease and 20 peptide substrates with single amino acid substitutions at positions from P4 to P3', where the scissile bond is between P1 and P1'. were built and compared with kinetic measurements. The unsubstituted peptide substrate. Pro-Ala-Val-Ser-Leu-Ala-Met-Thr, represents the NC-PR cleavage site of RSV protease. Models were built of two intermediates in the catalytic reaction, RSV protease with peptide substrate and with the tetrahedral intermediate. The energy minimization used an algorithm that increased the speed and eliminated a cutoff for nonbonded interactions. The calculated protease-substrate interaction energies showed correlation with the relative catalytic efficiency of peptide hydrolysis. The calculated interaction energies for the 8 RSV protease-substrate models with changes in P1 to P1' next to the scissile bond gave the highest correlation coefficient of 0.79 with the kinetic measurements, whereas all 20 substrates showed the lower, but still significant correlation of 0.46. Models of the tetrahedral reaction intermediates gave a correlation of 0.72 for the 8 substrates with changes next to the scissile bond, whereas a correlation coefficient of only 0.34 was observed for all 20 substrates. The differences between the energies calculated for the tetrahedral intermediate and the bound peptide gave the most significant correlation coefficients of 0.90 for models with changes in P1 and P1', and 0.56 for all substrates. These results are compared to those from similar calculations on HIV-1 protease and discussed in relation to the rate-limiting steps in the catalytic mechanism and the entropic contributions.