[gmx-users] AMBER and gromacs atoms renaming
Justin A. Lemkul
jalemkul at vt.edu
Sun Aug 10 18:56:47 CEST 2008
Ragnarok sdf wrote:
> Hi, I am running a simulation of a protein DNA complex, and for this i
> am using AMBER99 forcefield. The problem is that every complex I need
> to run, I have to rename several residues and atoms, especially
> hydrogen atoms one by one (well, actually using replace tool). And
> this takes hours to get done.
You could create (in a sense) an "alias" for the problematic residues. That is,
the AMBER equivalent of LYS is LYN, and for LYSH is LYP. What I've done with my
.rtp files is to copy the section that designates LYN and make a duplicate entry
that is headed [ LYS ]. Similarly, you can apply this type of manipulation to
ASP, GLU, CYS, HIS, and anything else you find that has different nomenclature.
You would then have to make corresponding changes in the .hdb file, as well.
Using pdb2gmx -ignh in conjunction with the .hdb file will eliminate the need to
manually re-name all hydrogen atoms. If you do this, the only names you will
have to change are the N- and C-termini.
If you don't want to do all that, find and replace within vi is quite simple,
will find and replace all LYS with LYP. If you need varying protonation states,
this can be more tedious. Otherwise, altering all the residue names should only
take a few minutes per file.
Another option is to use xleap to prepare your input. It is now freely
available as part of the AmberTools package. Using xleap should give all your
atoms the proper AMBER nomenclature.
> Is there any shortcut I can use, like a software that prepares my pdb
> files for an AMBER run, or something close to that?
> If not, Ill just keep doing things by hand.
> Thank you
> Fabrício Bracht
> gmx-users mailing list gmx-users at gromacs.org
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
jalemkul[at]vt.edu | (540) 231-9080
More information about the gromacs.org_gmx-users