> There is currently no working version of QM/MMGBSA to calculate the
> binding
> free energy in AMBER, although efforts are currently underway to
> include
> this implementation. A way for you to get around this would be to setup
> the
> unsolvated complex, receptor, and ligand mdcrds yourself, and then run
> a
> QM/MM using sander and set imin=5 to do a post-processing analysis on
> all
> three of those trajectories. You can grep out the energies and do the
> statistics analysis yourself, also.
What do you believe is the problem with using QM/MMGBSA? The QM/MM model in
AMBER supports Generalized Born it should be a simple matter to do MMGBSA
with it. Will the results be good? I have no idea but it should work right
now. You might have to hack the perl scripts a bit but that shouldn't be too
hard.
If you know of more fundamental reasons why it wouldn't work I'd be happy to
hear them and take a look.
All the best
Ross
/\
\/
|\oss Walker
| Assistant Research Professor |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- ross.rosswalker.co.uk |
|
http://www.rosswalker.co.uk |
http://www.wmd-lab.org/ |
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.
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Mar 14 2010 - 11:30:03 PDT
