This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.107.117770v1 94/4/1185    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Nymeyer, H. Articles by Zhou, H.-X. Search for Related Content PubMed PubMed Citation Articles by Nymeyer, H. Articles by Zhou, H.-X.

A Method to Determine Dielectric Constants in Nonhomogeneous Systems: Application to Biological Membranes

Hugh Nymeyer ^*  and Huan-Xiang Zhou  

^* Department of Chemistry and Biochemistry, Department of Physics, and School of Computational Science, Florida State University, Tallahassee, Florida 32306

Correspondence: Address reprint requests to Hugh Nymeyer, E-mail: hnymeyer{at}fsu.edu.

Continuum electrostatic models have had quantitative success in describing electrostatic-mediated phenomena on atomistic scales; however, there continues to be significant disagreement about how to assign dielectric constants in mixed, nonhomogeneous systems. We introduce a method for determining a position-dependent dielectric profile from molecular dynamics simulations. In this method, the free energy of introducing a test charge is computed two ways: from a free energy perturbation calculation and from a numerical solution to Poisson's Equation. The dielectric profile of the system is then determined by minimizing the discrepancy between these two calculations simultaneously for multiple positions of the test charge. We apply this method to determine the dielectric profile of a lipid bilayer surrounded by water. We find good agreement with dielectric models for lipid bilayers obtained by other approaches. The free energy of transferring an ion from bulk water to the lipid bilayer computed from the atomistic simulations indicates that large errors are introduced when the bilayer is represented as a single slab of low dielectric embedded in the higher-dielectric solvent. Significant improvement results from introducing an additional layer of intermediate dielectric (3) on each side of the low dielectric core extending from 12 Å to 18 Å. A small dip in transfer free energy just outside the lipid headgroups indicates the presence of a very high dielectric. These results have implications for the design of implicit membrane models and our understanding of protein-membrane interactions.