Water and hydrophobic gates in ion channels and nanopores

Mechanism and therapeutic targets of the involvement of a novel lysosomal proton channel TMEM175 in Parkinson's disease

Parkinson's disease (PD), recognized as the second most prevalent neurodegenerative disease in the aging population, presents a significant challenge due to the current lack of effective treatment methods to mitigate its progression. Many pathogenesis of PD are related to lysosomal dysfunction. Moreover, extensive genetic studies have shown a significant correlation between the lysosomal membrane protein TMEM175 and the risk of developing PD. Building on this discovery, TMEM175 has been identified as a novel potassium ion channel. Intriguingly, further investigations have found that potassium ion channels gradually close and transform into hydrion "excretion" channels in the microenvironment of lysosomes. This finding was further substantiated by studies on TMEM175 knockout mice, which exhibited pronounced motor dysfunction in pole climbing and suspension tests, alongside a notable reduction in dopamine neurons within the substantia nigra compacta. Despite these advancements, the current research landscape is not without its controversies. In light of this, the present review endeavors to methodically examine and consolidate a vast array of recent literature on TMEM175. This comprehensive analysis spans from the foundational research on the structure and function of TMEM175 to expansive population genetics studies and mechanism research utilizing cellular and animal models.A thorough understanding of the structure and function of TMEM175, coupled with insights into the intricate mechanisms underpinning lysosomal dysfunction in PD dopaminergic neurons, is imperative. Such knowledge is crucial for pinpointing precise intervention targets, thereby paving the way for novel therapeutic strategies that could potentially alter the neurodegenerative trajectory of PD.

Dynamic molecular portraits of ion-conducting pores characterize functional states of TRPV channels

Structural biology is solving an ever-increasing number of snapshots of ion channel conformational ensembles. Deciphering ion channel mechanisms, however, requires understanding the ensemble dynamics beyond the static structures. Here, we present a molecular modeling-based approach characterizing the ion channel structural intermediates, or their "dynamic molecular portraits", by assessing water and ion conductivity along with the detailed evaluation of pore hydrophobicity and residue packing. We illustrate the power of this approach by analyzing structures of few vanilloid-subfamily transient receptor potential (TRPV) channels. Based on the pore architecture, there are three major states that are common for TRPVs, which we call α-closed, π-closed, and π-open. We show that the pore hydrophobicity and residue packing for the open state is most favorable for the pore conductance. On the contrary, the α-closed state is the most hydrophobic and always non-conducting. Our approach can also be used for structural and functional classification of ion channels.

Ionic Liquid Gating Induces Anomalous Permeation through Membrane Channel Proteins

Membrane channel proteins (MCPs) play key roles in matter transport through cell membranes and act as major targets for vaccines and drugs. For emerging ionic liquid (IL) drugs, a rational understanding of how ILs affect the structure and transport function of MCP is crucial to their design. In this work, GPU-accelerated microsecond-long molecular dynamics simulations were employed to investigate the modulating mechanism of ILs on MCP. Interestingly, ILs prefer to insert into the lipid bilayer and channel of aquaporin-2 (AQP2) but adsorb on the entrance of voltage-gated sodium channels (Nav). Molecular trajectory and free energy analysis reflect that ILs have a minimal impact on the structure of MCPs but significantly influence MCP functions. It demonstrates that ILs can decrease the overall energy barrier for water through AQP2 by 1.88 kcal/mol, whereas that for Na+ through Nav is increased by 1.70 kcal/mol. Consequently, the permeation rates of water and Na+ can be enhanced and reduced by at least 1 order of magnitude, respectively. Furthermore, an abnormal IL gating mechanism was proposed by combining the hydrophobic nature of MCP and confined water/ion coordination effects. More importantly, we performed experiments to confirm the influence of ILs on AQP2 in human cells and found that treatment with ILs significantly accelerated the changes in cell volume in response to altered external osmotic pressure. Overall, these quantitative results will not only deepen the understanding of IL-cell interactions but may also shed light on the rational design of drugs and disease diagnosis.

Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants

SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at μ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.

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 chloride channel (hBest1) provide an opportunity to evaluate ion interactions predicted by molecular dynamics (MD) simulations against experimental observations. We implement the fully polarizable forcefield AMOEBA in MD simulations of open and partially-open states of the hBest1. The AMOEBA forcefield models multipole moments up to the quadrupole; therefore, it captures induced dipole and anion- π interactions. By including polarization we demonstrate the key role that aromatic residues play in ion permeation and the functional advantages of pore asymmetry within the highly conserved hydrophobic neck of the pore. We establish that these only arise when electronic polarization is included in the molecular models. We also show that Cl - permeation in this region can be achieved through hydrophobic solvation concomitant with partial ion dehydration, which is compensated for by the formation of contacts with the edge of the phenylalanine ring. Furthermore, we demonstrate how polarizable simulations can help determine the identity of ion-like densities within high-resolution cryo-EM structures. Crucially, neglecting polarization in simulation of these systems results in the localization of 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.

