Models of ion pores in N-type voltage-gated calcium channels

Two computer models of the outer vestibule of the pore of the N-type voltage-gated Ca²⁺ channel are predicted. The models are constructed from β-hairpin peptide segments in the S5–S6 loops of each of the four domains that produce the channel. These hairpins together are modeled to form a short eight-stranded β barrel. The models contain a ring of glutamates at the base of the barrel, which have been shown by mutagenesis experiments to function as a selectivity filter. These filters are suggested by the models to be of the correct dimensions to allow the permeation of a hydrated calcium ion, where the filter glutamates may substitute for molecules of water from the hydration shell of the ion. The models also suggest that a ring of threonines and an aspartate might be present between the mouth of the pore and the filter, and hence the models may prove useful in suggesting future mutagenesis experiments.     

Determining the conductance of the SecY protein translocation channel for small molecules

The channel formed by the SecY complex must maintain the membrane barrier for ions and other small molecules during the translocation of membrane or secretory proteins. We have tested the permeability of the channel by using planar bilayers containing reconstituted purified E. coli SecY complex. Wild-type SecY complex did not show any conductance for ions or water. Deletion of the "plug," a short helix normally located in the center of the SecY complex, or modification of a cysteine introduced into the plug resulted in transient channel openings; a similar effect was seen with a mutation in the pore ring, a constriction in the center of the channel. Permanent channel opening occurred when the plug was moved out of the way by disulfide-bridge formation. These data show that the resting channel on its own forms a barrier for small molecules, with both the pore ring and the plug required for the seal; channel opening requires movement of the plug.     

The Roles of Pore Ring and Plug in the SecY Protein-conducting Channel

The protein-conducting channel, or translocon, is an evolutionarily conserved complex that allows nascent proteins to cross a cellular membrane or integrate into it. The crystal structure of an archaeal translocon, the SecY complex, revealed that two elements contribute to sealing the channel: a small “plug” domain blocking the periplasmic region of the channel, and a pore ring composed of six hydrophobic residues acting as a constriction point at the channel's center. To determine the independent functions of these two elements, we have performed molecular dynamics simulations of the native channel as well as of two recently structurally resolved mutants in which portions of their plugs were deleted. We find that in the mutants, the instability in the plug region leads to a concomitant increase in flexibility of the pore ring. The instability is quantified by the rate of water permeation in each system as well as by the force required for oligopeptide translocation. Through a novel simulation in which the interactions between the plug and water were independently controlled, we find that the role of the plug in stabilizing the pore ring is significantly more important than its role as a purely steric barrier.     

A molecular model of the inner pore of the Ca channel in its open state

Structure of the Ca channel open pore is unlikely to be the same as that of the K channel because Ca channels do not contain the hinge residues Gly or Pro. The Ca channel does not have a wide entry into the inner pore, as is found in K channels. First we sought to simulate the open state of the Ca channel by modeling forced opening of the KcsA channel using a procedure of restrained minimization with distance constraints at the level of the α-helical bundle, corresponding to segments Thr-107-Val-115. This produced an intermediate open state, which was populated by amino acid residues of Ca channels and then successively optimized until the opening of the pore reached a diameter of about 10 Å, large enough to allow verapamil to enter and block the Ca channel from inside. Although this approach produced a sterically plausible structure, it was in significant disagreement with the MTSET accessibility data for single cysteine mutations of S6 segments of the P/Q channel¹ that do not fit with an α-helical pattern. Last we explored the idea that the four S6 segments of Ca channels may contain intra-molecular deformations that lead to reorientation of its side chains. After introduction of ≠-bulges, the model agreed with the MTSET accessibility data. MTSET modification of a cysteine at the C-end of only one S6 could produce physical occlusion and block of the inner pore of the open Ca channel, as observed experimentally, and as expected if the pore opening is narrower than that of K channels.     

Molecular dynamics study of water and Na+ ions in models of the pore region of the nicotinic acetylcholine receptor.

The nicotinic acetylcholine receptor (nAChR) is an integral membrane protein that forms ligand-gated and cation-selective channels. The central pore is lined by a bundle of five approximately parallel M2 helices, one from each subunit. Candidate model structures of the solvated pore region of a homopentameric (alpha7)5 nAChR channel in the open state, and in two possible forms of the closed state, have been studied using molecular dynamics simulations with restraining potentials. It is found that the mobility of the water is substantially lower within the pore than in bulk, and the water molecules become aligned with the M2 helix dipoles. Hydrogen-bonding patterns in the pore, especially around pore-lining charged and hydrophilic residues, and around exposed regions of the helix backbone, have been determined. Initial studies of systems containing both water and sodium ions together within the pore region have also been conducted. A sodium ion has been introduced into the solvated models at various points along the pore axis and its energy profile evaluated. It is found that the ion causes only a local perturbation of the water structure. The results of these calculations have been used to examine the effectiveness of the central ring of leucines as a component of a gate in the closed-channel model.     

