- Posts: 10
Parameterizing new molecule
- Prunotto
- Topic Author
- Offline
- Fresh Boarder
I am trying to parameterize a new lipid, which is not currently present in the MARTINI database.
For doing this, I am following the general outline that is provided in www.cgmartini.nl/index.php/tutorials-gen...rzining-new-molecule , where the toluene is taken into consideration as general example.
At step 1) ("Generate atomistic reference data"), you guys say that for toluene, 1 ns of simulation would be enough for sampling, as it is a quite simple molecule.
My molecule is quite more complicated: specifically, it is a lipid with 4 tails (therefore, the general structure is similar to cardiolipins) and it is composed by 288 atoms (hydrogens included). The MARTINI model that I created and that I am trying to parameterize is composed by 28 MARTINI beads.
Based on these data, can you estimate how long my simulation should be to be able to sample the bonds and angles accurately enough? Do you think that 40-50 ns would be enough?
Thanks a lot for the great job that you are doing! All the best!
Alessio
Please Log in or Create an account to join the conversation.
- helgi
- Offline
- Administrator
- Posts: 54
As a side node, you need to beware that your AA simulation only represents that specific condition (lipid phase, other lipid types in your simulation, temperature, etc) so if you want to make a general Martini lipid parameters it is a good idea not to over-fit, try to use similarity to already parameterized lipids, and/or do reference simulations under different conditions.
Hope this helps,
Cheers,
- Helgi
Please Log in or Create an account to join the conversation.
- Prunotto
- Topic Author
- Offline
- Fresh Boarder
- Posts: 10
I hope I am not bothering you too much, but I would still have a couple of questions:
1) I see from your tutorial that the atomistic simulations from which you sample the states, have a sampling frequency of about 1 frame every 1 ps. Is this necessary, or I can decrease a bit this sampling frequency? In general, is there a recommended sampling frequency or it depends by multiple factors?
2) The lipid that I have to parametrize in MARTINI will eventually be used to create a membrane. In this case, which solvent should I use to simulate my all-atom lipid?
As a matter of fact, it is true that the solvent that I will use in the MARTINI simulations is water, but since I will create a membrane patch, the lipid will spend all the time into a hydrophobic environment (the membrane itself).
In other words, if I simulate the all-atom lipid into water, will this simulation be able to capture the hydrophobic behaviour of the lipid tails?
I would like to thank you again and I'm looking forward to receiving your reply.
All the best,
Alessio
Please Log in or Create an account to join the conversation.
- helgi
- Offline
- Administrator
- Posts: 54
2) This will depend. Normally you could do a short AA simulation of one lipid in water just to see the “flexibility” of the lipids under those conditions. But your main AA reference simulation would be under “normal” conditions (as close as you can to those you want to use the lipid/molecule in. If this is a simple bilayer forming lipid e.g. POPC you can simulate a pure POPC bilayer. If this is a none bilayer forming lipid or a lipid that normally is in low quantity I would rather have simulate it in a reference bilayer (e.g. 5-20% of you lipid in a POPC bilayer or a bilayer similar to the one you want to simulate at the CG level).
Cheers,
- Helgi
Please Log in or Create an account to join the conversation.
- yogi@martini
- Offline
- Junior Boarder
- Posts: 31
I think there is a typo error in this tutorial, the statement "Remember that S-beads only have a small radius compared to other S-beads and will thus have the same partitioning free energy (between water and octanol) as normal beads." is confusing, did you mean to say "Remember that S-beads only have a small radius compared to other beads and will thus have the same partitioning free energy (between water and octanol) as normal beads."
What exactly the physical reasoning behind S-beads having small radius.
Thanks in advance.
Please Log in or Create an account to join the conversation.
- riccardo
- Offline
- Platinum Boarder
Regarding the physical reasoning behind the S-beads, the have a smaller radius as they describe 3 or 2 non-hydrogen atoms instead of 4.
Please make sure to check also the section about ring particles in S. J. Marrink et al., J. Phys. Chem. B, Vol. 111, No. 27, 2007.
Please Log in or Create an account to join the conversation.
- yogi@martini
- Offline
- Junior Boarder
- Posts: 31
helgi wrote: 2) This will depend. Normally you could do a short AA simulation of one lipid in water just to see the “flexibility” of the lipids under those conditions. But your main AA reference simulation would be under “normal” conditions (as close as you can to those you want to use the lipid/molecule in. If this is a simple bilayer forming lipid e.g. POPC you can simulate a pure POPC bilayer. If this is a none bilayer forming lipid or a lipid that normally is in low quantity I would rather have simulate it in a reference bilayer (e.g. 5-20% of you lipid in a POPC bilayer or a bilayer similar to the one you want to simulate at the CG level).
Cheers,
- Helgi
Hi Helgi. Thanks for the reply. But I could not get what you want to say. I am rephrasing the question. Suppose I have Atomistic simulation data for bilayer in water. And I want to build CG model for the lipid which forms the given bilayer, actually my plan is once I get correct CG model, I will self assemble it and see whether is forms the similar bilayer. To get bond and angle distribution should I do atomistic simulation of single lipid in water or not? OR should I analyze each and every lipid in given atomistic bilayer to get bond angle distribution. OR something else we need to do.
Thank you. Looking forward to get reply from you.
-Yogi
Please Log in or Create an account to join the conversation.
- helgi
- Offline
- Administrator
- Posts: 54
If you only run an AA bilayer simulation and parameterize your CG lipids to fit those, matching all angles / bond / dihedrals closely, you can tune them to get really good match between the lipid properties at both resolutions.
Now e.g. you raise the temperature (or take one lipid out into water) then your AA lipids will behave differently due to their new environment but your CG lipids not so much as you have looked them in place / over fit them to that specific AA reference environment. So it all depends on what you want to do and there is always a tradeoff.
Hope this helps,
- Helgi
Please Log in or Create an account to join the conversation.