Statement of Significance

Ion channels are nanoscale protein pores that enable the selective passage of charged ions across cell membranes. Understanding the underlying mechanisms for selective anion permeation through such pores remains a challenge. To simulate their behavior efficiently in silico , fixed charge models are typically employed. However, this approach is insufficient for the study of anions. Here, we use simulations with explicit treatment of electrostatics to investigate the interactions of chloride ions in the human bestrophin 1 channel. We find that electronic polarization tunes the state of the channel and affects the interactions of chloride ions thereby revealing a mechanism for permeation. Furthermore, these simulations can be used to resolve experimental ambiguity in ion-like densities from cryo-EM structures.

What keeps nanopores boiling

The liquid-to-vapor transition can occur under unexpected conditions in nanopores, opening the door to fundamental questions and new technologies. The physics of boiling in confinement is progressively introduced, starting from classical nucleation theory, passing through nanoscale effects, and terminating with the material and external parameters that affect the boiling conditions. The relevance of boiling in specific nanoconfined systems is discussed, focusing on heterogeneous lyophobic systems, chromatographic columns, and ion channels. The current level of control of boiling in nanopores enabled by microporous materials such as metal organic frameworks and biological nanopores paves the way to thrilling theoretical challenges and to new technological opportunities in the fields of energy, neuromorphic computing, and sensing.

Transmembrane Protein 175, a Lysosomal Ion Channel Related to Parkinson's Disease

Lysosomes are membrane-bound organelles with an acidic lumen and are traditionally characterized as a recycling center in cells. Lysosomal ion channels are integral membrane proteins that form pores in lysosomal membranes and allow the influx and efflux of essential ions. Transmembrane protein 175 (TMEM175) is a unique lysosomal potassium channel that shares little sequence similarity with other potassium channels. It is found in bacteria, archaea, and animals. The prokaryotic TMEM175 consists of one six-transmembrane domain that adopts a tetrameric architecture, while the mammalian TMEM175 is comprised of two six-transmembrane domains that function as a dimer in lysosomal membranes. Previous studies have demonstrated that the lysosomal K⁺ conductance mediated by TMEM175 is critical for setting membrane potential, maintaining pH stability, and regulating lysosome-autophagosome fusion. AKT and B-cell lymphoma 2 regulate TMEM175's channel activity through direct binding. Two recent studies reported that the human TMEM175 is also a proton-selective channel under normal lysosomal pH (4.5-5.5) as the K⁺ permeation dramatically decreased at low pH while the H⁺ current through TMEM175 greatly increased. Genome-wide association studies and functional studies in mouse models have established that TMEM175 is implicated in the pathogenesis of Parkinson's disease, which sparks more research interests in this lysosomal channel.

Structural Mechanism of Ionic Conductivity of the TRPV1 Channel

The so-called "hydrophobic gating" is widely discussed as a putative mechanism to control water and ion conduction via ion channels. This effect can occur in narrow areas of the channels pore lined by non-polar residues. In the closed state of the channel, such regions may spontaneously transit to a dehydrated state to block water and ions transport without full pore occlusion. In the open state, the hydrophobic gate is wide enough to provide sustainable hydration and conduction. Apparently, the transport through the open hydrophobic gate may by facilitated by some polar residues that assist polar/charged substances to overcome the energy barrier created by nonpolar environment. In this work, we investigated the behavior of Na⁺ ions and their hydration shells in the open pore of the rat TRPV1 ion channel by molecular dynamics simulations. We show that polar protein groups coordinate water molecules in such a way as to restore the hydration shell of ions in the hydrophobic gate that ensures ion transport through the gate in a fully hydrated state.

Molecular dynamics study of Cl- permeation through cystic fibrosis transmembrane conductance regulator (CFTR)

