Role of charged residues at the OmpF porin channel constriction probed by mutagenesis and simulation

The channel constriction of OmpF porin, a pore protein in the bacterial outer membrane, is highly charged due to the presence of three arginines (R42, R82, and R132) and two acidic residues (D113 and E117). The influence of these charges on ion conductance, ion selectivity, and voltage gating has been studied with mutants D113N/E117Q, R42A/R82A/R132A/D113N/E117Q, and V18K/G131K, which were designed to remove or add protein charge at the channel constriction. The crystal structures revealed no or only local changes compared to wild-type OmpF, thus allowing a comparative study. The single-channel conductance of the isosteric D113N/E117Q variant was found to be 2-fold reduced, and that of the pentuple mutant was 70% of the wild-type value, despite a considerably larger pore cross section. Ion selectivity was drastically altered by the mutations with cation/anion permeability ratios ranging from 1 to 12. Ion flow through these and eight other mutants, which have been characterized previously, was simulated by Brownian dynamics based on the detailed crystal structures. The calculated ion selectivity and relative channel conductance values agree well with the experimental data. This demonstrates that ion translocation through porin is mainly governed by pore geometry and charge, the two factors that are properly represented in the simulations.     

Ion permeation and selectivity of OmpF porin: a theoretical study based on molecular dynamics,Brownian dynamics,and continuum electrodiffusion theory

Three different theoretical approaches are used and compared to refine our understanding of ion permeation through the channel formed by OmpF porin from Escherichia coli. Those approaches are all-atom molecular dynamics (MD) in which ions, solvent, and lipids are represented explicitly, Brownian dynamics (BD) in which ions are represented explicitly, while solvent and lipids are represented as featureless dielectrics, and Poisson-Nernst-Planck (PNP) electrodiffusion theory in which both solvent and local ion concentrations are represented as a continuum. First, the ability of the different theoretical approaches in reproducing the equilibrium average ion density distribution in OmpF porin bathed by a 1M KCl symmetric salt solution is examined. Under those conditions the PNP theory is equivalent to the non-linear Poisson-Boltzmann (PB) theory. Analysis shows that all the three approaches are able to capture the important electrostatic interactions between ions and the charge distribution of the channel that govern ion permeation and selectivity in OmpF. The K(+) and Cl(-) density distributions obtained from the three approaches are very consistent with one another, which suggests that a treatment on the basis of a rigid protein and continuum dielectric solvent is valid in the case of OmpF. Interestingly, both BD and continuum electrostatics reproduce the distinct left-handed twisted ion pathways for K(+) and Cl(-) extending over the length of the pore which were observed previously in MD. Equilibrium BD simulations in the grand canonical ensemble indicate that the channel is very attractive for cations, particularly at low salt concentration. On an average there is 1.55 K(+) inside the pore in 10mM KCl. Remarkably, there is still 0.17 K(+) on average inside the pore even at a concentration as low as 1microM KCl. Secondly, non-equilibrium ion flow through OmpF is calculated using BD and PNP and compared with experimental data. The channel conductance in 0.2M and 1M KCl calculated using BD is in excellent accord with the experimental data. The calculations reproduce the experimentally well-known conductance-concentration relation and also reveal an asymmetry in the channel conductance (a larger conductance is observed under a positive transmembrane potential). Calculations of the channel conductance for three mutants (R168A, R132A, and K16A) in 1M KCl suggest that the asymmetry in the channel conductance arises mostly from the permanent charge distribution of the channel rather than the shape of the pore itself. Lastly, the calculated reversal potential in a tenfold salt gradient (0.1:1M KCl) is 27.4(+/-1.3)mV (BD) and 22.1(+/-0.6)mV (PNP), in excellent accord with the experimental value of 24.3mV. Although most of the results from PNP are qualitatively reasonable, the calculated channel conductance is about 50% higher than that calculated from BD probably because of a lack of some dynamical ion-ion correlations.     

Brownian dynamics simulation of ion flow through porin channels

