Atomistic molecular dynamics (MD) simulations of druglike molecules embedded in lipid bilayers are of considerable interest as models for drug penetration and positioning in biological membranes. Here we analyze partitioning of coumarin in dioleoylphosphatidylcholine (DOPC) bilayer, based on both multiple, unbiased 3 μs MD simulations (total length) and free energy profiles along the bilayer normal calculated by biased MD simulations (∼7 μs in total). The convergences in time of free energy profiles calculated by both umbrella sampling and z-constraint techniques are thoroughly analyzed. Two sets of starting structures are also considered, one from unbiased MD simulation and the other from "pulling" coumarin along the bilayer normal. The structures obtained by pulling simulation contain water defects on the lipid bilayer surface, while those acquired from unbiased simulation have no membrane defects. The free energy profiles converge more rapidly when starting frames from unbiased simulations are used. In addition, z-constraint simulation leads to more rapid convergence than umbrella sampling, due to quicker relaxation of membrane defects. Furthermore, we show that the choice of RESP, PRODRG, or Mulliken charges considerably affects the resulting free energy profile of our model drug along the bilayer normal. We recommend using z-constraint biased MD simulations based on starting geometries acquired from unbiased MD simulations for efficient calculation of convergent free energy profiles of druglike molecules along bilayer normals. The calculation of free energy profile should start with an unbiased simulation, though the polar molecules might need a slow pulling afterward. Results obtained with the recommended simulation protocol agree well with available experimental data for two coumarin derivatives.

Regional Centre of Advanced Technologies and Materials Department of Physical Chemistry Faculty of Science Palacky University tr 17 listopadu 12 771 46 Olomouc Czech Republic

Zobrazit více v PubMed

Orsi M.; Essex J. W.. Passive Permeation Across Lipid Bilayers: a Literature Review. In Molecular Simulations and Biomembranes, 1st ed.; Sansom M. S. P., Biggin P. C., Eds.; Royal Society of Chemistry: 2010; pp 76–90.

Balaz S. Modeling Kinetics of Subcellular Disposition of Chemicals. Chem. Rev. 2009, 109, 1793–1899. PubMed PMC

Afri M.; Gottlieb H. E.; Frimer A. A. Superoxide Organic Chemistry within the Liposomal Bilayer, part II: a Correlation between Location and Chemistry. Free Radical Biol. Med. 2002, 32, 605–618. PubMed

Xiang T.-X.; Anderson B. D. Liposomal Drug Transport: a Molecular Perspective from Molecular Dynamics Simulations in Lipid Bilayers. Adv. Drug Delivery Rev. 2006, 58, 1357–1378. PubMed

Berka K.; Hendrychová T.; Anzenbacher P.; Otyepka M. Membrane Position of Ibuprofen Agrees with Suggested Access Path Entrance to Cytochrome P450 2C9 Active Site. J. Phys. Chem. A 2011, 115, 11248–11255. PubMed PMC

Alberts B.; Johnson A.; Lewis J.; Raff M.; Roberts K.; Walter P.. Molecular Biology of the Cell, 4th ed.; Garland Science: New York, 2002.

Cooper G. M.The Cell: A Molecular Approach, 2nd ed.; Sinauer Associates: Sunderland MA, 2000.

van Meer G.; Voelker D. R.; Feigenson G. W. Membrane Lipids: Where They Are and How They Behave. Nat. Rev. Mol. Cell Biol. 2008, 9, 112–124. PubMed PMC

Black S. D. Membrane Topology of the Mammalian P450 Cytochromes. FASEB J. 1992, 6, 680–685. PubMed

Ishii Y.; Takeda S.; Yamada H. Modulation of UDP-glucuronosyltransferase Activity by Protein-Protein Association. Drug Metab. Rev. 2010, 42, 145–158. PubMed

Wu N. A. N.; Palczewski K.; Mu D. J. Vertebrate Membrane Proteins: Structure, Function, and Insights from Biophysical Approaches. Pharmacol. Rev. 2008, 60, 43–78. PubMed PMC

Marrink S.-J.; Berendsen H. J. C. Simulation of Water Transport through a Lipid Membrane. J. Phys. Chem. 1994, 98, 4155–4168.