The recent elucidation of atomistic structures of Cl^- channel CFTR provides opportunities for understanding the molecular basis of cystic fibrosis. Despite having been activated through phosphorylation and provided with ATP ligands, several near-atomistic cryo-EM structures of CFTR are in a closed state, as inferred from the lack of a continuous passage through a hydrophobic bottleneck region located in the extracellular portion of the pore. Here, we present repeated, microsecond-long molecular dynamics simulations of human CFTR solvated in a lipid bilayer and aqueous NaCl. At equilibrium, Cl^- ions enter the channel through a lateral intracellular portal and bind to two distinct cationic sites inside the channel pore but do not traverse the narrow, de-wetted bottleneck. Simulations conducted in the presence of a strong hyperpolarizing electric field led to spontaneous Cl^- translocation events through the bottleneck region of the channel, suggesting that the protein relaxed to a functionally open state. Conformational changes of small magnitude involving transmembrane helices 1 and 6 preceded ion permeation through diverging exit routes at the extracellular end of the pore. The pore bottleneck undergoes wetting prior to Cl^- translocation, suggesting that it acts as a hydrophobic gate. Although permeating Cl^- ions remain mostly hydrated, partial dehydration occurs at the binding sites and in the bottleneck. The observed Cl^- pathway is largely consistent with the loci of mutations that alter channel conductance, anion binding, and ion selectivity, supporting the model of the open state of CFTR obtained in the present study.

Quantum Mechanical Aspects in the Pathophysiology of Neuropathic Pain

Neuropathic pain is a challenging complaint for patients and clinicians since there are no effective agents available to get satisfactory outcomes even though the pharmacological agents target reasonable pathophysiological mechanisms. This may indicate that other aspects in these mechanisms should be unveiled to comprehend the pathogenesis of neuropathic pain and thus find more effective treatments. Therefore, in the present study, several mechanisms are chosen to be reconsidered in the pathophysiology of neuropathic pain from a quantum mechanical perspective. The mathematical model of the ions quantum tunneling model is used to provide quantum aspects in the pathophysiology of neuropathic pain. Three major pathophysiological mechanisms are revisited in the context of the quantum tunneling model. These include: (1) the depolarized membrane potential of neurons; (2) the cross-talk or the ephaptic coupling between the neurons; and (3) the spontaneous neuronal activity and the emergence of ectopic action potentials. We will show mathematically that the quantum tunneling model can predict the occurrence of neuronal membrane depolarization attributed to the quantum tunneling current of sodium ions. Moreover, the probability of inducing an ectopic action potential in the axons of neurons will be calculated and will be shown to be significant and influential. These ectopic action potentials are generated due to the formation of quantum synapses which are assumed to be the mechanism behind the ephaptic transmission. Furthermore, the spontaneous neuronal activity and the emergence of ectopic action potentials independently from any adjacent stimulated neurons are predicted to occur according to the quantum tunneling model. All these quantum mechanical aspects contribute to the overall hyperexcitability of the neurons and to the pathogenesis of neuropathic pain. Additionally, providing a new perspective in the pathophysiology of neuropathic pain may improve our understanding of how the neuropathic pain is generated and maintained and may offer new effective agents that can improve the overall clinical outcomes of the patients.

Opening of glutamate receptor channel to subconductance levels

Ionotropic glutamate receptors (iGluRs) are tetrameric ligand-gated ion channels that open their pores in response to binding of the agonist glutamate1-3. An ionic current through a single iGluR channel shows up to four discrete conductance levels (O1-O4)4-6. Higher conductance levels have been associated with an increased number of agonist molecules bound to four individual ligand-binding domains (LBDs)6-10. Here we determine structures of a synaptic complex of AMPA-subtype iGluR and the auxiliary subunit γ2 in non-desensitizing conditions with various occupancy of the LBDs by glutamate. We show that glutamate binds to LBDs of subunits B and D only after it is already bound to at least the same number of LBDs that belong to subunits A and C. Our structures combined with single-channel recordings, molecular dynamics simulations and machine-learning analysis suggest that channel opening requires agonist binding to at least two LBDs. Conversely, agonist binding to all four LBDs does not guarantee maximal channel conductance and favours subconductance states O1 and O2, with O3 and O4 being rare and not captured structurally. The lack of subunit independence and low efficiency coupling of glutamate binding to channel opening underlie the gating of synaptic complexes to submaximal conductance levels, which provide a potential for upregulation of synaptic activity.

The Quantum Tunneling of Ions Model Can Explain the Pathophysiology of Tinnitus

