;;;;;; Martini lipid topology for C16:2-C18:2 LPE (I0PE), generated using:
; The Martini lipid itp generator version 0.6  Args are: -o martini_v2.0_IPE_01.itp -alname I0PE-alhead 'E P' -allink 'G G' -altail '- CDDC'
; WARNING: Lipids topology was generated following the Martini 2.0 guidelines but this specific lipid type might
;          not have been tested and should therefore be used with care. 
;
; Description:
;   A general model lysophosphatidylethanolamine (LPE) lipid corresponding to atomistic e.g. C18:2(9c,12c) linoleic acid tails.
; Parameterization:
;   This topology follows the standard Martini 2.0 lipid definitions and building block rules.
; Reference(s): 
;   S.J. Marrink, A.H. de Vries, A.E. Mark. Coarse grained model for semi-quantitative lipid simulations. JPC-B, 108:750-760, 
;   2004. doi:10.1021/jp036508g 
;   S.J. Marrink, H.J. Risselada, S. Yefimov, D.P. Tieleman, A.H. de Vries. The MARTINI force field: coarse grained model for 
;   biomolecular simulations. JPC-B, 111:7812-7824, 2007. doi:10.1021/jp071097f 
;   T.A. Wassenaar, H.I. Ingolfsson, R.A. Bockmann, D.P. Tieleman, S.J. Marrink. Computational lipidomics with insane: a versatile 
;   tool for generating custom membranes for molecular simulations. JCTC, 150410125128004, 2015. doi:10.1021/acs.jctc.5b00209
; Created: 2016.09.27
;
;@INSANE alhead=E P, allink=G G, altail=- CDDC, alname=IPE, charge=0.0
;@RESNTEST IPE==I0PEif: atoms[0]==NH3
;@BEADS NH3 PO4 GL1 GL2 C1B D2B D3B C4B 
;@BONDS NH3-PO4 PO4-GL1 GL1-GL2 GL2-C1B C1B-D2B D2B-D3B D3B-C4B 
;

[moleculetype]
; molname      nrexcl
  I0PE         1

[atoms]
; id 	type 	resnr 	residu 	atom 	cgnr 	charge
   1 	Qd 	 1 	I0PE	NH3 	 1 	1.0 	
   2 	Qa 	 1 	I0PE	PO4 	 2 	-1.0 	
   3 	P1 	 1 	I0PE	GL1 	 3 	0 	
   4 	Na 	 1 	I0PE	GL2 	 4 	0 	
   5 	C1 	 1 	I0PE	C1B 	 5 	0 	
   6 	C4 	 1 	I0PE	D2B 	 6 	0 	
   7 	C4 	 1 	I0PE	D3B 	 7 	0 	
   8 	C1 	 1 	I0PE	C4B 	 8 	0 	

[bonds]
;  i  j 	funct 	length 	force.c.
   1  2 	1 	0.47 	1250 	
   2  3 	1 	0.47 	1250 	
   3  4 	1 	0.37 	1250 	
   4  5 	1 	0.47 	1250 	
   5  6 	1 	0.47 	1250 	
   6  7 	1 	0.47 	1250 	
   7  8 	1 	0.47 	1250 	

[angles]
;  i  j  k 	funct 	angle 	force.c.
   2  3  4 	2 	180.0 	25.0 	
   3  4  5 	2 	180.0 	25.0 	
   4  5  6 	2 	180.0 	25.0 	
   5  6  7 	2 	100.0 	10.0 	
   6  7  8 	2 	120.0 	45.0 	


; Include Position restraint file
#ifdef POSRES
#include "I0PE_posre.itp"
#endif