Bacterial porins, which allow the passage of solutes across the outer bacterial membrane, are structurally well characterized. They therefore lend themselves to detailed studies of the determinants of ion flow through transmembraneous channels. In a comparative study, we have performed Brownian dynamics simulations to obtain statistically significant transfer efficiencies for cations and anions through matrix porin OmpF, osmoporin OmpK36, phosphoporin PhoE and two OmpF charge mutants.The simulations show that the electrostatic potential at the highly charged channel constriction serves to enhance ion permeability of either cations or anions, dependent on the type of porin. At the same time translocation of counterions is not severely impeded. At the constriction, cations and anions follow distinct trajectories, due to the segregation of basic and acidic protein residues.Simulated ion selectivity and relative conductance agree well with experimental values, and are dependent crucially on the charge constellation at the pore constriction. The experimentally observed decrease in ion selectivity and single channel conductance with increasing ionic strength is well reproduced and can be attributed to electrostatic shielding of the pore lining.     

Free energy analysis of conductivity and charge selectivity of M2GlyR-derived synthetic channels

Significant progresses have been made in the design, synthesis, modeling and in vitro testing of channel-forming peptides derived from the second transmembrane domain of the α-subunit of the glycine receptor (GlyR). The latest designs, including p22 (KKKKP ARVGL GITTV LTMTT QS), are highly soluble in water with minimal aggregation propensity and insert efficiently into cell membranes to form highly conductive ion channels. The last obstacle to a potential lead sequence for channel replacement treatment of CF patients is achieving adequate chloride selectivity. We have performed free energy simulation to analyze the conductance and charge selectivity of M2GlyR-derived synthetic channels. The results reveal that the pentameric p22 pore is non-selective. Moderate barriers for permeation of both K⁺ and Cl⁻ are dominated by the desolvation cost. Despite previous evidence suggesting a potential role of threonine side chains in anion selectivity, the hydroxyl group is not a good surrogate of water for coordinating these ions. We have also tested initial ideas of introducing additional rings of positive changes to various positions along the pore to increase anion selectivity. The results support the feasibility of achieving anion selectivity by modifying the electrostatic properties of the pore, but at the same time suggest that the peptide assembly and pore topology may also be dramatically modified, which could abolish the effects of modified electrostatics on anion selectivity. This was confirmed by subsequent two-electrode voltage clamp measurements showing that none of the tested mono-, di- and tri-Dap substituted sequences was selective. The current study thus highlights the importance of controlling channel topology besides modifying pore electrostatics for achieving anion selectivity. Several strategies are now being explored in our continued efforts to design an anion selective peptide channel with suitable biophysical, physiological and pharmacological properties as a potential treatment modality for channel replacement therapy. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

Brownian dynamics simulations of ion transport through the VDAC

It is important to gain a physical understanding of ion transport through the voltage-dependent anion channel (VDAC) because this channel provides primary permeation pathways for metabolites and electrolytes between the cytosol and mitochondria. We performed grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) simulations to explore the ion transport properties of human VDAC isoform 1 (hVDAC1; PDB:2K4T) embedded in an implicit membrane. When the MD-derived, space-dependent diffusion constant was used in the GCMC/BD simulations, the current-voltage characteristics and ion number profiles inside the pore showed excellent agreement with those calculated from all-atom molecular-dynamics (MD) simulations, thereby validating the GCMC/BD approach. Of the 20 NMR models of hVDAC1 currently available, the third one (NMR03) best reproduces both experimental single-channel conductance and ion selectivity (i.e., the reversal potential). In addition, detailed analyses of the ion trajectories, one-dimensional multi-ion potential of mean force, and protein charge distribution reveal that electrostatic interactions play an important role in the channel structure and ion transport relationship. Finally, the GCMC/BD simulations of various mutants based on NMR03 show good agreement with experimental ion selectivity. The difference in ion selectivity between the wild-type and the mutants is the result of altered potential of mean force profiles that are dominated by the electrostatic interactions.     

