A method for producing enzyme mutants includes a generative model parameterized using homologous sequences of wild-type enzymes. Based on a distance between a mutated amino acid residue of each mutant and a corresponding substrate when the substrate is bound to the mutant being less than or equal to a first threshold, the processor may select a first candidate set of enzyme mutants and calculate a statistical energy of each mutant of the first candidate set, and identify, based on the calculated statistical energies, a first set of enzyme mutants among the first candidate set to improve enzyme efficiency. Based on the distance being larger than or equal to a second threshold, the processor may select a second candidate set of enzyme mutants and calculate a statistical energy of each mutant of the second candidate set, and identify, based on the calculated statistical energies, a second set of enzyme mutants among the second candidate set of enzyme mutants to improve protein stability.