How to query for all the nitrogen coordinates in pdb file?

2 次查看(过去 30 天)
0 down vote favorite
I was trying to extract Nitrogen coordinates from ubiquitin protein. I have the 1UBQ.pdb file from http://rcsb.org/pdb/home/home.do website. I have done the following.
pdb1 ='/home/devanandt/Documents/VMD/1UBQ.pdb';
x=pdbread(pdb1)
y=x.Model.Atom
'y' variable gives 1x602 struct array with many fields including the co-ordinates X,Y,Z. There are 76 residues in this protein and so 76 nitrogens. How to extract (X,Y,Z) data separately to an array?

采纳的回答

DEVANAND
DEVANAND 2014-10-7
编辑:DEVANAND 2014-10-7
pdb1 ='/home/devanandt/Documents/VMD/1UBQ.pdb';
x=pdbread(pdb1);
atom={x.Model.Atom.AtomName};
n_i=find(strcmp(atom,'N')); % Find indices of atoms
X = [x.Model.Atom.X];
Y = [x.Model.Atom.Y];
Z = [x.Model.Atom.Z];
X_n = X(n_i); % X Y Z coordinates of N atoms
Y_n = Y(n_i);
Z_n = Z(n_i);

更多回答(0 个)

类别

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