A method for writing enzyme rate equations: application to the estimation of the number and size of key proton families