From: Ritu Arora (ritu.arora_at_mail.concordia.ca)
Date: Sun Nov 08 2020 - 18:24:15 CST
Dear All,
I wish to run a simulation of 150 molecules of trehalose. I adapted a psfgen script from earlier NAMD threads, but I am not able to get the script working. The trial script for the carbohydrate section of my system follows:
package require psfgen
topology top_all36_carb.rtf
mol new carb_only.pdb
set sel [atomselect top 'all']
for {set i 0} {$i<150} {incr i} {
resetpsf
segment L$i {
residue 1 AGLC
residue 2 AGLC
}
pdbalias atom 'AGLC 1' O1 1O1
pdbalias atom 'AGLC 1' C1 1C1
pdbalias atom 'AGLC 2' C1 2C1
#
patch 11aa L$i:1 L$i:2
$sel set segname L$i
$set moveby [list 1 1 1]
$set writepdb tmp.pdb
regenerate angles dihederals
coordpdb L$i tmp.pdb
}
guesscoord
writepdb $name.pdb
writepsf $name.psf
The confusion is that there are two AGLCs in trehalose and each has either resid 1 or 2, rather than sequential increments for whole molecule of trehalose. The pdbalias command also seems to be wrong.
I also tried to split the two fragments of trehalose in two separate pdbs:
pdbalias atom 'AGLC 1' O1 1O1
pdbalias atom 'AGLC 2' C1 2C1
for {set i 0} {$i < 150} {incr i} {
segment TRE1$i {pdb output/carb1.pdb }
segment TRE2$i {pdb output/carb2.pdb }
patch 11aa TRE1$i:1 TRE2$i:2
coordpdb output/carb1.pdb TRE1$i
coordpdb output/carb2.pdb TRE2$i
}
I would be thankful for any suggestions to improvise the script.
Best,
Ritu
This archive was generated by hypermail 2.1.6 : Fri Dec 31 2021 - 23:17:10 CST