[gmx-users] Average area per lipid

[Stephen Dutz]

Hi all,

 

I'm using gromacs for MD simulations on membrane proteins.  As part of
the analysis I would like to calculate the average area per lipid.  I've
been going about this by pulling the coordinates from the trajectory in
VMD, then importing these to Mathematica and using it to plot a voronoi
tessellation.  The problem I run in to is I have to discard ~10% of the
data (when averaging over multiple frames) due to Mathematica not being
able to complete the tessellation.  I believe this is occurring as a
result of PBC and Mathematica plotting concave polygons.  That being
said, I have just about given up on my method.

 

I was wondering if anyone out there has a script or methodology they
would be willing to share with me in order to calculate the area per
lipid with a protein present (i.e using a voronoi tessellation).  I
would appreciate it tremendously!

 

Sincerely,

 

Stephen Dutz

Department of Chemistry &

Center for Macromolecular Modeling and Materials Design 

California State Polytechnic University, Pomona

http://www.csupomona.edu/~smdutz

 

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20080228/5b96584b/attachment.html>