A new realization of the constant pH molecular dynamics simulation is proposed in mean force potential of proton reservoir in constant titration conditions. The MD-pH-ET method is used for a protein in the most probable ionization microstate, which is the state with minimal energy of the current conformation taking into account the correcting ionization potential of mean force considering an equilibrium ensemble of ionization states. The new MD-pH-ET method allows one to carry out an optimization of protein structure and the total free energy of a protein in the aqueous solution at constant pH and to calculate the pH-dependent properties. The MD-pH-ET method possesses the following unique features: (1) it uses the most precise and computationally effective realization of calculation of electrostatic energy of a protein in aqueous solution, the model of continuous dielectric media with Poisson equation, and the generalized Born method with "ideal" Born atomic radii; (2) it uses the same model of the potential energy surface in the ionization-conformation phase space both for calculating the potential energy of the protein and atomic forces and for determining the most probable ionization states; (3) it calculates the total free energy of the protein in the aqueous solution in a proton reservoir under the conditions of equilibrium titration. The workability of the new method MD-pH-ET is demonstrated for a molecule of the BPTI protein.