Molecular study of the diffusional process of DMSO in double lipid bilayers

As a way to quantify the diffusion process of molecular compounds through biological membranes, we investigated in this study the dynamics of DMSO through an 1,2-Dipalmitoyl-sn-Glycero-3-Phosphocholine (DPPC) bilayer system. To properly account for the diffusion of DMSO due to a concentration gradient, a double DPPC bilayer was setup for our simulations. In such configuration, the aqueous phases can be explicitly associated with the extra and intracellular domains of the membrane, which is seldom the case in studies of single lipid bilayer due to the periodicity imposed by the simulations. DMSO molecules were initially contained in one of the aqueous phases (extracellular region) at a concentration of 5 wt.%. Molecular dynamics simulation was performed in this system for 95 ns at 350 K and 1 bar. The simulations showed that although many DMSO molecules penetrated the lipid bilayer, only about 10% of them crossed the bilayer to reach the other aqueous phase corresponding to the intracellular region of the membrane. The simulation time considered was insufficient to reach equilibrium of the DMSO concentration between the aqueous phases. However, the simulations provided sufficient information to estimate parameters to apply Fick's Law to model the diffusion process of the system. Using this model, we predicted that for the time considered in our simulation, the concentration of DMSO in the intracellular domain should have been about half of the actual value obtained. The model also predicted that equilibrium of the DMSO concentration in the system would be reached after about 2000 ns, approximately 20 times longer than the performed simulation.