Distribution of amino acids in a lipid bilayer from computer simulations.

We have calculated the distribution in a lipid bilayer of small molecules mimicking 17 natural amino acids in atomistic detail by molecular dynamics simulation. We considered both charged and uncharged forms for Lys, Arg, Glu, and Asp. The results give detailed insight in the molecular basis of the preferred location and orientation of each side chain as well the preferred charge state for ionizable residues. Partitioning of charged and polar side chains is accompanied by water defects connecting the side chains to bulk water. These water defects dominate the energetic of partitioning, rather than simple partitioning between water and a hydrophobic phase. Lys, Glu, and Asp become uncharged well before reaching the center of the membrane, but Arg may be either charged or uncharged at the center of the membrane. Phe has a broad distribution in the membrane but Trp and Tyr localize strongly to the interfacial region. The distributions are useful for the development of coarse-grained and implicit membrane potentials for simulation and structure prediction. We discuss the relationship between the distribution in membranes, bulk partitioning to cyclohexane, and several amino acid hydrophobicity scales.

[1]  D. Tieleman,et al.  The MARTINI force field: coarse grained model for biomolecular simulations. , 2007, The journal of physical chemistry. B.

[2]  Justin L. MacCallum,et al.  Calculation of the water–cyclohexane transfer free energies of neutral amino acid side‐chain analogs using the OPLS all‐atom force field , 2003, J. Comput. Chem..

[3]  Harry A. Stern,et al.  Calculation of the dielectric permittivity profile for a nonuniform system: Application to a lipid bilayer simulation , 2003 .

[4]  Alessandra Villa,et al.  Calculation of the free energy of solvation for neutral analogs of amino acid side chains , 2002, J. Comput. Chem..

[5]  Toby W Allen,et al.  On the thermodynamic stability of a charged arginine side chain in a transmembrane helix , 2007, Proceedings of the National Academy of Sciences.

[6]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[7]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[8]  A. Mark,et al.  Coarse grained model for semiquantitative lipid simulations , 2004 .

[9]  M S Sansom,et al.  Lipid properties and the orientation of aromatic residues in OmpF, influenza M2, and alamethicin systems: molecular dynamics simulations. , 1998, Biochemistry.

[10]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[11]  C. Brooks,et al.  An implicit membrane generalized born theory for the study of structure, stability, and interactions of membrane proteins. , 2003, Biophysical journal.

[12]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[13]  Michael R. Shirts,et al.  Solvation free energies of amino acid side chain analogs for common molecular mechanics water models. , 2005, The Journal of chemical physics.

[14]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[15]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[16]  J. Killian,et al.  How proteins adapt to a membrane-water interface. , 2000, Trends in biochemical sciences.

[17]  Justin L. MacCallum,et al.  Molecular dynamics simulations of pentapeptides at interfaces: Salt bridge and cation-π interactions , 2003 .

[18]  Richard Wolfenden,et al.  Comparing the polarities of the amino acids: side-chain distribution coefficients between the vapor phase, cyclohexane, 1-octanol, and neutral aqueous solution , 1988 .

[19]  M. Cadene,et al.  X-ray structure of a voltage-dependent K+ channel , 2003, Nature.

[20]  Francisco Bezanilla,et al.  Two atomic constraints unambiguously position the S4 segment relative to S1 and S2 segments in the closed state of Shaker K channel , 2007, Proceedings of the National Academy of Sciences.

[21]  M. Klein,et al.  Partitioning of anesthetics into a lipid bilayer and their interaction with membrane-bound peptide bundles. , 2006, Biophysical journal.

[22]  Richard Wolfenden,et al.  Experimental Measures of Amino Acid Hydrophobicity and the Topology of Transmembrane and Globular Proteins , 2007, The Journal of general physiology.

[23]  D. Peter Tieleman,et al.  Modifying the OPLS‐AA force field to improve hydration free energies for several amino acid side chains using new atomic charges and an off‐plane charge model for aromatic residues , 2007, J. Comput. Chem..

[24]  J. Killian.  Gramicidin and gramicidin-lipid interactions. , 1992, Biochimica et biophysica acta.

[25]  David Baker,et al.  Voltage sensor conformations in the open and closed states in ROSETTA structural models of K(+) channels. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[26]  H. Nymeyer,et al.  Indole localization in lipid membranes revealed by molecular simulation. , 2006, Biophysical journal.

[27]  Youxing Jiang,et al.  The principle of gating charge movement in a voltage-dependent K+ channel , 2003, Nature.

[28]  Stephen H. White,et al.  Experimentally determined hydrophobicity scale for proteins at membrane interfaces , 1996, Nature Structural Biology.

[29]  E. Campbell,et al.  Crystal Structure of a Mammalian Voltage-Dependent Shaker Family K+ Channel , 2005, Science.

