Ionic partition and transport in multi-ionic channels: a molecular dynamics simulation study of the OmpF bacterial porin

We performed all-atom molecular dynamics simulations studying the partition of ions and the ionic current through the bacterial porin OmpF and two selected mutants. The study is motivated by new, interesting experimental findings concerning their selectivity and conductance behavior at neutral pH. The mutations considered here are designed to study the effect of removal of negative charges present in the constriction zone of the wild-type OmpF channel (which contains, on one side, a cluster with three positive residues, and on the other side, two negatively charged residues). Our results show that these mutations induce an exclusion of cations from the constriction zone of the channel, substantially reducing the flow of cations. In fact, the partition of ions inside the mutant channels is strongly inhomogeneous, with regions containing an excess of cations and regions containing an excess of anions. Interestingly, the overall number of cations inside the channel is larger than the number of anions, this excess being different for each protein channel. We found that the differences in ionic charge inside these channels are justified by the differences in electric charge between the wild-type OmpF and the mutants, following an electroneutral balance.     

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.     

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.     

Probing tubulin-blocked state of VDAC by varying membrane surface charge

Reversible blockage of the voltage-dependent anion channel (VDAC) of the mitochondrial outer membrane by dimeric tubulin is being recognized as a potent regulator of mitochondrial respiration. The tubulin-blocked state of VDAC is impermeant for ATP but only partially closed for small ions. This residual conductance allows studying the nature of the tubulin-blocked state in single-channel reconstitution experiments. Here we probe this state by changing lipid bilayer charge from positive to neutral to negative. We find that voltage sensitivity of the tubulin-VDAC blockage practically does not depend on the lipid charge and salt concentration with the effective gating charge staying within the range of 10-14 elementary charges. At physiologically relevant low salt concentrations, the conductance of the tubulin-blocked state is decreased by positive and increased by negative charge of the lipids, whereas the conductance of the open channel is much less sensitive to this parameter. Such a behavior supports the model in which tubulin's negatively charged tail enters the VDAC pore, inverting its anionic selectivity to cationic and increasing proximity of ion pathways to the nearest lipid charges as compared with the open state of the channel.     

Topology of the selectivity filter of a TRPV channel: rapid accessibility of contiguous residues from the external medium

The transient receptor potential type V5 (TRPV5) channel is a six-transmembrane domain ion channel that is highly selective to Ca(2+). To study the topology of the selectivity filter using the substituted cysteine accessibility method (SCAM), cysteine mutants at positions 541-547 were studied as heterotetramers using dimeric constructs that couple the control channel in tandem with a cysteine-bearing subunit. Whole cell currents of dimeric constructs D542C, G543C, P544C, A545C, and Y547C were rapidly inhibited by positively charged 2-(trimethyl ammonium)methyl methane thiosulfonate bromide (MTSMT), 2-(aminoethyl)methane thiosulfonate bromide (MTSEA), and 2-(trimethyl ammonium)ethyl methane thiosulfonate bromide (MTSET) reagents, whereas D542C, P544C, and A545C were inhibited only by negatively charged sodium 2-(sulfonatoethyl)methane thiosulfonate (MTSES). In contrast, the I541C dimer remained insensitive to positive and negative reagents. However, I541C/D542G and I541C/D542N dimeric constructs were rapidly (<30 s) and strongly inhibited by positively and negatively charged methane thiosulfonate reagents, suggesting that removing two of the four carboxylate residues at position 542 disrupts a constriction point in the selectivity filter. Taken together, these results establish that the side chains of contiguous amino acids in the selectivity filter of TRPV5 are rapidly accessible from the external medium, in contrast to the three-dimensional structure of the selectivity filter in K(+) channels, where main chain carbonyls were shown to project toward a narrow permeation pathway. The I541C data further suggest that the selectivity filter of the TRPV5 channel espouses a specific conformation that restrains accessibility in the presence of four carboxylate residues at position 542.     

Site-directed mutagenesis of tyrosine 118 within the central constriction site of the LamB (Maltoporin) channel of Escherichia coli. I. Effect on ion transport

