Previous Thread
Next Thread
Print Thread
Solvent accessible area script and radial distribution script
#18819 09/03/08 02:51 PM
Joined: Sep 2008
Posts: 5
L
Forum Member
OP Offline
Forum Member
L
Joined: Sep 2008
Posts: 5
Dear all:

I am a new user. I have a finished simulation and I have all the output trjectory files (*.dcd). I am wondering, for the whole trajectory file, how I can calculate the total SAS for the selected atoms, and g(r) of water proton whithin certain distance away from selected atoms - selected atoms.

Re: Solvent accessible area script and radial distribution script
Lingling_Shen #18820 09/03/08 02:58 PM
Joined: Sep 2003
Posts: 4,793
Likes: 2
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 4,793
Likes: 2
See Sample scripts for basic tasks in CHARMM in the Script Archive for an example of SAS (ASA) calculation. RDFs are described with examples in corman.doc (coor anal ....).


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
Re: Solvent accessible area script and radial distribution script
lennart #18821 09/03/08 07:38 PM
Joined: Sep 2008
Posts: 5
L
Forum Member
OP Offline
Forum Member
L
Joined: Sep 2008
Posts: 5
Thanks.
I am using the MD output from Discovery Studio. I am running nucleotide simulation. I am trying to use your rgyr script to calculate radial distribution. This is my script:
-----------
stream "ModelBuild.stream"

STREAM "noneSolvation.stream"

STREAM "nonbond.stream"

read coor card name T10_3N_single_prod.crd

open unit 51 read unform name T10_3N_single_prod.dcd

open unit 21 read unform name @9pept500.cor
open unit 31 write form name @9pept500.gonh
coor anal select type HT end - ! Water oxygens
site select type NT end multi - ! and the amide hydrogens
firstu 1 nunit 1 skip 10 - ! trajectory specification
isdist 31 - ! do the g(r) (here solute-solvent)
mgn 100 dr 0.1 - ! comp. g(r) at MGN points separated by DR
rdsp 2.5 -
rsph 999.9 - ! we use ALL waters for the calculation
xbox @6 ybox @7 zbox @8 ! and we did use PBC


-------------
It does not work. I am wondering what is wrong?
Thanks,

-L

Attached Files
18989-test.txt (0 Bytes, 728 downloads)
Re: Solvent accessible area script and radial distribution script
Lingling_Shen #18822 09/03/08 08:25 PM
Joined: Sep 2003
Posts: 4,793
Likes: 2
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 4,793
Likes: 2

1. Read the error message - you have some kind of problem with your basic setup, in particlar with the topology/parameter files. It usually (read: always) is important to pay close attention to the first error message. When you want (free) help from others, do not hide the error message in an attachment.
You have not yet reached the part of the problem you asked about in the first place.
2. You have a commercial version of CHARMM, and thus you have the privilege of commercial support via Accelrys - use it, your company has probably paied real money for it!
3. The stream commands to unknown files makes it difficult to judge the quality of your input.
4. Still there are several inconsistencies in your input. Do you actually intend to read a trajectory contained in the file named pept500.cor? And what does the variable @9 represent in this command?
"open unit 21 read unform name @9pept500.cor"

!Innocent, but what makes you think that "type HT" refers to water oxygens?
coor anal select type HT end - ! Water oxygens
!Recheck "type NT"
site select type NT end multi - ! and the amide hydrogens
!What file did you open for reading on unit 1, referenced below?
firstu 1 nunit 1 skip ...
!RDSP 2.5???
!Did you use Periodic Boundary Conditions (of the simple type supported by this CHARMM command)? What about variables 6,7, and 8?

I hope the above helps you formulating good questions to Accelrys support!


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden

Moderated by  lennart, rmv 

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

PHP: 5.6.33-0+deb8u1 Page Time: 0.008s Queries: 23 (0.005s) Memory: 0.9144 MB (Peak: 1.0015 MB) Data Comp: Off Server Time: 2020-09-25 01:08:55 UTC
Valid HTML 5 and Valid CSS