Missing DNA bases with negative residue number

[mihir41@terpmail.umd.edu]

The PDB ID: 9GBV has 217 bp long dsDNA (Chain E & F). I plan to perform a mutation in some bases. But many of the bases are not displayed. Specifically, all the bases having a negative residue number are omitted by the portal.

[xiangjun]

Hi,

Thank you for using our tools and for reaching out on the forum.

To help us investigate this further, could you please provide a step-by-step description of your workflow? Specifically, please include the exact commands or options you used so that we can reproduce the issue on our end. This will help us determine why the residues with negative numbers are not being displayed as expected.

Best regards,

Xiang-Jun

[mihir41@terpmail.umd.edu]

Hello Xiang-Jun,

I used the online tool for DNA mutation.

On the page https://web.x3dna-dssr.org/mutation

I entered PDB ID: "9GBV" in the left box. Then clicked on the "List all bases" button.

On the next page, I see bases listed for chains E (from resid number 0 to 121) and F (from resid number 0 to 80).

However, if you check the sequence map on the RCSB (https://www.rcsb.org/3d-view/9GBV), you would see that chain E residues range from -37 to 121, and chain F residues range from -78 to 80.

Please let me know if you need more information.

Best,
Mihir

[xiangjun]

Hi Mihir,

Great catch! I have successfully reproduced the issue. The bug is due to a parsing error within the wDSSR web interface. As shown below, the core DSSR engine itself handles negative residue numbering without any issues.

Execute the following command to see the enumeration:

Code: [Select]

x3dna-dssr mutate -i=9GBV.pdb --enum
An excerpt of the output is as follows:

# For each entry to be mutated, remove the leading #, and add to=base

# To mutate A.A9 to G in PDB entry 1ehz (tRNA):
#          #  chain=A name=A num=9  #9  A.A9
#       --->  chain=A name=A num=9 to=G

# To mutate B.DT19 to DC in PDB entry 355d (B-DNA):
#          #  chain=B name=DT num=19  #19  B.DT19, pairedNt=A.DA6
#       --->  chain=B name=DT num=19 to=DC

# Empty or comment (starting with #) lines are ignored.

#  chain=E name=DT num=-37  #1  E.DT-37, pairedNt=F.DA80
#  chain=E name=DG num=-36  #2  E.DG-36, pairedNt=F.DC79
#  chain=E name=DC num=-35  #3  E.DC-35, pairedNt=F.DG78
#  chain=E name=DC num=-34  #4  E.DC-34, pairedNt=F.DG77
#  chain=E name=DA num=-33  #5  E.DA-33, pairedNt=F.DT76
#  chain=E name=DT num=-32  #6  E.DT-32, pairedNt=F.DA75
#  chain=E name=DC num=-31  #7  E.DC-31, pairedNt=F.DG74
#  chain=E name=DA num=-30  #8  E.DA-30, pairedNt=F.DT73
#  chain=E name=DG num=-29  #9  E.DG-29, pairedNt=F.DC72
#  chain=E name=DA num=-28  #10  E.DA-28, pairedNt=F.DT71
#  chain=E name=DA num=-27  #11  E.DA-27, pairedNt=F.DT70
#  chain=E name=DA num=-26  #12  E.DA-26, pairedNt=F.DT69
#  chain=E name=DA num=-25  #13  E.DA-25, pairedNt=F.DT68
#  chain=E name=DC num=-24  #14  E.DC-24, pairedNt=F.DG67
......

We will resolve this wDSSR error shortly. In the meantime, please continue using the tool and report any further issues you encounter.

If it is an option for your workflow, I recommend installing the standalone version of DSSR on your computer and using the command-line interface directly.

Best regards,

Xiang-Jun