One option may be considered if entropy is neglected:

=> psf2top can be use to convert psf file to gromacs format topology file. (https://github.com/resal81/PyTopol/wiki/psf2top-Usage)

=> catDCD can be use to convert dcd to gromacs format trr trrajectory file. (http://www.ks.uiuc.edu/Development/MDTools/catdcd/)

=>A GROMACS tpr file can be generated through grompp (A gromacs module)

=> Then g_mmpbsa can be use to calculate the binding energy through tpr and trr files. This tool does not include entropy calculation. (http://rashmikumari.github.io/g_mmpbsa/index.html)

I do not think you can do all necessary calculations using namd, you may need apbs for calculation of solvation energy, you may use vmd for mm calculations, and other software to do normal mode or quasi harmanic analysis for entropic calculations. In addition, you need other tools for sasa calculations.

But, if you convert your trajectory in a format acceptable for amber, you may use mm-pbsa tools of amber.

You need to download APBS software for calculating PB. APBS can be linked to NAMD during installation. A tutorial about how you can do it is found on installation documents for APBS and a tutorial about how to use APBS calculation with NAMD can be found on the NAMD tutorials.

One option may be considered if entropy is neglected:

=> psf2top can be use to convert psf file to gromacs format topology file. (https://github.com/resal81/PyTopol/wiki/psf2top-Usage)

=> catDCD can be use to convert dcd to gromacs format trr trrajectory file. (http://www.ks.uiuc.edu/Development/MDTools/catdcd/)

=>A GROMACS tpr file can be generated through grompp (A gromacs module)

=> Then g_mmpbsa can be use to calculate the binding energy through tpr and trr files. This tool does not include entropy calculation. (http://rashmikumari.github.io/g_mmpbsa/index.html)

Dear all,

I performed the g_mmpbsa analysis after great difficulty in installing the utility with the command:

g_mmpbsa -s md3.tpr -f md3.xtc

and the output was an energy file that I have attached herewith.

Is this the only output that the analysis gives?

I am apprehensive about this.

Kindly clarify.

@Revathi

The MM-PBSA method consists of three energy components: Molecular mechanics(MM), polar-solvation and non-polar solvation energy. You have calculated only MM energy.

To calculate other two components, you have to use option "-pbsa". You will need an input file for this calculation. For example:

g_mmpbsa -f md3.xtc -s md3.tpr -i mmpbsa.mdp -nomme -pbsa -pol -apol

"-nomme" suppresses the calculation of MM energy as you have already calculated it.

Details are provided in http://rashmikumari.github.io/g_mmpbsa/Usage.html#g_mmpbsa .

Keywords/parameters acceptable in mmpbsa.mdp file are described in http://rashmikumari.github.io/g_mmpbsa/Parameters.html

Once energy_MM.xvg, polar.xvg and apolar.xvg are obtained, "MmPbSaStat.py" script is required to calculate final binding energy.