Permeation across hydrated DPPC lipid bilayers: Simulation of the titrable amphiphilic drug valproic acid

Valproic acid is a short branched fatty acid used as an anticonvulsant drug whose therapeutic action has been proposed to arise from membrane-disordering properties. Static and kinetic properties of valproic acid interacting with fully hydrated dipalmitoyl phosphatidylcholine lipid bilayers are studied using molecular-dynamics simulations. We calculate spatially resolved free energy profiles and local diffusion coefficients using the distance between the bilayer and valproic acid respective centers-of-mass along the bilayer normal as reaction coordinate. To investigate the pH dependence, we calculate profiles for the neutral valproic acid as well as its water-soluble anionic conjugate base valproate. The local diffusion constants for valproate/valproic acid along the bilayer normal are found to be similar to10(-6) to 10(-5) cm(2) s(-1). Assuming protonation of valproic acid upon association with - or insertion into - the lipid bilayer, we calculate the permeation coefficient to be similar to2.0 10(-3) cm s(-1), consistent with recent experimental estimates of fast fatty acid transport. The ability of the lipid bilayer to sustain local defects such as water intrusions stresses the importance of going beyond mean field and taking into account correlation effects in theoretical descriptions of bilayer translocation processes.