A BEST example of channel structure annotation by molecular simulation

Exploring the Activation Process of the Glycine Receptor

Glycine receptors (GlyR) conduct inhibitory glycinergic neurotransmission in the spinal cord and the brainstem. They play an important role in muscle tone, motor coordination, respiration, and pain perception. However, the mechanism underlying GlyR activation remains unclear. There are five potential glycine binding sites in α1 GlyR, and different binding patterns may cause distinct activation or desensitization behaviors. In this study, we investigated the coupling of protein conformational changes and glycine binding events to elucidate the influence of binding patterns on the activation and desensitization processes of α1 GlyRs. Subsequently, we explored the energetic distinctions between the apical and lateral pathways during α1 GlyR conduction to identify the pivotal factors in the ion conduction pathway preference. Moreover, we predicted the mutational effects of the key residues and verified our predictions using electrophysiological experiments. For the mutants that can be activated by glycine, the predictions of the mutational directions were all correct. The strength of the mutational effects was assessed using Pearson's correlation coefficient, yielding a value of -0.77 between the calculated highest energy barriers and experimental maximum current amplitudes. These findings contribute to our understanding of GlyR activation, identify the key residues of GlyRs, and provide guidance for mechanistic studies on other pLGICs.

Electronic Polarizability Tunes the Function of the Human Bestrophin 1 Cl- Channel

Mechanisms of anion permeation within ion channels and nanopores remain poorly understood. Recent cryo-electron microscopy structures of the human bestrophin 1 Cl- channel (hBest1) provide an opportunity to evaluate ion interactions predicted by molecular dynamics (MD) simulations against experimental observations. Here, we implement the fully polarizable forcefield AMOEBA in MD simulations on different conformations of hBest1. This forcefield models multipole moments up to the quadrupole; therefore, it captures induced dipole and anion-π interactions. We show that key biophysical properties of the channel can only be simulated when electronic polarization is included in the molecular models and that Cl- permeation through the neck of the pore is achieved through hydrophobic solvation concomitant with partial ion dehydration. Furthermore, we demonstrate how such polarizable simulations can help determine the identity of ion-like densities within high-resolution cryo-EM structures and that neglecting polarization places Cl- at positions that do not correspond with their experimentally resolved location. Overall, our results demonstrate the importance of including electronic polarization in realistic and physically accurate models of biological systems, especially channels and pores that selectively permeate anions.

Predicting Mutational Effects on Ca2+-Activated Chloride Conduction of TMEM16A Based on a Simulation Study

The dysfunction and defects of ion channels are associated with many human diseases, especially for loss-of-function mutations in ion channels such as cystic fibrosis transmembrane conductance regulator mutations in cystic fibrosis. Understanding ion channels is of great current importance for both medical and fundamental purposes. Such an understanding should include the ability to predict mutational effects and describe functional and mechanistic effects. In this work, we introduce an approach to predict mutational effects based on kinetic information (including reaction barriers and transition state locations) obtained by studying the working mechanism of target proteins. Specifically, we take the Ca2+-activated chloride channel TMEM16A as an example and utilize the computational biology model to predict the mutational effects of key residues. Encouragingly, we verified our predictions through electrophysiological experiments, demonstrating a 94% prediction accuracy regarding mutational directions. The mutational strength assessed by Pearson's correlation coefficient is -0.80 between our calculations and the experimental results. These findings suggest that the proposed methodology is reliable and can provide valuable guidance for revealing functional mechanisms and identifying key residues of the TMEM16A channel. The proposed approach can be extended to a broad scope of biophysical systems.

When is a hydrophobic gate not a hydrophobic gate?

The flux of ions through a channel is most commonly regulated by changes that result in steric occlusion of its pore. However, ion permeation can also be prevented by formation of a desolvation barrier created by hydrophobic residues that line the pore. As a result of relatively minor structural changes, confined hydrophobic regions in channels may undergo transitions between wet and dry states to gate the pore closed without physical constriction of the permeation pathway. This concept is referred to as hydrophobic gating, and many examples of this process have been demonstrated. However, the term is also now being used in a much broader context that often deviates from its original meaning. In this Viewpoint, we explore the formal definition of a hydrophobic gate, discuss examples of this process compared with other gating mechanisms that simply exploit hydrophobic residues and/or lipids in steric closure of the pore, and describe the best practice for identification of a hydrophobic gate.

Using Metadynamics To Explore the Free Energy of Dewetting in Biologically Relevant Nanopores

