Preventing rotation of a protein

From: Thomas Albers (talbers_at_binghamton.edu)
Date: Wed Jul 17 2013 - 21:18:16 CDT

Hello!

I have got a system with a 400 amino acid dimer that fits into a
128*96*96 A oblong water box. That's about 100K atoms. A cubic
128*128*128 water box would be about 200K atoms.

If I restrain the rotation of the protein with an orientation colvar,
as hinted at on p. 123 of the NAMD 2.9 manual there's a *substantial*
speed penalty; IIRC it's so bad that I might as well drop the colvar
and use a cubic box instead. However, in the recent protein-ligand
binding tutorial (this one:
http://www.ks.uiuc.edu/Publications/Stories/plbtutorial/) the colvar
is defined only for the backbone of the protein, this is much faster,
but I fear I might introduce bias in the internal degrees of freedom.

Does anyone on the list have some input?

Regards,
Thomas

This archive was generated by hypermail 2.1.6 : Wed Dec 31 2014 - 23:21:26 CST