Molecular dynamics simulations of hydrophilic pores in lipid bilayers

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.     

An amphipathic peptide from the C-terminal region of the human immunodeficiency virus envelope glycoprotein causes pore formation in membranes.

The peptide fragment of the carboxy-terminal region of the human immunodeficiency virus (HIV) transmembrane protein (gp41) has been implicated in T-cell death. This positively charged, amphipathic helix (amino acids 828 to 848) of the envelope protein is located within virions or cytoplasm. We studied the interaction of the isolated, synthetic amphipathic helix of gp41 with planar phospholipid bilayer membranes and with Sf9 cells using voltage clamp, potentiodynamic, and single-cell recording techniques. We found that the peptide binds strongly to planar membranes, especially to the negatively charged phosphatidylserine bilayer. In the presence of micromolar concentrations of peptide sufficient to make its surface densities comparable with those of envelope glycoprotein molecules in HIV virions, an increase in bilayer conductance and a decrease in bilayer stability were observed, showing pore formation in the planar lipid bilayers. These pores were permeable to both monovalent and divalent cations, as well as to chloride. The exposure of the inner leaflet of cell membranes to even 25 nM peptide increased membrane conductance. We suggest that the carboxy-terminal fragment of the HIV type 1 envelope protein may interact with the cell membrane of infected T cells to create lipidic pores which increase membrane permeability, leading to sodium and calcium flux into cells, osmotic swelling, and T-cell necrosis or apoptosis.     

Barrel-stave model or toroidal model? A case study on melittin pores

Transmembrane pores induced by amphiphilic peptides, including melittin, are often modeled with the barrel-stave model after the alamethicin pore. We examine this assumption on melittin by using two methods, oriented circular dichroism (OCD) for detecting the orientation of melittin helix and neutron scattering for detecting transmembrane pores. OCD spectra of melittin were systematically measured. Melittin can orient either perpendicularly or parallel to a lipid bilayer, depending on the physical condition and the composition of the bilayer. Transmembrane pores were detected when the helices oriented perpendicularly to the plane of the bilayers, not when the helices oriented parallel to the bilayers. The evidence that led to the barrel-stave model for alamethicin and that to the toroidal model for magainin were reviewed. The properties of melittin pores are closely similar to that of magainin but unlike that of alamethicin. We conclude that, among naturally produced peptides that we have investigated, only alamethicin conforms to the barrel-stave model. Other peptides, including magainins, melittin and protegrins, all appear to induce transmembrane pores that conform to the toroidal model in which the lipid monolayer bends continuously through the pore so that the water core is lined by both the peptides and the lipid headgroups.     

Restricted diffusion of molecules in porous affinity chromatography adsorbents.

A restricted diffusion model is constructed and solved in order to study the permeability of large adsorbate molecules in the pores of affinity chromatography media, when the adsorbate molecules are adsorbed onto immobilized ligands. The combined effects of steric hindrance at the entrance to the pores and frictional resistance within the pores, as well as the effects of pore size distribution, pore connectivity of the adsorbent, molecular size of adsorbate and ligand, and the fractional saturation of adsorption sites (ligands), are considered. Affinity adsorbents with dilute and high ligand concentrations are examined, and the permeability of the adsorbate in porous networks of connectivity nT is studied by means of effective medium approximation (EMA) numerical solutions. As expected, the permeability of the adsorbate decreases as the size of the adsorbate and/or ligand molecule increases. The permeability also decreases when the fractional saturation of the ligands increases, as well as when the pore connectivity of the network decreases. The dependence of the permeability on the pore connectivity tends to be less marked in adsorbents with concentrated ligand than in porous media with dilute ligand concentration. The conditions are also presented for which the percolation threshold is attained in a number of different systems. The restricted diffusion model and results of this work may be of importance in studies involving the modeling, prediction of the dynamic behavior, design, and control of affinity chromatography (biospecific adsorption) systems employing porous adsorbents. The theoretical results may also have important implications in the selection of a ligand as well as in the selection and construction of an affinity porous matrix, so that the adsorbate of interest can be efficiently separated from a given solution. Furthermore, with appropriate modifications this restricted diffusion model may be used in studies involving the immobilization of ligands or enzymes in porous solids.     