Tinnitus is a well-known pathological entity in clinical practice. However, the pathophysiological mechanisms behind tinnitus seem to be elusive and cannot provide a comprehensive understanding of its pathogenesis and clinical manifestations. Hence, in the present study, we explore the mathematical model of ions’ quantum tunneling to propose an original pathophysiological mechanism for the sensation of tinnitus. The present model focuses on two major aspects: The first aspect is the ability of ions, including sodium, potassium, and calcium, to depolarize the membrane potential of inner hair cells and the neurons of the auditory pathway. This membrane depolarization is induced via the quantum tunneling of ions through closed voltage-gated channels. The state of membrane depolarization can be a state of hyper-excitability or hypo-excitability, depending on the degree of depolarization. Both of these states aid in understanding the pathophysiology of tinnitus. The second aspect is the quantum tunneling signals between the demyelinated neurons of the auditory pathway. These signals are mediated via the quantum tunneling of potassium ions, which exit to the extracellular fluid during an action potential event. These quantum signals can be viewed as a “quantum synapse” between neurons. The formation of quantum synapses results in hyper-excitability among the demyelinated neurons of the auditory pathway. Both of these aspects augment and amplify the electrical signals in the auditory pathway and result in a loss of the spatiotemporal fidelity of sound signals going to the brain centers. The brain interprets this hyper-excitability and loss of spatiotemporal fidelity as tinnitus. Herein, we show mathematically that the quantum tunneling of ions can depolarize the membrane potential of the inner hair cells and neurons of the auditory pathway. Moreover, we calculate the probability of action potential induction in the neurons of the auditory pathway generated by the quantum tunneling signals of potassium ions.

An open state of a voltage-gated sodium channel involving a π-helix and conserved pore-facing asparagine

Voltage-gated sodium (Nav) channels play critical roles in propagating action potentials and otherwise manipulating ionic gradients in excitable cells. These channels open in response to membrane depolarization, selectively permeating sodium ions until rapidly inactivating. Structural characterization of the gating cycle in this channel family has proved challenging, particularly due to the transient nature of the open state. A structure from the bacterium Magnetococcus marinus Nav (NavMs) was initially proposed to be open, based on its pore diameter and voltage-sensor conformation. However, the functional annotation of this model, and the structural details of the open state, remain disputed. In this work, we used molecular modeling and simulations to test possible open-state models of NavMs. The full-length experimental structure, termed here the α-model, was consistently dehydrated at the activation gate, indicating an inability to conduct ions. Based on a spontaneous transition observed in extended simulations, and sequence/structure comparison to other Nav channels, we built an alternative π-model featuring a helix transition and the rotation of a conserved asparagine residue into the activation gate. Pore hydration, ion permeation, and state-dependent drug binding in this model were consistent with an open functional state. This work thus offers both a functional annotation of the full-length NavMs structure and a detailed model for a stable Nav open state, with potential conservation in diverse ion-channel families.

Water Nanoconfined in a Hydrophobic Pore: Molecular Dynamics Simulations of Transmembrane Protein 175 and the Influence of Water Models

Water molecules within biological ion channels are in a nanoconfined environment and therefore exhibit behaviors which differ from that of bulk water. Here, we investigate the phenomenon of hydrophobic gating, the process by which a nanopore may spontaneously dewet to form a "vapor lock" if the pore is sufficiently hydrophobic and/or narrow. This occurs without steric occlusion of the pore. Using molecular dynamics simulations with both rigid fixed-charge and polarizable (AMOEBA) force fields, we investigate this wetting/dewetting behavior in the transmembrane protein 175 ion channel. We examine how a range of rigid fixed-charge and polarizable water models affect wetting/dewetting in both the wild-type structure and in mutants chosen to cover a range of nanopore radii and pore-lining hydrophobicities. Crucially, we find that the rigid fixed-charge water models lead to similar wetting/dewetting behaviors, but that the polarizable water model resulted in an increased wettability of the hydrophobic gating region of the pore. This has significant implications for molecular simulations of nanoconfined water, as it implies that polarizability may need to be included if we are to gain detailed mechanistic insights into wetting/dewetting processes. These findings are of importance for the design of functionalized biomimetic nanopores (e.g., sensing or desalination) as well as for furthering our understanding of the mechanistic processes underlying biological ion channel function.

Ion permeation controlled by hydrophobic residues and proton binding in the proton-activated chloride channel

Recently identified proton-activated chloride channel (PAC) contains two transmembrane helices (S1–S2) and is involved in lysosome function, hypoxia adaption, stroke, and carcinogenesis. Although a PAC structure was recently resolved, its gating and activation mechanisms remained largely unknown. By the two-electrode voltage clamp electrophysiology in Xenopus oocytes, we found that the hydrophobicity of site 304 at fenestrations, but not that of neighbor sites, is important for maintaining PAC at a closed state at pH 7.5. When activated at acidic pH, PAC activity significantly increased with the hydrophilicity of site 307 within S2, but not with that of neighbor sites, suggesting that 307 acts as an activation gate. We identified six conserved protonatable residues critical for proton-induced activation, consistent with structural studies. Our study depicted a scheme in which proton binding induces conformational changes from the W304-controlled closed state at fenestrations to an activated state controlled by activation gate I307 in helix S2.

