which might not be a good approach. Several computational studies have shown that the pKa of the aspartyl dyad changes when inhibitors are present in the binding pocket . Constant pH molecular dynamics might be able to capture the effect of ligand binding on the protonation states of the receptors, but is computationally expensive.