Mol2db

From DISI
Revision as of 19:01, 17 November 2009 by Rgc (talk | contribs) (saving realization about source of bugs in confhier.c)
Jump to navigation Jump to search

mol2db is a program developed to combine multimol2 files into the Flexibase format read by DOCK.

The purpose of this page is to document how mol2db works and any changes made to mol2db, and is not necessarily a guide to preparing files to run or running them through mol2db.

ui.c reads the inhier parameter file and sets variables

hiergen.c uses the multiple conformations to construct the hierarchy data structure by the following algorithm:

  1. find the first atom in the rigid component as the first atom with only 1 unique position
  2. expand through bonds to any atoms that have the same or fewer unique positions. bonded atoms with more unique positions are recursively added as branches
  3. for each branch, repeat the examining of bonded atoms and examine unique positions

this would appear to have flaws if the atom positions downstream overlap, especially since a distance tolerance is used to consider things overlapped.

confhier.c puts the conformations into the hierarchy. the potential problem here is that only one atom in each conformation in used to determine if the conformation differs from a previous one. when there are some atoms in a conformation that are the same while others are different, the results will be random based on atom ordering in the mol2 files, and can be incorrect.