Monte Carlo Energy Minimization Research Articles

Molecular Modeling of Cardiac Sodium Channel with Mexiletine.

A sodium channel blocker mexiletine (MEX) is used to treat chronic pain, myotonia and some arrhythmias. Mutations in the pore domain (PD) of voltage-gated sodium channels differently affect tonic block (TB) and use-dependent block (UDB) by MEX. Previous studies identified several MEX-sensing residues in the hNav1.5 channel and demonstrated that the channel block by MEX increases with activation of the voltage-sensing domain III (VSDIII), whereas MEX stabilizes the activated state of VSDIII. Structural rationales for these observations are unclear. Here, Monte Carlo (MC) energy minimizations were used to dock MEX and its more potent analog, Thio-Me2, into the hNav1.5 cryo-EM structure with activated VSDs and presumably inactivated PD. Computations yielded two ensembles of ligand binding poses in close contacts with known MEX-sensing residues in helices S6III, S6IV and P1IV. In both ensembles, the ligand NH3 group approached the cation-attractive site between backbone carbonyls at the outer-pore bottom, while the aromatic ring protruded ether into the inner pore (putative UDB pose) or into the III/IV fenestration (putative TB pose). In silico deactivation of VSDIII shifted helices S4-S5III, S5III, S6III and S6IV and tightened the TB site. In a model with activated VSDIII and three resting VSDs, MC-minimized energy profile of MEX pulled from the TB site towards lipids shows a deep local minimum due to interactions with 11 residues in S5III, P1III, S6III and S6IV. The minimum may correspond to an interim binding site for MEX in the hydrophobic path to the TB site along the lipid-exposed sides of repeats III and IV where 15 polar and aromatic residues would attract cationic blockers. The study explains numerous experimental data and suggests the mechanism of allosteric modification of the MEX binding site by VSDIII.

Pyrethroids in an AlphaFold2 Model of the Insect Sodium Channel.

Simple SummaryPyrethroids are a major class of insecticides for controlling arthropod pests and human disease transmitting vectors. Resistance to pyrethroid insecticides presents a serious obstacle to effective pest control. One major mechanism of pyrethroid resistance is caused by mutations in the voltage-gated sodium channel, the target of pyrethroid toxic action. Previous mutational analyses, coupled with homology modeling, predicted two pyrethroid receptor sites on insect sodium channels. Mutations that are located in or close to these receptor sites likely confer resistance by directly reducing pyrethroid binding affinity. However, many mutations appear to be far from the receptor sites, and how these mutations confer pyrethroid resistance remains unknown. In this study, we employed the AlphaFold2 neural network to generate new models of mosquito and cockroach sodium channels. We then used computational methods to dock pyrethroids in the open and closed states of sodium channels to understand the atomic mechanisms of interactions between pyrethroids and sodium channels. Our analysis suggests that mutations residing beyond the receptor sites alter the action of pyrethroids by allosterically affecting the geometry of the receptor sites. The information gained from this study could be valuable for a structure-based design of novel pyrethroids.Pyrethroid insecticides stabilize the open state of insect sodium channels. Previous mutational, electrophysiological, and computational analyses led to the development of homology models predicting two pyrethroid receptor sites, PyR1 and PyR2. Many of the naturally occurring sodium channel mutations, which confer knockdown resistance (kdr) to pyrethroids, are located within or close to these receptor sites, indicating that these mutations impair pyrethroid binding. However, the mechanism of the state-dependent action of pyrethroids and the mechanisms by which kdr mutations beyond the receptor sites confer resistance remain unclear. Recent advances in protein structure prediction using the AlphaFold2 (AF2) neural network allowed us to generate a new model of the mosquito sodium channel AaNav1-1, with the activated voltage-sensing domains (VSMs) and the presumably inactivated pore domain (PM). We further employed Monte Carlo energy minimizations to open PM and deactivate VSM-I and VSM-II to generate additional models. The docking of a Type II pyrethroid deltamethrin in the models predicted its interactions with many known pyrethroid-sensing residues in the PyR1 and PyR2 sites and revealed ligand-channel interactions that stabilized the open PM and activated VSMs. Our study confirms the predicted two pyrethroid receptor sites, explains the state-dependent action of pyrethroids, and proposes the mechanisms of the allosteric effects of various kdr mutations on pyrethroid action. The AF2-based models may assist in the structure-based design of new insecticides.

