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: [email protected] [[email protected]] on behalf of Mark Abraham [[email protected]] 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 <[email protected]>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 [email protected] 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 [email protected]. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- gmx-users mailing list [email protected] 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 [email protected]. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
