“…Self-consistent calculation of pK a values of surface and buried groups using a single value of ⑀ in is impossible with semi-macroscopic methods unless protein flexibility is taken into account explicitly. When conformational relaxation is treated explicitly, low ⑀ in can yield reasonable results (Langsetmo et al, 1991;Antosiewicz et al, 1994Antosiewicz et al, , 1996You and Bashford, 1995;Alexov and Gunner, 1997;Zhou and Vijayakumar, 1997;Rabenstein et al, 1998;van Vlijmen et al, 1998;Havranek and Harbury, 1999;Scharnagl et al, 1999;Ullmann and Knapp, 1999), especially when the neutral and the ionized states of the group of interest are treated separately to capture contributions by the relaxation of local dipoles upon ionization (Langen et al, 1992;Yang et al, 1993;Alexov and Gunner, 1997;Sham et al, 1997;Zhou and Vijayakumar, 1997;Rabenstein and Knapp, 2001). In microscopic methods, such as the protein dipole-Langevin dipole (PDLD) method developed by Warshel and co-workers, the relaxation of protein dipoles has always been included explicitly (Warshel and Russell, 1984).…”