The three-dimensional structure of the malto-oligosaccharide-specific LamB-channel of Escherichia coli (also called maltoporin) is known from x-ray crystallography. The central constriction of the channel formed by the external loop 3 is controlled by a tyrosine residue (Y118). Y118 was replaced by site-directed mutagenesis by ten other amino acids (alanine, isoleucine, asparagine, serine, cysteine, aspartic acid, arginine, histidine, phenylalanine, and tryptophane) including neutral ones, negatively and positively charged amino acids to study the effect of their size, hydrophobicity, and charge on ion transport through LamB. The mutant proteins were purified to homogeneity. They were reconstituted into lipid bilayer membranes and single-channel conductance and ion selectivity were measured to get insight into the mechanism of ion transport through LamB. The mutation of Y118 to any other nonaromatic amino acid led to a substantial increase of the single-channel conductance by more than a factor of six at maximum. The highest effect was observed for Y118D. Additionally, a nonlinear relationship between the salt concentration in the aqueous phase and the channel conductance was observed for this mutant, indicating strong discrete charge effects on ion conductance. For all other mutants, with the exception of Y118R, linear relationships were found between single-channel conductance and bulk aqueous concentration. The individual hydrophobicity indices of the amino acids introduced inside the central constriction of the LamB channel had a somewhat smaller effect on the single-channel conductance as compared with the effect of their size and charge.     

Ion permeation through the alpha-hemolysin channel: theoretical studies based on Brownian dynamics and Poisson-Nernst-Plank electrodiffusion theory

Identification of the molecular interaction governing ion conduction through biological pores is one of the most important goals of modern electrophysiology. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) and three-dimensional Poisson-Nernst-Plank (3d-PNP) electrodiffusion algorithms offer powerful and general approaches to study of ion permeation through wide molecular pores. A detailed analysis of ion flows through the staphylococcal alpha-hemolysin channel based on series of simulations at different concentrations and transmembrane potentials is presented. The position-dependent diffusion coefficient is approximated on the basis of a hydrodynamic model. The channel conductance calculated by GCMC/BD is approximately 10% higher than (electrophysiologically measured) experimental values, whereas results from 3d-PNP are always 30-50% larger. Both methods are able to capture all important electrostatic interactions in equilibrium conditions. The asymmetric conductance upon the polarity of the transmembrane potential observed experimentally is reproduced by GCMC/BD and 3d-PNP. The separation of geometrical and energetic influence of the channel on ion conduction reveals that such asymmetries arise from the permanent charge distribution inside the pore. The major determinant of the asymmetry is unbalanced charge in the triad of polar residues D127, D128, and K131. The GCMC/BD or 3d-PNP calculations reproduce also experimental reversal potentials and permeability rations in asymmetric ionic solutions. The weak anionic selectivity of the channel results from the presence of the salt bridge between E111 and K147 in the constriction zone. The calculations also reproduce the experimentally derived dependence of the reversible potential to the direction of the salt gradient. The origin of such effect arises from the asymmetrical distribution of energetic barriers along the channel axis, which modulates the preferential ion passage in different directions.     

Multistep mechanism of chloride translocation in a strongly anion-selective porin channel

The strongly anion-selective porin channel Omp32 from the bacterium Delftia acidovorans differs from other unspecific porins by its pronounced selectivity for anions and its particularly small channel cross-section. Multinanosecond molecular dynamics simulations of chloride ion movement in this pore protein suggest that translocated anions interact intimately with the charges of a "basic ladder", whose dynamics lead the anions in a stepwise manner through the constriction zone of the channel. The ladder-steps comprise the central clustered arginine groups and flanking basic residues at its exoplasmic and periplasmic sides. The computed free energy profile of ion movement in and around the constriction zone shows a corresponding succession of free energy minima and barriers. A number of polar atoms from other amino acids contribute to the coordination of Cl(-) at certain sites and to its temporary immobilization in the channel. A special binding site occurs at the transition of the constriction zone to the periplasmic funnel, binding the chloride ion over significant lengths of time. The results from our MD study offer a possible explanation for the nonlinear conductance properties and unusual salt-dependent characteristics of Omp32 observed earlier in experimental measurements.     