•

The hydrophobicity of site 304 is critical for maintaining PAC at a closed state

•

The function of activated PAC is modulated by the hydrophilicity of site 307

•

Six protonatable amino acids are involved in proton-induced PAC activation

•

H⁺ binding seem to change PAC from W304-controlled closed to I307-gated open state

Water-Induced Self-Assembly of Amphiphilic Discotic Molecules for Adaptive Artificial Water Channels

Inspired by the induced-fit mechanism in nature, we developed the process of water-induced self-assembly (WISA) to make water an active substrate that regulates the self-assembly and function of amphiphilic discotic molecules (ADMs). The ADM is an isotropic liquid that self-assembles only when in contact with water. Characterization results indicate that water fits into the hydrophilic core of the ADMs and induces the formation of a hexagonal columnar phase (Col_h), where each column contains a hydrated artificial water channel (AWC). The hydrated AWCs are adaptive rather than static; the dynamic incorporation/removal of water results in the reversible assembly/disassembly of the adaptive AWCs (aAWCs). Furthermore, its dynamic characteristics can enable water to act as an orientation-directional guest molecule that controls the growth direction of the aAWCs. Well-aligned aAWC arrays that showed the ability of water transport were obtained via a "directional WISA" method. In WISA, water thus governs the supramolecular chemistry and function of synthetic molecules as it does with natural materials. By making water an active component in adaptive chemistry and enabling host molecules to dynamically interact with water, this adaptive aquatic material may motivate the development of synthetic molecules further toward biomaterials.

Transmembrane Self-Assembled Cyclic Peptide Nanotubes Based on α-Residues and Cyclic δ-Amino Acids: A Computational Study

Self-assembling cyclic peptide nanotubes have been shown to function as synthetic, integral transmembrane channels. The combination of natural and nonnatural aminoacids in the sequence of cyclic peptides enables the control not only of their outer surface but also of the inner cavity behavior and properties, affecting, for instance, their permeability to different molecules including water and ions. Here, a thorough computational study on a new class of self-assembling peptide motifs, in which δ-aminocycloalkanecarboxylic acids are alternated with natural α-amino acids, is presented. The presence of synthetic δ-residues creates hydrophobic regions in these α,δ-SCPNs, which makes them especially attractive for their potential implementation in the design of new drug or diagnostic agent carrier systems. Using molecular dynamics simulations, the behavior of water molecules, different ions (Li+, Na+, K+, Cs+, and Ca2+), and their correspondent counter Cl- anions is extensively investigated in the nanoconfined environment. The structure and dynamics are mutually combined in a diving immersion inside these transmembrane channels to discover a fascinating submarine nanoworld where star-shaped water channels guide the passage of cations and anions therethrough.

Adaptable and Multifunctional Ion-Conducting Aquaporins

Aquaporins function as water and neutral solute channels, signaling hubs, disease virulence factors, and metabolon components. We consider plant aquaporins that transport ions compared to some animal counterparts. These are candidates for important, as yet unidentified, cation and anion channels in plasma, tonoplast, and symbiotic membranes. For those individual isoforms that transport ions, water, and gases, the permeability spans 12 orders of magnitude. This requires tight regulation of selectivity via protein interactions and posttranslational modifications. A phosphorylation-dependent switch between ion and water permeation in AtPIP2;1 might be explained by coupling between the gates of the four monomer water channels and the central pore of the tetramer. We consider the potential for coupling between ion and water fluxes that could form the basis of an electroosmotic transducer. A grand challenge in understanding the roles of ion transporting aquaporins is their multifunctional modes that are dependent on location, stress, time, and development.

Hydrophobic Gating and 1/f Noise of the Anthrax Toxin Channel

"Pink" or 1/f noise is a natural phenomenon omnipresent in physics, economics, astrophysics, biology, and even music and languages. In electrophysiology, the stochastic activity of a number of biological ion channels and artificial nanopores could be characterized by current noise with a 1/f power spectral density. In the anthrax toxin channel (PA₆₃), it appears as fast voltage-independent current interruptions between conducting and nonconducting states. This behavior hampers potential development of PA₆₃ as an ion-channel biosensor. On the bright side, the PA₆₃ flickering represents a mesmerizing phenomenon to investigate. Notably, similar 1/f fluctuations are observed in the channel-forming components of clostridial binary C2 and iota toxins, which share functional and structural similarities with the anthrax toxin channel. Similar to PA₆₃, they are evolved to translocate the enzymatic components of the toxins into the cytosol. Here, using high-resolution single-channel lipid bilayer experiments and all-atom molecular dynamic simulations, we suggest that the 1/f noise in PA₆₃ occurs as a result of "hydrophobic gating" at the ϕ-clamp region, the phenomenon earlier observed in several water-filled channels "fastened" inside by the hydrophobic belts. The ϕ-clamp is a narrow "hydrophobic ring" in the PA₆₃ lumen formed by seven or eight phenylalanine residues at position 427, conserved in the C2 and iota toxin channels, which catalyzes protein translocation. Notably, the 1/f noise remains undetected in the F427A PA₆₃ mutant. This finding can elucidate the functional purpose of 1/f noise and its possible role in the transport of the enzymatic components of binary toxins.

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.

