We generalize the frequently used model of isotropic reorientational Brownian motion via small but finite angular steps to the case in which molecular reorientations by different angles take place around arbitrary axes. A simple approximation allows to give analytical expressions for experimentally relevant correlation functions, including spin lattice relaxation rates. We compare the results of our model to earlier treatments of rotational Brownian motion.