A detailed protocol of molecular dynamics simulations of supercoiled DNA molecules that can be in addition knotted or catenated is described. We also describe how to model ongoing action of DNA gyrase that introduces negative supercoing into DNA molecules. The protocols provide detailed instructions about model parameters, equations of used potentials, simulation, and visualization. Implementation of the model into a frequently used molecular dynamics simulation environment, ESPResSo, is shown step by step.