[Biopython] The PDBParser Permissive setting

Mark Livingstone livingstonemark at gmail.com
Tue Jul 17 01:49:37 UTC 2012


Hi Guys,

In my code I am experimenting with different ways of doing RMSD
calculations. I have code which in addition to normal CA based RMSD
can do (CA & CB) RMSD and also sidechain RMSD. On a perfect PDB file
this works well. Unfortunately, the curation I have is fairly average
/ poor in quality :-( and I only find out when one of the liberal
number of Try/Except blocks falls over.

I need a better way to find out sooner if a PDB file is missing data.

I am wondering therefore is for PDBParser I set Permissive=0, and
after setting the relevant models and chains etc, I did


wt_atoms = Bio.PDB.Selection.unfold_entities(wtc, 'A')

If this successfully works without throwing an Exception, can I assume
that this unfolded chain is perfect, or are there ways that I could
still be tripped up?

Alternatively, can anyone suggest code that I can employ in my
curation process that will give me a decent sanity check of PDB
quality, so I can get on writing experimental code - and not
Try/Except blocks :-(

Thanks in advance,

MarkL



More information about the Biopython mailing list