From: amin_at_imtech.res.in
Date: Sat May 24 2014 - 14:32:55 CDT
Dear all,
I am trying to run a simulation of a protein in namd using amber ff. For
constraining the protein backbone, I generated a pdb file as follows.
Loaded prmtop file and then inpcrd file.
set all [atomselect top all]
set sel [atomselect top "protein and noh"]
$all set beta 0
$sel set beta 1
$all writepdb protein_cons.pdb
Now when I see this pdb file in vmd by coloring based on beta, it looks
good. But when I use this pdb for fixing or constraining atoms, it
doesn't work and all the atoms keep moving throughout the simulation.
The relevant section of my config file is
fixedAtoms on
fixedAtomsForces on
fixedAtomsFile protein_restrained.pdb
fixedAtomsCol B
constraints on
consRef protein_ca.pdb
consKFile protein_ca.pdb
consKCol B
This method works when i use a psf/pdb file which saves the pdb in a
slightly different format than on using top/crd files. For example
for psf/pdb
CRYST1 0.000 0.000 0.000 90.00 90.00 90.00 P 1 1
ATOM 1 N POPCL 1 -68.750 -75.621 20.887 1.00 1.00 L11 N
ATOM 2 C11 POPCL 1 -68.023 -74.839 22.005 1.00 1.00 L11 C
ATOM 3 C12 POPCL 1 -68.874 -74.862 19.662 1.00 1.00 L11 C
ATOM 4 C13 POPCL 1 -70.159 -76.034 21.305 1.00 1.00 L11 C
ATOM 5 C14 POPCL 1 -68.115 -76.902 20.570 1.00 1.00 L11 C
for top/crd
CRYST1 0.000 0.000 0.000 90.00 90.00 90.00 P 1 1
ATOM 1 N TYR X 1 19.125 9.627 -48.120 0.00 1.00
ATOM 2 H1 TYR X 1 18.710 10.461 -48.511 0.00 0.00
ATOM 3 H2 TYR X 1 19.507 9.072 -48.872 0.00 0.00
ATOM 4 H3 TYR X 1 19.709 9.891 -47.339 0.00 0.00
ATOM 5 CA TYR X 1 18.046 8.808 -47.540 0.00 1.00
Can this be the reason for not getting required constraints? I would be
really thankful for any help.
Regards.
Amin.
This archive was generated by hypermail 2.1.6 : Thu Dec 31 2015 - 23:20:49 CST