Our approach to labelling residues and linkages doesn't work for ano-ano linkages
- Dominant language
- C++
- Stars
- 7
- Forks
- 11
- PR merge metrics
- No merged PRs in 30d
Description
For non ano-ano linkages our current labeling system works fine:
IndexOrdered:
DNeu5Aca2-6DGalpb1-4DGlcpNAc[3S]b1-2DManpa1-3[DGlcpNAcb1-4][DManp[2S,3Me]a1-6DManpa1-6]DManpb1-4DGlcpNAc[6Me]b1-4DGlcpNAcb1-OH
IndexOrderedLabeled:
DNeu5Ac&Label=residue-14;a2-6&Label=link-12;DGalp&Label=residue-13;b1-4&Label=link-11;DGlcpNAc&Label=residue-12;[3&Label=link-13;S&Label=residue-15;]b1-2&Label=link-10;DManp&Label=residue-11;a1-3&Label=link-9;[DGlcpNAc&Label=residue-10;b1-4&Label=link-8;][DManp&Label=residue-7;[2&Label=link-7;S&Label=residue-9;,3&Label=link-6;Me&Label=residue-8;]a1-6&Label=link-5;DManp&Label=residue-6;a1-6&Label=link-4;]DManp&Label=residue-5;b1-4&Label=link-3;DGlcpNAc&Label=residue-3;[6&Label=link-2;Me&Label=residue-4;]b1-4&Label=link-1;DGlcpNAc&Label=residue-2;b1-&Label=link-0;OH&Label=residue-1;
However for ano-ano linkages, this means the "a" gets lost.
IndexOrdered:
DGlcpa1-1DGlcpa
IndexOrderedLabeled:
DGlcp&Label=residue-2;a1-1&Label=link-0;DGlcp&Label=residue-1;
Essentially it's because we put the alpha of DGlcpa into the linkage, and not the sugar residue. This follows the convention of what scientists do, i.e. we think like this:

But what's cleaner to program is that the linkage is e.g. "1-2" and the sugar should be "DGlcpa"
With the current approach I'm losing the a info for Glc

I'm not even sure where we would want the a to be in this representation?
Contributor guide
No contributing guide indexed for this repository
Assessment
This issue has not been assessed yet.