Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation

Ion channel selectivity through ion-modulated changes of selectivity filter pKa values

In bacterial voltage-gated sodium channels, the passage of ions through the pore is controlled by a selectivity filter (SF) composed of four glutamate residues. The mechanism of selectivity has been the subject of intense research, with suggested mechanisms based on steric effects, and ion-triggered conformational change. Here, we propose an alternative mechanism based on ion-triggered shifts in pKa values of SF glutamates. We study the NavMs channel for which the open channel structure is available. Our free-energy calculations based on molecular dynamics simulations suggest that pKa values of the four glutamates are higher in solution of K+ ions than in solution of Na+ ions. Higher pKa in the presence of K+ stems primarily from the higher population of dunked conformations of the protonated Glu sidechain, which exhibit a higher pKa shift. Since pKa values are close to the physiological pH, this results in predominant population of the fully deprotonated state of glutamates in Na+ solution, while protonated states are predominantly populated in K+ solution. Through molecular dynamics simulations we calculate that the deprotonated state is the most conductive, the singly protonated state is less conductive, and the doubly protonated state has significantly reduced conductance. Thus, we propose that a significant component of selectivity is achieved through ion-triggered shifts in the protonation state, which favors more conductive states for Na+ ions and less conductive states for K+ ions. This mechanism also suggests a strong pH dependence of selectivity, which has been experimentally observed in structurally similar NaChBac channels.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

Reconciling membrane protein simulations with experimental DEER spectroscopy data

Spectroscopy experiments are crucial to study membrane proteins for which traditional structure determination methods still prove challenging. Double electron-electron resonance (DEER) spectroscopy experiments provide protein residue-pair distance distributions that are indicative of their conformational heterogeneity. Atomistic molecular dynamics (MD) simulations are another tool that have been proven to be vital to study the structural dynamics of membrane proteins such as to identify inward-open, occluded, and outward-open conformations of transporter membrane proteins, among other partially open or closed states of the protein. Yet, studies have reported that there is no direct consensus between the distributional data from DEER experiments and MD simulations, which has challenged validation of structures obtained from long-timescale simulations and using simulations to design experiments. Current coping strategies for comparisons rely on heuristics, such as mapping the nearest matching peaks between two ensembles or biased simulations. Here we examine the differences in residue-pair distance distributions arising due to the choice of membranes around the protein and covalent modification of a pair of residues to nitroxide spin labels in DEER experiments. Through comparing MD simulations of two proteins, PepT_So and LeuT-both of which have been characterized using DEER experiments previously-we show that the proteins' dynamics are similar despite the choice of the detergent micelle as a membrane mimetic in DEER experiments. On the other hand, covalently modified residues show slight local differences in their dynamics and a huge divergence when the oxygen atom pair distances between spin labeled residues are measured rather than protein backbone distances. Given the computational expense associated with pairwise MTSSL labeled MD simulations, we examine the use of biased simulations to explore the conformational dynamics of the spin labels only to reveal that such simulations alter the underlying protein dynamics. Our study identifies the main cause for the mismatch between DEER experiments and MD simulations and will accelerate the development of potential mitigation strategies to improve the match.

Mechanism of Calcium Permeation in a Glutamate Receptor Ion Channel

The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) are neurotransmitter-activated cation channels ubiquitously expressed in vertebrate brains. The regulation of calcium flux through the channel pore by RNA-editing is linked to synaptic plasticity while excessive calcium influx poses a risk for neurodegeneration. Unfortunately, the molecular mechanisms underlying this key process are mostly unknown. Here, we investigated calcium conduction in calcium-permeable AMPAR using Molecular Dynamics (MD) simulations with recently introduced multisite force-field parameters for Ca²⁺. Our calculations are consistent with experiment and explain the distinct calcium permeability in different RNA-edited forms of GluA2. For one of the identified metal binding sites, multiscale Quantum Mechanics/Molecular Mechanics (QM/MM) simulations further validated the results from MD and revealed small but reproducible charge transfer between the metal ion and its first solvation shell. In addition, the ion occupancy derived from MD simulations independently reproduced the Ca²⁺ binding profile in an X-ray structure of an NaK channel mimicking the AMPAR selectivity filter. This integrated study comprising X-ray crystallography, multisite MD, and multiscale QM/MM simulations provides unprecedented insights into Ca²⁺ permeation mechanisms in AMPARs, and paves the way for studying other biological processes in which Ca²⁺ plays a pivotal role.