Neale C.; Bennett W. F. D.; Tieleman D. P.; Pomès R. Statistical Convergence of Equilibrium Properties in Simulations of Molecular Solutes Embedded in Lipid Bilayers. J. Chem. Theory Comput. 2011, 7, 4175–4188. PubMed

Orsi M.; Essex J. W. Permeability of Drugs and Hormones through a Lipid Bilayer: Insights from Dual-Resolution Molecular Dynamics. Soft Matter 2010, 6, 3797–3808.

Bemporad D.; Luttmann C.; Essex J. W. Computer Simulation of Small Molecule Permeation across a Lipid Bilayer: Dependence on Bilayer Properties and Solute Volume, Size, and Cross-Sectional Area. Biophys. J. 2004, 87, 1–13. PubMed PMC

Bemporad D.; Luttmann C.; Essex J. W. Behaviour of Small Solutes and Large Drugs in a Lipid Bilayer from Computer Simulations. Biochim. Biophys. Acta 2005, 1718, 1–21. PubMed

Orsi M.; Sanderson W. E.; Essex J. W. Permeability of Small Molecules through a Lipid Bilayer: a Multiscale Simulation Study. J. Phys. Chem. B 2009, 113, 12019–12029. PubMed

Boggara M. B.; Krishnamoorti R. Partitioning of Nonsteroidal Antiinflammatory Drugs in Lipid Membranes: a Molecular Dynamics Simulation Study. Biophys. J. 2010, 98, 586–595. PubMed PMC

MacCallum J. L.; Tieleman D. P. Computer Simulation of the Distribution of Hexane in a Lipid Bilayer: Spatially Resolved Free Energy, Entropy, and Enthalpy Profiles. J. Am. Chem. Soc. 2006, 128, 125–130. PubMed

Kumar S.; Rosenberg J.; Bouzida D.; Swensen R. H.; Kollman P. A. The Weighted Histogram Analysis Method for Free Energy Calculations on Biomolecules. I. The Method. J. Comput. Chem. 1992, 13, 1011–1021.

Torrie G. M.; Calleau J. P. Nonphysical Sampling Distribution in Monte Carlo Free Energy Estimation: Umbrella Sampling. J. Comput. Phys. 1997, 23, 187–199.

Marrink S. J.; Berendsen H. J. C. Permeation Process of Small Molecules across Lipid Membranes Studied by Molecular Dynamics Simulations. J. Phys. Chem. 1996, 100, 16729–16738.

Bemporad D.; Essex J. W.; Luttmann C. Permeation of Small Molecules through a Lipid Bilayer: A Computer Simulation Study. J. Phys. Chem. B 2004, 108, 4875–4884.

Zhang Y.; Voth G. A. Combined Metadynamics and Umbrella Sampling Method for the Calculation of Ion Permeation Free Energy Profiles. J. Chem. Theory Comput. 2011, 7, 2277–2283. PubMed PMC

Laio A.; Parrinello M. Escaping Free-Energy Minima. Proc. Natl. Acad. Sci. U. S. A. 2002, 99, 12562–12566. PubMed PMC

Wei C.; Pohorille A. Permeation of Membranes by Ribose and Its Diastereomers. J. Am. Chem. Soc. 2009, 131, 10237–10245. PubMed

Darve E.; Pohorille A. Calculating Free Energies Using Average Force. J. Chem. Phys. 2001, 115, 9169.

Tai K. Conformational Sampling for the Impatient. Biophys. Chem. 2004, 107, 213–220. PubMed

Roux B. The Calculation of the Potential of Mean Force Using Computer Simulations. Comput. Phys. Commun. 1995, 91, 275–282.

MacCallum J. L.; Bennett W. F. D.; Tieleman D. P. Distribution of Amino Acids in a Lipid Bilayer from Computer Simulations. Biophys. J. 2008, 94, 3393–3404. PubMed PMC

Kandt C.; Ash W. L.; Tieleman D. P. Setting up and Running Molecular Dynamics Simulations of Membrane Proteins. Methods 2007, 41, 475–488. PubMed

Humphrey W.; Dalke A.; Schulten K. VMD - Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33–38. PubMed

Tejwani R. W.; Davis M. E.; Anderson B. D.; Stouch T. R. Functional Group Dependence of Solute Partitioning to Various Locations within a DOPC Bilayer: A Comparison of Molecular Dynamics Simulations with Experiment. J. Pharm. Sci. 2011, 100, 2136–2146. PubMed

