AMBER: PMF calculation

From: Jena M <>
Date: Mon, 25 Jun 2007 20:59:02 +0100

I am having some difficulties in creating the umbrella sampling input file. I want to calculate PMF of a ligand moving from 8 to 2 of a residue, keeping the protein residues within 5 angstroms of the ligand and of the residue frozen. I want to decrease the distance at 0.2 intervals. So, for instance for the first window do I have to set r2=r3=7.8? Here is my complete input file:
  imin = 1,
  ntmin = 0,
  maxcyc = 3000,
  ntb = 1,
  ntr = 1,
  restraintmask='!:566,667 <:5 | (!:566,667)',
  cut = 10
 &wt type='DUMPFREQ', istep1=10 /
 &wt type='END' /
#distance restrainst
&rst iat=662,667, r1=0.0, r2=7.8, r3=7.8, r4=100.0, rk2 =100., rk3 = 100.,
As far as I understood we have two ways to prepare the input files for each window/simulation. One is to use the configuration of the previous window as input for the next window. The second one is to use always the same input file for each window. Is this correct? Which of the methods is better?
Any help is highly appreciated
Play free games, earn tickets, get cool prizes! Join Live Search Club.
The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed Jun 27 2007 - 06:07:24 PDT
Custom Search