In this work, we have performed a combination of MD and well-tempered metadynamics simulations on four representative PETases with different optimal temperatures and pH for hydrolytic activity, with the goal of providing an atomistic explanation on why the S/I substitution of two residues in the vicinity of the active site open cleft enhances the activity.