Hi,
I'm using CPPTRAJ to do a PCA analysis on a dimer. However, monomers are non-covalently attached therefore moving and rotating a lot from one another. I'm having trouble aligning both monomers at the same time because of this. If I align the entire protein due to movements of monomers some parts of the monomers are not aligning well and that gives me an incorrect output. is there a way to align these structures? I have also copied my commands below (which also includes correlation matrix and distance matrix commands)
parm nowater_glyco.psf
trajin nowater1.dcd
autoimage anchor :588-1174
createcrd deglyco-average-traj
run
parm ionizedph77.pdb [structX]
reference ionizedph77.pdb parm [structX]
crdaction deglyco-average-traj rmsd reference [structX] :1-1174.CA,C,O,N
run
crdaction deglyco-average-traj matrix covar \
name covar :1-1174.CA
runanalysis diagmatrix covar out evecs.dat \
vecs 10 name myEvecs \
nmwiz nmwizvecs 5 nmwizfile glyco2A.nmd nmwizmask :1-1174.CA
run
crdaction deglyco-average-traj matrix correl :1-1174.CA :1-1174.CA out glyco2correlmatA.dat
run
crdaction deglyco-average-traj matrix dist :1-1174.CA :1-1174.CA out glyco2distmatA.dat
Thank you !
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jan 21 2022 - 08:30:02 PST
