Coarse-Grained Model for Phospholipid/Cholestrol Bilayer
No Thumbnail Available
Journal Title
Journal ISSN
Volume Title
Helsinki University of Technology |
Diplomityö
Checking the digitized thesis and permission for publishing
Instructions for the author
Instructions for the author
Author
Date
2004
Department
Major/Subject
Fysiikka (laskennallinen fysiikka)
Mcode
Tfy-105
Degree programme
Language
en
Pages
81
Series
Abstract
Biological membranes are complex systems with a large number of degrees of freedom. One of the main structural components of these membranes is a lipid bilayer. Computer simulations of lipid bilayers are very challenging: many interesting phenomena occur at much longer length and time scales than those accessible with atomic-level simulations. For studies beyond the microscopic regime, simplified coarse-grained models are therefore needed. We have constructed a simple, two-dimensional coarse-grained model for a lipid bilayer consisting of dipalmitoylphosphatidylcholine (DPPC) and cholesterol. The model has been constructed using a so-called Inverse Monte Carlo method, using results from atomic-level molecular dynamics simulations as input. The main goals have been to apply the Inverse Monte Carlo technique to lipid bilayers and find how it is suited for studying the large-scale properties of bilayers. In addition, we have explored the limitations of this approach. In our model each molecule is described by a single point-like particle. To find effective interactions for the coarse-grained particles, we have first performed detailed molecular dynamics simulations of the system at various cholesterol concentrations over a wide concentration range. We have then used the Inverse Monte Carlo method to construct effective interaction potentials for the particles. These potentials are constructed such that the radial distribution functions calculated from the coarse-grained system are identical to those calculated from the molecular dynamics simulations. Therefore the model gives correct results for structural properties at short length scales. Using the coarse-grained model, we have performed Monte Carlo simulations to study the large-scale structural properties of lipid bilayers. These results seem promising. We see evidence of formation of cholesterol-rich and cholesterol-poor domains at intermediate cholesterol concentrations. This is in agreement with the phase diagram of the system. \\lork is still in progress to further understand the nature of this organisation. We have also found that the model is unable to describe the liquid crystalline to gel transition of pure DPPC bilayers. Additionally, we have learnt that the effective potentials are rather specific to both the concentration and temperature in which they were determined. We have also performed Brownian dynamics simulations to study general relationships between different diffusion coefficients. From these studies we can conclude that our approach for constructing coarse-grained models can be a useful tool in the future. It allows relatively easy adjustment of the level of coarse-graining, and it can be applied to various systems without major modifications. Also, our model could be improved in several ways to account for e.g. the ordering of the hydrocarbon tails of the DPPC molecules. Thus active research on the subject continues.Biologiset membraanit kuten solukalvot ovat monimutkaisia systeemejä. Yksi niiden keskeisimmistä rakenneosista on lipidikaksoiskerros. Lipidikaksoiskerrosten simulointi tietokoneella on haastavaa, sillä monet mielenkiintoiset ilmiöt tapahtuvat paljon suuremmassa mittakaavassa kuin mitä atomitason simulaatioilla voidaan saavuttaa. Tämän takia suuremman skaalan ilmiöiden tutkimiseen tarvitaan yksinkertaistettuja eli karkeistettuja malleja. Tässä työssä on muodostettu yksinkertainen, kaksiulotteinen malli dipalmitiinifosfatidyylikoliinista (DPPC) ja kolesterolista muodostuvalle lipidikaksoiskerrokselle. Mallin muodostamisessa on käytetty niin kutsuttua Inverse Monte Carlo (IMC) -tekniikkaa käyttäen hyväksi atomitason molekyylidynamiikkasimulaatioista laskettuja tuloksia. Tärkeimpinä tavoitteina on ollut IMC-tekniikan soveltaminen lipidikaksoiskerroksiin tämän lähestymistavan soveltuvuuden arvioimiseksi, kaksoiskerroksen suuren skaalan ominaisuuksien tutkiminen ja yksinkertaistetun mallin rajojen tutkiminen. Karkeistetussa mallissa jokainen molekyyli on kuvattu yksittäisellä pistemäisellä hiukkasella. Näiden hiukkasten väliset vuorovaikutukset on laskettu kahdessa vaiheessa: ensin on tehty atomitason molekyylidynamiikkasimulaatioita, joista hiukkasille on laskettu parikorrelaatiofunktiot. Tämän jälkeen on sovellettu IMC-tekniikkaa, joka etsii sellaiset vuorovaikutukset, joita käyttäen karkeistettu malli antaa samat parikorrelaatiofunktiot kuin molekyylidynamiikkasimulaatiot. Tämä takaa, että malli antaa oikeat rakenteelliset ominaisuudet pienessä skaalassa tarkasteltuna. Mallin muodostamisen lisäksi tässä työssä on tutkittu lipidikaksoiskerroksen suuren mittakaavan rakenteellisia ominaisuuksia käyttäen muodostettua karkeistettua mallia ja Monte Carlo -simulaatiotekniikkaa. Tulokset vaikuttavat lupaavilta. Pienillä ja suurilla kolesterolikonsentraatioilla kolesteroli vaikuttaa jakautuvan tasaisesti, mutta näiden konsentraatioiden välissä nähdään kolesterolirikkaiden alueiden muodostumista. Tämä on sopusoinnussa alkuperäisen systeemin faasikäyttäytymisen kanssa. Tutkimusta tämän järjestäytymisen luonteesta jatketaan edelleen. Simulaatioista havaitaan myös, että malli ei pysty kuvaamaan puhtaan DPPC-kaksoiskerroksen nestekiteisestä faasista geelifaasiin tapahtuvaa transitiota. Lisäksi havaitaan, että lasketut vuorovaikutuspotentiaalit eivät sovellu erityisen hyvin simulointiin kaukana siitä konsentraatiosta ja lämpötilasta, jossa ne on määritetty. Tässä työssä on tutkittu myös lyhyesti eri diffuusiotekijöiden välisiä yhteyksiä. Tulosten perusteella käytetty lähestymistapa voi olla tulevaisuudessa hyödyllinen työkalu karkeistettujen mallien rakentamisessa. Karkeistuksen taso on helposti säädettävissä, ja lähestymistapaa voi soveltaa erilaisiin systeemeihin. Tässä työssä käytettyä mallia voisi myös parantaa monin tavoin, esimerkiksi sisällyttämällä DPPC-molekyylien järjestäytyminen malliin. Tutkimuksia aihepiiristä jatketaan edelleen.Description
Supervisor
Ala-Nissilä, TapioThesis advisor
Vattulainen, IlpoFalck, Emma
Keywords
coarse-grained model, karkeistettu malli, lipid bilayer, lipidikaksoiskerros, cholesterol, kolesteroli