A molecular dynamics study of the pores formed by Escherichia coli OmpF porin in a fully hydrated palmitoyloleoylphosphatidylcholine bilayer

In this paper we study the properties of pores formed by OmpF porin from Escherichia coli, based on a molecular dynamics simulation of the OmpF trimer, 318 palmitoyl-oleoyl-phosphatidylethanolamine lipids, 27 Na+ ions, and 12,992 water molecules. After equilibration and a nanosecond production run, the OmpF trimer exhibits a C-alpha root mean square deviation from the crystal structure of 0.23 nm and a stable secondary structure. No evidence is found for large-scale motions of the L3 loop. We investigate the pore dimensions, conductance, and the properties of water inside the pore. This water forms a complicated pattern, even when averaged over 1 ns of simulation time. Around the pore constriction zone the water dipoles are highly structured in the plane of the membrane, oriented by the strong transversal electric field. In addition, there is a net orientation along the pore axis pointing from the extracellular to the intracellular side of the bilayer. The diffusion coefficients of water inside the pore are greatly reduced compared to bulk. We compare our results to results from model pores (Breed et at., 1996. Biophys. J. 70:1643-1661; Sansom et al. 1997. Biophys. J. 73:2404-2415) and discuss implications for further theoretical work.