This protocol allows running classical molecular dynamics on any potential provided by cuby interfaces. Verlet propagator is used and several thermostats are available. It is possible to run the MD only on a part of the system, keeping the rest frozen.