Dear Amber Users
I want to know. How we do the hydrogen bond analysis for Heme , NO like
hetroatoms. I ran the cpptraj command but it shows the hydrogen bond for
only residues not for the heme or NO. which is the part of my system.
Kindly guide If any tutorial is explained it.
Regards
script i used is this
# Hydrogen bond analysis with cpptraj
# Load topology and trajectory
parm noH_solv.prmtop
trajin noH_5.mdcrd
# All solute-solute and solute-solvent hydrogen bonds + water bridges.
hbond All out All.hbvtime.dat solventdonor :WAT solventacceptor :WAT.O \
avgout All.UU.avg.dat solvout All.UV.avg.dat bridgeout All.bridge.avg.dat
# Backbone hydrogen bonds only. Save time series data.
hbond Backbone :1-189.C,O,N,H avgout BB.avg.dat series uuseries bbhbond.gnu
# Create file containing only number of solute-solute and solute-solvent
# hydrogen bonds, and number of solute-solvent-solute bridges.
create nhbvtime.agr All[UU] Backbone[UU] All[UV] All[Bridge]
# Calculate RMSD to first frame.
rms BBrmsd :1-189.C,CA,N out BBrmsd.agr
# Run to actually find hydrogen bonds and generate data.
run
# Perform lifetime analysis on backbone hydrogen bonds
lifetime Backbone[solutehb] out backbone.lifetime.dat
runanalysis
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Nov 26 2017 - 06:30:02 PST