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 <charles at
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
<http://cake.cit.nih.gov/pipermail/xplor-nih/attachments/20140415/c4de5903/attachment.html>
