Re: Specific colvar component needed.

From: Aron Broom (broomsday_at_gmail.com)
Date: Mon May 07 2012 - 11:02:30 CDT

can I ask why you need either of these restrictions? The first restriction
at the very least seems pointless. For having the CA atom fixed in place,
I suspect your are trying to mimic AFM experiments or something, but
physics is all relative, if you simply have the colvar be the distance
between the C-terminal CA and the CA of the N-terminus, this is no
different, the relative distance between the two is all that matters
(saying one is fixed is meaningless I would think, maybe I'm missing the
boat though).

In terms of the second restriction, I can see that perhaps this is needed
to accurately recapitulate a pulling experiment, but given the rotational
degrees of freedom around all the residues, I don't think just pulling them
apart along a general distance colvar will be any different than
restricting that pulling along one axis in particular, maybe someone more
mathematically inclined can comment.

~Aron

On Mon, May 7, 2012 at 11:55 AM, Germán Andrés Miño <germino_at_u.uchile.cl>wrote:

> Hi namd2.8 users.
>
> I would like to do and protein unfolding ABF calculation in an water
> box with the following restrictions:
>
> CA atom of the C-terminal fixed in space.
> CA atom of the N-terminal moving along the (0,0,Z) axis.
>
> I have tried with distanceZ (see below) component but my system drift
> away from the Z axis.
>
> I need to minimize the size of the water box (actually a rectangle)
> because limitations in numbers of processors and computational time
> available.
>
> The colvar I am using now is:
>
>
> colvarsTrajFrequency 2000
> colvarsRestartFrequency 20000
>
> colvar {
> name d
>
> width 0.1
>
> lowerboundary 50.36
> upperboundary 380.00
>
> lowerwallconstant 10.0
> upperwallconstant 10.0
>
> distanceZ {
> main { atomNumbers 1921
> }
> ref {
> atomNumbers 959
> }
> axis (0.0, 0.0, 1.0)
>
> }
> }
>
> abf {
> colvars d
> fullSamples 500
> hideJacobian
> }
>
>
> Any suggestion is truly valuated.
> Thanks in advance.
>
> German Miño Galaz.
> germino_at_u.uchile.cl
> Postdoctoral Fellow of Grupo de Namomateriales (www.gnm.cl)
> Departamento de Fisica
> Universidad de Chile.
> Palmeras 3425, Nunoa.
> Santiago Chile
> Tel: 0056-2-978-7439
> Fax: 0056-2-2712973
>
>

-- 
Aron Broom M.Sc
PhD Student
Department of Chemistry
University of Waterloo

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