Verschueren K.Handbook of Environmental Data on Organic Chemicals, 3rd ed.; Van Nostrand Reinhold: New York,1996; pp 541–542.

Gutiérrez-Sánchez C.; Calvino-Casilda V.; Pérez-Mayoral E.; Martín-Aranda R. M.; López-Peinado A. J.; Bejblová M.; Čejka J. Coumarins Preparation by Pechmann Reaction under Ultrasound Irradiation. Synthesis of Hymecromone as Insecticide Intermediate. Cat. Lett. 2008, 128, 318–322.

Anzenbacher P.; Anzenbacherová E. Cellular and Molecular Life Sciences Cytochromes P450 and Metabolism of Xenobiotics. Cell. Mol. Life Sci. 2001, 58, 737–747. PubMed PMC

Pelkonen O.; Rautio A.; Raunio H.; Pasanen M. CYP2A6: a Human Coumarin 7-hydroxylase. Toxicology 2000, 144, 139–147. PubMed

Schüttelkopf A. W.; van Aalten D. M. F. PRODRG: a Tool for High-Throughput Crystallography of Protein-Ligand Complexes. Acta Crystallogr., Sect. D: Biol. Crystallogr. 2004, 60, 1355–1363. PubMed

Oostenbrink C.; Soares T. A.; van der Vegt N. F. A.; van Gunsteren W. F. Validation of the 53A6 GROMOS Force Field. Eur. Biophys. J. 2005, 34, 273–284. PubMed

Lemkul J. A.; Allen W. J.; Bevan D. R. Practical Considerations for Guilding GROMOS-Compatible Small-Molecule Topologies. J. Chem. Inf. Model. 2010, 50, 2221–2235. PubMed

Frisch M. J.; Trucks G. W.; Schlegel H. B.; Scuseria G. E.; Robb M. A.; Cheeseman J. R.; Montgomery J. A. Jr.; Vreven T.; Kudin K. N.; Burant J. C.; Millam J. M.; Iyengar S. S.; Tomasi J.; Barone V.; Mennucci B.; Cossi M.; Scalmani G.; Rega N.; Petersson G. A.; Nakatsuji H.; Hada M.; Ehara M.; Toyota K.; Fukuda R.; Hasegawa J.; Ishida M.; Nakajima T.; Honda Y.; Kitao O.; Nakai H.; Klene M.; Li X.; Knox J. E.; Hratchian H. P.; Cross J. B.; Bakken V.; Adamo C.; Jaramillo J.; Gomperts R.; Stratmann R. E.; Yazyev O.; Austin A. J.; Cammi R.; Pomelli C.; Ochterski J. W.; Ayala P. Y.; Morokuma K.; Voth G. A.; Salvador P.; Dannenberg J. J.; Zakrzewski V. G.; Dapprich S.; Daniels A. D.; Strain M. C.; Farkas O.; Malick D. K.; Rabuck A. D.; Raghavachari K.; Foresman J. B.; Ortiz J. V.; Cui Q.; Baboul A. G.; Clifford S.; Cioslowski J.; Stefanov B. B.; Liu G.; Liashenko A.; Piskorz P.; Komaromi I.; Martin R. L.; Fox D. J.; Keith T.; Al-Laham M. A.; Peng C. Y.; Nanayakkara A.; Challacombe M.; Gill P. M. W.; Johnson B.; Chen W.; Wong M. W.; Gonzalez C.; Pople J. A. Gaussian 03, Revision E.01; Gassuian, Inc.: Wallingford, CT, 2004.

Cieplak P.; Caldwell J.; Kollman P. Molecular Mechanical Models for Organic and Biological Systems Going Beyond the Atom Centered Two Body Additive Approximation: Aqueous Solution Free Energies of Methanol and N-Methyl Acetamide, Nucleic Acid Base, and Amide Hydrogen Bonding and Chloroform/ Water Partition Coefficients of the Nucleic Acid Bases. J. Comput. Chem. 2001, 22, 1048–1057.

