Role of Aromatic Localization in the Gating Process of a Potassium Channel

Position of the transmembrane aromatic residues of the KirBac1.1 potassium channel shifts from an even distribution in the closed state toward the membrane/solute interface in the open state model. This is the first example of an integral membrane protein making use of the observed preference for transmembrane aromatic residues to reside at the interfaces. The process of aromatic localization is proposed as a means of directing and stabilizing structural changes during conformational transitions within the transmembrane region of integral membrane proteins. All-atom molecular dynamics simulations of the open and closed conformers in a membrane environment have been carried out to take account of the interactions between the aromatic residues and the lipids, which may be involved in the conformational change, e.g., the gating of the channel.     

The closed structure of the MscS mechanosensitive channel. Cross-linking of single cysteine mutants

Mechanosensitive channels must make a large conformational change during the transition from the closed to the open state. The crystal structure of the open form of the Escherichia coli MscS channel was recently solved and depicts a homoheptamer (1). In this study, cross-linking of site-specific cysteine substitutions demonstrates that residues up to 10-33 A apart in the crystal structure readily form disulfide bridges in the closed form and can also be cross-linked by a 10-A linker. Cross-linking between adjacent subunits stabilizes the heptameric form of the channel providing biochemical evidence to support the crystal structure. The data are consistent with the published model (1) in that the membrane domain is highly flexible and that the closed to open transition may involve a significant displacement of transmembrane helices 1 and 2, possibly by as much as 30 A. The data are also consistent with significant flexibility of the cytoplasmic domain.     

Gating Transition of Pentameric Ligand-Gated Ion Channels

Pentameric ligand-gated ion channels are an important family of membrane proteins and play key roles in physiological processes, including signal transduction at chemical synapses. Here, we study the conformational changes associated with the opening and closing of the channel pore. Based on recent crystal structures of two prokaryotic members of the family in open and closed states, respectively, mixed elastic network models are constructed for the transmembrane domain. To explore the conformational changes in the gating transition, a coarse-grained transition path is computed that smoothly connects the closed and open conformations of the channel. We find that the conformational transition involves no major rotations of the transmembrane helices, and is instead characterized by a concerted tilting of helices M2 and M3. In addition, helix M2 changes its bending state, which results in an early closure of the pore during the open-to-closed transition.     

New insights of membrane environment effects on MscL channel mechanics from theoretical approaches

The prokaryotic mechanosensitive channel of large conductance (MscL) is a remarkable integral membrane protein. During hypo-osmotic shock, it responses to membrane tension through large conformational changes, that lead to an open state of the pore. The structure of the channel from Mycobacterium tuberculosis has been resolved in the closed state. Numerous experiments have attempted to trap the channel in its open state but they did not succeed in obtaining a structure. A gating mechanism has been proposed based on different experimental data but there is no experimental technique available to follow this process in atomic details. In addition, it has been shown that a decrease of the lipid bilayer thickness lowered MscL activation energy and stabilized a structurally distinct closed channel intermediate. Here, we use atomistic molecular dynamics simulations to investigate the effect of the lipid bilayer thinning on our model of the structure of the Escherichia coli. We thoroughly analyze simulations of the channel embedded in two pre-equilibrated membranes differing by their hydrophobic tail length (DMPE and POPE). The MscL structure remains stable in POPE, whereas a distinct structural state is obtained in DMPE in response to hydrophobic mismatch. This latter is obtained by tilts and kinks of the transmembrane helices, leading to a widening and a diminution of the channel height. Part of these motions is guided by a competition between solvent and lipids for the interaction with the periplasmic loops. We finally conduct a principal component analysis of the simulation and compare anharmonic motions with harmonic ones, previously obtained from a coarse-grained normal mode analysis performed on the same structural model. Significant similarities exist between low-frequency harmonic motions and those observed with essential dynamics in DMPE. In summary, change in membrane thickness permits to accelerate the conformational changes involved in the mechanics of the E. coli channel, providing a closed structural intermediate en route to the open state. These results give clues for better understanding why the channel activation energy is lowered in a thinner membrane.     

Conformational changes in a pore-lining helix coupled to cystic fibrosis transmembrane conductance regulator channel gating