Mutation-Specific Differences in Kv7.1 (KCNQ1) and Kv11.1 (KCNH2) Channel Dysfunction and Long QT Syndrome Phenotypes

The electrocardiogram (ECG) empowered clinician scientists to measure the electrical activity of the heart noninvasively to identify arrhythmias and heart disease. Shortly after the standardization of the 12-lead ECG for the diagnosis of heart disease, several families with autosomal recessive (Jervell and Lange-Nielsen Syndrome) and dominant (Romano-Ward Syndrome) forms of long QT syndrome (LQTS) were identified. An abnormally long heart rate-corrected QT-interval was established as a biomarker for the risk of sudden cardiac death. Since then, the International LQTS Registry was established; a phenotypic scoring system to identify LQTS patients was developed; the major genes that associate with typical forms of LQTS were identified; and guidelines for the successful management of patients advanced. In this review, we discuss the molecular and cellular mechanisms for LQTS associated with missense variants in KCNQ1 (LQT1) and KCNH2 (LQT2). We move beyond the "benign" to a "pathogenic" binary classification scheme for different KCNQ1 and KCNH2 missense variants and discuss gene- and mutation-specific differences in K+ channel dysfunction, which can predispose people to distinct clinical phenotypes (e.g., concealed, pleiotropic, severe, etc.). We conclude by discussing the emerging computational structural modeling strategies that will distinguish between dysfunctional subtypes of KCNQ1 and KCNH2 variants, with the goal of realizing a layered precision medicine approach focused on individuals.

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.

Reorganization Free Energy of Copper Proteins in Solution, in Vacuum, and on Metal Surfaces

Metalloproteins, known to efficiently transfer electronic charge in biological systems, recently found their utilization in nanobiotechnological devices where the protein is placed into direct contact with metal surfaces. The feasibility of oxidation/reduction of the protein redox sites is affected by the reorganization free energies, one of the key parameters determining the transfer rates. While their values have been measured and computed for proteins in their native environments, i.e., in aqueous solution, the reorganization free energies of dry proteins or proteins adsorbed to metal surfaces remain unknown. Here, we investigate the redox properties of blue copper protein azurin, a prototypical redox-active metalloprotein previously probed by various experimental techniques both in solution and on metal/vacuum interfaces. We used a hybrid QM/MM computational technique based on DFT to explore protein dynamics, flexibility, and corresponding reorganization free energies in aqueous solution, vacuum, and on vacuum gold interfaces. Somewhat surprisingly, the reorganization free energy only slightly decreases when azurin is dried because the loss of the hydration shell leads to larger flexibility of the protein near its redox site. At the vacuum gold surfaces, the energetics of the structure relaxation depends on the adsorption geometry, however, significant reduction of the reorganization free energy was not observed. These findings have important consequences for the charge transport mechanism in vacuum devices, showing that the free energy barriers for protein oxidation remain significant even under ultra-high vacuum conditions.

What have molecular simulations contributed to understanding of Gram-negative bacterial cell envelopes?

Bacterial cell envelopes are compositionally complex and crowded and while highly dynamic in some areas, their molecular motion is very limited, to the point of being almost static in others. Therefore, it is no real surprise that studying them at high resolution across a range of temporal and spatial scales requires a number of different techniques. Details at atomistic to molecular scales for up to tens of microseconds are now within range for molecular dynamics simulations. Here we review how such simulations have contributed to our current understanding of the cell envelopes of Gram-negative bacteria.

Molecular Insights Into Binding and Activation of the Human KCNQ2 Channel by Retigabine

Voltage-gated potassium channels of the Kv7.x family are involved in a plethora of biological processes across many tissues in animals, and their misfunctioning could lead to several pathologies ranging from diseases caused by neuronal hyperexcitability, such as epilepsy, or traumatic injuries and painful diabetic neuropathy to autoimmune disorders. Among the members of this family, the Kv7.2 channel can form hetero-tetramers together with Kv7.3, forming the so-called M-channels, which are primary regulators of intrinsic electrical properties of neurons and of their responsiveness to synaptic inputs. Here, prompted by the similarity between the M-current and that in Kv7.2 alone, we perform a computational-based characterization of this channel in its different conformational states and in complex with the modulator retigabine. After validation of the structural models of the channel by comparison with experimental data, we investigate the effect of retigabine binding on the two extreme states of Kv7.2 (resting-closed and activated-open). Our results suggest that binding, so far structurally characterized only in the intermediate activated-closed state, is possible also in the other two functional states. Moreover, we show that some effects of this binding, such as increased flexibility of voltage sensing domains and propensity of the pore for open conformations, are virtually independent on the conformational state of the protein. Overall, our results provide new structural and dynamic insights into the functioning and the modulation of Kv7.2 and related channels.