Efficient d-dimensional molecular dynamics simulations for studies of the glass-jamming transition.

We develop an algorithm suitable for parallel molecular dynamics simulations in d spatial dimensions and describe its implementation in C++. All routines work in arbitrary d; the maximum simulated d is limited only by available computing resources. These routines include several that are particularly useful for studies of the glass-jamming transition, such as SWAP Monte Carlo and FIRE energy minimization. The scalings of simulation runtimes with the number of particles N and number of simulation threads n_{threads} are comparable to popular molecular dynamics codes such as LAMMPS. The efficient parallel implementation allows simulation of systems that are much larger than those employed in previous high-dimensional glass-transition studies. As a demonstration of the code's capabilities, we show that supercooled d=6 liquids can possess dynamics that are substantially more heterogeneous and experience a breakdown of the Stokes-Einstein relation that is substantially stronger than previously reported, owing at least in part to the much smaller system sizes employed in earlier simulations.

Molecular structure characterization of middle-high rank coal via 13C NMR, XPS and FTIR spectroscopy

Elemental analysis, 13C nuclear magnetic resonance (NMR) spectroscopy, X-ray photoelectron spectroscopy (XPS), Fourier transform infrared spectroscopy (FTIR), combined with Grand Canonical ensemble Monte Carlo (GCMC) simulations and energy minimizations were used to establish the molecular structure model of middle-high coal rank accurately. The three-dimensional molecular structure models of middle-high coal rank (from the No. 6 Mine of Pingdingshan (PDS) Mining Area and the No. 4 Mine of Hebi (HB) Mining Area in Henan Province, Xiegou (XG) Mine, Changping (CP) Mine, and Sihe (SH) Mine in Shanxi Province, China) were optimized and constructed. The results showed that the structural arrangement is more orderly and compact with the increase of coal rank. The established molecular formulae for PDS, HB, XG, CP, and SH are respectively C69H43NO3, C166H93N3O8, C191H125N3O7, C169H105N3O7, and C221H123N3O7.

Possible Mechanism of Ion Selectivity in Eukaryotic Voltage-Gated Sodium Channels

The outer pore of Nav1.x channels is lined by the selectivity-filter ring Asp-Glu-Lys-Ala (DEKA), an outer ring of carboxylates, and two inner rings of backbone carbonyls. A key role of Lys in the Na+/K+ selectivity is known, but the mechanism is unclear. Here, contacts involving DEKA residues in 15 cryo-EM structures of Nav1.x channels were analyzed and Monte Carlo (MC) energy minimizations of models with the DEKA residues in different protonation states, with or without Na+ or K+, were performed. In MC-minimized structures, protonated Lys+ was salt-bridged with Glu, whereas deprotonated Lys•• "dunked" to the inner rings. When Na+ was pulled through the outer pore, it was inevitably chelated by Glu and Lys•• at the narrow pore levels. Lys•• further escorted Na+ to the inner rings and in several steps mutual dispositions of the DEKA residues are similar to those seen in cryo-EM structures. Analogous results were obtained in models with DEKA mutants, which have high, but not low Na+/K+ selectivity. When K+ was pulled through the pore, it was also chelated between Glu and Lys••, but respective distances were larger and K+ energy was higher than in models with Na+. The computations suggest that salt-bridged Lys+ and Glu block the pore. Approaching Na+ would knock out H+, squeeze between Glu and Lys••, and move down escorted by Lys••, whereas the displaced H+ would stay nearby in a H-bond involving Glu or/and Asp. When Na+ leaves the outer pore, reprotonated Lys•• would rejoin Glu to complete the permeation cycle.

