AMBER: How to analyze the trajectory generated by LES simulation

From: pang zhao <zhao.pang.gmail.com>
Date: Fri, 2 Jun 2006 15:08:03 +0800

Dear All,
I performe a LES simulation for 500 ps with 5 copys of the ligand.
I first run 2 ns NPT MD with explicit solvent and then with use the restrt
file "2ns.restrt" and prmtop file "npt.prmtop"
to generate the two files for LES, les.top, les.crd. I successfully run 500
ps LES simulation.

But I could not know how to analyze the trajectory generated by LES
simulation. Because the LES trajectory is different
from the normal trajectory.

when I try to do analyze hydrogen bond using ptraj with the file les.top,
the error occurs: ERROR in readParm: ...failed to find SOLVENT_POINTERS

When I use npt.prmtop for the ptraj analysis, it works.
But I find the analysis result is too different from that with the first 2
ns NPT trajectory.

Such as the RMSD analysis for the LES trajectory, for the protein the C
alpha RMSD is too big.

the picutre about the RMSD analysis is attached, one is the C alpha RMSD,
the second is that of the ligand.

the ptraj input :
trajin lesnpt.mdcrd
center :1-343 origin mass
image origin center familiar
rms first out protein_les_rms.out :1-342.CA
rms first out lig_les_rms.out :343


For the hydrogen analysis, I could not find the hydrogen bond of the ligand
with the LES trajectory,
but I can find the two hydrogen bonds of the ligand with the 2 ns NPT
trajectory.



Because the LES trajectory is different from the normal trajectory. I think
there the LES analysis may be diffrent from
the Normal MD simulations. But How to correctly analyze the LES trajectory?

Can I use the LES trajectory for MM-PBSA?

Great thanks!



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

C_alpha_RMSD.jpg
(image/jpeg attachment: C_alpha_RMSD.jpg)

lig_rmsd.jpg
(image/jpeg attachment: lig_rmsd.jpg)

Received on Sun Jun 04 2006 - 06:07:14 PDT
Custom Search