- Posts: 10
Something about bilayer spontaneous-assembly
- Ray
- Topic Author
- Offline
- Fresh Boarder
My problem is same to the previous topic "512 DPPG bilayer", but i didn't use the polarized water.
I‘m looking forward very much to your kind reply.
Best regards.
Please Log in or Create an account to join the conversation.
- djurre
- Offline
- Admin
- Posts: 272
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
Yes,I have tried that way.I made the 512 lipid bilayer by expanding a 128 bilayer in the x and y direction using "genconf".The problem is the bilayer is not a symmetrical system after 40ns equilibrium in which water is probably in one side and some lipids run out of the bilayer.What's the problem?
Please Log in or Create an account to join the conversation.
- djurre
- Offline
- Admin
- Posts: 272
In the original 128 DPPC bilayer there was no water in the bilayer. Are you solvating the bilayer after you have multiplied it? Because you shouldn't: the system already contains water.
Please Log in or Create an account to join the conversation.
- siewert
- Offline
- Admin
- Posts: 61
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
Firstly I will still show my deep gratitude for your concern to my problem.
I'm afraid I have not made myself clear. Can you see the figure attach to the email ? I mean the lipid bilayer is probably OK in my system except some lipids run out of the bilayer. The main question is the water. Water is not on either hand of bilayer . They are almost on one hand of the bilayer! There is no water in the bilayer.
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
Thank you very much!
I also suspect the periodic boundary conditions in the beginning. The periodic images is not good enough because there are some lipids out of bilayer. What do you mean by "shift the bilayer in the box"?
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
Please Log in or Create an account to join the conversation.
- richor
- Offline
- Fresh Boarder
- Posts: 19
I have a question here.
Would there be a artificial periodic if we multiply the 128 system to 512?
For larger systems like 6000,will the effect be smaller?
Appreciate for your reply!
Please Log in or Create an account to join the conversation.
- djurre
- Offline
- Admin
- Posts: 272
Please Log in or Create an account to join the conversation.
- richor
- Offline
- Fresh Boarder
- Posts: 19
As we are using the Periodic Boundary Conditions,I cannot figure out why the patterns will be
lost.I have looked through the Gromacs tutorial,no clue was found.
Can you help me?
Best regards!
Please Log in or Create an account to join the conversation.
- djurre
- Offline
- Admin
- Posts: 272
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
I have recently got some problems when I tried to establish a relative large lipid bilayer system. So could I ask you a favor?
I wanted to build the system including 1152 lipid molecules. But I failed to get a plane lipid bilayer. In the first case, there are 1152lipids and 20000 water,the bilayer is not a single plane,it is just like a irregular lipid system mixed up with two bilayers. I then added more water molecules to 60000. I surprised to get a cylindric bilayer. It so confused me. Could you please help me check this question whether there are some details I have not noticed.
I'm looking forward very much to your kind reply.
Best regards.
Sincerely yours,
Ray
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
I'm sorry to reply later. You still remember my "cool" lipid bilayer? I have done that by enlarging the system "128 lipids" with "genconf" to "1152 lipids" system as the method you mentioned. However, I still can't obtain a plain bilayer. I think the problem exisited after energy minimization. After energy minimization, the lipid molecules moved to the surface of the box, which might lead to a cylindric bilayer. What's your opinion?
Please Log in or Create an account to join the conversation.
- djurre
- Offline
- Admin
- Posts: 272
100 W
100 DPPC
100 W
100 DPPC
100 W
100 DPPC
100 W
100 DPPC
and should NOT look like:
400 W
400 DPPC
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
I sent three pictures: "image after -nbox 3 3 1", "image after energy minimization", "image after MD for 90ns" to your email box and hope you could help me check the question.
Thanks a lot for your help.
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
I have sent these files to your email box. Thanks again!
Please Log in or Create an account to join the conversation.
- djurre
- Offline
- Admin
- Posts: 272
As I wrote in my previous post, your .top files has the wrong order. If you look in your starting.gro there are first 128 DPPC molecules, then comes W and then again 128 DPPC molecules, etc. However, in your top file you start with all the DPPC (1152) and then all the W. That is wrong. Read my last post again and change it.
Please Log in or Create an account to join the conversation.
- Ray
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
Yes, it is the mistake in the top file.I was so careless about your suggestion in your last post. And again, I'm very grateful for your time and kind help!
Best,
Ray
Please Log in or Create an account to join the conversation.