[AMBER] MMPBSA.py.MPI TypeError when using &pb

From: Eiros Zamora, Juan <j.eiros-zamora14.imperial.ac.uk>
Date: Mon, 30 Apr 2018 17:49:48 +0000

Dear Amber Experts,

I’m trying to run parallel MMPBSA calculations on a ligand binding trajectory.

Here’s the input file:

Input file for running PB and GB
&general
    verbose=2,
    strip_mask=:Na+,Cl-,
    ligand_mask=:EGC,
    netcdf=1,
    interval=10,
    entropy=1
/
&pb
    istrng=0.15,
    fillratio=4.0
/
&nmode
    nmode_istrng=0.15,
/

And this is the command I issue:

mpirun -np 4 MMPBSA.py.MPI -O -i mmpbsa.in \
    -sp stripped.prmtop \
    -cp complex.prmtop -rp receptor.prmtop \
    -lp ligand.prmtop -eo energies.csv -y Production.nc

This is the output that I get:

Loading and checking parameter files for compatibility...
cpptraj found! Using /usr/local/amber16/bin/cpptraj
mmpbsa_py_energy found! Using /usr/local/amber16/bin/mmpbsa_py_energy
mmpbsa_py_nabnmode found! Using /usr/local/amber16/bin/mmpbsa_py_nabnmode
Preparing trajectories for simulation...
About to call autoimage
Autoimage should have been called
125 frames were processed by cpptraj for use in calculation.
125 frames were processed by cpptraj for nmode calculations.

Running calculations on normal system...

Beginning PB calculations with /usr/local/amber16/bin/mmpbsa_py_energy
  calculating complex contribution...
/usr/local/amber16/lib/python3.6/site-packages/MMPBSA_mods/main.py:604: UserWarning: ligand_mask overwritten with default

  warnings.warn('ligand_mask overwritten with default\n')
  File "/usr/local/amber16/bin/MMPBSA.py.MPI", line 100, in <module>
    app.run_mmpbsa()
  File "/usr/local/amber16/lib/python3.6/site-packages/MMPBSA_mods/main.py", line 218, in run_mmpbsa
    self.calc_list.run(rank, self.stdout)
  File "/usr/local/amber16/lib/python3.6/site-packages/MMPBSA_mods/calculation.py", line 82, in run
    calc.run(rank, stdout=stdout, stderr=stderr)
  File "/usr/local/amber16/lib/python3.6/site-packages/MMPBSA_mods/calculation.py", line 428, in run
    error_list = [s.strip() for s in out.split('\n')
TypeError: a bytes-like object is required, not 'str'
Error occured on rank 0.
Exiting. All files have been retained.
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

Am I doing anything wrong? What’s weird is that an MMGBSA calculation with the following inputs runs fine:

Input file for running PB and GB
&general
    verbose=2,
    strip_mask=:Na+,Cl-,
    ligand_mask=:EGC,
    netcdf=1,
    interval=10
/
&gb
    igb=5,
    saltcon=0.150
/

My machine is a macOS High Sierra v. 10.13.4.

Note: I have not updated to ambertools 18 since I’m using a manually modified and re-compiled from source version of the MMPBSA.py program that avoids calling the autoimage command: see this thread of mine from awhile ago (http://archive.ambermd.org/201702/0101.html) for an explanation.

Thanks for any help!

Juan
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Apr 30 2018 - 11:00:03 PDT
Custom Search