VMD-L Mailing List
From: John Stone (johns_at_ks.uiuc.edu)
Date: Thu Dec 07 2006 - 11:04:39 CST
- Next message: Nels Thorsteinson: "problem/bug? calling hbonds within a loop over the residues (tcl)"
- Previous message: Streutker, J: "drawing helices from partial aminoacids"
- In reply to: Streutker, J: "drawing helices from partial aminoacids"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
Hi,
In order to draw ribbons, your structure needs to contain
both CA and O atoms, for cartoons you'll need CA and P atoms,
for tubes you just need CA atoms. STRIDE probably won't like the
coarse-grained structure very much, but it can often be convinced
to work by atom renaming etc. I'm not sure what the tube error
you're getting is, but if you send me a copy of the structure(s) you're
trying to load, I'd be happy to have a look and give suggestions.
John Stone
vmd_at_ks.uiuc.edu
On Thu, Dec 07, 2006 at 11:40:24AM +0100, Streutker, J wrote:
> Hi everyone,
>
> I`m having trouble with VMD. I`m using trajectory and PDB file, so simulate a lipid bilayer. All of this is done coarse-grained, so I only have partial amino acids (# atoms ranges from 1 to 3), instead of complete amino acids. A part of my trajectory contains 5 helices, but I cant get VMD/stride te draw them as (new)ribbons (or some other kind of helix shape). If I`m not mistaken, ribbons are drawn through the CA atoms of each aminoacid. And my file contains all CA atoms, from every molecule from the helix. But nothing seems to appear. When I select Tubes, I get an error: Internal error in draw_tube for CA atoms: last_resid = 2 but no first atom was added. So I thought: I cheat and rename every molecule to CA. Unfortunately that didnt work either. Well it works partially. I get some lines and bends, but no real helix shape. After this I didnt really know what to do next, so I used PSFgen to complete the aminoacids. This worked, I see helices. There is only one major problem: some H atoms and every O and C
atom at the end of a aminoacid, get the coordinates: 0.0, 0.0, 0.0. Adding to that, they also have an atomic charge of -1.0 in the psf file. So it is happing every time, with the same kind of atom, but why? perhaps it has something to do with the fact that the are added by PSFgen, although most of the are placed correctly by PSFgen...
>
>
> Kind regards,
>
> J. Streutker
> University of Groningen
> Netherlands
-- NIH Resource for Macromolecular Modeling and Bioinformatics Beckman Institute for Advanced Science and Technology University of Illinois, 405 N. Mathews Ave, Urbana, IL 61801 Email: johns_at_ks.uiuc.edu Phone: 217-244-3349 WWW: http://www.ks.uiuc.edu/~johns/ Fax: 217-244-6078
- Next message: Nels Thorsteinson: "problem/bug? calling hbonds within a loop over the residues (tcl)"
- Previous message: Streutker, J: "drawing helices from partial aminoacids"
- In reply to: Streutker, J: "drawing helices from partial aminoacids"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]