From: Ajasja Ljubetič (
Date: Wed Feb 22 2012 - 11:51:51 CST

> >c) decide that you don't want the calcium to leave the binding site and
> lock the distances and orientation using additional >bonds.

> do u mean to say that , show i fix these calcium atoms using "fixedatoms"
> options [setting Bfactor column to 1] or apply tclforces to control the
> distances and orientations ..?
> No, I mean that you change the force field parameters to
define additional covalent bonds between the Ca and the atoms to which the
Ca is coordinated.
This is probably quite a bit of work, since you have to estimate the
bond constants and angle bending constants. (Or just set them to high
numbers to have the Ca bound rigidly)

I also have another idea: perhaps you could use the colvar
between the Ca and the coordinating atoms) and lock it to the desired value
using harmonic restraints<>.
This is probably simpler, it it works.

Best regards,

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:14 CST