How to use pdbsuperpose without .Sequence field?

1 Ansicht (letzte 30 Tage)
Steven
Steven am 3 Dez. 2015
I created a pdb from XYZ coordinates using the function mat2pdb (<http://www.mathworks.com/matlabcentral/fileexchange/42957-read-and-write-pdb-files-using-matlab/content/mat2pdb.m)>. This function does not create a .Sequence field. The matlab function pdbsuperpose.m seems to require this field to do the alignment. I was wondering if there was a solution where I could conduct superposition without the sequence field. I am not working with a protein structure.
Error: Reference to non-existent field 'Sequence'. Error in pdbsuperpose (line 133) chainList1 = {pdbStruct1.Sequence.ChainID}; % all chains

Antworten (2)

Yuebin Zhou
Yuebin Zhou am 7 Dez. 2015
You may want to use getpdb function to first retrieve protein structure data from Protein Data Bank (PDB) database, then you should be able to use pdbsuperpose to process it.
Please check the links below:

Steven
Steven am 7 Dez. 2015
I believe I found a solution. In my case I only had x, y, z coordinates for the structure. It was not a protein structure. However, I wanted to view the x,y,z coordinates with molviewer. Thus there was no corresponding pdb accession associated with the molecule I was interested in. In fact, the structure I was interested in was a chromosome 3D model. The answer was to conduct the alignment on the x, y, z using procrustes analysis (<http://www.mathworks.com/help/stats/procrustes.html>) then convert that aligned structures to pdb using mat2pdb (<http://www.mathworks.com/matlabcentral/fileexchange/42957-read-and-write-pdb-files-using-matlab/content/mat2pdb.m)>. Finally the pdb formatted x,y,z coordinates could be viewed with molviewer from within matlab.
  1 Kommentar
OLUWATOSIN OLUWADARE
OLUWATOSIN OLUWADARE am 25 Aug. 2016
Bearbeitet: OLUWATOSIN OLUWADARE am 25 Aug. 2016
Can you please explain what you mean by "conduct the alignment on the x, y, z using procrustes analysis ". I know procustes function only take 2 arguments as input. I am also tring to do this but I get the error too. Thanks.

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Genomics and Next Generation Sequencing finden Sie in Help Center und File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by