Conformation and hydrogen ion titration of proteins: a continuum electrostatic model with conformational flexibility

A new method for including local conformational flexibility in calculations of the hydrogen ion titration of proteins using macroscopic electrostatic models is presented. Intrinsic pKa values and electrostatic interactions between titrating sites are calculated from an ensemble of conformers in which the positions of titrating side chains are systematically varied. The method is applied to the Asp, Glu, and Tyr residues of hen lysozyme. The effects of different minimization and/or sampling protocols for both single-conformer and multi-conformer calculations are studied. For single-conformer calculations it is found that the results are sensitive to the choice of all-hydrogen versus polar-hydrogen-only atomic models and to the minimization protocol chosen. The best overall agreement of single-conformer calculations with experiment is obtained with an all-hydrogen model and either a two-step minimization process or minimization using a high dielectric constant. Multi-conformational calculations give significantly improved agreement with experiment, slightly smaller shifts between model compound pKa values and calculated intrinsic pKa values, and reduced sensitivity of the intrinsic pKa calculations to the initial details of the structure compared to single-conformer calculations. The extent of these improvements depends on the type of minimization used during the generation of conformers, with more extensive minimization giving greater improvements. The ordering of the titrations of the active-site residues, Glu-35 and Asp-52, is particularly sensitive to the minimization and sampling protocols used. The balance of strong site-site interactions in the active site suggests a need for including site-site conformational correlations.