Dear Mark and David,

Thanks for your reply. Dimension of 128-lipid bilayer system is 6.32x6.43x8.57nm, which leads to the projected area/lipid of ~63.5 A^2, and the box size is consistent during whole simulation time. In order to measure total SASA, I typed "g_sas -f **.xtc -s **.tpr -o **.xvg -solsize **" with different values of -solsize. Outputs are shown below. Here, my questions is 1) With bigger solvent probes, why does the total SASA decrease ? I thought that SASA should increase because bigger solvents will not insert into the lipid headgroups.

2) with -solsize 1.9, area/lipid can be obtained from 43.0275/2/64=33.6 A^2. How can this be smaller than the projected area/lipid, ~63.5 ?

For -solsize 1.5,   -----------------------------
# This file was created Tue Jul 10 10:25:24 2007
# by the following command:
# g_sas -f test.xtc -s md.tpr -o area0.xvg -solsize 1.5 -b 100000 -e 200000
#
# g_sas is part of G R O M A C S:
#
# Gnomes, ROck Monsters And Chili Sauce
#
@    title "Solvent Accessible Surface"
@    xaxis  label "Time (ps)"
@    yaxis  label "Area (nm\S2\N)"
@TYPE xy
@ view 0.15, 0.15, 0.75, 0.85
@ legend on
@ legend box on
@ legend loctype view
@ legend 0.78, 0.8
@ legend length 2
@ s0 legend "Hydrophobic"
@ s1 legend "Hydrophilic"
@ s2 legend "Total"
@ s3 legend "D Gsolv"
   100000           0     85.8878     85.8878           0
   100050           0     88.9415     88.9415           0
   100100           0     88.0002     88.0002           0
-----------------------------------------------

For, -solsize 1.7,
-----------------------------------------------
# This file was created Tue Jul 10 10:27:27 2007
# by the following command:
# g_sas -f test.xtc -s md.tpr -o area1.xvg -solsize 1.7 -b 100000 -e 200000
#
# g_sas is part of G R O M A C S:
#
# GRowing Old MAkes el Chrono Sweat
#
@    title "Solvent Accessible Surface"
@    xaxis  label "Time (ps)"
@    yaxis  label "Area (nm\S2\N)"
@TYPE xy
@ view 0.15, 0.15, 0.75, 0.85
@ legend on
@ legend box on
@ legend loctype view
@ legend 0.78, 0.8
@ legend length 2
@ s0 legend "Hydrophobic"
@ s1 legend "Hydrophilic"
@ s2 legend "Total"
@ s3 legend "D Gsolv"
   100000           0     82.7674     82.7674           0
   100050           0     84.0109     84.0109           0
   100100           0     85.2114     85.2114           0
-----------------------------------------

For -solsize 1.9
-----------------------------------------
This file was created Tue Jul 10 10:28:40 2007
# by the following command:
# g_sas -f test.xtc -s md.tpr -o area2.xvg -solsize 1.9 -b 100000 -e 200000
#
# g_sas is part of G R O M A C S:
#
# Good gRace! Old Maple Actually Chews Slate
#
@    title "Solvent Accessible Surface"
@    xaxis  label "Time (ps)"
@    yaxis  label "Area (nm\S2\N)"
@TYPE xy
@ view 0.15, 0.15, 0.75, 0.85
@ legend on
@ legend box on
@ legend loctype view
@ legend 0.78, 0.8
@ legend length 2
@ s0 legend "Hydrophobic"
@ s1 legend "Hydrophilic"
@ s2 legend "Total"
@ s3 legend "D Gsolv"
   100000           0     43.0275     43.0275           0
   100050           0      47.978      47.978           0
   100100           0     38.3633     38.3633           0
----------------------------------------------


Thanks,

best,
Hwankyu.

Message: 5
Date: Tue, 10 Jul 2007 15:55:29 +1000
From: Mark Abraham <[EMAIL PROTECTED]>
Subject: Re: [gmx-users] Re: Calculations of area per lipid using
       g_sas
To: Discussion list for GROMACS users <gmx-users@gromacs.org>
Message-ID: <[EMAIL PROTECTED]>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed

Hwankyu Lee wrote:
[Hide Quoted Text]
Dear Mark,

I just measured the projected XY-area (simply XY-size of box), and used
typical periodic boundary condition (In g_sas, default considers PBC).
In the small system (128 lipid bilayer), there is no undulation, so
projected area per lipid may be almost the same as non-projected area
per lipid.  So, total solvent accessible surface area (the third column)
should be always larger than the projected XY-area, but I observed that
total SASA become smaller than the projected XY-area at some point.
Could you give me some suggestions?

Show us your box dimensions at a point when you observe this to occur -
or preferably a relevant chunk of both g_sas and g_energy output.

Mark
_______________________________________________
gmx-users mailing list    gmx-users@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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to