AMBER: Reproducing constant pH results !!

From: Pradipta Bandyopadhyay <>
Date: Thu, 27 Sep 2007 10:46:49 +0530

Dear All,

  I am trying to reproduce the results of constant pH simulation in the JCC paper (Mongan, Case, McCammon).
I am using AMBER 8. I tried the 1AKI structure of Lysozyme. I followed the protocol given in the paper (an input file is given later), except initial structure preparation by optimizing H-bond network was not done. All simulations were run for 1ns as done in the JCC paper.

The results I am getting is OK for some residues but very different for others. For instance, at pH=4 the pKas are

residue my result JCC paper (from Table 2)
ASP18 2.05 1.98
ASP48 0.93 1.94 ( I got only 8 MC transitions - perhaps this is one reason for this difference)
ASP87 1.98 3.12

One typical MD input file is like this

 md simulation
 ntt=1,temp0=300.0 tempi=300.0,tautp=2.0

My questions are (a) What is causing this large difference for some residues - one is the initial structure is not prepared
                              using optimizing H-bond network? Anything else?
                         (b) Can anybody give me a structure of 1AKI where H-bond network is optimized?
                          (c) Is there anything wrong I am doing in MD (as seen from the input file).
Any other help will be appreciated as well.



Dr. Pradipta Bandyopadhyay
Associate Professor
Centre for computational biology and bioinformatics
School of Information Technology
Jawaharlal Nehru University, New Delhi, INDIA

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Sun Sep 30 2007 - 06:07:20 PDT
Custom Search