Dear Justin Thanks for your answer. I want to use GROMOS96 43a1 force field for my protein and PRODRG parameters for my ligand. I want to use leap-frog algorithm for integrator. Based on your opinion, I don't use QM charges. In these conditions, I will use my first way mentioned in my previous email. Is my manner true?