AMBER Archive (2007)

Subject: AMBER: md potential energy

From: Francesco Pietra (chiendarret_at_yahoo.com)
Date: Mon Jul 23 2007 - 11:14:48 CDT


With a non polymeric large new molecule (parameters from antechamber/gaff) I
run (Amber9) a simple GB MD in analogy with 4.2 of tutorial B1 (igb=1 standard,
ie water).

The molecule represents the less strained conformer obtained from a
global-space conformational search with a molecular mechanics package (MM3-like
force field in vacuum). I.e., I am using a "discrete conformational sampling"
approach, though I hope also to find later the way to a suitable simulated
annealing with Amber9.

The potential energy soon becomes stable, just as in the xmgrace plot of at
4.3 of tutorial B1. Though, it starts at ca +30, attaining constancy at ca
+100.

What these positive values of the potential energy mean? Examination of the
mdcr file with VMD shows the molecule grossly in order. Nothing grossly wrong
occurred. I have not calculated rms.

Although all that is encouraging (final task: interaction with receptors), I am
worried by the positive values.

Thanks

francesco pietra

       
____________________________________________________________________________________
Take the Internet to Go: Yahoo!Go puts the Internet in your pocket: mail, news, photos & more.
http://mobile.yahoo.com/go?refer=1GNXIC
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu