How pH opens a H+ channel: the gating mechanism of influenza A M2

The tetrameric M2 protein from influenza A is one of the simplest pH-gated H+ channels known, offering the potential of structurally characterizing its gating mechanism. Since the only ionizable groups in the pore are four histidines, we investigated the stability and dynamics of all six possible protonation states of the protein by using molecular dynamics. We show that while all channel protonation states are surprisingly stable, only systems with two or more charged histidines are appreciably conductive. The structural switch, from a uniprotonated to a biprotonated channel, causes an electrostatic repulsion between the charged histidines that pushes the helices apart. This results in the formation of a continuous water file that conducts protons via a H+ wire. pKa calculations place this transition at a pH of 5.6, in remarkable agreement with the experimental value. Since the conversion from uniprotonation to biprotonation occurs during endosome acidification, this explains how M2 is activated in vivo.     

Charge dissociation and the rising phase of gating currents

Gating currents from voltage-sensitive channels are generally attributed to the translocation or redistribution of ionic charge associated with the channel molecule. Such charge moves in the direction of the applied field to produce a decreasing current in the external circuit. An early rising phase for the gating current is observed for a number of channel systems and might be either some special kinetic redistribution of charge or an experimental artifact. A model that produces net charge in the channel through sequential molecular dissociation of a charged channel segment gives a rising phase for the gating current. Translocation of the charged segment produces the decay phase for a biphasic gating current. This kinetic molecular explanation constitutes a physical explanation for the biphasic gating currents that is consistent with present views of channel structure.     

Molecular dynamics simulations investigate the mechanism of Psalmotoxin 1 regulating gating process of an acid-sensing ion channel 1a at pH 5.5

Acid-sensing ion channel 1a (ASIC1a) is a cation channel activated by protons and causes neuronal death through central nervous system. Psalmotoxin1 (PcTx1) is a gating modifier for ASIC1a. The process of PcTx1 regulating the channel gating from the extracellular domain to the transmembrane domain is unclear. Here we used molecular dynamics (MD) simulations method to investigate how PcTx1 regulates the gating of the ASIC1a. Our results indicated that PcTx1can mainly regulate ASIC1a gating process through hydrogen bonds, which can affect their relative positions of several key domains in ASIC1a, further, a long-range conformational changes path was determined, which is composed of β1, β2, β10, α6, α7, β11, and β12 in ASIC1a.     

Electrophoretic mobility of DNA in gels. I. New biased reptation theory including fluctuations

In weak electric fields, the motion of DNA molecules undergoing gel electrophoresis may be described by biased reptation. We argue that the degree of molecular orientation induced by the field is determined by a competition between longitudinal fluctuations and drift of the molecule along the tube. A self-consistent calculation shows that the end-to-end vector of long molecules varies with the square root of the field strength, and not linearly as previously supposed. This leads to a number of new predictions about the field dependence of the molecular mobility and the size limit of resolution. We present the results of computer simulations that support the predictions of the theory of biased reptation including fluctuations. Finally, we discuss the correspondence with experimental data and the implications of our findings for the optimization of DNA electrophoresis. © 1994 John Wiley & Sons, Inc.     

Intrinsic aqueduct orifices facilitate K+ channel gating

The ion-conducting pore of potassium channels, which can open and close to regulate ion passage, was at long thought to be a one-dimensional pore structure with a water-filled central cavity. Here, we find four orifices in the KcsA potassium channel, which are perpendicular to the pore and stretch out from the cavity. Equilibrium molecular dynamics simulations show that water molecules can flow between the cavity and orifices. Targeted molecular dynamics simulations show that during the opening process, water molecules can move into the cavity through the orifices to facilitate channel gating, whereas blocking the aqueduct orifices makes the channel difficult to open.     

A combined coarse-grained and all-atom simulation of TRPV1 channel gating and heat activation

