Dear Amber users,
I am trying to calculate the intramolecular pairwise-residue energy
decomposition for a small protein. The protein is an apo monomer (no
ligand). I have tried the MMPBSA.py method as explained in the manual.
Since the system is not a complex, I used a single monomer for the complex
file with the command:
MMPBSA.py -O -i mmpbsa_decomp.in -cp receptor.top -y receptor.crd
In the output, the total energy for Res# n - Res# n is not null between
identical residues as I would be expected.
I tried also to create a 'complex' pdb by putting the same protein as both
receptor and ligand (superimposed) with exactly the same coordinates for
all the atoms. So I have 3 top files - com.top (2xMonomer), rec.top
(Monomer), lig.top (Monomer). In this case, the calculation exits with the
following error:
"Beginning GB calculations with /opt/amber12/bin/sander
calculating complex contribution...
CalcError: /opt/amber12/bin/sander failed with prmtop com.top! "
So I would like to know if it is possible to calculate, using MMPBSA.py,
the intramolecular pairwise-residue energy decomposition for each
residue-pair for a protein which is not complexed with any ligand or
another protein.
Thank you
Best regards
Rohit
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Sep 09 2013 - 01:00:02 PDT
