findHBonds finds all H bonds in the given models (though the donors/acceptors can be restricted). The kind of restrictions that you can specify in the hbonds command are built on top of the findHBonds call. There is another function in the FIndHBonds module that can help you here: filterHBondsBySel from FindHBonds import findHBonds, filterHBondsBySel, recDistSlop, recAngleSlop allHBonds = findHBonds(chimera.selection.currentMolecules(), distSlop=recDistSlop, angleSlop=recAngleSlop) loopHBonds = filterHBondsBySel(allHBonds, loop, "cross") --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu On Jun 22, 2015, at 2:47 AM, Ahir Pushpanath <ahir29@gmail.com> wrote:
Dear All, I am confused about what input I need to provide to the findHbond command in the FindHBond midas module. This is the help() on the command.
findHBonds(models, intermodel=True, intramodel=True, donors=None, acceptors=None, distSlop=0.0, angleSlop=0.0, interSubmodel=False, cacheDA=False).
For example if I want to find all H bonds between a pre-selected group of residues in loop, saved as:
loop = chimera.selection.currentAtoms() and the rest of the protein molecule, what command would I have to run? I am assuming the output(such as the one given in the findHbonds runcommand) would be saved into a variable if I provide one?
Ahir
-- Dr. Ahir Pushpanath PhD. Senior Biologist, Johnson Matthey.
_______________________________________________ Chimera-users mailing list Chimera-users@cgl.ucsf.edu http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users