Also, I received an error when trying to run ante-MMPBSA.py to generate the new prmtop files to exclude residues 40 through 51 from my solvated complex.
This is the command I used:
ante-MMPBSA.py -p complex3solvated.prmtop -c newcomplex3.prmtop -r newReceptor3.prmtop -l newLigand3.prmtop -s :40-51 -n "!(:208)"
This is the error I received:
mtop -l newLigand3.prmtop -s :40-51 -n "!(:208)"
Error: parmed.py cannot be found!
I'm not sure if my syntax is wrong or if the issue is something else.
Thanks!
On Aug 23, 2019 11:22 AM, A K <adk194.outlook.com> wrote:
Hi,
I just wanted to ask this question again regarding how to run MMGBSA calculations to excluding certain residues since I haven't heard back.
I ran pmemd with a protein-ligand complex, but my prepared pdb still had a few residues missing in the middle of one of its two chains, so there's an unattached part of my protein that's moving around a lot when I view it in VMD and also affecting the backbone RMSD graphs I created with cpptraj. So, I'm trying to re-calculate the free binding energy since I think these unattached residues would probably be affecting the overall free binding energy calculation as well, and might make the default masks not accurate? My protein has three separated parts: residues 1 to 39 are chain A, residues 40 to 51 are chain B, and residues 52 to 207 are also chain B. There's around 6 missing residues between residues 51 and 52 in chain B that are not in my prepared complex, so residues 40 to 51 are moving around a lot when I view the VMD trajectory. My ligands are residue 208. How would I run MMPBSA.py to exclude residues 40 to 51? I'm not sure why the mask I tried isn't working.
The error I got is this:
Loading and checking parameter files for compatibility...
File "/amber14/bin/MMPBSA.py", line 94, in <module>
app.loadcheck_prmtops()
File "/amber14/bin/MMPBSA_mods/main.py", line 608, in loadcheck_prmtops
self.normal_system.Map(INPUT['receptor_mask'], INPUT['ligand_mask'])
File "/amber14/bin/MMPBSA_mods/parm_setup.py", line 130, in Map
"every atom in the complex topology!")
PrmtopError: provided receptor/ligand masks don't select every atom in the complex topology!
Exiting. All files have been retained.
tar: _MMPBSA*: Cannot stat: No such file or directory
tar: Exiting with failure status due to previous errors
rm: cannot remove `_MMPBSA*': No such file or directory
Job complete
My input file, mmGBSA.in, is below:
&general
receptor_mask=':1-39,52-207',
ligand_mask=':208',
keep_files=2,interval=1,
/
&gb
igb=${i}, saltcon=0.100
/
The command I used is "MMPBSA.py -O -i mmGBSA.in -o complex1_output.dat -sp complex1solvated.prmtop -cp complex1.prmtop -rp receptor.prmtop -lp ligand1.prmtop -eo complex1_output.csv -y prod1_*.mdcrd"
Would I need to recreate the prmtop files again to exclude the residues as well to run MMPBSA.py or is specifying the masks enough but my syntax is just wrong? Also, the new prmtop files excluding residues 40 to 51 could be created either through leap or ante-mmpbsa.py, right? Would recreating the prmtops lead to an incompatibility with the coordinate/trajectory files though? Since my mdcrd files had these residues included when I ran the MD. I just want to make sure that my estimated binding energies would still be valid.
Thanks!
_______________________________________________
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 Fri Aug 23 2019 - 14:30:02 PDT
