Elucidating Enzymatic Catalysis Using Fast Quantum Chemical Descriptors.

In general, computational simulations of enzymatic catalysis processes are thermodynamic and structural surveys to complement experimental studies, requiring high level computational methods to match accurate energy values. In the present work, we propose the usage of reactivity descriptors, theoretical quantities calculated from the electronic structure, to characterize enzymatic catalysis outlining its reaction profile using low-level computational methods, such as semiempirical Hamiltonians. We simulate three enzymatic reactions paths, one containing two reaction coordinates and without prior computational study performed, and calculate the reactivity descriptors for all obtained structures. We observed that the active site local hardness does not change substantially, even more so for the amino-acid residues that are said to stabilize the reaction structures. This corroborates with the theory that activation energy lowering is caused by the electrostatic environment of the active sites. Also, for the quantities describing the atom electrophilicity and nucleophilicity, we observed abrupt changes along the reaction coordinates, which also shows the enzyme participation as a reactant in the catalyzed reaction. We expect that such electronic structure analysis allows the expedient proposition and/or prediction of new mechanisms, providing chemical characterization of the enzyme active sites, thus hastening the process of transforming the resolved protein three-dimensional structures in catalytic information.