Re: [AMBER] Issue with MMGBSA Calculation for Protein–siRNA Complex

From: Bill Miller III via AMBER <amber.ambermd.org>
Date: Fri, 19 Dec 2025 12:33:30 -0600

When MMPBSA.py says it "mmpbsa_py_energy failed with prmtop", I typically
look in the mdout files created by MMPBSA.py for the error message. If
there are any other std out/err files, check those, too. But the error
messages shown are not very indicative of the cause of the error.

-Bill

On Fri, Dec 19, 2025 at 10:04 AM YASHIKA . via AMBER <amber.ambermd.org>
wrote:

> Here is command that I used:
> MMPBSA.py -O -i mmpbsa.in -o mmpbsa_result_wt.dat -cp
> complex_wt_nosolvent.parm7 -rp protein_wt.parm7 -lp sirna_wt.parm7 -y
> traj_wt.mdcrd
> I also tried with the full location of working directory like
> MMPBSA.py -O -i /home/username/Download/184-22/mmpbsa.in....
>
> On Fri, 19 Dec 2025, 20:01 Tomás Benjamín Cáceres Madariaga, <
> caceres.tomas.uc.cl> wrote:
>
> > Hi, may you share the command that you are using to call MMPBSA.py?
> >
> > Tomás Cáceres
> > Bioquímico UC
> > Asistente de Investigación | Laboratorio de Simulación Computacional y
> > Diseño de Fármacos
> > ------------------------------
> > *From:* YASHIKA . via AMBER <amber.ambermd.org>
> > *Sent:* Friday, December 19, 2025 5:30:26 AM
> > *To:* amber.ambermd.org <amber.ambermd.org>
> > *Subject:* [AMBER] Issue with MMGBSA Calculation for Protein–siRNA
> Complex
> >
> > Dear Team,
> >
> > I am performing MMGBSA calculations on a protein–siRNA complex and
> > encountering an error. Below are the details:
> >
> > *System and files:*
> >
> > -
> >
> > *Complex topology (complex_wt_nosolvent.parm7):* 8691 atoms, 502
> > residues, 3 molecules
> > -
> >
> > *Protein topology (protein_wt.parm7):* 7352 atoms, 460 residues, 1
> > molecule
> > -
> >
> > *siRNA topology (sirna_wt.parm7):* 1339 atoms, 42 residues, 2
> molecules
> > -
> >
> > *Trajectory (traj.mdcrd):* 8691 atoms, 502 residues, 3 molecules;
> > stripped of water and ions
> >
> > *MMGBSA input:*
> >
> > &general
> > startframe=1
> > endframe=5000
> > interval=5
> > verbose=2
> > keep_files=0
> > /
> >
> > &gb
> > igb=5
> > saltcon=0.150
> > radiopt=0
> > /
> > *Observed error:*
> > Loading and checking parameter files for compatibility...
> > cpptraj found! Using /home/leena/miniconda3/envs/amber_fix/bin/cpptraj
> > mmpbsa_py_energy found! Using
> > /home/leena/miniconda3/envs/amber_fix/bin/mmpbsa_py_energy
> > Preparing trajectories for simulation...
> > 1000 frames were processed by cpptraj for use in calculation.
> >
> > Running calculations on normal system...
> >
> > Beginning GB calculations with
> > /home/leena/miniconda3/envs/amber_fix/bin/mmpbsa_py_energy
> > calculating complex contribution...
> > File "/home/leena/miniconda3/envs/amber_fix/bin/MMPBSA.py", line 100,
> in
> > <module>
> > app.run_mmpbsa()
> > File
> >
> >
> "/home/leena/miniconda3/envs/amber_fix/lib/python3.10/site-packages/MMPBSA_mods/main.py",
> > line 224, in run_mmpbsa
> > self.calc_list.run(rank, self.stdout)
> > File
> >
> >
> "/home/leena/miniconda3/envs/amber_fix/lib/python3.10/site-packages/MMPBSA_mods/calculation.py",
> > line 82, in run
> > calc.run(rank, stdout=stdout, stderr=stderr)
> > File
> >
> >
> "/home/leena/miniconda3/envs/amber_fix/lib/python3.10/site-packages/MMPBSA_mods/calculation.py",
> > line 156, in run
> > raise CalcError('%s failed with prmtop %s!' % (self.program,
> > CalcError: /home/leena/miniconda3/envs/amber_fix/bin/mmpbsa_py_energy
> > failed with prmtop
> >
> >
> /home/leena/Downloads/scmv1/simulation/184_22_out/complex_wt_nosolvent.parm7!
> >
> > Fatal Error!
> > All files have been retained for your error investigation:
> > You should begin by examining the output files of the first failed
> > calculation.
> > Consult the "Temporary Files" subsection of the MMPBSA.py chapter in the
> > manual for file naming conventions.
> > Could anyone please advise on what might be causing this failure, or
> > suggest steps to resolve it?
> > Thank You
> > With regards
> > Yashika
> > NSUT, Delhi
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> >
> >
> https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=05%7C02%7Ccaceres.tomas%40uc.cl%7C6270c62ecee2411aedf708de3ed96722%7C5ff5d9faf83f4ac1a4d2eb48ea0a00d2%7C0%7C0%7C639017300624802474%7CUnknown%7CTWFpbGZsb3d8eyJFbXB0eU1hcGkiOnRydWUsIlYiOiIwLjAuMDAwMCIsIlAiOiJXaW4zMiIsIkFOIjoiTWFpbCIsIldUIjoyfQ%3D%3D%7C0%7C%7C%7C&sdata=fEr%2FOP63kJRwpma4qBY6Pjw2Bu4SKIE4TSKA%2Bp3plM4%3D&reserved=0
> > <http://lists.ambermd.org/mailman/listinfo/amber>
> >
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>


-- 
Bill Miller III
Associate Professor of Biochemistry
A.T. Still University
417-549-0952
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Dec 19 2025 - 11:00:02 PST
Custom Search