Any one using OpenMM dynamic simulation should help with the steps and python script for integrating ligand to the protein in PDB.

Also, I find that Charmm parameterization for ligand has no template for NH3 which is in my ligand. Is there any way out?

Thanks