# Best Practices for Simulating Small Ligand-Protein Complexes in GROMACS?

**URL:** <https://gromacs.bioexcel.eu/t/best-practices-for-simulating-small-ligand-protein-complexes-in-gromacs/11793>\
**Category:** User discussions\
**Created:** [March 23, 2025, 10:24pm UTC](https://gromacs.bioexcel.eu/t/best-practices-for-simulating-small-ligand-protein-complexes-in-gromacs/11793 "2025-03-23T22:24:04Z")\
**Posts on this page:** 1\
**Page:** 1

<div class="post-metadata">

**Author:** ![Zepry](https://avatars.discourse-cdn.com/v4/letter/z/f9ae1b/32.png) [@Zepry](https://gromacs.bioexcel.eu/u/Zepry)\
**Post date:** [March 23, 2025, 10:24pm UTC](https://gromacs.bioexcel.eu/t/best-practices-for-simulating-small-ligand-protein-complexes-in-gromacs/11793/1 "2025-03-23T22:24:04Z")

</div>

Hi everyone,

I’m relatively new to GROMACS and I’ve been exploring its capabilities for simulating ligand-protein systems. I’ve mostly worked on protein-only simulations so far, but now I want to move into protein-ligand complexes (small molecules docked into active sites). I’m wondering:

- What are the best practices for setting up such systems in GROMACS?
- Are there any specific force fields or tools recommended for generating ligand topologies (besides PRODRG, which seems a bit outdated)?
- How do people usually validate or check the stability of such complexes over long simulations?

Would really appreciate any insights or shared experiences—especially any gotchas I should be aware of during preprocessing or production runs. Thanks in advance!