Methodological issues in molecular dynamics (MD) simulations, such as the treatment of long-range electrostatic interactions or the type of pressure coupling, have important consequences for the equilibrium properties observed. We report a series of long (up to 150 ns) MD simulations of dipalmitoylphosphatidylcholine (DPPC) bilayers in which the methodology of simulation is systematically varied. Comparisons of simulations with truncation schemes, Ewald summations, and modified Coulomb interactions, either by shift functions or reaction field models, to describe long-range electrostatics point out the artifacts inherent in each of these methods and above all those of straight cutoff methods. We further show that bilayer properties are less sensitive to the details of the pressure-coupling algorithm and that an increased integration time step of 5 fs can be safely used in simulations of phosphatidylcholine lipid bilayers.

The compulsive nature of tobacco use is attributable to nicotine addiction. Nicotine is eliminated by metabolism through the cytochrome P450 2A6 (CYP2A6) enzyme in liver. Inhibition of CYP2A6 by chemical compounds may represent a potential supplement to anti-smoking therapy. The purpose of this study was to rationally design potent inhibitors of CYP2A6. 3D-QSAR models were constructed to find out which structural characteristics are important for inhibition potency. Specifically located hydrophobic and hydrogen donor features were found to affect inhibition potency. These features were used in virtual screening of over 60 000 compounds in the Maybridge chemical database. A total of 22 candidate molecules were selected and tested for inhibition potency. Four of these were potent and selective CYP2A6 inhibitors with IC 50 values lower than 1 mM. They represent novel structures of CYP2A6 inhibitors, especially N1-(4-fluorophenyl)cyclopropane-1-carboxamide. This compound can be used as a lead in the design of CYP2A6 inhibitor drugs to combat nicotine addiction.