Cryo-EM structure of the inner ring from the Xenopus laevis nuclear pore complex

Nuclear pore complex (NPC) mediates nucleocytoplasmic shuttling. Here we present single-particle cryo-electron microscopy structure of the inner ring (IR) subunit from the Xenopus laevis NPC at an average resolution of 4.2 Å. A homo-dimer of Nup205 resides at the center of the IR subunit, flanked by two molecules of Nup188. Four molecules of Nup93 each places an extended helix into the axial groove of Nup205 or Nup188, together constituting the central scaffold. The channel nucleoporin hetero-trimer of Nup62/58/54 is anchored on the central scaffold. Six Nup155 molecules interact with the central scaffold and together with the NDC1–ALADIN hetero-dimers anchor the IR subunit to the nuclear envelope and to outer rings. The scarce inter-subunit contacts may allow sufficient latitude in conformation and diameter of the IR. Our structure reveals the molecular basis for the IR subunit assembly of a vertebrate NPC.Subject terms: Cryoelectron microscopy, Nuclear organization     

A molecular model of the inner pore of the Ca channel in its open state

Structure of the Ca channel open pore is unlikely to be the same as that of the K channel because Ca channels do not contain the hinge residues Gly or Pro. The Ca channel does not have a wide entry into the inner pore, as is found in K channels. First we sought to simulate the open state of the Ca channel by modeling forced opening of the KcsA channel using a procedure of restrained minimization with distance constraints at the level of the α-helical bundle, corresponding to segments Thr-107-Val-115. This produced an intermediate open state, which was populated by amino acid residues of Ca channels and then successively optimized until the opening of the pore reached a diameter of about 10 Å, large enough to allow verapamil to enter and block the Ca channel from inside. Although this approach produced a sterically plausible structure, it was in significant disagreement with the MTSET accessibility data for single cysteine mutations of S6 segments of the P/Q channel¹ that do not fit with an α-helical pattern. Last we explored the idea that the four S6 segments of Ca channels may contain intra-molecular deformations that lead to reorientation of its side chains. After introduction of π-bulges, the model agreed with the MTSET accessibility data. MTSET modification of a cysteine at the C-end of only one S6 could produce physical occlusion and block of the inner pore of the open Ca channel, as observed experimentally, and as expected if the pore opening is narrower than that of K channels.Key words: calcium channels, homology modeling, π-bulges, restrained minimization     

Exchange of conductance and gating properties between gap junction hemichannels.

Gap junction channels span the membranes of two adjacent cells and allow the gated transit of molecules as large as second messengers from cell to cell. The structure of the gap junction channel pore is not resolved. For identification of pore determinants we used a chimera of two connexins, cx46 and cx32E(1)43, that form membrane channels with distinct unit conductances and channel kinetics. Exchange of the first transmembrane segment (M1) between these connexins resulted in a chimera that exhibited most of the channel properties of the M1 donor, including single channel conductance, channel kinetics, and the preference to dwell at a subconductance level. The M1 segment thus appears to be an important determinant of conductance and gating properties of connexin channels.     

Voltage-dependent hydration and conduction properties of the hydrophobic pore of the mechanosensitive channel of small conductance

A detailed picture of water and ion properties in small pores is important for understanding the behavior of biological ion channels. Several recent modeling studies have shown that small, hydrophobic pores exclude water and ions even if they are physically large enough to accommodate them, a mechanism called hydrophobic gating. This mechanism has been implicated in the gating of several channels, including the mechanosensitive channel of small conductance (MscS). Although the pore in the crystal structure of MscS is wide and was initially hypothesized to be open, it is lined by hydrophobic residues and may represent a nonconducting state. Molecular dynamics simulations were performed on MscS to determine whether or not the structure can conduct ions. Unlike previous simulations of hydrophobic nanopores, electric fields were applied to this system to model the transmembrane potential, which proved to be important. Although simulations without a potential resulted in a dehydrated, occluded pore, the application of a potential increased the hydration of the pore and resulted in current flow through the channel. The calculated channel conductance was in good agreement with experiment. Therefore, it is likely that the MscS crystal structure is closer to a conducting than a nonconducting state.     