Side-chain charge effects and conductance determinants in the pore of ClC-0 chloride channels

The charge on the side chain of the internal pore residue lysine 519 (K519) of the Torpedo ClC-0 chloride (Cl-) channel affects channel conductance. Experiments that replace wild-type (WT) lysine with neutral or negatively charged residues or that modify the K519C mutant with various methane thiosulfonate (MTS) reagents show that the conductance of the channel decreases when the charge at position 519 is made more negative. This charge effect on the channel conductance diminishes in the presence of a high intracellular Cl- concentration ([Cl-]i). However, the application of high concentrations of nonpermeant ions, such as glutamate or sulfate (SO42-), does not change the conductance, suggesting that the electrostatic effects created by the charge at position 519 are unlikely due to a surface charge mechanism. Another pore residue, glutamate 127 (E127), plays an even more critical role in controlling channel conductance. This negatively charged residue, based on the structures of the homologous bacterial ClC channels, lies 4-5 A from K519. Altering the charge of this residue can influence the apparent Cl- affinity as well as the saturated pore conductance in the conductance-Cl- activity curve. Amino acid residues at the selectivity filter also control the pore conductance but mutating these residues mainly affects the maximal pore conductance. These results suggest at least two different conductance determinants in the pore of ClC-0, consistent with the most recent crystal structure of the bacterial ClC channel solved to 2.5 A, in which multiple Cl--binding sites were identified in the pore. Thus, we suggest that the occupancy of the internal Cl--binding site is directly controlled by the charged residues located at the inner pore mouth. On the other hand, the Cl--binding site at the selectivity filter controls the exit rate of Cl- and therefore determines the maximal channel conductance.     

Comparing the Temperature-Dependent Conductance of the Two Structurally Similar E. coli Porins OmpC and OmpF

The temperature-dependent ion conductance of OmpC, a major outer membrane channel of Escherichia coli, is predicted using all-atom molecular dynamics simulations and experimentally verified. To generalize previous results, OmpC is compared to its structural homolog OmpF at different KCl concentrations, pH values, and a broad temperature range. At low salt concentrations and up to room temperature, the molecular modeling predicts the experimental conductance accurately. At high salt concentrations above 1 M KCl and above room temperature, the simulations underestimate the conductance. Moreover, the temperature dependence of the channel conductance is different from that of the bulk, both in experiment and simulation, indicating a strong contribution of surface effects to the ion conductance. With respect to OmpC, subconductance levels can be observed in experiments only. Subconductance and gating levels can be clearly distinguished by their differences in conductance values and temperature-dependent behavior. With increasing temperature, the probability of a subconductance state to occur, increases, while the dwell time is decreased. The open probability, frequency, and dwell time of such states is largely pH- and KCl concentration-independent, while their amplitudes show a lower increase with increasing salt concentration than gating amplitudes. Voltage dependence of subconductance has been found to be negligible within the uncertainty of the measurements.     

Glutamate substitution in repeat IV alters divalent and monovalent cation permeation in the heart Ca2+ channel.

