A computer simulation of functional group contributions to free energy in water and a DPPC lipid bilayer.

A series of all-atom molecular dynamics simulations has been performed to evaluate the contributions of various functional groups to the free energy of solvation in water and a dipalmitoylphospatidylcholine lipid bilayer membrane and to the free energies of solute transfer (Delta(DeltaG(o))X) from water into the ordered-chain interior of the bilayer. Free energies for mutations of the alpha-H atom in p-toluic acid to six different substituents (-CH3, -Cl, -OCH3, -CN, -OH, -COOH) were calculated by a combined thermodynamic integration and perturbation method and compared to literature results from vapor pressure measurements, partition coefficients, and membrane transport experiments. Convergence of the calculated free energies was indicated by substantial declines in standard deviations for the calculated free energies with increased simulation length, by the independence of the ensemble-averaged Boltzmann factors to simulation length, and the weak dependence of hysteresis effects on simulation length over two different simulation lengths and starting from different initial configurations. Calculated values of Delta(DeltaG(o))X correlate linearly with corresponding values obtained from lipid bilayer transport experiments with a slope of 1.1 and from measurements of partition coefficients between water and hexadecane or decadiene, with slopes of 1.1 and 0.9, respectively. Van der Waals interactions between the functional group of interest and the acyl chains in the ordered chain region account for more than 95% of the overall potential energy of interaction. These results support the view that the ordered chain region within the bilayer interior is the barrier domain for transport and that solvation interactions within this region resemble those occurring in a nonpolar hydrocarbon.

[1]  Peter A. Kollman,et al.  Alternative approaches to potential of mean force calculations: Free energy perturbation versus thermodynamic integration. Case study of some representative nonpolar interactions , 1996, Journal of Computational Chemistry.

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

[3]  Terry R. Stouch,et al.  Characterization of force fields for lipid molecules: Applications to crystal structures , 1993, J. Comput. Chem..

[4]  M. Klein,et al.  Constant pressure and temperature molecular dynamics simulation of a fully hydrated liquid crystal phase dipalmitoylphosphatidylcholine bilayer. , 1995, Biophysical journal.

[5]  H. Berendsen,et al.  Molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine bilayer with different macroscopic boundary conditions and parameters , 1996 .

[6]  R. Suter,et al.  X-ray structure determination of fully hydrated L alpha phase dipalmitoylphosphatidylcholine bilayers. , 1996, Biophysical journal.

[7]  Ken A. Dill,et al.  Solute partitioning into chain molecule interphases: Monolayers, bilayer membranes, and micelles , 1986 .

[8]  J. Martín,et al.  Solution Thermodynamics from Gas-Liquid Chromatography. , 1966 .

[9]  T. Xiang,et al.  Substituent contributions to the transport of substituted p-toluic acids across lipid bilayer membranes. , 1994, Journal of pharmaceutical sciences.

[10]  L. E. Chirlian,et al.  Atomic charges derived from electrostatic potentials: A detailed study , 1987 .

[11]  Peter A. Kollman,et al.  LIGAND BINDING IN THE CATALYTIC ANTIBODY 17E8. A FREE ENERGY PERTURBATION CALCULATION STUDY , 1997 .

[12]  Keith B. Ward,et al.  Simulations of lipid crystals: Characterization of potential energy functions and parameters for lecithin molecules , 1991 .

[13]  S. Davis,et al.  Determination of thermodynamics of the methylene group in solutions of drug molecules. , 1972, The Journal of pharmacy and pharmacology.

[14]  A. Leo.  CALCULATING LOG POCT FROM STRUCTURES , 1993 .

[15]  A. Pohorille,et al.  Excess chemical potential of small solutes across water--membrane and water--hexane interfaces. , 1996, The Journal of chemical physics.

[16]  M. Mezei.  The finite difference thermodynamic integration, tested on calculating the hydration free energy difference between acetone and dimethylamine in water , 1987 .

[17]  C. Hansch,et al.  Linear free-energy relationship between partition coefficients and the aqueous solubility of organic liquids , 1968 .

[18]  J. Seelig,et al.  Deuterium order parameters in relation to thermodynamic properties of a phospholiped bilayer. A statistical mechanical interpretation. , 1975, Biochemistry.

[19]  J. Prausnitz,et al.  Activity Coefficients of Aromatic Solutes in Dilute Aqueous Solutions , 1971 .

