Toward decrypting the allosteric mechanism of the ryanodine receptor based on coarse‐grained structural and dynamic modeling

The ryanodine receptors (RyRs) are a family of calcium (Ca) channels that regulate Ca release by undergoing a closed‐to‐open gating transition in response to action potential or Ca binding. The allosteric mechanism of RyRs gating, which is activated/regulated by ligand/protein binding >200 Å away from the channel gate, remains elusive for the lack of high‐resolution structures. Recent solution of the closed‐form structures of the RyR1 isoform by cryo‐electron microscopy has paved the way for detailed structure‐driven studies of RyRs functions. Toward elucidating the allosteric mechanism of RyRs gating, we performed coarse‐grained modeling based on the newly solved closed‐form structures of RyR1. Our normal mode analysis captured a key mode of collective motions dominating the observed structural variations in RyR1, which features large outward and downward movements of the peripheral domains with the channel remaining closed, and involves hotspot residues that overlap well with key functional sites and disease mutations. In particular, we found a key interaction between a peripheral domain and the Ca‐binding EF hand domain, which may allow for direct coupling of Ca binding to the collective motions as captured by the above mode. This key mode was robustly reproduced by the normal mode analysis of the other two closed‐form structures of RyR1 solved independently. To elucidate the closed‐to‐open conformational changes in RyR1 with amino‐acid level of details, we flexibly fitted the closed‐form structures of RyR1 into a 10‐Å cryo‐electron microscopy map of the open state. We observed extensive structural changes involving the peripheral domains and the central domains, resulting in the channel pore opening. In sum, our findings have offered unprecedented structural and dynamic insights to the allosteric mechanism of RyR1 via modulation of the key collective motions involved in RyR1 gating. The predicted hotspot residues and open‐form conformation of RyR1 will guide future mutational and functional studies. Proteins 2015; 83:2307–2318. © 2015 Wiley Periodicals, Inc.     

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.     

Pharmacological properties of native CaCCs and TMEM16A

We previously reported that TREK-1 gating by internal pH and pressure occurs close to or within the selectivity filter. These conclusions were based upon kinetic measurements of high-affinity block by quaternary ammonium (QA) ions that appeared to exhibit state-independent accessibility to their binding site within the pore. Intriguingly, recent crystal structures of two related K2P potassium channels were also both found to be open at the helix bundle crossing. However, this did not exclude the possibility of gating at the bundle crossing and it was suggested that side-fenestrations within these structures might allow state-independent access of QA ions to their binding site. In this addendum to our original study we demonstrate that even hydrophobic QA ions do not access the TREK-1 pore via these fenestrations. Furthermore, by using a chemically reactive QA ion immobilized within the pore via covalent cysteine modification we provide additional evidence that the QA binding site remains accessible to the cytoplasm in the closed state. These results support models of K2P channel gating which occur close to or within the selectivity filter and do not involve closure at the helix bundle crossing.     

State-independent intracellular access of quaternary ammonium blockers to the pore of TREK-1

We previously reported that TREK-1 gating by internal pH and pressure occurs close to or within the selectivity filter. These conclusions were based upon kinetic measurements of high-affinity block by quaternary ammonium (QA) ions that appeared to exhibit state-independent accessibility to their binding site within the pore. Intriguingly, recent crystal structures of two related K2P potassium channels were also both found to be open at the helix bundle crossing. However, this did not exclude the possibility of gating at the bundle crossing and it was suggested that side-fenestrations within these structures might allow state-independent access of QA ions to their binding site. In this addendum to our original study we demonstrate that even hydrophobic QA ions do not access the TREK-1 pore via these fenestrations. Furthermore, by using a chemically reactive QA ion immobilized within the pore via covalent cysteine modification we provide additional evidence that the QA binding site remains accessible to the cytoplasm in the closed state. These results support models of K2P channel gating which occur close to or within the selectivity filter and do not involve closure at the helix bundle crossing.     

Probing Conformational Changes during the Gating Cycle of a Potassium Channel in Lipid Bilayers

