[gmx-users] How to proceed with molecule.itp file

[Cara Kreck]

Hi Sotirios,

You need to either edit the .top file produced by pdb2gmx or make your own from scratch. You only need the pdb2gmx output if your system includes standard amino acids.

If you have both a standard protein and a non-standard molecule, you should make a .pdb file with just the protein to use with pdb2gmx. I would then recommend you remove all the .top specific parts of this file and convert it to a .itp file for just the protein. Then make your own .top file using the bits you removed and link everything together. Alternatively you can just add to the include and molecule sections at the end of pdb2gmx's .top file.

If you don't have any protein, don't bother with pdb2gmx at all.

Here's an example DIY .top file:

; Forcefield files
#include "gromos53a6.ff/forcefield.itp"
#include "./lipids_vdw.itp"

; Parameter set for DOPC
#include "./dopc.itp"

; Include DMSO topology
#include "./dmso.itp"

; Include water topology
#include "spc.itp"

[ system ]
; Name
128 DOPC in DMSO solution

[ molecules ]
; name     number
DOPC          128
DMSO          496
SOL          7573


Make sure you use the molecule names defined in your .itp files and list the molecules in the same order that they appear in your coordinate file. 

Hope that helps,

Cara


> From: si.papadatos at edu.cut.ac.cy
> To: gromacs.org_gmx-users at maillist.sys.kth.se
> Date: Thu, 21 May 2015 07:13:00 +0000
> Subject: [gmx-users] How to proceed with molecule.itp file
> 
> Hi,
> 
> I have made a molecule.itp file according to the chapter 5 of Gromacs Manual and existing entries in the ff directory I want to use.
> 
> My problem is that I don't know how to continue, in theory I have to include it in the .top file but a .top file is produced by pdb2gmx which I cannot use because there isn't an entry in .rtp.
> 
> Please help
> -- 
> Gromacs Users mailing list
> 
> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
> 
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> 
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.