Department of Biochemistry, Virginia Tech
This tutorial will guide the user through the process of setting up and running pulling simulations necessary to calculate binding energy between two species. The tutorial assumes the user has already successfully completed the Lysozyme tutorial, some other tutorial, or is otherwise well-versed in basic GROMACS simulation methods and topology organization. Special attention will be paid to the methods for properly building the system and settings for the pull code itself.
The binding energy (ΔGbind) is derived from the potential of mean force (PMF), extracted from a series of umbrella sampling simulations. A series of initial configurations is generated, each corresponding to a location wherein the molecule of interest (generally referred to as a "ligand") is harmonically restrained at increasing center-of-mass (COM) distance from a reference molecule using an umbrella biasing potential. This restraint allows the ligand to sample the configurational space in a defined region along a reaction coordinate between it and its reference molecule or binding partner. The windows must allow for slight overlap of the ligand positions for proper reconstruction of the PMF curve.
The steps for such a procedure (and the ones utilized in this tutorial) are as follows:
The tutorial assumes that the reader is using GROMACS version 4.5.3 or later. My original work (from which this workflow was derived) was conducted with version 4.0.5, but in principle can be applied to any version in the 4.0.x or 4.5.x series. The pull code was completely re-written after version 3.3.3, such that none of the information contained herein (beyond the basic theory of the technique) is applicable to any GROMACS version prior to 4.0.
Bevan Lab Homepage
Virginia Tech Homepage
Virginia Tech Biochemistry
Site design copyright 2008-2012 by Justin Lemkul
Problems with the site? Send them to the Webmaster