The Langevin equation is a stochastic differential equation frequently used in molecular dynamics for simulating systems with a constant temperature. Recent developments have given rise to wide uses of Langevin dynamics at different levels of spatial resolution, which necessitate time step and friction parameter choices outside of the range for which many existing temporal discretization methods were originally developed. We first study the GJ--F, BAOAB and BBK numerical algorithms, originally developed for atomistic simulations, on a coarse-grained polymer melt, paying close attention to the large time step regime. The results of this study then inspire our search for new algorithms and lead to a general class of velocity Verlet-based time-stepping schemes designed to perform well for all parameter regions, by ensuring that they faithfully reproduce statistical quantities for the case of a free particle and harmonic oscillator. This family of methods depends on the choice of a single free parameter function and we explore some of the methods defined for certain choices of this parameter on realistic coarse-grained and atomistic molecular systems relevant in material and bio-molecular science. In addition, we provide an equivalent splitting formulation of this one-parameter family which allows for enhanced insight into the hidden time scaling induced by the choice of the free parameter in the Hamiltonian and stochastic time scales. / Mathematics