A broad array of canonical sampling methods are available for molecular simulation based on stochastic-dynamical perturbation of Newtonian dynamics, including Langevin dynamics, Stochastic Velocity Rescaling, and methods that combine Nos\'{e}-Hoover dynamics with stochastic perturbation. In this article we discuss several stochastic-dynamical thermostats in the setting of simulating systems with holonomic constraints. The approaches described are easily implemented and facilitate the recovery of correct canonical averages with minimal disturbance of the underlying dynamics. For the purpose of illustrating our results, we examine the numerical application of these methods to a simple atomic chain, where a Fixman term is required to correct the thermodynamic ensemble.