AMBER: umbrella sampling for proton transfer

From: <cniu.ualberta.ca>
Date: Tue, 25 Mar 2008 18:30:33 -0600

Hi, all,

As a new user for amber, I am trying to run some umbrella sampling to
study the catalytic reaction of an enzyme using AMBER9. The reaction
involves a proton being transferred from S-H of a cys to N of a His.
Following is the input file:

&cntrl
   imin = 0,
   ntx = 7, irest = 1,
   ntb = 2, pres0 =1.0, ntp = 1,
   taup = 2,
   ntr = 0,
   cut = 20,
   ntc = 2, ntf = 2,
   tempi = 300.0, temp0 = 300.0,
   ntt = 3,
   gamma_ln = 1.0,
   nstlim = 5000, dt = 0.001,
   ntpr = 1, ntwx = 500, ntwr = 1000,
   nmropt = 1
  /
  &wt type = 'DUMPFREQ', istep1 = 10 /
  &wt type = 'END' /
  DISANG = s_h_n.RST
  DUMPAVE = s_h_n_vs_t.1

and s_h_n.RST is
&rst iat = 2238,2239, r1=0.0, r2=1.5, r3=1.5, r4=90.0, rk2=50., rk3=50., /
&rst iat = 615,2239, r1=0.0, r2=2.5, r3=2.5 r4=99.0, rk2=50., rk3=50., /

and part of result of s_h_n_vs_t.1 is
        0 1.336 2.705
       10 1.336 2.486
       20 1.336 2.407
       30 1.336 2.567
       40 1.336 2.481
       50 1.336 2.413
       60 1.336 2.494

The hydrogen atom in question is initially bonded to S . 1.3 ang, I
would like to simulate this proton leaving the sulfur atom and moving
closer to a nearby nitrogen (initially at 2.7 ang).

My questions are:
1. why the distance between S and H stays at the same value? I tried
to increase the value of rk2 and rk3 to something like 500, but S-H
still keeps at 1.336. I made a test run to increase the distance
between convalently bonded S and C, and the value of S-C distance did
change. So I guess there must be something special to deal with
hydrogen in Amber. Am I right? What should I do to deal with this
problem?
2. Am I right to use more than one restraints in s_h_n.RST to run
umbrella sampling (like reaction coordinate R=r1-r2)? In fact, my
"real" mechanism will likely involve more than 3 restraints.

Thanks in advance to anyone that will give any suggestion.


Amara



-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Apr 18 2008 - 21:14:03 PDT
Custom Search