Computation of the electrophoretic mobility of proteins

A scheme is presented for computing the electrophoretic mobility of proteins in free solution, accounting for the details of the protein shape and charge distribution. The method of Teubner is implemented using a boundary integral formulation within which the velocity distribution, the equilibrium electrical potential around the molecule, and the potential distribution due to the applied field are solved for numerically using the boundary element method. Good agreement of the numerical result is obtained for spheres with the corresponding semi-analytical specialization of Henry's analysis. For protein systems, the method is applied to lysozyme and ribonuclease A. In both cases, the predicted mobility tensors are fairly isotropic, with the resulting scalar mobilities being significantly smaller than for spheres of equal volume and net charge. Comparisons with previously published experimental results for ribonuclease show agreement to be excellent in the presence of a net charge, but poorer at the point of zero charge. The approach may be useful for evaluating approximate methods for estimating protein electrophoretic mobilities and for using electrophoretic measurements to obtain insight into charge distributions on proteins.