The transient receptor potential (TRP) channels act as key sensors of various chemical and physical stimuli in eukaryotic cells. Despite years of study, the molecular mechanisms of TRP channel activation remain unclear. To elucidate the structural, dynamic, and energetic basis of gating in TRPV1 (a founding member of the TRPV subfamily), we performed coarse-grained modeling and all-atom molecular dynamics (MD) simulation based on the recently solved high resolution structures of the open and closed form of TRPV1. Our coarse-grained normal mode analysis captures two key modes of collective motions involved in the TRPV1 gating transition, featuring a quaternary twist motion of the transmembrane domains (TMDs) relative to the intracellular domains (ICDs). Our transition pathway modeling predicts a sequence of structural movements that propagate from the ICDs to the TMDs via key interface domains (including the membrane proximal domain and the C-terminal domain), leading to sequential opening of the selectivity filter followed by the lower gate in the channel pore (confirmed by modeling conformational changes induced by the activation of ICDs). The above findings of coarse-grained modeling are robust to perturbation by lipids. Finally, our MD simulation of the ICD identifies key residues that contribute differently to the nonpolar energy of the open and closed state, and these residues are predicted to control the temperature sensitivity of TRPV1 gating. These computational predictions offer new insights to the mechanism for heat activation of TRPV1 gating, and will guide our future electrophysiology and mutagenesis studies.     

Inherent Dynamics of the Acid-Sensing Ion Channel 1 Correlates with the Gating Mechanism

The acid-sensing ion channel 1 (ASIC1) is a key receptor for extracellular protons. Although numerous structural and functional studies have been performed on this channel, the structural dynamics underlying the gating mechanism remains unknown. We used normal mode analysis, mutagenesis, and electrophysiological methods to explore the relationship between the inherent dynamics of ASIC1 and its gating mechanism. Here we show that a series of collective motions among the domains and subdomains of ASIC1 correlate with its acid-sensing function. The normal mode analysis result reveals that the intrinsic rotation of the extracellular domain and the collective motions between the thumb and finger induced by proton binding drive the receptor to experience a deformation from the extracellular domain to the transmembrane domain, triggering the channel pore to undergo “twist-to-open” motions. The movements in the transmembrane domain indicate that the likely position of the channel gate is around Leu440. These motion modes are compatible with a wide body of our complementary mutations and electrophysiological data. This study provides the dynamic fundamentals of ASIC1 gating.     

Conformational dynamics of helix S6 from Shaker potassium channel: simulation studies

Prolines in transmembrane (TM) alpha-helices are believed to play an important structural and/or functional role in membrane proteins. At a structural level a proline residue distorts alpha-helical structure due to the loss of at least one stabilizing backbone hydrogen bond, and introduces flexibility in the helix that may result in substantial kink and swivel motions about the effective "hinge." At a functional level, for example in Kv channels, it is believed that proline-induced molecular hinges may have a direct role in gating, i.e., the conformational change linked to opening/closing the channel to movement of ions. In this article we study the conformational dynamics of the S6 TM helix from of the Kv channel Shaker, which possesses the motif PVP--a motif that is conserved in Kv channels. We perform multiple molecular dynamics simulations of single S6 helices in a membrane-mimetic environment in order to effectively map the kink-swivel conformational space of the protein, exploiting the ability of multiple simulations to achieve greater sampling. We show that the presence of proline locally perturbs the helix, disrupting local dihedral angles and producing local twist and unwinding in the region of the hinge--an effect that is relaxed with distance from the PVP motif. We furthermore show that motions about the hinge are highly anisotropic, reflecting a preferred region of kink-swivel conformation space that may have implications for the gating process.     

Mechanics of channel gating of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a key molecule involved in the propagation of signals in the central nervous system and peripheral synapses. Although numerous computational and experimental studies have been performed on this receptor, the structural dynamics of the receptor underlying the gating mechanism is still unclear. To address the mechanical fundamentals of nAChR gating, both conventional molecular dynamics (CMD) and steered rotation molecular dynamics (SRMD) simulations have been conducted on the cryo-electron microscopy (cryo-EM) structure of nAChR embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer and water molecules. A 30-ns CMD simulation revealed a collective motion amongst C-loops, M1, and M2 helices. The inward movement of C-loops accompanying the shrinking of acetylcholine (ACh) binding pockets induced an inward and upward motion of the outer β-sheet composed of β9 and β10 strands, which in turn causes M1 and M2 to undergo anticlockwise motions around the pore axis. Rotational motion of the entire receptor around the pore axis and twisting motions among extracellular (EC), transmembrane (TM), and intracellular MA domains were also detected by the CMD simulation. Moreover, M2 helices undergo a local twisting motion synthesized by their bending vibration and rotation. The hinge of either twisting motion or bending vibration is located at the middle of M2, possibly the gate of the receptor. A complementary twisting-to-open motion throughout the receptor was detected by a normal mode analysis (NMA). To mimic the pulsive action of ACh binding, nonequilibrium MD simulations were performed by using the SRMD method developed in one of our laboratories. The result confirmed all the motions derived from the CMD simulation and NMA. In addition, the SRMD simulation indicated that the channel may undergo an open-close (O ↔ C) motion. The present MD simulations explore the structural dynamics of the receptor under its gating process and provide a new insight into the gating mechanism of nAChR at the atomic level.     

