MD trajectory hbond analysis
Hi everybody, I wonder is there's a way that Chimera produces hbond occupancy statistics from amber trajectory files. I can produce hbond information between a ligand and a receptor molecule frame by frame. As the ligand forms hbonds with different amino acid residues in the binding site during the trajectory, I'd like to work out which ones are the most frequently formed. Thanks in advance for your advice Best regards George
Hi George, There is a FindHBond tool and findhbond command -- when you are showing a trajectory in the MD Movie tool (could be from Amber or any of several other programs), you can use the "per-frame script" feature in that tool to run the calculation at each step and save the results. Please see these previous posts for more details on doing that: <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2010-January/004736.html> <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2010-January/004751.html> It sounds like you would need to save a file listing the hbonds at each step, and then do some of your own postprocessing or analyses on these output files to determine how many times each hbond is found. Currently there is no built-in function to do that part. I hope this helps, Elaine ---------- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco On Sep 13, 2010, at 7:02 AM, George Tzotzos wrote:
Hi everybody, I wonder is there's a way that Chimera produces hbond occupancy statistics from amber trajectory files. I can produce hbond information between a ligand and a receptor molecule frame by frame. As the ligand forms hbonds with different amino acid residues in the binding site during the trajectory, I'd like to work out which ones are the most frequently formed. Thanks in advance for your advice Best regards George
Also, this could be accomplished more directly by using an MD Movie per-frame Python script, like so: if not hasattr(chimera, 'hbCounts'): # each execution of script is in new namespace, # so store counts somewhere permanent chimera.hbCounts = {} from FindHBond import findHBonds, recDistSlop, recAngleSlop hbs = findHBonds([mdInfo['mol']], distSlop=recDistSlop, angleSlop=recAngleSlop, cacheDA=True) selected = set(chimera.selection.currentAtoms()) for hb in hbs: if len([da for da in hb if da in selected]) == 1: # exactly one end selected if hb[0] in selected: ra = hb[1] else: ra = hb[0] chimera.hbCounts[ra] = chimera.hbCounts.get(ra, 0) + 1 To use the above, you would delete all atoms except for the ligand and receptor from your trajectory, select the ligand, and then loop through your trajectory once. This will populate the chimera.hbCounts dictionary. To dig out the data from that dictionary you will need to open the IDLE environment (General Controls->IDLE) and type: for a, count in chimera.hbCounts.items(): print a, count If you are at all conversant with Python, it is also trivial to sort the counts or to write the hbCounts info to a file. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu On Sep 13, 2010, at 9:13 AM, Elaine Meng wrote:
Hi George, There is a FindHBond tool and findhbond command -- when you are showing a trajectory in the MD Movie tool (could be from Amber or any of several other programs), you can use the "per-frame script" feature in that tool to run the calculation at each step and save the results.
Please see these previous posts for more details on doing that: <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2010-January/004736.html
<http://plato.cgl.ucsf.edu/pipermail/chimera-users/2010-January/004751.html
It sounds like you would need to save a file listing the hbonds at each step, and then do some of your own postprocessing or analyses on these output files to determine how many times each hbond is found. Currently there is no built-in function to do that part.
I hope this helps, Elaine ---------- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco
On Sep 13, 2010, at 7:02 AM, George Tzotzos wrote:
Hi everybody, I wonder is there's a way that Chimera produces hbond occupancy statistics from amber trajectory files. I can produce hbond information between a ligand and a receptor molecule frame by frame. As the ligand forms hbonds with different amino acid residues in the binding site during the trajectory, I'd like to work out which ones are the most frequently formed. Thanks in advance for your advice Best regards George
_______________________________________________ Chimera-users mailing list Chimera-users@cgl.ucsf.edu http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users
participants (3)
-
Elaine Meng -
Eric Pettersen -
George Tzotzos