γ-Aminobutyric acid type A (GABAA) receptors are ligand-gated ion channels in the central nervous system with largely inhibitory function. Despite being a target for drugs including general anesthetics and benzodiazepines, experimental structures have yet to capture an open state of classical synaptic α1β2γ2 GABAA receptors. Here, we use a goal-oriented adaptive sampling strategy in molecular dynamics simulations followed by Markov state modeling to capture an energetically stable putative open state of the receptor. The model conducts chloride ions with comparable conductance as in electrophysiology measurements. Relative to experimental structures, our open model is relatively expanded at both the cytoplasmic (−2′) and central (9′) gates, coordinated with distinctive rearrangements at the transmembrane αβ subunit interface. Consistent with previous experiments, targeted substitutions disrupting interactions at this interface slowed the open-to-desensitized transition rate. This work demonstrates the capacity of advanced simulation techniques to investigate a computationally and experimentally plausible functionally critical of a complex membrane protein yet to be resolved by experimental methods.

A diverse set of modulators, including stimulants and anesthetics, regulates ion channel function in our nervous system. However, structures of ligand-bound complexes can be difficult to capture by experimental methods, particularly when binding is dynamic. Here, we used computational methods and electrophysiology to identify a possible bound state of a modulatory stimulant derivative in a cryptic vestibular pocket of a mammalian serotonin-3 receptor. We first applied a molecular dynamics simulation–based goal-oriented adaptive sampling method to identify possible open-pocket conformations, followed by Boltzmann docking that combines traditional docking with Markov state modeling. Clustering and analysis of stability and accessibility of docked poses supported a preferred binding site; we further validated this site by mutagenesis and electrophysiology, suggesting a mechanism of potentiation by stabilizing intersubunit contacts. Given the pharmaceutical relevance of serotonin-3 receptors in emesis, psychiatric, and gastrointestinal diseases, characterizing relatively unexplored modulatory sites such as these could open valuable avenues to understanding conformational cycling and designing state-dependent drugs.

In the Big Data era, a change of paradigm in the use of molecular dynamics is required. Trajectories should be stored under FAIR (findable, accessible, interoperable and reusable) requirements to favor its reuse by the community under an open science paradigm.

Background and Purpose: Allosteric modulation of pentameric ligand-gated ion channels (pLGICs) are critical for the action of neurotransmitters and many psychoactive drugs. However, details of their modulatory mechanisms remain unclear, especially beyond the orthosteric neurotransmitter-binding sites. The recently reported prokaryotic symbiont of Tevnia jerichonana ligand-gated ion channel (sTeLIC), a pH-gated homologue of eukaryotic receptors in the pLGIC family, is thought to be modulated by aromatic compounds via a relatively uncharacterised modulatory site in the extracellular vestibule. Experimental Approach: We have characterised the effects of psychostimulant derivatives on sTeLIC using two-electrode voltage-clamp electrophysiology in the presence and absence of engineered mutations, and determined X-ray and cryo-EM structures of the channel in both closed and open states. Key Results: We have shown that sTeLIC is sensitive to potentiation by several amphiphilic compounds, which preferentially bind to a vestibular pocket in the contracted open-state extracellular domain. Conclusions and Implications: This work provides a detailed structure–function mechanism for allosteric potentiation via a noncanonical ligand site, with potential conservation of the eukaryotic pentameric ligand-gated ion channels.

Glycine receptors are pentameric ligand-gated ion channels that conduct chloride ions across postsynaptic membranes to facilitate fast inhibitory neurotransmission. In addition to gating by the glycine agonist, interactions with lipids and other compounds in the surrounding membrane environment modulate their function, but molecular details of these interactions remain unclear, in particular, for cholesterol. Here, we report coarse-grained simulations in a model neuronal membrane for three zebrafish glycine receptor structures representing apparent resting, open, and desensitized states. We then converted the systems to all-atom models to examine detailed lipid interactions. Cholesterol bound to the receptor at an outer-leaflet intersubunit site, with a preference for the open and desensitized versus resting states, indicating that it can bias receptor function. Finally, we used short atomistic simulations and iterative amino acid perturbations to identify residues that may mediate allosteric gating transitions. Frequent cholesterol contacts in atomistic simulations clustered with residues identified by perturbation analysis and overlapped with mutations influencing channel function and pathology. Cholesterol binding at this site was also observed in a recently reported pig heteromeric glycine receptor. These results indicate state-dependent lipid interactions relevant to allosteric transitions of glycine receptors, including specific amino acid contacts applicable to biophysical modeling and pharmaceutical design.

