No BSD License  

Highlights from
PDBRead

4.4

4.4 | 5 ratings Rate this file 5 Downloads (last 30 days) File Size: 28.86 KB File ID: #1288

PDBRead

by Sachin Nikumbh

 

05 Feb 2002 (Updated 05 Feb 2002)

Read a PDB file into a MATLAB structure.

| Watch this File

File Information
Description

PDBRead is a function that will read a Protein Data Bank file into a MATLAB structure.

PDB_struct = PDBRead(filename) reads the file corresponding to filename and stores the information contained in this file in the PDB_struct. e.g. PDB_struct = PDBRead('PDBSilk.txt')

The file that is being read by PDBRead should be compatible to the PDB file format described at the RCSB (Research Collaboratory for Structural Bioinformatics) web site. For more information about this format, please visit the following URL:

http://www.rcsb.org/pdb/docs/format/pdbguide2.2/guide2.2_frame.html

The file "PDBSilk.txt" mentioned in the sample code above was obtained from the following web site:

http://cmm.info.nih.gov/modeling/pdb_at_a_glance.html

The current version of this function can recognize the standard record types mentioned at the above URL. The information stored in the output structure can be extracted by using the . (dot) operator similar to a normal structure in MATLAB. e.g PDB_struct.HEADER will give the information about the HEADER record type.
  
Each record type has format fields. For more information about the fields,please visit the above mentioned URL. Information for each field for a record type can be accessed by using the . (dot) operator. e.g. PDB_struct.HEADER.depDate gives the deposition date for the particular molecule structure. The field names used for a particular record type are according to the information contained at the RCSB web site.

MATLAB release MATLAB 6.1 (R12.1)
Tags for This File  
Everyone's Tags
Tags I've Applied
Add New Tags Please login to tag files.
Comments and Ratings (7)
28 Oct 2004 Ping-Hung Chen  
01 Feb 2005 Joe Cereal

well, i needed an mfile example of file io with structures.. and this fit the bill nicely.

30 Jul 2005 Oliver Lange

Syntax Error when i try to read a pdb-file:
Reference to non-existent field 'REMARK'

06 Aug 2006 ashiwn rao

excellent

16 Jan 2008 Rishi Mukhopadhyay

Possibly its an issue specific to my setup (Kubuntu linux, matlab 7), but when I run the program on the included PDBSilk.txt file, I get an error message:
??? Cell contents indices must be greater than 0.

Error in ==> PDBRead at 348
                    PDB_struct.SEQRES.ResChainDetail(PDB_struct.SEQRES.NoOfResChain).AminoAcids = strcat(PDB_struct.SEQRES.ResChainDetail(PDB_struct.SEQRES.NoOfResChain).AminoAcids,GetAminoAcids(

30 Apr 2008 Celso Araujo

I ger almost the same error as the last poster, under matlab 7.4.0 and linux:

>> p = PDBRead('PDBSilk.txt');
Extracting the information............
??? Subscript indices must either be real positive integers or logicals.
Error in ==> PDBRead at 348
                    PDB_struct.SEQRES.ResChainDetail(PDB_struct.SEQRES.NoOfResChain).AminoAcids = strcat(PDB_struct.SEQRES.ResChainDetail(PDB_struct.SEQRES.NoOfResChain).AminoAcids,GetAminoAcids(

15 Jun 2011 Laurent Guerin

There is a small bug in the code.
On line 341, replace :
if CurRes ~= PrevRes
by
if ~strcmp(CurRes,PrevRes)

Please login to add a comment or rating.
Tag Activity for this File
Tag Applied By Date/Time
biotech Sachin Nikumbh 22 Oct 2008 06:40:43
pharmaceutical Sachin Nikumbh 22 Oct 2008 06:40:43
protein data bank Sachin Nikumbh 22 Oct 2008 06:40:43
pdb Sachin Nikumbh 22 Oct 2008 06:40:43
rcsb Sachin Nikumbh 22 Oct 2008 06:40:43
bioinformatics Sachin Nikumbh 22 Oct 2008 06:40:43

Contact us at files@mathworks.com