<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
<br>
Hi,<br>
<br>
It is not a pdb2gmx feature, but a global one.<br>
pdb2gmx is only affected in the sense that it retains the residue numbers from the input pdb.<br>
I assume those will be different for most lipid pdb files (if you use pdb2gmx for those).<br>
This renumbering is done by any program that needs to output or select global residue numbers.<br>
Currently this is switchable with the env.var. I mentioned.<br>
If you put this env.var. in your GMXRC, you can get things how you want them.<br>
Adding an option to all programs is not a good idea.<br>
<br>
But I am open for any suggestions on this issue.<br>
<br>
There is something I don't get though.<br>
In the problematic output the lipids can not be one single residue, but should be two or more residues.<br>
<br>
Berk<br>
<br>
> Date: Mon, 18 Jan 2010 10:31:19 -0500<br>
> From: <a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a><br>
> To: <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
> Subject: Re: [gmx-users] trjconv problem<br>
><br>
><br>
> Could renumbering be a switchable feature? For instance, pdb2gmx -[no]renumber,<br>
> with "no" being default? Otherwise, could you provide an example of how to<br>
> properly use the environment variable you posted, since this would seem to<br>
> affect all of us in the membrane protein world quite distinctly (i.e., a<br>
> singly-numbered lipid, as reported, kills many of the programs we have written<br>
> for lipid analysis).<br>
><br>
> -Justin<br>
><br>
> Berk Hess wrote:<br>
> > Hi,<br>
> ><br>
> > This is a feature.<br>
> > For a long time users have been complaining that pdb2gmx renumbers the<br>
> > residues in a protein.<br>
> > I have now changed this such that the residue numbers in the pdb are<br>
> > retained.<br>
> > But for for instance solvent you would not like to have this behavior.<br>
> > So I decided to keep numbering single-residue molecules.<br>
> > If you also want you lipids to continue numbering, you'll have to set<br>
> > the env.var. GMX_MAXRESRENUM<br>
> > to the number of residues in a lipid.<br>
> ><br>
> > Berk<br>
> ><br>
> > ------------------------------------------------------------------------<br>
> > Date: Mon, 18 Jan 2010 13:05:42 -0200<br>
> > From: <a href="mailto:stefhoor@gmail.com">stefhoor@gmail.com</a><br>
> > To: <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
> > Subject: [gmx-users] trjconv problem<br>
> ><br>
> > I am trying to use trjconv to extract frames from my protein/membrane<br>
> > system in order to analyse with gdmat. Before using git's latest gromacs<br>
> > version, everything worked just fine. But now, every coordinate file<br>
> > trjconv gives me has the DPPC lipid molecules without numbering. It is<br>
> > like if my membrane was formed of a single DPPC (huge) residue. So<br>
> > instead of generating DPPC residue 1, 2, 3, 4 ....etc, trjconv gives me<br>
> > DPPC residue 1 with 5000 atoms.<br>
> > Some light on the matter would be great.<br>
> > Thanks<br>
> ><br></blockquote></div><br>I thank you all for the contributions, but the problem is that the command "trjconv" is making funny things and not pdb2gmx. My residues are all numbered correctly. Even my coordinate.gro file that is generated at the end of the simulation has the correct numbering. The problem is specifically with "trjconv". <br>
The output from trjconv comes out like this:<br>"...<br>66ASN O 671 3.212 3.554 5.248<br>66ASN HO 672 3.258 3.625 5.302<br> 1DPP C33 673 0.623 5.221 1.344<br> 1DPP C34 674 0.461 5.360 1.434<br>
1DPP C35 675 0.697 5.438 1.420<br> 1DPP N 676 0.604 5.327 1.444<br> 1DPP C32 677 0.607 5.278 1.583<br> 1DPP C31 678 0.722 5.193 1.637<br> 1DPP O32 679 0.860 5.229 1.645<br>
1DPP P 680 0.954 5.102 1.677<br> 1DPP O33 681 0.882 4.986 1.620<br> 1DPP O34 682 1.094 5.137 1.647<br> 1DPP O31 683 0.929 5.088 1.836<br> 1DPP C3 684 1.003 5.184 1.912<br>
1DPP C2 685 0.946 5.209 2.052<br> 1DPP O21 686 0.988 5.332 2.111<br> 1DPP C21 687 0.939 5.454 2.087<br> 1DPP O22 688 0.893 5.479 1.976<br> 1DPP C22 689 0.952 5.554 2.195<br>
1DPP C23 690 1.086 5.626 2.185<br> 1DPP C24 691 1.132 -0.046 2.310<br> 1DPP C25 692 1.029 0.049 2.372<br> 1DPP C26 693 1.080 0.108 2.504<br> 1DPP C27 694 1.100 0.015 2.624<br>
1DPP C28 695 1.147 0.089 2.750<br> 1DPP C29 696 1.047 0.189 2.809<br> 1DPP C210 697 1.127 0.248 2.925<br> 1DPP C211 698 1.046 0.363 2.985<br> 1DPP C212 699 1.127 0.411 3.105<br>
1DPP C213 700 1.093 0.323 3.226<br> 1DPP C214 701 1.141 0.402 3.348<br> 1DPP C215 702 1.107 0.339 3.483<br> 1DPP C216 703 1.135 0.428 3.604<br> 1DPP C1 704 0.997 5.089 2.132<br>
1DPP O11 705 0.921 5.069 2.251<br> 1DPP C11 706 0.969 4.974 2.332<br> 1DPP O12 707 1.070 4.914 2.295<br> 1DPP C12 708 0.886 4.938 2.449<br> 1DPP C13 709 0.889 4.790 2.488<br>
1DPP C14 710 0.996 4.756 2.592<br> 1DPP C15 711 0.985 4.603 2.597<br> 1DPP C16 712 1.098 4.556 2.689<br> 1DPP C17 713 1.061 4.414 2.730<br> 1DPP C18 714 1.174 4.351 2.813<br>
1DPP C19 715 1.189 4.364 2.965<br> 1DPP C110 716 1.330 4.335 3.017<br> 1DPP C111 717 1.367 4.338 3.165<br> 1DPP C112 718 1.514 4.304 3.191<br> 1DPP C113 719 1.542 4.329 3.339<br>
1DPP C114 720 1.691 4.329 3.372<br> 1DPP C115 721 1.711 4.346 3.523<br> 1DPP C116 722 1.861 4.344 3.549<br> 1DPP C33 723 1.394 3.827 1.098<br> 1DPP C34 724 1.402 3.868 1.324<br>
1DPP C35 725 1.337 3.651 1.249<br> 1DPP N 726 1.324 3.795 1.223<br> 1DPP C32 727 1.185 3.841 1.213<br> 1DPP C31 728 1.074 3.797 1.309<br> 1DPP O32 729 1.100 3.828 1.446<br>
1DPP P 730 1.029 3.739 1.561<br> 1DPP O33 731 0.892 3.693 1.529<br> 1DPP O34 732 1.132 3.642 1.606<br> 1DPP O31 733 1.020 3.846 1.680<br> 1DPP C3 734 1.118 3.824 1.782<br>
1DPP C2 735 1.123 3.924 1.898<br> 1DPP O21 736 1.254 3.939 1.955<br> 1DPP C21 737 1.326 4.049 1.942...."<br><br>You see that starting from the second nitrogen atom of my DPPC molecules, the residue number should be 68 and not 2. So this is what trjconv is actually giving me. It is renumbering the solvent molecules but it is not recognising different DPPC molecules. Instead it is writing my DPPC membrane as one single huge molecule.<br>