Ion conduction across the cellular membrane requires the simultaneous opening of activation and inactivation gates of the K⁺ channel pore. The bacterial KcsA channel has served as a powerful system for dissecting the structural changes that are related to four major functional states associated with K⁺ gating. Yet, the direct observation of the full gating cycle of KcsA has remained structurally elusive, and crystal structures mimicking these gating events require mutations in or stabilization of functionally relevant channel segments. Here, we found that changes in lipid composition strongly increased the KcsA open probability. This enabled us to probe all four major gating states in native-like membranes by combining electrophysiological and solid-state NMR experiments. In contrast to previous crystallographic views, we found that the selectivity filter and turret region, coupled to the surrounding bilayer, were actively involved in channel gating. The increase in overall steady-state open probability was accompanied by a reduction in activation-gate opening, underscoring the important role of the surrounding lipid bilayer in the delicate conformational coupling of the inactivation and activation gates.     

pH-Induced Conformational Change of the β-Barrel-Forming Protein OmpG Reconstituted into Native E. coli Lipids

A gating mechanism of the β-barrel-forming outer membrane protein G (OmpG) from Escherichia coli was recently presented. The mechanism was based on X-ray structures revealed from crystals grown from solubilized OmpG at both neutral pH and acidic pH. To investigate whether these conformations represent the naturally occurring gating mechanism, we reconstituted OmpG in native E. coli lipids and applied high-resolution atomic force microscopy. The reconstituted OmpG molecules assembled into both monomers and dimers. Single monomeric and dimeric OmpG molecules showed open channel entrances at pH 7.5 and at room temperature. The extracellular loops connecting the β-strands that form the transmembrane β-barrel pore exhibited elevated structural flexibility. Upon lowering the pH to 5.0, the conformation of OmpG molecules changed to close the extracellular entrance of their channel. It appears that one or more of the extracellular loops collapsed onto the channel entrance. This conformational change was fully reversible. Our data confirm that the previously reported gating mechanism of OmpG occurs at physiological conditions in E. coli lipid membranes.     

The structure of MgtE in the absence of magnesium provides new insights into channel gating

MgtE is a Mg²⁺ channel conserved in organisms ranging from prokaryotes to eukaryotes, including humans, and plays an important role in Mg²⁺ homeostasis. The previously determined MgtE structures in the Mg²⁺-bound, closed-state, and structure-based functional analyses of MgtE revealed that the binding of Mg²⁺ ions to the MgtE cytoplasmic domain induces channel inactivation to maintain Mg²⁺ homeostasis. There are no structures of the transmembrane (TM) domain for MgtE in Mg²⁺-free conditions, and the pore-opening mechanism has thus remained unclear.Here, we determined the cryo-electron microscopy (cryo-EM) structure of the MgtE-Fab complex in the absence of Mg²⁺ ions. The Mg²⁺-free MgtE TM domain structure and its comparison with the Mg²⁺-bound, closed-state structure, together with functional analyses, showed the Mg²⁺-dependent pore opening of MgtE on the cytoplasmic side and revealed the kink motions of the TM2 and TM5 helices at the glycine residues, which are important for channel activity. Overall, our work provides structure-based mechanistic insights into the channel gating of MgtE.

MgtE is a magnesium-selective ion channel whose gating is regulated by cytoplasmic magnesium concentration; this cryo-EM study reveals how MgtE undergoes magnesium-dependent structural changes to open the pore on the cytoplasmic side.     

Sensing and Responding to Membrane Tension: The Bacterial MscL Channel as a Model System

Mechanosensors are important for many life functions, including the senses of touch, balance, and proprioception; cardiovascular regulation; kidney function; and osmoregulation. Many channels from an assortment of families are now candidates for eukaryotic mechanosensors and proprioception, as well as cardiovascular regulation, kidney function, and osmoregulation. Bacteria also possess two families of mechanosensitive channels, termed MscL and MscS, that function as osmotic emergency release valves. Of the two channels, MscL is the most conserved, most streamlined in structure, and largest in conductance at 3.6 nS with a pore diameter in excess of 30 Å; hence, the structural changes required for gating are exaggerated and perhaps more easily defined. Because of these properties, as well as its tractable nature, MscL represents a excellent model for studying how a channel can sense and respond to biophysical changes of a lipid bilayer. Many of the properties of the MscL channel, such as the sensitivity to amphipaths, a helix that runs along the membrane surface and is connected to the pore via a glycine, a twisting and turning of the transmembrane domains upon gating, and the dynamic changes in membrane interactions, may be common to other candidate mechanosensors. Here we review many of these properties and discuss their structural and functional implications.     

