|
|
Joined: Mar 2006
Posts: 18
Forum Member
|
OP
Forum Member
Joined: Mar 2006
Posts: 18 |
I plan to put membrane, protein, water together and run simulations. I was suggested to run NAPzT simulaton (pressure is coupled only to z direction and the xy-corss area of the simulation box is a constant), because the force filed (c27) was optimized in NAPzT simulaton. I am wondering if i can run NPT simulation, where pressure is also coupled to xy direction with the newly developed charmm force field, such as c36? is c36 well developed so that i can use this force field just like c27? thank you.
Last edited by parachuter; 02/05/10 10:20 PM.
|
|
|
|
Joined: Dec 2005
Posts: 1,535
Forum Member
|
Forum Member
Joined: Dec 2005
Posts: 1,535 |
Ummm... there is a c36 force field for carbohydrates and the CHARMM general force field (CGenFF) was also released with c36, but you're only simulating "membrane, protein, water". Are you talking about the c32 lipid force field?
|
|
|
|
Joined: Sep 2003
Posts: 8,596 Likes: 12
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 8,596 Likes: 12 |
The C36 lipid force field, which is suitable for NPT simulations, has not been released yet; the manuscript should be submitted for publication in the next couple of weeks.
Rick Venable computational chemist
|
|
|
|
Joined: Mar 2006
Posts: 18
Forum Member
|
OP
Forum Member
Joined: Mar 2006
Posts: 18 |
so, only c36 force field can be used to do NPT simulations? and it is well developed for NPT simulations with different kinds of lipids and amino acid (as long as these molecules are included in the force field)? thanks. "The C36 lipid force field, which is suitable for NPT simulations, has not been released yet"
|
|
|
|
Joined: Sep 2003
Posts: 8,596 Likes: 12
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 8,596 Likes: 12 |
NPT simulations of a range of PC lipid bilayers stay close to their expt area/lipid values for extended periods with the new c36 force field.
With the earlier c27 (and c32) lipid force fields, NPT simulations tend to condense in the xy plane, so that the area/lipid shrinks. To compensate, one can use the NPgammmaT ensemble with an applied surface tension (gamma) of about 15 dyn/cm.
Rick Venable computational chemist
|
|
|
|
Joined: Jul 2004
Posts: 165
Forum Member
|
Forum Member
Joined: Jul 2004
Posts: 165 |
Dear Rick,
Is there any update for the C36 lipid force field?
Where can we get that?
Thank you very much
|
|
|
|
Joined: Sep 2003
Posts: 8,596 Likes: 12
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 8,596 Likes: 12 |
Since the paper is now available online (J. Phys. Chem. B, ASAP), the files are undergoing final review, and are about to be released on Alex MacKerell's site, maybe next week.
Rick Venable computational chemist
|
|
|
|
|
|