Example: Relative Binding Affinity

From AlchemistryWiki
Jump to navigation Jump to search



This example will be a bit more complex than the solvation example since we are now transforming one molecule to another within an environment. Again, the purpose of this page is not to take you through step-by-step the process of simulation, but instead guide you through the logical decisions that can be applied to any simulation package of your choosing. We will work through the same questions that were asked on the previous example as these should be asked for any free energy simulation.

For this particular example, we will look at the difference in the relative binding free energies of toluene and phenol in the apolar cavity of T4 Lysozyme.

Toluene in the pocket...
Doublearrow.png
... becomes phenol in the pocket.

What is the Thermodynamic Cycle?

Example of using a thermodynamic path to find which molecule has a stronger binding affinity. We could either remove both molecules from the pocket which is time prohibitive, or we could simply transform one molecule to another both in and out of pocket to find the difference that way.

If you have already read through the thermodynamic cycle page, then you will have already seen this picture on the left. For our simulation, we can define the thermodynamic cycle either by directly converting toluene to phenol (which would be very costly due to the large number of changing interactions and box size), or we can recognize that free energy is a state function and we can use this relation to solve for it as well:

[math]\displaystyle{ \Delta \Delta A_{\mathrm{bind}} = \Delta A_{\mathrm{bind}}^B - \Delta A_{\mathrm{bind}}^A = \Delta A_{A\rightarrow B}^{\mathrm{bound}} - \Delta A_{A\rightarrow B}^{\mathrm{unbound}} }[/math]

which does not require an excessively large box and can be done relatively simply by choosing the correct intermediate states to transform the ligands.

To provide the visuals for this particular system, please see the rendered images below; note that the solvent can be either implicit or explicit and the cycle will still be valid.

T4-Toluene.png A
[math]\displaystyle{ \Delta A_{\mathrm{bind}} }[/math]
Doublearrow.png
T4-lysozyme.png Plus sign.jpg Toluene.png
Doublearrowupdown.png Bound
[math]\displaystyle{ \displaystyle \Delta A_{A\rightarrow B} }[/math]
Doublearrowupdown.png Unound
[math]\displaystyle{ \displaystyle \Delta A_{A\rightarrow B} }[/math]
T4-Phenol.png B
[math]\displaystyle{ \Delta A_{\mathrm{bind}} }[/math]
Doublearrow.png
T4-lysozyme.png Plus sign.jpg Phenol.png

What are the End States?

What are the Intermediate States?

What Simulations to Run?

What Analysis do we Perform?

Catches, Traps, and Anything to Watch out for

References