Synthetic Macrocycle Nanopore for Potassium-Selective Transmembrane Transport

Reproducing the structure and function of biological membrane channels, synthetic nanopores have been developed for applications in membrane filtration technologies and biomolecular sensing. Stable stand-alone synthetic nanopores have been created from a variety of materials, including peptides, nucleic acids, synthetic polymers, and solid-state membranes. In contrast to biological nanopores, however, furnishing such synthetic nanopores with an atomically defined shape, including deliberate placement of each and every chemical group, remains a major challenge. Here, we introduce a chemosynthetic macromolecule-extended pillararene macrocycle (EPM)-as a chemically defined transmembrane nanopore that exhibits selective transmembrane transport. Our ionic current measurements reveal stable insertion of individual EPM nanopores into a lipid bilayer membrane and remarkable cation type-selective transport, with up to a 21-fold selectivity for potassium over sodium ions. Taken together, direct chemical synthesis offers a path to de novo design of a new class of synthetic nanopores with custom transport functionality imprinted in their atomically defined chemical structure.

Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators

Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.

Conductance selectivity of Na+ across the K+ channel via Na+ trapped in a tortuous trajectory

Ion selectivity of the potassium channel is crucial for regulating electrical activity in living cells; however, the mechanism underlying the potassium channel selectivity that favors large K⁺ over small Na⁺ remains unclear. Generally, Na⁺ is not completely excluded from permeation through potassium channels. Herein, the distinct nature of Na⁺ conduction through the prototypical KcsA potassium channel was examined. Single-channel current recordings revealed that, at a high Na⁺ concentration (200 mM), the channel was blocked by Na⁺, and this blocking was relieved at high membrane potentials, suggesting the passage of Na⁺ across the channel. At a 2,000 mM Na⁺ concentration, single-channel Na⁺ conductance was measured as one-eightieth of the K⁺ conductance, indicating that the selectivity filter allows substantial conduit of Na⁺ Molecular dynamics simulations revealed unprecedented atomic trajectories of Na⁺ permeation. In the selectivity filter having a series of carbonyl oxygen rings, a smaller Na⁺ was distributed off-center in eight carbonyl oxygen-coordinated sites as well as on-center in four carbonyl oxygen-coordinated sites. This amphipathic nature of Na⁺ coordination yielded a continuous but tortuous path along the filter. Trapping of Na⁺ in many deep free energy wells in the filter caused slow elution. Conversely, K⁺ is conducted via a straight path, and as the number of occupied K⁺ ions increased to three, the concerted conduction was accelerated dramatically, generating the conductance selectivity ratio of up to 80. The selectivity filter allows accommodation of different ion species, but the ion coordination and interactions between ions render contrast conduction rates, constituting the potassium channel conductance selectivity.

A double bilayer to study the nonequilibrium environmental response of GIRK2 in complex states

G protein-gated inwardly rectifying potassium (GIRK) channels play essential roles in electrical signaling in neurons and muscle cells. Nonequilibrium environments provide crucial driving forces behind many cellular events. Here, we apply the antiparallel alignment double bilayer model to study GIRK2 in response to the time-dependent membrane potential. Using molecular dynamics and umbrella sampling, we examined the time-dependent environmental impact on the ion conduction, energy basis, and primary motions of GIRK2 in different complex states with phosphatidylinositol-4,5-bisphosphate (PIP2) and G-protein βγ subunits (Gβγ). The antiparallel alignment double bilayer model enables us to study the transport performance in inward and outward K+ and mixed K+ and Na+. We obtained the recoverable discharge process of GIRK2 complexed with both PIP2 and Gβγ, compared with occasional conduction under PIP2-only regulation. Calculations of potential of mean force suggest different regulation by the helix bundle crossing (HBC) gate and G-loop gate regarding different complex states and under a membrane potential. In a nonequilibrium environment, distinct functional rocking motions of GIRK2 were identified under strengthened correlations between the transmembrane helices and downstream cytoplasmic domains with binding of PIP2, cations, and Gβγ. The findings suggest the potential domain motions and dynamics associated with a nonequilibrium environment and highlight the application of the antiparallel alignment double bilayer model to investigate factors in an asymmetric environment.

