Dear all,
In the pmx tutorial, it is described a workflow for thermodynamic integration of ligand. For the generation of the double topologies a gromacs file *.itp is required: is there a way to write it? I previously parametrized my system with tleap in order to use Amber force field. So my gromacs topology file has all the solute atoms in a single group since I used the Acpype parser to obtain input file for MD simulations. Do I have to select by hand all the ligand parameters or is there a simpler way?
Thanks for the help
I recommend using pdb2gmx to generate topologies for protein. The ligand topologies then can be created with acpype+antechamber combination.
I can generate the protein as well through the Amber tleap module. Instead, with the combination acype+antechamber I can generate a *top and *gro files with just the ligand. How the *itp with just ligand parameters is generated?
.itp file is similar to .top, it just contains entries relevant for the given molecule, e.g. moleculetype, atoms, bonds, angles, dihedrals, etc.
You can find more information in the Gromacs manual: File formats — GROMACS 2021.3 documentation