The solubilization of a small unilamellar vesicle by bile salts is studied using Monte Carlo and Langevin dynamics simulation with an implicit solvent. By using coarse-grained models of dipalmitoylphosphatidylcholine and sodium cholate, it has been shown that the mechanism of transition from vesicle to mixedmicelle consists of the following stages: (a) separation of the bile salts from the solution phase and adsorption on the vesicle membrane, (b) enhanced pore formation with increasing bile salt concentration in the membrane, (c) detachment of worm-like mixed micelles from the pores, and (d) decay of the worm-like mixed micelles into smaller mixed micelles which have a reduced ratio of lipid to cholate molecules.The solubilization of lipid vesicles by bile salts is a crucial step in lipid digestion in vertebrates.1 There is a huge body of experimental research that focuses on the association behavior of bile salts and their biological relevance to lipid transport. Over the last decades, much experimental effort has been dedicated to find a detailed mechanism that can explain the effect of bile salts on the transition of lipid vesicles into mixed micelles. On the basis of these researches, transition mechanisms comprising different stages, such as adhesion and incorporation of bile salts into liposomes followed by the release of micelles and unfolding of liposomes into bilayer patches, have been proposed.
2The particular molecular architecture of bile salts, which consists of a convex hydrophobic side and a concave hydrophilic side, differs from the head/tail structure of common surfactants. As a consequence, the bile salts show a unique association behavior. They do not have a sharp critical micelle concentration, and the number of molecules in the aggregates is relatively small.Computer simulations can help to clarify the peculiarities of the vesicle-to-micelle transition induced by bile salts. Simulations exploring the formation of bile salt micelles have suggested several micelle structures which are in agreement with some features of the experimental data. The direct simulation of the vesicle-to-micelle transition is much more demanding. To our knowledge, no simulation of the solubilization of whole vesicles by bile salts is available. The main limitations are the number of molecules required for the vesicle to be sufficiently large, and the high bile salt concentration needed in the simulation box. In our simulations, we have used coarse-grained models with an implicit solvent. The solvent effect is included in the nonbonded interaction scheme. We have used two different simulation methods: Metropolis Monte Carlo and Langevin dynamics simulations. The lipid, dipalmitoylphosphatidylcholine (DPPC) was used for building the lipid vesicle, and sodium cholate was used as an example of the bile salts. The coarse-grained models of DPPC and sodium cholate are shown in Figure 1 (left) and Figure 2. The geometry of the coarse-grained lipid model has been adapted from the MARTINI model.3 All spherical segments in the...