In voltage-gated ion channels, residues responsible for ion selectivity were identified in the pore-lining SS1-SS2 segments. Negatively charged glutamate residues (E393, E736, E1145, and E1446) found in each of the four repeats of the alpha 1C subunit were identified as the major determinant of selectivity in Ca2+ channels. Neutralization of glutamate residues by glutamine in repeat I (E393Q), repeat III (E1145Q), and repeat IV (E1446Q) decreased the channel affinity for calcium ions 10-fold from the wild-type channel. In contrast, neutralization of glutamate residues in repeat II failed to significantly alter Ca2+ affinity. Likewise, mutation of neighboring residues in E1149K and D1450N did not affect the channel affinity, further supporting the unique role of glutamate residues E1145 in repeat III and E1446 in repeat IV in determining Ca2+ selectivity. Conservative mutations E1145D and E1446D preserved high-affinity Ca2+ binding, which suggests that the interaction between Ca2+ and the pore ligand sites is predominantly electrostatic and involves charge neutralization. Mutational analysis of E1446 showed additionally that polar residues could achieve higher Ca2+ affinity than small hydrophobic residues could. The role of high-affinity calcium binding sites in channel permeation was investigated at the single-channel level. Neutralization of glutamate residue in repeats I, II, and III did not affect single-channel properties measured with 115 mM BaCl2. However, mutation of the high-affinity binding site E1446 was found to significantly affect the single-channel conductance for Ba2+ and Li+, providing strong evidence that E1446 is located in the narrow region of the channel outer mouth. Side-chain substitutions at 1446 in repeat IV were used to probe the nature of divalent cation-ligand interaction and monovalent cation-ligand interaction in the calcium channel pore. Monovalent permeation was found to be inversely proportional to the volume of the side chain at position 1446, with small neutral residues such as alanine and glycine producing higher Li+ currents than the wild-type channel. This suggests that steric hindrance is a major determinant for monovalent cation conductance. Divalent permeation was more complex. Ba2+ single-channel conductance decreased when small neutral residues such as glycine were replaced by bulkier ones such as glutamine. However, negatively charged amino acids produced single-channel conductance higher than predicted from the size of their side chain. Hence, negatively charged residues at position 1446 in repeat IV are required for divalent cation permeation.     

Probing the role of negatively charged amino acid residues in ion permeation of skeletal muscle ryanodine receptor

Sequence comparison suggests that the ryanodine receptors (RyRs) have pore architecture similar to that of the bacterial K+ channel KcsA. The lumenal loop linking the two most C-terminal transmembrane spanning segments in the RyRs has a predicted pore helix and an amino acid motif (GGGIG) similar to the selectivity filter (TVGYG) of KcsA identified by x-ray analysis. The RyRs have many negatively charged amino acid residues in the two regions linking the GGGIG motif and predicted pore helix with the two most C-terminal transmembrane spanning segments. We tested the role of these residues by generating single-site mutants, focusing on amino acid residues conserved among the mammalian RyRs. Replacement of two acidic residues immediately after the GGGIG motif in skeletal muscle ryanodine receptor (RyR1-D4899 and -E4900) with asparagine and glutamine profoundly affected ion permeation and selectivity. By comparison, mutagenesis of aspartate and glutamate residues in the putative linker regions showed a K+ conductance and selectivity for Ca2+ compared to K+ (P(Ca)/P(K)) close to wild-type. The results show that the negatively charged carboxyl oxygens of D4899 and E4900 side chains are major determinants of RyR ion conductance and selectivity.     

Positive charges at the intracellular mouth of the pore regulate anion conduction in the CFTR chloride channel

Many different ion channel pores are thought to have charged amino acid residues clustered around their entrances. The so-called surface charges contributed by these residues can play important roles in attracting oppositely charged ions from the bulk solution on one side of the membrane, increasing effective local counterion concentration and favoring rapid ion movement through the channel. Here we use site-directed mutagenesis to identify arginine residues contributing important surface charges in the intracellular mouth of the cystic fibrosis transmembrane conductance regulator (CFTR) Cl(-) channel pore. While wild-type CFTR was associated with a linear current-voltage relationship with symmetrical solutions, strong outward rectification was observed after mutagenesis of two arginine residues (R303 and R352) located near the intracellular ends of the fifth and sixth transmembrane regions. Current rectification was dependent on the charge present at these positions, consistent with an electrostatic effect. Furthermore, mutagenesis-induced rectification was more pronounced at lower Cl(-) concentrations, suggesting that these mutants had a reduced ability to concentrate Cl(-) ions near the inner pore mouth. R303 and R352 mutants exhibited reduced single channel conductance, especially at negative membrane potentials, that was dependent on the charge of the amino acid residue present at these positions. However, the very low conductance of both R303E and R352E-CFTR could be greatly increased by elevating intracellular Cl(-) concentration. Modification of an introduced cysteine residue at position 303 by charged methanethiosulfonate reagents reproduced charge-dependent effects on current rectification. Mutagenesis of arginine residues in the second and tenth transmembrane regions also altered channel permeation properties, however these effects were not consistent with changes in channel surface charges. These results suggest that positively charged arginine residues act to concentrate Cl(-) ions at the inner mouth of the CFTR pore, and that this contributes to maximization of the rate of Cl(-) ion permeation through the pore.     

