Previous Thread
Next Thread
Print Thread
Joined: Sep 2021
Posts: 5
D
Dani Offline OP
Forum Member
OP Offline
Forum Member
D
Joined: Sep 2021
Posts: 5
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,

Attached Images
error1.zip (16.16 KB, 356 downloads)
SHA1: 5fdca146c3156ad73a58bdbf762d3503a879e174
Joined: Sep 2003
Posts: 8,623
Likes: 24
rmv Offline
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 8,623
Likes: 24
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.


Rick Venable
computational chemist

Joined: Sep 2021
Posts: 5
D
Dani Offline OP
Forum Member
OP Offline
Forum Member
D
Joined: Sep 2021
Posts: 5
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-
Thanks

Joined: Sep 2003
Posts: 8,623
Likes: 24
rmv Offline
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 8,623
Likes: 24
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.


Rick Venable
computational chemist

Joined: Sep 2021
Posts: 5
D
Dani Offline OP
Forum Member
OP Offline
Forum Member
D
Joined: Sep 2021
Posts: 5
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.

Joined: Sep 2003
Posts: 8,623
Likes: 24
rmv Offline
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 8,623
Likes: 24
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.


Rick Venable
computational chemist


Moderated by  lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u1 Page Time: 0.011s Queries: 27 (0.007s) Memory: 0.7537 MB (Peak: 0.8273 MB) Data Comp: Off Server Time: 2022-12-08 16:01:08 UTC
Valid HTML 5 and Valid CSS