Although it is well-known that limited local mutations of enzymes, such as matrix metalloproteinases (MMPs), may change enzyme activity by orders of magnitude as well as its stability, the completely rational design of proteins is still challenging. These local changes alter the electrostatic potential and thus local electrostatic fields, which impacts the dynamics of water molecules close the protein surface. Here we show by a combined computational design, experimental, and molecular dynamics (MD) study that local mutations have not only a local but also a global effect on the solvent: In the specific case of the matrix metalloprotease MMP14, we found that the nature of local mutations, coupled with surface morphology, have the ability to influence large patches of the water hydrogen-bonding network at the protein surface, which is correlated with stability. The solvent contribution can be experimentally probed via terahertz (THz) spectroscopy, thus opening the door to the exciting perspective of rational protein design in which a systematic tuning of hydration water properties allows manipulation of protein stability and enzymatic activity.