Hi Xiangjun,
Thank you so much for your reply.
(1) Format:
I am thinking the similar format for DSSR should work great. If I remember it correctly, in DSSR json, you have an object called "hbonds" which you store all the hbonds pair right? For example:
[yuze@summer:dssr_learn] jq .hbonds[1] 6tna.json
{
"index": 2,
"atom1_serNum": 20,
"atom2_serNum": 1564,
"donAcc_type": "standard",
"distance": 2.663,
"atom1_id": "N1@1:A.G1",
"atom2_id": "N3@1:A.C72",
"atom_pair": "N:N",
"residue_pair": "nt:nt"
}
I think similar structure is fine. You can also add stacking interaction as a separate object in the json file.
By the way, in the introduction, you said
SNAP recognizes and outputs base-amino-acid H-bonding and stacking interactions.
Does that mean SNAP cannot output nucleic acid backbone interaction (only nucleic acid base)?
(2) Feedback:
I am currently training an undergraduate to analyze the interaction mode of RNA and protein. We will do the analyze for all the deposited protein RNA structures in RCSB. So we will definitely provide useful feedback when we encounter some wired cases.
Best,
Honglue