The velocity fluctuations present in macroscopically homogeneous suspensions of neutrally buoyant, non-Brownian spheres undergoing simple shear flow, and their dependence on the microstructure developed by the suspensions, are investigated in the limit of vanishingly small Reynolds numbers using Stokesian dynamics simulations. We show that, in the dilute limit, the standard deviation of the velocity fluctuations (the so-called suspension temperature) is proportional to the volume fraction, in both the transverse and the flow directions, and that a theoretical prediction, which considers only the hydrodynamic interactions between isolated pairs of spheres, is in good agreement with the numerical results even up to volume fractions of the order of 10%.