In order to obtain binding energies between two proteins I have run a dynamics of the complex and I would like to extract the component proteins in order to obtain the energy differences.
The script that I generated produces an error for a missing unit number.
I include the input and the ouput file.
Thanks in advance,
The input file is corrupted and has garbage text in it, so you'll need to clean it up or get a new copy of the file.
Dear Dr. Venable,
Thanks for your answer.
Anyhow I would to know if the commands I used in the script to separate the complex into its components are the appropiate ones-
It's not completely clear what you're trying to do, and why you want to separate the components. Binding energy can be difficult to compute; getting the potential energy of an interaction is fairly easy in CHARMM, and does not require separating the components. The difficult part is estimating the entropy contribution, the TdS term in the standard free energy expression. Part of the problem is sampling; most simulations are too short to explore a major fraction of the phase space of a binding interaction for two large, flexible molecules.
I am trying to estimate theses energies without the entropy component as a zero order approximation.
I suppose (in the rigid body approxumation) I could also select the atoms of both proteins and then ask for an interaction energy. This would avoid having to split the system and generating the energy difference a priori.
You'd probably need to run separate simulations of the two proteins isolated from each other, and include the solvation energies in the calculations at the least. Because of the likely conformational changes of surface groups (and more) involved in the binding, the entropy effects could be substantial.