Brownian dynamics simulations of shear flows are carried out for various suspensions of ellipsoids interacting via the Gay-Berne potential. In this simulation all the systems of the suspension are in a liquid crystalline phase at rest. In a shear flow they exhibit various motions of the director depending on the shear rate: the continuous rotation, the intermittent rotation, the wagging-like oscillation, and the aligning. The director is almost always out of the vorticity plane when it rotates, that is the kayaking. The number density of the system and the inter-particle potential intensity significantly affect the shear rate dependence of orientation. In particular, the continuous rotation of director is maintained to higher shear rates for the system with a stronger potential. Furthermore, the rheological properties are examined. The shear-thinning in viscosity is observed, but the negative first normal difference is not obtained.