[20]  J. Kh,et al.  Advances in Pharmaceutical Sciences , 1982 .

[21]  A. Leo,et al.  Partition coefficients and their uses , 1971 .

[22]  T. Xiang,et al.  Transport methods for probing the barrier domain of lipid bilayer membranes. , 1992, Biophysical journal.

[23]  T. Xiang,et al.  Influence of chain ordering on the selectivity of dipalmitoylphosphatidylcholine bilayer membranes for permeant size and shape. , 1998, Biophysical journal.

[24]  David J.W. Grant,et al.  Solubility behavior of organic compounds , 1990 .

[25]  Peter A. Kollman,et al.  The overlooked bond‐stretching contribution in free energy perturbation calculations , 1991 .

[26]  M. Rami Reddy,et al.  Calculation of relative solvation free energy differences by thermodynamic perturbation method: Dependence of free energy results on simulation length , 1999, J. Comput. Chem..

[27]  A. Smondyrev,et al.  Molecular dynamics simulation of DPPC bilayer in DMSO. , 1999, Biophysical journal.

[28]  T. Xiang,et al.  Influence of a Transmembrane Protein on the Permeability of Small Molecules Across Lipid Membranes , 2000, The Journal of Membrane Biology.

[29]  M. Mezei,et al.  Calculation of the Free Energy Profile of H2O, O2, CO, CO2, NO, and CHCl3 in a Lipid Bilayer with a Cavity Insertion Variant of the Widom Method , 2000 .

[30]  Arieh Ben-Naim,et al.  Solvation thermodynamics of nonionic solutes , 1984 .

[31]  P. Kollman,et al.  Biomolecular simulations: recent developments in force fields, simulations of enzyme catalysis, protein-ligand, protein-protein, and protein-nucleic acid noncovalent interactions. , 2001, Annual review of biophysics and biomolecular structure.

[32]  M. Erion,et al.  Calculation of Relative Hydration Free Energy Differences for Heteroaromatic Compounds: Use in the Design of Adenosine Deaminase and Cytidine Deaminase Inhibitors , 1998 .

[33]  R. Pastor,et al.  Computer simulation of liquid/liquid interfaces. II. Surface tension-area dependence of a bilayer and monolayer , 1995 .

[34]  E. D. Cyan.  Handbook of Chemistry and Physics , 1970 .

[35]  Ronald M. Levy,et al.  SOLVATION FREE ENERGIES OF SMALL AMIDES AND AMINES FROM MOLECULAR DYNAMICS/FREE ENERGY PERTURBATION SIMULATIONS USING PAIRWISE ADDITIVE AND MANY-BODY POLARIZABLE POTENTIALS , 1995 .

[36]  K. Miller,et al.  The solubility of hydrocarbon gases in lipid bilayers , 1977 .

[37]  J. Seelig,et al.  The dynamic structure of fatty acyl chains in a phospholipid bilayer measured by deuterium magnetic resonance. , 1974, Biochemistry.

[38]  Terry R. Stouch,et al.  Effects of Switching Functions on the Behavior of Liquid Water in Molecular Dynamics Simulations , 1994 .

[39]  M. Rami Reddy,et al.  A molecular dynamics study of the structure and dynamics of water between dilauroylphosphatidylethanolamine bilayers , 1992 .

[40]  J. Butler.  The energy and entropy of hydration of organic compounds , 1937 .

[41]  E. Wright,et al.  Biological membranes: the physical basis of ion and nonelectrolyte selectivity. , 1969, Annual review of physiology.

[42]  J. Hine,et al.  Structural effects on rates and equilibriums. XIX. Intrinsic hydrophilic character of organic compounds. Correlations in terms of structural contributions , 1975 .

[43]  William H. Press,et al.  Numerical recipes , 1990 .

[44]  P A Kollman,et al.  Ranking ligand binding affinities with avidin: a molecular dynamics‐based interaction energy study , 1999, Proteins.

[45]  T. Xiang,et al.  Molecular distributions in interphases: statistical mechanical theory combined with molecular dynamics simulation of a model lipid bilayer. , 1994, Biophysical journal.

[46]  T. Higuchi,et al.  Thermodynamic group contributions for hydroxyl, amino, and methylene groups. , 1978, Journal of pharmaceutical sciences.

