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 - 07:00:03 PST