Re: [AMBER] MMPBSA.py PB calculation not working

From: Ray Luo <rluo.uci.edu>
Date: Wed, 6 Apr 2022 10:13:25 -0700

Richard,

Okay, so something is wrong with the prmtop/inpcrd files. Could you share
with me one of the pb.out files after the failure?

Make sure you are running your job in a brand new folder without any old
files.

Ray
--
Ray Luo, Ph.D.
Professor of Structural Biology/Biochemistry/Biophysics,
Chemical and Materials Physics, Chemical and Biomolecular Engineering,
Biomedical Engineering, and Materials Science and Engineering
Department of Molecular Biology and Biochemistry
University of California, Irvine, CA 92697-3900
On Wed, Apr 6, 2022 at 2:45 AM Richard Kullmann <
Richard.Kullmann.mpikg.mpg.de> wrote:
> Hello again,
>
> I tried to run the test cases with my files. For that, I just replaced
> the top and inpcrd files with my corresponding files (and commented out
> the dacdif part). For GB, it ran through again without issue. For PB, I
> got the same error as before. So maybe there is something wrong with my
> files? What part could let GB run but not PB?
> An important note might be that I created the top files for MMPBSA AFTER
> I had already done the simulation, and I only added the PBRadii option
> then. Furthermore I use tip5p, but that should not be relevant for
> implicit solvent calculations, right? Finally I am also loading custom
> parameters for special residues, but this has also not been a problem
> before.
> Does any of this give you an idea of whats going on?
> Thank you very much again!
>
> Best regards,
>
> Richard
>
> On 4/5/22 20:42, Ray Luo wrote:
> > Great! Could you next use an input file in one of the test cases that is
> > closest to your own job without change and run your job again? Let's see
> > what happens.
> >
> > Ray
> > --
> > Ray Luo, Ph.D.
> > Professor of Structural Biology/Biochemistry/Biophysics,
> > Chemical and Materials Physics, Chemical and Biomolecular Engineering,
> > Biomedical Engineering, and Materials Science and Engineering
> > Department of Molecular Biology and Biochemistry
> > University of California, Irvine, CA 92697-3900
> >
> >
> > On Tue, Apr 5, 2022 at 8:38 AM Richard Kullmann <
> > Richard.Kullmann.mpikg.mpg.de> wrote:
> >
> >> Dear Ray,
> >>
> >> both GB and PB tests run without issues.
> >>
> >> Best,
> >>
> >> Richard
> >>
> >>
> >> -------- Forwarded Message --------
> >> Subject:        Re: [AMBER] MMPBSA.py PB calculation not working
> >> Date:   Tue, 5 Apr 2022 07:48:50 -0700
> >> From:   Ray Luo <rluo.uci.edu>
> >> To:     Richard Kullmann <richard.kullmann.mpikg.mpg.de>
> >>
> >>
> >>
> >> Richard,
> >>
> >> Do you see any problem running mmpbsa.py tests under
> >> $AMBERHOME/AmberTools/test/mmpbsa_py?
> >>
> >> All the best,
> >> Ray
> >> --
> >> Ray Luo, Ph.D.
> >> Professor of Structural Biology/Biochemistry/Biophysics,
> >> Chemical and Materials Physics, Chemical and Biomolecular Engineering,
> >> Biomedical Engineering, and Materials Science and Engineering
> >> Department of Molecular Biology and Biochemistry
> >> University of California, Irvine, CA 92697-3900
> >>
> >>
> >> On Tue, Apr 5, 2022 at 1:23 AM Richard Kullmann
> >> <Richard.Kullmann.mpikg.mpg.de <mailto:Richard.Kullmann.mpikg.mpg.de>>
> >> wrote:
> >>
> >>       Hello Ray,
> >>
> >>       I had actually tried that as well, and got exactly the same error
> >>       message...Considering that GB worked fine, I would also assume
> that
> >> the
> >>       input prmtop and inpcrd files are fine. Is that correct?
> >>       Best regards,
> >>
> >>       Richard
> >>
> >>       On 4/4/22 17:33, Ray Luo wrote:
> >>        > Richard,
> >>        >
> >>        > How about you provide at least one option, such as "ipb=2," for
> >>       the &pb
> >>        > section, as in the &gb section of the input?
> >>        >
> >>        > Ray
> >>        > --
> >>        > Ray Luo, Ph.D.
> >>        > Professor of Structural Biology/Biochemistry/Biophysics,
> >>        > Chemical and Materials Physics, Chemical and Biomolecular
> >>       Engineering,
> >>        > Biomedical Engineering, and Materials Science and Engineering
> >>        > Department of Molecular Biology and Biochemistry
> >>        > University of California, Irvine, CA 92697-3900
> >>        >
> >>        >
> >>        > On Mon, Apr 4, 2022 at 3:02 AM Richard Kullmann <
> >>        > Richard.Kullmann.mpikg.mpg.de
> >>       <mailto:Richard.Kullmann.mpikg.mpg.de>> wrote:
> >>        >
> >>        >> Hello,
> >>        >>
> >>        >> I am trying to run MMPBSA.py on my trajectories doing PB and
> GB
> >>        >> calculations. Here is the input file I used for trying out the
> >>       script:
> >>        >>
> >>        >> Input file for running PB and GB
> >>        >> &general
> >>        >>      verbose=1,
> >>        >> /
> >>        >> &gb
> >>        >>     igb=2,
> >>        >> /
> >>        >> &pb
> >>        >> /
> >>        >>
> >>        >> When calling it with a command of the form
> >>        >>
> >>        >> MMPBSA.py -O -i mmpbsa.in <http://mmpbsa.in> -cp com.top -rp
> >>       rec.top -lp lig.top -y
> >>        >> traj.crd  ,
> >>        >>
> >>        >> this is the final part of the output I get:
> >>        >>
> >>        >> Running calculations on normal system...
> >>        >>
> >>        >> Beginning GB calculations with
> >>        >> /cluster/apps/amber20/gnu/amber20_src/bin/mmpbsa_py_energy
> >>        >>     calculating complex contribution...
> >>        >>     calculating receptor contribution...
> >>        >>     calculating ligand contribution...
> >>        >>
> >>        >> Beginning PB calculations with
> >>        >> /cluster/apps/amber20/gnu/amber20_src/bin/mmpbsa_py_energy
> >>        >>     calculating complex contribution...
> >>        >>     File
> "/cluster/apps/amber20/gnu/amber20_src/bin/MMPBSA.py",
> >>       line 100,
> >>        >> in <module>
> >>        >>       app.run_mmpbsa()
> >>        >>     File
> >>        >>
> >>
> >>
> >>
> "/cluster/apps/amber20/gnu/amber20_src/lib/python3.8/site-packages/MMPBSA_mods/main.py",
> >>        >>
> >>        >> line 218, in run_mmpbsa
> >>        >>       self.calc_list.run(rank, self.stdout)
> >>        >>     File
> >>        >>
> >>
> >>
> >>
> "/cluster/apps/amber20/gnu/amber20_src/lib/python3.8/site-packages/MMPBSA_mods/calculation.py",
> >>        >>
> >>        >> line 82, in run
> >>        >>       calc.run(rank, stdout=stdout, stderr=stderr)
> >>        >>     File
> >>        >>
> >>
> >>
> >>
> "/cluster/apps/amber20/gnu/amber20_src/lib/python3.8/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'
> >>        >> Exiting. All files have been retained.
> >>        >>
> >>        >> It seems that GB runs smoothly but PB doesn't.
> >>        >> Do you have any idea what could be the issue here?
> >>        >> Thank you very much in advance and best regards,
> >>        >>
> >>        >> Richard
> >>        >>
> >>        >> _______________________________________________
> >>        >> AMBER mailing list
> >>        >> AMBER.ambermd.org <mailto:AMBER.ambermd.org>
> >>        >> http://lists.ambermd.org/mailman/listinfo/amber
> >>       <http://lists.ambermd.org/mailman/listinfo/amber>
> >>        >>
> >>        > _______________________________________________
> >>        > AMBER mailing list
> >>        > AMBER.ambermd.org <mailto:AMBER.ambermd.org>
> >>        > http://lists.ambermd.org/mailman/listinfo/amber
> >>       <http://lists.ambermd.org/mailman/listinfo/amber>
> >>
> >>
> >> _______________________________________________
> >> 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
>
> _______________________________________________
> 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 Wed Apr 06 2022 - 10:30:03 PDT
Custom Search