Computational Study of Ion Selectivity in Eukaryotic Voltage-Gated Sodium Channels

The outer pore of Nav1.x channels is lined by the selectivity-filter ring Asp-Glu-Lys-Ala (DEKA), an outer ring of carboxylates, and two inner rings of backbone carbonyls. A key role of Lys in the Na/K selectivity is known, but the mechanism is unclear. The DEKA residues may exist in different protonation forms. Here contacts involving DEKA residues in 15 cryo-EM structures of Nav1.x channels were analyzed. The analysis suggests that either Asp or Glu are protonated in the majority of the structures. Monte Carlo (MC) energy minimizations of models with the DEKA residues in different protonation states, with or without Na+ or K+ were further performed. In MC-minimized structures protonated Lys+ was salt bridged with Glu-, whereas deprotonated Lys·· “dunked” to the inner rings. When Na+ was pulled through the outer pore, it was inevitably chelated by Glu- and Lys·· at the narrow pore levels. Lys·· further escorted Na+ to the inner rings and in several steps mutual dispositions of the DEKA residues are similar to those seen in cryo-EM structures. Analogous results were obtained in models with DEKA mutants, which have high, but not low Na/K selectivity. When K+ was pulled through the pore, it was also chelated between Glu- and Lys··, but respective distances were expectedly larger and K+ energy was higher than in models with Na+. The computations suggest that salt-bridged Lys+ and Glu- block the pore. Approaching Na+ would knock-out H+, squeeze between Glu- and Lys·· and move down escorted by Lys··, whereas the displaced H+ would stay nearby in an H-bond involving Glu or/and Asp. When Na+ leaves the outer pore, Lys·· would be re-protonated and rejoin Glu- to complete the permeation cycle. Supported by NSERC (RGPIN-2020-07100) and RSF (17-15-01292-П).

Phenylalkylamines in calcium channels: computational analysis of experimental structures.

Experimental 3D structures of calcium channels with phenylalkylamines (PAAs) provide basis for further analysis of atomic mechanisms of these important cardiovascular drugs. In the crystal structure of the engineered calcium channel CavAb with Br-verapamil and in the cryo-EM structure of the Cav1.1 channel with verapamil, the ligands bind in the inner pore. However, there are significant differences between these structures. In the crystal structure the ligand ammonium group is much closer to the ion in the selectivity-filter region Site 3, which is most proximal to the inner pore, than in the cryo-EM structure. Here we used Monte Carlo energy minimizations to dock PAAs in calcium channels. Our computations suggest that in the crystal structure Site 3 is occupied by a water molecule rather than by a calcium ion. Analysis of the published electron density map does not rule out this possibility. In the cryo-EM structures the ammonium group of verapamil is shifted from the calcium ion in Site 3 either along the pore axis, towards the cytoplasm or away from the axis. Our unbiased docking reproduced these binding modes. However, in the cryo-EM structures detergent and lipid molecules interact with verapamil. When we removed these molecules, the nitrile group of verapamil bound to the calcium ion in Site 3. Models of Cav1.2 with different PAAs suggest similar binding modes and direct contacts of the ligands electronegative atoms with the calcium ion in Site 3. Such interactions explain paradoxes in structure-activity relationships of PAAs.

Identification of New Templates for the Synthesis of BEA, BEC, and ISV Zeolites Using Molecular Topology and Monte Carlo Techniques.