Dynamical properties of the MscL of Escherichia coli: a normal mode analysis

The mechanosensitive channel (MscL) is an integral membrane protein which gates in response to membrane tension. Physiological data have shown that the gating transition involves a very large change in the conformation, and that the open state of the channel forms a large non-specific pore with a high conductance. The Escherichia coli channel structure was first modeled by homology modeling, starting with the X-ray structure of the homologous from Mycobacterium tuberculosis. Then, the dynamical and conformational properties of the channel were explored, using normal mode analysis. Such an analysis was also performed with the different structures proposed recently by Sukharev and co-workers. Similar dynamical behaviors are observed, which are characteristic of the channel architecture, subtle differences being due to the different relative positioning of the structural elements. The ability of particular regions of the channel to deform is discussed with respect to the functional and structural properties, implied in the gating process. Our results show that the first step of the gating mechanism can be described with three low-frequency modes only. The movement associated to these modes is clearly an iris-like movement involving both tilt and twist rotation.     

Affixing N-terminal α-helix to the wall of the voltage-dependent anion channel does not prevent its voltage gating

The voltage-dependent anion channel (VDAC) governs the free exchange of ions and metabolites between the mitochondria and the rest of the cell. The three-dimensional structure of VDAC1 reveals a channel formed by 19 β-strands and an N-terminal α-helix located near the midpoint of the pore. The position of this α-helix causes a narrowing of the cavity, but ample space for metabolite passage remains. The participation of the N-terminus of VDAC1 in the voltage-gating process has been well established, but the molecular mechanism continues to be debated; however, the majority of models entail large conformational changes of this N-terminal segment. Here we report that the pore-lining N-terminal α-helix does not undergo independent structural rearrangements during channel gating. We engineered a double Cys mutant in murine VDAC1 that cross-links the α-helix to the wall of the β-barrel pore and reconstituted the modified protein into planar lipid bilayers. The modified murine VDAC1 exhibited typical voltage gating. These results suggest that the N-terminal α-helix is located inside the pore of VDAC in the open state and remains associated with β-strand 11 of the pore wall during voltage gating.     

The pore structure and gating mechanism of K2P channels

Two-pore domain (K2P) potassium channels are important regulators of cellular electrical excitability. However, the structure of these channels and their gating mechanism, in particular the role of the bundle-crossing gate, are not well understood. Here, we report that quaternary ammonium (QA) ions bind with high-affinity deep within the pore of TREK-1 and have free access to their binding site before channel activation by intracellular pH or pressure. This demonstrates that, unlike most other K(+) channels, the bundle-crossing gate in this K2P channel is constitutively open. Furthermore, we used QA ions to probe the pore structure of TREK-1 by systematic scanning mutagenesis and comparison of these results with different possible structural models. This revealed that the TREK-1 pore most closely resembles the open-state structure of KvAP. We also found that mutations close to the selectivity filter and the nature of the permeant ion profoundly influence TREK-1 channel gating. These results demonstrate that the primary activation mechanisms in TREK-1 reside close to, or within the selectivity filter and do not involve gating at the cytoplasmic bundle crossing.     

Conserved gating hinge in ligand- and voltage-dependent K+ channels

Ion channels open and close their pore in a process called gating. On the basis of crystal structures of two voltage-independent K(+) channels, KcsA and MthK, a conformational change for gating has been proposed whereby the inner helix bends at a glycine hinge point (gating hinge) to open the pore and straightens to close it. Here we ask if a similar gating hinge conformational change underlies the mechanics of pore opening of two eukaryotic voltage-dependent K(+) channels, Shaker and BK channels. In the Shaker channel, substitution of the gating hinge glycine with alanine and several other amino acids prevents pore opening, but the ability to open is recovered if a secondary glycine is introduced at an adjacent position. A proline at the gating hinge favors the open state of the Shaker channel as if by preventing inner helix straightening. In BK channels, which have two adjacent glycine residues, opening is significantly hindered in a graded manner with single and double mutations to alanine. These results suggest that K(+) channels, whether ligand- or voltage-dependent, open when the inner helix bends at a conserved glycine gating hinge.     