Gas-Induced Drying of Nanopores

Here, we investigate the role of a dilute hydrophobic gas on the phase behavior of water confined in hydrophobic nanopores. Molecular dynamics showed that gas atoms are hydrophobically attracted within the pores, where even a single particle is able to induce spontaneous drying of the whole pore. The drying process is rationalized in terms of its free-energy landscape, revealing that the penetration of a gas atom is able to suppress the drying free-energy barriers of hydrophobic pores. Results provide insights into the role of gases on the wettability of nanopores and evidence of a possibile physical mechanism for the action of volatile anesthetics on some kinds of ion channels. Results also indicate a novel, bioinspired strategy for controlling bubble formation in nanopores for sensing and energy applications.

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.

Multiscale Simulation Reveals Passive Proton Transport Through SERCA on the Microsecond Timescale

The sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) transports two Ca²⁺ ions from the cytoplasm to the reticulum lumen at the expense of ATP hydrolysis. In addition to transporting Ca²⁺, SERCA facilitates bidirectional proton transport across the sarcoplasmic reticulum to maintain the charge balance of the transport sites and to balance the charge deficit generated by the exchange of Ca²⁺. Previous studies have shown the existence of a transient water-filled pore in SERCA that connects the Ca²⁺ binding sites with the lumen, but the capacity of this pathway to sustain passive proton transport has remained unknown. In this study, we used the multiscale reactive molecular dynamics method and free energy sampling to quantify the free energy profile and timescale of the proton transport across this pathway while also explicitly accounting for the dynamically coupled hydration changes of the pore. We find that proton transport from the central binding site to the lumen has a microsecond timescale, revealing a novel passive cytoplasm-to-lumen proton flow beside the well-known inverse proton countertransport occurring in active Ca²⁺ transport. We propose that this proton transport mechanism is operational and serves as a functional conduit for passive proton transport across the sarcoplasmic reticulum.

The mechanism and energetics of a ligand-controlled hydrophobic gate in a mammalian two pore channel

In the last decade two-pore intracellular channels (TPCs) attracted the interest of researchers, still some key questions remain open. Their importance for vacuolar (plants) and endo-lysosomal (animals) function highlights them as a very attractive system to study, both theoretically and experimentally. Indicated as key players in the trafficking of the cell, today they are considered a new potential target for avoiding virus infections, including those from coronaviruses. A particular boost for theoretical examinations has been made with recent high-resolution X-ray and cryo-EM structures. These findings have opened the way for efficient and precise computational studies at the atomistic level. Here we report a set of multiscale-calculations performed on the mTPC1, a ligand- and voltage-gated sodium selective channel. The molecular dynamics and enhanced molecular dynamics simulations were used for a thorough analysis of the mammalian TPC1 behaviour in the presence and absence of the ligand molecule, with a special accent on the supposed bottleneck, the hydrophobic gate. Moreover, from the reconstructed free energy obtained from enhanced simulations, we have calculated the macroscopic conductance of sodium ions through the mTPC1, which we compared with measured single-channel conductance values. The hydrophobic gate works as a steric barrier and the key parameters are its flexibility and the dimension of the sodium first hydration shell.

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.

Conformational dynamics and interfacial interactions of peptide-appended pillar[5]arene water channels in biomimetic membranes