An energy-efficient gating mechanism in the acetylcholine receptor channel suggested by molecular and Brownian dynamics

Acetylcholine receptors mediate electrical signaling between nerve and muscle by opening and closing a transmembrane ion conductive pore. Molecular and Brownian dynamics simulations are used to shed light on the location and mechanism of the channel gate. Four separate 5 ns molecular dynamics simulations are carried out on the imaged structure of the channel, a hypothetical open structure with a slightly wider pore and a mutant structure in which a central ring of hydrophobic residues is replaced by polar groups. Water is found to partially evacuate the pore during molecular simulations of the imaged structure, whereas ions face a large energy barrier and do not conduct through the channel in Brownian dynamics simulations. The pore appears to be in a closed configuration despite containing an unobstructed pathway across the membrane as a series of hydrophobic residues in the center of the channel provide an unfavorable home to water and ions. When the channel is widened slightly, water floods into the channel and ions conduct at a rate comparable to the currents measured experimentally in open channels. The pore remains permeable to ions provided the extracellular end of the pore-lining helix is restrained near the putative open configuration to mimic the presence of the ligand binding domain. Replacing some of the hydrophobic residues with polar ones decreases the barrier for ion permeation but does not result in significant currents. The channel is posited to utilize an energy efficient gating mechanism in which only minor conformational changes of the hydrophobic region of the pore are required to create macroscopic changes in conductance.     

Molecular mechanism of H+ conduction in the single-file water chain of the gramicidin channel

The conduction of protons in the hydrogen-bonded chain of water molecules (or "proton wire") embedded in the lumen of gramicidin A is studied with molecular dynamics free energy simulations. The process may be described as a "hop-and-turn" or Grotthuss mechanism involving the chemical exchange (hop) of hydrogen nuclei between hydrogen-bonded water molecules arranged in single file in the lumen of the pore, and the subsequent reorganization (turn) of the hydrogen-bonded network. Accordingly, the conduction cycle is modeled by two complementary steps corresponding respectively to the translocation 1) of an ionic defect (H+) and 2) of a bonding defect along the hydrogen-bonded chain of water molecules in the pore interior. The molecular mechanism and the potential of mean force are analyzed for each of these two translocation steps. It is found that the mobility of protons in gramicidin A is essentially determined by the fine structure and the dynamic fluctuations of the hydrogen-bonded network. The translocation of H+ is mediated by spontaneous (thermal) fluctuations in the relative positions of oxygen atoms in the wire. In this diffusive mechanism, a shallow free-energy well slightly favors the presence of the excess proton near the middle of the channel. In the absence of H+, the water chain adopts either one of two polarized configurations, each of which corresponds to an oriented donor-acceptor hydrogen-bond pattern along the channel axis. Interconversion between these two conformations is an activated process that occurs through the sequential and directional reorientation of water molecules of the wire. The effect of hydrogen-bonding interactions between channel and water on proton translocation is analyzed from a comparison to the results obtained previously in a study of model nonpolar channels, in which such interactions were missing. Hydrogen-bond donation from water to the backbone carbonyl oxygen atoms lining the pore interior has a dual effect: it provides a coordination of water molecules well suited both to proton hydration and to high proton mobility, and it facilitates the slower reorientation or turn step of the Grotthuss mechanism by stabilizing intermediate configurations of the hydrogen-bonded network in which water molecules are in the process of flipping between their two preferred, polarized states. This mechanism offers a detailed molecular model for the rapid transport of protons in channels, in energy-transducing membrane proteins, and in enzymes.     

Ion coordination in the amphotericin B channel.

The antifungal polyene antibiotic amphotericin B forms channels in lipid membranes that are permeable to ions, water, and nonelectrolytes. Anion, cation, and ion pair coordination in the water-filled pore of the "barrel" unit of the channels was studied by molecular dynamics simulations. Unlike the case of the gramicidin A channel, the water molecules do not create a single-file configuration in the pore, and some cross sections of the channel contain three or four water molecules. Both the anion and cation are strongly bound to ligand groups and water molecules located in the channel. The coordination number of the ions is about six. The chloride has two binding sites in the pore. The binding with water is dominant; more than four water molecules are localized in the anion coordination sphere. Three motifs of the ion coordination were monitored. The dominant motif occurs when the anion is bound to one ligand group. The ion is bound to two or three ligand groups in the less favorable configurations. The strong affinity of cations to the channel is determined by the negatively charged ligand oxygens, whose electrostatic field dominates over the field of the hydrogens. The ligand contribution to the coordination number of the sodium ion is noticeably higher than in the case of the anion. As in the case of the anion, there are three motifs of the cation coordination. The favorable one occurs when the cation is bound to two ligand oxygens. In the less favorable cases, the cation is bound to three or four oxygens. In the contact ion pair, the cation and anion are bound to two ligand oxygens and one ligand hydrogen, respectively. There exist intermediate solvent-shared states of the ion pair. The average distances between ions in these states are twice as large as that of the contact ion pair. The stability of the solvent-shared state is defined by the water molecule oriented along the electrostatic field of both ions.     

