[gmx-users] how to creat new.rtp for acetylated lysine (ALY)

Justin Lemkul jalemkul at vt.edu
Mon Dec 4 15:25:14 CET 2017



On 12/4/17 5:27 AM, Mehreen Jan wrote:
> HI!
> Hope you all will be fine.
> Version : Gromacs5.1.4
> Thank you for your reply. According to your advice I have pasted parameters for acetylated lysine (ALY) from “stream/prot/toppar_all36_prot_modify_res.str (RESI ALY)” in the Charmm36ff (merged.rtp). As you suggested that just create a new .rtp entry based on LYS and add it to residuetypes.dat.  Can we rename merged.rtp with new .rtp? How to create new.rtp? Can we add ALY as protein in   residuetypes.dat?   When we past ALY parameters in merged.rtp gives the following error
>                                                                  :at line ALY error *.
> I have removed these asterisk (*) but it gives the following error:
>                                                      : in .rtp file in residue ALY at line
>                                                            : Atom N NH1 -0.47 !.
> Would you like to please tell me how these parameters should be integrated in any force field or in Charmm36 that

You can't simply copy and paste from a CHARMM RTF into a GROMACS .rtp 
file and expect it to work. The format is entirely different.

Look at standard LYS and see how it is formatted. Copy its entry, modify 
as needed, and keep it in the same .rtp file (no need for more .rtp 
files, just write a new residue entry in the existing file).

-Justin

-- 
==================================================

Justin A. Lemkul, Ph.D.
Assistant Professor
Virginia Tech Department of Biochemistry

303 Engel Hall
340 West Campus Dr.
Blacksburg, VA 24061

jalemkul at vt.edu | (540) 231-3129
http://www.biochem.vt.edu/people/faculty/JustinLemkul.html

==================================================



More information about the gromacs.org_gmx-users mailing list