Cystic fibrosis transmembrane conductance regulator (CFTR), the protein dysfunctional in cystic fibrosis, is unique among ATP-binding cassette transporters in that it functions as an ion channel. In CFTR, ATP binding opens the channel, and its subsequent hydrolysis causes channel closure. We studied the conformational changes in the pore-lining sixth transmembrane segment upon ATP binding by measuring state-dependent changes in accessibility of substituted cysteines to methanethiosulfonate reagents. Modification rates of three residues (resides 331, 333, and 335) near the extracellular side were 10-1000-fold slower in the open state than in the closed state. Introduction of a charged residue by chemical modification at two of these positions (resides 331 and 333) affected CFTR single-channel gating. In contrast, modifications of pore-lining residues 334 and 338 were not state-dependent. Our results suggest that ATP binding induces a modest conformational change in the sixth transmembrane segment, and this conformational change is coupled to the gating mechanism that regulates ion conduction. These results may establish a structural basis of gating involving the dynamic rearrangement of transmembrane domains necessary for vectorial transport of substrates in ATP-binding cassette transporters.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

Toward elucidating the heat activation mechanism of the TRPV1 channel gating by molecular dynamics simulation

As a key cellular sensor, the TRPV1 cation channel undergoes a gating transition from a closed state to an open state in response to various physical and chemical stimuli including noxious heat. Despite years of study, the heat activation mechanism of TRPV1 gating remains enigmatic at the molecular level. Toward elucidating the structural and energetic basis of TRPV1 gating, we have performed molecular dynamics (MD) simulations (with cumulative simulation time of 3 μs), starting from the high‐resolution closed and open structures of TRPV1 solved by cryo‐electron microscopy. In the closed‐state simulations at 30°C, we observed a stably closed channel constricted at the lower gate (near residue I679), while the upper gate (near residues G643 and M644) is dynamic and undergoes flickery opening/closing. In the open‐state simulations at 60°C, we found higher conformational variation consistent with a large entropy increase required for the heat activation, and both the lower and upper gates are dynamic with transient opening/closing. Through ensemble‐based structural analyses of the closed state versus the open state, we revealed pronounced closed‐to‐open conformational changes involving the membrane proximal domain (MPD) linker, the outer pore, and the TRP helix, which are accompanied by breaking/forming of a network of closed/open‐state specific hydrogen bonds. By comparing the closed‐state simulations at 30°C and 60°C, we observed heat‐activated conformational changes in the MPD linker, the outer pore, and the TRP helix that resemble the closed‐to‐open conformational changes, along with partial formation of the open‐state specific hydrogen bonds. Some of the residues involved in the above key hydrogen bonds were validated by previous mutational studies. Taken together, our MD simulations have offered rich structural and dynamic details beyond the static structures of TRPV1, and promising targets for future mutagenesis and functional studies of the TRPV1 channel. Proteins 2016; 84:1938–1949. © 2016 Wiley Periodicals, Inc.     

Hinge-bending motions in the pore domain of a bacterial voltage-gated sodium channel

Computational methods and experimental data are used to provide structural models for NaChBac, the homo-tetrameric voltage-gated sodium channel from the bacterium Bacillus halodurans, with a closed and partially open pore domain. Molecular dynamic (MD) simulations on membrane-bound homo-tetrameric NaChBac structures, each comprising six helical transmembrane segments (labeled S1 through S6), reveal that the shape of the lumen, which is defined by the bundle of four alpha-helical S6 segments, is modulated by hinge bending motions around the S6 glycine residues. Mutation of these glycine residues into proline and alanine affects, respectively, the structure and conformational flexibility of the S6 bundle. In the closed channel conformation, a cluster of stacked phenylalanine residues from the four S6 helices hinders diffusion of water molecules and Na(+) ions. Activation of the voltage sensor domains causes destabilization of the aforementioned cluster of phenylalanines, leading to a more open structure. The conformational change involving the phenylalanine cluster promotes a kink in S6, suggesting that channel gating likely results from the combined action of hinge-bending motions of the S6 bundle and concerted reorientation of the aromatic phenylalanine side-chains.     

Defining the Role of the Tension Sensor in the Mechanosensitive Channel of Small Conductance

