Accurate charge representation is essential for correctly modeling the interaction between the enzyme and substrate, particularly in systems where electrostatics play a pivotal role in reaction mechanisms. As we mentioned before, for technical reasons, we chose the Mulliken charge model, but depending on the future applications of the semiempirical Hamiltonian, more accurate charge models or electronic density-derived properties maybe implemented.