From: Vermaas, Joshua (Joshua.Vermaas_at_nrel.gov)
Date: Sat Apr 07 2018 - 10:01:45 CDT
See http://www.ks.uiuc.edu/Research/namd/2.12/ug/node57.html#SECTION0001341240000000000000. This is an example from a system I was working on, where I wanted to change the rotation of a helix by 30 degrees. The helix was tagged with a 1 in the beta column of spin.pdb.
colvar {
name spin
spinAngle {
atoms {
atomsFile spin.pdb
atomsCol B
atomsColValue 1
}
refPositionsFile spin.pdb
refPositionsCol B
refPositionsColValue 1
axis (-0.05879103049526372, 0.37592884870839527, 0.9247816582534958)
}
}
harmonic {
name harm
colvars spin
forceConstant 1000
centers 0
targetCenters -30
targetNumSteps 50000
}
-Josh
On 2018-04-07 00:06:39-06:00 owner-namd-l_at_ks.uiuc.edu wrote:
Hi
I want to calculate helical axis of an alpha-helix of my protein along which I want to rotate that particular helix with respect to rest of the protein. is there any build in function in namd that I can use to calculate it or do I need to write the script to do it? also during the simulation, the atoms in alpha helix will move, changing the helical axis. so do I need to keeps calculating and updating the axis of rotation or will it be safe to calculate it at the beginning of the simulation and don't change it during the simulation? if I am calculating the helical axis on the the fly during simulation, will it slow down the simulation significantly?
Thank you
Dhiraj
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:20:58 CST