Re: [AMBER] problem with extraction of potential energy for the ligand

From: Bill Ross <ross.cgl.ucsf.edu>
Date: Sun, 07 Jun 2015 14:20:28 -0700

Interestingly it is having a problem reading your lig_en_extract.in -
was there no other info after the stack trace to indicate what line of
the file is involved? Could there be Windows line endings in the file?

Aside from that issue, I see nothing that directs sander to read the
trajectory file - note that with -O plus -c lig.traj you could expect
the file to be truncated, if -c is still an output file specifier. I
believe sander may have an option to read an mdcrd and calculate
energies - note that it would require a prmtop for the ligand rather
than the complex, in case que.prmtop refers to complex below.

Bill

On 6/7/2015 10:14 AM, Maryam Hamzehee wrote:
> Dear Amber List
> I am trying to extract potential energy for my ligand from the ligand-receptor complex during the MD simulation (1ns). I have extracted the trajectory file of the ligand from ligand-receptor complex using ptraj, the input file is as follows:
> [maryam.cluster quercetin]$ more ptraj.in
> trajin md.binpos
> strip_mask :1-1330
> trajout lig.traj nobox
> #rms first out ligandrmsd.txt :1331
> #average ligaverage.pdb pdb :1331
>
> For extraction of potential energy I used the following input file:
> [maryam.cluster quercetin]$ more lig_en_extract.in
>
> &cntrl
> imin=5,maxcyc=1,
> cut=8.0,ntb=1,
> ntc=2,ntf=2,
> ntpr=0,
> /
> But the following error has been observed:
>
> [maryam.cluster quercetin]$ /share/apps/amber11/exe/sander -O -i lig_en_extract.in -o lig_en_extract.out -p que.prmtop -c lig.traj -y que.inpcrd
> forrtl: severe (64): input conversion error, unit -5, file Internal Formatted Read
> Image PC Routine Line Source
> libintlc.so.5 00002B1E31584AB1 Unknown Unknown Unknown
> libintlc.so.5 00002B1E31583A85 Unknown Unknown Unknown
> libifcore.so.5 00002B1E3107FB43 Unknown Unknown Unknown
> libifcore.so.5 00002B1E30FF4ADE Unknown Unknown Unknown
> libifcore.so.5 00002B1E30FF4286 Unknown Unknown Unknown
> libifcore.so.5 00002B1E31023576 Unknown Unknown Unknown
> libifcore.so.5 00002B1E31021085 Unknown Unknown Unknown
> sander 00000000004EFF34 Unknown Unknown Unknown
> sander 00000000004AF671 Unknown Unknown Unknown
> sander 000000000048C7B5 Unknown Unknown Unknown
> sander 000000000048C60C Unknown Unknown Unknown
> sander 0000000000402DEC Unknown Unknown Unknown
> libc.so.6 0000003C0021D994 Unknown Unknown Unknown
> sander 0000000000402CF9 Unknown Unknown Unknown
> How can I extract the potential energy for the ligand. I would be grateful if you could kindly help me in this regard.
>
> Any help in this regard would be highly appreciated.
> All the bestMaryam
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Jun 07 2015 - 14:30:02 PDT
Custom Search