logo logo logo
Barnase barstar association rate calculation

This script uses all pdb structures and grids in data_grid/. You need to run at least one time the script in prepare_grid_and_ecm/ to generate the grids and effective charges in this folder.

If you want more information about this step, refer to the tutorial in examples/Tutorial_SDA7/.

Running the example

cd bnbs_assoc/

You can check that the output files are similar to the ones generated at HITS in bnbs_assoc_hits/

Details of the script

The script will run the calculation of the association rate between barnase and barstar in two different ways :

To draw the association rates you can use xmgrace : xmgrace -settype xydy rates_combined
To vizualise the film : vmd ../p2_noh.pdb trajectory1_apolar.dcd, and then load p1_noh.pdb (p1 never moves)

Description of the results

In SDA 7, rate constants are calculated using the Northrup, Allison and McCammon method (see Northrup, Allison McCammon, (1984), J. Chem. Phys.).
Simulations with only electrostatic terms reproduce the experimental association rate constant when the formation of 2 contacts (red curve) at ~6Å is used as the criterion for formation of an encounter complex.
Namely, three or four contacts are formed at the same rate with contact distances from 6 to 10 Å. 

This graph has been drawn from the results in bnbs_assoc_hits/, from the simulation with bootstrap.
It is only 2 500 trajectories. Usually, we run 25 000-100 000 trajectories for producing reliable results.

[Back to Main Documentation]