Normal Mode Entropy calculations¶

Requirements¶

In this case, gmx_MMPBSA requires:

-> Must be defined -- -> Optional, but recommended -- -> Optional

See a detailed list of all the flags in gmx_MMPBSA command line here

Command-line¶

That being said, once you are in the folder containing all files, the command-line will be as follows:

gmx_MMPBSA -O -i mmpbsa.in -cs com.tpr -ct com_traj.xtc -ci index.ndx -cg 3 4 -cp topol.top -o FINAL_RESULTS_MMPBSA.dat -eo FINAL_RESULTS_MMPBSA.csv

mpirun -np 2 gmx_MMPBSA -O -i mmpbsa.in -cs com.tpr -ct com_traj.xtc -ci index.ndx -cg 3 4 -cp topol.top -o FINAL_RESULTS_MMPBSA.dat -eo FINAL_RESULTS_MMPBSA.csv

gmx_MMPBSA_test -t 17

where the mmpbsa.in input file, is a text file containing the following lines:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28

Sample input file for entropy calculations (nmode)
This input file is meant to show only that gmx_MMPBSA works. Althought,
we tried to used the input files as recommended in the Amber manual,
some parameters have been changed to perform more expensive calculations
in a reasonable amount of time. Feel free to change the parameters 
according to what is better for your system.

&general
sys_name="NMode",
startframe=1,
endframe=10,
forcefields="leaprc.protein.ff14SB",
PBRadii=2,
/
&gb
igb=1, saltcon=0.150,
/
#note that nmode will use only a fraction of the no. of frames selected in 
#&general variable. This way, nmode will only process 8 frames (3rd to 10th frames)
#note also that some parameters have been change to perform the calculation 
#faster (maxcyc=5, drms=2). The typical values for these parameters 
#are (maxcyc=50000, drms=0.001)
&nmode
nmstartframe=3,
nmendframe=10,
nminterval=1,
maxcyc=5, drms=2,
/

## Considerations In this case, a single trajectory (ST) approximation is followed, which means the receptor and ligand (in this case, the ligand is also another protein) amber format topologies and trajectories will be obtained from that of the complex. To do so, an MD Structure+mass(db) file (com.tpr), an index file (index.ndx), a trajectory file (com_traj.xtc), and both the receptor and ligand group numbers in the index file (3 4) are needed. The mmpbsa.in input file will contain all the parameters needed for the MM/PB(GB)SA calculation. In this case, 10 frames are going to be used when performing the MM/PB(GB)SA calculation with the igb1 (GB-HCT) model and a salt concentration = 0.15M.

nmode will be used for estimating the entropic contribution, though it's way more expensive in computation as compared with IE or C2 methods. Note that in this case we have "relaxed" the convergence criteria to speed up the calculation.

A plain text output file with all the statistics (default: FINAL_RESULTS_MMPBSA.dat) and a CSV-format output file containing all energy terms for every frame in every calculation will be saved. The file name in '-eo' flag will be forced to end in [.csv] (FINAL_RESULTS_MMPBSA.csv in this case). This file is only written when specified on the command-line.