A model of voltage gating developed using the KvAP channel crystal structure

Having inspected the crystal structure of the complete KvAP channel protein, we suspect that the voltage-sensing domain is too distorted to provide reliable information about its native tertiary structure or its interactions with the central pore-forming domain. On the other hand, a second crystal structure of the isolated voltage-sensing domain may well correspond to a native open conformation. We also observe that the paddle model of gating developed from these two structures is inconsistent with many experimental results, and suspect it to be energetically unrealistic. Here we show that the isolated voltage-sensing domain crystal structure can be docked onto the pore domain portion of the full-length KvAP crystal structure in an energetically favorable way to create a model of the open conformation. Using this as a starting point, we have developed rather conventional models of resting and transition conformations based on the helical screw mechanism for the transition from the open to the resting conformation. Our models are consistent with both theoretical considerations and experimental results.     

Crystal structures of a ligand-free MthK gating ring: insights into the ligand gating mechanism of K+ channels

MthK is a prokaryotic Ca(2+)-gated K(+) channel that, like other ligand-gated channels, converts the chemical energy of ligand binding to the mechanical force of channel opening. The channel's eight ligand-binding domains, the RCK domains, form an octameric gating ring in which Ca(2+) binding induces conformational changes that open the channel. Here we present the crystal structures of the MthK gating ring in closed and partially open states at 2.8 A, both obtained from the same crystal grown in the absence of Ca(2+). Furthermore, our biochemical and electrophysiological analyses demonstrate that MthK is regulated by both Ca(2+) and pH. Ca(2+) regulates the channel by changing the equilibrium of the gating ring between closed and open states, while pH regulates channel gating by affecting gating-ring stability. Our findings, along with the previously determined open MthK structure, allow us to elucidate the ligand gating mechanism of RCK-regulated K(+) channels.     

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.     

ClpP: a structurally dynamic protease regulated by AAA+ proteins

Proteolysis is an important process for many aspects of bacterial physiology. Clp proteases carry out a large proportion of protein degradation in bacteria. These enzymes assemble in complexes that combine the protease ClpP and the unfoldase, ClpA or ClpX. ClpP oligomerizes as two stacked heptameric rings enclosing a central chamber containing the proteolytic sites. ClpX and ClpA assemble into hexameric rings that bind both axial surfaces of the ClpP tetradecamer forming a barrel-like complex. ClpP requires association with ClpA or ClpX to unfold and thread protein substrates through the axial pore into the inner chamber where degradation occurs. A gating mechanism regulated by the ATPase exists at the entry of the ClpP axial pore and involves the N-terminal regions of the ClpP protomers. These gating motifs are located at the axial regions of the tetradecamer but in most crystal structures they are not visible. We also lack structural information about the ClpAP or ClpXP complexes. Therefore, the structural details of how the axial gate in ClpP is regulated by the ATPases are unknown. Here, we review our current understanding of the conformational changes that ClpA or ClpX induce in ClpP to open the axial gate and increase substrate accessibility into the degradation chamber. Most of this knowledge comes from the recent crystal structures of ClpP in complex with acyldepsipeptides (ADEP) antibiotics. These small molecules are providing new insights into the gating mechanism of this protease because they imitate the interaction of ClpA/ClpX with ClpP and activate its protease activity.     

Structures of the MthK RCK domain and the effect of Ca2+ on gating ring stability

