Re: AW: Energy minimization of membrane

From: Chris Harrison (charris5_at_gmail.com)
Date: Thu Feb 09 2012 - 01:28:36 CST

And to answer your other question, using constantArea is one way to
maintain the Area per Lipid.

Best,
Chris

On Fri, Feb 3, 2012 at 8:53 AM, R. Charbel Maroun
<charbel.maroun_at_inserm.fr> wrote:
> Hi,
>
> No, I haven't run the equilibration in CHARMM 'cause I'm using NAMD. In
> C-GUI's step5_assembly.str, I see :
>
>  SET BOXTYPE  = RECT
>  SET XTLTYPE  = TETRAGONAL
>  SET A        = 85.3711
>  SET B        = 85.3711
>  SET C        = 108.095
>  SET ALPHA    = 90.0
>  SET BETA     = 90.0
>  SET GAMMA    = 90.0
>  SET ZCEN     = 0.0
>  SET NLIPTOP  = 94
>  SET NLIPBOT  = 94
>  SET NWATER   = 15800
>  SET NPWATER  = 3
>  SET POSID    = POT
>  SET NEGID    = CLA
>  SET NPOS     = 45
>  SET NNEG     = 45
>
> So, no problem for A, B, C. Is ZCEN the cellOrigin? If I don't use Charmm,
> where do I get the grid size ? I found the following in the namd wiki :
> "The PME grid size should be about one cell per Angstrom in each dimension,
> with factors of 2, 3, and 5 to be sure that the FFT is efficient (i.e.,
> order N log N complexity). Newer versions of NAMD will pick this
> automatically if you set "pmeGridSpacing 1.0", but the grid may change when
> a constant pressure simulation is restarted." So, something like 90, 90 and
> 111 should be suitable.
>
> On another hand, using
>
> useFlexibleCell       yes
> useConstantRatio      yes
>
> I'm afraid that the cell will get distorted and not respect anymore the
> starting experimental cross-sectional area (63A² for DPPC). Wouldn't it be
> better to use useConstant Area ?
>
> Cheers,
>
> CM
>
>
>
> Ajasja Ljubetič <ajasja.ljubetic_at_gmail.com> a écrit :
>
>
>> Hi,
>>
>>
>>> I used CHARMM-GUI for building my membrane system.
>>
>>
>>
>> Did you run the equlibration in CHARMM as well? If so, just open the
>> file step6.6_equilibration.str, where you have the box size (A, B, C) and
>> the FFT grid size information.
>> If you didn't run the equilibration in charm then look
>> into step5_assembly.str.
>>
>> Usually I ignore the grid used by charmm and just set PMEGridSpacing to
>> 1.0.
>>
>> You might also want to use
>>
>> useFlexibleCell       yes
>> useConstantRatio      yes
>>
>> Best regards,
>> Ajasja
>>
>
>
>
> ----------------------------------------------------------------
> This message was sent using IMP, the Internet Messaging Program.
>
>
>

-- 
Chris Harrison, Ph.D.
Theoretical and Computational Biophysics Group
NIH Resource for Macromolecular Modeling and Bioinformatics
Beckman Institute for Advanced Science and Technology
University of Illinois, 405 N. Mathews Ave., Urbana, IL 61801
char_at_ks.uiuc.edu                            Voice: 773-570-6078
http://www.ks.uiuc.edu/~char               Fax: 217-244-6078

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 05:24:43 CST