Direct proof of soft knock-on mechanism of ion permeation in a voltage gated sodium channel

Sodium channels selectively conduct Na⁺ ions across cellular membrane with extraordinary efficiency, which is essential for initiating action potentials. However, how Na⁺ ions permeate the ionic channels remains obscure and ambiguous. With more than 40 conductance events from microsecond molecular dynamics simulation, the soft knock-on ion permeation mediated by water molecules was observed and confirmed by the free energy profile and electrostatic potential calculation in this study. During the soft knock-on process, the change of average distance between four oxygen atoms in Glu177-Glu177 plays a very important role for the permeation of Na⁺ ion. Exploration of the ionic conductance mechanism could provide a guideline for designing ion channel targeted drug.

Determinants of conductance of a bacterial voltage-gated sodium channel

Through molecular dynamics (MD) and free energy simulations in electric fields, we examine the factors influencing conductance of bacterial voltage-gated sodium channel Na_vMs. The channel utilizes four glutamic acid residues in the selectivity filter (SF). Previously, we have shown, through constant pH and free energy calculations of pKa values, that fully deprotonated, singly protonated, and doubly protonated states are all feasible at physiological pH, depending on how many ions are bound in the SF. With 173 MD simulations of 450 or 500 ns and additional free energy simulations, we determine that the conductance is highest for the deprotonated state and decreases with each additional proton bound. We also determine that the pKa value of the four glutamic residues for the transition between deprotonated and singly protonated states is close to the physiological pH and that there is a small voltage dependence. The pKa value and conductance trends are in agreement with experimental work on bacterial Na_v channels, which show a decrease in maximal conductance with lowering of pH, with pKa in the physiological range. We examine binding sites for Na⁺ in the SF, compare with previous work, and note a dependence on starting structures. We find that narrowing of the gate backbone to values lower than the crystal structure's backbone radius reduces the conductance, whereas increasing the gate radius further does not affect the conductance. Simulations with some amount of negatively charged lipids as opposed to purely neutral lipids increases the conductance, as do simulations at higher voltages.

Molecular Insights Into the Gating Kinetics of the Cardiac hERG Channel, Illuminated by Structure and Molecular Dynamics

The rapidly activating delayed rectifier K⁺ current generated by the cardiac hERG potassium channel encoded by KCNH2 is the most important reserve current for cardiac repolarization. The unique inward rectification characteristics of the hERG channel depend on the gating regulation, which involves crucial structural domains and key single amino acid residues in the full-length hERG channel. Identifying critical molecules involved in the regulation of gating kinetics for the hERG channel requires high-resolution structures and molecular dynamics simulation models. Based on the latest progress in hERG structure and molecular dynamics simulation research, summarizing the molecules involved in the changes in the channel state helps to elucidate the unique gating characteristics of the channel and the reason for its high affinity to cardiotoxic drugs. In this review, we aim to summarize the significant advances in understanding the voltage gating regulation of the hERG channel based on its structure obtained from cryo-electron microscopy and computer simulations, which reveal the critical roles of several specific structural domains and amino acid residues.

Potassium and sodium ion complexes with a partial peptide of the selectivity filter in K+ channels studied by cold ion trap infrared spectroscopy: the effect of hydration

Potassium channels allow K+ to rapidly diffuse, while the selectivity filter (SF) actively blocks Na+. The presence of water in the SF during ion translocation remains under debate due to the experimental and computational challenges in characterizing the interactions between water, ions, and the SF. Our bottom-up approach has been applied to a system composed of a partial peptide of the SF (Ac-tyrosine-NHMe) with a metal ion and a single water molecule to probe these interactions. The IR photodissociation spectra of M+Ac-tyrosine-NHMe(H2O) (M = Na, K) combined with quantum chemical calculations revealed that the water molecule binding sites are ion-dependent. In addition, the ion-peptide distances are elongated significantly for the K+ complex in comparison to the Na+ complex by the addition of a single water molecule. This striking structural difference with the water molecule is discussed in relation to ion selectivity and translocation within the K+ channel.

