VMD-L Mailing List
From: Irene Newhouse (einew_at_hotmail.com)
Date: Tue Mar 02 2010 - 13:51:25 CST
- Next message: Edelmiro Moman: "Re: measure command problem"
- Previous message: radak004_at_umn.edu: "Running Python within VMD"
- Next in thread: Irene Newhouse: "RE: wrapping dimer units in a Desmond traj"
- Reply: Irene Newhouse: "RE: wrapping dimer units in a Desmond traj"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
I've got Desmond trajectories in which one of the monomer units of a dimeric protein wraps. I can still align on each monomer & get backbone RMSD values for each monomer, but I'd also like to analyze whether the dimer is stable or drifting apart.
So I found the web references to the pbctools plugin. I used the initial frame to determine values for a, b, & c, using the measure minmax command, and set the cell with pbc set "64.4 96.3 99.2" -all . Note that the use of "" is not shown in the documentation, but not using "" resulted in pbctools errors. Then I pbc wrap -sel "protein or resname HEM or resname ACT" -all . The result is disappointing - VMD does a lot of work, but the dimer is not reconstituted. Instead, there are anywhere from 2 to 4 segments of protein. They are distributed differently from the original, but they're still not in closest proximity. Frame 0 contains the unwrapped protein, while most of the other frames are wrapped; I don't know if this might be contributing.
Thanks for your help!
Irene Newhouse
_________________________________________________________________
Hotmail: Powerful Free email with security by Microsoft.
http://clk.atdmt.com/GBL/go/201469230/direct/01/
- Next message: Edelmiro Moman: "Re: measure command problem"
- Previous message: radak004_at_umn.edu: "Running Python within VMD"
- Next in thread: Irene Newhouse: "RE: wrapping dimer units in a Desmond traj"
- Reply: Irene Newhouse: "RE: wrapping dimer units in a Desmond traj"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]