Loss-of-function mutations of the CFTR gene cause the life-shortening genetic disease cystic fibrosis (CF), whereas overactivity of CFTR may lead to secretory diarrhea and polycystic kidney disease. While effective drugs targeting the CFTR protein have been developed for the treatment of CF, little progress has been made for diseases caused by hyper-activated CFTR. Here, we solve the cryo-EM structure of CFTR in complex with CFTRinh-172 (Inh-172), a CFTR gating inhibitor with promising potency and efficacy. We find that Inh-172 binds inside the pore of CFTR, interacting with amino acid residues from transmembrane segments (TMs) 1, 6, 8, 9, and 12 through mostly hydrophobic interactions and a salt bridge. Substitution of these residues lowers the apparent affinity of Inh-172. The inhibitor-bound structure reveals re-orientations of the extracellular segment of TMs 1, 8, and 12, supporting an allosteric modulation mechanism involving post-binding conformational changes. This allosteric inhibitory mechanism readily explains our observations that pig CFTR, which preserves all the amino acid residues involved in Inh-172 binding, exhibits a much-reduced sensitivity to Inh-172 and that the apparent affinity of Inh-172 is altered by the CF drug ivacaftor (i.e., VX-770) which enhances CFTR’s activity through binding to a site also comprising TM8.

ρ-type γ-aminobutyric acid-A (GABAA) receptors are widely distributed in the retina and brain, and are potential drug targets for the treatment of visual, sleep and cognitive disorders. Endogenous neuroactive steroids including β-estradiol and pregnenolone sulfate negatively modulate the function of ρ1 GABAA receptors, but their inhibitory mechanisms are not clear. By combining five cryo-EM structures with electrophysiology and molecular dynamics simulations, we characterize binding sites and negative modulation mechanisms of β-estradiol and pregnenolone sulfate at the human ρ1 GABAA receptor. β-estradiol binds in a pocket at the interface between extracellular and transmembrane domains, apparently specific to the ρ subfamily, and disturbs allosteric conformational transitions linking GABA binding to pore opening. In contrast, pregnenolone sulfate binds inside the pore to block ion permeation, with a preference for activated structures. These results illuminate contrasting mechanisms of ρ1 inhibition by two different neuroactive steroids, with potential implications for subtype-specific gating and pharmacological design.

Transient receptor potential (TRP) ion channels are involved in the surveillance or regulation of the acid-base balance. Here, we demonstrate that weak carbonic acids, including acetic acid, lactic acid, and CO₂ activate and sensitize TRPV2 through a mechanism requiring permeation through the cell membrane. TRPV2 channels in cell-free inside-out patches maintain weak acid-sensitivity, but protons applied on either side of the membrane do not induce channel activation or sensitization. The involvement of proton modulation sites for weak acid-sensitivity was supported by the identification of titratable extracellular (Glu495, Glu561) and intracellular (His521) residues on a cryo-EM structure of rat TRPV2 (rTRPV2) treated with acetic acid. Molecular dynamics simulations as well as patch clamp experiments on mutant rTRPV2 constructs confirmed that these residues are critical for weak acid-sensitivity. We also demonstrate that the pore residue Glu609 dictates an inhibition of weak acid-induced currents by extracellular calcium. Finally, TRPV2-expression in HEK293 cells is associated with an increased weak acid-induced cytotoxicity. Together, our data provide new insights into weak acids as endogenous modulators of TRPV2.

We developed fluorescent protein sensors for nicotine with improved sensitivity. For iNicSnFR12 at pH 7.4, the proportionality constant for ∆F/F₀vs [nicotine] (δ-slope, 2.7 μM⁻¹) is 6.1-fold higher than the previously reported iNicSnFR3a. The activated state of iNicSnFR12 has a fluorescence quantum yield of at least 0.6. We measured similar dose-response relations for the nicotine-induced absorbance increase and fluorescence increase, suggesting that the absorbance increase leads to the fluorescence increase via the previously described nicotine-induced conformational change, the ‘candle snuffer’ mechanism. Molecular dynamics (MD) simulations identified a binding pose for nicotine, previously indeterminate from experimental data. MD simulations also showed that Helix 4 of the periplasmic binding protein (PBP) domain appears tilted in iNicSnFR12 relative to iNicSnFR3a, likely altering allosteric network(s) that link the ligand binding site to the fluorophore. In thermal melt experiments, nicotine stabilized the PBP of the tested iNicSnFR variants. iNicSnFR12 resolved nicotine in diluted mouse and human serum at 100 nM, the peak [nicotine] that occurs during smoking or vaping, and possibly at the decreasing levels during intervals between sessions. NicSnFR12 was also partially activated by unidentified endogenous ligand(s) in biofluids. Improved iNicSnFR12 variants could become the molecular sensors in continuous nicotine monitors for animal and human biofluids.

It is commonly assumed that ionizable molecules, such as drugs, permeate through the skin barrier in their neutral form. By using molecular dynamics simulations of the charged and neutral states separately, we can study the dynamic protonation behavior during the permeation process. We have studied three weak acids and three weak bases and conclude that the acids are ionized to a larger extent than the bases, when passing through the headgroup region of the lipid barrier structure, at pH values close to their pKa. It can also be observed that even if these dynamic protonation simulations are informative, in the cases studied herein they are not necessary for the calculation of permeability coefficients. It is sufficient to base the calculations only on the neutral form, as is commonly done.