Hi Matthew,
In order to do a perturbation from H --> NH2, you need to use 2 extra
dummy atoms in the initial (starting) structure. According to AMBER 8
Manual, this isnt possible. The reason is because of 'mixing' function
used in pp. 133. This function, (1-lambda)^k, does not have the first
derivative equal to zero at lambda = 0. As a result, it is not a smooth
transition function at the beginning (around lambda = 0). I have been
working on this problem for 1-2 month, and I have changed the source code
of amber, and replaced this function with something else, which will be a
smooth function at both ends; at lambda=0 and lambda=1. I have a test case
where I am transforming a methane molecule (CH4) to ethane (C2H6). As far
as now, the simulations are giving reasonable results. By that I mean the
following: I am doing a perturbation from eth --> meth and use the
trapezoidal rule to calculate the integral (pp.132l Eq.1). I do the same
thing for meth -> eth perturbation, and calculate the integral.
Theoretically, they should yield the same absolute value with a sign
difference. I am getting some expected results, but I am doing some other
simulations to see more about it.
I will write more about my results later, hopefully soon. Best,
On Mon, 25 Apr 2005 mrbroad.ilstu.edu wrote:
>
>
> hey guys~
>
> I am using AMBER 8 and trying to preform some FEP calculation on guanosine
> triphosphate and its analogs. I'm not sure however to make atoms appear.
> What I am specifically trying to do is change a (H) into an (NH2). I was
> guessing that I might have to make some dummy atoms, but I am not sure the
> proper way to do this. Any help is greatly appreciated.
>
> Matthew Broadus
>
> ------------------------------------------------------------
> Illinois State University Webmail https://webmail2.ilstu.edu
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
>
>
--
Ilyas Yildirim
---------------------------------------------------------------
- Department of Chemisty - -
- University of Rochester - -
- Hutchison Hall, # B10 - -
- Rochester, NY 14627-0216 - Ph.:(585) 275 67 66 (Office) -
- http://www.pas.rochester.edu/~yildirim/ -
---------------------------------------------------------------
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Mon Apr 25 2005 - 21:53:02 PDT
