Dear Grpmac's Users!<br><br><br>Recently, I've done Protein-ligand tutorial consisted of simulation of the lysozyme with the ligand.<br><br>I have some questions about analysisng of the results, in particular I want to analyse all non-covalent interactions in that system<br>
<br><br>1) First off all I've tried to analyse potential H-bonds in that system. I've analysed very short 0.5nm trajectory since my simulation have not been completed yet so I've obtained<br><br>Found 256 donors and 481 acceptors<br>
Will do grid-seach on 18x18x13 grid, rcut=0.35<br>Reading frame 200 time 400.000 <br>Found 2 different hydrogen bonds in trajectory<br>Found 2 different atom-pairs within hydrogen bonding distance<br>Merging hbonds with Acceptor and Donor swapped<br>
256/256<br>- Reduced number of hbonds from 2 to 2<br>- Reduced number of distances from 2 to 2<br>HB lifetime = 4.03 ps<br><br>From that information I've found that 2 h-bonds have been emerged during 0.5 simulation with their lifetime 4ps.<br>
<br>Firstly, during analysis of the <b><tt>-num</tt></b> <tt><a href="http://manual.gromacs.org/current/online/xvg.html"> hbnum.xvg</a> file I have not found presence of the two above bonds (this graphs was consisted of the time on X and strange scale measured from 0 to 1 on Y). How that graphs might be interpreted?</tt> By the way what are the typical time of the arising of the h-bonds in such systems and how it could be dependent on the background factors ?<br>
<br>Than I've found that all hbonds are calculated in the scope of the cutt-off distance ( from above log I've found that such distance is the 0.35 isnt it? ) How I can define such cutoff rad? In my MDF file I've found only parametries for the cuttoffs for the Electrosatics as well as for the vdv contacts. Does the cut off for the H-bond was calculated on the combination of the above values?<br>
<br><br>2) Also I want to investigate posible contributions of the others non-covalent contacts in my protein-ligand system. What gromacs programs might be used for the calculation of the electrostatics and the hydrophobic contacts by the analogy of the ussage of the g_hbond?<br>
<br><br>Thanks for your help,<br><br>James<br><br>