We propose an integrated MOZYME + DFT method for calculation of the pKa val
ues of ionizable residues in a protein. The method is applied to two protei
n systems, bacteriorhodopsin (bR) and ribonuclease T1. It is shown that the
present method reproduces the experimental pKa values of several targeted
residues better than other conventional methods. Thus, it is available for
pKa calculation for large molecular systems. (C) 2001 Elsevier Science B.V.
All rights reserved.