Potential enhancements to 'hbond'

[drroe]

• nointrares keyword?
• bridging time series?
• Ignore UA - V -UA bridges with nointramol hbond
• allow bridges within one residue
• hbond bridging track atom numbers
• Add option to 'hbond' to only keep series above certain cutoff; Something like seriescut 0.2 which will keep all series with occupancy >= 0.2.
• Water-water hbond tracking.
• Record standard deviation for distance/angle.
• 'byresidue' option
• Be able to generate a sparse 2D matrix (by atom or by residue) containing fraction hydrogen bonds present for atom/residue pairs.
• Option to report/calculate hydrogen bonds based on H-A distance.
• Auto-detect solvent hydrogen bond sites.
• Option to print detected donors/acceptors to a file.
• Track bridging interactions of higher orders (e.g. U-V-V-U).

[hainm]

@drroe

per solvent bridge hbond ID:

354(11+12+) means solute 354 is bridge between res 11 and 12? is index 0-based or 1-based?
what does the extra '+' in 12+ mean?

[drroe]

is index 0-based or 1-based?

1-based.

what does the extra '+' in 12+ mean?

Nothing special - every residue number has a + after it.

[envirodug]

@drroe, for the bridgebyatom new feature in hbond (thanks for this!), is it possible to add a few new features (or is it already possible to do this)?

1. An option to only report heavy atoms in the output (e.g. instead of detailing hbonds to Lys HZ1, HZ2, and HZ3, just report to NZ).

2. Map the atom numbers to the residue number and name and print to the output file (e.g. including 1:GLU, 2:ASP, etc. as in the hbond calculation without specifying bridgebyatom in the options).

[drroe]

An option to only report heavy atoms in the output (e.g. instead of detailing hbonds to Lys HZ1, HZ2, and HZ3, just report to NZ).

You mean in the average output right? An option doesn't exist for doing this yet - I'll put it on the to-do list.

Map the atom numbers to the residue number and name and print to the output file (e.g. including 1:GLU, 2:ASP, etc. as in the hbond calculation without specifying bridgebyatom in the options).

Does the printatomnum keyword do what you want?