VMD-L Mailing List
From: Norman Geist (norman.geist_at_uni-greifswald.de)
Date: Tue Jan 31 2017 - 09:35:48 CST
- Next message: Ruggero Cortini: "Re: tkcon cursors not working"
- Previous message: Norman Geist: "AW: SASA area"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
This looks ok so far, but inside VMD you use the „-restrict“ flag to obtain the surface residue based. How do you do the same when using the external tool? So it’s reasonable you get different results.
Von: Anna Battisti [mailto:battistianna_at_googlemail.com]
Gesendet: Dienstag, 31. Januar 2017 12:31
An: Norman Geist <norman.geist_at_uni-greifswald.de>
Betreff: Re: vmd-l: SASA area
Thanks a lot for replying!
Sorry, what do you mean? maybe the flag I used?
./surf_LINUXAMD64 -W 1 -R 1.4 fort.1000
Anyway, step by steps, I have done the calculation in the following way:
1) starting from pdb file I generated the fort.1000 file:
awk 'BEGIN{while((getline < "output_VMDradius.dat")>0){n++; a[n]=$1}}{printf "%4i%10.4f%10.4f%10.4f%10.4f\n",NR,a[NR],$6,$7,$8}' temp2 > fort.1000
the beginning of the file is:
1 1.4000 31.3700 63.3600 63.0500
2 1.5000 32.1400 62.9100 61.8800
3 1.5000 31.2900 61.9600 61.0200
4 1.5000 32.1000 61.0600 60.1400
5 1.3000 30.4700 61.2400 61.9500
6 1.5000 32.6000 64.1100 61.0600
7 1.3000 31.9900 65.1700 61.1500
8 1.4000 33.6700 63.9300 60.2900
9 1.5000 34.2400 65.0600 59.5500
10 1.5000 35.6200 65.3500 60.1900
11 1.5000 36.6100 65.9800 59.2900
12 1.3000 35.4200 66.2000 61.3300
--------
The first column is the number of the line, the second is the VDW radius (calculated with VMD) then the coordinates (x, y, z)
2) Then I processed the fort.1000, and obtained fort.1000.tri in this way:
./surf_LINUXAMD64 -W 1 -R 1.4 fort.1000
3) I calculated the area of the triangles of the triangles.
Is this procedure correct?
On Tue, Jan 31, 2017 at 12:05 PM, Norman Geist <norman.geist_at_uni-greifswald.de <mailto:norman.geist_at_uni-greifswald.de> > wrote:
How did you tell the surf_LINUXAMD64 the restrict settings that you used in the TCL version?
Von: owner-vmd-l_at_ks.uiuc.edu <mailto:owner-vmd-l_at_ks.uiuc.edu> [mailto:owner-vmd-l_at_ks.uiuc.edu <mailto:owner-vmd-l_at_ks.uiuc.edu> ] Im Auftrag von Anna Battisti
Gesendet: Dienstag, 31. Januar 2017 08:28
An: vmd-l_at_ks.uiuc.edu <mailto:vmd-l_at_ks.uiuc.edu>
Betreff: vmd-l: SASA area
Dear all,
I'm doing the sasa area calculation but something strange is happening: I'm seeing different results by comparing the calculation made with
set output [open "output_SASA.dat" w]
for {set x 0} {$x <= 2084} {incr x} {
set sel [atomselect top "index $x"]
set protein [atomselect top "protein"]
set sasa [measure sasa 1.4 $protein -restrict $sel]
puts $output "$sasa"
}
close $output
and that I obtained calculating the area of each atom as sum of the associated triangles area, when the triangles are obtained with:
./surf_LINUXAMD64 -W 1 -R 1.4 file > file.tri
I have supposed that the procedure to calculate the sasa with VMD is the same , namely surf_LINUXAMD64 is the VMD program.
Kindly, could you give me any kind of suggestion or tell me if there is some specific approximation?
Thanks
Anna
- Next message: Ruggero Cortini: "Re: tkcon cursors not working"
- Previous message: Norman Geist: "AW: SASA area"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]