Mutations that alter the phenotypic behavior of the Escherichia coli mechanosensitive channel of small conductance (MscS) have been identified; however, most of these residues play critical roles in the transition between the closed and open states of the channel and are not directly involved in lipid interactions that transduce the tension response. In this study, we use molecular dynamic simulations to predict critical lipid interacting residues in the closed state of MscS. The physiological role of these residues was then investigated by performing osmotic downshock assays on MscS mutants where the lipid interacting residues were mutated to alanine. These experiments identified seven residues in the first and second transmembrane helices as lipid-sensing residues. The majority of these residues are hydrophobic amino acids located near the extracellular interface of the membrane. All of these residues interact strongly with the lipid bilayer in the closed state of MscS, but do not face the bilayer directly in structures associated with the open and desensitized states of the channel. Thus, the position of these residues relative to the lipid membrane appears related to the ability of the channel to sense tension in its different physiological states.     

Loss-of-function mutations at the rim of the funnel of mechanosensitive channel MscL

MscL is a bacterial mechanosensitive channel that is activated directly by membrane stretch. Although the gene has been cloned and the crystal structure of the closed channel has been defined, how membrane tension causes conformational changes in MscL remains largely unknown. To identify the site where MscL senses membrane tension, we examined the function of the mutants generated by random and scanning mutagenesis. In vitro (patch-clamp) and in vivo (hypoosmotic-shock) experiments showed that when a hydrophilic amino acid replaces one of the hydrophobic residues that are thought to make contact with the membrane lipid near the periplasmic end of the M1 or M2 transmembrane domain, MscL loses the ability to open in response to membrane tension. Hydrophilic (asparagine) substitution of the other residues in the lipid-protein interface did not impair the channel's mechanosensitivity. These observations suggest that the disturbance of the hydrophobic interaction between the membrane lipid and the periplasmic rim of the channel's funnel impairs the function of MscL.     

Structure of the KcsA channel intracellular gate in the open state

Ion channels catalyze the selective transfer of ions across the membrane in response to a variety of stimuli. These channels gate by controlling the access of ions to a centrally located water-filled pore. The crystal structure of the Streptomyces lividans potassium channel (KcsA) has allowed a molecular exploration of this mechanism. Electron paramagnetic resonance (EPR) studies have uncovered significant conformational changes at the intracellular end of the second transmembrane helix (TM2) upon gating. We have used site-directed spin labeling (SDSL) and EPR spectroscopy in an attempt to quantify the structural rearrangements of the KcsA TM2 bundle underlying the transition from the closed to the open state. Under conditions favoring the closed and open conformations, 10 intersubunit distances were obtained across TM2 segments from tandem dimer constructs. Analysis of these data points to a mechanism in which each TM2 helix tilts away from the permeation pathway, towards the membrane plane, and rotates about its helical axis, supporting a scissoring-type motion with a pivot point near residues 107-108. These movements are accompanied by a large increase in the diameter of the vestibule below the central water-filled cavity.     

Changes in local S4 environment provide a voltage-sensing mechanism for mammalian hyperpolarization-activated HCN channels

The positively charged S4 transmembrane segment of voltage-gated channels is thought to function as the voltage sensor by moving charge through the membrane electric field in response to depolarization. Here we studied S4 movements in the mammalian HCN pacemaker channels. Unlike most voltage-gated channel family members that are activated by depolarization, HCN channels are activated by hyperpolarization. We determined the reactivity of the charged sulfhydryl-modifying reagent, MTSET, with substituted cysteine (Cys) residues along the HCN1 S4 segment. Using an HCN1 channel engineered to be MTS resistant except for the chosen S4 Cys substitution, we determined the reactivity of 12 S4 residues to external or internal MTSET application in either the closed or open state of the channel. Cys substitutions in the NH2-terminal half of S4 only reacted with external MTSET; the rates of reactivity were rapid, regardless of whether the channel was open or closed. In contrast, Cys substitutions in the COOH-terminal half of S4 selectively reacted with internal MTSET when the channel was open. In the open state, the boundary between externally and internally accessible residues was remarkably narrow (approximately 3 residues). This suggests that S4 lies in a water-filled gating canal with a very narrow barrier between the external and internal solutions, similar to depolarization-gated channels. However, the pattern of reactivity is incompatible with either classical gating models, which postulate a large translational or rotational movement of S4 within a gating canal, or with a recent model in which S4 forms a peripheral voltage-sensing paddle (with S3b) that moves within the lipid bilayer (the KvAP model). Rather, we suggest that voltage sensing is due to a rearrangement in transmembrane segments surrounding S4, leading to a collapse of an internal gating canal upon channel closure that alters the shape of the membrane field around a relatively static S4 segment.     

