  An interfacial regularized Stokeslet scheme is presented to predict the
motion of solid bodies (e.g. proteins or gel-phase domains) embedded within
flowing lipid bilayer membranes. The approach provides a numerical route to
calculate velocities and angular velocities in complex flow fields that are not
amenable to simple Fax\'en-like approximations. Additionally, when applied to
shearing motions, the calculations yield predictions for the effective surface
viscosity of dilute rigid body-laden membranes. In the case of cylindrical
proteins, effective viscosity calculations are compared to two prior analytical
predictions from the literature. Effective viscosity predictions for a dilute
suspension of rod-shaped objects in the membrane are also presented.
