An algorithm is developed to determine the electrophoretic mobility of a rigid polyion modeled as a low dielectric volume element of arbitrary shape containing an arbitrary charge distribution. The solvent is modeled as a high dielectric continuum with salt distributed according to the linearized Poisson Boltzmann equation. Account is also taken of a Stern layer that separates the molecular surface and the surface of hydrodynamic shear, or Stern surface. Relaxation of the ion atmosphere because of the presence of the external field is ignored. The electrostatic and hydrodynamic problems are both solved by boundary element methods, The procedure is first applied to spherical polyions containing monopolar, dipolar, and quadrupolar charge distributions, and calculated mobilities are found to be in excellent agreement with the theory of Yoon and Kim. It is then applied to lysozyme by using models that account for the detailed shape and charge distribution of the enzyme. For reasonable choices of the molecular and Stern surfaces, calculated and experimental mobilities are found to be in fair agreement with each other, However, if a pH independent Stern layer (or, equivalently, translational diffusion constant, D-t) is assumed, the calculated mobilities exhibit a stronger pH dependence than is observed experimentally. A small increase in D-t with increasing pH could correct this discrepancy.