[gmx-users] renumbering residue numbers for version 4.0.2

[Nash, Anthony]

Thanks for that, I've renumbered the residues in the .gro file, and even so far as in the respective protein .itp files. However, grompp is still throwing out the same error. Not sure where it is holding onto the number 652 (which is the start amino acid index of my protein) as I have reset the residue numbers to begin from 1. 

Residue numbers in the .top are not numbered consecutively from 1 (rather, resnumber = 652 and resnr_diff = -1)

My top file is:

#include "forcefield_popc.itp"
#include "mut_helixA_renumber.itp"
#include "mut_helixB_renumber.itp"

#include "popc.itp"
#ifdef FLEX_SPC
#include "flexspc.itp"
#else
#include "spc.itp"
#endif
#include "ions.itp"

[ system ]
Motif A oncogene in water

[ molecules ]
helixA          1
helixB          1
POPC            395
SOL         23765
CL              10

The grompp command is:
~/gromacs-localpressure/bin/grompp_d -c starting_config.gro -n lateral.ndx -p ../../mut_system_renumber.top -o local_pressure -f local_pressure.mdp

Any thoughts?
Many thanks

________________________________________
From: gmx-users-bounces at gromacs.org [gmx-users-bounces at gromacs.org] on behalf of Mark Abraham [mark.j.abraham at gmail.com]
Sent: 29 March 2013 14:17
To: Discussion list for GROMACS users
Subject: Re: [gmx-users] renumbering residue numbers for version 4.0.2

On Fri, Mar 29, 2013 at 2:10 PM, Nash, Anthony
<Anthony.Nash at warwick.ac.uk>wrote:

> Apologies if this has been covered once before, I've had a hearty search
> and I've yet to find anything.
>
> I have just installed a modified verison of Gromacs 4.0.2 which calculate
> lateral pressure. However, when I run grommp I get:
>
> Residue numbers in the .top are not numbered consecutively from 1 (rather,
> resnumber = 652 and resnr_diff = -1)
>
> Yet, running grompp on the same files with version 4.5.1 everything is
> perfectly fine. I am guessing somewhere in between the two versions
> something was changed.
>
> So, does this sound about right? Does the order of residue numbers no
> longer mater with the latest versions of grompp? And ultimately, does
> anyone know how I can go through a .gro file and renumber every residue
> (short of writing a perl script)?
>

Use genconf -renumber from probably any version (and curse whoever thought
that was a good place to put it :-))

Mark
--
gmx-users mailing list    gmx-users at gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists