I would suggest you use multiple analysis approach.
1. Firstly you can calculated distance among residues during the simulation
time to capture the shift. Say for example distance between CA of Asp120
and another neighbouring group which can capture the shift. Use CPPTRAJ
distance script for that.
2. Dump PDBs at different interval using CPPTRAJ script to check evolution
of Asp120 and its shift in real time. Use this script.
trajin md2_dry.mdcrd
autoimage
outtraj snapshot10.pdb onlyframes 10000 (this implies sanpshot10.pdb
corresponds to 10ns and you can dump at any interval if you save
trajectories in every 1ps)
3. Do a Residue Interaction Analysis otherwise known as RIN on snapshots
dumped. Say for example in 10ns time theere is connection between Asp120
and X21 and you get the edge value which signifies the power of connection.
Now in 30ns that connection value is less that signifies there is a shift
which will be backed up the distance calculation. Use Cytoscape and
RINAnalyzer plugin for that.
HTH!!
SB
On Wed, Jul 16, 2014 at 2:00 PM, anu chandra <anu80125.gmail.com> wrote:
> Dear amber users,
>
> I am working with MD simulation of protein-ligand system. From literature,
> I came to know that a residues, say ASP120, in the binding site shift about
> 5 angstroms during ligand binding. I would like to capture this shift from
> my simulation. Is there a quantitative analysis tool in PTRAJ/CPPTRAJ,
> which can explain the shift in angstroms itself?. I can try with evolution
> of dihedral angle about CA-CB plane, but it wont help me to see the extend
> of shift?. Waiting for your valuable reply
>
>
>
> Thanks in advance
>
> Sincerely
> Anu
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
Thanks & Regards;
Soumendranath Bhakat
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jul 16 2014 - 05:30:03 PDT