Free-energy landscapes of ion-channel gating are malleable: changes in the number of bound ligands are accompanied by changes in the location of the transition state in acetylcholine-receptor channels

Acetylcholine-receptor channels (AChRs) are allosteric membrane proteins that mediate synaptic transmission by alternatively opening and closing ("gating") a cation-selective transmembrane pore. Although ligand binding is not required for the channel to open, the binding of agonists (for example, acetylcholine) increases the closed right harpoon over left harpoon open equilibrium constant because the ion-impermeable --> ion-permeable transition of the ion pathway is accompanied by a low-affinity --> high-affinity change at the agonist-binding sites. The fact that the gating conformational change of muscle AChRs can be kinetically modeled as a two-state reaction has paved the way to the experimental characterization of the corresponding transition state, which represents a snapshot of the continuous sequence of molecular events separating the closed and open states. Previous studies of fully (di) liganded AChRs, combining single-channel kinetic measurements, site-directed mutagenesis, and data analysis in the framework of the linear free-energy relationships of physical organic chemistry, have suggested a transition-state structure that is consistent with channel opening being an asynchronous conformational change that starts at the extracellular agonist-binding sites and propagates toward the intracellular end of the pore. In this paper, I characterize the gating transition state of unliganded AChRs, and report a remarkable difference: unlike that of diliganded gating, the unliganded transition state is not a hybrid of the closed- and open-state structures but, rather, is almost indistinguishable from the open state itself. This displacement of the transition state along the reaction coordinate obscures the mechanism underlying the unliganded closed right harpoon over left harpoon open reaction but brings to light the malleable nature of free-energy landscapes of ion-channel gating.     

Cystic Fibrosis Transmembrane Conductance Regulator (CFTR): CLOSED AND OPEN STATE CHANNEL MODELS*

The cystic fibrosis transmembrane conductance regulator (CFTR) is a member of the ATP-binding cassette (ABC) transporter superfamily. CFTR controls the flow of anions through the apical membrane of epithelia. Dysfunctional CFTR causes the common lethal genetic disease cystic fibrosis. Transitions between open and closed states of CFTR are regulated by ATP binding and hydrolysis on the cytosolic nucleotide binding domains, which are coupled with the transmembrane (TM) domains forming the pathway for anion permeation. Lack of structural data hampers a global understanding of CFTR and thus the development of “rational” approaches directly targeting defective CFTR. In this work, we explored possible conformational states of the CFTR gating cycle by means of homology modeling. As templates, we used structures of homologous ABC transporters, namely TM(287–288), ABC-B10, McjD, and Sav1866. In the light of published experimental results, structural analysis of the transmembrane cavity suggests that the TM(287–288)-based CFTR model could correspond to a commonly occupied closed state, whereas the McjD-based model could represent an open state. The models capture the important role played by Phe-337 as a filter/gating residue and provide structural information on the conformational transition from closed to open channel.     

Solution X-ray scattering evidence for agonist- and antagonist-induced modulation of cleft closure in a glutamate receptor ligand-binding domain

Agonist-induced conformational changes in the ligand-binding domains (LBD) of glutamate receptor ion channels provide the driving force for molecular rearrangements that mediate channel opening and subsequent desensitization. The resulting regulated transmembrane ion fluxes form the basis for most excitatory neuronal signaling in the brain. Crystallographic analysis of the GluR2 LBD core has revealed a ligand-binding cleft located between two lobes. Channel antagonists stabilize an open cleft, whereas agonists stabilize a closed cleft. The crystal structure of the apo form is similar to the antagonist-bound, open state. To understand the conformational behavior of the LBD in the absence of crystal lattice constraints, and thus better to appreciate the thermodynamic constraints on ligand binding, we have undertaken a solution x-ray scattering study using two different constructs encoding either the core or an extended LBD. In agreement with the GluR2 crystal structures, the LBD is more compact in the presence of agonist than it is in the presence of antagonist. However, the time-averaged conformation of the ligand-free core in solution is intermediate between the open, antagonist-bound state and the closed, agonist-bound state, suggesting a conformational equilibrium. Addition of peptide moieties that connect the core domain to the other functional domains in each channel subunit appears to constrain the conformational equilibrium in favor of the open state.     

Multiscale simulation of transmembrane proteins

