The effect of hydrodynamic torques on the association rate constants for enzyme-ligand complexation is investigated by Brownian dynamics simulations. Our hydrodynamic models of the enzyme and ligand are composed of spherical elements with friction forces acting at their centers. A quantitative measure of hydrodynamic torque orientational effects is introduced by choosing, as a reference system, an enzyme-ligand model with the same average hydrodynamic interactions but without orientational dependence. Our simple models show a 15% increase in the rate constant caused by hydrodynamic torques at physiological ionic strength. For more realistic hydrodynamic models, which are not computationally feasible at present, this effect is probably higher. The most important finding of this work is that hydrodynamic complementarity in shape (i.e. like the fitting together of pieces of a puzzle) is most effective for interactions between molecules at physiological ionic strength.