Thank you, Dr. Scwhieters. I'm yet to get this working though. If you would be
kind enough to show it from the loading of a PDB file to calculating most
pairwise interaction energies (intra-chain ) primarily HBond, VdW & may be
electrostatics, it would be great. As am new to applying Xplor for such
purposes, this will be a good start. Moreover, as I found out there are not
great many tools that does this simple energy evaluation that is so regularly
done and hence this example could serve a general purpose as well.
Thanks again for your time.
On Tuesday, April 15, 2014 1:34 AM, Charles Schwieters <char...@schwieters.org>
wrote:
Hello Nahsivar--
> How to effectively calculate particular (pairwise) energy terms such as VdW,
> H-bond & electrostatic energies
> for a given model/conformation or even in general a PDB structure (a chain
> from a simple PDB file) ?
>
> Using MD programs like Gromacs or CHARMM or AMBER one could do it. But since
> X-plor also includes the corresponding FF terms it would be possible to
> easily evaluate such pairwise terms with more control (in the scripting)
>
> Simply, for interacting residue/atomic pairs, how to evaluate corresponding
> energies (from a given FF) particularly H-Bond energies & VdW (L-J potential
> energies). If this is not straightforward, where does one find the
> appropriate parameters ?
>
You can accomplish this task using the XPLOR CONStraints INTEraction
statement. From Python this would be e.g.
xplor.command("constr inte (name N and resid 1) (name HN and resid 1) end")
xplor.command("flags exclude * include BOND end")
bondEnergy = xplor.command("energy end",("ENER",))
(or something like this)
I hope this helps--
Charles_______________________________________________
Xplor-nih mailing list
Xplor-nih@cake.cit.nih.gov
http://cake.cit.nih.gov/mailman/listinfo/xplor-nih
