Thanks for your answer but I guess the links that you have sent is more for the prediction for protein-ligand complex protonation state rather than just the protein. I am not able to get the point whether the affinity would be same for predicting the protonation state correctly for my purpose. Specifically for Histidine protonation state.
Propka uses the PDB2PQR which I already have used and the other server preddimer and might not be suitable for my purpose. proteinplus might be one of the useful for me out of all your suggestions.
Hi Vishal, take a look at the ProteinPrepare web application (www.playmolecule.org). It uses PROPKA3.1 for protonation prediction plus PDB2PQR for H-bond network optimization. You can then download results and optimized structure for simulation.