Methodological aspects of QM/MM calculations: A case study on matrix metalloproteinase-2

We address methodological issues in quantum mechanics/molecular mechanics (QM/MM) calculations on a zinc-dependent enzyme. We focus on the first stage of peptide bond cleavage by matrix metalloproteinase-2 (MMP-2), that is, the nucleophilic attack of the zinc-coordinating water molecule on the carbonyl carbon atom of the scissile fragment of the substrate. This step is accompanied by significant charge redistribution around the zinc cation, bond cleavage, and bond formation. We vary the size and initial geometry of the model system as well as the computational protocol to demonstrate the influence of these choices on the results obtained. We present QM/MM potential energy profiles for a set of snapshots randomly selected from QM/MM-based molecular dynamics simulations and analyze the differences in the computed profiles in structural terms. Since the substrate in MMP-2 is located on the protein surface, we investigate the influence of the thickness of the water layer around the enzyme on the QM/MM energy profile. Thin water layers (0-2 Å) give unrealistic results because of structural reorganizations in the active-site region at the protein surface. A 12 Å water layer appears to be sufficient to capture the effect of the solvent; the corresponding QM/MM energy profile is very close to that obtained from QM/MM/SMBP calculations using the solvent macromolecular boundary potential (SMBP). We apply the optimized computational protocol to explain the origin of the different catalytic activity of the Glu116Asp mutant: the energy barrier for the first step is higher, which is rationalized on structural grounds. © 2016 Wiley Periodicals, Inc.