AMBER Archive (2008)

Subject: AMBER: ptraj, hbond and non-bonded contacts

From: Samuel Genheden (a03samge) (a03samge_at_student.his.se)
Date: Thu Mar 13 2008 - 09:54:03 CDT

Hello, Amber users
 
I've have modelled a protein complex using a docking algorithm and now I'm performing an MD simulation of the docked complex. From the docking algorithm I have got a list of possible hydrogen bonds and a list of possible non-bonded contacts (carbon-carbon interactions). Now I want to track these contacts during the MD simulation...
 
For the hydrogen contacts it is easy, I just use ptraj and the hbond facilities and it works fine.
 
But how about the non-bonded contacts. Is it possible to use the hbond facility even for these contacts? As I see is the only problem how to define the acceptor atom. The acceptor command requires an hydrogen and the heavy atom which it is bonded to. Now I have carbons contacts that I want to keep track of, so which bonded atom should I use? My idea is just to take the heavy atom previous in the chain, e.g. a CB for a CG and CA for the CB. Is this a plausible solution?
 
Best regards, Samuel
 

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu