- Posts: 3
developping ornithin lipid parameters
- vanliefferinge
- Topic Author
- Offline
- Fresh Boarder
Less
More
5 years 6 months ago #7793
by vanliefferinge
developping ornithin lipid parameters was created by vanliefferinge
Hi,
I would like to run CG simulations with a ornithine lipid bilayer, but ornithin lipid is not in the Martini force field yet.
I wanted to developped it by analogy with existing parameters. Ornithine is quite similar to lysine (minus a CH2 ) and the carbon chain is classical.
In this picture I present the assignment of the beads type for ornithine lipids with lysine as reference.
https://www.dropbox.com/home/ornithin?preview=ornithine.jpg
I have mainly two questions:
1) concerning the side chain of ornithine. It is one CH2 shorter than the one for lysine. Will the C3 bead be affected? Should I change the type of the bead or remove it?
2) I'm not sure about the P1 bead. It is not a classical lipid junction between the polar head and the carbon chains, so by comparison I decide for a P1. Do you have any advises regarding this?
thank you in advance
François
I would like to run CG simulations with a ornithine lipid bilayer, but ornithin lipid is not in the Martini force field yet.
I wanted to developped it by analogy with existing parameters. Ornithine is quite similar to lysine (minus a CH2 ) and the carbon chain is classical.
In this picture I present the assignment of the beads type for ornithine lipids with lysine as reference.
https://www.dropbox.com/home/ornithin?preview=ornithine.jpg
I have mainly two questions:
1) concerning the side chain of ornithine. It is one CH2 shorter than the one for lysine. Will the C3 bead be affected? Should I change the type of the bead or remove it?
2) I'm not sure about the P1 bead. It is not a classical lipid junction between the polar head and the carbon chains, so by comparison I decide for a P1. Do you have any advises regarding this?
thank you in advance
François
Please Log in or Create an account to join the conversation.
- riccardo
- Offline
- Platinum Boarder
5 years 5 months ago #7821
by riccardo
Replied by riccardo on topic developping ornithin lipid parameters
Sorry for the late reply. I am not able to see the figure you linked via dropbox.
Maybe you can post it in some other way if you are still interested.
Maybe you can post it in some other way if you are still interested.
Please Log in or Create an account to join the conversation.
- vanliefferinge
- Topic Author
- Offline
- Fresh Boarder
Less
More
- Posts: 3
5 years 5 months ago #7833
by vanliefferinge
Replied by vanliefferinge on topic developping ornithin lipid parameters
Hello,
thank you for your reply.
You can find the picture in the following link:
ibb.co/k2Azaf
regards
francois
thank you for your reply.
You can find the picture in the following link:
ibb.co/k2Azaf
regards
francois
Please Log in or Create an account to join the conversation.
- peterkroon
- Offline
- Gold Boarder
Less
More
- Posts: 210
5 years 5 months ago #7834
by peterkroon
Replied by peterkroon on topic developping ornithin lipid parameters
Hey Francois,
As far as molecules and mappings go, this is a hard one. If you want to keep your mapping mostly as is, I'd drop the C3 bead. Also, I expect P1 to be too polar for a ketone group, but see the original Martini paper.
As an alternative, you can deviate from the Lysine mapping (it's not Lysine anyway, with the missing CH2 and all). In that case, I propose something like the following:
H2N-CH2-CH2-CH2;
HC-COOH;
HN-(C=O)-CH2;
rest.
This has the advantage that it keeps functional group together better, and it very nicely follows 4-1 atoms per bead. The disadvantage is that it looks less like Lysine, so you might need to do some extra validation.
I'm also interested to hear what Riccardo has to add to this :)
As far as molecules and mappings go, this is a hard one. If you want to keep your mapping mostly as is, I'd drop the C3 bead. Also, I expect P1 to be too polar for a ketone group, but see the original Martini paper.
As an alternative, you can deviate from the Lysine mapping (it's not Lysine anyway, with the missing CH2 and all). In that case, I propose something like the following:
H2N-CH2-CH2-CH2;
HC-COOH;
HN-(C=O)-CH2;
rest.
This has the advantage that it keeps functional group together better, and it very nicely follows 4-1 atoms per bead. The disadvantage is that it looks less like Lysine, so you might need to do some extra validation.
I'm also interested to hear what Riccardo has to add to this :)
Please Log in or Create an account to join the conversation.
- vanliefferinge
- Topic Author
- Offline
- Fresh Boarder
Less
More
- Posts: 3
5 years 5 months ago #7836
by vanliefferinge
Replied by vanliefferinge on topic developping ornithin lipid parameters
Thank you Peterkroon for your reply,
Your mapping look much better than mine but which type of bead would you assign to these 3 beads?
H2N-CH2-CH2-CH2;
HC-COOH;
HN-(C=O)-CH2;
thank you
Your mapping look much better than mine but which type of bead would you assign to these 3 beads?
H2N-CH2-CH2-CH2;
HC-COOH;
HN-(C=O)-CH2;
thank you
Please Log in or Create an account to join the conversation.
- riccardo
- Offline
- Platinum Boarder
5 years 5 months ago #7837
by riccardo
Replied by riccardo on topic developping ornithin lipid parameters
I would go for the mapping proposed by Peter:
H2N-CH2-CH2-CH2; (Qd as in Lys)
HC-COOH; (P3?)
HN-(C=O)-CH2; (P4? P5?)
rest.
as it nicely follows the 4-1 mapping scheme + keeps functional groups together.
Regarding the bead types, I put some guesses based on Table 3 of Marrink et al 2007 JPCB. Ideally, if you had partitioning data on the fragment you are attaching to the lipid chains you could validate these bead choices by computing free energies of transfer for the fragment and comparing them to the experimental ones. Another option is to look in the literature for other Martini lipids which have a similar head - say, e.g., without the Lys-like piece but with the other two groups - to check how other people mapped those groups (and perhaps you'll find already validated choices).
H2N-CH2-CH2-CH2; (Qd as in Lys)
HC-COOH; (P3?)
HN-(C=O)-CH2; (P4? P5?)
rest.
as it nicely follows the 4-1 mapping scheme + keeps functional groups together.
Regarding the bead types, I put some guesses based on Table 3 of Marrink et al 2007 JPCB. Ideally, if you had partitioning data on the fragment you are attaching to the lipid chains you could validate these bead choices by computing free energies of transfer for the fragment and comparing them to the experimental ones. Another option is to look in the literature for other Martini lipids which have a similar head - say, e.g., without the Lys-like piece but with the other two groups - to check how other people mapped those groups (and perhaps you'll find already validated choices).
Please Log in or Create an account to join the conversation.
- peterkroon
- Offline
- Gold Boarder
Less
More
- Posts: 210
5 years 5 months ago #7839
by peterkroon
Replied by peterkroon on topic developping ornithin lipid parameters
HC-COOH is going to be deprotonated in pH 7, so I'd make it a Qa as first guess.
Please Log in or Create an account to join the conversation.
Time to create page: 0.112 seconds