[gmx-users] Qestion about how to define groups by chain id in a protein ?

Liu Shiyong liushiyong at gmail.com
Tue Jan 20 21:17:40 CET 2009


Thanks.

I used the PDB file with H .  It works though it's still a mystery for me
why PDB file without H couldn't work.



On Tue, Jan 20, 2009 at 1:53 PM, Justin A. Lemkul <jalemkul at vt.edu> wrote:

>
>
> Liu Shiyong wrote:
>
> <snip>
>
> The problem lies here:
>
> Including chain 1 in system: 1296 atoms 125 residues
> Including chain 2 in system: 1274 atoms 123 residues
> Including chain 3 in system: 2085 atoms 201 residues
>
> This suggests that chain 2 (Protein B) should contain numbers up to about
> 2500.
>
>         ATOM   1996  O   ASN B 248      49.634   9.874  85.195  1.00  0.00
>>        ATOM   1997  OXT ASN B 248      50.217  10.536  83.158  1.00  0.00
>>        TER
>>        ATOM   1998  N   GLY C 249      70.273  30.186  73.098  1.00  0.00
>>        ATOM   1999  CA  GLY C 249      68.973  30.327  72.421  1.00  0.00
>>
>>
> This is the original .pdb file, then?  The hydrogens will be missing from
> the appropriate groups in the pdb2gmx-processed output structures.  You can
> check the numbering (and pertinent charge groups) in the topology for each
> chain, to be sure.
>
>                make_ndx -f ${df}.pdb  -o ${file}.ndx >
>>        ${file}.output.make_ndx
>>               << _EOF_
>>               del 0-9
>>               chain A and B
>>               chain C
>>               q
>>               _EOF_
>>
>
> This should work, with the right input :)
>
> -Justin
>
>
> --
> ========================================
>
> Justin A. Lemkul
> Graduate Research Assistant
> Department of Biochemistry
> Virginia Tech
> Blacksburg, VA
> jalemkul[at]vt.edu | (540) 231-9080
> http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
>
> ========================================
>



-- 
Shiyong Liu
Postdoc
center for bioinformatics in the university of kansas
Lab: (785)864-1962
Email: syliu at ku.edu (shiyongliu at ku.edu or liushiyong at ku.edu)
Homepage: http://www.people.ku.edu/~syliu
Lab:    http://vakser.bioinformatics.ku.edu/people
Phone:      (785) 864-1962
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090120/32c4182c/attachment.html>


More information about the gromacs.org_gmx-users mailing list