logo logo logo
Barnase-barstar koff calculation

This example (which is still in the testing phase) shows how to compute the dissociation rate constant for diffusional dissociation of barstar from a complex with barnase. For the input, it uses all pdb structures and grids precomputed in data_grid/. You need to run the script in prepare_grid_and_ecm/ at least once to generate the grids and effective charges in this folder.

Running the example

cd bnbs_koff/

Details of the script

This input will run a sda_koff computathe dissociation rate constant for diffusional dissociation of a solute from a complex. It is similar to a normal docking simulation, but the initial configuration should be read from a complexes file rather than starting randomly at a b-surface distance.
In addition, account can be taken of the occurrence of the configurations using the keyword account_occurency in the input file.

If a restart file is not provided, the simulation will start from the absolute position of solute 2, taken from the pdb input file.
The option is similar to the SDA 6 implementation of sda_koff

sda_koff computes the average time that the reaction criteria are satisfied before unbinding (in the file residence_time), and may be used to provide an estimation of k_off (results not validated).

Description of the results

The residence time is plotted on the figure below

  snap docking

At 6 Å, for 2 contacts, it is approximately 1.8 ns.
It is more than was predicted by SDA 6. This option is still in the testing phase.

[Back to Main Documentation]