From: Vermaas, Joshua (Joshua.Vermaas_at_nrel.gov)
Date: Wed Jun 07 2017 - 12:15:45 CDT
Its n/m. Technically it would be 0/0 when d=d0, but it doesn't evaluate that way within NAMD, which effectively takes the limit at that point, which is n/m. I found it handy to plot the function for myself, especially when I'm toying around with using different values for n, m, and d0.
-Josh
On 06/07/2017 12:38 AM, Souvik Sinha wrote:
Hi,
There I have a query regarding the description of the switching function. Its written that at d = d0 , the function will have a value of n/m. But looking at the function, it seems to be of an indeterminate form at d = d0.
Souvik
On Tue, Jun 6, 2017 at 10:29 PM, Vermaas, Joshua <Joshua.Vermaas_at_nrel.gov<mailto:Joshua.Vermaas_at_nrel.gov>> wrote:
Hi Charles,
If you only care about keeping things out of the pore, but not really the direction, I'd actually use a coordnum collective variable. You can do the same sort of thing with the tcl forces interface, but there the scripting is a bit less straightforward, since instead of forcing the atom purely along z, technically you may want to force things towards the outside. Performancewise tclforces and collective variables are about equal if you choose sane reaction coordinates. If you use coordnum, this is what you'd do. First, you'd define the K ions as well as what you want to keep them away from:
set allsel [atomselect top "all"]
$allsel set beta 0
set ksel [atomselect top "name K"]
set psel [atosmleect top "noh and protein"]
$ksel set beta 1
$psel set beta 2
$allsel writepdb tags.pdb
Then you'd define a collective variable to drive the number of contacts between the ksel and psel groups to 0 (this goes in a seperate file. colvars.conf or something like that).
colvar {
Then you'd add the required elements to your NAMD configuration file to use this new coordination number collective variable to eliminate contacts.
colvars on
-Josh
On 06/05/2017 10:04 PM, Charles Hoying wrote:
I have a membrane protein that I want to equilibrate in NAMD without allowing K ions from the solvation/ionization to get close to it. Specifically, it has a pore with stable K+ binding sites on the outside, so during equilibration before I run simulations I want to keep K+ ions away from that pore.
I was thinking of modifying the "keep_water_out.tcl" script from the NAMD Membrane Proteins Tutorial (http://www.ks.uiuc.edu/Research/namd/mailing_list/ I am open to any and all suggestions whether the script or some other approach can help me solve the problem! Please let me know if there is anything I can clarify to help answer my question.
Best,
This archive was generated by hypermail 2.1.6
: Sun Dec 31 2017 - 23:21:21 CST
name contacts
upperWall 0
upperWallConstant 1
coordNum {
group1 {
atomsFile tags.pdb
atomsCol B
atomsColValue 1
}
group2 {
atomsFile tags.pdb
atomsCol B
atomsColValue 2
}
#These last three parameters are up to you. See http://www.ks.uiuc.edu/Research/namd/2.12/ug/node57.html#SECTION0001341120000000000000
cutoff 6
expNumer 8
expDenom 16
}
}
colvarsConfig colvars.conf
Hi,
Charles
--
Souvik Sinha
Research Fellow
Bioinformatics Centre (SGD LAB)
Bose Institute
Contact: 033 25693275