30.08.2013 Views

Mesoscopic models of lipid bilayers and bilayers with embedded ...

Mesoscopic models of lipid bilayers and bilayers with embedded ...

Mesoscopic models of lipid bilayers and bilayers with embedded ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

6.2 Computational details 81<br />

For reference purpose we also report the interaction parameters between water <strong>and</strong><br />

<strong>lipid</strong>s. By choosing these three different types <strong>of</strong> solute molecules we can discriminate<br />

between the chemical affinity <strong>and</strong> the mechanic effects on the distribution <strong>of</strong><br />

the solutes in the bilayer.<br />

dbA dbH dbN<br />

a (db1,w) 25 80 25 a (w,w) 25<br />

a (db1,h) 25 80 35 a (w,h) 15<br />

a (db1,t) 50 25 25 a (w,t) 80<br />

a (db2,w) 80 80 25 a (h,h) 35<br />

a (db2,h) 80 80 25 a (h,t) 80<br />

a (db2,t) 25 25 25 a (t,t) 25<br />

a (db1,db1) 25 25 25<br />

a (db2,db2) 25 25 25<br />

a (db1,db2) 25 25 25<br />

Table 6.1: Repulsion parameters aij (for the force <strong>of</strong> equation F C ij = aij(1−rij/Rc)^rij) between<br />

dumb-bells, water, <strong>and</strong> <strong>lipid</strong>s. Water beads are indicated as w, <strong>lipid</strong> hydrophilic head-beads as<br />

h, <strong>lipid</strong> hydrophobic tail-beads as t, <strong>and</strong> the two beads <strong>of</strong> a dumb-bell as db1 <strong>and</strong> db2 respectively.<br />

Three types <strong>of</strong> dumb-bells are considered: hydrophobic (dbH), amphiphilic (dbA) <strong>and</strong><br />

neutral (dbN).<br />

To investigate the influence <strong>of</strong> the <strong>lipid</strong> characteristics on the distribution <strong>of</strong> the<br />

solutes, we consider three <strong>lipid</strong> types whose density distribution <strong>and</strong> lateral pressure<br />

pr<strong>of</strong>ile show different shapes. All the considered <strong>lipid</strong>s have a single tail <strong>of</strong> 7<br />

hydrophobic beads attached to one hydrophilic head-bead. The <strong>lipid</strong> tail can have<br />

different stiffness. We consider <strong>bilayers</strong> formed by fully flexible <strong>lipid</strong>s (ht7), by linear,<br />

t), <strong>and</strong> by poly-unsaturated <strong>lipid</strong>s, where the unsaturation points are<br />

stiff <strong>lipid</strong>s (ht (L)<br />

6<br />

represented by kinks along the chain (ht (L) t (K)<br />

4 t(L) t). The <strong>lipid</strong>s nomenclature follows<br />

the convention introduced in section 4.3 <strong>of</strong> Chapter 4, <strong>and</strong> the different <strong>lipid</strong> topologies<br />

used in this work are shown in figure 6.1. In all cases, the solutes were added<br />

to equilibrated, tensionless <strong>bilayers</strong> consisting <strong>of</strong> 200 <strong>lipid</strong>s (100 in each monolayer),<br />

<strong>and</strong> approximately 2000 water beads, at an overall bead density <strong>of</strong> 3. The dumb-bells<br />

mole fraction (respect to the number <strong>of</strong> <strong>lipid</strong>s) was 25% (the same used in [119]),<br />

which results in 50 dumb-bells for the <strong>bilayers</strong> considered here.<br />

The initial positions <strong>of</strong> the dumb-bells were chosen at r<strong>and</strong>om <strong>with</strong>in the hydrophobic<br />

part <strong>of</strong> the bilayer. The system was then equilibrated for 5000 DPD steps<br />

to allow the added molecules to diffuse in the bilayer. After this equilibration period<br />

the molecules were already located in the bilayer regions where they would reside<br />

throughout the rest <strong>of</strong> the simulation. Further 20000 MC-DPD cycles <strong>with</strong> imposed<br />

zero surface tension concluded the equilibration. Density <strong>and</strong> pressure pr<strong>of</strong>iles, <strong>and</strong><br />

bilayer structural characteristics were measured over 50000 MC-DPD cycles (at zero<br />

surface tension) <strong>with</strong> sampling frequency <strong>of</strong> 10 cycles.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!