How to use pdbsuperpose without .Sequence field?

1 view (last 30 days)
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

Answers (2)

Yuebin Zhou
Yuebin Zhou on 7 Dec 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 on 7 Dec 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 Comment
OLUWATOSIN OLUWADARE
OLUWATOSIN OLUWADARE on 25 Aug 2016
Edited: OLUWATOSIN OLUWADARE on 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.

Sign in to comment.

Categories

Find more on Genomics and Next Generation Sequencing in Help Center and File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!