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

From: YASHIKA . via AMBER <amber.ambermd.org>
Date: Fri, 19 Dec 2025 21:33:39 +0530

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
Received on Fri Dec 19 2025 - 08:30:02 PST
Custom Search