Hi Matt,

this was a bug with the dynamic_measures setting. It was fixed in PyMOL 1.7.0. 
So please upgrade PyMOL or read this old post for workarounds:

http://sourceforge.net/p/pymol/mailman/pymol-users/thread/501949e8.2080...@users.sourceforge.net/

Cheers,
  Thomas

On 27 Aug 2014, at 10:30, Matthew Baumgartner <mp...@pitt.edu> wrote:

> Hi,
> I ran into a bug with the find polar contacts function (distance with 
> mode=2) today. The problem arises when I was trying to find the hydrogen 
> bonds between a ligand and a receptor that had been processed with 
> reduce (http://kinemage.biochem.duke.edu/software/reduce.php). A 
> bug/feature of reduce is that when it adds hydrogens, it gives all of 
> the newly added hydrogens a 0 as the atom index in the pdb file.
> So what happens is when you try to find the polar contacts between the 
> ligand and the receptor, it draws a bond between a ligand atom and the 
> last hydrogen (which happens to be non-polar) which is 12 A away. I'm 
> guessing it is because pymol is using the atom index from the pdb.
> One possible alternative would be to index based on the set of the atom 
> id, atom name, residue name, chain id, and residue ID, but I don't know 
> how feasible that is.
> 
> The work around that I have found is to pass the receptor through babel 
> after processing with reduce. babel prot.pdb prot2.pdb , which 
> re-indexes the atom ids.
> 
> Matt Baumgartner

-- 
Thomas Holder
PyMOL Developer
Schrödinger, Inc.


------------------------------------------------------------------------------
Slashdot TV.  
Video for Nerds.  Stuff that matters.
http://tv.slashdot.org/
_______________________________________________
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Reply via email to