In this work, we examine the viscosity of a dilute suspension of irregularly shaped particles at low shear. A particle is modeled as a rigid array of nonoverlapping beads of variable size and geometry. Starting from a boundary element formalism, approximate account is taken of the variation in hydrodynamic stress over the surface of the individual beads. For a touching dimer of two identical beads, the predicted viscosity is lower than the exact value by 5.2%. The methodology is then applied to several other model systems including tetramers of variable conformation and linear strings of touching beads. An analysis is also carried out of the viscosity and translational diffusion of several dilute amino acids and diglycine in water. It is concluded that continuum hydrodynamic modeling with stick boundary conditions is unable to account for the experimental viscosity and diffusion data simultaneously. A model intermediate between "stick" and "slip" could possibly reconcile theory and experiment.