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.