[gmx-users] choosing appropriate force field
elham tazikeh
elham.tazikeh at gmail.com
Thu Jan 1 11:13:47 CET 2015
Dear Justin
thanks for your help
i used of *ATB server* for making *.itp* file for SAL (salicylic acid), as
below:
[ moleculetype ]
; Name nrexcl
_N4Y 3
[ atoms ]
; nr type resnr resid atom cgnr charge mass total_charge
1 OA 1 _N4Y O10 1 -0.546 15.9994
2 H 1 _N4Y H16 1 0.469 1.0080
3 C 1 _N4Y C8 1 0.643 12.0110
4 O 1 _N4Y O9 1 -0.566 15.9994 ; 0.000
5 C 1 _N4Y C3 2 -0.221 12.0110
6 C 1 _N4Y C2 2 0.388 12.0110
7 OA 1 _N4Y O7 2 -0.604 15.9994
8 H 1 _N4Y H15 2 0.437 1.0080 ; 0.000
9 C 1 _N4Y C1 3 -0.266 12.0110
10 HC 1 _N4Y H11 3 0.175 1.0080
11 C 1 _N4Y C6 3 -0.049 12.0110
12 HC 1 _N4Y H14 3 0.140 1.0080 ; 0.000
13 C 1 _N4Y C5 4 -0.201 12.0110
14 HC 1 _N4Y H13 4 0.141 1.0080
15 C 1 _N4Y C4 4 -0.073 12.0110
16 HC 1 _N4Y H12 4 0.133 1.0080 ; 0.000
; total charge of the molecule: 0.000
here, i have a question:
in ATB server , there are 2 force field as default:* gromos 54A7 , gromos
53A6*
i did not see other force field,right?
for my simulation, can i choose *gromos43a1* or i have to choosing the same
force field (as ATB), i mean i must choose gromos 54A7?
on theother hand, Formats in ATB,were:
Gromacs
Gromos 11
Gromos 96
for my study, i choosed *Gromacs* as format in ATB , is it correct?
finally, for SAL molecule ; there is SAL in ATB as default, and i used of
it without loading of anywhere (same PRODRG)
is this correct?
at the end, shoul i change *pdb* format (produced by ATB) to *gro* with
*editconf* tool in gromacs and then *copy* et the end of my protein
coordinate file?
or
change my *complex.pdb* (get from RCSB) format to *coordinate.gro?*
i really appriciate for your advise
cheers
More information about the gromacs.org_gmx-users
mailing list