Peptide appended pillar[5]arene (PAP) is an artificial water channel resembling biological water channel proteins, which has shown a significant potential for designing bioinspired water purification systems. Given that PAP channels need to be incorporated at a high density in membrane matrices, it is critical to examine the role of channel-channel and channel-membrane interactions in governing the structural and functional characteristics of channels. To resolve the atomic-scale details of these interactions, we have carried out atomistic molecular dynamics (MD) simulations of multiple PAP channels inserted in a lipid or a block-copolymer (BCP) membrane matrix. Classical MD simulations on a sub-microsecond timescale showed clustering of channels only in the lipid membrane, but enhanced sampling MD simulations showed thermodynamically-favorable dimerized states of channels in both lipid and BCP membranes. The dimerized configurations of channels, with an extensive buried surface area, were stabilized via interactions between the aromatic groups in the peptide arms of neighboring channels. The conformational metrics characterizing the orientational and structural changes in channels revealed a higher flexibility in the lipid membrane as opposed to the BCP membrane although hydrogen bonds between the channel and the membrane molecules were not a major contributor to the stability of channels in the BCP membrane. We also found that the channels undergo wetting/dewetting transitions in both lipid and BCP membranes with a marginally higher probability of undergoing a dewetting transition in the BCP membrane. Collectively, these results highlight the role of channel dynamics in governing channel-channel and channel-membrane interfacial interactions, and provide atomic-scale insights needed to design stable and functional biomimetic membranes for efficient separations.

Functionalization of the TMEM175 p.M393T variant as a risk factor for Parkinson disease

Multiple genome-wide association studies (GWAS) in Parkinson disease (PD) have identified a signal at chromosome 4p16.3; however, the causal variant has not been established for this locus. Deep investigation of the region resulted in one identified variant, the rs34311866 missense SNP (p.M393T) in TMEM175, which is 20 orders of magnitude more significant than any other SNP in the region. Because TMEM175 is a lysosomal gene that has been shown to influence α-synuclein phosphorylation and autophagy, the p.M393T variant is an attractive candidate, and we have examined its effect on TMEM175 protein and PD-related biology. After knocking down each of the genes located under the GWAS peak via multiple shRNAs, only TMEM175 was found to consistently influence accumulation of phosphorylated α-synuclein (p-α-syn). Examination of the p.M393T variant showed effects on TMEM175 function that were intermediate between the wild-type (WT) and knockout phenotypes, with reduced regulation of lysosomal pH in response to starvation and minor changes in clearance of autophagy substrates, reduced lysosomal localization, and increased accumulation of p-α-syn. Finally, overexpression of WT TMEM175 protein reduced p-α-syn, while overexpression of the p.M393T variant resulted in no change in α-synuclein phosphorylation. These results suggest that the main signal in the chromosome 4p16.3 PD risk locus is driven by the TMEM175 p.M393T variant. Modulation of TMEM175 may impact α-synuclein biology and therefore may be a rational therapeutic strategy for PD.

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.

Confined Dynamics of Water in Transmembrane Pore of TRPV1 Ion Channel

Solvation effects play a key role in chemical and biological processes. The microscopic properties of water near molecular surfaces are radically different from those in the bulk. Furthermore, the behavior of water in confined volumes of a nanometer scale, including transmembrane pores of ion channels, is especially nontrivial. Knowledge at the molecular level of structural and dynamic parameters of water in such systems is necessary to understand the mechanisms of ion channels functioning. In this work, the results of molecular dynamics (MD) simulations of water in the pore and selectivity filter domains of TRPV1 (Transient Receptor Potential Vanilloid type 1) membrane channel are considered. These domains represent nanoscale volumes with strongly amphiphilic walls, where physical behavior of water radically differs from that of free hydration (e.g., at protein interfaces) or in the bulk. Inside the pore and filter domains, water reveals a very heterogeneous spatial distribution and unusual dynamics: It forms compact areas localized near polar groups of particular residues. Residence time of water molecules in such areas is at least 1.5 to 3 times larger than that observed for similar groups at the protein surface. Presumably, these water “blobs” play an important role in the functional activity of TRPV1. In particular, they take part in hydration of the hydrophobic TRPV1 pore by localizing up to six waters near the so-called “lower gate” of the channel and reducing by this way the free energy barrier for ion and water transport. Although the channel is formed by four identical protein subunits, which are symmetrically packed in the initial experimental 3D structure, in the course of MD simulations, hydration of the same amino acid residues of individual subunits may differ significantly. This greatly affects the microscopic picture of the distribution of water in the channel and, potentially, the mechanism of its functioning. Therefore, reconstruction of the full picture of TRPV1 channel solvation requires thorough atomistic simulations and analysis. It is important that the naturally occurring porous volumes, like ion-conducting protein domains, reveal much more sophisticated and fine-tuned regulation of solvation than, e.g., artificially designed carbon nanotubes.

Hydrophobic gasket mutation produces gating pore currents in closed human voltage-gated proton channels

