How similar are enzyme active site geometries derived from quantum mechanical theozymes to crystal structures of enzyme‐inhibitor complexes? Implications for enzyme design

Quantum mechanical optimizations of theoretical enzymes (theozymes), which are predicted catalytic arrays of biological functionalities stabilizing a transition state, have been carried out for a set of nine diverse enzyme active sites. For each enzyme, the theozyme for the rate‐determining transition state plus the catalytic groups modeled by side‐chain mimics was optimized using B3LYP/6–31G(d) or, in one case, HF/3–21G(d) quantum mechanical calculations. To determine if the theozyme can reproduce the natural evolutionary catalytic geometry, the positions of optimized catalytic atoms, i.e., covalent, partial covalent, or stabilizing interactions with transition state atoms, are compared to the positions of the atoms in the X‐ray crystal structure with a bound inhibitor. These structure comparisons are contrasted to computed substrate–active site structures surrounded by the same theozyme residues. The theozyme/transition structure is shown to predict geometries of active sites with an average RMSD of 0.64 Å from the crystal structure, while the RMSD for the bound intermediate complexes are significantly higher at 1.42 Å. The implications for computational enzyme design are discussed.