Rob Carlson wrote:

An Appeal to Gromacs Users,

I have been learning computer modeling for a while, on my own, and without any formal training. Having said that, I am a little confused about some concepts concerning free energies of binding. Therefore, I would like to propose a scenario and hopefully you can tell me if I am right, completely wrong, or on the proper track. Any response is sincerely appreciated. The scenario is as follows:

Suppose we have a ligand and a receptor. The ligand has 4 conformations in solution. The receptor has 3 conformations in solution. The complex (receptor+ligand) has 2 conformations in solution.

I will now define some terms. The binding free energy is G(bind). The complex free energy is G(complex). The ligand free energy is G(ligand). The receptor free energy is G(receptor).

Thus, G(bind) = G(complex) - G(receptor) - G(ligand)

To calculate G(ligand), I will use the following equation: G(ligand) = -RT ln Z(ligand), where Z(ligand) is the configuration integral of the ligand.

Now here is where I probably start misunderstanding everything.

In theory Z = integral over "r" of [ integral over "r" of ( EXP(-E(r)/RT) ] ........ where E(r) is the internal energy (often called U) of the considered species (ligand, receptor, or complex) plus the solvation energy (unless explicit water molecules are used). "r" is the atomic coordinates of the species under consideration.

You are about right but here Z would be the Partition function. The following simplification that would be really too easy if it was thrue. We all would be doing very different things. The idea is correct expect that you have to have the total partition function to be able to do so and the complete partition function implies that you sample the full
conformational space !!
In principle if you know the proportion of the different states you could weight them ubt in practice this will not
give a realistic result, or may be you are super lucky.

I am not a specialist of this so you could llok up some books or paper that would talk about binding affinity and molecular dynamics. Alan Mark, Wvan Gunsteren have done a lot of work on this subject.

So, in my case, since I have 4, 3, and 2 conformations for the ligand, receptor, and complex respectively, I am thinking that the above equation can be simplified.

For the case of the ligand, Z(ligand) = EXP( -E(r1) / RT ) + EXP( -E(r2) / RT ) + EXP( -E(r2) / RT ) + EXP( -E(r2) / RT ), where r1 to r4 are the coordinates of the 4 ligand conformations. In practice, this means that for each of the 4 ligand conformations I will run a 0 step dynamics simulation at a desired temperature to obtain the energy E of each conformation. The E is plugged back into that equation to get Z(ligand). The reason I run a 0 step dynamics is that I don't want the conformation to change and I need to get the kinetic and potential energies. A minimization would only give me potential energies. Solvation energies I can get from GB or PB. If I use explicit water molecules, however, I don't need to do this because the energy E obtained from the 0 step dynamics will contain the contributions from the explicit waters.

I will then perform similar calculations to get Z(complex) and Z(receptor). Once I get all the Z values I can get the individual G values and then the G(bind) from the equations above.

So, hopefully my understanding is sound. But as I mentioned above, any comments or criticisms is appreciated, as is any recommendations for books, publications, or webpages that can help me better understand these concepts. Thank you.

Best continuation
XAvier

--
----------------------------------
Xavier Periole - Ph.D.

Dept. of Biophysical Chemistry / MD Group Univ. of Groningen
Nijenborgh 4
9747 AG Groningen
The Netherlands

Tel: +31-503634329
Fax: +31-503634398
email: [EMAIL PROTECTED]
web-page: http://md.chem.rug.nl/~periole
----------------------------------

_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to