Water confined within hydrophobic spaces can undergo cooperative dewetting transitions due to slight changes in water density and pressure that push water toward the vapor phase. Many transmembrane protein ion channels contain nanoscale hydrophobic pores that could undergo dewetting transitions, sometimes blocking the flow of ions without physical blockages. Standard molecular dynamics simulations have been extensively applied to study the behavior of water in nanoscale pores, but the large free energy barriers of dewetting often prevent direct sampling of both wet and dry states and quantitative studies of the hydration thermodynamics of biologically relevant pores. Here, we describe a metadynamics protocol that uses the number of waters within the pore as the collective variable to drive many reversible transitions between relevant hydration states and calculate well-converged free energy profiles of pore hydration. By creating model nanopore systems and changing their radius and morphology and including various cosolvents, we quantify how these pore properties and cosolvents affect the dewetting transition. The results reveal that the dewetting free energy of nanoscale pores is determined by two key thermodynamic parameters, namely, the effective surface tension coefficients of water-air and water-pore interfaces. Importantly, while the effect of salt can be fully captured in the water activity dependence, amphipathic cosolvents such as alcohols modify both dry and wet states of the pore and dramatically shift the wet-dry equilibrium. The metadynamics approach could be applied to studies of dewetting transitions within nanoscale pores of proteins and provide new insights into why different pore properties evolved in biological systems.

Structure and Function of the Bestrophin family of calcium-activated chloride channels

Bestrophins are a family of calcium-activated chloride channels (CaCCs) with relevance to human physiology and a myriad of eye diseases termed "bestrophinopathies". Since the identification of bestrophins as CaCCs nearly two decades ago, extensive studies from electrophysiological and structural biology perspectives have sought to define their key channel features including calcium sensing, gating, inactivation, and anion selectivity. The initial X-ray crystallography studies on the prokaryotic homolog of Best1, Klebsiella pneumoniae (KpBest), and the Best1 homolog from Gallus gallus (chicken Best1, cBest1), laid the foundational groundwork for establishing the architecture of Best1. Recent progress utilizing single-particle cryogenic electron microscopy has further elucidated the molecular mechanism of gating in cBest1 and, separately, the structure of Best2 from Bos taurus (bovine Best2, bBest2). Meanwhile, whole-cell patch clamp, planar lipid bilayer, and other electrophysiologic analyses using these models as well as the human Best1 (hBest1) have provided ample evidence describing the functional properties of the bestrophin channels. This review seeks to consolidate these structural and functional results to paint a broad picture of the underlying mechanisms comprising the bestrophin family's structure-function relationship.

Hydrophobic dewetting in gating and regulation of transmembrane protein ion channels

Water is at the heart of almost all biological phenomena, without which no life that we know of would have been possible. It is a misleadingly complex liquid that exists in near coexistence with the vapor phase under ambient conditions. Confinement within a hydrophobic cavity can tip this balance enough to drive a cooperative dewetting transition. For a nanometer-scale pore, the dewetting transition leads to a stable dry state that is physically open but impermeable to ions. This phenomenon is often referred to as hydrophobic gating. Numerous transmembrane protein ion channels have now been observed to utilize hydrophobic gating in their activation and regulation. Here, we review recent theoretical, simulation, and experimental studies that together have started to establish the principles of hydrophobic gating and discuss how channels of various sizes, topologies, and biological functions can utilize these principles to control the thermodynamic properties of water within their interior pores for gating and regulation. Exciting opportunities remain in multiple areas, particularly on direct experimental detection of hydrophobic dewetting in biological channels and on understanding how the cell may control the hydrophobic gating in regulation of ion channels.

Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions

Ion channels are proteins which form gated nanopores in biological membranes. Many channels exhibit hydrophobic gating, whereby functional closure of a pore occurs by local dewetting. The pentameric ligand gated ion channels (pLGICs) provide a biologically important example of hydrophobic gating. Molecular simulation studies comparing additive vs polarizable models indicate predictions of hydrophobic gating are robust to the model employed. However, polarizable models suggest favorable interactions of hydrophobic pore-lining regions with chloride ions, of relevance to both synthetic carriers and channel proteins. Electrowetting of a closed pLGIC hydrophobic gate requires too high a voltage to occur physiologically but may inform designs for switchable nanopores. Global analysis of ∼200 channels yields a simple heuristic for structure-based prediction of (closed) hydrophobic gates. Simulation-based analysis is shown to provide an aid to interpretation of functional states of new channel structures. These studies indicate the importance of understanding the behavior of water and ions within the nanoconfined environment presented by ion channels.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

Structural basis for ion selectivity in TMEM175 K+ channels