Bioinformatic and mutational analysis of channelrhodopsin-2 protein cation-conducting pathway

Channelrhodopsin-2 (ChR2) is a light-gated cation channel widely used as a biotechnological tool to control membrane depolarization in various cell types and tissues. Although several ChR2 variants with modified properties have been generated, the structural determinants of the protein function are largely unresolved. We used bioinformatic modeling of the ChR2 structure to identify the putative cationic pathway within the channel, which is formed by a system of inner cavities that are uniquely present in this microbial rhodopsin. Site-directed mutagenesis combined with patch clamp analysis in HeLa cells was used to determine key residues involved in ChR2 conductance and selectivity. Among them, Gln-56 is important for ion conductance, whereas Ser-63, Thr-250, and Asn-258 are previously unrecognized residues involved in ion selectivity and photocurrent kinetics. This study widens the current structural information on ChR2 and can assist in the design of new improved variants for specific biological applications.     

X-ray crystallographic and mass spectrometric structure determination and functional characterization of succinylated porin from Rhodobacter capsulatus: implications for ion selectivity and single-channel conductance.

The role of charges near the pore mouth has been discussed in theoretical work about ion channels. To introduce new negative charges in a channel protein, amino groups of porin from Rhodobacter capsulatus 37b4 were succinylated with succinic anhydride, and the precise extent and sites of succinylations and structures of the succinylporins determined by mass spectrometry and X-ray crystallography. Molecular weight and peptide mapping analyses using matrix-assisted laser desorption-ionization mass spectrometry identified selective succinylation of three lysine-epsilon-amino groups (Lys-46, Lys-298, Lys-300) and the N-terminal alpha-amino group. The structure of a tetra-succinylated porin (TS-porin) was determined to 2.4 A and was generally found unchanged in comparison to native porin to form a trimeric complex. All succinylated amino groups found in a mono/di-succinylated porin (MS-porin) and a TS-porin are localized at the inner channel surface and are solvent-accessible: Lys-46 is located at the channel constriction site, whereas Lys-298, Lys-300, and the N-terminus are all near the periplasmic entrance of the channel. The Lys-46 residue at the central constriction loop was modeled as succinyl-lysine from the electron density data and shown to bend toward the periplasmic pore mouth. The electrical properties of the MS-and TS-porins were determined by reconstitution into black lipid membranes, and showed a negative charge effect on ion transport and an increased cation selectivity through the porin channel. The properties of a typical general diffusion porin changed to those of a channel that contains point charges near the pore mouth. The single-channel conductance was no longer a linear function of the bulk aqueous salt concentration. The substantially higher cation selectivity of the succinylated porins compared with the native protein is consistent with the increase of negatively charged groups introduced. These results show tertiary structure-selective modification of charged residues as an efficient approach in the structure-function evaluation of ion channels, and X-ray crystallography and mass spectrometry as complementary analytical tools for defining precisely the chemically modified structures.     

Conducting-state properties of the KcsA potassium channel from molecular and Brownian dynamics simulations.

The mechanisms underlying transport of ions across the potassium channel are examined using electrostatic calculations and three-dimensional Brownian dynamics simulations. We first build open-state configurations of the channel with molecular dynamics simulations, by pulling the transmembrane helices outward until the channel attains the desired interior radius. To gain insights into ion permeation, we construct potential energy profiles experienced by an ion traversing the channel in the presence of other resident ions. These profiles reveal that in the absence of an applied field the channel accommodates three potassium ions in a stable equilibrium, two in the selectivity filter and one in the central cavity. In the presence of a driving potential, this three-ion state becomes unstable, and ion permeation across the channel is observed. These qualitative explanations are confirmed by the results of three-dimensional Brownian dynamics simulations. We find that the channel conducts when the ionizable residues near the extracellular entrance are fully charged and those near the intracellular side are partially charged. The conductance increases steeply as the radius of the intracellular mouth of the channel is increased from 2 A to 5 A. Our simulation results reproduce several experimental observations, including the current-voltage curves, conductance-concentration relationships, and outward rectification of currents.     