The effect of hexadecaprenol on molecular organisation and transport properties of model membranes

The Langmuir monolayer technique and voltammetric analysis were used to investigate the properties of model lipid membranes prepared from dioleoylphosphatidylcholine (DOPC), hexadecaprenol (C80), and their mixtures. Surface pressure-molecular area isotherms, current-voltage characteristics, and membrane conductance-temperature were measured. Molecular area isobars, specific molecular areas, excess free energy of mixing, collapse pressure and collapse area were determined for lipid monolayers. Membrane conductance, activation energy of ion migration across the membrane, and membrane permeability coefficient for chloride ions were determined for lipid bilayers. Hexadecaprenol decreases the activation energy and increases membrane conductance and membrane permeability coefficient. The results of monolayer and bilayer investigations show that some electrical, transport and packing properties of lipid membranes change under the influence of hexadecaprenol. The results indicate that hexadecaprenol modulates the molecular organisation of the membrane and that the specific molecular area of polyprenol molecules depends on the relative concentration of polyprenols in membranes. We suggest that hexadecaprenol modifies lipid membranes by the formation of fluid microdomains. The results also indicate that electrical transmembrane potential can accelerate the formation of pores in lipid bilayers modified by long chain polyprenols.     

Synthetic multifunctional pores that open and close in response to chemical stimulation

Studies on synthetic multifunctional pores with external and internal active sites for ligand gating and noncompetitive blockage are presented, with emphasis on the contribution of external ligands to the characteristics of pore. A comparison between different synthetic multifunctional pores reveals that the location of functional groups in rigid-rod beta-barrel pores is precisely reflected in the function: molecular recognition at the outer barrel surface results in pore opening, while molecular recognition at the inner barrel surface results in pore closing. Negligible nonspecific leakage, disappearance of pH gating, inhibition of intervesicular pore transfer, and maybe also the flickering of currents of single open pores characterize external ligands as adhesive cushions that liberate the pore from lateral pressure exerted by the surrounding membrane. Refined molecular models show good agreement with pore design and experimental facts with regard to function.     

Organizations of melittin peptides after spontaneous penetration into cell membranes

The antimicrobial peptide, melittin, is a potential next-generation antibiotic because melittin can spontaneously form pores in bacterial cell membranes and cause cytoplasm leakage. However, the organizations of melittin peptides in cell membranes remain elusive, which impedes the understanding of the poration mechanism. In this work, we use coarse-grained and all-atom molecular dynamics (MD) simulations to investigate the organizations of melittin peptides during and after spontaneous penetration into DPPC/POPG lipid bilayers. We find that the peptides in lipid bilayers adopt either a transmembrane conformation or a U-shaped conformation, which are referred to as T- and U-peptides, respectively. Several U-peptides and/or T-peptides aggregate to form stable pores. We analyze a T-pore consisting of four T-peptides and a U-pore consisting of three U-peptides and one T-peptide. In both pores, peptides are organized in a manner such that polar residues face inward and hydrophobic residues face outward, which stabilizes the pores and produces water channels. Compared with the U-pore, the T-pore has lower energy, larger pore diameter, and higher permeability. However, the T-pore occurs less frequently than the U-pore in our simulations, probably because the formation of the T-pore is kinetically slower than the U-pore. The stability and permeability of both pores are confirmed by 300 ns all-atom MD simulations. The peptide organizations obtained in this work should deepen the understanding of the stability, poration mechanism, and permeability of melittin, and facilitate the optimization of melittin to enhance the antibacterial ability.     

Role of cholesterol flip-flop in oxidized lipid bilayers

