RE: AMBER: question about QM/MM in amber9

From: Ross Walker <>
Date: Thu, 25 May 2006 08:33:17 -0700

Hi Aixiao

> I have a question about the QM/MM module in Amber9, I
> don't know if QM/MM can define the QM center dynamically, did
> not find it in manual. For example, for receptor+ligand
> simulation in Waterbox, first I define ligand and one water
> near to it as QM center, but later, after running MD for tens
> or hundreds ps, this water molecule may go far from the
> ligand and have no interaction with it, while another water
> come to the ligand, so, is it possible to change QM center
> dynamically during MD? Based on distance from ligand or some
> other criteria instead of using residue number?

At present one cannot do this. Computationally it is not simple as there is
a very large amount of bookkeeping involved. This feature may be available
in Amber 10 but that is some ways off. You could try and work around the
problem in two ways, e.g. apply very very weak constraints on the QM water
molecule so that it stays in the vicinity of the QM region. Although this
may affect the overall dynamics. You could also script things to run say
25ps of simulation, have ptraj find you the atom numbers of the closest
water and then use these to write your input file for the 25ps. In this way
you could 'dynamically' alter the QM region in 25 ps steps. If you script it
carefully you should be able to make it entirely automatic. Again not ideal
but these are the only options at present.

All the best

|\oss Walker

| HPC Consultant and Staff Scientist |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- |
| | PGP Key available on request |

Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Fri May 26 2006 - 21:05:21 PDT
Custom Search