Determination of pK(a) values of titrating residues in proteins provides a direct means of studying electrostatic coupling as well as pH-dependent stability. The B1 domain of protein G provides an excellent model system for such investigations. In this work, we analyze the observed pK(a) values of all carboxyl groups in a variant of PGB1 (T2Q, N8D, N37D) at low and high ionic strength as determined using (1)H-(13)C heteronuclear NMR in a structural context. The pK(a) values are used to calculate the pH-dependent stability in low and high salt and to investigate electrostatic coupling in the system. The observed pK(a) values can explain the pH dependence of protein stability but require pK(a) shifts relative to model values in the unfolded state, consistent with persistent residual structure in the denatured state. In particular, we find that most of the deviations from the expected random coil values can be explained by a significantly upshifted pK(a) value. We show also that (13)C backbone carbonyl data can be used to study electrostatic coupling in proteins and provide specific information on hydrogen bonding and electrostatic potential at nontitrating sites.