Minimize solvent
 &cntrl

  ntx    = 1,       irest  = 0,
  ntpr   = 10,

  ntb    = 1,
  cut    = 9.0,     nsnb   = 10,   

  imin   = 1,       

  ## run minimization 1000 cycles
  maxcyc = 1000,
  
  ## first steepest decent than conjugate gradient
  ntmin  = 1,

  ## reduce initial step size by 10
  dx0    = 0.001,

  ## cartesian restraints
  ntr=1,
  restraint_wt=500.0,
  restraintmask=' :1-%(lenres)i',
 &end

