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.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Mar 16 2008 - 06:07:23 PDT
