James Starlight wrote:
Dear Justin,

Indeed I didnt take into account square degree so all values presented in the tutorial exactly true so this was exactly my mistake. :(

After some iterations I've obtained value for S per lipid beetween 0.65-0.71 A^2. This is the example of the system wich I've obtained after solvation and futher minimization. http://www.sendspace.com/file/h6urhq It looks good. :) By the way I'd like to know about correct value for the S per protein ? Where I can obtain such information for my protein ? ( e.g for big proteins such receptors and channels )


Such information depends on whether there is experimental evidence for the protein you're looking for.

The error with GenBox was due to the big vdv R for the carbons. It's intresting that I couldnt increase this value mush as 1.2. But my finally system didnt consist of some water within the membrane so I supose that it's good.

Now my system on the stage of the equilibtation so I have not any questions yet but I only wounder to know about teqnique of the force fied preparation for the membrane protein simulation.

In this work I've done this step in accordance with the tutorial. On the last stage I've copied gromos53a6_lipid.ff directory in $GMXLIB so now I can use it in pdb2gmx for instance. But when I've tried to make gro and top files for my lipid.pdb via pdb2gmx I've obtained errors that Residue 'DPP' not found in residue topology database so I think that in addition to the operations wich I've already done I should also make some changes in the
/gromos53a6_lipid.ff/aminoacids.hdb
 ./gromos53a6_lipid.ff/aminoacids.n.tdb
 ./gromos53a6_lipid.ff/aminoacids.c.tdb
shouldn't I?


pdb2gmx is not suited for building topologies of multi-molecule, heterogeneous systems like solvated membranes. You may be able to get a topology, but it will be a redundant mess. Follow the protocol in the tutorial; it avoids these issues in the most streamlined way possible. If you need a topology of a new lipid, create an .rtp entry (if one does not exist already) and run pdb2gmx on a single lipid molecule. You can then #include the new .itp file in your system topology.

Finally, I've heard that there are already pre-build force fields for the membrane smulation. E.g I've heard that there is some modified CHARMM ff wich gave good results in simulation. Do you know something about this?


Various CHARMM parameter sets exist. CHARMM27 is built into the Gromacs 4.5.x series, and CHARMM36 is available from the User Contributions section of the website.

-Justin


Thank you for your help,

James


--
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to