Hi List,
I have a system which has two antiparallel monomers. I m using the following
commands below using Amber 9. I can get the correct distances between Ca-Ca
but the values get messed up with helix.rmsd and equil.rmsd after 90 ns. Any
suggestions?
....
....
trajin prod197.mdcrd.gz 1 50 1
trajin prod198.mdcrd.gz 1 50 1
trajin prod199.mdcrd.gz 1 50 1
trajin prod200.mdcrd.gz 1 50 1
trajin prod201.mdcrd.gz 1 50 1
center :1-144 origin mass
image origin center
distance CALPHA :62.CA :137.CA out dist_CALPHA.list
distance HIPGLN :7.H :32.NE2 out dist_HIPGLN.list
reference pro_solvated.inpcrd
rms reference out equil.rmsd .CA,C,N
rms reference out helix.rmsd :54-71.N,CA,C,:126-143.N,CA,C
rms reference out helix1.rmsd :54-71.N,CA,C
rms reference out helix2.rmsd :126-143.N,CA,C
Many thanks,
Regards,
Neha Gandhi Gajaria,
School of Biomedical Sciences,
Curtin University of Technology,
GPO Box U1987 Perth,
Western Australia 6845
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Dec 08 2009 - 22:30:02 PST
