Molecular dynamics simulation of the partitioning of benzocaine and phenytoin into a lipid bilayer

Molecular dynamics simulations were used to examine the partitioning behaviour of the local anaesthetic benzocaine and the anti-epileptic phenytoin into lipid bilayers, a factor that is critical to their mode of action. Free energy methods are used to quantify the thermodynamics of drug movement between water and octanol as well as for permeation across a POPC membrane. Both drugs are shown to favourably partition into the lipid bilayer from water and are likely to accumulate just inside the lipid headgroups where they may alter bilayer properties or interact with target proteins. Phenytoin experiences a large barrier to cross the centre of the bilayer due to less favourable energetic interactions in this less dense region of the bilayer. Remarkably, in our simulations both drugs are able to pull water into the bilayer, creating water chains that extend back to bulk, and which may modify the local bilayer properties. We find that the choice of atomic partial charges can have a significant impact on the quantitative results, meaning that careful validation of parameters for new drugs, such as performed here, should be performed prior to their use in biomolecular simulations. (C) 2013 Elsevier B.V. All rights reserved.