[30]  R. Larson,et al.  The MARTINI Coarse-Grained Force Field: Extension to Proteins. , 2008, Journal of chemical theory and computation.

[31]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[32]  Michael Feig,et al.  Molecular dynamics simulations of large integral membrane proteins with an implicit membrane model. , 2006, The journal of physical chemistry. B.

[33]  Herman J. C. Berendsen,et al.  Permeation Process of Small Molecules across Lipid Membranes Studied by Molecular Dynamics Simulations , 1996 .

[34]  S. Feller,et al.  Nuclear Overhauser enhancement spectroscopy cross-relaxation rates and ethanol distribution across membranes. , 2002, Biophysical journal.

[35]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[36]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[37]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[38]  Gunnar von Heijne,et al.  Formation of Transmembrane Helices In Vivo—Is Hydrophobicity All that Matters? , 2007 .

[39]  Justin L. MacCallum,et al.  Partitioning of Amino Acid Side Chains into Lipid Bilayers: Results from Computer Simulations and Comparison to Experiment , 2007, The Journal of general physiology.

[40]  S. Sukharev,et al.  Characterization of the resting MscS: modeling and analysis of the closed bacterial mechanosensitive channel of small conductance. , 2008, Biophysical journal.

[41]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[42]  Christian Kandt,et al.  Membrane protein simulations with a united-atom lipid and all-atom protein model: lipid–protein interactions, side chain transfer free energies and model proteins , 2006, Journal of physics. Condensed matter : an Institute of Physics journal.

[43]  G. Heijne,et al.  Recognition of transmembrane helices by the endoplasmic reticulum translocon , 2005, Nature.

[44]  D Peter Tieleman,et al.  Computer simulation of partitioning of ten pentapeptides Ace-WLXLL at the cyclohexane/water and phospholipid/water interfaces , 2005 .

[45]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[46]  Stephen H. White,et al.  Membrane Protein Insertion: The Biology–Physics Nexus , 2007, The Journal of general physiology.

[47]  M. Feig,et al.  A generalized Born formalism for heterogeneous dielectric environments: application to the implicit modeling of biological membranes. , 2005, The Journal of chemical physics.

[48]  Jonathan W. Essex,et al.  Permeation of small molecules through a lipid bilayer: a computer simulation study , 2004 .

[49]  Herman J. C. Berendsen,et al.  Simulation of Water Transport through a Lipid Membrane , 1994 .

[50]  H. Berendsen,et al.  Proton transport across transient single-file water pores in a lipid membrane studied by molecular dynamics simulations. , 1996, Biophysical journal.

[51]  Martin B Ulmschneider,et al.  Properties of integral membrane protein structures: Derivation of an implicit membrane potential , 2005, Proteins.

[52]  Justin L MacCallum,et al.  Computer simulation of the distribution of hexane in a lipid bilayer: spatially resolved free energy, entropy, and enthalpy profiles. , 2006, Journal of the American Chemical Society.

[53]  J. Essex,et al.  Computer simulation of small molecule permeation across a lipid bilayer: dependence on bilayer properties and solute volume, size, and cross-sectional area. , 2004, Biophysical journal.

[54]  Eduardo Perozo,et al.  Structure and mechanism in prokaryotic mechanosensitive channels. , 2003, Current opinion in structural biology.

[55]  Y. Jan,et al.  A quantitative assessment of models for voltage-dependent gating of ion channels. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[56]  Gunnar von Heijne,et al.  Membrane Insertion of a Potassium-Channel Voltage Sensor , 2005, Science.

[57]  T. Lazaridis.  Effective energy function for proteins in lipid membranes , 2003, Proteins.

[58]  O. Berger,et al.  Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. , 1997, Biophysical journal.

[59]  T. Creamer,et al.  Solvation energies of amino acid side chains and backbone in a family of host-guest pentapeptides. , 1996, Biochemistry.

[60]  G. von Heijne,et al.  Interface connections of a transmembrane voltage sensor. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[61]  S. White,et al.  Membrane protein folding and stability: physical principles. , 1999, Annual review of biophysics and biomolecular structure.

[62]  S. White,et al.  Location of hexane in lipid bilayers determined by neutron diffraction , 1981, Nature.

[63]  S. White,et al.  The preference of tryptophan for membrane interfaces. , 1998, Biochemistry.

[64]  Justin L MacCallum,et al.  Structures of neat and hydrated 1-octanol from computer simulations. , 2002, Journal of the American Chemical Society.

[65]  E. Campbell,et al.  Voltage Sensor of Kv1.2: Structural Basis of Electromechanical Coupling , 2005, Science.

[66]  D. Baker,et al.  Multipass membrane protein structure prediction using Rosetta , 2005, Proteins.

[67]  E. Lindahl,et al.  Amino-acid solvation structure in transmembrane helices from molecular dynamics simulations. , 2006, Biophysical journal.