H bonding at the helix-bundle crossing controls gating in Kir potassium channels

Specific stimuli such as intracellular H+ and phosphoinositides (e.g., PIP2) gate inwardly rectifying potassium (Kir) channels by controlling the reversible transition between the closed and open states. This gating mechanism underlies many aspects of Kir channel physiology and pathophysiology; however, its structural basis is not well understood. Here, we demonstrate that H+ and PIP2 use a conserved gating mechanism defined by similar structural changes in the transmembrane (TM) helices and the selectivity filter. Our data support a model in which the gating motion of the TM helices is controlled by an intrasubunit hydrogen bond between TM1 and TM2 at the helix-bundle crossing, and we show that this defines a common gating motif in the Kir channel superfamily. Furthermore, we show that this proposed H-bonding interaction determines Kir channel pH sensitivity, pH and PIP2 gating kinetics, as well as a K+-dependent inactivation process at the selectivity filter and therefore many of the key regulatory mechanisms of Kir channel physiology.     

Decrypting the Heat Activation Mechanism of TRPV1 Channel by Molecular Dynamics Simulation

As a prototype cellular sensor, the TRPV1 cation channel undergoes a closed-to-open gating transition in response to various physical and chemical stimuli including noxious heat. Despite recent progress, the molecular mechanism of heat activation of TRPV1 gating remains enigmatic. Toward decrypting the structural basis of TRPV1 heat activation, we performed extensive molecular dynamics simulations (with cumulative simulation time of ～11 μs) for the wild-type channel and a constitutively active double mutant at different temperatures (30, 60, and 72°C), starting from a high-resolution closed-channel structure of TRPV1 solved by cryo-electron microscopy. In the wild-type simulations, we observed heat-activated conformational changes (e.g., expansion or contraction) in various key domains of TRPV1 (e.g., the S2-S3 and S4-S5 linkers) to prime the channel for gating. These conformational changes involve a number of dynamic hydrogen-bond interactions that were validated with previous mutational studies. Next, our mutant simulations observed channel opening after a series of conformational changes that propagate from the channel periphery to the channel pore via key intermediate domains (including the S2-S3 and S4-S5 linkers). The gating transition is accompanied by a large increase in the protein-water electrostatic interaction energy, which supports the contribution of desolvation of polar/charged residues to the temperature-sensitive TRPV1 gating. Taken together, our molecular dynamics simulations and analyses offered, to our knowledge, new structural, dynamic, and energetic information to guide future mutagenesis and functional studies of the TRPV1 channels and development of TRPV1-targeting drugs.     

Bacterial mechanosensitive channels as a paradigm for mechanosensory transduction

Research on bacterial mechanosensitive (MS) channels has since their discovery been at the forefront of the MS channel field due to extensive studies of the structure and function of MscL and MscS, two of the several different types of MS channels found in bacteria. Just a few years after these two MS channels were cloned their 3D structure was solved by X-ray crystallography. Today, the repertoire of multidisciplinary approaches used in experimental and theoretical studies following the cloning and crystallographic determination of the MscL and MscS structure has expanded by including electronparamagnetic resonance (EPR) and F?rster resonance energy transfer (FRET) spectroscopy aided by computational modelling employing molecular dynamics as well as Brownian dynamics simulations, which significantly advanced the understanding of structural determinants of the gating and conduction properties of these two MS channels. These extensive multidisciplinary studies of MscL and MscS have greatly contributed to elucidation of the basic physical principles of MS channel gating by mechanical force. This review summarizes briefly the major experimental and conceptual advancements, which helped in establishing MscL and MscS as a major paradigm of mechanosensory transduction in living cells.     

Inhibition mechanism of the acetylcholine receptor by alpha-neurotoxins as revealed by normal-mode dynamics

The nicotinic acetylcholine receptor (AChR) is the prototype of ligand-gated ion channels. Here, we calculate the dynamics of the muscle AChR using normal modes. The calculations reveal a twist-like gating motion responsible for channel opening. The ion channel diameter is shown to increase with this twist motion. Strikingly, the twist motion and the increase in channel diameter are not observed for the AChR in complex with two alpha-bungarotoxin (alphaBTX) molecules. The toxins seems to lock together neighboring receptor subunits, thereby inhibiting channel opening. Interestingly, one alphaBTX molecule suffices to prevent the twist motion. These results shed light on the gating mechanism of the AChR and present a complementary inhibition mechanism by snake-venom-derived alpha-neurotoxins.     

