VMD-L Mailing List
From: Martin, Erik W (Erik.Martin_at_stjude.org)
Date: Tue Feb 18 2014 - 14:50:23 CST
- Next message: Carlos Simmerling: "ACS COMP Division student and junior faculty awards for Fall 2014 meeting in San Francisco"
- Previous message: KK R: "Extracting complex of protein-ligand+water from every frame of DCD trajectory"
- In reply to: KK R: "Extracting complex of protein-ligand+water from every frame of DCD trajectory"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
Add "same residue as" before your water selection
-Erik
Sent from my iPhone
On Feb 18, 2014, at 12:44 PM, "KK R" <kkres14_at_gmail.com<mailto:kkres14_at_gmail.com>> wrote:
Hi,
I have a result from NAMD simulation. I noticed that there were two water clusters around the ligand within the active site of the protein. Therefore, I tried to extract a set of protein+ligand+water (within 10A of ligand) from every frame of my DCD trajectory. I could accomplish it using the following script in the Tk console.
--------------------------------------------------------------------------------------------------------------------------------
set sel [atomselect top "{protein} or {resname LIG} or {water within 10 of (resname LIG)}"]
set n [molinfo top get numframes]
set outname "protligwater_"
for { set i 0 } { $i < $n } { incr i } {
$sel frame $i
$sel update
set name [format "%s%.9d.pdb" $outname $i]
$sel writepdb $name
}
---------------------------------------------------------------------------------------------------------------------------------
where LIG is the segid for ligand.
But, the above script writes protein+ligand+(water residues as well as only part of water residues (like H or OH etc) lying within that radius around ligand).
Can anybody suggest me about how to specify in the command line to select complete water residues lying within any particular radius, even if only part of some water resides with that radius?
I hope that my question is clear to you all. If not, please let me know.
Thanks in advance,
Best regards,
Roy
________________________________
Email Disclaimer: www.stjude.org/emaildisclaimer
Consultation Disclaimer: www.stjude.org/consultationdisclaimer
- Next message: Carlos Simmerling: "ACS COMP Division student and junior faculty awards for Fall 2014 meeting in San Francisco"
- Previous message: KK R: "Extracting complex of protein-ligand+water from every frame of DCD trajectory"
- In reply to: KK R: "Extracting complex of protein-ligand+water from every frame of DCD trajectory"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]