Multiscale simulation is employed to examine changes in atomistic-level protein structure due to long wavelength membrane undulations and plane stress fields. An ensemble of atomistic-level simulations of a model of a transmembrane influenza A virus M2 proton channel in a dimyristoylphosphatidylcholine (DMPC) bilayer is coupled to a corresponding mesoscopic model of a DMPC bilayer in an explicit mesoscopic solvent. Structural variations in the key proton gating His37 residues of the M2 channel are examined. Small, but distinct variations in the structure of the His37 residues are observed in both the open and closed states of the channel as a result of the coupling to mesoscopic-level membrane motions.     

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.     

Interaction of anesthetics with open and closed conformations of a potassium channel studied via molecular dynamics and normal mode analysis

A variety of experiments suggest that membrane proteins are important targets of anesthetic molecules, and that ion channels interact differently with anesthetics in their open and closed conformations. The availability of an open and a closed structural model for the KirBac1.1 potassium channel has made it possible to perform a comparative analysis of the interactions of anesthetics with the same channel in its open and closed states. To this end, all-atom molecular dynamics simulations supplemented by normal mode analysis have been employed to probe the interactions of the inhalational anesthetic halothane with both an open and closed conformer of KirBac1.1 embedded in a lipid bilayer. Normal mode analysis on the closed and open channel, in the presence and absence of halothane, reveals that the anesthetic modulates the global as well as the local dynamics of both conformations differently. In the case of the open channel, the observed reduction of flexibility of residues in the inner helices suggests a functional modification action of anesthetics on ion channels. In this context, preferential quenching of the aromatic residue motion and modulation of global dynamics by halothane may be seen as steps toward potentiating or favoring open state conformations. These molecular dynamics simulations provide the first insights into possible specific interactions between anesthetic molecules and ion channels in different conformations.     

Mutations in the S4 region isolate the final voltage-dependent cooperative step in potassium channel activation

Charged residues in the S4 transmembrane segment play a key role in determining the sensitivity of voltage-gated ion channels to changes in voltage across the cell membrane. However, cooperative interactions between subunits also affect the voltage dependence of channel opening, and these interactions can be altered by making substitutions at uncharged residues in the S4 region. We have studied the activation of two mutant Shaker channels that have different S4 amino acid sequences, ILT (V369I, I372L, and S376T) and Shaw S4 (the S4 of Drosophila Shaw substituted into Shaker), and yet have very similar ionic current properties. Both mutations affect cooperativity, making a cooperative transition in the activation pathway rate limiting and shifting it to very positive voltages, but analysis of gating and ionic current recordings reveals that the ILT and Shaw S4 mutant channels have different activation pathways. Analysis of gating currents suggests that the dominant effect of the ILT mutation is to make the final cooperative transition to the open state of the channel rate limiting in an activation pathway that otherwise resembles that of Shaker. The charge movement associated with the final gating transition in ILT activation can be measured as an isolated component of charge movement in the voltage range of channel opening and accounts for 13% ( approximately 1.8 e0) of the total charge moved in the ILT activation pathway. The remainder of the ILT gating charge (87%) moves at negative voltages, where channels do not open, and confirms the presence of Shaker-like conformational changes between closed states in the activation pathway. In contrast to ILT, the activation pathway of Shaw S4 seems to involve a single cooperative charge-moving step between a closed and an open state. We cannot detect any voltage-dependent transitions between closed states for Shaw S4. Restoring basic residues that are missing in Shaw S4 (R1, R2, and K7) rescues charge movement between closed states in the activation pathway, but does not alter the voltage dependence of the rate-limiting transition in activation.     

Conformational changes in S6 coupled to the opening of cyclic nucleotide-gated channels.

In cyclic nucleotide-gated channels (CNG), direct binding of cyclic nucleotides in the carboxy-terminal region is allosterically coupled to opening of the pore. A CNG1 channel pore was probed using site-directed cysteine substitution to elucidate conformational changes associated with channel opening. The effects of cysteine modification on permeation suggest a structural homology between CNG and KcsA pores. We found that intersubunit disulfide bonds form spontaneously between S399C residues in the helix bundle when channels are in the closed but not in the open state. While MTSET modification of pore-lining residues was state dependent, Ag(+) modification of V391C, in the inner vestibule, occurred at the same diffusion-limited rate in both open and closed states. Our results suggest that the helix bundle undergoes a conformational change associated with gating but is not the activation gate for CNG channels.