Most of the simulations are carried out in neutral pH. http://www.gromacs.org/Documentation/How-tos/Constant_pH_Simulation you can refer this paper which for different protonation state. http://pubs.acs.org/doi/abs/10.1021/ct401042b