Rather than increasing the timescale of the MD simulations to obtain converged dihedral and protonation state distributions, or introducing enhanced sampling techniques, we propose to reduce the barriers for dihedral rotations in a systematic way. We will demonstrate that such optimized dihedral force field parameters improve pK a estimates of amino acids, without compromising the overall conformational sampling of the protein.