Colloid-osmotic lysis is a general feature of the mechanism of action of Bacillus thuringiensis δ-endotoxins with different insect specificity

The mechanism of action of Bacillus thuringiensis insecticidal δ-endotoxins has long been the subject of controversy. As our working hypothesis we propose a two-step model in which, after binding a specific plasma membrane receptor, the action of all the δ-endotoxins studied here is to generate small pores in the plasma membrane, either directly by inserting into the membrane, or indirectly by perturbing resident plasma membrane molecules. The creation of these pores will lead to colloid-osmotic lysis, i.e., an equilibration of ions through the pore resulting in a net inflow of ions, an accompanying influx of water, cell swelling and eventual lysis. Our observations that cell swelling precedes lysis, that small molecules leak out of the cell before large ones, that osmotic protectants inhibit or delay cytolysis, and that the toxin-induced pore of 0.5–1.0 nm radius will allow equilibration of ions but not leakage of cytoplasmic macromolecules, are in full agreement with the predictions of this hypothesis. To explain the specificity of the δ-endotoxin-induced lytic pore formation, we propose that prior interaction between the toxin and cell-specific plasma membrane recpetors is necessary before these toxins can insert into, or interact with, the membrane.     

The Gating Mechanism of the Human Aquaporin 5 Revealed by Molecular Dynamics Simulations

Aquaporins are protein channels located across the cell membrane with the role of conducting water or other small sugar alcohol molecules (aquaglyceroporins). The high-resolution X-ray structure of the human aquaporin 5 (HsAQP5) shows that HsAQP5, as all the other known aquaporins, exhibits tetrameric structure. By means of molecular dynamics simulations we analyzed the role of spontaneous fluctuations on the structural behavior of the human AQP5. We found that different conformations within the tetramer lead to a distribution of monomeric channel structures, which can be characterized as open or closed. The switch between the two states of a channel is a tap-like mechanism at the cytoplasmic end which regulates the water passage through the pore. The channel is closed by a translation of the His67 residue inside the pore. Moreover, water permeation rate calculations revealed that the selectivity filter, located at the other end of the channel, regulates the flow rate of water molecules when the channel is open, by locally modifying the orientation of His173. Furthermore, the calculated permeation rates of a fully open channel are in good agreement with the reported experimental value.     

Quantum calculations on water in the KcsA channel cavity with permeant and non-permeant ions

Different ions in the pore of the KcsA channel behave differently, and we relate this to their solvation. We show that the selectivity is dependent, in part, on the solvation in the cavity (sometimes referred to as the vestibule, it is the region containing water molecules between the intracellular gate and the selectivity filter at the extracellular end of the pore). We have shown earlier that potassium is more dependent at the upper end of the cavity region on solvation by the threonines there, while sodium ion has more water molecules as ligands. In addition, sodium ion is placed asymmetrically, while potassium is nearly exactly symmetric with respect to the four-fold symmetry of the channel. We have now extended these calculations to rubidium and cesium ions, and find that rubidium solvation resembles that of potassium (and both are permeant ions), while cesium resembles sodium (and both are non-permeant), in terms of the geometry of up to eight hydrating, and four non-hydrating, water molecules. In each case, a maximum of 12 water molecules are relevant to the calculation. The placement of the water molecules in the two cases is essentially the same as found from the electron density in the X-ray structure of Zhou and MacKinnon. For Na⁺ and K⁺, we show that energy decreases from bulk to the cavity to the lowest position in the selectivity filter (accurate energy could not be calculated for the heavier ions). A separate calculation shows that fixing the Na⁺ ion at the position of the K⁺ minimum, followed by re-optimization produced a significantly modified system, not something that could be produced by thermal fluctuations. Moving the K⁺ into the Na⁺ position in the upper cavity led to a small increase in energy, ≈ 3 k_BT, but was accompanied by large shifts in the positions of hydrating waters, which would create a major kinetic barrier. Therefore, thermal fluctuations could not invalidate the conclusions of the main calculations.     

