Electrostatics of cysteine residues in proteins: Parameterization and validation of a simple model

One of the most popular and simple models for the calculation of pK_as from a protein structure is the semi‐macroscopic electrostatic model MEAD. This model requires empirical parameters for each residue to calculate pK_as. Analysis of current, widely used empirical parameters for cysteine residues showed that they did not reproduce expected cysteine pK_as; thus, we set out to identify parameters consistent with the CHARMM27 force field that capture both the behavior of typical cysteines in proteins and the behavior of cysteines which have perturbed pK_as. The new parameters were validated in three ways: (1) calculation across a large set of typical cysteines in proteins (where the calculations are expected to reproduce expected ensemble behavior); (2) calculation across a set of perturbed cysteines in proteins (where the calculations are expected to reproduce the shifted ensemble behavior); and (3) comparison to experimentally determined pK_a values (where the calculation should reproduce the pK_a within experimental error). Both the general behavior of cysteines in proteins and the perturbed pK_a in some proteins can be predicted reasonably well using the newly determined empirical parameters within the MEAD model for protein electrostatics. This study provides the first general analysis of the electrostatics of cysteines in proteins, with specific attention paid to capturing both the behavior of typical cysteines in a protein and the behavior of cysteines whose pK_a should be shifted, and validation of force field parameters for cysteine residues. Proteins 2012. © 2012 Wiley Periodicals, Inc.