Computer simulation of the KvAP voltage-gated potassium channel: steered molecular dynamics of the voltage sensor

The recent crystal structures of the voltage-gated potassium channel KvAP and its isolated voltage-sensing 'paddle' (composed of segments S1-S4) challenge existing models of voltage gating and raise a number of questions about the structure of the physiologically relevant state. We investigate a possible gating mechanism based on the crystal structures in a 10 ns steered molecular dynamics simulation of KvAP in a membrane-mimetic octane layer. The structure of the full KvAP protein has been modified by restraining the S2-S4 domain to the conformation of the isolated high-resolution paddle structure. After an initial relaxation, the paddle tips are pulled through the membrane from the intracellular to the extracellular side, corresponding to a putative change from closed to open. We describe the effect of this large-scale motion on the central pore domain, which remains largely unchanged, on the protein hydrogen-bonding network and on solvent. We analyze the motion of the S3b-S4 portion of the protein and propose a possible coupling mechanism between the paddle motion and the opening of the channel. Interactions between the arginine residues in S4, solvent and chloride ions are likely to play a role in the gating charge.     

Computational study of the transmembrane domain of the acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein whose transmembrane domain (TM-domain) is believed to be responsible for channel gating via a hydrophobic effect. In this work, we perform molecular dynamics and Brownian dynamics simulations to investigate the effect of transmembrane potential on the conformation and water occupancy of TM-domain, and the resulting ion permeation events. The results show that the behavior of the hydrophobic gate is voltage-dependent. Large hyperpolarized membrane potential can change the conformation of TM-domain and water occupancy in this region, which may enable ion conduction. An electrostatic gating mechanism is also proposed from our simulations, which seems to play a role in addition to the well-known hydrophobic effect.     

Multiscale modeling shows that dielectric differences make NaV channels faster than KV channels

The generation of action potentials in excitable cells requires different activation kinetics of voltage-gated Na (Na_V) and K (K_V) channels. Na_V channels activate much faster and allow the initial Na⁺ influx that generates the depolarizing phase and propagates the signal. Recent experimental results suggest that the molecular basis for this kinetic difference is an amino acid side chain located in the gating pore of the voltage sensor domain, which is a highly conserved isoleucine in K_V channels but an equally highly conserved threonine in Na_V channels. Mutagenesis suggests that the hydrophobicity of this side chain in Shaker K_V channels regulates the energetic barrier that gating charges cross as they move through the gating pore and control the rate of channel opening. We use a multiscale modeling approach to test this hypothesis. We use high-resolution molecular dynamics to study the effect of the mutation on polarization charge within the gating pore. We then incorporate these results in a lower-resolution model of voltage gating to predict the effect of the mutation on the movement of gating charges. The predictions of our hierarchical model are fully consistent with the tested hypothesis, thus suggesting that the faster activation kinetics of Na_V channels comes from a stronger dielectric polarization by threonine (Na_V channel) produced as the first gating charge enters the gating pore compared with isoleucine (K_V channel).     

Statistical properties of ion channel records. Part II: estimation from the macroscopic current

Macroscopic ion channel current is the summation of the stochastic records of individual channel currents and therefore relates to their statistical properties. As a consequence of this relationship, it may be possible to derive certain statistical properties of single channel records or even generate some estimates of the records themselves from the macroscopic current when the direct measurement of single channel currents is not applicable. We present a procedure for generating the single channel records of an ion channel from its macroscopic current when the stochastic process of channel gating has the following two properties: (I) the open duration is independent of the time of opening event and has a single exponential probability density function (pdf), (II) all the channels have the same probability to open at time t. The application of this procedure is considered for cases where direct measurement of single channel records is difficult or impossible. First, the probability density function (pdf) of opening events, a statistical property of single channel records, is derived from the normalized macroscopic current and mean channel open duration. Second, it is shown that under the conditions (I) and (II), a non-stationary Markov model can represent the stochastic process of channel gating. Third, the non-stationary Markov model is calibrated using the results of the first step. The non-stationary formulation increases the model ability to generate a variety of different single channel records compared to common stationary Markov models. The model is then used to generate single channel records and to obtain other statistical properties of the records. Experimental single channel records of inactivating BK potassium channels are used to evaluate how accurately this procedure reconstructs measured single channel sweeps.