A large family of membrane proteins, voltage-gated ion channels, regulate a vast array of physiological functions in essentially all life forms. How these molecules sense membrane potential and respond by creating ionic conduction is incompletely understood. The voltage sensors of these channels contain a “hydrophobic gasket,” a ring of hydrophobic amino acids near the center of the membrane, separating internal and external aqueous solutions. Although voltage-gated proton channels, H_V1, resemble voltage-sensing domains of other channels, they differ fundamentally. On depolarization, H_V1 conducts protons, whereas other voltage sensors open a physically distinct pore. We identify Val¹⁰⁹, Phe¹⁵⁰, Val¹⁷⁷, and Val¹⁷⁸ as the hH_V1 hydrophobic gasket. Replacement with less hydrophobic amino acids accelerated channel opening and caused proton-selective leak through closed channels.

The hydrophobic gasket (HG), a ring of hydrophobic amino acids in the voltage-sensing domain of most voltage-gated ion channels, forms a constriction between internal and external aqueous vestibules. Cationic Arg or Lys side chains lining the S4 helix move through this “gating pore” when the channel opens. S4 movement may occur during gating of the human voltage-gated proton channel, hH_V1, but proton current flows through the same pore in open channels. Here, we replaced putative HG residues with less hydrophobic residues or acidic Asp. Substitution of individuals, pairs, or all 3 HG positions did not impair proton selectivity. Evidently, the HG does not act as a secondary selectivity filter. However, 2 unexpected functions of the HG in H_V1 were discovered. Mutating HG residues independently accelerated channel opening and compromised the closed state. Mutants exhibited open–closed gating, but strikingly, at negative voltages where “normal” gating produces a nonconducting closed state, the channel leaked protons. Closed-channel proton current was smaller than open-channel current and was inhibited by 10 μM Zn²⁺. Extreme hyperpolarization produced a deeper closed state through a weakly voltage-dependent transition. We functionally identify the HG as Val¹⁰⁹, Phe¹⁵⁰, Val¹⁷⁷, and Val¹⁷⁸, which play a critical and exclusive role in preventing H⁺ influx through closed channels. Molecular dynamics simulations revealed enhanced mobility of Arg²⁰⁸ in mutants exhibiting H⁺ leak. Mutation of HG residues produces gating pore currents reminiscent of several channelopathies.

A hydrophobic gate in the inner pore helix is the major determinant of inactivation in mechanosensitive Piezo channels

Piezo1 and Piezo2 belong to a family of mechanically-activated ion channels implicated in a wide range of physiological processes. Mechanical stimulation triggers Piezo channels to open, but their characteristic fast inactivation process results in rapid closure. Several disease-causing mutations in Piezo1 alter the rate of inactivation, highlighting the importance of inactivation to the normal function of this channel. However, despite the structural identification of two physical constrictions within the closed pore, the mechanism of inactivation remains unknown. Here we identify a functionally conserved inactivation gate in the pore-lining inner helix of mouse Piezo1 and Piezo2 that is distinct from the two constrictions. We show that this gate controls the majority of Piezo1 inactivation via a hydrophobic mechanism and that one of the physical constrictions acts as a secondary gate. Our results suggest that, unlike other rapidly inactivating ion channels, a hydrophobic barrier gives rise to fast inactivation in Piezo channels.

Discrete and Continuum Models for the Salt in Crowded Environments of Suspended Charged Particles

Electrostatic forces greatly affect the overall dynamics and diffusional activities of suspended charged particles in crowded environments. Accordingly, the concentration of counter- or co-ions in a fluid-''the salt"-determines the range, strength, and order of electrostatic interactions between particles. This environment fosters engineering routes for controlling directed assembly of particles at both the micro- and nanoscale. Here, we analyzed two computational modeling schemes that considered salt within suspensions of charged particles, or polyelectrolytes: discrete and continuum. Electrostatic interactions were included through a Green's function formalism, where the confined fundamental solution for Poisson's equation is resolved by the general geometry Ewald-like method. For the discrete model, the salt was considered as regularized point-charges with a specific valence and size, while concentration fields were defined for each ionic species for the continuum model. These considerations were evolved using Brownian dynamics of the suspended charged particles and the discrete salt ions, while a convection-diffusion transport equation, including the Nernst-Planck diffusion mechanism, accounted for the dynamics of the concentration fields. The salt/particle models were considered as suspensions under slit-confinement conditions for creating crowded "macro-ions", where density distributions and radial distribution functions were used to compare and differentiate computational models. Importantly, our analysis shows that disparate length scales or increased system size presented by the salt and suspended particles are best dealt with using concentration fields to model the ions. These findings were then validated by novel simulations of a semipermeable polyelectrolyte membrane, at the mesoscale, from which ionic channels emerged and enable ion conduction.

Artificial water channels: inspiration, progress, and challenges

Developing bioinspired artificial water channels may lead to the next-generation filtration membranes with ultra-high pore density and exclusive water permeability.