Permeation properties of an engineered bacterial OmpF porin containing the EEEE-locus of Ca2+ channels

The selectivity filter of the bacterial porin OmpF carries a small net charge close to -1 e and is therefore only slightly cation-selective. Calcium channels, on the other hand, contain four negatively charged glutamates, the EEEE-locus, and are among the most selective cation channels known. We aimed to turn the essentially nonselective OmpF into a Ca2+-selective channel. To that end, two additional glutamates (R42E and R132E) were introduced in the OmpF constriction zone that already contains D113 and E117. Mutant OmpF containing this DEEE-locus has a high Ca2+ over Cl- selectivity and a Na+ current with a strongly increased sensitivity to 1 mM Ca2+. The charge/space competition model, initially applied to the L-type Ca2+ channel, identifies the fixed charge and filter volume as key determinants of ion selectivity, with the precise atomic arrangement having only second-order effects. By implication, the reproduction of fixed charge and filter volume should transform two channels into channels of similar selectivity, even if the two belong to entirely different ion channel families, as is the case for OmpF and the L-type Ca2+ channel. The results presented here fit quite well in the framework of charge/space competition theory.     

Single channel analysis of conductance and rectification in cation-selective, mutant glycine receptor channels

Members of the ligand-gated ion channel superfamily mediate fast synaptic transmission in the nervous system. In this study, we investigate the molecular determinants and mechanisms of ion permeation and ion charge selectivity in this family of channels by characterizing the single channel conductance and rectification of alpha1 homomeric human glycine receptor channels (GlyRs) containing pore mutations that impart cation selectivity. The A-1'E mutant GlyR and the selectivity double mutant ([SDM], A-1'E, P-2' Delta) GlyR, had mean inward chord conductances (at -60 mV) of 7 pS and mean outward conductances of 11 and 12 pS (60 mV), respectively. This indicates that the mutations have not simply reduced anion permeability, but have replaced the previous anion conductance with a cation one. An additional mutation to neutralize the ring of positive charge at the extracellular mouth of the channel (SDM+R19'A GlyR) made the conductance-voltage relationship linear (14 pS at both 60 and -60 mV). When this external charged ring was made negative (SDM+R19'E GlyR), the inward conductance was further increased (to 22 pS) and now became sensitive to external divalent cations (being 32 pS in their absence). The effects of the mutations to the external ring of charge on conductance and rectification could be fit to a model where only the main external energy barrier height for permeation was changed. Mean outward conductances in the SDM+R19'A and SDM+R19'E GlyRs were increased when internal divalent cations were absent, consistent with the intracellular end of the pore being flanked by fixed negative charges. This supports our hypothesis that the ion charge selectivity mutations have inverted the electrostatic profile of the pore by introducing a negatively charged ring at the putative selectivity filter. These results also further confirm the role of external pore vestibule electrostatics in determining the conductance and rectification properties of the ligand-gated ion channels.     

Cation-selective mutations in the M2 domain of the inhibitory glycine receptor channel reveal determinants of ion-charge selectivity

Ligand-gated ion channel receptors mediate neuronal inhibition or excitation depending on their ion charge selectivity. An investigation into the determinants of ion charge selectivity of the anion-selective alpha1 homomeric glycine receptor (alpha1 glycine receptor [GlyR]) was undertaken using point mutations to residues lining the extra- and intracellular ends of the ion channel. Five mutant GlyRs were studied. A single substitution at the intracellular mouth of the channel (A-1'E GlyR) was sufficient to convert the channels to select cations over anions with P(Cl)/P(Na) = 0.34. This result delimits the selectivity filter and provides evidence that electrostatic interactions between permeating ions and pore residues are a critical factor in ion charge selectivity. The P-2'Delta mutant GlyR retained its anion selectivity (P(Cl)/P(Na) = 3.81), but it was much reduced compared with the wild-type (WT) GlyR (P(Cl)/P(Na) = 27.9). When the A-1'E and the P-2'Delta mutations were combined (selectivity double mutant [SDM] GlyR), the relative cation permeability was enhanced (P(Cl)/P(Na) = 0.13). The SDM GlyR was also Ca(2+) permeable (P(Ca)/P(Na) = 0.29). Neutralizing the extracellular mouth of the SDM GlyR ion channel (SDM+R19'A GlyR) produced a more Ca(2+)-permeable channel (P(Ca)/P(Na) = 0.73), without drastically altering monovalent charge selectivity (P(Cl)/P(Na) = 0.23). The SDM+R19'E GlyR, which introduces a negatively charged ring at the extracellular mouth of the channel, further enhanced Ca(2+) permeability (P(Ca)/P(Na) = 0.92), with little effect on monovalent selectivity (P(Cl)/P(Na) = 0.19). Estimates of the minimum pore diameter of the A-1'E, SDM, SDM+R19'A, and SDM+R19'E GlyRs revealed that these pores are larger than the alpha1 GlyR, with the SDM-based GlyRs being comparable in diameter to the cation-selective nicotinic acetylcholine receptors. This result provides evidence that the diameter of the ion channel is also an important factor in ion charge selectivity.     

Conductance and selectivity fluctuations in D127 mutants of the bacterial porin OmpF

A recent molecular dynamics study questioned the protonation state and physiological role of aspartate 127 (D127) of E. coli porin OmpF. To address that question we isolated two OmpF mutants with D127 either neutralized (D127N) or replaced by a positively charged lysine (D127K). The charge state of the residue at position 127 has clear effects on both conductance and selectivity. The D127K but not the D127N mutant expresses resilient conductance and selectivity fluctuations. These fluctuations reflect, we think, either changes in the ionization state of K127 and/or transitions between unstable subconformations as induced by the electrostatic repulsion between two positively charged residues, K127 and the nearby R167. Our results slightly favor the view that in WT OmpF residue D127 is deprotonated. As for the role of D127 in OmpF functionality, the gating of both mutants shows very similar sensitivity toward voltage as WT OmpF. Moreover, the current fluctuations of the D127K mutant were observed also in the absence of an applied electric field. We therefore dismiss D127 as a key residue in the control mechanism of the voltage-dependent gating of OmpF.     

Charge selectivity of the Cys-loop family of ligand-gated ion channels

The determinants of charge selectivity of the Cys-loop family of ligand-gated ion channels have been studied for more than a decade. The investigations have mainly covered homomeric receptors e.g. the nicotinic acetylcholine receptor alpha7, the glycine receptor alpha1 and the serotonin receptor 5-HT(3A). Only recently, the determinants of charge selectivity of heteromeric receptors have been addressed for the GABA(A) receptor alpha2beta3gamma2. For all receptor subtypes, the selectivity determinants have been located to an intracellular linker between transmembrane domains M1 and M2. Two features of the M1-M2 linker appear to control ion selectivity. A central role for charged amino acid residues in selectivity has been almost universally observed. Furthermore, recent studies point to an important role of the size of the narrowest constriction in the pore. In the present review, these determinants of charge selectivity of the Cys-loop family of ligand-gated ion channels will be discussed in detail.     

The Binding of Antibiotics in OmpF Porin

1. Download : Download high-res image (445KB)
2. Download : Download full-size image

Highlights? OmpF porin binds antibiotics in the extracellular and periplasmic pore vestibules ? MD simulations show that ionic current through OmpF is blocked by bound ampicillin ? Carbenicillin binding alters ion selectivity of OmpF but not total ionic current ? Disruption of binding increases the susceptibility of E. coli for antibiotics     

The ionization state of D37 in E. coli porin OmpF and the nature of conductance fluctuations in D37 mutants

Understanding the flow of ions through E. coli porin outer membrane protein F (OmpF) requires knowledge of the charge state of all titratable residues located along the permeation pathway. Earlier theoretical studies proved successful in the calculation of the pK values of most residues. The (apparent) pK of Asp37 (D37), on the other hand, appeared rather sensitive to the (unknown) protein dielectric used. We addressed the protonation state of D37 experimentally by replacing D37 with a (neutral) valine. This D37V mutant expressed reduced cation selectivity, in agreement with the view that D37 in wild-type (WT) OmpF is fully ionized, i.e., deprotonated. The introduction of a (positively charged) arginine at position 37 evoked current fluctuations. Similar behavior was observed in the D37K mutant and the cysteine mutants D37C-MTSEA and D37C-MTSET. Nontitratable [2-(trimethylammonium)ethyl]-methanethiosulfonate (MTSET) carries a permanent and pH-independent charge of 1e, implying that the fluctuations of the D37C-MTSET mutant do not represent (de)protonation reactions of MTSET. We therefore conclude that these fluctuations reflect transitions between conformational substates evoked by structural instabilities due to the positive charge at that particular position in the pore lumen. Based on the similarities between D37C-MTSET fluctuations and those seen in the other mutants, notably D37K, the underlying mechanism of these fluctuations may be (essentially) the same in all four mutants studied.     

Mechanism of Cl- selection by a glutamate-gated chloride (GluCl) receptor revealed through mutations in the selectivity filter

To learn about the mechanism of ion charge selectivity by invertebrate glutamate-gated chloride (GluCl) channels, we swapped segments between the GluClbeta receptor of Caenorhabditis elegans and the vertebrate cationic alpha7-acetylcholine receptor and monitored anionic/cationic permeability ratios. Complete conversion of the ion charge selectivity in a set of receptor microchimeras indicates that the selectivity filter of the GluClbeta receptor is created by a sequence connecting the first with the second transmembrane segments. A single substitution of a negatively charged residue within this sequence converted the selectivity of the GluClbeta receptor's pore from anionic to cationic. Unexpectedly, elimination of the charge of each basic residue of the selectivity filter, one at a time or concomitantly, moderately reduced the P(Cl)/P(Na) ratios, but the GluClbeta receptor's mutants retained high capacity to select Cl(-) over Na(+). These results indicate that, unlike the proposed case of anionic Gly- and gamma-aminobutyric acid-gated ion channels, positively charged residues do not play the key role in the selection of ionic charge by the GluClbeta receptor. Taken together with measurements of the effective open pore diameter and with structural modeling, the study presented here collectively indicates that in the most constricted part of the open GluClbeta receptor's channel, Cl(-) interacts with backbone amides, where it undergoes partial dehydration necessary for traversing the pore.     

Crystal structure of Omp32, the anion-selective porin from Comamonas acidovorans, in complex with a periplasmic peptide at 2.1 A resolution

BACKGROUND: Porins provide diffusion channels for salts and small organic molecules in the outer membrane of bacteria. In OmpF from Escherichia coli and related porins, an electrostatic field across the channel and a potential, originating from a surplus of negative charges, create moderate cation selectivity. Here, we investigate the strongly anion-selective porin Omp32 from Comamonas acidovorans, which is closely homologous to the porins of pathogenic Bordetella and Neisseria species. RESULTS: The crystal structure of Omp32 was determined to a resolution of 2.1 A using single isomorphous replacement with anomalous scattering (SIRAS). The porin consists of a 16-stranded beta barrel with eight external loops and seven periplasmic turns. Loops 3 and 8, together with a protrusion located within beta-strand 2, narrow the cross-section of the pore considerably. Arginine residues create a charge filter in the constriction zone and a positive surface potential at the external and periplasmic faces. One sulfate ion was bound to Arg38 in the channel constriction zone. A peptide of 5.8 kDa appeared bound to Omp32 in a 1:1 stoichiometry on the periplasmic side close to the symmetry axis of the trimer. Eight amino acids of this peptide could be identified, revealing specific interactions with beta-strand 1 of the porin. CONCLUSIONS: The Omp32 structure explains the strong anion selectivity of this porin. Selectivity is conferred by a positive potential, which is not attenuated by negative charges inside the channel, and by an extremely narrow constriction zone. Moreover, Omp32 represents the anchor molecule for a peptide which is homologous to proteins that link the outer membrane to the cell wall peptidoglycan.     

Anion pathway and potential energy profiles along curvilinear bacterial ClC Cl- pores: electrostatic effects of charged residues

X-ray structures permit theoretical study of Cl(-) permeation along bacterial ClC Cl(-) pores. We determined the lowest energy curvilinear pathway, identified anion-coordinating amino acids, and calculated the electrostatic potential energy profiles. We find that all four bacterial ClC Cl(-) crystal structures correspond to closed states. E148 and S107 side chains form steric barriers on both sides of the crystal binding site in the StClC wild-type and EcClC wild-type crystals; both the EcClC(E148A) and EcClC(E148Q) mutants are blocked at the S107 site. We studied the effect that mutating the charge of some strongly conserved pore-lining amino acids has on the electrostatic potential energy profiles. When E148 is neutralized, it creates an electrostatic trap, binding the ion near midmembrane. This suggests a possible electrostatic mechanism for controlling anion flow: neutralize E148, displace the side chain of E148 from the pore pathway to relieve the steric barrier, then trap the anion at midmembrane, and finally either deprotonate E148 and block the pore (pore closure) or bring a second Cl(-) into the pore to promote anion flow (pore conductance). Side-chain displacement may arise by competition for the binding site between the oxygens of E148 and the anion moving down the electrostatic energy gradient. We also find that the charge state of E111 and E113 may electrostatically control anion conductance and occupancy of the binding site within the cytoplasmic pore.     

Site selectivity in the binding of inorganic anions to serum transferrin

Equilibrium constants for the sequential binding of two anions at the specific metal-binding sites of apotransferrin have been measured by difference ultraviolet spectroscopy in 0.1 M N-2-hydroxyethylpiperazine-N'-2-ethanesulfonic acid (Hepes) at pH 7.4 and 25 degrees C. Log K1 values for phosphate, phosphite, sulfate, and arsenate fall in the narrow range of 3.5-4.0, while the log K1 for bicarbonate is 2.73. No binding is observed for nitrate, perchlorate, or borate. A dinegative charge appears to be the most important criterion for anion binding. Equilibrium constants have also been measured for binding of anions to both forms of mono(ferric)transferrin. There appears to be a very small site selectivity (0.2 to 0.4 log units) for phosphate, arsenate, and phosphite that favors binding to the N-terminal site, but there is no detectable selectivity for binding of sulfate or bicarbonate. Comparison of the binding affinities and anion selectivity with literature data on anion-binding to protonated macrocyles and cryptates strongly supports the existence of specific anion-binding sites on the protein. Binding constants were also measured in 0.01 M Hepes. The anionic sulfonate group of the buffer appears to have a small effect on anion binding.     

Distinct sensitivities of OmpF and PhoE porins to charged modulators

The inhibition of the anion-selective PhoE porin by ATP and of the cation-selective OmpF porin by polyamines has been previously documented. In the present study, we have extended the comparison of the inhibitor-porin pairs by investigating the effect of anions (ATP and aspartate) and positively charged polyamines (spermine and cadaverine) on both OmpF and PhoE with the patch-clamp technique, and by comparing directly the gating kinetics of the channels modulated by their respective substrates. The novel findings reported here are (1) that the activity of PhoE is completely unaffected by polyamines, and (2) that the kinetic changes induced by ATP on PhoE or polyamines on OmpF suggest different mechanisms of inhibition. ATP induces a high degree of flickering in the PhoE-mediated current and appears to behave as a blocker of ion flow during its presumed transport through PhoE. Polyamines modulate the kinetics of openings and closings of OmpF, in addition to promoting a blocker-like flickering activity. The strong correlation between sensitivity to inhibitors and ion selectivity suggests that some common molecular determinants are involved in these two properties and is in agreement with the hypothesis that polyamines bind inside the pore of cationic porins.     

Chemical rescue of histidine selectivity filter mutants of the M2 ion channel of influenza A virus

The influenza virus M2 proton-selective ion channel activity facilitates virus uncoating, a process that occurs in the acidic environment of the endosome. The M2 channel causes acidification of the interior of the virus particle, which results in viral protein-protein dissociation. The M2 protein is a homotetramer that contains in its aqueous pore a histidine residue (His-37) that acts as a selectivity filter and a tryptophan residue (Trp-41) that acts as a channel gate. Substitution of His-37 modifies M2 ion channel properties drastically. However, the results of such experiments are difficult to interpret because substitution of His-37 could cause gross structural changes to the channel pore. We described here experiments in which partial or, in some cases, full rescue of specific M2 ion channel properties of His-37 substitution mutants was achieved by addition of imidazole to the bathing medium. Chemical rescue was demonstrated for three histidine substitution mutant ion channels (M2-H37G, M2-H37S, and M2-H37T) and for two double mutants in which the Trp-41 channel gate was also mutated (H37G/W41Y and H37G/W41A). Currents of the M2-H37G mutant ion channel were inhibited by Cu(II), which has been shown to coordinate with His-37 in the wild-type channel. Chemical rescue was very specific for imidazole. Buffer molecules that were neutral when protonated (4-morpholineethanesulfonic acid and 3-morpholino-2-hydroxypropanesulfonic acid) did not rescue ion channel activity of the M2-H37G mutant ion channel, but 1-methylimidazole did provide partial rescue of function. These results were consistent with a model for proton transport through the pore of the wild-type channel in which the imidazole side chain of His-37 acted as an intermediate proton acceptor/donor group.     

Electrostatics and the ion selectivity of ligand-gated channels.

The nicotinic acetylcholine receptor (nAChR) is a cation-selective ion channel that opens in response to acetylcholine binding. The related glycine receptor (GlyR) is anion selective. The pore-lining domain of each protein may be modeled as a bundle of five parallel M2 helices. Models of the pore-lining domains of homopentameric nAChR and GlyR have been used in continuum electrostatics calculations to probe the origins of ion selectivity. Calculated pKA values suggest that "rings" of acidic or basic side chains at the mouths of the nAChR or GlyR M2 helix bundles, respectively, may not be fully ionized. In particular, for the nAChR the ring of glutamate side chains at the extracellular mouth of the pore is predicted to be largely protonated at neutral pH, whereas those glutamate side chains in the intracellular and intermediate rings (at the opposite mouth of the pore) are predicted to be fully ionized. Inclusion of the other domains of each protein represented as an irregular cylindrical tube in which the M2 bundles are embedded suggests that both the M2 helices and the extramembrane domains play significant roles in determining ion selectivity.     

Role of lysines in ion selectivity of bacterial outer membrane porins

The epsilon-amino groups of available lysine residues of the OmpC, OmpF and PhoE porin proteins of Escherichia coli and of the protein P porin of Pseudomonas aeruginosa, were modified by the bulky reagent trinitrobenzenesulphonic acid. Approximately 78% of the lysines of the anion-selective protein P and PhoE porins were modified whereas only 40-50% of the lysines of the cation selective OmpF and OmpC porins were altered. After modification, the three E. coli porins had very similar high selectivities for cations over anions, in contrast to the native porins which varied 86-fold in ion selectivity. Despite the large size of the trinitrophenyl group attached to modified lysines (i.e., a disc of approx. 0.86 nm diameter X 0.36 nm high) relative to the reported size of the constrictions of the E. coli porins (1.0-1.2 nm diameter), only the anion-selective PhoE porin was substantially blocked after trinitrophenylation. The protein P porin channel was relatively unaffected by trinitrophenylation, in contrast to previous data showing dramatic effects of acetylation of lysines on protein P conductance and selectivity. This favoured a model in which the critical lysines involved in anion binding by protein P were present in a constriction of the channel that was too small for trinitrobenzenesulphonic acid to enter. Overall, the data suggest that both the number and relative position of charged lysines are major determinants of ion selectivity.