atom selection for nucleic acid residues

Issue #771 closed
Former user created an issue

Hii all, i am new to bio3d getting thus error. my nucleic acid residues are 48 while with atoms it shows 1522. from galaxy i am able to generate the graph for all 1522 atoms but i want graph of only 48 nucleotide. please help me with this. i just want the C alpha atoms and nucleic acid residues DCCM graph. in galaxy if i want to do which kind of residues i have to select? as i select residues like residue id A, G , C, U it give error. i had attached by pdb file for the same.

read.pdb(wtpdb , multi=TRUE)

Call: read.pdb(file = wtpdb, multi = TRUE)

Total Models#: 1 Total Atoms#: 2984, XYZs#: 8952 Chains#: 1 (values: A)

 Protein Atoms#: 1462  (residues/Calpha atoms#: 94)
 Nucleic acid Atoms#: 1522  (residues/phosphate atoms#: 48)

 Non-protein/nucleic Atoms#: 0  (residues: 0)
 Non-protein/nucleic resid values: [ none ]

Protein sequence: GSSGSSGMSGGLAPSKSTVYVSNLPFSLTNNDLYRIFSKYGKVVKVTIMKDKDTRKSKGV AFILFLDKDSAQNCTRAINNKQLFGRVIKASIAI

Nucleic acid sequence: AUUAAAGACCAACCAACUUCUGUUCUCUAAACGAACUUUAAAUUAAUA

  • attr: atom, xyz, calpha, call

    p.inds <- atom.select(wtpdb, "phosphates") Error in UseMethod("atom.select") : no applicable method for 'atom.select' applied to an object of class "character" sele = atom.select(wtpdb, "protein") Error in UseMethod("atom.select") : no applicable method for 'atom.select' applied to an object of class "character" sele = atom.select(wtpdb, "protein") Error in UseMethod("atom.select") : no applicable method for 'atom.select' applied to an object of class "character" pdb <- read.pdb('wtpdb') Error in read.pdb("wtpdb") : No input PDB file found: check filename

Comments (3)

  1. Xinqiu Yao

    Hi,

    It is not very clear what you intended to do exactly. Note that this is not the place to ask questions about galaxy, which has no relationship with bio3d.

    Regarding the ‘atom.select()’ errors, instead of using the unsupported string 'phosphates', use a combination of the type of polymer and element. For example, p.inds ← atom.select(wtpdb,"nucleic", elety="P"). See ?atom.select for more detail.

    Hope it helps.

  2. Log in to comment