Dear Amber Community,
I need to quantify the total number of H-bonds for each pair of residues in my protein. That would mean a 1321 x 1321 matrix. I am having trouble in writing a bash script to process the CPPTRAJ hbond output (e.g. All.UU.avg.dat) to compute each matrix element. I’ve thought of running CPPTRAJ for every possible pair of donormask and acceptormask, but this is truly an ugly brute-force approach. Does anyone have suggestions on how I can compute the total number of H-bonds between each pair of residues?
Best,
Matthew
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Mar 16 2022 - 07:00:02 PDT