The TMEM175 family constitutes recently discovered K+channels that are important for autophagosome turnover and lysosomal pH regulation and are associated with the early onset of Parkinson Disease. TMEM175 channels lack a P-loop selectivity filter, a hallmark of all known K+ channels, raising the question how selectivity is achieved. Here, we report the X-ray structure of a closed bacterial TMEM175 channel in complex with a nanobody fusion-protein disclosing bound K+ ions. Our analysis revealed that a highly conserved layer of threonine residues in the pore conveys a basal K+ selectivity. An additional layer comprising two serines in human TMEM175 increases selectivity further and renders this channel sensitive to 4-aminopyridine and Zn2+. Our findings suggest that large hydrophobic side chains occlude the pore, forming a physical gate, and that channel opening by iris-like motions simultaneously relocates the gate and exposes the otherwise concealed selectivity filter to the pore lumen.

Conduction through a narrow inward-rectifier K+ channel pore

G-protein–gated inwardly rectifying potassium channels are important mediators of inhibitory neurotransmission. Based on microsecond-scale molecular dynamics simulations, Bernsteiner et al. propose novel gating details that may enable K⁺ flux via a direct knock-on mechanism.

Inwardly rectifying potassium (Kir) channels play a key role in controlling membrane potentials in excitable and unexcitable cells, thereby regulating a plethora of physiological processes. G-protein–gated Kir channels control heart rate and neuronal excitability via small hyperpolarizing outward K⁺ currents near the resting membrane potential. Despite recent breakthroughs in x-ray crystallography and cryo-EM, the gating and conduction mechanisms of these channels are poorly understood. MD simulations have provided unprecedented details concerning the gating and conduction mechanisms of voltage-gated K⁺ and Na⁺ channels. Here, we use multi-microsecond–timescale MD simulations based on the crystal structures of GIRK2 (Kir3.2) bound to phosphatidylinositol-4,5-bisphosphate to provide detailed insights into the channel’s gating dynamics, including insights into the behavior of the G-loop gate. The simulations also elucidate the elementary steps that underlie the movement of K⁺ ions through an inward-rectifier K⁺ channel under an applied electric field. Our simulations suggest that K⁺ permeation might occur via direct knock-on, similar to the mechanism recently shown for K_v channels.

Water and hydrophobic gates in ion channels and nanopores

Simulations of water behaviour have been used to probe hydrophobic gates in BEST1 and TMEM175, which can reveal important design principles for the engineering of gates in novel biomimetic nanopores.

Ion channel proteins form nanopores in biological membranes which allow the passage of ions and water molecules. Hydrophobic constrictions in such pores can form gates, i.e. energetic barriers to water and ion permeation. Molecular dynamics simulations of water in ion channels may be used to assess whether a hydrophobic gate is closed (i.e. impermeable to ions) or open. If there is an energetic barrier to water permeation then it is likely that a gate will also be impermeable to ions. Simulations of water behaviour have been used to probe hydrophobic gates in two recently reported ion channel structures: BEST1 and TMEM175. In each of these channels a narrow region is formed by three consecutive rings of hydrophobic sidechains and in both cases such analysis demonstrates that the crystal structures correspond to a closed state of the channel. In silico mutations of BEST1 have also been used to explore the effect of changes in the hydrophobicity of the gating constriction, demonstrating that substitution of hydrophobic sidechains with more polar sidechains results in an open gate which allows water permeation. A possible open state of the TMEM175 channel was modelled by the in silico expansion of the hydrophobic gate resulting in the wetting of the pore and free permeation of potassium ions through the channel. Finally, a preliminary study suggests that a hydrophobic gate motif can be transplanted in silico from the BEST1 channel into a simple β-barrel pore template. Overall, these results suggest that simulations of the behaviour of water in hydrophobic gates can reveal important design principles for the engineering of gates in novel biomimetic nanopores.

A heuristic derived from analysis of the ion channel structural proteome permits the rapid identification of hydrophobic gates

Ion channels are nanoscale protein pores in cell membranes. An exponentially increasing number of structures for channels means that computational methods for predicting their functional state are needed. Hydrophobic gates in ion channels result in local dewetting of pores, which functionally closes them to water and ion permeation. We use simulations of water behavior within nearly 200 different ion channel structures to explore how the radius and hydrophobicity of pores determine their hydration vs. dewetting behavior. Machine learning-assisted analysis of these simulations allowed us to propose a simple model for this relationship and present an easy method for rapidly predicting the functional state of new channel structures as they emerge.

Ion channel proteins control ionic flux across biological membranes through conformational changes in their transmembrane pores. An exponentially increasing number of channel structures captured in different conformational states are now being determined; however, these newly resolved structures are commonly classified as either open or closed based solely on the physical dimensions of their pore, and it is now known that more accurate annotation of their conductive state requires additional assessment of the effect of pore hydrophobicity. A narrow hydrophobic gate region may disfavor liquid-phase water, leading to local dewetting, which will form an energetic barrier to water and ion permeation without steric occlusion of the pore. Here we quantify the combined influence of radius and hydrophobicity on pore dewetting by applying molecular dynamics simulations and machine learning to nearly 200 ion channel structures. This allows us to propose a simple simulation-free heuristic model that rapidly and accurately predicts the presence of hydrophobic gates. This not only enables the functional annotation of new channel structures as soon as they are determined, but also may facilitate the design of novel nanopores controlled by hydrophobic gates.