The presence of organic structure directing agents (templates) in the synthesis of zeolites allows the synthesis to be directed, in many cases, toward structures in which there is a large stabilization between the template and the zeolite micropore due to dispersion interactions. Although other factors are also important (temperature, pH, Si/Al ratio, etc.), systems with strong zeolite-template interactions are good candidates for an application of new computational algorithms, for instance those based in molecular topology (MT), that can be used in combination with large databases of organic molecules. Computational design of new templates allows the synthesis of existing and new zeolites to be expanded and refined. Three zeolites with similar 3-D large pore systems, BEA, BEC, and ISV, were selected with the aim of finding new templates for their selective syntheses. Using a training set of active and inactive templates (obtained from the literature) for the synthesis of target zeolites, it was possible to select chemical descriptors related to activity, meaning a good candidate template. With a discriminant function defined upon MT, the screening through a database of organic molecules led to a small subset (preselection) of candidate templates for the synthesis of BEA, BEC, and ISV. As far as we know, this is the first time that topological/topochemical descriptors, which do not consider 3-D information on the molecules, have been used to predict the activity of zeolite structure directing agents (SDAs). Following the prediction of SDAs using MT, an automated approach of sequential template filling of micropores based on a combination of Monte Carlo and lattice energy minimization was applied for all the candidate templates in the three zeolites. Two results can be obtained from this: an evaluation of the quality of the molecular topology QSAR models leading to the preselection of templates, and a final selection of candidate templates for the selective synthesis of BEA, BEC, and ISV. Regarding the latter, a good template will be that which maximizes the zeolite-template dispersion interactions with one, and only one, of the three zeolites. The presented methodology can be used to find alternative (maybe cheaper or perhaps more selective) templates than those already known.

Improving ligand 3D shape similarity-based pose prediction with a continuum solvent model.

In order to improve the pose prediction performance of docking methods, we have previously developed the pose prediction using shape similarity (PoPSS) method. It identifies a ligand conformation of the highest shape similarity with target protein crystal ligands. The identified ligand conformation is then placed into the target protein binding pocket and refined using side-chain repacking and Monte Carlo energy minimization. Subsequently, we have reported a modification to PoPSS, named as PoPSS-Lite, using a simple grid-based energy minimization for side-chain repacking and Tversky correlation coefficient as the similarity metric. This modification has improved the pose prediction performance and PoPSS-Lite was one of the top performers in D3R GC3. Here we report a further modification to PoPSS that utilizes a continuum solvent model to account for water mediated protein ligand interactions. In this approach, named as PoPSS-PB, the ligand conformation of the highest shape similarity with crystal ligands is refined along with the target protein binding site by incorporating the Poisson-Boltzmann electrostatics. The performance of PoPSS-PB along with PoPSS and PoPSS-Lite was prospectively evaluated in D3R GC4. PoPSS-PB not only demonstrated excellent performance with mean and median RMSDs of 1.20 and 1.13Å but also achieved improved performance over PoPSS and PoPSS-Lite. Furthermore, the comparison with other D3R GC4 pose prediction submissions revealed admirable performance. Our results showed that the binding poses of ligands with unknown binding modes can be successfully predicted by utilizing ligand 3D shape similarity with known crystallographic ligands and that taking the solvation into consideration improves pose prediction.

Discovery of a Novel Series of Tricyclic Oxadiazine 4a-Methyl Esters Based on Indoxacarb as Potential Sodium Channel Blocker/Modulator Insecticides.

Indoxacarb, a commercialized oxadiazine insecticide, nearly irreversibly blocks open/inactivated, but not resting sodium channels. The structure-activity relationships showed that the substituents at the position of the chiral atom in the oxadiazine ring are very important to the biological activity of oxadiazine insecticide. Here we synthesized a series of tricyclic oxadiazine 4a-methyl ester derivatives. The chiral atom in the oxadiazine ring has been epimerized and substituted with either pyrethric acid or cinnamic acid derivatives. Benzene ring in the tricyclic moiety was substituted with a chlorine, fluorine, or bromine atom, and nitrogen-linked benzene ring was substituted with a trifluoromethyl or trifluoromethoxy group. Toxicity of these compounds against Spodoptera litura F. was evaluated. Diastereoisomers of most toxic compounds J7 and J9 with pyrethric acid moiety were separated by flash column chromatography. The more polar diastereoisomers, J7-L-Rf and J9-L-Rf, and compounds J24 and J26 with cinnamic acid moiety exhibited highest insecticidal activities. We further used Monte Carlo energy minimizations to dock compound J7 and J24 in the NavMs-based homology model of the open cockroach sodium channel. In the low-energy binding modes, the compound interacted with residues in the inner pore and domain interfaces, which previously were proposed to contribute to receptors of pyrethroids and sodium channel blocker insecticides. Our results define compound J7 and J24 as a potentially useful optimized hit for the development of multiple sites sodium channel blocker or modulator.

