next up previous contents index
Next: General parameters Up: Collective Variable-based Calculations (Colvars) Previous: Overview   Contents   Index


A crash course

Suppose that we want to run a steered MD experiment where a small molecule is pulled away from a protein binding site. In Colvars terms, this is done by applying a moving restraint to the distance between the two objects. The configuration will contain two blocks, one defining the distance variable (see section 9.3 and 9.3.1), and the other the moving harmonic restraint (9.5.4).


colvar {
  name dist
  distance {
    group1 { atomNumbersRange 42-55 }
    group2 {
      psfSegID PR
      atomNameResidueRange CA 15-30 }
    }
  }
}

harmonic {
  colvars dist
  forceConstant 20.0
  centers 4.          # initial distance
  targetCenters 15.   # final distance
  targetNumSteps 500000
}

Reading this input in plain English: the variable here named dist consists in a distance function between the centers of two groups: the ligand (atoms 42 to 55) and the alpha carbon atoms (CA) of residues 15 to 30 in the protein (segment name PR). The atom selection syntax is detailed in section 9.4.

To the ``dist'' variable, we apply a harmonic potential of force constant 20 kcal/mol/Å$ ^2$ , initially centered around a value of 4 Å, which will increase to 15 Å over 500,000 simulation steps.


next up previous contents index
Next: General parameters Up: Collective Variable-based Calculations (Colvars) Previous: Overview   Contents   Index
http://www.ks.uiuc.edu/Research/namd/