Previous Thread
Next Thread
Print Thread
#216 10/23/03 09:26 AM
Joined: Sep 2003
Posts: 4,850
Likes: 10
lennart Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 4,850
Likes: 10
*FILENAME: asa-protein.str
*PURPOSE: compute Lee-Richards type accessible surface area of a protein probe radius 1.4A
*AUTHOR: Lennart Nilsson, Karolinska Institutet, October 2003
*
!RETURNS:
! ASA computed for the selection defined as SOLUTE (which also gives the context to be used)
! is returned in variable ASA
! One could extend this to subdivide the calculation based on atom or aa type
!ASSUMPTIONS:
! psf and coordinates in MAIN coordinate set
! selection SOLUTE defined
!USAGE EXAMPLE
! read rtf,para,psf, coor
! define SOLUTE sele segid prot end
! stream asa-protein.str

coor surf rprobe 1.4 select SOLUTE end
! ASA per atom is now in WMAIN array, sum up over selected atoms
scalar wmain statistics select SOLUTE end
set ASA ?STOT
!scalar wmain statistics select "NONPOLAR" end
!set ASANP ?STOT

return


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
lennart #217 03/18/09 01:51 PM
Joined: Feb 2009
Posts: 40
N
nqn Offline
Forum Member
Offline
Forum Member
N
Joined: Feb 2009
Posts: 40
Hi,
I would like to know which residues are accessible by the solvent. Is it possible with this script (if I understand correctly, this script returns only area of solute accessible, isn'it?)?

Thanks for advice!

nqn #218 03/18/09 01:57 PM
Joined: Sep 2003
Posts: 4,850
Likes: 10
lennart Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 4,850
Likes: 10
Yes. The key is here:
coor surf rprobe 1.4 select SOLUTE end
! ASA per atom is now in WMAIN array,


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
lennart #219 03/18/09 02:12 PM
Joined: Feb 2009
Posts: 40
N
nqn Offline
Forum Member
Offline
Forum Member
N
Joined: Feb 2009
Posts: 40
yes but how may I see to WMAIN?

Sorry...

nqn #220 03/18/09 02:43 PM
Joined: Sep 2003
Posts: 4,850
Likes: 10
lennart Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 4,850
Likes: 10
write/print coor ...
or
scalar ...


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
lennart #221 03/18/09 02:53 PM
Joined: Feb 2009
Posts: 40
N
nqn Offline
Forum Member
Offline
Forum Member
N
Joined: Feb 2009
Posts: 40
thanks a lot Lennart! Area per atom correspond at the last column? Unit is A²?

nqn #222 03/18/09 02:55 PM
Joined: Sep 2003
Posts: 4,850
Likes: 10
lennart Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 4,850
Likes: 10
Correct


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden

Moderated by  lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u1 Page Time: 0.007s Queries: 28 (0.003s) Memory: 0.7559 MB (Peak: 0.8173 MB) Data Comp: Off Server Time: 2022-10-06 07:54:54 UTC
Valid HTML 5 and Valid CSS