Interactions of Sea Anemone Toxins with Insect Sodium Channel-Insights from Electrophysiology and Molecular Docking Studies

Animal venoms are considered as a promising source of new drugs. Sea anemones release polypeptides that affect electrical activity of neurons of their prey. Voltage dependent sodium (Nav) channels are the common targets of Av1, Av2, and Av3 toxins from Anemonia viridis and CgNa from Condylactis gigantea. The toxins bind to the extracellular side of a channel and slow its fast inactivation, but molecular details of the binding modes are not known. Electrophysiological measurements on Periplaneta americana neuronal preparation revealed differences in potency of these toxins to increase nerve activity. Av1 and CgNa exhibit the strongest effects, while Av2 the weakest effect. Extensive molecular docking using a modern SMINA computer method revealed only partial overlap among the sets of toxins’ and channel’s amino acid residues responsible for the selectivity and binding modes. Docking positions support earlier supposition that the higher neuronal activity observed in electrophysiology should be attributed to hampering the fast inactivation gate by interactions of an anemone toxin with the voltage driven S4 helix from domain IV of cockroach Nav channel (NavPaS). Our modelling provides new data linking activity of toxins with their mode of binding in site 3 of NavPaS channel.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.

On the relationship between anion binding and chloride conductance in the CFTR anion channel

Mutations at many sites within the cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel pore region result in changes in chloride conductance. Although chloride binding in the pore - as well as interactions between concurrently bound chloride ions - are thought to be important facets of the chloride permeation mechanism, little is known about the relationship between anion binding and chloride conductance. The present work presents a comprehensive investigation of a number of anion binding properties in different pore mutants with differential effects on chloride conductance. When multiple pore mutants are compared, conductance appears best correlated with the ability of anions to bind to the pore when it is already occupied by chloride ions. In contrast, conductance was not correlated with biophysical measures of anion:anion interactions inside the pore. Although these findings suggest anion binding is required for high conductance, mutations that strengthened anion binding had very little effect on conductance, especially at high chloride concentrations, suggesting that the wild-type CFTR pore is already close to saturated with chloride ions. These results are used to support a revised model of chloride permeation in CFTR in which the overall chloride occupancy of multiple loosely-defined chloride binding sites results in high chloride conductance through the pore.

Unraveling of a Strongly Correlated Dynamical Network of Residues Controlling the Permeation of Potassium in KcsA Ion Channel

The complicated patterns of the single-channel currents in potassium ion channel KcsA are governed by the structural variability of the selectivity filter. A comparative analysis of the dynamics of the wild type KcsA channel and several of its mutants showing different conducting patterns was performed. A strongly correlated dynamical network of interacting residues is found to play a key role in regulating the state of the wild type channel. The network is centered on the aspartate D80 which plays the role of a hub by strong interacting via hydrogen bonds with residues E71, R64, R89, and W67. Residue D80 also affects the selectivity filter via its backbones. This network further compromises ions and water molecules located inside the channel that results in the mutual influence: the permeation depends on the configuration of residues in the network, and the dynamics of network's residues depends on locations of ions and water molecules inside the selectivity filter. Some features of the network provide a further understanding of experimental results describing the KcsA activity. In particular, the necessity of anionic lipids to be present for functioning the channel is explained by the interaction between the lipids and the arginine residues R64 and R89 that prevents destabilizing the structure of the selectivity filter.

Targeting thermoTRP ion channels: in silico preclinical approaches and opportunities

INTRODUCTION

A myriad of cellular pathophysiological responses are mediated by polymodal ion channels that respond to chemical and physical stimuli such as thermoTRP channels. Intriguingly, these channels are pivotal therapeutic targets with limited clinical pharmacology. In silico methods offer an unprecedented opportunity for discovering new lead compounds targeting thermoTRP channels with improved pharmacological activity and therapeutic index.

AREAS COVERED

This article reviews the progress on thermoTRP channel pharmacology because of (i) advances in solving their atomic structure using cryo-electron microscopy and, (ii) progress on computational techniques including homology modeling, molecular docking, virtual screening, molecular dynamics, ADME/Tox and artificial intelligence. Together, they have increased the number of lead compounds with clinical potential to treat a variety of pathologies. We used original and review articles from Pubmed (1997-2020), as well as the clinicaltrials.gov database, containing the terms thermoTRP, artificial intelligence, docking, and molecular dynamics.

EXPERT OPINION