MthK is a Ca2+-gated K+ channel from Methanobacterium autotrophicum. The crystal structure of the MthK channel in a Ca2+-bound open state was previously determined at 3.3 A and revealed an octameric gating ring composed of eight intracellular ligand-binding RCK (regulate the conductance of K+) domains. It was suggested that Ca2+ binding regulates the gating ring conformation, which in turn leads to the opening and closing of the channel. However, at 3.3 AA resolution, the molecular details of the structure are not well defined, and many of the conclusions drawn from that structure were hypothetical. Here we have presented high resolution structures of the MthK RCK domain with and without Ca2+ bound from three different crystals. These structures revealed a dimeric architecture of the RCK domain and allowed us to visualize the Ca2+ binding and protein-protein contacts at atomic detail. The dimerization of RCK domains is also conserved in other RCK-regulated K+ channels and transporters, suggesting that the RCK dimer serves as a basic unit in the gating ring assembly. A comparison of these dimer structures confirmed that the dimer interface is indeed flexible as suggested previously. However, the conformational change at the flexible interface is of an extent smaller than the previously hypothesized gating ring movement, and a reconstruction of these dimers into octamers by applying protein-protein contacts at the fixed interface did not generate enclosed gating rings. This indicated that there is a high probability that the previously defined fixed interface may not be fixed during channel gating. In addition to the structural studies, we have also carried out biochemical analyses and have shown that near physiological pH, isolated RCK domains form a stable octamer in solution, supporting the notion that the formation of octameric gating ring is a functionally relevant event in MthK gating. Additionally, our stability studies indicated that Ca2+ binding stabilizes the RCK domains in this octameric state.     

Molecular dynamics simulations of Piezo1 channel opening by increases in membrane tension

Piezo1 is a mechanosensitive channel involved in many cellular functions and responsible for sensing shear stress and pressure forces in cells. Piezo1 has a unique trilobed topology with a curved membrane region in the closed state. It has been suggested that upon activation Piezo1 adopts a flattened conformation, but the molecular and structural changes underpinning the Piezo1 gating and opening mechanisms and how the channel senses forces in the membrane remain elusive. Here, we used molecular dynamics simulations to reveal the structural rearrangements that occur when Piezo1 moves from a closed to an open state in response to increased mechanical tension applied to a model membrane. We find that membrane stretching causes Piezo1 to flatten and expand its blade region, resulting in tilting and lateral movement of the pore-lining transmembrane helices 37 and 38. This is associated with the opening of the channel and movement of lipids out of the pore region. Our results reveal that because of the rather loose packing of Piezo1 pore region, movement of the lipids outside the pore region is critical for the opening of the pore. Our simulations also suggest synchronous flattening of the Piezo1 blades during Piezo1 activation. The flattened structure lifts the C-terminal extracellular domain up, exposing it more to the extracellular space. Our studies support the idea that it is the blade region of Piezo1 that senses tension in the membrane because pore opening failed in the absence of the blades. Additionally, our simulations reveal that upon opening, water molecules occupy lateral fenestrations in the cytosolic region of Piezo1, which might be likely paths for ion permeation. Our results provide a model for how mechanical force opens the Piezo1 channel and thus how it might couple mechanical force to biological response.     

Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors

The initial coupling between ligand binding and channel gating in the human α7 nicotinic acetylcholine receptor (nAChR) has been investigated with targeted molecular dynamics (TMD) simulation. During the simulation, eight residues at the tip of the C-loop in two alternating subunits were forced to move toward a ligand-bound conformation as captured in the crystallographic structure of acetylcholine binding protein (AChBP) in complex with carbamoylcholine. Comparison of apo- and ligand-bound AChBP structures shows only minor rearrangements distal from the ligand-binding site. In contrast, comparison of apo and TMD simulation structures of the nAChR reveals significant changes toward the bottom of the ligand-binding domain. These structural rearrangements are subsequently translated to the pore domain, leading to a partly open channel within 4 ns of TMD simulation. Furthermore, we confirmed that two highly conserved residue pairs, one located near the ligand-binding pocket (Lys145 and Tyr188), and the other located toward the bottom of the ligand-binding domain (Arg206 and Glu45), are likely to play important roles in coupling agonist binding to channel gating. Overall, our simulations suggest that gating movements of the α7 receptor may involve relatively small structural changes within the ligand-binding domain, implying that the gating transition is energy-efficient and can be easily modulated by agonist binding/unbinding.     

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.