We study the geodesic motion of test particles in the space-time of non-compact boson stars. These objects are made of a self-interacting scalar field and -- depending on the scalar field's mass -- can be as dense as neutron stars or even black holes. In contrast to the former these objects do not contain a well-defined surface, while in contrast to the latter the space-time of boson stars is globally regular, can -- however -- only be given numerically. Hence, the geodesic equation also has to be studied numerically. We discuss the possible orbits for massive and massless test particles and classify them according to the particle's energy and angular momentum. The space-time of a boson star approaches the Schwarzschild space-time asymptotically, however deviates strongly from it close to the center of the star. As a consequence, we find additional bound orbits of massive test particles close to the center of the star that are not present in the Schwarzschild case. Our results can be used to make predictions about extreme-mass-ratio inspirals (EMRIs) and we hence compare our results to recent observational data of the stars orbiting Sagittarius A* - the radiosource at the center of our own galaxy.