GROMACS simulation of dimers using the Amber FF14SB forcefield

GROMACS version: 2018.4
GROMACS modification: No

Hi there,

I am attempting to perform MD simulations on a functional dimer Using the Amber FF14SB forcefield set.
The current workflow has me importing the PDB file containing the dimer to amber, generating the box, solvation and ion topology and perameters using tleap, then converting the resultant files to GROMACS .top and .gro formats using parmed. As this process dosent generate positional restraints in the same way that pdb2gmx does, I then generate positional restraints using gmx genrestr.

The issue I am encountering is in the way the topology file is generated in this method. It essentially lists the topology for the monomer in this case termed ‘system1’ and then simply lists ‘system 1’ twice in the [ molecules ] section of the topology file. This becomes an issue as I wish to be able to restrain one half of the dimer pair in order to perform umbrella sampling, as the atom numbering in any POSRE.itp file, generated using gmx genrestr from index groups, wont match the numbering in the general topology file as it essentially only defines topolgy for half the atoms and then duplicates it.

Any recomendations on how to get around it? I havent been able to idenitfy if this is due to the way that tleap or parmed treats protein chains? I am of course happy to provide any of the files I am working with.


An option is to generate manually posre.itp file, putting the desired atoms number. Below an example of position restrains file

; In this topology include file, you will find position restraint
; entries for all the heavy atoms in your original pdb file.
; This means that all the protons which were added by pdb2gmx are
; not restrained.

[ position_restraints ]
; atom type fx fy fz
1 1 1000 1000 1000
5 1 1000 1000 1000
7 1 1000 1000 1000

Best regards

Hi, Alessandra

Thanks this did the job!

All the best,