On 11/12/12 7:03 AM, jonas87 wrote:
Ok, what i've currently done:

I start with a pdb file of a single molecule Lipid_A and its topology and
put it in a box:

editconf -f Lipid_A.pdb -o Lipid_A-box.gro -bt triclinic -d 1.0

I want a total of 128 molecules of this lipid in my system so I use genbox
to add 127 more:

genbox -cp Lipid_A-box.gro -ci Lipid_A.pdb -p Lipid_A.top -o 128_Lipid_A.pdb
-nmol 127

This only ads 8 molecules, I'm going to guess because more don't fit in my
box. How do I deal with this? How can i know in advance how big i have to
make the box?


You would need to know the volume of an individual lipid in the configuration that you are supplying, plus a volume that will allow for any other species in the system (like water) to yield the desired composition of the system. The above series of commands yields a very small box, as you have found. Most 128-lipid membranes are on the order of 6-8 nm cubes, depending on the specific properties of the lipids being studied.

-Justin

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

Justin A. Lemkul, Ph.D.
Research Scientist
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