It is known that protonation states of the receptor and the ligand before and after the binding are a very important factor for the binding, as indicated by both experimental data and computational modeling in literature.
I want to perform a molecular docking for protein-ligand varying pH, that is, changing the protonation states of the protein.
My approach is to use the Chimera software to modify the protonation state of the protein and carry out a docking for each molecule (with and without protonation states). However, I do not know if this approach is complete, and although the importance of the protonation state is highlighted in the literature, I have not found a working protocol for this.