Case D. A.; Darden T. A.; Cheatham T. E. III; Simmerling C. L.; Wang J.; Duke R. E.; Luo R.; Walker R. C.; Zhang W.; Merz K. M.; Roberts B.; Wang B.; Hayik S.; Roitberg A.; Seabra G.; Kolossváry I.; Wong K. F.; Paesani F.; Vanicek J.; Liu J.; Wu X.; Brozell S. R.; Steinbrecher T.; Gohlke H.; Cai Q.; Ye X.; Hsieh M.-J.; Cui G.; Roe D. R.; Mathews D. H.; Seetin M. G.; Sagui C.; Babin V.; Luchko T.; Gusarov S.; Kovalenko A.; Kollman P. A.. AMBER 11; University of California: San Francisco, 2010.

Berger O.; Edholm O.; Jahnig F. Molecular Dynamics Simulations of a Fluid Bilayer of Dipalmitoylphosphatidylcholine at Full Hydration, Constant Pressure, and Constant Temperature. Biophys. J. 1997, 72, 2002–2013. PubMed PMC

Siu S.; Vácha R.; Jungwirt P.; Böckmann R. A. Biomolecular Simulation of Membranes: Physical Properties from Different Force Fields. J. Chem. Phys. 2008, 128, 125103. PubMed

Domański J.; Stansfeld P. J.; Sansom M. S. P.; Beckstein O. Lipidbook: a Public Repository for Force-Field Parameters Used in Membrane Simulations. J. Membr. Biol. 2010, 236, 255–258. PubMed

Berendsen H. J. C.; Postma J.P.M.; Gunsteren W. F. van; Hermans J. Interaction Models for Water in Relation to Protein Hydration. In Intermol. Forces; Pullman B, Ed.; Reidel Publishing Company: 1981; pp 331–338.

Hess B.; Kutzner C.; van der Spoel D.; Lindahl E. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435–447. PubMed

Shinoda W.; Mikami M.; Baba T.; Hato M. Molecular Dynamics Study on the Effects of Chain Branching on the Physical Properties of Lipid Bilayers: 2. Permeability. J. Phys. Chem. B 2004, 108, 9346–9356.

Chiu S. W.; Clark M.; Balaji V.; Subramaniam S.; Scott H. L.; Jakobsson E. Incorporation of Surface Tension into Molecular Dynamics Simulation of an Interface: a Fluid Phase Lipid Bilayer Membrane. Biophys. J. 1995, 69, 1230–1245. PubMed PMC

Darden T.; York D.; Pedersen L. Particle Mesh Ewald: An N.log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089–10092.

Hess B.; Bekker H.; Berendsen H. J. C.; Fraaije J. G. E. M. LINCS: A Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 1463–1472.

Bussi G.; Donadio D.; Parrinello M. Canonical Sampling Through Velocity Rescaling. J. Chem. Phys. 2007, 126, 014101. PubMed

Berendsen H.; Postma J.; Vangunsteren W.; Dinola A.; Haak J. Molecular-Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 3684–3690.

Van Der Spoel D.; Lindahl E.; Hess B.; Groenhof G.; Mark A. E.; Berendsen H. J. C. GROMACS: Fast, Flexible, and Free. J. Comput. Chem. 2005, 26, 1701–1718. PubMed

Hub J. S.; Groot B. L. D.; Spoel D. V. D. g_wham-A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates. J. Chem. Theory Comput. 2010, 6, 3713–3720.

Eriksson E. S. E.; Eriksson L. A. The Influence of Cholesterol on the Properties and Permeability of Hypericin Derivatives in Lipid Membranes. J. Chem. Theory Comput. 2011, 7, 560–574. PubMed

Biocomputing at the University Of Calgary. http://people.ucalgary.ca/∼tieleman/download.html (accessed Oct. 12, 2011).

Cohen Y.; Afri M.; Frimer A. A. NMR-Based Molecular Ruler for Determining the Depth of Intercalants within the Lipid Bilayer Part II. The Preparation of a Molecular Ruler. Chem. Phys. Lipids 2008, 155, 114–119. PubMed

Cohen Y.; Bodner E.; Richman M.; Afri M.; Frimer A. A. NMR-Based Molecular Ruler for Determining the Depth of Intercalants within the Lipid Bilayer Part I. Discovering the Guidelines. Chem. Phys. Lipids 2008, 155, 98–113. PubMed

Trzesniak D.; Kunz A.-P. E.; van Gunsteren W. F. A Comparison of Methods to Compute the Potential of Mean Force. ChemPhysChem 2007, 8, 162–169. PubMed