[47]  J. Hine,et al.  THE INTRINSIC HYDROPHILIC CHARACTER OF ORGANIC COMPOUNDS, CORRELATIONS IN TERMS OF STRUCTURAL CONTRIBUTIONS , 1975 .

[48]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[49]  R. C. Weast.  CRC Handbook of Chemistry and Physics , 1973 .

[50]  Interaction of K+ with a Phospholipid Bilayer: A Molecular Dynamics Study , 1997 .

[51]  T. Xiang,et al.  The Barrier Domain for Solute Permeation Varies With Lipid Bilayer Phase Structure , 1998, The Journal of Membrane Biology.

[52]  T. Xiang,et al.  Translational Diffusion in Lipid Bilayers: Dynamic Free-Volume Theory and Molecular Dynamics Simulation , 1999 .

[53]  K. James.  Drug design , 1992, Nature.

[54]  R. Prosser,et al.  Deuterium NMR relaxation studies of peptide-lipid interactions. , 1992, Biochemistry.

[55]  Towards the development of a second-order approximation in activity coefficient models based on group contributions , 1996 .

[56]  T. Xiang,et al.  MEAN MOLECULAR POTENTIALS IN A MODEL LIPID BILAYER: A MOLECULAR DYNAMICS SIMULATION , 1995 .

[57]  C. Tanford,et al.  The solubility of amino acids and two glycine peptides in aqueous ethanol and dioxane solutions. Establishment of a hydrophobicity scale. , 1971, The Journal of biological chemistry.

[58]  Volkhard Helms,et al.  Free energies of hydration from thermodynamic integration: Comparison of molecular mechanics force fields and evaluation of calculation accuracy , 1997, J. Comput. Chem..

[59]  J. Seelig,et al.  Two pictures of a lipid bilayer. A comparison between deuterium label and spin-label experiments. , 1974, Biochemistry.

[60]  C. Hansch,et al.  SUBSTITUENT CONSTANTS FOR ALIPHATIC FUNCTIONS OBTAINED FROM PARTITION COEFFICIENTS. , 1965, Journal of medicinal chemistry.

[61]  S H White,et al.  Hydrophobic interactions of peptides with membrane interfaces. , 1998, Biochimica et biophysica acta.

[62]  J. Andrew McCammon,et al.  Free energy difference calculations by thermodynamic integration: Difficulties in obtaining a precise value , 1991 .

[63]  R. Collander,et al.  The Permeability of Nitella Cells to Non‐Eleetrolytes , 1954 .

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

[65]  W. Richards,et al.  Calculation of the tautomer ratio of histamine in aqueous solution using free energy perturbation methods : an in-depth study , 1994 .

[66]  William I. Higuchi,et al.  Heterogeneity Effects on Permeability–Partition Coefficient Relationships in Human Stratum Corneum , 1988, Pharmaceutical Research.

[67]  T. Xiang,et al.  Molecular dissolution processes in lipid bilayers: A molecular dynamics simulation , 1999 .

[68]  Heather A. Carlson,et al.  Accuracy of free energies of hydration for organic molecules from 6‐31g*‐derived partial charges , 1993, J. Comput. Chem..

[69]  A Kitao,et al.  Improved protein free energy calculation by more accurate treatment of nonbonded energy: Application to chymotrypsin inhibitor 2, V57A , 1998, Proteins.

[70]  Kazuhisa Ohtaguchi,et al.  Application of Principal Components Analysis to Calculate Henry's Constant from Molecular Structure , 1992, Comput. Chem..

[71]  Terry R. Stouch,et al.  Computer simulation of a phospholipid monolayer‐water system: The influence of long range forces on water structure and dynamics , 1993 .

[72]  M Mezei,et al.  Free Energy Simulations a , 1986, Annals of the New York Academy of Sciences.

[73]  P. Kollman,et al.  Molecular Dynamics Free Energy Perturbation Calculations: Influence of Nonbonded Parameters on the Free Energy of Hydration of Charged and Neutral Species , 1994 .

[74]  E. Jakobsson,et al.  Incorporation of surface tension into molecular dynamics simulation of an interface: a fluid phase lipid bilayer membrane. , 1995, Biophysical journal.

[75]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

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

[77]  T R Stouch,et al.  Solute diffusion in lipid bilayer membranes: an atomic level study by molecular dynamics simulation. , 1993, Biochemistry.