pKa calculations based on advanced MD simulations