[gmx-users] Calculating bond energy

[Yulia Einav]

Erik and Mark, thank you for the reply.
Well, I suppose there is no direct measure to do. But I will definitely try to perform some indirect analysis, such as distances, angles, torsion, and maybe FEP.
Thanks!
 

________________________________

From: gmx-users-bounces at gromacs.org on behalf of Mark Abraham
Sent: Thu 5/29/2008 12:10
To: Discussion list for GROMACS users
Subject: Re: [gmx-users] Calculating bond energy



Yulia Einav wrote:
> Thank you for your reply, Mark.
> My protein has four disulfide bonds. I did in silico mutation of one of them and, then, MD simulations. Now I want to know whether breaking one bond influence the others. I can measure the distance between the centers of mass of two sulfur atoms, but is it sufficient to supply the information about bond strength of each disulfide? I guess I mean the bond dissociation energy (not the average one, but rather for each bond individually) of the remaining S-S bonds in the protein.

In principle, I guess you could measure such a bond dissociation free
energy with free energy perturbation (see
http://en.wikipedia.org/wiki/Free_energy_perturbation in the first
instance) but offhand I'd guess that a state with two bound sulfur atoms
is going to be too far removed from two -SH groups that are free to
dissociate for FEP methods to work in finite time.

It's definitely true that any S-S atom distance doesn't measure anything
relevant. Even if you haven't constrained your bond lengths, there's no
reason to suppose that the contribution to the total potential energy of
the "bond stretching" term has any strong relationship with bond strain
or such.

This is a good lesson in why one needs to design a computational study
with a specific goal in mind, and then to choose tools that can do the
job - e.g. measure that observable.

Mark
_______________________________________________
gmx-users mailing list    gmx-users at gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


-------------- next part --------------
A non-text attachment was scrubbed...
Name: winmail.dat
Type: application/ms-tnef
Size: 5128 bytes
Desc: not available
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20080529/98530832/attachment.bin>