Intersegment contacts determine geometry of the open and closed states in P-loop channels

Despite impressive progress in experimental studies of ion channels, determinants of their state-dependent geometry are not completely understood. Previous studies of P-loop channels suggested that the gating mechanism involves coupled movement of the S4-S5 cuff and S6 bundle and emphasized importance of specific intersegment contacts in stabilizing different states. However, it is unclear whether or not such contacts are sufficient to computationally reproduce gating rearrangements. Here we analyzed X-ray and cryo-EM structures of several channels in different functional states and selected structures with the wide cuff (open-state Kv1.2) and narrow cuff (closed-state MlotiK1) for detailed analysis. We revealed three categories of inter-residue contacts within the pore domain: (i) state-dependent and state-independent contacts between helices S4-S5 and S5, which provide integrity and state-dependent dimensions of the S4-S5 cuff; (ii) state-independent contacts between helices S4-S5 and S6 that enable their coupled movement during gating and (iii) state-dependent contacts between S6s that stabilize the open and/or closed activation gate. We imposed these contacts to transform the channels from the open to closed state and vice versa using Monte Carlo energy minimizations. In all cases, the target structures were reached with a good precision. Thus, a limited set of inter-residue contacts can be used to predict computationally state-dependent geometry of the pore domain in P-loop channels. Effects of various engineered and naturally occurring mutations (channelopathies) on the channel gating can be rationalized in view of the contacts.Communicated by Ramaswamy H. Sarma

Atomic Mechanisms of Timothy Syndrome-Associated Mutations in Calcium Channel Cav1.2.

Timothy syndrome (TS) is a very rare multisystem disorder almost exclusively associated with mutations G402S and G406R in helix IS6 of Cav1.2. Recently, mutations R518C/H in helix IIS0 of the voltage sensing domain II (VSD-II) were described as a cause of cardiac-only TS. The three mutations are known to decelerate voltage-dependent inactivation (VDI). Here, we report a case of cardiac-only TS caused by mutation R518C. To explore possible impact of the three mutations on interdomain contacts, we modeled channel Cav1.2 using as templates Class Ia and Class II cryo-EM structures of presumably inactivated channel Cav1.1. In both models, R518 and several other residues in VSD-II donated H-bonds to the IS6-linked α1-interaction domain (AID). We further employed steered Monte Carlo energy minimizations to move helices S4–S5, S5, and S6 from the inactivated-state positions to those seen in the X-ray structures of the open and closed NavAb channel. In the open-state models, positions of AID and VSD-II were similar to those in Cav1.1. In the closed-state models, AID moved along the β subunit (Cavβ) toward the pore axis and shifted AID-bound VSD-II. In all the models R518 retained strong contacts with AID. Our calculations suggest that conformational changes in VSD-II upon its deactivation would shift AID along Cavβ toward the pore axis. The AID-linked IS6 would bend at flexible G402 and G406, facilitating the activation gate closure. Mutations R518C/H weakened the IIS0-AID contacts and would retard the AID shift. Mutations G406R and G402S stabilized the open state and would resist the pore closure. Several Cav1.2 mutations associated with long QT syndromes are consistent with this proposition. Our results provide a mechanistic rationale for the VDI deceleration caused by TS-associated mutations and suggest targets for further studies of calcium channelopathies.

Timothy Syndrome-Associated Mutations Affect State-Dependent Contacts in L-Type Calcium Channel

