atom selection for nucleic acid residues
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)
-
-
Not a bug
-
- changed status to closed
Closed - no reply
- Log in to comment
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.