[gmx-users] RE: Insertion of protein in lipid bilayer
graham.smith at biiuk.com
Tue Aug 2 22:40:23 CEST 2005
> I did follow the position restraint step in my run by restricting P8
> z direction. What happened was that the 2 popc layers flipped right
> bat in the md so that the head groups and water from each layer were
> each other while the hydrophobic tails from each layer were situated
> away from each other as possible. It was really quite bizzare and
> getting a few chuckles from it, I'd like to find out what I did wrong.
> is what I did:
> 1) got the 128 popc pdb and itp from Sir Tieleman, used
I think they haven't "really" done that, you've probably got an illusion
due to periodic boundary conditions (sketch your system replicated along
z and you'll see what I mean). The simulation you've got isn't wrong but
you may want to change the coord file and repeat it so things look more
normal. See note (c) below ...
> 2) followed the how-to html to create holes using msms.
> A few questions here:
> a) should I minimize the new membrane with hole?
It won't do any harm, but it might well work without it anyway.
> b) how critical is removing all the possible lipids that may clash
Not at all. The whole idea is that it will push the remaining ones out.
> c) I simply used vmd to manually put my proteins into the membrane
> saved the coordinate of the protein in a seperate file. Then I copied
> pasted the simulation box size from the membrane pdb file to the top
> protein file. Is this acceptable? What would be the most bullet-proof
Moving files around with vmd could perhaps be the source of the apparent
"flipping" of the lipids - e.g. if it has centered something at (0,0,0);
gromacs avoids -ve coordinates and wraps them round into the next
periodic image. I generally use editconf -center [rvec] or editconf -c
to move things around, then subsequently concatenate the pdb files and
view the result to check I've got it right.
> 3) grompp and mdrun using the following parameters.
> A few questions here too:
> a) what would typical values be for the tau_p and tau_t?
I use tau-t = 0.1 and tau-p = 1.0 with Berendsen coupling but I don't
think it's critical, especially in an unphysical simulation like this.
> 4) I havn't got to this step yet but after the hole has been generated
> successfully generated in the membrane, do I append the protein.pdb
> the end of the membrane_hole.pdb file and run a normal em/md to the
Yes. Well, I normally put the protein first, but the order of the
components (protein, lipid, water, ions if you have them) doesn't
matter, as long as they're the same in the pdb file and the topology.
Dr Graham R. Smith
Biosystems Informatics Institute
Newcastle upon Tyne NE1 4EP
The information contained in this e-mail and any files transmitted with it are confidential
and intended solely for the use of the individual or entity to whom they are addressed. If
you are not the intended recipient or have received this communication in error, please
notify the sender by e-mail or telephone (+44 (0)191 211 2560) and delete it from your
system and destroy any copies of it. Any copying, retransmission, distribution, action
taken, or omitted to be taken, in reliance upon it is strictly prohibited and may be unlawful.
Please note that whilst steps have been taken by BII to ensure that this e-mail and
attachments are virus free, neither BII nor the sender accepts responsibility for viruses
and the recipient should ensure that the e-mail and attachments (if any) are actually
More information about the gromacs.org_gmx-users