[gmx-users] Membrane protein simulation
yz30 at nyu.edu
Sat Dec 22 21:07:59 CET 2007
> Dear Yanzi,
> Of course the system charge is 30+ so I have to add 30 Cl- . So I put “nn 30” to the command like this:
> genion -s protein_pope_em.tpr -o protein_pope_ion.pdb -nname CL- -nn 30 -g protein_pope.log
> But it seems only one water has been replaced with Cl-. Do I have to repeat this step for 30 times?!! or something wrong in this command.
Maybe something is wrong. Generally, you only need to do this for one
time. The command is right. Check your pdb file. One more thing, you
should add water first, then add ions.
> The protein is larger than box in the Z axis. I tried to add more water in the Z axis by genbox but it just added a few water in one side despite it needs more. Could you please tell me what I should do to fill the empty regions of the box with water. What do you mean about the pope may not be suitable for this protein and how do you understand it?
When adding water, before using genbox, you can use editconf at first.
For Gromacs, The origin of a box is always at (0, 0, 0). When you modify
the box size by genbox, the coordinates of the system wouldn't change.
By editconf, you can center molecule in box, translate or rotate the
molecule. When you use RASMOL to view the image, type "unitcell true",
you can see the box. Maybe VMD has the same function, but I just don't
What membrane is the protein found from? What does the membrane consist
of? Is the protein active in POPE? Read literature carefully, and choose
the suitable lipids.
> Thank you very much in advance for your kindness and expertize.
> Best regards,
> gmx-users mailing list gmx-users at gromacs.org
> 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 gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
More information about the gromacs.org_gmx-users