A potential slice of SOD passing through the two binding sites. The linearized and nonlinear Poisson-Boltzmann models yield substantially different electrostatic potential values near the binding sites, leading to correspondingly different reaction rates predicted by Brownian dynamics simulations. This implies that the full nonlinear model is important for certain modeling situations.