We performed a series of molecular dynamics simulations of cholesterol (Chol) in nonoxidized 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) bilayer and in binary mixtures of PLPC-oxidized-lipid-bilayers with 0–50% Chol concentration and oxidized lipids with hydroperoxide and aldehyde oxidized functional groups. From the 60 unbiased molecular dynamics simulations (total of 161 μs), we found that Chol inhibited pore formation in the aldehyde-containing oxidized lipid bilayers at concentrations greater than 11%. For both pure PLPC bilayer and bilayers with hydroperoxide lipids, no pores were observed at any Chol concentration. Furthermore, increasing cholesterol concentration led to a change of phase state from the liquid-disordered to the liquid-ordered phase. This condensing effect of Chol was observed in all systems. Data analysis shows that the addition of Chol results in an increase in bilayer thickness. Interestingly, we observed Chol flip-flop only in the aldehyde-containing lipid bilayer but neither in the PLPC nor the hydroperoxide bilayers. Umbrella-sampling simulations were performed to calculate the translocation free energies and the Chol flip-flop rates. The results show that Chol’s flip-flop rate depends on the lipid bilayer type, and the highest rate are found in aldehyde bilayers. As the main finding, we shown that Chol stabilizes the oxidized lipid bilayer by confining the distribution of the oxidized functional groups.     

Permeation of halide anions through phospholipid bilayers occurs by the solubility-diffusion mechanism.

Two alternative mechanisms are frequently used to describe ionic permeation of lipid bilayers. In the first, ions partition into the hydrophobic phase and then diffuse across (the solubility-diffusion mechanism). The second mechanism assumes that ions traverse the bilayer through transient hydrophilic defects caused by thermal fluctuations (the pore mechanism). The theoretical predictions made by both models were tested for halide anions by measuring the permeability coefficients for chloride, bromide, and iodide as a function of bilayer thickness, ionic radius, and sign of charge. To vary the bilayer thickness systematically, liposomes were prepared from monounsaturated phosphatidylcholines (PC) with chain lengths between 16 and 24 carbon atoms. The fluorescent dye MQAE (N-(ethoxycarbonylmethyl)-6-methoxyquinolinium bromide) served as an indicator for halide concentration inside the liposomes and was used to follow the kinetics of halide flux across the bilayer membranes. The observed permeability coefficients ranged from 10(-9) to 10(-7) cm/s and increased as the bilayer thickness was reduced. Bromide was found to permeate approximately six times faster than chloride through bilayers of identical thickness, and iodide permeated three to four times faster than bromide. The dependence of the halide permeability coefficients on bilayer thickness and on ionic size were consistent with permeation of hydrated ions by a solubility-diffusion mechanism rather than through transient pores. Halide permeation therefore differs from that of a monovalent cation such as potassium, which has been accounted for by a combination of the two mechanisms depending on bilayer thickness.     

Interactions of bovine lactoferricin with acidic phospholipid bilayers and its antimicrobial activity as studied by solid-state NMR

Bovine lactoferricin (LfcinB) is an antimicrobial peptide released by pepsin cleavage of lactoferrin. In this work, the interaction between LfcinB and acidic phospholipid bilayers with the weight percentage of 65% dimyristoylphosphatidylglycerol (DMPG), 10% cardiolipin (CL) and 25% dimyristoylphosphatidylcholine (DMPC) was investigated as a mimic of cell membrane of Staphylococcus aureus by means of quartz crystal microbalance (QCM) and solid-state (31)P and (1)H NMR spectroscopy. Moreover, we elucidated a molecular mechanism of the antimicrobial activity of LfcinB by means of potassium ion selective electrode (ISE). It turned out that affinity of LfcinB for acidic phospholipid bilayers was higher than that for neutral phospholipid bilayers. It was also revealed that the association constant of LfcinB was larger than that of lactoferrin as a result of QCM measurements. (31)P DD-static NMR spectra indicated that LfcinB interacted with acidic phospholipid bilayers and bilayer defects were observed in the bilayer systems because isotropic peaks were clearly appeared. Gel-to-liquid crystalline phase transition temperatures (Tc) in the mixed bilayer systems were determined by measuring the temperature variation of relative intensities of acyl chains in (1)H MAS NMR spectra. Tc values of the acidic phospholipid and LfcinB-acidic phospholipid bilayer systems were 21.5 degrees C and 24.0 degrees C, respectively. To characterize the bilayer defects, potassium ion permeation across the membrane was observed by ISE measurements. The experimental results suggest that LfcinB caused pores in the acidic phospholipid bilayers. Because these pores lead the permeability across the membrane, the molecular mechanism of the antimicrobial activity could be attributed to the pore formation in the bacterial membrane induced by LfcinB.     