Timothy syndrome (TS) is a very rare disorder almost exclusively associated with mutations G402S and G406R in IS6 of Cav1.2. Recently, mutations R518C/H in IIS0 of the voltage sensing domain II (VSD-II) were described as a cause of cardiac-only Timothy syndrome. The three mutations decelerate the voltage-dependent inactivation (VDI). We describe a new case of cardiac-only Timothy syndrome caused by mutation R518C. We used Class Ia and Class II cryo-EM structures of presumably inactivated-state channel Cav1.1 to model Cav1.2 and explore possible impact of the three mutations on interdomain contacts. In both models, R518 and several other residues in VSD-II donated H-bonds to the IS6-linked α1-interaction domain (AID). We further employed steered Monte Carlo energy minimizations (MCM) to move helices S4-S5, S5 and S6 to positions seen in the X-ray structures of open and closed sodium channel NavAb. In the open-state models, positions of AID and VSD-II were similar to those in Cav1.1. In the closed-state models, AID moved along the β subunit (Cavβ) towards the pore axis and shifted AID-bound VSD-II. In all the models, R518 retained H-bonds with AID. Our models suggest that VSD-II rearrangement upon its activation would shift AID along Cavβ towards the pore axis. The AID-linked IS6 would bend at flexible G402 and G406, facilitating closure of the activation gate. Mutations R518C/H, which weakened the IIS0-AID contacts in our models, would retard the AID shift. Mutations G406R and G402S stabilized the open state and would resist the pore closure. Several Cav1.2 mutations associated with long QT syndromes are consistent with this scenario. Our results provide a mechanistic rationale for the VDI deceleration caused by TS-associated mutations. They also suggest targets for further studies of calcium channelopathies. Supported by grants to BSZ from RSF and NSERC.

Shape similarity guided pose prediction: lessons from D3R Grand Challenge 3.

To extend the utility of ligand 3D shape similarity into pose prediction and virtual screening, we have previously developed CDVS and PoPSS methods. Both of them utilize ligand 3D shape similarity with the crystallographic ligands to improve pose prediction. While CDVS utilizes shape similarity to select suitable receptor structures for molecular docking, PoPSS places a ligand conformation of the highest shape similarity with crystal ligands into the target protein binding pocket which is then refined by side-chain repacking and Monte Carlo energy minimization. Analyses of PoPSS revealed some drawbacks in ligand conformation generation and the scoring scheme used. Moreover, as PoPSS does not sample the ligand conformation after placing it in the binding pocket, it relies solely on conformation generation methods to produce native like conformations. To address these limitations of PoPSS method, we report here a modified approach named as PoPSS-Lite, where side-chain repacking was replaced by a simple grid-based energy minimization. This modification also allowed the sampling of terminal functional groups while keeping the core scaffold fixed. Furthermore, shape similarity calculations were improved by increasing the number of ligand conformations and using a different similarity metric. The performance of PoPSS-Lite was prospectively evaluated in D3R GC3. Comparison of PoPSS-Lite demonstrated superior performance over PoPSS and CDVS with lower mean and median RMSDs. Furthermore, comparison with other D3R GC3 pose prediction submissions revealed top performance for PoPSS-Lite. Our D3R GC3 result extends our perspective that ligand 3D shape similarity with known crystallographic information can be successfully used to predict the binding pose of ligands with unknown binding modes. Our D3R GC3 results further highlight the necessity for improvement in conformer generation methods in order to improve shape similarity guided pose prediction.

Homology modeling of the transmembrane domain of the GABAA receptor

