На пути к валидированию модели четырехкомпонентной липидной мембраны
firstname.lastname@example.org, email@example.com МГУ им. Ломоносова, Биологический Факультет, Кафедра Биофизики
In order to construct molecular model of reaction center of purple bacteria Rhodobacter sphaeroides we try to reproduce it's native lipid membrane environment. It was shown that there are lipid-binding sites on RC's hydrophobic part, at least some of which bind particular lipid molecules and can influence RC functionality [1, 2]. Here we try to recreate angular lipid environment of RC, including in our simulation all four of it's predominant lipid molecules, namely phosphatydylglycerol (PG), cardiolipin (CL), phosphatydylethanolamine (PE) and phosphatydylcholine (PC).
This work is dedicated solely to modeling of phospolipid bilayers of complex composition, including negatively charged PG and CL. Other two lipids, PE and PC, are zwitterionic and we expect them to play more passive role in intramembrane processes. Our ultimate goal is to find a set of parameters, that would allow us to obtain a reasonable area per lipid (ApL) values in GROMACS OPLS-AA force field simulations of pure and mixed lipid bilayers of different composition. Experimental data, we use for validation of our models, we get from classical Langmuir experiments with lipid monolayers [3, 4, 5] as well as other, like NMR spectroscopy .
It was shown that free energy of lipid mixtures usually depends nonlinearly on bilayer composition. This could be partly explained by tendency of charged lipids to assemble in ordered domains, thus decreasing their effective concentration in mixture . Hydrogen bonds between ions and water molecules and negatively charged phosphorous groups of PG, CL may account for strong dependency of free energy on ion strength . Upon wide range of concentrations the same nonlinear behavior is observed even in mixtures of zwitterionic lipids like PE and PC, comprising the same set of acyl chains [5,6]. Acyl chains themselves, their length and presence of double bonds, influence on lipid bilayer thickness in a very complicated manner, with interactions of additional methylen groups simply adding to net Van-der-Waals forces but also changing order parameter. Although each of mentioned factors of interaction has it's own distinct nature, all of them affect macroscopic thermodynamic parameters, including ApL. By comparison of bicomponent lipid mixtures of different compositions we hope to find some general dependencies of ApL on molar fraction. Ideally we hope to find a way to combine those dependencies to get some insight upon how ApL behaves in ternary and quaternary mixtures of these lipids.
1. Roszak AW, Gardiner AT, et Al. Biochemistry. 2007;46(11):2909-2916. doi:10.1021/bi062154i.
2. Giustini M, Castelli F, et Al. J Phys Chem B. 2005;109(44):21187-21196. doi:10.1021/jp054104d.
3. Sennato S, Bordi F, et Al. J Phys Chem B. 2005;109(33):15950-15957. doi:10.1021/jp051893q
4. Roach C, Feller SE, et Al. C Biochemistry. 2004;43(20):6344-6351. doi:10.1021/bi049917r.
5. Poger D, Mark AE. Journal of Chemical Theory and Computation. 2009 Nov 23;6(1):325-36.
6. Thurmond RL, Dodd SW, Brown MF. Biophysical J . 1991 Jan;59(1):108. doi: 10.1016/S0006-3495(91)82203-2
7. Murzyn K, Róg T, et Al. Biophysical J. 2005;88(2):1091-1103. doi:10.1529/biophysj.104.048835.