CHAP: A Versatile Tool for the Structural and Functional Annotation of Ion Channel Pores

The control of ion channel permeation requires the modulation of energetic barriers or “gates” within their pores. However, such barriers are often simply identified from the physical dimensions of the pore. Such approaches have worked well in the past, but there is now evidence that the unusual behavior of water within narrow hydrophobic pores can produce an energetic barrier to permeation without requiring steric occlusion of the pathway. Many different ion channels have now been shown to exploit “hydrophobic gating” to regulate ion flow, and it is clear that new tools are required for more accurate functional annotation of the increasing number of ion channel structures becoming available. We have previously shown how molecular dynamics simulations of water can be used as a proxy to predict hydrophobic gates, and we now present a new and highly versatile computational tool, the Channel Annotation Package (CHAP) that implements this methodology.

•

Some ion channels exhibit hydrophobic gating via dewetting of the central pore.

•

This cannot be predicted from their pore radius alone.

•

The hydrophobicity of a pore contributes to this dewetting effect.

•

We have developed a new tool (CHAP) that combines these measurements.

•

CHAP now enables the rapid functional annotation of ion channel structures.

Automated cryo-EM structure refinement using correlation-driven molecular dynamics

We present a correlation-driven molecular dynamics (CDMD) method for automated refinement of atomistic models into cryo-electron microscopy (cryo-EM) maps at resolutions ranging from near-atomic to subnanometer. It utilizes a chemically accurate force field and thermodynamic sampling to improve the real-space correlation between the modeled structure and the cryo-EM map. Our framework employs a gradual increase in resolution and map-model agreement as well as simulated annealing, and allows fully automated refinement without manual intervention or any additional rotamer- and backbone-specific restraints. Using multiple challenging systems covering a wide range of map resolutions, system sizes, starting model geometries and distances from the target state, we assess the quality of generated models in terms of both model accuracy and potential of overfitting. To provide an objective comparison, we apply several well-established methods across all examples and demonstrate that CDMD performs best in most cases.

Molecular mechanisms of gating in the calcium-activated chloride channel bestrophin

Bestrophin (BEST1-4) ligand-gated chloride (Cl^-) channels are activated by calcium (Ca²⁺). Mutation of BEST1 causes retinal disease. Partly because bestrophin channels have no sequence or structural similarity to other ion channels, the molecular mechanisms underlying gating are unknown. Here, we present a series of cryo-electron microscopy structures of chicken BEST1, determined at 3.1 Å resolution or better, that represent the channel’s principal gating states. Unlike other channels, opening of the pore is due to the repositioning of tethered pore-lining helices within a surrounding protein shell that dramatically widens a neck of the pore through a concertina of amino acid rearrangements. The neck serves as both the activation and the inactivation gate. Ca²⁺ binding instigates opening of the neck through allosteric means whereas inactivation peptide binding induces closing. An aperture within the otherwise wide pore controls anion permeability. The studies define a new molecular paradigm for gating among ligand-gated ion channels.

An allosteric mechanism of inactivation in the calcium-dependent chloride channel BEST1

Bestrophin proteins are calcium (Ca²⁺)-activated chloride channels. Mutations in bestrophin 1 (BEST1) cause macular degenerative disorders. Whole-cell recordings show that ionic currents through BEST1 run down over time, but it is unclear whether this behavior is intrinsic to the channel or the result of cellular factors. Here, using planar lipid bilayer recordings of purified BEST1, we show that current rundown is an inherent property of the channel that can now be characterized as inactivation. Inactivation depends on the cytosolic concentration of Ca²⁺, such that higher concentrations stimulate inactivation. We identify a C-terminal inactivation peptide that is necessary for inactivation and dynamically interacts with a receptor site on the channel. Alterations of the peptide or its receptor dramatically reduce inactivation. Unlike inactivation peptides of voltage-gated channels that bind within the ion pore, the receptor for the inactivation peptide is on the cytosolic surface of the channel and separated from the pore. Biochemical, structural, and electrophysiological analyses indicate that binding of the peptide to its receptor promotes inactivation, whereas dissociation prevents it. Using additional mutational studies we find that the "neck" constriction of the pore, which we have previously shown to act as the Ca²⁺-dependent activation gate, also functions as the inactivation gate. Our results indicate that unlike a ball-and-chain inactivation mechanism involving physical occlusion of the pore, inactivation in BEST1 occurs through an allosteric mechanism wherein binding of a peptide to a surface-exposed receptor controls a structurally distant gate.