Energetics of pore formation induced by membrane active peptides

Antimicrobial peptides are known to form pores in cell membranes. We study this process in model bilayers of various lipid compositions. We use two of the best-studied peptides, alamethicin and melittin, to represent peptides making two types of pores, that is, barrel-stave pores and toroidal pores. In both cases, the key control variable is the concentration of the bound peptides in the lipid bilayers (expressed in the peptide-lipid molar ratio, P/L). The method of oriented circular dichroism (OCD) was used to monitor the peptide orientation in bilayers as a function of P/L. The same samples were scanned by X-ray diffraction to measure the bilayer thickness. In all cases, the bilayer thickness decreases linearly with P/L and then levels off after P/L exceeds a lipid-dependent critical value, (P/L)*. OCD spectra showed that the helical peptides are oriented parallel to the bilayers as long as P/L < (P/L)*, but as P/L increases over (P/L)*, an increasing fraction of peptides changed orientation to become perpendicular to the bilayer. We analyzed the data by assuming an internal membrane tension associated with the membrane thinning. The free energy containing this tension term leads to a relation explaining the P/L-dependence observed in the OCD and X-ray diffraction measurements. We extracted the experimental parameters from this thermodynamic relation. We believe that they are the quantities that characterize the peptide-lipid interactions related to the mechanism of pore formation. We discuss the meaning of these parameters and compare their values for different lipids and for the two different types of pores. These experimental parameters are useful for further molecular analysis and are excellent targets for molecular dynamic simulation studies.     

The internal cavity of the staphylococcal alpha-hemolysin pore accommodates approximately 175 exogenous amino acid residues