The SecY complex forms a channel capable of ionic discrimination

Protein translocation across the bacterial membrane occurs at the SecY complex or channel. The resting SecY channel is impermeable to small molecules owing to a plug domain that creates a seal. Here, we report that a channel loosely sealed, or with a plug locked open, does not, however, lead to general membrane permeability. Instead, strong selectivity towards small monovalent anions, especially chloride, is observed. Mutations in the pore ring‐structure increase both the translocation activity of the channel and its ionic conductance, however the selectivity is maintained. The same ionic specificity also occurs at the onset of protein translocation and across the archaeal SecY complex. Thus, the ion‐conducting characteristic of the channel seems to be conserved as a normal consequence of protein translocation. We propose that the pore ring‐structure forms a selectivity filter, allowing cells to tolerate channels with imperfect plugs.     

Water transport in human aquaporin-4: molecular dynamics (MD) simulations

Aquaporin-4 (AQP4) is the predominant water channel in the central nervous system, where it has been reported to be involved in many pathophysiological roles including water transport. In this paper, the AQP4 tetramer was modeled from its PDB structure file, embedded in a palmitoyl-oleoyl-phosphatidyl-choline (POPC) lipid bilayer, solvated in water, then minimized and equilibrated by means of molecular dynamics simulations. Analysis of the equilibrated structure showed that the central pore along the fourfold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-195, Leu-191 and Leu-75, form the narrowest part of the pore. Therefore water molecules are not expected to transport through the central pore, which was confirmed by MD simulations. Each monomer of the AQP4 tetramers forms a channel whose walls consist mostly of hydrophilic residues. There are eight water molecules in single file observed in each of the four channels, transporting through the selectivity filter containing Arg-216, His-201, Phe-77, Ala-210, and the two conserved Asn-Pro-Ala (NPA) motifs containing Asn-213 and Asn-97. By using Brownian dynamics fluctuation–dissipation-theorem (BD-FDT), the overall free-energy profile was obtained for water transporting through AQP4 for the first time, which gives a complete map of the entire channel of water permeation.     

A computational study of the closed and open states of the influenza a M2 proton channel

In this study, four possible conformations of the His-37 and Trp-41 residues for the closed state of the influenza M2 ion channel were identified by a conformation scan based on a solid-state NMR restraint. In the four conformations, the His-37 residue can be of either the t-160 or t60 rotamer, whereas Trp-41 can be of either the t-105 or t90 rotamer. These conformations were further analyzed by density functional theory calculations and molecular dynamics simulations, and the data indicate that the His-37 residue most likely adopts the t60 rotamer and should be monoprotonated at the delta-nitrogen site, whereas Trp-41 adopts the t90 rotamer. This result is consistent with published experimental data and points to a simple gating mechanism: in the closed state, the His-37 and Trp-41 residues adopt the (t60, t90) conformation, which nearly occludes the pore, preventing nonproton ions from passing through due to the steric and desolvation effects. Moreover, the His-37 tetrad interrupts the otherwise continuous hydrogen-bonding network of the pore water by forcing the water molecules above and below it to adopt opposite orientations, thus adding to the blockage of proton shuttling. The channel can be easily opened by rotating the His-37 chi2 angle from 60 to 0 degrees . This open structure allows pore water to penetrate the constrictive region and to form a continuous water wire for protons to shuttle through, while being still narrow enough to exclude other ions.     

Molecular dynamics studies of the archaeal translocon

The translocon is a protein-conducting channel conserved over all domains of life that serves to translocate proteins across or into membranes. Although this channel has been well studied for many years, the recent discovery of a high-resolution crystal structure opens up new avenues of exploration. Taking advantage of this, we performed molecular dynamics simulations of the translocon in a fully solvated lipid bilayer, examining the translocation abilities of monomeric SecYEbeta by forcing two helices comprised of different amino acid sequences to cross the channel. The simulations revealed that the so-called plug of SecYEbeta swings open during translocation, closing thereafter. Likewise, it was established that the so-called pore ring region of SecYEbeta forms an elastic, yet tight, seal around the translocating oligopeptides. The closed state of the channel was found to block permeation of all ions and water molecules; in the open state, ions were blocked. Our results suggest that the SecYEbeta monomer is capable of forming an active channel.