This problem cannot have an ultimate solution unless your chains are of the
same length and you postulate that residues should be corresponded one-to one.
But if this is the case, LSQKAB is the tool.
If correspondence between C-alphas is ambiguous, you need to align structures
in order to find the best set of corresponding atoms. This is what SSM (or
GESAMT) does and it cannot be guaranteed that a an atom in one chain will have
a match in another chain. Note that alignment != best rms superposition.
Eugene
On 12 Sep 2014, at 14:45, KAUSHIK H.S. wrote:
Hi,
I want to calculate main-chain and side-chain deviations between sets of
identical proteins modeled into different 3D structures. While SSM superposes
well, it calculates RMS only for the residues considered for superposing.
Though LSQKAB program calculates RMS for all positions, fails when the topology
is not closely related.
For now, I do SSM superpose, identify the residues, explicitly mention them in
LSQKAB program and then calculate RMS for each position.
However, is there a way we can make SSM calculate RMS for all positions?
Or,
Is there a program which will give all atom RMS for a previously superposed
protein structure ensemble?
Thanks,
Kaushik,
Molecular Biophysics Unit,
IISc, India
--
Scanned by iCritical.