The cavity within the cap domain of the transmembrane staphylococcal alpha-hemolysin (alphaHL) pore is roughly a sphere of diameter approximately 45 A (molecular surface volume approximately 39,500 A(3)). We tested the ability of the cavity to accommodate exogenous polypeptide chains. Concatemerized Gly/Ser-containing sequences ("loops", L; number of repeats = n; number of residues = 10n + 5, n = 0-21) were inserted at a position located within the cavity of the fully assembled heptameric alphaHL pore. Homomeric pores containing 25 or less residues in each loop (n or= 7, only one L subunit was incorporated. As the inserted loop was lengthened, transient closures were observed in planar bilayer experiments with single pores. However, L(1)W(6) pores with very long loops (n = 14 and 21) had unitary conductance values close to those of W(7), suggesting that the loop is extruded through the opening in the cap of the pore into the external medium. Further analysis of bilayer recordings and electrophoretic migration patterns indicates that the upper capacity of the cavity is approximately 175 amino acids. The findings suggest that small functional peptides or proteins might be assembled within the alphaHL pore.     

Life Cycle of an Electropore: Field-Dependent and Field-Independent Steps in Pore Creation and Annihilation

Electropermeabilization, an electric field-induced modification of the barrier functions of the cell membrane, is widely used in laboratories and increasingly in the clinic; but the mechanisms and physical structures associated with the electromanipulation of membrane permeability have not been definitively characterized. Indirect experimental observations of electrical conductance and small molecule transport as well as molecular dynamics simulations have led to models in which hydrophilic pores form in phospholipid bilayers with increased probability in the presence of an electric field. Presently available methods do not permit the direct, nanoscale examination of electroporated membranes that would confirm the existence of these structures. To facilitate the reconciliation of poration models with the observed properties of electropermeabilized lipid bilayers and cell membranes, we propose a scheme for characterizing the stages of electropore formation and resealing. This electropore life cycle, based on molecular dynamics simulations of phospholipid bilayers, defines a sequence of discrete steps in the electric field-driven restructuring of the membrane that leads to the formation of a head group-lined, aqueous pore and then, after the field is removed, to the dismantling of the pore and reassembly of the intact bilayer. Utilizing this scheme we can systematically analyze the interactions between the electric field and the bilayer components involved in pore initiation, construction and resealing. We find that the pore creation time depends strongly on the electric field gradient across the membrane interface and that the pore annihilation time is at least weakly dependent on the magnitude of the pore-initiating electric field and, in general, much longer than the pore creation time.     

Exploring gas permeability of cellular membranes and membrane channels with molecular dynamics

Aquaporins are a family of membrane proteins specialized in rapid water conduction across biological membranes. Whether these channels also conduct gas molecules and the physiological significance of this potential function have not been well understood. Here we report 140 ns of molecular dynamics simulations of membrane-embedded AQP1 and of a pure POPE bilayer addressing these questions. The permeability of AQP1 to two types of gas molecules, O2 and CO2, was investigated using two complementary methods, namely, explicit gas diffusion simulation and implicit ligand sampling. The simulations show that the central (tetrameric) pore of AQP1 can be readily used by either gas molecule to permeate the channel. The two approaches produced similar free energy profiles associated with gas permeation through the central pore: a -0.4 to -1.7 kcal/mol energy well in the middle, and a 3.6-4.6 kcal/mol energy barrier in the periplasmic vestibule. The barrier appears to be mainly due to a dense cluster of water molecules anchored in the periplasmic mouth of the central pore by four aspartate residues. Water pores show a very low permeability to O2, but may contribute to the overall permeation of CO2 due to its more hydrophilic nature. Although the central pore of AQP1 is found to be gas permeable, the pure POPE bilayer provides a much larger cross-sectional area, thus exhibiting a much lower free energy barrier for CO2 and O2 permeation. As such, gas conduction through AQP1 may only be physiologically relevant either in membranes of low gas permeability, or in cells where a major fraction of the cellular membrane is occupied by AQPs.     

Microfabricated teflon membranes for low-noise recordings of ion channels in planar lipid bilayers

We present a straightforward, accessible method for the fabrication of micropores with diameters from 2 to 800 micro m in films of amorphous Teflon (Teflon AF). Pores with diameters 相似文献   

Comparative study of molecular dynamics,diffusion, and permeability for ligands in biomembranes of different lipid composition

A comparative study of several model lipid bilayers of different composition, which included analysis of kinetic parameters of model lipid bilayers and permeability of bilayer membranes for small molecules, has been carried out. The conformity of results of numeric experiments to experimental data (structure of membrane lipid bilayers, lateral diffusion coefficients, and relative permeability of biomembranes for ligands) is discussed in the framework of a standard molecular dynamics protocol.     

Permeation of protons, potassium ions, and small polar molecules through phospholipid bilayers as a function of membrane thickness.

Two mechanisms have been proposed to account for solute permeation of lipid bilayers. Partitioning into the hydrophobic phase of the bilayer, followed by diffusion, is accepted by many for the permeation of water and other small neutral solutes, but transient pores have also been proposed to account for both water and ionic solute permeation. These two mechanisms make distinctively different predictions about the permeability coefficient as a function of bilayer thickness. Whereas the solubility-diffusion mechanism predicts only a modest variation related to bilayer thickness, the pore model predicts an exponential relationship. To test these models, we measured the permeability of phospholipid bilayers to protons, potassium ions, water, urea, and glycerol. Bilayers were prepared as liposomes, and thickness was varied systematically by using unsaturated lipids with chain lengths ranging from 14 to 24 carbon atoms. The permeability coefficient of water and neutral polar solutes displayed a modest dependence on bilayer thickness, with an approximately linear fivefold decrease as the carbon number varied from 14 to 24 atoms. In contrast, the permeability to protons and potassium ions decreased sharply by two orders of magnitude between 14 and 18 carbon atoms, and leveled off, when the chain length was further extended to 24 carbon atoms. The results for water and the neutral permeating solutes are best explained by the solubility-diffusion mechanism. The results for protons and potassium ions in shorter-chain lipids are consistent with the transient pore model, but better fit the theoretical line predicted by the solubility-diffusion model at longer chain lengths.     

Interactions of bovine lactoferricin with acidic phospholipid bilayers and its antimicrobial activity as studied by solid-state NMR

Bovine lactoferricin (LfcinB) is an antimicrobial peptide released by pepsin cleavage of lactoferrin. In this work, the interaction between LfcinB and acidic phospholipid bilayers with the weight percentage of 65% dimyristoylphosphatidylglycerol (DMPG), 10% cardiolipin (CL) and 25% dimyristoylphosphatidylcholine (DMPC) was investigated as a mimic of cell membrane of Staphylococcus aureus by means of quartz crystal microbalance (QCM) and solid-state ³¹P and ¹H NMR spectroscopy. Moreover, we elucidated a molecular mechanism of the antimicrobial activity of LfcinB by means of potassium ion selective electrode (ISE). It turned out that affinity of LfcinB for acidic phospholipid bilayers was higher than that for neutral phospholipid bilayers. It was also revealed that the association constant of LfcinB was larger than that of lactoferrin as a result of QCM measurements. ³¹P DD-static NMR spectra indicated that LfcinB interacted with acidic phospholipid bilayers and bilayer defects were observed in the bilayer systems because isotropic peaks were clearly appeared. Gel-to-liquid crystalline phase transition temperatures (Tc) in the mixed bilayer systems were determined by measuring the temperature variation of relative intensities of acyl chains in ¹H MAS NMR spectra. Tc values of the acidic phospholipid and LfcinB-acidic phospholipid bilayer systems were 21.5 °C and 24.0 °C, respectively. To characterize the bilayer defects, potassium ion permeation across the membrane was observed by ISE measurements. The experimental results suggest that LfcinB caused pores in the acidic phospholipid bilayers. Because these pores lead the permeability across the membrane, the molecular mechanism of the antimicrobial activity could be attributed to the pore formation in the bacterial membrane induced by LfcinB.     

Melittin-induced bilayer leakage depends on lipid material properties: evidence for toroidal pores

The membrane-lytic peptide melittin has previously been shown to form pores in lipid bilayers that have been described in terms of two different structural models. In the "barrel stave" model the bilayer remains more or less flat, with the peptides penetrating across the bilayer hydrocarbon region and aggregating to form a pore, whereas in the "toroidal pore" melittin induces defects in the bilayer such that the bilayer bends sharply inward to form a pore lined by both peptides and lipid headgroups. Here we test these models by measuring both the free energy of melittin transfer (DeltaG degrees ) and melittin-induced leakage as a function of bilayer elastic (material) properties that determine the energetics of bilayer bending, including the area compressibility modulus (K(a)), bilayer bending modulus (k(c)), and monolayer spontaneous curvature (R(o)). The addition of cholesterol to phosphatidylcholine (PC) bilayers, which increases K(a) and k(c), decreases both DeltaG degrees and the melittin-induced vesicle leakage. In contrast, the addition to PC bilayers of molecules with either positive R(o), such as lysoPC, or negative R(o), such as dioleoylglycerol, has little effect on DeltaG degrees , but produces large changes in melittin-induced leakage, from 86% for 8:2 PC/lysoPC to 18% for 8:2 PC/dioleoylglycerol. We observe linear relationships between melittin-induced leakage and both K(a) and 1/R(o)(2). However, in contrast to what would be expected for a barrel stave model, there is no correlation between observed leakage and bilayer hydrocarbon thickness. All of these results demonstrate the importance of bilayer material properties on melittin-induced leakage and indicate that the melittin-induced pores are defects in the bilayer lined in part by lipid molecules.     

Ion transport across transmembrane pores

To study the pore-mediated transport of ionic species across a lipid membrane, a series of molecular dynamics simulations have been performed of a dipalmitoyl-phosphatidyl-choline bilayer containing a preformed water pore in the presence of sodium and chloride ions. It is found that the stability of the transient water pores is greatly reduced in the presence of the ions. Specifically, the binding of sodium cations at the lipid/water interface increases the pore line tension, resulting in a destabilization of the pore. However, the application of mechanical stress opposes this effect. The flux of ions through these mechanically stabilized pores has been analyzed. Simulations indicate that the transport of the ions through the pores depends strongly on the size of the water channel. In the presence of small pores (radius <1.5 nm) permeation is slow, with both sodium and chloride permeating at similar rates. In the case in which the pores are larger (radius >1.5 nm), a crossover is observed to a regime where the anion flux is greatly enhanced. Based on these observations, a mechanism for the basal membrane permeability of ions is discussed.