- 1Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
- 2Beijing Institute of Big Data Research, Beijing, China
- 3DP Technology, Beijing, China
- 4School of Pharmaceutical Sciences, Peking University, Beijing, China
- 5National Engineering Research Center of Visual Technology, Peking University, Beijing, China
Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.
Introduction
Ion channels are a group of pore-forming proteins located at the membrane of cells or intracellular organelles, controlling the flow of ions across the membranes. Such transportation of ions plays a critical role in the physiology of all living cells (Hille, 1978; Hille, 1986), thus the malfunction of ion channels leads to numerous diseases at the fundamental level. Although spanning a broad spectrum of families, ion channels can be roughly categorized into the voltage-gated and the ligand-gated based on activation mechanisms and structural similarities (Figure 1). Most voltage-gated ion channels contain four individual subunits or four repeats within one continuous polypeptide, each includes a similar six-transmembrane-helix core architecture, which can be further divided into a four-helix voltage sensor domain S1-S4 and the pore domain S5-S6 connected by the pore loop (Gulbis et al., 1999). The typical voltage-gated ion channels, mainly the sodium (Catterall, 2000), potassium (Wulff et al., 2009), and calcium channels (Catterall, 2011), are mostly depolarization activated and responsible for fast reactions upon voltage changes at the cell membrane, thereby playing a key part in electric signaling (Catterall, 2000; Sands et al., 2005; Catterall, 2010). While some members of the transient receptor potential (TRP) family, mainly facilitate their biological function via other stimuli including temperature, force, and chemical compounds. They serve as the vanguard of the sensory system (Clapham, 2003). As for ligand-gated ion channels, they are mostly activated by ligands like neurotransmitters and are crucial for nervous activities. There are three families in mammals: Cys-loop receptor ion channels (Sine and Engel, 2006), Glutamate receptor ion channels (Traynelis et al., 2010), and ATP gated ion channels (Schmid and Evans, 2019), they are typically composed of a transmembrane domain forming the pore and an extracellular domain carrying the ligand-binding sites. Besides, some other ion channels are beyond these two types, for example, the ubiquitously expressed trimeric Piezo 1/2 channels that respond to mechanical stimuli (Coste et al., 2010).
FIGURE 1. Classification of ion channels based on activation mechanisms and structural similarities. The evolutionary relationships among these ion channels are measured by the phylogenetic tree of ion channels (Alexander et al., 2019) with representative channel structures (center). Representative ion channels selected from each class are shown with their 3D molecular structures.
Due to the vital function of ion channels in physiology, mutations of ion channels were identified to be responsible for various diseases: inherited neuronal diseases like epilepsy (Escayg et al., 2000), cardiac diseases like Long-QT syndrome (Keating, 1996), muscle diseases like Hyperkalemic periodic paralysis (Brugnoni et al., 2022) and many others (Ackerman and Clapham, 1997). Also, the blockage of ion channels can directly regulate sensory perception such as pain (Tominaga and Caterina, 2004; Bourinet et al., 2014; Bennett et al., 2019). Consequently, ion channels are among the most pursued drug targets in the past decades (Bagal et al., 2013) and several drugs were developed against a wide range of diseases (Clare, 2010). To be specific, for the most related central nervous system (CNS) diseases: Diazepam, the famous anxiolytic launched in 1963 and the best-seller in the United States between 1968 and 1982 (Calcaterra and Barrow, 2014), was identified as a positive allosteric modulator of the gamma-aminobutyric acid (GABA) type-A receptors in 1977 (Braestrup and Squires, 1977; Mohler and Okada, 1977). A detoured but successful development of ketamine, the N-methyl-D-aspartate (NMDA) receptor antagonist (Alessandri et al., 1989; Zhang et al., 2021) into an anti-depressant drug is even considered a sign of a renaissance in the psychiatric drug industry (Daly et al., 2018; Reardon, 2018). Recent studies also revealed more potential applications of channel modulators in Parkinson’s (Daniel et al., 2021) and Alzheimer’s diseases (Tan et al., 2012; Gonzales and Sumien, 2017; Lu et al., 2017). Additionally, the prevalent expression makes ion channels adequate targets for other diseases as well. Varenicline is derived from the modification of (-)-cytisine, an alkaloid that could be found in nature (Coe et al., 2005). As a partial agonist of nicotinic receptors with the beta 2 subunit (Papke and Heinemann, 1994), Chantix (varenicline tartrate) was launched by Pfizer in 2006 for smoking cessation (Coe et al., 2005) and is still one of the best-selling drugs in the world with annual sales of 919 million dollars in the year of 2020 (McGrath et al., 2010; Njardarson et al., 2020). In addition, cystic fibrosis, a heritable disease, is caused by a defect in the ATP-activated chloride channel, cystic fibrosis transmembrane conductance regulator (CFTR) which is involved in water flow control during the production of sweat, digestive fluids, and mucus (Bobadilla et al., 2002). The latest treatment for cystic fibrosis, Trikafta, was developed by Vertex Pharmaceuticals and hit the market in 2019, and was considered at the top of 2019's new approvals by expected revenue (Mullard, 2020). Another case for diabetes treatment is the blocking of ATP-sensitive potassium channel KATP by Glipizide leads to depolarization of beta cells, resulting in the opening of voltage-gated calcium channels which encourages insulin release (Shuman, 1983). Ion channels were also reported to be potential anti-cancer targets. Their role in proliferation, migration, and metastasis of cancer was well-explored (Litan and Langhans, 2015) and the over-expression of TRPA1 and TRPV1 channels was observed in several tumors suitable for topical administration, allowing for a safe and efficient therapy (Kiss et al., 2020).
Despite the vast opportunities, the discovery of drugs targeting ion channels faces unique obstacles. To begin with, the development of biological assays for ion channels is usually challenging, and the low expression rates further hindered the high-throughput screening. Ever since Hodgkin and Huxley firstly described the movement of ions in nerve cells of squid axons during an action potential (Hodgkin and Huxley, 1952), electrophysiology has been used to study the function of ion channels. Yet, such methods require considerable expertise and are labor-intensive and time-consuming, therefore not realistic for large-scale screening, even with the great advancement in automated patch clamp devices (Dunlop et al., 2008; Kodandaramaiah et al., 2012). Another major stumbling block is ligand selectivity because channels from the same or close families usually share great similarities in sequence and structure yet bear diverse functions. Lead compounds, especially channel blockers, often suffer from side effects due to poor selectivity (Vandenberg et al., 2017; Flood et al., 2019). As the “low-hanging fruit” is becoming rare in modern drug discovery, cases like carbamazepine, a pan-ion-channel inhibitor with broad pharmacological properties (Bagal et al., 2013), may be very rare and not welcomed nowadays, for which mechanism studies came several decades later than discovery (Schindler and Hafliger, 1954; Grant and Faulds, 1992).
Because of the aforementioned problems, the rational design of drugs is critical for the efficiency and efficacy of ion channel drug discovery. Such designs require precise structural information, which is largely missing until the rise of single-particle cryo-electron microscopy (cryo-EM) (Li et al., 2013; Liao et al., 2013). Starting from the high-resolution structures, computer-aided drug design (CADD) can be applied to perform large scale virtual screenings before the costly wet-lab experiments to greatly narrow down the size of the library to an affordable range (Lyne, 2002); molecular dynamics (MD) simulations can be performed to further model the structures and study the mechanism at atomistic details (Maffeo et al., 2012). On the other hand, recent advancement in machine learning begins to demonstrate their power in cryo-EM and drug development. Towards accelerating the ion channel drug discovery, in this review, we will provide a broad overview of the current machine-learning and simulation-based techniques in structural biology and drug design, and discuss how they can be applied to ion channel research.
Computational Methods for Single-Particle Cryo-EM
While being essential for mechanistic studies and rational drug design, molecular structures for ion channels are more difficult to purify and crystallize than soluble proteins, mainly due to the necessity and difficulty of preserving the membrane-like environment. Detergents, amphipols, and nanodiscs have been commonly used to extract membrane proteins such as ion channels, and also serve as substitutes for the local membrane to stabilize the transmembrane domains (Zampieri et al., 2021). In recent years, the single-particle cryo-EM technique has rapidly evolved as a powerful method for structure determination for various ion channels (Martin et al., 2017; Basak et al., 2019; Dang et al., 2019; Masiulis et al., 2019; Zhao et al., 2019; Wang Q. et al., 2020; Sun and MacKinnon, 2020; Lin et al., 2021; Song et al., 2021; Yu et al., 2021). Cryo-EM eliminates the need for protein crystallization and is capable of resolving multiple protein conformations within a single cryo-EM dataset (Frank and Ourmazd, 2016; Frank, 2017; Cheng, 2018; Wu and Lander, 2020). The structural resolution of cryo-EM experiments has also been improved over the years, with the highest resolution reaching ∼1.2 Å (Nakane et al., 2020; Yip et al., 2020). High-resolution cryo-EM enables its use in target identification and lead-optimization for developing drugs targeting ion channels.
A typical single-particle cryo-EM experiment includes sample purification, vitrification via plunge freezing, electron microscopy data collection, image processing, and structural model building (Figure 2) (Lyumkis, 2019; Wu and Lander, 2020). We will focus on the last two steps in this review, which rely heavily on computational techniques.
FIGURE 2. Singe-particle cryo-EM workflow and relevant computational methods. Traditional methods for image processing and structural model building are shown in grey and ML-based methods in blue. Electron micrographs and particle images were visualized using cryoSPARC (Punjani et al., 2017). The density map and structural model were generated using Mol* viewer (Sehnal et al., 2021).
Cryo-EM Image Analysis
The goal of cryo-EM image analysis is to reconstruct the three-dimensional (3D) structure of the target proteins represented as statistically estimated density maps. A large collection of two-dimensional (2D) projection maps are subjected to image preprocessing, particle picking, 2D image classifications, 3D classifications, and 3D reconstructions (Singer and Sigworth, 2021). Many software programs are available for building cryo-EM density maps from 2D projection maps, such as RELION (Scheres, 2012; Kimanius et al., 2021), cryoSPARC (Punjani et al., 2017), Xmipp (Sorzano et al., 2004), Warp (Tegunov and Cramer, 2019), CisTEM (Grant et al., 2018), EMAN2 (Tang et al., 2007). Although these programs are sufficient for most cryo-EM experiments, their output density maps are mostly limited to a few representative classes of particles in the projection map dataset. Recent machine-learning-based methods focus on characterizing the continuous distribution of particle distributions within a cryo-EM experiment, e.g., cryoDRGN (Zhong et al., 2020; Zhong ED. et al., 2021; Zhong E. D. et al., 2021) and e2gmm (Chen and Ludtke, 2021). CryoDRGN is a method for heterogeneous cryo-EM reconstruction based on deep neural networks. It is powerful and general for analyzing structural heterogeneity in macromolecular complexes of various sizes and degrees of heterogeneity. The cryoDRGN model consists of two neural networks structured in an image-encoder-volume-decoder architecture with a continuous latent variable representation for describing the sample heterogeneity. Although the initial version of cryoDRGN (Zhong et al., 2020) can optimize image poses together with image reconstruction/classification, for better efficiency and accuracy, the pose estimation component was dropped and pre-computed image poses estimated from an upstream homogeneous reconstruction were used (Zhong ED. et al., 2021). CryoDRGN2 combines a traditional cryo-EM image pose search algorithm with deep neural networks, which eliminates the need for pre-computed image poses (Zhong E. D. et al., 2021). Compared with the initial version of cryoDRGN, both the speed and accuracy are increased significantly in pose estimation. The e2gmm algorithm can also find the structural variability of target macromolecules captured by cryo-EM experiments. In this method, each structural domain is modeled by a 3-dimensional (3D) Gaussian function, parameterized by its center coordinates, amplitude, and width. Autoencoder-like neural network architecture is used to generate a latent vector representation of each projection image through an encoder. Then the decoder component converts each latent vector into a set of 3D Gaussian function parameters which represents a Gaussian mixture model (GMM). This GMM representation is converted back to a 2D projection image through a projection operation. In e2gmm, training was first done for the decoder only, then for the encoder and decoder together. Unlike the standard autoencoder, the input to the encoder is not the projection image, but the gradient of loss function score with respect to the GMM parameters (Chen and Ludtke, 2021). The projection orientation of each image must be given as inputs which can be refined later once the network has been trained. The conformational variability of the cryo-EM dataset can be visualized through dimension reduction of the 4D latent space based on the latent vectors from the encoder. The cryoDRGN and e2gmm algorithms are similar in terms of the use of latent vectors for characterizing the structure variability, although latent vectors are obtained differently.
Reconstructed cryo-EM density maps usually suffer from loss of contrast at high resolution due to factors such as specimen movement, radiation damage, particle flexibility/heterogeneity, and deficiencies in extracting/averaging micrograph signals in the reconstruction algorithms (Rosenthal and Henderson, 2003). Algorithms, known as sharpening, have been designed to reduce this type of contrast loss. A commonly used sharpening method is based on global B-factor correction (Rosenthal and Henderson, 2003). Other methods, such as LocScale (Jakobi et al., 2017), LocalDeblur (Ramirez-Aportela et al., 2020), and LocSpiral (Kaur et al., 2021), improved upon the previous method by considering local variations in B-factors. Instead of explicitly performing sharpening calculations, deep-learning-based methods can implicitly encode the sharpening transformation. One example is DeepSharpen (Zehni et al., 2020), which uses a Convolutional Neural Network (CNN) trained on ∼15,000 pairs of low and high-resolution synthetic cryo-EM density maps. The benefit of this approach is that sharpening on a new density map only requires inference through trained CNN parameters which can be very fast. Deep cryo-EM Map Enhancer (DeepEMhancer) is another deep-learning-based sharpening method (Sanchez-Garcia et al., 2021). It uses a 3D U-net architecture (Ronneberger et al., 2015) to learn the local sharpening effect of the LocScale algorithm by training on paired map datasets from the Electron Microscopy Data Bank (EMDB). Each map pair consists of an experimental density map and a post-processed map using LocScale. DeepSharpen and DeepEMhancer share similar deep learning designs, although their training strategies and training targets differ. More realistic and better-labeled training datasets can further improve these supervised-learning-based methods.
Cryo-EM Structural Model Building
In the final step of cryo-EM experiments, 3D atomic structural models are built based on the reconstructed 3D density maps. The goal is to search for the most probable 3D atomic structure that can best interpret the reconstructed density map. Traditionally, these atomic structural models are built semi-automatically using software like Coot (Casañal et al., 2020) and Phenix (Liebschner et al., 2019). These types of methods require the users to have extensive knowledge of the biochemical compositions of the target biological entities. More automated methods can be roughly divided into three categories. The first class of automated methods performs rigid-body fitting of initial structures into the density maps. Examples of such methods include EMfit (Rossmann et al., 2001), HermiteFit (Derevyanko and Grudinin, 2014), GTMC (Wu et al., 2003), ADP_EM (Garzón et al., 2007), and Situs-colores (Chacón and Wriggers, 2002). These methods solve the problem of placing an atomic structure model into a target cryo-EM density map as a rigid body and they differ in their specific implementation of translational/rotational search. Most of these methods are developed for low-resolution maps before the so-called “resolution revolution” during the development of cryo-EM (Crowther, 2016; Frank, 2017; Alnabati and Kihara, 2019).
As the average resolution of cryo-EM increased, the second class of automated methods emerged which allows flexible fitting of initial structures into density maps. One type of flexible fitting method combines molecular dynamics (MD) simulations and the experimental density maps to automatically fit some initial structural models into density maps. Examples of such methods include molecular dynamics flexible fitting (MDFF) and correlation-driven molecular dynamics (CDMD). In MDFF (Trabuco et al., 2008; Chan et al., 2011; Singharoy et al., 2016), an additional potential energy term derived from the density map was added to the simulation force field. The derivative of this potential generates forces that drive the initial structural model into the target density map. CDMD acts similarly (Igaev et al., 2019). Instead of directly adding a density-map-derived potential energy term, the cross-correlation coefficient (ccc) between the intermediate structures and the target density map is calculated. Forces derived from ccc drive the initial structure into the density map. Although both methods are equivalent, CDMD is more sensitive to the quality of initial structures and noises in the density maps compared to MDFF in practice. ISOLDE is a hybrid method that combines user interactions and MDFF to achieve a semi-automated structure modeling (Croll, 2018). The same hybrid modeling can also be performed using interactive MDFF (iMDFF) (McGreevy et al., 2016). NMMD is another MD-based cryo-EM structural fitting method (Vuillemot et al., 2022). Compared to CDMD and MDFF, NMMD combines normal mode analysis (NMA) and MD to perform global structural fitting with NMA and local structural fitting with MD. Other NMA-based flexible fitting methods include MDeNM-EMFit (Costa et al., 2020), NMFF (Tama et al., 2004a; b), NORMA (Suhre et al., 2006), and a coarse-grained fitting method with a modified elastic network model (Zheng, 2011). Besides MD, Monte Carlo (MC) simulations together with MD can also be used to fit structures into density maps (Topf et al., 2008).
Another approach to building cryo-EM structures is de novo modeling, which doesn’t require any initial structural models. This class of methods includes RosettaES (Frenz et al., 2017), EM-fold (Lindert et al., 2009; Lindert et al., 2012a; Lindert et al., 2012b), Gorgon (Baker et al., 2011), Pathwalking (Baker et al., 2012; Chen et al., 2016), Phenix (Terwilliger et al., 2018; Terwilliger et al., 2020), MAINMAST (Terashi and Kihara, 2018b; a), etc. (see the recent review (Alnabati and Kihara, 2019)). More recently, deep-learning-based de novo cryo-EM structural model building has been gaining much ground, such as DeepTracer (Pfab et al., 2021), CR-I-TASSER (Zhang et al., 2022), DeepMM (He and Huang, 2021), and A2-Net (Xu et al., 2019). DeepTracer (Pfab et al., 2021) uses four U-Net-based neural networks to extract structural information from the input density amp. Each network can classify each density map voxel into a certain category. The Atoms U-Net predicts four classes: C-alpha, nitrogen, carbon (non-C-alpha), and non-atom. The Backbone U-Net predicts three classes: backbone, side chain, and neither. The Secondary Structure U-Net predicts four classes: loop, helix, sheet, and none of these. The Amino Acid Type U-Net predicts 21 classes: 20 types of amino acids and non-protein density. Then, protein backbone tracing is done in three steps: identifying protein chains based on Backbone U-Net outputs, estimating C-alpha atom coordinates, and connecting C-alpha atoms into chains using a modified traveling salesman algorithm. The assignment of amino acid residue types to each C-alpha atom was done through a customized sequence alignment algorithm based on the U-Net-predicted and input protein sequences. The all-atom backbone structure is subsequently reconstructed based on information from the Atoms U-Net predictions and some assumptions on the peptide bond geometries. One of the major limitations of DeepTracer is its inability to predict the side-chain conformations, which is delegated to an external program called SCWRL4 (Krivov et al., 2009). CR-I-TASSER (cryo-EM iterative threading assembly refinement) is a hybrid method that uses information from both homology modeling and density-map-based de novo C-alpha tracing (Zhang et al., 2022). Like DeepTracer, C-alpha atom positions are predicted based on density maps using machine learning techniques, which is a deep convolutional network (3D-CNN) in this case. Then, CR-I-TASSER solves the problem of C-alpha tracing using a local meta-threading server (LOMETS). After that, the atom structure models are built using the iterative threading assembly refinement method (I-TASSER) driven by model-map correlation under deep-learning boosted template restraints. The novelty of this method is the C-alpha position prediction using 3D-CNN and the use of C-alpha locations for better template reselection in LOMETS. CR-I-TASSER shows good performance in terms of TM scores compared to other methods such as MAINMAST and MDFF. However, its accuracy is limited for low-resolution cryo-EM data, and its C-alpha tracing method is also problematic in long-loop/tail or disordered regions. Another major limitation of CR-I-TASSER is that it can only handle single-chain proteins, thus the density maps have to be manually segmented first. DeepMM takes a workflow design similar to DeepTracer (He and Huang, 2021). Based on cryo-EM density maps, it uses one Densely Connected Convolutional Network (DenseNet) to predict the main chain (backbone N, C, and C-alpha atoms) and C-alpha probabilities at each density map voxel. Using the MAINMAST method, the main-chain paths are determined based on the main-chain probability map. Then, another DenseNet is used to predict the residue identities and secondary structure types for each main chain point. Smith-Waterman dynamic programming is used to align the target sequence to the main-chain paths. Finally, all-atom structure models are based on the top-10 C-alpha models using the ctrip program (Xiang and Honig, 2001; Petrey et al., 2003) and refined by energy minimization using Amber (Case et al., 2005). Like CR-I-TASSER, DeepMM can only handle building single protein chains and segmentation tools, such as Segger (Pintilie et al., 2010) are needed to preprocess multi-chain density maps into individual chain map segments. Compared to DeepTracer, DeepMM’s main feature is the ability to build full-length protein structures from high-resolution EM maps, while DeepTracer only builds backbone segments in high-resolution map regions. All these methods only apply machine learning techniques to convert cryo-EM maps into basic protein structural models without side chains. A2-Net (Xu et al., 2019) is another de novo structure building method based on a deep Convolutional Neural Network (Srivastava et al., 2015). Like the methods mentioned above, it contains neural networks for amino acid residue identity detection. The main difference is that A2-net also contains a poseNet based on the 3D stacked hourglass network (Newell et al., 2016) for estimating the coordinates of each residue (side chain included). Then, these residues are connected to form the target protein using a Monte Carlo tree search algorithm. In the advent of accurate protein structure prediction methods such as AlphaFold2 (Jumper et al., 2021) and RoseTTAFold (Baek et al., 2021), new possibilities are brought to cryo-EM structure building. For example, ChimeraX added a new feature to allow users to start cryo-EM structure modeling from AlphaFold predictions1. Phenix also added a new workflow2 (Terwilliger et al., 2022), to iteratively refine cryo-EM structural models by integrating AlphaFold predictions using ColabFold (Mirdita et al., 2022). Recently, Yan and Shen reported that AlphaFold2-predicted structure could help model the density map regions with moderate resolutions in voltage-gated sodium channel Nav1.7 (Huang et al., 2022). With the rapid progress in the field of machine learning, more accurate and efficient de novo cryo-EM structure modeling methods are expected in the near future.
Evaluation of the Cryo-EM Structural Model Quality
Evaluation of the quality of the density-map-derived molecular structural models is yet another challenge. Global metrics (such as map-model cross-correlation coefficient and FSC curve) between the experimental density map and a simulated map converted from the structural model are often insensitive to local mis-fittings (Hryc et al., 2017). Moreover, these metrics usually do not account for densities not modeled by the atomic structural models such as lipids and detergents, and masks are usually needed to access local features (Pintilie et al., 2016). EMRinger and Z-scores are metrics defined based on how well the structure models interpret the experimental density maps. EMRinger score is based on the expectation that the map density values near the beta-carbon atoms should peak at certain angles when rotating the side-chain chi-1 dihedral angle (Barad et al., 2015). Although the EMRinger score correlates well with the protein backbone modeling quality and map density quality near the beta-carbons, it is suspectable to noise and lacks assessment of the rest of the side chains. Z-score has been proposed to quantify the density map quality at the secondary structure element level and side-chain level (Pintilie and Chiu, 2018). Its definition is based on the cross-correlation coefficients between the experimental and simulated maps of the secondary structure elements or side chains with and without some geometric displacements (Pintilie and Chiu, 2018). To make even finer-scale assessments, Q-score is established to quantify the resolvability of individual atoms in cryo-EM maps (Pintilie et al., 2020). The basic idea of the Q-score is to measure the cross-correlation coefficient between experimental density map values for individual atoms and ideal densities represented by a three-dimensional (3D) Gaussian-like function centered at each atom in the structure model (Pintilie et al., 2020). FSC-Q is another method that can also assess structure modeling quality on a per-atom basis (Ramírez-Aportela et al., 2021). FSC-Q calculates the difference between two local resolution maps calculated using the blocres program. The first local resolution map is calculated using the full density map and the simulated map based on the structure model, while the second is calculated using the half maps. Then the FSC-Q difference map can be projected onto atomic models to obtain an FSC-Q score for each atom. Alternatively, the FSC-Q map can be normalized by dividing by the local resolution map based on the half maps for better cross-comparison between experimental density maps of different resolutions.
Challenges in Cryo-EM Structure Determination for Ion Channels
Ion channels naturally reside in biological membranes which means preserving the amphiphilic local environment for the channel transmembrane domains is essential during cryo-EM sample preparation. Techniques such as nanodiscs can serve as good replacement for the membrane (Zampieri et al., 2021). However, the disordered nature of these membrane memetics also add difficulties to the cryo-EM structural analysis. In cryo-EM image analysis, algorithms such as non-uniform refinement can reduce the negative impact of the disordered detergent/lipid molecules (Punjani et al., 2020). In addition, certain ion channels, such as mechanosensitive channels, have flexible transmembrane domains and often require specific protein-lipid interactions for structural stability (Kefauver et al., 2020). Preserving their local membrane environment can be very difficult. For example, it has been reported that mechanosensitive channel YnaI’s transmembrane domain couldn’t be resolved using cryo-EM with membrane-active polymer SMA2000 except for core helices TM4 and TM5 (Catalano et al., 2021). Computational techniques such as normal model analysis and molecular dynamics can offer additional insights into these protein-lipid interactions required for the structural stability and functions of mechanosensitive channels (Jojoa-Cruz et al., 2018; Argudo et al., 2019). In addition, the application of toxins, truncation of flexible regions and mutations are often introduced in cryo-EM studies of ion channels (Jiang D. et al., 2021). Although sometimes intended by the researchers for their specific projects, these conditions have also been used to artificially stabilize ion channels to capture a certain conformation during structure determination (Rook et al., 2021). Computational techniques can reverse these effects if needed and can also be used to evaluate the impact of such conditions (Gupta and Vadde, 2021; Saikia et al., 2021).
Computer-Aided Drug Design Approaches Targeting Ion Channels
In the glorious history of the era driven by occasions and classic synthetic chemistry, numerous ion-channel modulators were discovered. For example, the series of “-caine” analogues was inspired and chemically modified from cocaine, an alkaloid derived from plants (Plowman, 1982), including benzocaine, amylocaine, procaine, procainamide, and lignocaine (lidocaine), which showed potency as sodium channel inhibitors and were approved for medical use as anesthetics and/or antiarrhythmics (Cox, 2015). Typical cases also include the previously mentioned varenicline and carbamazepine, and maybe even ethanol and neurotoxic peptides in venoms. Ethanol was proved to be able to interact with various post-synaptic ion channel receptors (Howard et al., 2011; Murail et al., 2012; Zhang et al., 2016), while neurotoxic peptides in venoms were sometimes regarded as good starting points for modification as modulators targeting ion channels such as voltage-gated sodium, potassium, and calcium channels as well as some ligand-gated channels (Clark et al., 2010; Dutertre and Lewis, 2010; Bagal et al., 2013; Wulff et al., 2019).
In recent years, under the concept of rational drug design, considering the features of ligands, the structures of receptors, and the motions and dynamics of the system becomes more and more important. Accordingly, numerous computational methods and machine-learning algorithms have been developed (Figure 3).
FIGURE 3. Computational approaches in the structure modeling, mechanistic study and drug discovery of ion channels. Several subplots were collected using Hermite (https://hermite.dp.tech), Mol* viewer (Sehnal et al., 2021), PyMOL (https://github.com/schrodinger/pymol-open-source) and the Pred-hERG web-server (Braga et al., 2015).
Ligand-Based Approaches
Benefit from these data collected in the early era, ligand-based approaches could be applied to build relatively reliable models to describe the relationship between modulators and the corresponding ion channels, especially when the experimentally-resolved protein structures are absent. One of the typical cases should be the development of ligand-based models for the prediction of human ether-a-go-go related gene hERG potassium ion channel blockage, which is believed to be important in avoiding sudden cardiac death considering its relationship with drug-induced long QT syndrome. In the past few years, a considerable amount of associated data was collected (Konda et al., 2019) and numerous quantitative structure-activity relationship (QSAR) or classification models (Cavalli et al., 2002; Ekins et al., 2002; Roche et al., 2002; Keseru, 2003; Nisius and Goller, 2009; Braga et al., 2015; Siramshetty et al., 2018; Konda et al., 2019) were built for hEGR blockage prediction based on the traditional regression models such as partial least squares, machine learning algorithms including support vector machines (SVM), random forest (RF), gradient boosting machine (GBM) or neural network models. These models have been systematically discussed in the previous reviews (Aronov, 2005; Taboureau and Jorgensen, 2011; Villoutreix and Taboureau, 2015; Vandenberg et al., 2017; Menke et al., 2021).
Actually, besides hERG, the application of the mentioned models is also a common strategy in ligand-based drug design of ion channel modulators when facing abundant research data. For example, in the study of structural-activity relationships of Nav1.4 blockers, both traditional two-dimensional and advanced three-dimensional QSAR models were applied by Carrieri et al. to reveal the structural characteristics of bioactive molecules (Carrieri et al., 2009). Further design of sodium channel blockers was thus inspired. Consequently, a series of tocainide analogues were designed and one of them was verified to have a favorable pharmacodynamic profile and was proposed to be a valid Nav1.4 blocker (Muraglia et al., 2014). Based on the data collected from ChEMBL, BindingDB, and in-house databases, various machine learning-based QSAR models were built by Kristam et al. to predict blockers of the voltage-gated sodium ion channel Nav1.5, with the balanced accuracy of 0.88 (at the threshold of 1 μM) and predicted R2 of 0.71 (RMSE = 0.73 for pIC50) for the classification and regression models, respectively (Khalifa et al., 2020). Similarly, Huang and Xie et al. trained and filtered a classification model based on the data from ChEMBL and BindingDB for the discovery of Nav1.7 blockers (Kong et al., 2020). The Grammar Variational Autoencoder, the trained classification model, and simulated annealing were combined to conduct the molecular optimization and an active compound was found and identified experimentally (Kong et al., 2020). Li et al. studied the pharmacophore hypothesis of ligands that bind at the benzodiazepine site of GABAA receptors based on ligand-based pharmacophore, 3D-QSAR analysis, and Bayesian models, which might provide useful viewpoints in the discovery of GABAA modulators (Yang Y. et al., 2013). Because of the potentially limited chemical space the pore-blockers could adapt, approaches based on the three-dimensional similarity between the ligands and reference molecules that take the spatial molecular shape and electrostatic features into consideration (Cleves et al., 2019; Jiang Z. et al., 2021) are possibly another practical ligand-based method in discovering ion channel modulators. Additionally, Bahar et al. proposed a probabilistic matrix factorization (PMF) based scheme to predict unknown drug-target interactions, which focused on the phenotypic similarity of the drugs by grouping the drugs according to the corresponding therapeutic effects (Cobanoglu et al., 2013). The obtained model showed good performance when applied to large data sets of ion-channel-drug pairs (Cobanoglu et al., 2013).
Structure-Based Approaches
In terms of the structural-based approaches, recent computational and machine-learning technologies have empowered the exploration of ion channel modulators mainly in three aspects, viz., structures with higher quality, binding affinity prediction with higher accuracy, and virtual screening with higher throughput (Figure 3).
The recent rapid development of protein structure modeling, as well as structural biology, possibly enables the rational design and exploration of ion channel modulators with atomic-resolution structures. Taking again the case of hERG as an example, in the past, owing to the lack of solved hERG channel structure, structural-based exploration mainly relies on structures from homology modeling. With the modeled structure and docked model, Vaz et al. illustrated that the π-stacking formed between the blockers and Phe656 in the protein as well as that of the cation-π interaction involving Tyr652 might be the potential key elements of the pharmacophore (Pearlstein et al., 2003). Similar conclusions associated with the important character of Phe656 were also drawn by the work of Åqvist et, al., which explored the binding of potential hERG blockers using docking and molecular dynamics combined with the linear interaction energy method to evaluate the binding affinities of the potential blockers (Osterberg and Aqvist, 2005). However, numerous experimental ion channel structures have been determined with the cryo-EM technique nowadays, including the human hERG structures that were resolved in 2017 (Wang and MacKinnon, 2017), which might provide a significant cornerstone of the structure-based prediction. Beneficial from the high-resolution structures of the TRPM8 channel in various states (Yin et al., 2018; Diver et al., 2019; Yin et al., 2019), Yang et al. conducted a rational design of the TRPM8 channel and developed a modality-specific inhibitor DeC-1.2, which showed the potential as a novel analgesic against oxaliplatin-induced neuropathic pain (Aierken et al., 2021). Additionally, various protein-folding and modeling tools were also developed (Baek et al., 2021; Jumper et al., 2021), which provide an alternative way for protein structure modeling besides the classic homology modeling approach (Huang et al., 2022), and may also contribute as a start point for the construction of single-particle cryo-EM structures (Figure 3). The application of molecular dynamics may be another powerful approach to structural exploration and refinement (Hollingsworth and Dror, 2018; Carnevale et al., 2021). For example, following a multistep protocol, Noskov et al. modeled the open, closed and open-inactive states of the hERG channel to explore the molecular mechanisms behind the state-dependent binding of potential hERG blockers (Durdagi et al., 2012). The missing elements in the homology model were complemented with Rosetta, following all-atom molecular dynamics simulation to refine the obtained structure. The Poisson-Boltzmann calculations conducted based on the ligand-channel complex obtained by the Glide induced-fit docking (IFD) protocol showed a good correlation with the experimental results (Durdagi et al., 2012). Also, starting from the modeled structure, Barakat et al. performed the structural refinement of the Nav1.5 ion channel based on a long-time molecular dynamics simulation of 680 ns in the lipid membrane bilayer (Ahmed et al., 2017). The binding affinities between the reported molecules and the obtained model of the channel were evaluated with AMBER-MM/GBSA. The results demonstrated a good computational-experimental correlation (Pearson coefficient = 0.7) to distinguishing drugs that block or not block the Nav1.5 channel (Ahmed et al., 2017).
Binding affinities of the ion channel modulators are usually evaluated after the construction of reliable modulator-channel complexes. As discussed above, numerous algorithms are available in the historical explorations for binding affinity prediction, including LIE, MM/PBSA, and MM/GBSA, and showed relatively satisfactory results (Osterberg and Aqvist, 2005; Durdagi et al., 2012; Ahmed et al., 2017). Nevertheless, to derive the theoretically more accurate results, the more rigorous defined alchemical methods such as free energy perturbation (FEP) could also be applied to access the binding potency of the channel modulators (Limongelli, 2020). Based on the extensive search for the isoflurane binding sites in the proton gated Gloeobacter violaceus ligand-gated ion channel (GLIC) in their previous study (Brannigan et al., 2010), Klein and Brannigan et al. performed free energy perturbations to evaluate the binding affinities of two anesthetics against GLIC, viz., isoflurane and propofol (LeBard et al., 2012). The evaluated affinities demonstrated good correlations with the previous electrophysiology data (LeBard et al., 2012). In the work of Rempe and Ren et al., FEP was used to determine accurately the thermodynamic stabilities of ion configurations for the potassium channel of streptomyces A (KcsA), by comparing the relative free energy difference between various configurations (Jing et al., 2021). The computation results showed good agreement with evidence derived from crystal structures, ion binding experiments, and mutagenesis experiments (Jing et al., 2021). Moreover, machine learning techniques could be introduced in large-scale screening considering the large computational cost of alchemical methods. Bhat et al. showed that it is possible and acceptable to perform more than 5000 FEP calculations in one lead optimization task with the acceleration of active learning (Konze et al., 2019). System-specific scoring functions might provide another efficient way for discovering ion channel modulators. These scoring functions could be constructed either with deep learning or machine learning methods alone or combined with physical-based models (Ain et al., 2015; Wang et al., 2019; Guedes et al., 2021). Numerous studies have demonstrated that the application of target-specific scoring functions outperforms the universal scoring functions on the validation set of DUD-E, which usually covers the major druggable target classes and ion channels are also included (Ain et al., 2015; Wojcikowski et al., 2017; Wang et al., 2019; Guedes et al., 2021).
The development of ultra-high-throughput virtual screening might also be a chance to accelerate the discovery of ion channel modulators. Though various experimental techniques including ion-flux or fluorescence-based assays and electrophysiology measurements are developed for discovering ion channel modulators, however, most of the options are time-consuming and costly for high-throughput screening (Zhang et al., 2016). Under the logic that screening of larger libraries could generally lead to higher quality results (Gorgulla et al., 2020), a variety of programs and algorithms have been developed for ultra-high throughput virtual screening. Generally, the ideas for these approaches include optimization in algorithm engineering to enable high-performance computing (Alhossary et al., 2015; Hassan et al., 2017) and the usage of graphics processing units (Santos-Martins et al., 2021; Shidi et al., 2022), the development of highly-organized workflows (Gorgulla et al., 2020), the application of deep learning which trained QSAR models on docking scores of subsets of the library (Gentile et al., 2020), and the improvement in searching efficiency in a fragment-based manner with the idea partial similar with dynamic programming (Sadybekov et al., 2022). These methods have pushed the limit of virtual screening throughput to libraries with more than tens of billion compounds, and may potentially help the drug discovery of various target classes including ion channels.
Computational Approaches for Functional Analysis of Ion Channels
Recent advances in structural biology such as high-resolution cryo-EM have enabled the in-depth study of the functions and mechanisms of ion channels at the molecular scale. However, protein dynamics is essential for the functions of ion channels, which is difficult to infer from static structures. Therefore, computational approaches have been widely applied to explore the dynamic responses and microscopic interactions of macromolecules (Figure 3).
Simulation of Ion Conductance
Since the first ion channel structure was published in 1998 by Doyle et al. in their study of the KcsA channel (Doyle et al., 1998), scientists have been using molecular dynamics (MD) simulations to study the interactions between ion channels and their solvent environment. In 2000, Berneche et al. simulated the KcsA K+ channel in dipalmitoylphosphatidylcholine (DPPC) phospholipid bilayer with KCl aqueous salt solution (Bernèche and Roux, 2000). Their work showed that the translocation of K+ ion along the ion conductance pathway is facilitated by the structural movements of the extracellular residues Glu 71, Asp80, Arg89, and the selectivity filter residues Val76, Gly77. The dynamics of water and K+ ions were reported in several simulation studies, where a single file of potassium with water molecules interspersed in between passes along the channel pore (Guidoni et al., 2000; Shrivastava and Sansom, 2000). The observed interactions with the fluctuating selectivity filter laid the foundation for understanding the basis for K+ selectivity over Na+. Later, Shrivastava et al. showed that K+ ions are preferentially coordinated by eight carbonyl oxygens of the filter (Shrivastava et al., 2002), which mimics the salvation shell of hydrated K+, while the smaller Na+ ion interacts with four carbonyl organs and two water molecules, which exerts larger distortion on the selectivity filter. Further simulation studies revealed that to accommodate the smaller ionic radius of Na+ and provide the expected coordination, the dynamic carbonyl groups need to adapt to a “collapsed” state, which builds up unfavorable strain energy. The local interaction mediated by the filter residues and coordinating waters provides a molecular mechanism for ion selectivity in KcsA (Noskov et al., 2004; Noskov and Roux, 2006; Bostick and Brooks, 2007; Noskov and Roux, 2007). After the initial surge of both experimental and simulation studies of the KcsA channel with the availability of its structural information, research has been continuously growing in studying the structures and ion conductance of potassium channels. Long before the crystal structure of the inward rectifier potassium channel Kir6.2 was solved, MD simulation studies had already been conducted on its homology model, which represents another field of computational applications in pursuing structure-based approaches to research (John and Sali, 2003). Capener et al. showed that the Kir6.2 channel model has a similar ion passing motion with KcsA (Capener et al., 2000).
Despite the success of these early studies, a complete permeation event was not described until Khalili-Araghi et al., who applied an electric field of 1 V across the membrane bilayer for the open state of Kv1.2 channel during 25 ns of MD simulations, and observed permeation events that are consistent with previous studies, where K+ ions pass the channel pore through the water-mediated knock-on mechanism (Khalili-Araghi et al., 2006). With increasing computation power, more detailed studies of the ion conductance mechanisms have been carried out. In 2010, Jensen et al. performed simulations on Kv1.2 with the voltage ranging between −180 and 180 mV over a total of ∼30 μs(Jensen et al., 2010). In this work, not only do they show that the rate-limiting step of K+ conductance is the formation of the knock-on intermediate, where two ions form direct contact but also revealed a then-novel gating mechanism that at reverse or zero voltage, the channel undergoes a dewetting transition into the intrinsically more stable closed state.
In addition to the direct application of electric field across the membrane, computational electrophysiology has also been used to study ion permeation. This approach uses a double membrane to separate the system into two isolated compartments, and different ion concentrations are created in each compartment to directly form an electrochemical gradient across the membrane. Köpfer et al. used this method to study the K+ conductance of the archaeal MthK channel from Methanobacterium thermoautotrophicum and the eukaryotic Kv1.2-Kv2.1 chimeric channel in the physiological voltage range (Köpfer et al., 2014). Their results suggest that direct ion-ion contacts instead of co-translocation of alternating ion and water is the key to highly efficient K+ conductance. Several subsequent studies reported this direct knock-on mechanism for different potassium channels under a variety of simulation conditions (Schewe et al., 2016; Kopec et al., 2018; Kopec et al., 2019; Lolicato et al., 2020). These studies proposed that strong electrostatic repulsion between ions is the main driving force for the fasting permeation of K+ ions, and the higher energetic penalty for complete desolvation of Na+ over K+ contributes to the ion selectivity. The two permeation mechanisms are still under debate. Both experimental and computational research has been providing various support for either of the mechanisms. It cannot be excluded that a mixture of these mechanisms exists depending on the channel and its physiological conditions.
In contrast to K+ channels, the studies of ion conductance of sodium selective channels have been hindered by relatively fewer structural information. The publication of the prokaryotic voltage-gated sodium channel NavAb from Arcobacter butzleri with an activated voltage sensor and a closed pore first provided a structural basis for investigations of the selectivity and transport of Na+ ions (Payandeh et al., 2011). A notable difference between Kv and Nav channels lies in the selectivity filter. The backbone carbonyl oxygens of the Kv channel selectivity filter loop form the cation binding sites to coordinate and pass on K+ ions. However, Nav channel selectivity filters are lined by the side chains, which contribute to the relatively wider and more diversely shaped ion pathways in Nav channels. This difference puzzled scientists as to how the wider filter provides selectivity for the ion with a smaller atomic radius and prevents K+ and Ca2+ to permeate. Corry and Thomas (Corry and Thomas, 2012) calculated the energetics of ion permeation in NavAb via MD simulations, and showed that a plane formed by four of the filter lining glutamate residues prohibits K+ to be rightly coordinated with water molecule bridging between the cation and the carboxylate groups. Similar binding geometry of partially-hydrated Na+ has also been reported in the open state sodium channel from the marine bacterium Magnetococcus sp. (NavMs) (McCusker et al., 2012), showing that the glutamate side chains form the high field strength ion binding site, which has a lower permeation barrier for Na+ over K+ (Ulmschneider et al., 2013). However, the timescale of these studies hampered the observation of concerted movement between ions and proteins. Simulations on the microsecond scale showed that the selectivity filter fluctuates in a coordinated way with ion translocation (Chakrabarti et al., 2013; Boiteux et al., 2014; Ke et al., 2014; Furini and Domene, 2018). With the ion moving across the narrowest part of the pore, the GLU side chains break the symmetric arrangement to coordinate multiple ions to facilitate the passage of the ion into the cavity (Chakrabarti et al., 2013; Boiteux et al., 2014; Guardiani et al., 2017; Callahan and Roux, 2018; Chen et al., 2021). Although simulation studies of bacterial Nav channels have provided valuable insights into the selectivity and ion conductance of Nav channels, structures of eukaryotic Nav channels were not available until 2017 (Shen et al., 2017; Yan et al., 2017; Pan et al., 2018; Shen et al., 2018; Pan et al., 2019; Shen et al., 2019; Jiang et al., 2020). Instead of the homo-tetrameric configuration of the bacterial Nav channels, the four domains of eukaryotic Nav channels are formed by a single polypeptide chain, which gives rise to an asymmetrical pore structure. Compared to the EEEE motif of the prokaryotic SF, eukaryotic SF is constituted by aspartic acid, glutamic acid, lysine, and alanine (DEKA), each from one channel domain. Simulation studies suggested that despite this striking structural difference, the asymmetrical SF functions similarly to the symmetrical prokaryotic counterpart. The ASP and GLU residues coordinate the Na+ ion, and the protonated LYS has been proposed to act as another cation in the knock-on or pass-by mechanism (Zhang J. et al., 2018; Flood et al., 2018).
Analysis of Functional States
Ion channels control the flux of ions by changing between their active and inactive states, which relies on changes between different conformations. Such transitions have proven to be difficult to capture by experimental methods, and even when structures of different states are available, the process of induction and transformation of conformational changes is still hard to explain with the static pictures, which might be important for rational drug design. As a result, computational techniques, especially MD simulations, become useful in studying different function states of ion channels.
As an example, Wen and Zheng performed extensive simulations on the wild-type closed state and the constitutively active mutant of the heat-sensitive cation channel TRPV1 at different temperatures (Wen and Zheng, 2018). They observed a range of hydrogen-bonding network rearrangement between domains in the WT simulations, which are consistent with previous mutational studies, providing a mechanistic explanation of the transition between a closed to pre-open state. Simulations of the gain-of-function mutant complete the process from pre-open to an open state. Similar propagation of movements in the S2-S3 and S4-S5 linkers were observed, suggesting that the dynamic motions in these domains are the key to the channel opening. To study the gating process, Guardiani et al. used targeted molecular dynamics to simulate the transition between closed and open states of the calcium release-activated calcium channel (CRAC) (Guardiani et al., 2021). They reported that upon binding of the activator protein STIM1, a propagation of conformation change occurs from the TM4 to TM1 helices. The extension of TM4 pulls TM3 outward, which in turn moves the lower part of TM1 backward and opens the hydrophobic region of the pore. This result is consistent with patch-clamp experiments on a series of mutants, which showed that disruption of the TM11-TM3 interactions reduces Ca2+ influx (Liu et al., 2019).
For mechanosensitive ion channel TREK-2, Aryal et al. simulated changes in membrane tension by changing the area per lipid, and thus studied the effect of lateral pressures on the conformation of the channel (Aryal et al., 2017). They observed that increase in membrane tension induced a transition from the downstate to the upstate. The expansion of the bilayer increases the cross-sectional area of the lower half of the protein, while the upper half remains unchanged to maintain the integrity of the selectivity filter. Interestingly, such membrane stretch does not induce a conformational change in the homologous non-mechanosensitive TWIK-1 channel, indicating that the membrane tension-dependent conformational change is specific to mechanosensitive K2P channels. NOMPC is also a mechanosensitive ion channel, it has, however, been shown to be activated by compression of the intracellular ankyrin repeat domain on the normal plane of the membrane instead of a stretch along the bilayer. In MD simulations, Wang et al. observed that when a pushing force is applied, the TRP domain undergoes an upward movement and clockwise rotation, which induced a rotation in the S6 helices leading to the opening of the channel (Wang et al., 2021). They validated this using patch-clamp experiments and showed that clear electrical signals could be detected when a compressing force was applied, which were abolished with the addition of a NOMPC blocker. Yet, a pulling force was not able to induce a clear NOMPC-dependent current. Another well-known mechanosensitive channel Piezo1 has also attracted a large amount of research attentions. Due to the local bilayer convex curvature imposed by the bowl-like shape of Piezo1, it has been proposed that as membrane tension increases, the protein flattens, which leads to opening of the channel (Guo and MacKinnon, 2017). Since the membrane is a crowded environment, any effects of alteration of the membrane are likely to propagate beyond a single channel. Jiang et al. used a hyperbolic tangent model to study the membrane topology and channel opening as a result of overlapping neighboring Piezo1 membrane footprints. Together with atomistic MD simulations, they found that this overlap decreased bilayer curvature, creating a tension-free opening of Piezo1 (Jiang W. et al., 2021).
Not only has computational methods showed their broad applicability in studying single protein properties, but they also have been used to investigate interactions between channels and their auxiliary proteins. Catte et al. simulated the human Kv4.3 channel in complex with one of its auxiliary β-subunits, Kv channel-interacting protein 1 (KChIP1), and detailed the structural and energetic changes of the complex upon mutations in the interfaces of the complex (Catte et al., 2019). By combining computational protein-protein docking, MD simulations, and electrophysiology, Kuenze and co-workers showed that the S1, S4, S5 helices of voltage-gated KCNQ1 potassium channel interacts with a three-amino-acid motif (F57-T58-L59) in the auxiliary protein KCNE1, which modulates the channel gating through an allosteric network with S5-S6 pore helices (Kuenze et al., 2020). With more structures of channel-auxiliary protein complexes available, such as the sodium leak channel NALCN in complex with auxiliary subunit FAM155A (Xie et al., 2020) or the AMPA-subtype ionotropic glutamate receptor with the auxiliary subunit γ2, more computational investigations can be conducted to reveal the underlying mechanistic details of such interactions (Yelshanskaya et al., 2022).
Identification of Binding Site and Allosteric Site
Identifying potential binding sites has become a key step in rational drug discovery. Computational methods for identifying protein-ligand binding sites can be roughly classified into sequence-based and structure-based methods. The main assumption for sequence-based approaches is that functionally important binding sites are highly conserved for proteins to maintain their functionalities. Therefore, these methods scan protein sequences for conserved residues as potential binding sites (Capra and Singh, 2007). However, the lack of spatial and physicochemical information in such predictions limits the ability of these methods in providing reliable binding sites. Consequently, structure-based methods have become the mainstream methods for predicting potential binding sites with the advances in obtaining reliable structures both experimentally and computationally. One of the most common approaches in this category is searching for cavities on protein surfaces. Since such methods were pioneered in the 1990s (Levitt and Banaszak, 1992; Laskowski, 1995), developments in this field have not only incorporated both structural and conservation information (Huang and Schroeder, 2006; Capra et al., 2009), but also allowed predictions regarding properties that are highly relevant to drug discovery such as binding affinities (Laurie and Jackson, 2005; Ngan et al., 2012; Liu et al., 2020) and druggability (Le Guilloux et al., 2009; Volkamer et al., 2012; Yuan et al., 2013; Hussein et al., 2015). Methods that infer binding sites from structurally similar template proteins have also been developed (Brylinski and Skolnick, 2008; Wass et al., 2010; Roy and Zhang, 2012; Yang J. et al., 2013). They take advantage of the quick accumulation of protein structures, and provide relatively reliable binding sites for proteins with templates of high structural similarities. In recent years, artificial intelligence methods have also been applied in the prediction of protein-ligand binding sites, which learned the structural or interaction patterns of the binding sites using schemes such as 3D convolutional neural networks (Jiménez et al., 2017) or grid-based approaches (Pu et al., 2019), and may provide satisfactory results when compared with traditional structure-based methods (Capra et al., 2009; Le Guilloux et al., 2009). Additionally, to consider the dynamic nature of ligand binding, MD simulations are also frequently used besides or together with static pocket detection algorithms (Schmidtke et al., 2011; Ung et al., 2016; Chen et al., 2019).
These methods have been widely applied to investigate ligand binding in ion channels. To validate their deep learning-based binding site detection algorithm BiteNet, Igor and Popov applied it to the ATP-gated cation channel P2X3 (Kozlovskii and Popov, 2020). They were able to successfully identify the binding sites for the endogenous agonist ATP and the antagonist AF-219. It is worth noting that such success depends on the protein conformational state. The difference between the agonist-bound and the antagonist-bound conformations prevents the correct prediction of the agonist binding site in the antagonist-bound conformation and vice versa. This is a common limitation of pocket prediction methods using static structures. Nguyen et al. constructed a homology model of human Nav1.5 based on electric eel Nav1.4 (Yan et al., 2017), and used molecular docking to identify binding locations of antiarrhythmic and local anesthetic drugs, lidocaine, flecainide, and ranolazine (Nguyen et al., 2019). Starting from five different initial locations, they were able to locate binding hot spots of each drug, which are consistent with regions suggested by experimental data. Further multi-microsecond MD simulations revealed detailed mechanistic insights into the dynamic passage and binding of the ligands. Faulkner et al. performed atomistic MD simulations to probe the binding sites of the opioid analgesic fentanyl on the GLIC channel (Faulkner et al., 2019). Fentanyl molecules were placed in a bulk solution to allow free exploration of potential binding sites on the protein, and binding affinities of each site were assessed by molecular mechanics Poisson−Boltzmann surface area (MM/PBSA) calculations (Miller et al., 2012). The newly identified stable binding sites in GLIC are different from previously observed sites for other general anesthetics. They also reported that the binding of fentanyl to one of these novel sites leads to conformational changes, which results in the formation of a hydrophobic gate inhibiting ion conductance through the channel.
Targeting allosteric binding sites of ion channels has aroused more and more attention in recent years considering that the emblematic pore-blocking strategy that targets the highly conserved region might suffer from poor subtype specificity (Marzian et al., 2013). The allosteric effect refers to the phenomenon that the modulation of a binding site (named as allosteric site) differed from the classic active site (referred to as the orthosteric site) resulting in the functional change of the protein. A variety of allosteric sites and the corresponding modulators have been discovered in the field of ion channels, according to that organized in the Allosteric Database (Huang et al., 2011; Huang et al., 2014). Using the comprehensive alanine-scanning mutagenesis, patch-clamp electrophysiological recordings as well as molecular docking and molecular dynamic simulation approaches, Decher et al. showed that the Kv1 inhibitor Psora-4 could bind to another less conserved site of the channel in addition to the central pore cavity, which demonstrated a new allosteric site of the Kv1.x channels and thus provided the molecular basis for the development of novel selective voltage-gated channel inhibitors (Marzian et al., 2013). Under the similar strategies, with TKDC, an inhibitor of the TREK subfamily, as the chemical probe, Yang et al. conducted a study combining molecular simulation, mutagenesis, and electrophysiology, and revealed an allosteric site of the two-pore domain potassium (K2P) channel that located in the extracellular cap of the channel (Luo et al., 2017). Accordingly, virtual screening was performed and a series of new inhibitors were identified (Luo et al., 2017). With the combination of molecular dynamics and oocyte electrophysiology studies, Lindahl et al. revealed the special motion and modulation with the modification of a transmembrane binding site within each subunit of the GLIC channel, which thus supports the multisite model of transmembrane allosteric modulation of the channel (Heusser et al., 2018). Using microseconds long atomistic MD simulations, Botello-Smith et al. identified the binding site of the allosteric agonist Yoda1 of Piezo1 (Botello-Smith et al., 2019). Their study indicated that Yoda1 binds at a site approximately 40 A aways from the channel pore, and modulates the channel activity by facilitating force-induced motions of the bound subunit. Similar research was also conducted in various ion channel systems by different groups, including the voltage-gated KCNQ1 potassium channel (Kuenze et al., 2020) and the large-conductance mechanosensitive channel MscL (Kapsalis et al., 2019).
The prediction of potential allosteric sites might also be empowered with the combination of machine learning or deep learning approaches and molecular simulation strategies. The general idea of these algorithms is to analyze the correlation between the orthogonal site and the allosteric site. The prediction could be conducted based on either a single protein structure or combined with trajectories from molecular simulation. In the former case, the correlation might be analyzed based on perturbation exerted on the orthogonal sites using a network model including an elastic network model, anisotropic network model, and graph-based model (Amor et al., 2016; Dokholyan, 2016; Zhang and Nussinov, 2019; Wang J. et al., 2020; Huang et al., 2021; Mersmann et al., 2021; Zheng, 2021; Wu et al., 2022). In the latter case, as described in the cases mentioned above, the results of the simulation-based approaches rely heavily on sampling efficiency. Recently, the reinforced dynamics scheme was reported by Zhang and Wang et al. by constructing the free energy surface with deep neural network models, which turn out to be highly efficient in sampling high-dimensional free energy landscapes (Zhang L. et al., 2018; Wang et al., 2022), and thus might provide new insights into the allosteric modulation of ion channels based on the comprehensive exploration of the protein structural dynamics. Additionally, the evolutionary-based approach could predict the allosteric communication in proteins as well by analyzing the evolutionarily conserved networks of residues (Suel et al., 2003). Case studies also showed that the evolutionary insights could contribute to the understanding of allosteric regulation of the TRP Ion channels (Hilton et al., 2019). However, an in-depth discussion of the principle of the algorithms is beyond the scope of this paper. Details of the algorithms are referred to in the previous literature (Dokholyan, 2016; Liu and Nussinov, 2016; Zhang and Nussinov, 2019).
Conclusion
Ion channels have long been concerned as medically important drug targets. Increasing the ligand screening throughput and ligand binding selectivity are the main challenges in the search for highly potent ion-channel modulators. As the potential key to meeting these challenges, the practice of rational drug design has long been expected.
Resolving the molecular structures of ion channels is getting easier owing to the rapid development of single-particle cryo-EM. Also, the impressive progress in protein structure prediction is bringing fundamental changes to structural biology. Although accurate predictions of multimeric structures, protein-protein interactions and multiple conformations remain challenging, we are already seeing progress in these directions such as the structural prediction of core eukaryotic protein complexes (Humphreys et al., 2021), heterodimeric protein complexes (Bryant et al., 2022) and alternative conformations for receptors and transporters (Del Alamo et al., 2022). Innovations and applications of protein folding algorithms have been accelerated through open-sourcing AlphaFold2-comparable high-performance training code3, and online tools such as ColabFold (Mirdita et al., 2022) and Uni-Fold4. It is expected that protein folding algorithms will soon be used to study ion channel conformational states and channel-modulator interactions, even for new channel discoveries. Results from the exploration in the drug discovery era driven by occasions and classic synthetic chemistry enabled reliable computational models through ligand-based approaches. More abundant and higher quality ion channel structures (obtained either experimentally or computationally), as well as accurate binding affinity predictions and ultra-high-throughput virtual screening, have pushed forward structure-based rational drug design. Besides techniques based on static protein structures, simulation and machine-learning-based methods also allow us to explore the dynamic nature of ion channels and direct the rational drug design process. Various computational studies have revealed the molecular mechanisms of channel functions and potential ligand-binding sites, especially allosteric binding sites. With help from these computational methods, we expect the implementation of a more effective drug discovery paradigm in the coming future (de Oliveira et al., 2021; Lees et al., 2021; Robertson et al., 2022).
Author Contributions
ZZ, ZD and HW conceptualized and finalized the manuscript; ZZ, ZD, QW, YW, DZ, RX, LG, HW drafted the manuscript.
Conflict of Interest
ZD, QW, YW, DZ, RX, LG, HW were employed by the company DP Technology and all authors have participated in the development of the Hermite platform or the Uni-Fold project.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
ZZ thanks Prof. Weinan E and Prof. Jianfeng Pei for the sponsorship of his postdoctoral research.
Footnotes
1https://www.rbvi.ucsf.edu/chimerax/data/stanford-sep2021/alphafold.html
2https://phenix-online.org/documentation/reference/alphafold_in_colab.html
3https://github.com/dptech-corp/Uni-Fold
4Hermite platform: https://hermite.dp.tech
References
Ackerman, M. J., and Clapham, D. E. (1997). Ion Channels-Bbasic Science and Clinical Disease. N. Engl. J. Med. 336, 1575–1586. doi:10.1056/nejm199705293362207
Ahmed, M., Jalily Hasani, H., Ganesan, A., Houghton, M., and Barakat, K. (2017). Modeling the Human Nav1.5 Sodium Channel: Structural and Mechanistic Insights of Ion Permeation and Drug Blockade. Drug Des. Devel Ther. 11, 2301–2324. doi:10.2147/dddt.S133944
Aierken, A., Xie, Y. K., Dong, W., Apaer, A., Lin, J. J., Zhao, Z., et al. (2021). Rational Design of a Modality‐Specific Inhibitor of TRPM8 Channel against Oxaliplatin‐Induced Cold Allodynia. Adv. Sci. 8, 2101717. doi:10.1002/advs.202101717
Ain, Q. U., Aleksandrova, A., Roessler, F. D., and Ballester, P. J. (2015). Machine-learning Scoring Functions to Improve Structure-Based Binding Affinity Prediction and Virtual Screening. Wiley Interdiscip. Rev. Comput. Mol. Sci. 5, 405–424. doi:10.1002/wcms.1225
Alessandri, B., Bättig, K., and Welzl, H. (1989). Effects of Ketamine on Tunnel Maze and Water Maze Performance in the Rat. Behav. Neural Biol. 52, 194–212. doi:10.1016/S0163-1047(89)90313-0
Alexander, S. P. H., Christopoulos, A., Davenport, A. P., Kelly, E., Mathie, A., Peters, J. A., et al. (2019). The Concise Guide to Pharmacology 2019/20: G Protein-Coupled Receptors. Br. J. Pharmacol. 176, S21–S141. doi:10.1111/bph.14748
Alhossary, A., Handoko, S. D., Mu, Y., and Kwoh, C. K. (2015). Fast, Accurate, and Reliable Molecular Docking with QuickVina 2. Bioinformatics 31, 2214–2216. doi:10.1093/bioinformatics/btv082
Alnabati, E., and Kihara, D. (2019). Advances in Structure Modeling Methods for Cryo-Electron Microscopy Maps. Molecules 25, 82. doi:10.3390/molecules25010082
Amor, B. R., Schaub, M. T., Yaliraki, S. N., and Barahona, M. (2016). Prediction of Allosteric Sites and Mediating Interactions through Bond-To-Bond Propensities. Nat. Commun. 7, 12477. doi:10.1038/ncomms12477
Argudo, D., Capponi, S., Bethel, N. P., and Grabe, M. (2019). A Multiscale Model of Mechanotransduction by the Ankyrin Chains of the NOMPC Channel. J. Gen. Physiol. 151, 316–327. doi:10.1085/jgp.201812266
Aronov, A. M. (2005). Predictive In Silico Modeling for hERG Channel Blockers. Drug Discov. Today 10, 149–155. doi:10.1016/s1359-6446(04)03278-7
Aryal, P., Jarerattanachat, V., Clausen, M. V., Schewe, M., McClenaghan, C., Argent, L., et al. (2017). Bilayer-Mediated Structural Transitions Control Mechanosensitivity of the TREK-2 K2P Channel. Structure 25, 708–e2. doi:10.1016/j.str.2017.03.006
Baek, M., DiMaio, F., Anishchenko, I., Dauparas, J., Ovchinnikov, S., Lee, G. R., et al. (2021). Accurate Prediction of Protein Structures and Interactions Using a Three-Track Neural Network. Science 373, 871–876. doi:10.1126/science.abj8754
Bagal, S. K., Brown, A. D., Cox, P. J., Omoto, K., Owen, R. M., Pryde, D. C., et al. (2013). Ion Channels as Therapeutic Targets: A Drug Discovery Perspective. J. Med. Chem. 56, 593–624. doi:10.1021/jm3011433
Baker, M. L., Abeysinghe, S. S., Schuh, S., Coleman, R. A., Abrams, A., Marsh, M. P., et al. (2011). Modeling Protein Structure at Near Atomic Resolutions with Gorgon. J. Struct. Biol. 174, 360–373. doi:10.1016/j.jsb.2011.01.015
Baker, M. R., Rees, I., Ludtke, S. J., Chiu, W., and Baker, M. L. (2012). Constructing and Validating Initial Cα Models from Subnanometer Resolution Density Maps with Pathwalking. Structure 20, 450–463. doi:10.1016/j.str.2012.01.008
Barad, B. A., Echols, N., Wang, R. Y., Cheng, Y., DiMaio, F., Adams, P. D., et al. (2015). EMRinger: Side Chain-Directed Model and Map Validation for 3D Cryo-Electron Microscopy. Nat. Methods 12, 943–946. doi:10.1038/nmeth.3541
Basak, S., Gicheru, Y., Kapoor, A., Mayer, M. L., Filizola, M., and Chakrapani, S. (2019). Molecular Mechanism of Setron-Mediated Inhibition of Full-Length 5-HT3A Receptor. Nat. Commun. 10, 3225. doi:10.1038/s41467-019-11142-8
Bennett, D. L., Clark, A. J., Huang, J., Waxman, S. G., and Dib-Hajj, S. D. (2019). The Role of Voltage-Gated Sodium Channels in Pain Signaling. Physiol. Rev. 99, 1079–1151. doi:10.1152/physrev.00052.2017
Bernèche, S., and Roux, B. (2000). Molecular Dynamics of the KcsA K(+) Channel in a Bilayer Membrane. Biophys. J. 78, 2900–2917. doi:10.1016/s0006-3495(00)76831-7
Bobadilla, J. L., Macek, M., Fine, J. P., and Farrell, P. M. (2002). Cystic Fibrosis: a Worldwide Analysis of CFTR Mutations-Ccorrelation with Incidence Data and Application to Screening. Hum. Mutat. 19, 575–606. doi:10.1002/humu.10041
Boiteux, C., Vorobyov, I., and Allen, T. W. (2014). Ion Conduction and Conformational Flexibility of a Bacterial Voltage-Gated Sodium Channel. Proc. Natl. Acad. Sci. U. S. A. 111, 3454–3459. doi:10.1073/pnas.1320907111
Bostick, D. L., and Brooks, C. L. (2007). Selectivity in K+ Channels Is Due to Topological Control of the Permeant Ion's Coordinated State. Proc. Natl. Acad. Sci. U. S. A. 104, 9260–9265. doi:10.1073/pnas.0700554104
Botello-Smith, W. M., Jiang, W., Zhang, H., Ozkan, A. D., Lin, Y. C., Pham, C. N., et al. (2019). A Mechanism for the Activation of the Mechanosensitive Piezo1 Channel by the Small Molecule Yoda1. Nat. Commun. 10, 4503. doi:10.1038/s41467-019-12501-1
Bourinet, E., Altier, C., Hildebrand, M. E., Trang, T., Salter, M. W., and Zamponi, G. W. (2014). Calcium-Permeable Ion Channels in Pain Signaling. Physiol. Rev. 94, 81–140. doi:10.1152/physrev.00023.2013
Braestrup, C., and Squires, R. F. (1977). Specific Benzodiazepine Receptors in Rat Brain Characterized by High-Affinity (3H)diazepam Binding. Proc. Natl. Acad. Sci. U. S. A. 74, 3805–3809. doi:10.1073/pnas.74.9.3805
Braga, R. C., Alves, V. M., Silva, M. F., Muratov, E., Fourches, D., Lião, L. M., et al. (2015). Pred-hERG: A Novel Web-Accessible Computational Tool for Predicting Cardiac Toxicity. Mol. Inf. 34, 698–701. doi:10.1002/minf.201500040
Brannigan, G., LeBard, D. N., Hénin, J., Eckenhoff, R. G., and Klein, M. L. (2010). Multiple Binding Sites for the General Anesthetic Isoflurane Identified in the Nicotinic Acetylcholine Receptor Transmembrane Domain. Proc. Natl. Acad. Sci. U. S. A. 107, 14122–14127. doi:10.1073/pnas.1008534107
Brugnoni, R., Canioni, E., Filosto, M., Pini, A., Tonin, P., Rossi, T., et al. (2022). Mutations Associated with Hypokalemic Periodic Paralysis: from Hotspot Regions to Complete Analysis of CACNA1S and SCN4A Genes. Neurogenetics 23, 19–25. doi:10.1007/s10048-021-00673-2
Bryant, P., Pozzati, G., and Elofsson, A. (2022). Improved Prediction of Protein-Protein Interactions Using AlphaFold2. Nat. Commun. 13. doi:10.1038/s41467-022-28865-w
Brylinski, M., and Skolnick, J. (2008). A Threading-Based Method (FINDSITE) for Ligand-Binding Site Prediction and Functional Annotation. Proc. Natl. Acad. Sci. U. S. A. 105, 129–134. doi:10.1073/pnas.0707684105
Calcaterra, N. E., and Barrow, J. C. (2014). Classics in Chemical Neuroscience: Diazepam (Valium). ACS Chem. Neurosci. 5, 253–260. doi:10.1021/cn5000056
Callahan, K. M., and Roux, B. (2018). Molecular Dynamics of Ion Conduction through the Selectivity Filter of the NaVAb Sodium Channel. J. Phys. Chem. B 122, 10126–10142. doi:10.1021/acs.jpcb.8b09678
Capener, C. E., Shrivastava, I. H., Ranatunga, K. M., Forrest, L. R., Smith, G. R., and Sansom, M. S. (2000). Homology Modeling and Molecular Dynamics Simulation Studies of an Inward Rectifier Potassium Channel. Biophys. J. 78, 2929–2942. doi:10.1016/s0006-3495(00)76833-0
Capra, J. A., Laskowski, R. A., Thornton, J. M., Singh, M., and Funkhouser, T. A. (2009). Predicting Protein Ligand Binding Sites by Combining Evolutionary Sequence Conservation and 3D Structure. PLoS. Comput. Biol. 5, e1000585. doi:10.1371/journal.pcbi.1000585
Capra, J. A., and Singh, M. (2007). Predicting Functionally Important Residues from Sequence Conservation. Bioinformatics 23, 1875–1882. doi:10.1093/bioinformatics/btm270
Carnevale, V., Delemotte, L., and Howard, R. J. (2021). Molecular Dynamics Simulations of Ion Channels. Trends biochem. Sci. 46, 621–622. doi:10.1016/j.tibs.2021.04.005
Carrieri, A., Muraglia, M., Corbo, F., and Pacifico, C. (2009). 2D- and 3D-QSAR of Tocainide and Mexiletine Analogues Acting as Na(v)1.4 Channel Blockers. Eur. J. Med. Chem. 44, 1477–1485. doi:10.1016/j.ejmech.2008.10.005
Casañal, A., Lohkamp, B., and Emsley, P. (2020). Current Developments in Coot for Macromolecular Model Building of Electron Cryo-Microscopy and Crystallographic Data. Protein Sci. 29, 1069–1078. doi:10.1002/pro.3791
Case, D. A., Cheatham, T. E., Darden, T., Gohlke, H., Luo, R., Merz, K. M., et al. (2005). The Amber Biomolecular Simulation Programs. J. Comput. Chem. 26, 1668–1688. doi:10.1002/jcc.20290
Catalano, C., Ben-Hail, D., Qiu, W., Blount, P., des Georges, A., and Guo, Y. (2021). Cryo-EM Structure of Mechanosensitive Channel YnaI Using SMA2000: Challenges and Opportunities. Membr. (Basel) 11, 849. doi:10.3390/membranes11110849
Catte, A., Ferbel, L., Bhattacharjee, N., Jan Akhunzada, M., D'Agostino, T., and Brancato, G. (2019). In Silico investigation of the Interaction between the Voltage-Gated Potassium Channel Kv4.3 and its Auxiliary Protein KChIP1. Phys. Chem. Chem. Phys. 21, 25290–25301. doi:10.1039/c9cp04082j
Catterall, W. A. (2000). From Ionic Currents to Molecular Mechanisms: the Structure and Function of Voltage-Gated Sodium Channels. Neuron 26, 13–25. doi:10.1016/s0896-6273(00)81133-2
Catterall, W. A. (2010). Ion Channel Voltage Sensors: Structure, Function, and Pathophysiology. Neuron 67, 915–928. doi:10.1016/j.neuron.2010.08.021
Catterall, W. A. (2011). Voltage-gated Calcium Channels. Cold Spring Harb. Perspect. Biol. 3, a003947. doi:10.1101/cshperspect.a003947
Cavalli, A., Poluzzi, E., De Ponti, F., and Recanatini, M. (2002). Toward a Pharmacophore for Drugs Inducing the Long QT Syndrome: Insights from a CoMFA Study of HERG K(+) Channel Blockers. J. Med. Chem. 45, 3844–3853. doi:10.1021/jm0208875
Chacón, P., and Wriggers, W. (2002). Multi-resolution Contour-Based Fitting of Macromolecular Structures. J. Mol. Biol. 317, 375–384. doi:10.1006/jmbi.2002.5438
Chakrabarti, N., Ing, C., Payandeh, J., Zheng, N., Catterall, W. A., and Pomès, R. (2013). Catalysis of Na+ Permeation in the Bacterial Sodium Channel Na(V)Ab. Proc. Natl. Acad. Sci. U. S. A. 110, 11331–11336. doi:10.1073/pnas.1309452110
Chan, K. Y., Gumbart, J., McGreevy, R., Watermeyer, J. M., Sewell, B. T., and Schulten, K. (2011). Symmetry-restrained Flexible Fitting for Symmetric EM Maps. Structure 19, 1211–1218. doi:10.1016/j.str.2011.07.017
Chen, A. Y., Brooks, B. R., and Damjanovic, A. (2021). Determinants of Conductance of a Bacterial Voltage-Gated Sodium Channel. Biophys. J. 120, 3050–3069. doi:10.1016/j.bpj.2021.06.013
Chen, M., Baldwin, P. R., Ludtke, S. J., and Baker, M. L. (2016). De Novo modeling in Cryo-EM Density Maps with Pathwalking. J. Struct. Biol. 196, 289–298. doi:10.1016/j.jsb.2016.06.004
Chen, M., and Ludtke, S. J. (2021). Deep Learning-Based Mixed-Dimensional Gaussian Mixture Model for Characterizing Variability in Cryo-EM. Nat. Methods 18, 930–936. doi:10.1038/s41592-021-01220-5
Chen, Z., Zhang, X., Peng, C., Wang, J., Xu, Z., Chen, K., et al. (2019). D3Pockets: A Method and Web Server for Systematic Analysis of Protein Pocket Dynamics. J. Chem. Inf. Model. 59, 3353–3358. doi:10.1021/acs.jcim.9b00332
Cheng, Y. (2018). Membrane Protein Structural Biology in the Era of Single Particle Cryo-EM. Curr. Opin. Struct. Biol. 52, 58–63. doi:10.1016/j.sbi.2018.08.008
Clapham, D. E. (2003). TRP Channels as Cellular Sensors. Nature 426, 517–524. doi:10.1038/nature02196
Clark, R. J., Jensen, J., Nevin, S. T., Callaghan, B. P., Adams, D. J., and Craik, D. J. (2010). The Engineering of an Orally Active Conotoxin for the Treatment of Neuropathic Pain. Angew. Chem. Int. Ed. Engl. 49, 6545–6548. doi:10.1002/anie.201000620
Cleves, A. E., Johnson, S. R., and Jain, A. N. (2019). Electrostatic-field and Surface-Shape Similarity for Virtual Screening and Pose Prediction. J. Comput. Aided Mol. Des. 33, 865–886. doi:10.1007/s10822-019-00236-6
Cobanoglu, M. C., Liu, C., Hu, F., Oltvai, Z. N., and Bahar, I. (2013). Predicting Drug-Target Interactions Using Probabilistic Matrix Factorization. J. Chem. Inf. Model. 53, 3399–3409. doi:10.1021/ci400219z
Coe, J. W., Brooks, P. R., Vetelino, M. G., Wirtz, M. C., Arnold, E. P., Huang, J., et al. (2005). Varenicline: an Alpha4beta2 Nicotinic Receptor Partial Agonist for Smoking Cessation. J. Med. Chem. 48, 3474–3477. doi:10.1021/jm050069n
Corry, B., and Thomas, M. (2012). Mechanism of Ion Permeation and Selectivity in a Voltage Gated Sodium Channel. J. Am. Chem. Soc. 134, 1840–1846. doi:10.1021/ja210020h
Costa, M. G. S., Fagnen, C., Vénien-Bryan, C., and Perahia, D. (2020). A New Strategy for Atomic Flexible Fitting in Cryo-EM Maps by Molecular Dynamics with Excited Normal Modes (MDeNM-EMfit). J. Chem. Inf. Model. 60, 2419–2423. doi:10.1021/acs.jcim.9b01148
Coste, B., Mathur, J., Schmidt, M., Earley, T. J., Ranade, S., Petrus, M. J., et al. (2010). Piezo1 and Piezo2 Are Essential Components of Distinct Mechanically Activated Cation Channels. Science 330, 55–60. doi:10.1126/science.1193270
Cox, B. (2015). “Chapter 1 Ion Channel Drug Discovery: a Historical Perspective,” in Ion Channel Drug Discovery (Cambridge, United Kingdom: The Royal Society of Chemistry), 1–15.
Croll, T. I. (2018). ISOLDE: a Physically Realistic Environment for Model Building into Low-Resolution Electron-Density Maps. Acta Crystallogr. D. Struct. Biol. 74, 519–530. doi:10.1107/s2059798318002425
Daly, E. J., Singh, J. B., Fedgchin, M., Cooper, K., Lim, P., Shelton, R. C., et al. (2018). Efficacy and Safety of Intranasal Esketamine Adjunctive to Oral Antidepressant Therapy in Treatment-Resistant Depression: A Randomized Clinical Trial. JAMA Psychiatry 75, 139–148. doi:10.1001/jamapsychiatry.2017.3739
Dang, S., van Goor, M. K., Asarnow, D., Wang, Y., Julius, D., Cheng, Y., et al. (2019). Structural Insight into TRPV5 Channel Function and Modulation. Proc. Natl. Acad. Sci. U. S. A. 116, 8869–8878. doi:10.1073/pnas.1820323116
Daniel, N. H., Aravind, A., and Thakur, P. (2021). Are Ion Channels Potential Therapeutic Targets for Parkinson's Disease? Neurotoxicology 87, 243–257. doi:10.1016/j.neuro.2021.10.008
de Oliveira, T. M., van Beek, L., Shilliday, F., Debreczeni, J. É., and Phillips, C. (2021). Cryo-EM: The Resolution Revolution and Drug Discovery. SLAS Discov. 26, 17–31. doi:10.1177/2472555220960401
Del Alamo, D., Sala, D., McHaourab, H. S., and Meiler, J. (2022). Sampling Alternative Conformational States of Transporters and Receptors with AlphaFold2. eLife 11, e75751. doi:10.7554/eLife.75751
Derevyanko, G., and Grudinin, S. (2014). HermiteFit: Fast-Fitting Atomic Structures into a Low-Resolution Density Map Using Three-Dimensional Orthogonal Hermite Functions. Acta Crystallogr. D. Biol. Crystallogr. 70, 2069–2084. doi:10.1107/s1399004714011493
Diver, M. M., Cheng, Y., and Julius, D. (2019). Structural Insights into TRPM8 Inhibition and Desensitization. Science 365, 1434–1440. doi:10.1126/science.aax6672
Dokholyan, N. V. (2016). Controlling Allosteric Networks in Proteins. Chem. Rev. 116, 6463–6487. doi:10.1021/acs.chemrev.5b00544
Doyle, D. A., Morais Cabral, J., Pfuetzner, R. A., Kuo, A., Gulbis, J. M., Cohen, S. L., et al. (1998). The Structure of the Potassium Channel: Molecular Basis of K+ Conduction and Selectivity. Science 280, 69–77. doi:10.1126/science.280.5360.69
Dunlop, J., Bowlby, M., Peri, R., Vasilyev, D., and Arias, R. (2008). High-throughput Electrophysiology: an Emerging Paradigm for Ion-Channel Screening and Physiology. Nat. Rev. Drug Discov. 7, 358–368. doi:10.1038/nrd2552
Durdagi, S., Deshpande, S., Duff, H. J., and Noskov, S. Y. (2012). Modeling of Open, Closed, and Open-Inactivated States of the hERG1 Channel: Structural Mechanisms of the State-dependent Drug Binding. J. Chem. Inf. Model. 52, 2760–2774. doi:10.1021/ci300353u
Dutertre, S., and Lewis, R. J. (2010). Use of Venom Peptides to Probe Ion Channel Structure and Function. J. Biol. Chem. 285, 13315–13320. doi:10.1074/jbc.R109.076596
Ekins, S., Crumb, W. J., Sarazan, R. D., Wikel, J. H., and Wrighton, S. A. (2002). Three-dimensional Quantitative Structure-Activity Relationship for Inhibition of Human Ether-A-Go-Go-Related Gene Potassium Channel. J. Pharmacol. Exp. Ther. 301, 427–434. doi:10.1124/jpet.301.2.427
Escayg, A., MacDonald, B. T., Meisler, M. H., Baulac, S., Huberfeld, G., An-Gourfinkel, I., et al. (2000). Mutations of SCN1A, Encoding a Neuronal Sodium Channel, in Two Families with GEFS+2. Nat. Genet. 24, 343–345. doi:10.1038/74159
Faulkner, C., Plant, D. F., and De Leeuw, N. H. (2019). Modulation of the Gloeobacter Violaceus Ion Channel by Fentanyl: A Molecular Dynamics Study. Biochemistry 58, 4804–4808. doi:10.1021/acs.biochem.9b00881
Flood, E., Boiteux, C., and Allen, T. W. (2018). Selective Ion Permeation Involves Complexation with Carboxylates and Lysine in a Model Human Sodium Channel. PLoS. Comput. Biol. 14, e1006398. doi:10.1371/journal.pcbi.1006398
Flood, E., Boiteux, C., Lev, B., Vorobyov, I., and Allen, T. W. (2019). Atomistic Simulations of Membrane Ion Channel Conduction, Gating, and Modulation. Chem. Rev. 119, 7737–7832. doi:10.1021/acs.chemrev.8b00630
Frank, J. (2017). Advances in the Field of Single-Particle Cryo-Electron Microscopy over the Last Decade. Nat. Protoc. 12, 209–212. doi:10.1038/nprot.2017.004
Frank, J., and Ourmazd, A. (2016). Continuous Changes in Structure Mapped by Manifold Embedding of Single-Particle Data in Cryo-EM. Methods 100, 61–67. doi:10.1016/j.ymeth.2016.02.007
Frenz, B., Walls, A. C., Egelman, E. H., Veesler, D., and DiMaio, F. (2017). RosettaES: a Sampling Strategy Enabling Automated Interpretation of Difficult Cryo-EM Maps. Nat. Methods 14, 797–800. doi:10.1038/nmeth.4340
Furini, S., and Domene, C. (2018). Ion-triggered Selectivity in Bacterial Sodium Channels. Proc. Natl. Acad. Sci. U. S. A. 115, 5450–5455. doi:10.1073/pnas.1722516115
Garzón, J. I., Kovacs, J., Abagyan, R., and Chacón, P. (2007). ADP_EM: Fast Exhaustive Multi-Resolution Docking for High-Throughput Coverage. Bioinformatics 23, 427–433. doi:10.1093/bioinformatics/btl625
Gentile, F., Agrawal, V., Hsing, M., Ton, A. T., Ban, F., Norinder, U., et al. (2020). Deep Docking: A Deep Learning Platform for Augmentation of Structure Based Drug Discovery. ACS Cent. Sci. 6, 939–949. doi:10.1021/acscentsci.0c00229
Gonzales, E. B., and Sumien, N. (2017). Acidity and Acid-Sensing Ion Channels in the Normal and Alzheimer's Disease Brain. J. Alzheimers Dis. 57, 1137–1144. doi:10.3233/jad-161131
Gorgulla, C., Boeszoermenyi, A., Wang, Z. F., Fischer, P. D., Coote, P. W., Padmanabha Das, K. M., et al. (2020). An Open-Source Drug Discovery Platform Enables Ultra-large Virtual Screens. Nature 580, 663–668. doi:10.1038/s41586-020-2117-z
Grant, S. M., and Faulds, D. (1992). Oxcarbazepine. A Review of its Pharmacology and Therapeutic Potential in Epilepsy, Trigeminal Neuralgia and Affective Disorders. Drugs 43, 873–888. doi:10.2165/00003495-199243060-00007
Grant, T., Rohou, A., and Grigorieff, N. (2018). cisTEM, User-Friendly Software for Single-Particle Image Processing. eLife 7, e35383. doi:10.7554/eLife.35383
Guardiani, C., Rodger, P. M., Fedorenko, O. A., Roberts, S. K., and Khovanov, I. A. (2017). Sodium Binding Sites and Permeation Mechanism in the NaChBac Channel: A Molecular Dynamics Study. J. Chem. Theory Comput. 13, 1389–1400. doi:10.1021/acs.jctc.6b01035
Guardiani, C., Sun, D., and Giacomello, A. (2021). Unveiling the Gating Mechanism of CRAC Channel: A Computational Study. Front. Mol. Biosci. 8. doi:10.3389/fmolb.2021.773388
Guedes, I. A., Barreto, A. M. S., Marinho, D., Krempser, E., Kuenemann, M. A., Sperandio, O., et al. (2021). New Machine Learning and Physics-Based Scoring Functions for Drug Discovery. Sci. Rep. 11, 3198. doi:10.1038/s41598-021-82410-1
Guidoni, L., Torre, V., and Carloni, P. (2000). Water and Potassium Dynamics inside the KcsA K(+) Channel. FEBS Lett. 477, 37–42. doi:10.1016/s0014-5793(00)01712-9
Gulbis, J. M., Mann, S., and MacKinnon, R. (1999). Structure of a Voltage-dependent K+ Channel Beta Subunit. Cell 97, 943–952. doi:10.1016/s0092-8674(00)80805-3
Guo, Y. R., and MacKinnon, R. (2017). Structure-based Membrane Dome Mechanism for Piezo Mechanosensitivity. eLife 6, e33660. doi:10.7554/eLife.33660
Gupta, M. K., and Vadde, R. (2021). A Computational Structural Biology Study to Understand the Impact of Mutation on Structure-Function Relationship of Inward-Rectifier Potassium Ion Channel Kir6.2 in Human. J. Biomol. Struct. Dyn. 39, 1447–1460. doi:10.1080/07391102.2020.1733666
Hassan, N. M., Alhossary, A. A., Mu, Y., and Kwoh, C. K. (2017). Protein-Ligand Blind Docking Using QuickVina-W with Inter-process Spatio-Temporal Integration. Sci. Rep. 7, 15451. doi:10.1038/s41598-017-15571-7
He, J., and Huang, S.-Y. (2021). Full-length De Novo Protein Structure Determination from Cryo-EM Maps Using Deep Learning. Bioinformatics 37, 3480–3490. doi:10.1093/bioinformatics/btab357
Heusser, S. A., Lycksell, M., Wang, X., McComas, S. E., Howard, R. J., and Lindahl, E. (2018). Allosteric Potentiation of a Ligand-Gated Ion Channel Is Mediated by Access to a Deep Membrane-Facing Cavity. Proc. Natl. Acad. Sci. U. S. A. 115, 10672–10677. doi:10.1073/pnas.1809650115
Hille, B. (1978). Ionic Channels in Excitable Membranes. Current Problems and Biophysical Approaches. Biophys. J. 22, 283–294. doi:10.1016/s0006-3495(78)85489-7
Hilton, J. K., Kim, M., and Van Horn, W. D. (2019). Structural and Evolutionary Insights Point to Allosteric Regulation of TRP Ion Channels. Acc. Chem. Res. 52, 1643–1652. doi:10.1021/acs.accounts.9b00075
Hodgkin, A. L., and Huxley, A. F. (1952). Currents Carried by Sodium and Potassium Ions through the Membrane of the Giant Axon of Loligo. J. Physiol. 116, 449–472. doi:10.1113/jphysiol.1952.sp004717
Hollingsworth, S. A., and Dror, R. O. (2018). Molecular Dynamics Simulation for All. Neuron 99, 1129–1143. doi:10.1016/j.neuron.2018.08.011
Howard, R. J., Murail, S., Ondricek, K. E., Corringer, P. J., Lindahl, E., Trudell, J. R., et al. (2011). Structural Basis for Alcohol Modulation of a Pentameric Ligand-Gated Ion Channel. Proc. Natl. Acad. Sci. U. S. A. 108, 12149–12154. doi:10.1073/pnas.1104480108
Hryc, C. F., Chen, D. H., Afonine, P. V., Jakana, J., Wang, Z., Haase-Pettingell, C., et al. (2017). Accurate Model Annotation of a Near-Atomic Resolution Cryo-EM Map. Proc. Natl. Acad. Sci. U. S. A. 114, 3103–3108. doi:10.1073/pnas.1621152114
Huang, B., and Schroeder, M. (2006). LIGSITEcsc: Predicting Ligand Binding Sites Using the Connolly Surface and Degree of Conservation. BMC Struct. Biol. 6, 19. doi:10.1186/1472-6807-6-19
Huang, G., Liu, D., Wang, W., Wu, Q., Chen, J., Pan, X., et al. (2022). High-resolution Structures of Human Nav1.7 Reveal Gating Modulation through α-π Helical Transition of S6IV. Cell Rep. 39, 110735. doi:10.1016/j.celrep.2022.110735
Huang, Q., Song, P., Chen, Y., Liu, Z., and Lai, L. (2021). Allosteric Type and Pathways Are Governed by the Forces of Protein-Ligand Binding. J. Phys. Chem. Lett. 12, 5404–5412. doi:10.1021/acs.jpclett.1c01253
Huang, Z., Mou, L., Shen, Q., Lu, S., Li, C., Liu, X., et al. (2014). ASD v2.0: Updated Content and Novel Features Focusing on Allosteric Regulation. Nucleic Acids Res. 42, D510–D516. doi:10.1093/nar/gkt1247
Huang, Z., Zhu, L., Cao, Y., Wu, G., Liu, X., Chen, Y., et al. (2011). ASD: a Comprehensive Database of Allosteric Proteins and Modulators. Nucleic Acids Res. 39, D663–D669. doi:10.1093/nar/gkq1022
Humphreys, I. R., Pei, J., Baek, M., Krishnakumar, A., Anishchenko, I., Ovchinnikov, S., et al. (2021). Computed Structures of Core Eukaryotic Protein Complexes. Science 374, 1340. doi:10.1126/science.abm4805
Hussein, H. A., Borrel, A., Geneix, C., Petitjean, M., Regad, L., and Camproux, A. C. (2015). PockDrug-Server: a New Web Server for Predicting Pocket Druggability on Holo and Apo Proteins. Nucleic Acids Res. 43, W436–W442. doi:10.1093/nar/gkv462
Igaev, M., Kutzner, C., Bock, L. V., Vaiana, A. C., and Grubmüller, H. (2019). Automated Cryo-EM Structure Refinement Using Correlation-Driven Molecular Dynamics. eLife 8. doi:10.7554/eLife.43542
Jakobi, A. J., Wilmanns, M., and Sachse, C. (2017). Model-based Local Density Sharpening of Cryo-EM Maps. eLife 6, e27131. doi:10.7554/eLife.27131
Jensen, M. Ø., Borhani, D. W., Lindorff-Larsen, K., Maragakis, P., Jogini, V., Eastwood, M. P., et al. (2010). Principles of Conduction and Hydrophobic Gating in K+ Channels. Proc. Natl. Acad. Sci. U. S. A. 107, 5833–5838. doi:10.1073/pnas.0911691107
Jiang, D., Gamal El-Din, T., Zheng, N., and Catterall, W. A. (2021a). Expression and Purification of the Cardiac Sodium Channel NaV1.5 for Cryo-EM Structure Determination. Methods Enzymol. 653, 89–101. doi:10.1016/bs.mie.2021.01.030
Jiang, D., Shi, H., Tonggu, L., Gamal El-Din, T. M., Lenaeus, M. J., Zhao, Y., et al. (2020). Structure of the Cardiac Sodium Channel. Cell 180, 122–e10. e110. doi:10.1016/j.cell.2019.11.041
Jiang, W., Del Rosario, J. S., Botello-Smith, W., Zhao, S., Lin, Y. C., Zhang, H., et al. (2021b). Crowding-induced Opening of the Mechanosensitive Piezo1 Channel In Silico. Commun. Biol. 4, 84. doi:10.1038/s42003-020-01600-1
Jiang, Z., Xu, J., Yan, A., and Wang, L. (2021c). A Comprehensive Comparative Assessment of 3D Molecular Similarity Tools in Ligand-Based Virtual Screening. Brief. Bioinform 22. doi:10.1093/bib/bbab231
Jiménez, J., Doerr, S., Martínez-Rosell, G., Rose, A. S., and De Fabritiis, G. (2017). DeepSite: Protein-Binding Site Predictor Using 3D-Convolutional Neural Networks. Bioinformatics 33, 3036–3042. doi:10.1093/bioinformatics/btx350
Jing, Z., Rackers, J. A., Pratt, L. R., Liu, C., Rempe, S. B., and Ren, P. (2021). Thermodynamics of Ion Binding and Occupancy in Potassium Channels. Chem. Sci. 12, 8920–8930. doi:10.1039/d1sc01887f
John, B., and Sali, A. (2003). Comparative Protein Structure Modeling by Iterative Alignment, Model Building and Model Assessment. Nucleic Acids Res. 31, 3982–3992. doi:10.1093/nar/gkg460
Jojoa-Cruz, S., Saotome, K., Murthy, S. E., Tsui, C. C. A., Sansom, M. S., Patapoutian, A., et al. (2018). Cryo-EM Structure of the Mechanically Activated Ion Channel OSCA1.2. eLife 7, e41845. doi:10.7554/eLife.41845
Jumper, J., Evans, R., Pritzel, A., Green, T., Figurnov, M., Ronneberger, O., et al. (2021). Highly Accurate Protein Structure Prediction with AlphaFold. Nature 596, 583–589. doi:10.1038/s41586-021-03819-2
Kapsalis, C., Wang, B., El Mkami, H., Pitt, S. J., Schnell, J. R., Smith, T. K., et al. (2019). Allosteric Activation of an Ion Channel Triggered by Modification of Mechanosensitive Nano-Pockets. Nat. Commun. 10, 4619. doi:10.1038/s41467-019-12591-x
Kaur, S., Gomez-Blanco, J., Khalifa, A. A. Z., Adinarayanan, S., Sanchez-Garcia, R., Wrapp, D., et al. (2021). Local Computational Methods to Improve the Interpretability and Analysis of Cryo-EM Maps. Nat. Commun. 12, 1240. doi:10.1038/s41467-021-21509-5
Ke, S., Timin, E. N., and Stary-Weinzinger, A. (2014). Different Inward and Outward Conduction Mechanisms in NaVMs Suggested by Molecular Dynamics Simulations. PLoS. Comput. Biol. 10, e1003746. doi:10.1371/journal.pcbi.1003746
Keating, M. T. (1996). The Long QT Syndrome. A Review of Recent Molecular Genetic and Physiologic Discoveries. Med. Baltim. 75, 1–5. doi:10.1097/00005792-199601000-00001
Kefauver, J. M., Ward, A. B., and Patapoutian, A. (2020). Discoveries in Structure and Physiology of Mechanically Activated Ion Channels. Nature 587, 567–576. doi:10.1038/s41586-020-2933-1
Keserü, G. M. (2003). Prediction of hERG Potassium Channel Affinity by Traditional and Hologram QSAR Methods. Bioorg. Med. Chem. Lett. 13, 2773–2775. doi:10.1016/s0960-894x(03)00492-x
Khalifa, N., Kumar Konda, L. S., and Kristam, R. (2020). Machine Learning-Based QSAR Models to Predict Sodium Ion Channel (Nav 1.5) Blockers. Future Med. Chem. 12, 1829–1843. doi:10.4155/fmc-2020-0156
Khalili-Araghi, F., Tajkhorshid, E., and Schulten, K. (2006). Dynamics of K+ Ion Conduction through Kv1.2. Biophys. J. 91, L72–L74. doi:10.1529/biophysj.106.091926
Kimanius, D., Dong, L., Sharov, G., Nakane, T., and Scheres, S. H. W. (2021). New Tools for Automated Cryo-EM Single-Particle Analysis in RELION-4.0. Biochem. J. 478, 4169–4185. doi:10.1042/bcj20210708
Kiss, F., Pohóczky, K., Szállási, A., and Helyes, Z. (2020). Transient Receptor Potential (TRP) Channels in Head-And-Neck Squamous Cell Carcinomas: Diagnostic, Prognostic, and Therapeutic Potentials. Int. J. Mol. Sci. 21. doi:10.3390/ijms21176374
Kodandaramaiah, S. B., Franzesi, G. T., Chow, B. Y., Boyden, E. S., and Forest, C. R. (2012). Automated Whole-Cell Patch-Clamp Electrophysiology of Neurons In Vivo. Nat. Methods 9, 585–587. doi:10.1038/nmeth.1993
Konda, L. S. K., Keerthi Praba, S., and Kristam, R. (2019). hERG Liability Classification Models Using Machine Learning Techniques. Comput. Toxicol. 12, 100089. doi:10.1016/j.comtox.2019.100089
Kong, W., Tu, X., Huang, W., Yang, Y., Xie, Z., and Huang, Z. (2020). Prediction and Optimization of NaV1.7 Sodium Channel Inhibitors Based on Machine Learning and Simulated Annealing. J. Chem. Inf. Model. 60, 2739–2753. doi:10.1021/acs.jcim.9b01180
Konze, K. D., Bos, P. H., Dahlgren, M. K., Leswing, K., Tubert-Brohman, I., Bortolato, A., et al. (2019). Reaction-Based Enumeration, Active Learning, and Free Energy Calculations to Rapidly Explore Synthetically Tractable Chemical Space and Optimize Potency of Cyclin-dependent Kinase 2 Inhibitors. J. Chem. Inf. Model. 59, 3782–3793. doi:10.1021/acs.jcim.9b00367
Kopec, W., Köpfer, D. A., Vickery, O. N., Bondarenko, A. S., Jansen, T. L. C., de Groot, B. L., et al. (2018). Direct Knock-On of Desolvated Ions Governs Strict Ion Selectivity in K+ Channels. Nat. Chem. 10, 813–820. doi:10.1038/s41557-018-0105-9
Kopec, W., Rothberg, B. S., and de Groot, B. L. (2019). Molecular Mechanism of a Potassium Channel Gating through Activation Gate-Selectivity Filter Coupling. Nat. Commun. 10, 5366. doi:10.1038/S41467-019-13227-w
Köpfer, D. A., Song, C., Gruene, T., Sheldrick, G. M., Zachariae, U., and De Groot, B. L. (2014). Ion Permeation in K⁺ Channels Occurs by Direct Coulomb Knock-On. Science 346, 352–355. doi:10.1126/science.1254840
Kozlovskii, I., and Popov, P. (2020). Spatiotemporal Identification of Druggable Binding Sites Using Deep Learning. Commun. Biol. 3, 618. doi:10.1038/s42003-020-01350-0
Krivov, G. G., Shapovalov, M. V., and Dunbrack, R. L. (2009). Improved Prediction of Protein Side-Chain Conformations with SCWRL4. Proteins 77, 778–795. doi:10.1002/prot.22488
Kuenze, G., Vanoye, C. G., Desai, R. R., Adusumilli, S., Brewer, K. R., Woods, H., et al. (2020). Allosteric Mechanism for KCNE1 Modulation of KCNQ1 Potassium Channel Activation. eLife 9, e57680. doi:10.7554/eLife.57680
Laskowski, R. A. (1995). SURFNET: a Program for Visualizing Molecular Surfaces, Cavities, and Intermolecular Interactions. J. Mol. Graph. 13, 323–328. doi:10.1016/0263-7855(95)00073-9
Laurie, A. T., and Jackson, R. M. (2005). Q-SiteFinder: an Energy-Based Method for the Prediction of Protein-Ligand Binding Sites. Bioinformatics 21, 1908–1916. doi:10.1093/bioinformatics/bti315
Le Guilloux, V., Schmidtke, P., and Tuffery, P. (2009). Fpocket: An Open Source Platform for Ligand Pocket Detection. BMC Bioinforma. 10, 168. doi:10.1186/1471-2105-10-168
LeBard, D. N., Hénin, J., Eckenhoff, R. G., Klein, M. L., and Brannigan, G. (2012). General Anesthetics Predicted to Block the GLIC Pore with Micromolar Affinity. PLoS. Comput. Biol. 8, e1002532. doi:10.1371/journal.pcbi.1002532
Lees, J. A., Dias, J. M., and Han, S. (2021). Applications of Cryo-EM in Small Molecule and Biologics Drug Design. Biochem. Soc. Trans. 49, 2627–2638. doi:10.1042/bst20210444
Levitt, D. G., and Banaszak, L. J. (1992). POCKET: a Computer Graphics Method for Identifying and Displaying Protein Cavities and Their Surrounding Amino Acids. J. Mol. Graph. 10, 229–234. doi:10.1016/0263-7855(92)80074-N
Li, X., Mooney, P., Zheng, S., Booth, C. R., Braunfeld, M. B., Gubbens, S., et al. (2013). Electron Counting and Beam-Induced Motion Correction Enable Near-Atomic-Resolution Single-Particle Cryo-EM. Nat. Methods 10, 584–590. doi:10.1038/nmeth.2472
Liao, M., Cao, E., Julius, D., and Cheng, Y. (2013). Structure of the TRPV1 Ion Channel Determined by Electron Cryo-Microscopy. Nature 504, 107–112. doi:10.1038/nature12822
Liebschner, D., Afonine, P. V., Baker, M. L., Bunkóczi, G., Chen, V. B., Croll, T. I., et al. (2019). Macromolecular Structure Determination Using X-Rays, Neutrons and Electrons: Recent Developments in Phenix. Acta Crystallogr. D. Struct. Biol. 75, 861–877. doi:10.1107/s2059798319011471
Limongelli, V. (2020). Ligand Binding Free Energy and Kinetics Calculation in 2020. WIREs Comput. Mol. Sci. 10. doi:10.1002/wcms.1455
Lin, S., Han, S., Cai, X., Tan, Q., Zhou, K., Wang, D., et al. (2021). Structures of Gi-Bound Metabotropic Glutamate Receptors mGlu2 and mGlu4. Nature 594, 583–588. doi:10.1038/s41586-021-03495-2
Lindert, S., Alexander, N., Wötzel, N., Karakaş, M., Stewart, P. L., and Meiler, J. (2012a). EM-fold: De Novo Atomic-Detail Protein Structure Determination from Medium-Resolution Density Maps. Structure 20, 464–478. doi:10.1016/j.str.2012.01.023
Lindert, S., Hofmann, T., Wötzel, N., Karakaş, M., Stewart, P. L., and Meiler, J. (2012b). Ab Initio protein Modeling into cryoEM Density Maps Using EM-fold. Biopolymers 97, 669–677. doi:10.1002/bip.22027
Lindert, S., Staritzbichler, R., Wötzel, N., Karakaş, M., Stewart, P. L., and Meiler, J. (2009). EM-fold: De Novo Folding of Alpha-Helical Proteins Guided by Intermediate-Resolution Electron Microscopy Density Maps. Structure 17, 990–1003. doi:10.1016/j.str.2009.06.001
Litan, A., and Langhans, S. A. (2015). Cancer as a Channelopathy: Ion Channels and Pumps in Tumor Development and Progression. Front. Cell. Neurosci. 9, 86. doi:10.3389/fncel.2015.00086
Liu, J., and Nussinov, R. (2016). Allostery: An Overview of its History, Concepts, Methods, and Applications. PLoS. Comput. Biol. 12, e1004966. doi:10.1371/journal.pcbi.1004966
Liu, X., Wu, G., Yu, Y., Chen, X., Ji, R., Lu, J., et al. (2019). Molecular Understanding of Calcium Permeation through the Open Orai Channel. PLoS. Biol. 17, e3000096. doi:10.1371/journal.pbio.3000096
Liu, Y., Grimm, M., Dai, W. T., Hou, M. C., Xiao, Z. X., and Cao, Y. (2020). CB-dock: a Web Server for Cavity Detection-Guided Protein-Ligand Blind Docking. Acta Pharmacol. Sin. 41, 138–144. doi:10.1038/s41401-019-0228-6
Lolicato, M., Natale, A. M., Abderemane-Ali, F., Crottès, D., Capponi, S., Duman, R., et al. (2020). K2P Channel C-type Gating Involves Asymmetric Selectivity Filter Order-Disorder Transitions. Sci. Adv. 6. doi:10.1126/sciadv.abc9174
Lu, R., He, Q., and Wang, J. (2017). “TRPC Channels and Alzheimer's Disease,” in Transient Receptor Potential Canonical Channels and Brain Diseases. Editor Y. Wang, 73–83. doi:10.1007/978-94-024-1088-4_7
Luo, Q., Chen, L., Cheng, X., Ma, Y., Li, X., Zhang, B., et al. (2017). An Allosteric Ligand-Binding Site in the Extracellular Cap of K2P Channels. Nat. Commun. 8, 378. doi:10.1038/s41467-017-00499-3
Lyne, P. D. (2002). Structure-based Virtual Screening: an Overview. Drug Discov. Today 7, 1047–1055. doi:10.1016/s1359-6446(02)02483-2
Lyumkis, D. (2019). Challenges and Opportunities in Cryo-EM Single-Particle Analysis. J. Biol. Chem. 294, 5181–5197. doi:10.1074/jbc.REV118.005602
Maffeo, C., Bhattacharya, S., Yoo, J., Wells, D., and Aksimentiev, A. (2012). Modeling and Simulation of Ion Channels. Chem. Rev. 112, 6250–6284. doi:10.1021/cr3002609
Martin, G. M., Kandasamy, B., DiMaio, F., Yoshioka, C., and Shyng, S. L. (2017). Anti-diabetic Drug Binding Site in a Mammalian KATP Channel Revealed by Cryo-EM. eLife 6, e31054. doi:10.7554/eLife.31054
Marzian, S., Stansfeld, P. J., Rapedius, M., Rinné, S., Nematian-Ardestani, E., Abbruzzese, J. L., et al. (2013). Side Pockets Provide the Basis for a New Mechanism of Kv Channel-specific Inhibition. Nat. Chem. Biol. 9, 507–513. doi:10.1038/nchembio.1271
Masiulis, S., Desai, R., Uchański, T., Serna Martin, I., Laverty, D., Karia, D., et al. (2019). GABAA Receptor Signalling Mechanisms Revealed by Structural Pharmacology. Nature 565, 454–459. doi:10.1038/s41586-018-0832-5
McCusker, E. C., Bagnéris, C., Naylor, C. E., Cole, A. R., D'Avanzo, N., Nichols, C. G., et al. (2012). Structure of a Bacterial Voltage-Gated Sodium Channel Pore Reveals Mechanisms of Opening and Closing. Nat. Commun. 3, 1102. doi:10.1038/ncomms2077
McGrath, N. A., Brichacek, M., and Njardarson, J. T. (2010). A Graphical Journey of Innovative Organic Architectures that Have Improved Our Lives. J. Chem. Educ. 87, 1348–1349. doi:10.1021/ed1003806
McGreevy, R., Teo, I., Singharoy, A., and Schulten, K. (2016). Advances in the Molecular Dynamics Flexible Fitting Method for Cryo-EM Modeling. Methods 100, 50–60. doi:10.1016/j.ymeth.2016.01.009
Menke, J., Maskri, S., and Koch, O. (2021). Computational Ion Channel Research: from the Application of Artificial Intelligence to Molecular Dynamics Simulations. Cell. Physiol. biochem. 55, 14–45. doi:10.33594/000000336
Mersmann, S. F., Strömich, L., Song, F. J., Wu, N., Vianello, F., Barahona, M., et al. (2021). ProteinLens: a Web-Based Application for the Analysis of Allosteric Signalling on Atomistic Graphs of Biomolecules. Nucleic Acids Res. 49, W551–W558. doi:10.1093/nar/gkab350
Miller, B. R., McGee, T. D., Swails, J. M., Homeyer, N., Gohlke, H., and Roitberg, A. E. (2012). MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. J. Chem. Theory Comput. 8, 3314–3321. doi:10.1021/ct300418h
Mirdita, M., Schütze, K., Moriwaki, Y., Heo, L., Ovchinnikov, S., and Steinegger, M. (2022). ColabFold: Making Protein Folding Accessible to All. Nat. Methods. doi:10.1038/s41592-022-01488-1
Mohler, H., and Okada, T. (1977). GABA Receptor Binding with 3H (+) Bicuculline-Methiodide in Rat CNS. Nature 267, 65–67. doi:10.1038/267065a0
Mullard, A. (2020). 2019 FDA Drug Approvals. Nat. Rev. Drug Discov. 19, 79–84. doi:10.1038/d41573-020-00001-7
Muraglia, M., De Bellis, M., Catalano, A., Carocci, A., Franchini, C., Carrieri, A., et al. (2014). N-Aryl-2,6-dimethylbenzamides, a New Generation of Tocainide Analogues as Blockers of Skeletal Muscle Voltage-Gated Sodium Channels. J. Med. Chem. 57, 2589–2600. doi:10.1021/jm401864b
Murail, S., Howard, R. J., Broemstrup, T., Bertaccini, E. J., Harris, R. A., Trudell, J. R., et al. (2012). Molecular Mechanism for the Dual Alcohol Modulation of Cys-Loop Receptors. PLoS. Comput. Biol. 8, e1002710. doi:10.1371/journal.pcbi.1002710
Nakane, T., Kotecha, A., Sente, A., McMullan, G., Masiulis, S., Brown, P. M. G. E., et al. (2020). Single-particle Cryo-EM at Atomic Resolution. Nature 587, 152–156. doi:10.1038/s41586-020-2829-0
Newell, A., Yang, K., and Deng, J. (2016). in Stacked Hourglass Networks for Human Pose Estimation. Editors B. Leibe, J. Matas, N. Sebe, and M. Welling (Cham: Springer International Publishing), 483–499.
Ngan, C. H., Hall, D. R., Zerbe, B., Grove, L. E., Kozakov, D., and Vajda, S. (2012). FTSite: High Accuracy Detection of Ligand Binding Sites on Unbound Protein Structures. Bioinformatics 28, 286–287. doi:10.1093/bioinformatics/btr651
Nguyen, P. T., DeMarco, K. R., Vorobyov, I., Clancy, C. E., and Yarov-Yarovoy, V. (2019). Structural Basis for Antiarrhythmic Drug Interactions with the Human Cardiac Sodium Channel. Proc. Natl. Acad. Sci. U. S. A. 116, 2945–2954. doi:10.1073/pnas.1817446116
Nisius, B., and Göller, A. H. (2009). Similarity-Based Classifier Using Topomers to Provide a Knowledge Base for hERG Channel Inhibition. J. Chem. Inf. Model. 49, 247–256. doi:10.1021/ci800304t
Njardarson, T. G. O. J. T. (2020). Top 200 Brand Name Drugs by Retails Sales in 2020. https://njardarson.lab.arizona.edu/content/top-pharmaceuticals-poster (Accessed March, 2022).
Noskov, S. Y., Bernèche, S., and Roux, B. (2004). Control of Ion Selectivity in Potassium Channels by Electrostatic and Dynamic Properties of Carbonyl Ligands. Nature 431, 830–834. doi:10.1038/nature02943
Noskov, S. Y., and Roux, B. (2007). Importance of Hydration and Dynamics on the Selectivity of the KcsA and NaK Channels. J. Gen. Physiol. 129, 135–143. doi:10.1085/jgp.200609633
Noskov, S. Y., and Roux, B. (2006). Ion Selectivity in Potassium Channels. Biophys. Chem. 124, 279–291. doi:10.1016/j.bpc.2006.05.033
Osterberg, F., and Aqvist, J. (2005). Exploring Blocker Binding to a Homology Model of the Open hERG K+ Channel Using Docking and Molecular Dynamics Methods. FEBS Lett. 579, 2939–2944. doi:10.1016/j.febslet.2005.04.039
Pan, X., Li, Z., Huang, X., Huang, G., Gao, S., Shen, H., et al. (2019). Molecular Basis for Pore Blockade of Human Na+ Channel Nav1.2 by the μ-conotoxin KIIIA. Science 363, 1309–1313. doi:10.1126/science.aaw2999
Pan, X., Li, Z., Zhou, Q., Shen, H., Wu, K., Huang, X., et al. (2018). Structure of the Human Voltage-Gated Sodium Channel Nav1.4 in Complex with β1. Science 362. doi:10.1126/science.aau2486
Papke, R. L., and Heinemann, S. F. (1994). Partial Agonist Properties of Cytisine on Neuronal Nicotinic Receptors Containing the Beta 2 Subunit. Mol. Pharmacol. 45, 142–149.
Payandeh, J., Scheuer, T., Zheng, N., and Catterall, W. A. (2011). The Crystal Structure of a Voltage-Gated Sodium Channel. Nature 475, 353–358. doi:10.1038/nature10238
Pearlstein, R. A., Vaz, R. J., Kang, J., Chen, X. L., Preobrazhenskaya, M., Shchekotikhin, A. E., et al. (2003). Characterization of HERG Potassium Channel Inhibition Using CoMSiA 3D QSAR and Homology Modeling Approaches. Bioorg. Med. Chem. Lett. 13, 1829–1835. doi:10.1016/s0960-894x(03)00196-3
Petrey, D., Xiang, Z., Tang, C. L., Xie, L., Gimpelev, M., Mitros, T., et al. (2003). Using Multiple Structure Alignments, Fast Model Building, and Energetic Analysis in Fold Recognition and Homology Modeling. Proteins 53, 430–435. doi:10.1002/prot.10550
Pfab, J., Phan, N. M., and Si, D. (2021). DeepTracer for Fast De Novo Cryo-EM Protein Structure Modeling and Special Studies on CoV-Related Complexes. Proc. Natl. Acad. Sci. U. S. A. 118. doi:10.1073/pnas.2017525118
Pintilie, G., Chen, D. H., Haase-Pettingell, C. A., King, J. A., and Chiu, W. (2016). Resolution and Probabilistic Models of Components in CryoEM Maps of Mature P22 Bacteriophage. Biophys. J. 110, 827–839. doi:10.1016/j.bpj.2015.11.3522
Pintilie, G., and Chiu, W. (2018). Assessment of Structural Features in Cryo-EM Density Maps Using SSE and Side Chain Z-Scores. J. Struct. Biol. 204, 564–571. doi:10.1016/j.jsb.2018.08.015
Pintilie, G., Zhang, K., Su, Z., Li, S., Schmid, M. F., and Chiu, W. (2020). Measurement of Atom Resolvability in Cryo-EM Maps with Q-Scores. Nat. Methods 17, 328–334. doi:10.1038/s41592-020-0731-1
Pintilie, G. D., Zhang, J., Goddard, T. D., Chiu, W., and Gossard, D. C. (2010). Quantitative Analysis of Cryo-EM Density Map Segmentation by Watershed and Scale-Space Filtering, and Fitting of Structures by Alignment to Regions. J. Struct. Biol. 170, 427–438. doi:10.1016/j.jsb.2010.03.007
Plowman, T. (1982). The Identification of Coca (Erythroxylum Species): 1860-1910. Bot. J. Linn. Soc. 84, 329–353. doi:10.1111/j.1095-8339.1982.tb00368.x
Pu, L., Govindaraj, R. G., Lemoine, J. M., Wu, H. C., and Brylinski, M. (2019). DeepDrug3D: Classification of Ligand-Binding Pockets in Proteins with a Convolutional Neural Network. PLoS. Comput. Biol. 15, e1006718. doi:10.1371/journal.pcbi.1006718
Punjani, A., Rubinstein, J. L., Fleet, D. J., and Brubaker, M. A. (2017). cryoSPARC: Algorithms for Rapid Unsupervised Cryo-EM Structure Determination. Nat. Methods 14, 290–296. doi:10.1038/nmeth.4169
Punjani, A., Zhang, H., and Fleet, D. J. (2020). Non-uniform Refinement: Adaptive Regularization Improves Single-Particle Cryo-EM Reconstruction. Nat. Methods 17, 1214–1221. doi:10.1038/s41592-020-00990-8
Ramírez-Aportela, E., Maluenda, D., Fonseca, Y. C., Conesa, P., Marabini, R., Heymann, J. B., et al. (2021). FSC-Q: a CryoEM Map-To-Atomic Model Quality Validation Based on the Local Fourier Shell Correlation. Nat. Commun. 12, 42. doi:10.1038/s41467-020-20295-w
Ramírez-Aportela, E., Vilas, J. L., Glukhova, A., Melero, R., Conesa, P., Martínez, M., et al. (2020). Automatic Local Resolution-Based Sharpening of Cryo-EM Maps. Bioinformatics 36, 765–772. doi:10.1093/bioinformatics/btz671
Reardon, S. (2018). 'Party Drug' Turned Antidepressant Approaches Approval. Nat. Rev. Drug Discov. 17, 773–775. doi:10.1038/nrd.2018.187
Robertson, M. J., Meyerowitz, J. G., and Skiniotis, G. (2022). Drug Discovery in the Era of Cryo-Electron Microscopy. Trends biochem. Sci. 47, 124–135. doi:10.1016/j.tibs.2021.06.008
Roche, O., Trube, G., Zuegge, J., Pflimlin, P., Alanine, A., and Schneider, G. (2002). A Virtual Screening Method for Prediction of the hERG Potassium Channel Liability of Compound Libraries. ChemBioChem 3, 455–459. doi:10.1002/1439-7633(20020503)3:5<455:AID-CBIC455>3.0.CO;2-L
Ronneberger, O., Fischer, P., and Brox, T. (2015). in U-net: Convolutional Networks for Biomedical Image Segmentation. Editors N. Navab, J. Hornegger, W. M. Wells, and A. F. Frangi (Cham: Springer International Publishing), 234–241.
Rook, M. L., Musgaard, M., and MacLean, D. M. (2021). Coupling Structure with Function in Acid-Sensing Ion Channels: Challenges in Pursuit of Proton Sensors. J. Physiol. 599, 417–430. doi:10.1113/JP278707
Rosenthal, P. B., and Henderson, R. (2003). Optimal Determination of Particle Orientation, Absolute Hand, and Contrast Loss in Single-Particle Electron Cryomicroscopy. J. Mol. Biol. 333, 721–745. doi:10.1016/j.jmb.2003.07.013
Rossmann, M. G., Bernal, R., and Pletnev, S. V. (2001). Combining Electron Microscopic with X-Ray Crystallographic Structures. J. Struct. Biol. 136, 190–200. doi:10.1006/jsbi.2002.4435
Roy, A., and Zhang, Y. (2012). Recognizing Protein-Ligand Binding Sites by Global Structural Alignment and Local Geometry Refinement. Structure 20, 987–997. doi:10.1016/j.str.2012.03.009
Sadybekov, A. A., Sadybekov, A. V., Liu, Y., Iliopoulos-Tsoutsouvas, C., Huang, X. P., Pickett, J., et al. (2022). Synthon-based Ligand Discovery in Virtual Libraries of over 11 Billion Compounds. Nature 601, 452–459. doi:10.1038/s41586-021-04220-9
Saikia, C., Ben-Nissan, G., Reuveny, E., and Karbat, I. (2021). Production of Recombinant Venom Peptides as Tools for Ion Channel Research. Methods Enzymol. 654, 169–201. doi:10.1016/bs.mie.2021.01.029
Sanchez-Garcia, R., Gomez-Blanco, J., Cuervo, A., Carazo, J. M., Sorzano, C. O. S., and Vargas, J. (2021). DeepEMhancer: a Deep Learning Solution for Cryo-EM Volume Post-processing. Commun. Biol. 4, 874. doi:10.1038/s42003-021-02399-1
Sands, Z., Grottesi, A., and Sansom, M. S. (2005). Voltage-gated Ion Channels. Curr. Biol. 15, R44–R47. doi:10.1016/j.cub.2004.12.050
Santos-Martins, D., Solis-Vasquez, L., Tillack, A. F., Sanner, M. F., Koch, A., and Forli, S. (2021). Accelerating AUTODOCK4 with GPUs and Gradient-Based Local Search. J. Chem. Theory Comput. 17, 1060–1073. doi:10.1021/acs.jctc.0c01006
Scheres, S. H. (2012). RELION: Implementation of a Bayesian Approach to Cryo-EM Structure Determination. J. Struct. Biol. 180, 519–530. doi:10.1016/j.jsb.2012.09.006
Schewe, M., Nematian-Ardestani, E., Sun, H., Musinszki, M., Cordeiro, S., Bucci, G., et al. (2016). A Non-canonical Voltage-Sensing Mechanism Controls Gating in K2P K(+) Channels. Cell 164, 937–949. doi:10.1016/j.cell.2016.02.002
Schindler, W., and Häfliger, F. (1954). Über Derivate des Iminodibenzyls. Hca 37, 472–483. doi:10.1002/hlca.19540370211
Schmid, R., and Evans, R. J. (2019). ATP-gated P2X Receptor Channels: Molecular Insights into Functional Roles. Annu. Rev. Physiol. 81, 43–62. doi:10.1146/annurev-physiol-020518-114259
Schmidtke, P., Bidon-Chanal, A., Luque, F. J., and Barril, X. (2011). MDpocket: Open-Source Cavity Detection and Characterization on Molecular Dynamics Trajectories. Bioinformatics 27, 3276–3285. doi:10.1093/bioinformatics/btr550
Sehnal, D., Bittrich, S., Deshpande, M., Svobodová, R., Berka, K., Bazgier, V., et al. (2021). Mol* Viewer: Modern Web App for 3D Visualization and Analysis of Large Biomolecular Structures. Nucleic Acids Res. 49, W431–W437. doi:10.1093/nar/gkab314
Shen, H., Li, Z., Jiang, Y., Pan, X., Wu, J., Cristofori-Armstrong, B., et al. (2018). Structural Basis for the Modulation of Voltage-Gated Sodium Channels by Animal Toxins. Science 362. doi:10.1126/science.aau2596
Shen, H., Liu, D., Wu, K., Lei, J., and Yan, N. (2019). Structures of Human Nav1.7 Channel in Complex with Auxiliary Subunits and Animal Toxins. Science 363, 1303–1308. doi:10.1126/science.aaw2493
Shen, H., Zhou, Q., Pan, X., Li, Z., Wu, J., and Yan, N. (2017). Structure of a Eukaryotic Voltage-Gated Sodium Channel at Near-Atomic Resolution. Science 355. doi:10.1126/science.aal4326
Shidi, T., Ruiqi, C., Mengru, L., Qingde, L., Yanxiang, Z., Ji, D., et al. (2022). Accelerating AutoDock VINA with GPUs. ChemRxiv. doi:10.26434/chemrxiv-2021-3qvn2
Shrivastava, I. H., and Sansom, M. S. (2000). Simulations of Ion Permeation through a Potassium Channel: Molecular Dynamics of KcsA in a Phospholipid Bilayer. Biophys. J. 78, 557–570. doi:10.1016/s0006-3495(00)76616-1
Shrivastava, I. H., Tieleman, D. P., Biggin, P. C., and Sansom, M. S. (2002). K(+) versus Na(+) Ions in a K Channel Selectivity Filter: a Simulation Study. Biophys. J. 83, 633–645. doi:10.1016/s0006-3495(02)75197-7
Shuman, C. R. (1983). Glipizide: an Overview. Am. J. Med. 75, 55–59. doi:10.1016/0002-9343(83)90254-1
Sine, S. M., and Engel, A. G. (2006). Recent Advances in Cys-Loop Receptor Structure and Function. Nature 440, 448–455. doi:10.1038/nature04708
Singer, A., and Sigworth, F. J. (2020). Computational Methods for Single-Particle Electron Cryomicroscopy. Annu. Rev. Biomed. Data Sci. 3, 163–190. doi:10.1146/annurev-biodatasci-021020-093826
Singharoy, A., Teo, I., McGreevy, R., Stone, J. E., Zhao, J., and Schulten, K. (2016). Molecular Dynamics-Based Refinement and Validation for Sub-5 Å Cryo-Electron Microscopy Maps. eLife 5, e16105. doi:10.7554/eLife.16105
Siramshetty, V. B., Chen, Q., Devarakonda, P., and Preissner, R. (2018). The Catch-22 of Predicting hERG Blockade Using Publicly Accessible Bioactivity Data. J. Chem. Inf. Model. 58, 1224–1233. doi:10.1021/acs.jcim.8b00150
Song, K., Wei, M., Guo, W., Quan, L., Kang, Y., Wu, J. X., et al. (2021). Structural Basis for Human TRPC5 Channel Inhibition by Two Distinct Inhibitors. eLife 10. doi:10.7554/eLife.63429
Sorzano, C. O., Marabini, R., Velázquez-Muriel, J., Bilbao-Castro, J. R., Scheres, S. H., Carazo, J. M., et al. (2004). XMIPP: a New Generation of an Open-Source Image Processing Package for Electron Microscopy. J. Struct. Biol. 148, 194–204. doi:10.1016/j.jsb.2004.06.006
Srivastava, R. K., Greff, K., and Schmidhuber, J. (2015). “Training Very Deep Networks,” in Advances in Neural Information Processing Systems (Montreal, Canada: Curran Associates, Inc.), 28.
Süel, G. M., Lockless, S. W., Wall, M. A., and Ranganathan, R. (2003). Evolutionarily Conserved Networks of Residues Mediate Allosteric Communication in Proteins. Nat. Struct. Biol. 10, 59–69. doi:10.1038/nsb881
Suhre, K., Navaza, J., and Sanejouand, Y. H. (2006). NORMA: a Tool for Flexible Fitting of High-Resolution Protein Structures into Low-Resolution Electron-Microscopy-Derived Density Maps. Acta Crystallogr. D. Biol. Crystallogr. 62, 1098–1100. doi:10.1107/s090744490602244x
Sun, J., and MacKinnon, R. (2020). Structural Basis of Human KCNQ1 Modulation and Gating. Cell 180, 340–e9. doi:10.1016/j.cell.2019.12.003
Taboureau, O., and Jørgensen, F. S. (2011). In Silico predictions of hERG Channel Blockers in Drug Discovery: from Ligand-Based and Target-Based Approaches to Systems Chemical Biology. Comb. Chem. High. Throughput Screen 14, 375–387. doi:10.2174/138620711795508322
Tama, F., Miyashita, O., and Brooks, C. L. (2004a). Flexible Multi-Scale Fitting of Atomic Structures into Low-Resolution Electron Density Maps with Elastic Network Normal Mode Analysis. J. Mol. Biol. 337, 985–999. doi:10.1016/j.jmb.2004.01.048
Tama, F., Miyashita, O., and Brooks, C. L. (2004b). Normal Mode Based Flexible Fitting of High-Resolution Structure into Low-Resolution Experimental Data from Cryo-EM. J. Struct. Biol. 147, 315–326. doi:10.1016/j.jsb.2004.03.002
Tan, Y., Deng, Y., and Qing, H. (2012). Calcium Channel Blockers and Alzheimer's Disease. Neural Regen. Res. 7, 137–140. doi:10.3969/j.issn.1673-5374.2012.02.010
Tang, G., Peng, L., Baldwin, P. R., Mann, D. S., Jiang, W., Rees, I., et al. (2007). EMAN2: An Extensible Image Processing Suite for Electron Microscopy. J. Struct. Biol. 157, 38–46. doi:10.1016/j.jsb.2006.05.009
Tegunov, D., and Cramer, P. (2019). Real-time Cryo-Electron Microscopy Data Preprocessing with Warp. Nat. Methods 16, 1146–1152. doi:10.1038/s41592-019-0580-y
Terashi, G., and Kihara, D. (2018a). De Novo main-chain Modeling for EM Maps Using MAINMAST. Nat. Commun. 9, 1618. doi:10.1038/s41467-018-04053-7
Terashi, G., and Kihara, D. (2018b). De Novo main-chain Modeling with MAINMAST in 2015/2016 EM Model Challenge. J. Struct. Biol. 204, 351–359. doi:10.1016/j.jsb.2018.07.013
Terwilliger, T. C., Adams, P. D., Afonine, P. V., and Sobolev, O. V. (2018). A Fully Automatic Method Yielding Initial Models from High-Resolution Cryo-Electron Microscopy Maps. Nat. Methods 15, 905–908. doi:10.1038/s41592-018-0173-1
Terwilliger, T. C., Adams, P. D., Afonine, P. V., and Sobolev, O. V. (2020). Cryo-EM Map Interpretation and Protein Model-Building Using Iterative Map Segmentation. Protein Sci. 29, 87–99. doi:10.1002/pro.3740
Terwilliger, T. C., Poon, B. K., Afonine, P. V., Schlicksup, C. J., Croll, T. I., Millán, C., et al. (2022). Improved AlphaFold Modeling with Implicit Experimental Information. bioRxiv. doi:10.1101/2022.01.07.475350
Tominaga, M., and Caterina, M. J. (2004). Thermosensation and Pain. J. Neurobiol. 61, 3–12. doi:10.1002/neu.20079
Topf, M., Lasker, K., Webb, B., Wolfson, H., Chiu, W., and Sali, A. (2008). Protein Structure Fitting and Refinement Guided by Cryo-EM Density. Structure 16, 295–307. doi:10.1016/j.str.2007.11.016
Trabuco, L. G., Villa, E., Mitra, K., Frank, J., and Schulten, K. (2008). Flexible Fitting of Atomic Structures into Electron Microscopy Maps Using Molecular Dynamics. Structure 16, 673–683. doi:10.1016/j.str.2008.03.005
Traynelis, S. F., Wollmuth, L. P., McBain, C. J., Menniti, F. S., Vance, K. M., Ogden, K. K., et al. (2010). Glutamate Receptor Ion Channels: Structure, Regulation, and Function. Pharmacol. Rev. 62, 405–496. doi:10.1124/pr.109.002451
Ulmschneider, M. B., Bagnéris, C., McCusker, E. C., DeCaen, P. G., Delling, M., Clapham, D. E., et al. (2013). Molecular Dynamics of Ion Transport through the Open Conformation of a Bacterial Voltage-Gated Sodium Channel. Proc. Natl. Acad. Sci. U. S. A. 110, 6364–6369. doi:10.1073/pnas.1214667110
Ung, P. M., Ghanakota, P., Graham, S. E., Lexa, K. W., and Carlson, H. A. (2016). Identifying Binding Hot Spots on Protein Surfaces by Mixed-Solvent Molecular Dynamics: HIV-1 Protease as a Test Case. Biopolymers 105, 21–34. doi:10.1002/bip.22742
Vandenberg, J. I., Perozo, E., and Allen, T. W. (2017). Towards a Structural View of Drug Binding to hERG K+ Channels. Trends Pharmacol. Sci. 38, 899–907. doi:10.1016/j.tips.2017.06.004
Villoutreix, B. O., and Taboureau, O. (2015). Computational Investigations of hERG Channel Blockers: New Insights and Current Predictive Models. Adv. Drug Deliv. Rev. 86, 72–82. doi:10.1016/j.addr.2015.03.003
Volkamer, A., Kuhn, D., Rippmann, F., and Rarey, M. (2012). DoGSiteScorer: a Web Server for Automatic Binding Site Prediction, Analysis and Druggability Assessment. Bioinformatics 28, 2074–2075. doi:10.1093/bioinformatics/bts310
Vuillemot, R., Miyashita, O., Tama, F., Rouiller, I., and Jonic, S. (2022). NMMD: Efficient Cryo-EM Flexible Fitting Based on Simultaneous Normal Mode and Molecular Dynamics Atomic Displacements. J. Mol. Biol. 434, 167483. doi:10.1016/j.jmb.2022.167483
Wang, D., Cui, C., Ding, X., Xiong, Z., Zheng, M., Luo, X., et al. (2019). Improving the Virtual Screening Ability of Target-specific Scoring Functions Using Deep Learning Methods. Front. Pharmacol. 10, 924. doi:10.3389/fphar.2019.00924
Wang, D., Wang, Y., Chang, J., Zhang, L., Wang, H., and E., W. (2022). Efficient Sampling of High-Dimensional Free Energy Landscapes Using Adaptive Reinforced Dynamics. Nat. Comput. Sci. 2, 20–29. doi:10.1038/s43588-021-00173-1
Wang, J., Jain, A., McDonald, L. R., Gambogi, C., Lee, A. L., and Dokholyan, N. V. (2020a). Mapping Allosteric Communications within Individual Proteins. Nat. Commun. 11, 3862. doi:10.1038/s41467-020-17618-2
Wang, Q., Corey, R. A., Hedger, G., Aryal, P., Grieben, M., Nasrallah, C., et al. (2020b). Lipid Interactions of a Ciliary Membrane TRP Channel: Simulation and Structural Studies of Polycystin-2. Structure 28, 169–e5. e165. doi:10.1016/j.str.2019.11.005
Wang, W., and MacKinnon, R. (2017). Cryo-EM Structure of the Open Human Ether-À-Go-Go-Related K+ Channel hERG. Cell 169, 422–e10. doi:10.1016/j.cell.2017.03.048
Wang, Y., Guo, Y., Li, G., Liu, C., Wang, L., Zhang, A., et al. (2021). The Push-To-Open Mechanism of the Tethered Mechanosensitive Ion Channel NompC. eLife 10. doi:10.7554/eLife.58388
Wass, M. N., Kelley, L. A., and Sternberg, M. J. (2010). 3DLigandSite: Predicting Ligand-Binding Sites Using Similar Structures. Nucleic Acids Res. 38, W469–W473. doi:10.1093/nar/gkq406
Wen, H., and Zheng, W. (2018). Decrypting the Heat Activation Mechanism of TRPV1 Channel by Molecular Dynamics Simulation. Biophys. J. 114, 40–52. doi:10.1016/j.bpj.2017.10.034
Wójcikowski, M., Ballester, P. J., and Siedlecki, P. (2017). Performance of Machine-Learning Scoring Functions in Structure-Based Virtual Screening. Sci. Rep. 7, 46710. doi:10.1038/srep46710
Wu, M., and Lander, G. C. (2020). Present and Emerging Methodologies in Cryo-EM Single-Particle Analysis. Biophys. J. 119, 1281–1289. doi:10.1016/j.bpj.2020.08.027
Wu, N., Strömich, L., and Yaliraki, S. N. (2022). Prediction of Allosteric Sites and Signaling: Insights from Benchmarking Datasets. Patterns (N Y) 3, 100408. doi:10.1016/j.patter.2021.100408
Wu, X., Milne, J. L., Borgnia, M. J., Rostapshov, A. V., Subramaniam, S., and Brooks, B. R. (2003). A Core-Weighted Fitting Method for Docking Atomic Structures into Low-Resolution Maps: Application to Cryo-Electron Microscopy. J. Struct. Biol. 141, 63–76. doi:10.1016/s1047-8477(02)00570-1
Wulff, H., Castle, N. A., and Pardo, L. A. (2009). Voltage-gated Potassium Channels as Therapeutic Targets. Nat. Rev. Drug Discov. 8, 982–1001. doi:10.1038/nrd2983
Wulff, H., Christophersen, P., Colussi, P., Chandy, K. G., and Yarov-Yarovoy, V. (2019). Antibodies and Venom Peptides: New Modalities for Ion Channels. Nat. Rev. Drug Discov. 18, 339–357. doi:10.1038/s41573-019-0013-8
Xiang, Z., and Honig, B. (2001). Extending the Accuracy Limits of Prediction for Side-Chain Conformations. J. Mol. Biol. 311, 421–430. doi:10.1006/jmbi.2001.4865
Xie, J., Ke, M., Xu, L., Lin, S., Huang, J., Zhang, J., et al. (2020). Structure of the Human Sodium Leak Channel NALCN in Complex with FAM155A. Nat. Commun. 11, 5831. doi:10.1038/s41467-020-19667-z
Xu, K., Wang, Z., Shi, J., Li, H., and Zhang, Q. C. (2019). A2-Net: Molecular Structure Estimation from Cryo-EM Density Volumes. Aaai 33, 1230–1237. doi:10.1609/aaai.v33i01.33011230
Yan, Z., Zhou, Q., Wang, L., Wu, J., Zhao, Y., Huang, G., et al. (2017). Structure of the Nav1.4-β1 Complex from Electric Eel. Cell 170, 470–e11. e411. doi:10.1016/j.cell.2017.06.039
Yang, J., Roy, A., and Zhang, Y. (2013a). Protein-ligand Binding Site Recognition Using Complementary Binding-specific Substructure Comparison and Sequence Profile Alignment. Bioinformatics 29, 2588–2595. doi:10.1093/bioinformatics/btt447
Yang, Y., Zhang, W., Cheng, J., Tang, Y., Peng, Y., and Li, Z. (2013b). Pharmacophore, 3D-QSAR, and Bayesian Model Analysis for Ligands Binding at the Benzodiazepine Site of GABAA Receptors: the Key Roles of Amino Group and Hydrophobic Sites. Chem. Biol. Drug Des. 81, 583–590. doi:10.1111/cbdd.12100
Yelshanskaya, M. V., Patel, D. S., Kottke, C. M., Kurnikova, M. G., and Sobolevsky, A. I. (2022). Opening of Glutamate Receptor Channel to Subconductance Levels. Nature 605, 172–178. doi:10.1038/s41586-022-04637-w
Yin, Y., Le, S. C., Hsu, A. L., Borgnia, M. J., Yang, H., and Lee, S. Y. (2019). Structural Basis of Cooling Agent and Lipid Sensing by the Cold-Activated TRPM8 Channel. Science 363, 945. doi:10.1126/science.aav9334
Yin, Y., Wu, M., Zubcevic, L., Borschel, W. F., Lander, G. C., and Lee, S. Y. (2018). Structure of the Cold- and Menthol-Sensing Ion Channel TRPM8. Science 359, 237–241. doi:10.1126/science.aan4325
Yip, K. M., Fischer, N., Paknia, E., Chari, A., and Stark, H. (2020). Atomic-resolution Protein Structure Determination by Cryo-EM. Nature 587, 157–161. doi:10.1038/s41586-020-2833-4
Yu, H., Bai, X. C., and Wang, W. (2021). Characterization of the Subunit Composition and Structure of Adult Human glycine Receptors. Neuron 109, 2707–e6. doi:10.1016/j.neuron.2021.08.019
Yuan, Y., Pei, J., and Lai, L. (2013). Binding Site Detection and Druggability Prediction of Protein Targets for Structure-Based Drug Design. Curr. Pharm. Des. 19, 2326–2333. doi:10.2174/1381612811319120019
Zampieri, V., Gobet, A., Robert, X., Falson, P., and Chaptal, V. (2021). CryoEM Reconstructions of Membrane Proteins Solved in Several Amphipathic Solvents, Nanodisc, Amphipol and Detergents, Yield Amphipathic Belts of Similar Sizes Corresponding to a Common Ordered Solvent Layer. Biochim. Biophys. Acta Biomembr. 1863, 183693. doi:10.1016/j.bbamem.2021.183693
Zehni, M., Do, M. N., and Zhao, Z. (2020). “Deepsharpen: Deep-Learning Based Sharpening of 3D Reconstruction Map from Cryo-Electron Microscopy,” in 2020 IEEE 17th International Symposium on Biomedical Imaging Workshops (Iowa City: ISBI Workshops), 1–4.
Zhang, J., Mao, W., Ren, Y., Sun, R. N., Yan, N., and Gong, H. (2018a). Simulating the Ion Permeation and Ion Selection for a Eukaryotic Voltage-Gated Sodium Channel NaVPaS. Protein Cell 9, 580–585. doi:10.1007/s13238-018-0522-y
Zhang, L., Wang, H., and E, W. (2018b). Reinforced Dynamics for Enhanced Sampling in Large Atomic and Molecular Systems. J. Chem. Phys. 148, 124113. doi:10.1063/1.5019675
Zhang, Q., Gao, Z., and Yang, H. (2016). Computer-Aided Drug Discovery and Design Targeting Ion Channels. Curr. Top. Med. Chem. 16, 1819–1829. doi:10.2174/1568026616666160315142346
Zhang, X., Zhang, B., Freddolino, P. L., and Zhang, Y. (2022). CR-I-TASSER: Assemble Protein Structures from Cryo-EM Density Maps Using Deep Convolutional Neural Networks. Nat. Methods 19, 195–204. doi:10.1038/s41592-021-01389-9
Zhang, Y., Ye, F., Zhang, T., Lv, S., Zhou, L., Du, D., et al. (2021). Structural Basis of Ketamine Action on Human NMDA Receptors. Nature 596, 301–305. doi:10.1038/s41586-021-03769-9
Zhao, Y., Huang, G., Wu, J., Wu, Q., Gao, S., Yan, Z., et al. (2019). Molecular Basis for Ligand Modulation of a Mammalian Voltage-Gated Ca2+ Channel. Cell 177, 1495–e12. doi:10.1016/j.cell.2019.04.043
Zheng, W. (2011). Accurate Flexible Fitting of High-Resolution Protein Structures into Cryo-Electron Microscopy Maps Using Coarse-Grained Pseudo-energy Minimization. Biophys. J. 100, 478–488. doi:10.1016/j.bpj.2010.12.3680
Zheng, W. (2021). Predicting Cryptic Ligand Binding Sites Based on Normal Modes Guided Conformational Sampling. Proteins 89, 416–426. doi:10.1002/prot.26027
Zhong, E. D., Bepler, T., Berger, B., and Davis, J. H. (2021a). CryoDRGN: Reconstruction of Heterogeneous Cryo-EM Structures Using Neural Networks. Nat. Methods 18, 176–185. doi:10.1038/s41592-020-01049-4
Zhong, E. D., Bepler, T., Davis, J. H., and Berger, B. 2021b. "Reconstructing Continuous Distributions of 3D Protein Structure from Cryo-EM Images", arXiv. doi:10.48550/arXiv.1909.05215
Keywords: ion channel, cryo-EM, machine learning, molecular dynamics, computer-aided drug design
Citation: Zhu Z, Deng Z, Wang Q, Wang Y, Zhang D, Xu R, Guo L and Wen H (2022) Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design. Front. Pharmacol. 13:939555. doi: 10.3389/fphar.2022.939555
Received: 09 May 2022; Accepted: 07 June 2022;
Published: 28 June 2022.
Edited by:
Chia-Hsueh Lee, St. Jude Children’s Research Hospital, United StatesReviewed by:
Hongtu Zhao, St. Jude Children’s Research Hospital, United StatesJing Huang, Westlake University, China
Copyright © 2022 Zhu, Deng, Wang, Wang, Zhang, Xu, Guo and Wen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Han Wen, d2VuaEBkcC50ZWNo
†These authors have contributed equally to this work and share the first authorship
‡These authors have contributed equally to this work and share senior authorship