We present a formulation of quantized Hamiltonian dynamics (QHD) using variables that arise naturally from the Heisenberg equation of motion. The QHD equations are obtained and solved either directly in terms of these generalized variables, or by employing a wavefunction ansatz. The approach avoids a Taylor expansion and other approximations to the potential, leading to more stable dynamics and a higher precision of the calculated quantities. The proposed formulation is also amenable to for analytic and numerical implementations, thus facilitating its use in molecular dynamics simulation.