The atomic structure of thermoTRP channels along with computational methods constitute a realistic first line strategy for designing drug candidates with improved pharmacology and clinical translation. In silico approaches can also help predict potential side-effects that can limit clinical development of drug candidates. Together, they should provide drug candidates with upgraded therapeutic properties.

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.

Long QT Syndrome Type 2: Emerging Strategies for Correcting Class 2 KCNH2 (hERG) Mutations and Identifying New Patients

Significant advances in our understanding of the molecular mechanisms that cause congenital long QT syndrome (LQTS) have been made. A wide variety of experimental approaches, including heterologous expression of mutant ion channel proteins and the use of inducible pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) from LQTS patients offer insights into etiology and new therapeutic strategies. This review briefly discusses the major molecular mechanisms underlying LQTS type 2 (LQT2), which is caused by loss-of-function (LOF) mutations in the KCNH2 gene (also known as the human ether-à-go-go-related gene or hERG). Almost half of suspected LQT2-causing mutations are missense mutations, and functional studies suggest that about 90% of these mutations disrupt the intracellular transport, or trafficking, of the KCNH2-encoded Kv11.1 channel protein to the cell surface membrane. In this review, we discuss emerging strategies that improve the trafficking and functional expression of trafficking-deficient LQT2 Kv11.1 channel proteins to the cell surface membrane and how new insights into the structure of the Kv11.1 channel protein will lead to computational approaches that identify which KCNH2 missense variants confer a high-risk for LQT2.

Long QT Syndrome Type 2: Emerging Strategies for Correcting Class 2 KCNH2 (hERG) Mutations and Identifying New Patients

Significant advances in our understanding of the molecular mechanisms that cause congenital long QT syndrome (LQTS) have been made. A wide variety of experimental approaches, including heterologous expression of mutant ion channel proteins and the use of inducible pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) from LQTS patients offer insights into etiology and new therapeutic strategies. This review briefly discusses the major molecular mechanisms underlying LQTS type 2 (LQT2), which is caused by loss-of-function (LOF) mutations in the KCNH2 gene (also known as the human ether-à-go-go-related gene or hERG). Almost half of suspected LQT2-causing mutations are missense mutations, and functional studies suggest that about 90% of these mutations disrupt the intracellular transport, or trafficking, of the KCNH2-encoded Kv11.1 channel protein to the cell surface membrane. In this review, we discuss emerging strategies that improve the trafficking and functional expression of trafficking-deficient LQT2 Kv11.1 channel proteins to the cell surface membrane and how new insights into the structure of the Kv11.1 channel protein will lead to computational approaches that identify which KCNH2 missense variants confer a high-risk for LQT2.

Selectivity filter modalities and rapid inactivation of the hERG1 channel

The human ether-á-go-go-related gene (hERG1) channel conducts small outward K⁺ currents that are critical for cardiomyocyte membrane repolarization. The gain-of-function mutation N629D at the outer mouth of the selectivity filter (SF) disrupts inactivation and K⁺-selective transport in hERG1, leading to arrhythmogenic phenotypes associated with long-QT syndrome. Here, we combined computational electrophysiology with Markov state model analysis to investigate how SF-level gating modalities control selective cation transport in wild-type (WT) and mutant (N629D) hERG1 variants. Starting from the recently reported cryogenic electron microscopy (cryo-EM) open-state channel structure, multiple microseconds-long molecular-dynamics (MD) trajectories were generated using different cation configurations at the filter, voltages, electrolyte concentrations, and force-field parameters. Most of the K⁺ permeation events observed in hERG1-WT simulations occurred at microsecond timescales, influenced by the spontaneous dehydration/rehydration dynamics at the filter. The SF region displayed conductive, constricted, occluded, and dilated states, in qualitative agreement with the well-documented flickering conductance of hERG1. In line with mutagenesis studies, these gating modalities resulted from dynamic interaction networks involving residues from the SF, outer-mouth vestibule, P-helices, and S5-P segments. We found that N629D mutation significantly stabilizes the SF in a state that is permeable to both K⁺ and Na⁺, which is reminiscent of the SF in the nonselective bacterial NaK channel. Increasing the external K⁺ concentration induced "WT-like" SF dynamics in N629D, in qualitative agreement with the recovery of flickering currents in experiments. Overall, our findings provide an understanding of the molecular mechanisms controlling selective transport in K⁺ channels with a nonconventional SF sequence.