The permeability of ion channels for ions and substances that bind inside the pore depends on the cross-sectional area of the pore. We have constructed models of the closed, open, and desensitized α1β2γ2 GABAA receptor on the basis of known structures of both prokaryotic and eukaryotic ligand-gated channels. We employed Monte Carlo energy minimization to optimize the model structures. We have found significant pore constrictions, whose diameter depends on the functional state of the receptor in the cytoplasmic, middle, and extracellular parts of the pore-forming M2 segments. It is known that the constrictions in the middle (the 9' ring of residues) and cytoplasmic (the 2' ring of residues) parts of the M2 helices form the activation and desensitization gates of the GABAA receptor. Our computations predict that the constriction in the extracellular part of the M2 helices (the 20' ring of residues) may also serve as a gate in the GABAA receptor, whose physiological role is still unclear. Our results imply that the structures of a number of prokaryotic and eukaryotic ligand-gated channels that have been found in bacteria and lower organisms can be used for homology modeling of the pore region of the human GABAA receptor.

Voltage-Dependent Inhibition of Glycine Receptor Channels by Niflumic Acid

Niflumic acid (NFA) is a member of the fenamate class of nonsteroidal anti-inflammatory drugs. This compound and its derivatives are used worldwide clinically for the relief of chronic and acute pain. NFA is also a commonly used blocker of voltage-gated chloride channels. Here we present evidence that NFA is an efficient blocker of chloride-permeable glycine receptors (GlyRs) with subunit heterogeneity of action. Using the whole-cell configuration of patch-clamp recordings and molecular modeling, we analyzed the action of NFA on homomeric α1ΔIns, α2B, α3L, and heteromeric α1β and α2β GlyRs expressed in CHO cells. NFA inhibited glycine-induced currents in a voltage-dependent manner and its blocking potency in α2 and α3 GlyRs was higher than that in α1 GlyR. The Woodhull analysis suggests that NFA blocks α1 and α2 GlyRs at the fractional electrical distances of 0.16 and 0.65 from the external membrane surface, respectively. Thus, NFA binding site in α1 GlyR is closer to the external part of the membrane, while in α2 GlyR it is significantly deeper in the pore. Mutation G254A at the cytoplasmic part of the α1 GlyR pore-lining TM2 helix (level 2′) increased the NFA blocking potency, while incorporation of the β subunit did not have a significant effect. The Hill plot analysis suggests that α1 and α2 GlyRs are preferably blocked by two and one NFA molecules, respectively. Molecular modeling using Monte Carlo energy minimizations provides the structural rationale for the experimental data and proposes more than one interaction site along the pore where NFA can suppress the ion permeation.

Side chain flexibility and the pore dimensions in the GABAA receptor.

Permeation of ions through open channels and their accessibility to pore-targeting drugs depend on the pore cross-sectional dimensions, which are known only for static X-ray and cryo-EM structures. Here, we have built homology models of the closed, open and desensitized α1β2γ2 GABAA receptor (GABAAR). The models are based, respectively, on the X-ray structure of α3 glycine receptor (α3 GlyR), cryo-EM structure of α1 GlyR and X-ray structure of β3 GABAAR. We employed Monte Carlo energy minimizations to explore how the pore lumen may increase due to repulsions of flexible side chains from a variable-diameter electroneutral atom (an expanding sphere) pulled through the pore. The expanding sphere computations predicted that the pore diameter averaged along the permeation pathway is larger by approximately 3Å than that computed for the models with fixed sidechains. Our models predict three major pore constrictions located at the levels of -2', 9' and 20' residues. Residues around the -2' and 9' rings are known to form the desensitization and activation gates of GABAAR. Our computations predict that the 20' ring may also serve as GABAAR gate whose physiological role is unclear. The side chain flexibility of residues -2', 9' and 20' and hence the dimensions of the constrictions depend on the GABAAR functional state.

Modeling interactions between blocking and permeant cations in the NavMs channel

Mechanisms of sodium channel block by local anesthetics (LAs) are still a matter of intensive studies. In the absence of high-resolution structures of eukaryotic channels, atomic details of LA-channel interactions are analyzed using homology modeling. LAs are predicted to access the closed channel through a sidewalk (fenestration) between the channel repeats, bind in a horizontal orientation, and leave its aromatic moiety in the interface. Recent X-ray structure of a bacterial sodium channel NavMs with a cationic molecule Pl1, which is structurally similar to LAs, has confirmed this theoretical prediction and demonstrated a reduced selectivity filter occupancy by the permeant ions in the Pl1-bound channel. However, the nature of the antagonism between LAs and permeant ions is still unclear. Here we used the NavMs structure and Monte Carlo energy minimizations to model Pl1 binding. Our computations predict that Pl1 can displace permeant ion(s) from the selectivity filter by both steric and electrostatic mechanisms. We hypothesize that the electrostatic mechanism is more general, because it is applicable to many LAs and related drugs, which lack a moiety capable to enter the selectivity filter and sterically displace the permeant ion. The electrostatic mechanism is also consistent with the data that various cationic blockers of potassium channels bind in the inner pore without entering the selectivity filter.

Block of GABA(A) receptor ion channel by penicillin: electrophysiological and modeling insights toward the mechanism.

GABA(A) receptors (GABA(A)R) mainly mediate fast inhibitory neurotransmission in the central nervous system. Different classes of modulators target GABA(A)R properties. Penicillin G (PNG) belongs to the class of noncompetitive antagonists blocking the open GABA(A)R and is a prototype of β-lactam antibiotics. In this study, we combined electrophysiological and modeling approaches to investigate the peculiarities of PNG blockade of GABA-activated currents recorded from isolated rat Purkinje cells and to predict the PNG binding site. Whole-cell patch-сlamp recording and fast application system was used in the electrophysiological experiments. PNG block developed after channel activation and increased with membrane depolarization suggesting that the ligand binds within the open channel pore. PNG blocked stationary component of GABA-activated currents in a concentration-dependent manner with IC50 value of 1.12mM at -70mV. The termination of GABA and PNG co-application was followed by a transient tail current. Protection of the tail current from bicuculline block and dependence of its kinetic parameters on agonist affinity suggest that PNG acts as a sequential open channel blocker that prevents agonist dissociation while the channel remains blocked. We built the GABA(A)R models based on nAChR and GLIC structures and performed an unbiased systematic search of the PNG binding site. Monte-Carlo energy minimization was used to find the lowest energy binding modes. We have shown that PNG binds close to the intracellular vestibule. In both models the maximum contribution to the energy of ligand-receptor interactions revealed residues located on the level of 2', 6' and 9' rings formed by a bundle of M2 transmembrane segments, indicating that these residues most likely participate in PNG binding. The predicted structural models support the described mechanism of PNG block.

Folding similarity of the outer pore region in prokaryotic and eukaryotic sodium channels revealed by docking of conotoxins GIIIA, PIIIA, and KIIIA in a NavAb-based model of Nav1.4.

Voltage-gated sodium channels are targets for many drugs and toxins. However, the rational design of medically relevant channel modulators is hampered by the lack of x-ray structures of eukaryotic channels. Here, we used a homology model based on the x-ray structure of the NavAb prokaryotic sodium channel together with published experimental data to analyze interactions of the μ-conotoxins GIIIA, PIIIA, and KIIIA with the Nav1.4 eukaryotic channel. Using Monte Carlo energy minimizations and published experimentally defined pairwise contacts as distance constraints, we developed a model in which specific contacts between GIIIA and Nav1.4 were readily reproduced without deformation of the channel or toxin backbones. Computed energies of specific interactions between individual residues of GIIIA and the channel correlated with experimental estimates. The predicted complexes of PIIIA and KIIIA with Nav1.4 are consistent with a large body of experimental data. In particular, a model of Nav1.4 interactions with KIIIA and tetrodotoxin (TTX) indicated that TTX can pass between Nav1.4 and channel-bound KIIIA to reach its binding site at the selectivity filter. Our models also allowed us to explain experimental data that currently lack structural interpretations. For instance, consistent with the incomplete block observed with KIIIA and some GIIIA and PIIIA mutants, our computations predict an uninterrupted pathway for sodium ions between the extracellular space and the selectivity filter if at least one of the four outer carboxylates is not bound to the toxin. We found a good correlation between computational and experimental data on complete and incomplete channel block by native and mutant toxins. Thus, our study suggests similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.