An efficient simulation method is presented for Brownian fiber suspensions, which includes both uncrossability of the fibers and hydrodynamic interactions between the fibers mediated by a mesoscopic solvent. To conserve hydrodynamics, collisions between the fibers are treated such that momentum and energy are conserved locally. The choice of simulation parameters is rationalized on the basis of dimensionless numbers expressing the relative strength of different physical processes. The method is applied to suspensions of semiflexible fibers with a contour length equal to the persistence length, and a mesh size to contour length ratio ranging from 0.055 to 0.32. For such fibers the effects of hydrodynamic interactions are observable, but relatively small. The noncrossing constraint, on the other hand, is very important and leads to hindered displacements of the fibers, with an effective tube diameter in agreement with recent theoretical predictions. The simulation technique opens the way to study the effect of viscous effects and hydrodynamic interactions in microrheology experiments where the response of an actively driven probe bead in a fiber suspension is measured.