[Biojava-l] Superimposing structure pieces

Daniel Asarnow dasarnow at gmail.com
Sun Oct 31 23:56:05 UTC 2010


I've been trying to pull out pieces of protein chains and superimpose
them...my current code (as generic-ified code snips below) works, but
I wonder if it couldn't be faster.
Has anyone worked on similar methods?  Any other advice?

Best regards everyone,
da

Getting residue CA's as Atom[]:

for (int i; i < length; i++) {
    someAtoms[i] = someChain.getSeqResGroup(start + i).getAtom("CA");
}

Superimposing/aligning:

SVDSuperimposer svds = new SVDSuperimposer(someAtoms1, someAtoms2);
Matrix rot = svds.getRotation();
Atom trans = svds.getTranslation();
for (int i = 0; i < length; i++) {
    Calc.rotate(someAtoms1[i], rot);
    Calc.shift(someAtoms1[i], trans);
}
SVDSuperimposer.getRmsd(someAtoms1, someAtoms2);



More information about the Biojava-l mailing list