Acidification Asymmetrically Affects Voltage-dependent Anion Channel Implicating the Involvement of Salt Bridges

The voltage-dependent anion channel (VDAC) is the major pathway for ATP, ADP, and other respiratory substrates through the mitochondrial outer membrane, constituting a crucial point of mitochondrial metabolism regulation. VDAC is characterized by its ability to “gate” between an open and several “closed” states under applied voltage. In the early stages of tumorigenesis or during ischemia, partial or total absence of oxygen supply to cells results in cytosolic acidification. Motivated by these facts, we investigated the effects of pH variations on VDAC gating properties. We reconstituted VDAC into planar lipid membranes and found that acidification reversibly increases its voltage-dependent gating. Furthermore, both VDAC anion selectivity and single channel conductance increased with acidification, in agreement with the titration of the negatively charged VDAC residues at low pH values. Analysis of the pH dependences of the gating and open channel parameters yielded similar pK_a values close to 4.0. We also found that the response of VDAC gating to acidification was highly asymmetric. The presumably cytosolic (cis) side of the channel was the most sensitive to acidification, whereas the mitochondrial intermembrane space (trans) side barely responded to pH changes. Molecular dynamic simulations suggested that stable salt bridges at the cis side, which are susceptible to disruption upon acidification, contribute to this asymmetry. The pronounced sensitivity of the cis side to pH variations found here in vitro might provide helpful insights into the regulatory role of VDAC in the protective effect of cytosolic acidification during ischemia in vivo.     

Molecular dynamics investigation of the ω-current in the Kv1.2 voltage sensor domains

Voltage sensor domains (VSD) are transmembrane proteins that respond to changes in membrane voltage and modulate the activity of ion channels, enzymes, or in the case of proton channels allow permeation of protons across the cell membrane. VSDs consist of four transmembrane segments, S1-S4, forming an antiparallel helical bundle. The S4 segment contains several positively charged residues, mainly arginines, located at every third position along the helix. In the voltage-gated Shaker K(+) channel, the mutation of the first arginine of S4 to a smaller uncharged amino acid allows permeation of cations through the VSD. These currents, known as ω-currents, pass through the VSD and are distinct from K(+) currents passing through the main ion conduction pore. Here we report molecular dynamics simulations of the ω-current in the resting-state conformation for Kv1.2 and for four of its mutants. The four tested mutants exhibit various degrees of conductivity for K(+) and Cl(-) ions, with a slight selectivity for K(+) over Cl(-). Analysis of the ion permeation pathway, in the case of a highly conductive mutant, reveals a negatively charged constriction region near the center of the membrane that might act as a selectivity filter to prevent permeation of anions through the pore. The residues R1 in S4 and E1 in S2 are located at the narrowest region of the ω-pore for the resting state conformation of the VSD, in agreement with experiments showing that the largest increase in current is produced by the double mutation E1D and R1S.     

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.     

Binding modes of μ-conotoxin to the bacterial sodium channel (NaVAb)

Polypeptide toxins isolated from the venom of cone snails, known as μ-conotoxins, block voltage-gated sodium channels by physically occluding the ion-conducting pathway. Using molecular dynamics, we show that one subtype of μ-conotoxins, PIIIA, effectively blocks the bacterial voltage-gated sodium channel Na_VAb, whose crystal structure has recently been elucidated. The spherically shaped toxin, carrying a net charge of +6 e with six basic residues protruding from its surface, is attracted by the negatively charged residues on the vestibular wall and the selectivity filter of the channel. The side chain of each of these six arginine and lysine residues can wedge into the selectivity filter, whereas the side chains of other basic residues form electrostatic complexes with two acidic residues on the channel. We construct the profile of potential of mean force for the unbinding of PIIIA from the channel, and predict that PIIIA blocks the bacterial sodium channel with subnanomolar affinity.