Applying weak harmonic force between two proteins

Date: Wed May 21 2008 - 11:12:18 CDT

Hello all,
I am interested in setting up a dimer simulation. The simulation cell has two
proteins and I am interested in studying how each protein responds to the
presence of the other protein. To efficiently do this I want to apply a weak
harmonic force between the center of mass of the two proteins to prevent them
from flying apart. I can write a detailed Tcl force script (Identifying all the
protein coordinates, calculating the center of mass, then calculating the force
on each and every atom) but I am afraid that it might be very inefficient. Is there
any way to do this succinctly and efficiently using NAMD's inbuilt commands
and functions.
Thank You

G. Reddy
Institute for Physical Science and Technology
University of Maryland, College Park
MD 20742

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:49:29 CST