After the construction of the systems, the potential energy of each system was minimized using the steepest descent algorithm, followed by 10 ps NVT equilibration and 10 ps NPT equilibration. During these simulations, the λ-coordinates of all titratable groups were kept constant. Then, a 10 ns constant pH simulation was performed for cardiotoxin V from a