Abstract
Using the potential of mean constrained force method, molecular dynamics simulations with atomistic details were performed to examine the partitioning and nature of interactions of two nonsteroidal antiinflammatory drugs, namely aspirin and ibuprofen, in bilayers of dipalmitoylphosphatidylcholine. Two charge states (neutral and anionic) of the drugs were simulated to understand the effect of protonation or pH on drug partitioning. Both drugs, irrespective of their charge state, were found to have high partition coefficients in the lipid bilayer from water. However, the values and trends of the free energy change and the location of the minima in the bilayer are seen to be influenced by the drug structure and charge state. In the context of the transport of the drugs through the bilayer, the charged forms were found to permeate fully hydrated in contrast to the neutral forms that permeate unhydrated.