How to use pdbsuperpose without .Sequence field?

조회 수: 1 (최근 30일)
Steven
Steven 2015년 12월 3일
편집: OLUWATOSIN OLUWADARE 2016년 8월 25일
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

답변 (2개)

Yuebin Zhou
Yuebin Zhou 2015년 12월 7일
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 2015년 12월 7일
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
OLUWATOSIN OLUWADARE
OLUWATOSIN OLUWADARE 2016년 8월 25일
편집: OLUWATOSIN OLUWADARE 2016년 8월 25일
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.

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Genomics and Next Generation Sequencing에 대해 자세히 알아보기

태그

Community Treasure Hunt

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

Start Hunting!

Translated by