Modeling of serine protease prototype reactions with the flexible effective fragment potential quantum mechanical/molecular mechanical method

A complete cycle of chemical transformations for the serine protease prototype reaction is modeled following calculations with the flexible effective fragment quantum mechanical/molecular mechanical (QM/MM) method. The initial molecular model is based on the crystal structure of the trypsin–bovine pancreatic trypsin inhibitor complex including all atoms of the enzyme within approximately 15–18 Å of the oxygen center O γ of the catalytic serine residue. Several selections of the QM/MM partitioning are considered. Fractions of the side chains of the residues from the catalytic triad (serine, histidine and aspartic acid) and a central part of a model substrate around the C–N bond to be cleaved are included into the QM subsystem. The remaining part, or the MM subsystem, is represented by flexible chains of small effective fragments, whose potentials explicitly contribute to the Hamiltonian of the QM part, but the corresponding fragment–fragment interactions are described by the MM force fields. The QM/MM boundaries are extended over the C α –C β bonds of the peptides assigned to the QM subsystem in the enzyme, C–C and C–N bonds in model substrates. Multiple geometry optimizations have been performed by using the RHF/6-31G method in the QM part and OPLSAA or AMBER sets of MM parameters, resulting in a series of stationary points on the complex potential-energy surfaces. All structures generally accepted for the serine protease catalytic cycle have been located. Energies at the stationary points found have been recomputed at the MP2/6-31+G* level for the QM part in the protein environment. Structural changes along the reaction path are analyzed with special attention to hydrogen-bonding networks. In the case of a model substrate selected as a short peptide CH3(NHCO-CH2)2 – HN–CO–(CH2–NHCO)CH3 the computed energy profile for the acylation step shows too high activation energy barriers. The energetics of this rate-limiting step is considerably improved, if more realistic model for the substrate is considered, following the motifs of the ThrI11–GlyI12–ProI13-–CysI14–LysI15–AlaI16–ArgI17–IleI18–IleI19 sequence of the bovine pancreatic trypsin inhibitor.