
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

51904
10.1038/s41467-024-51904-7
Article
Divergent mechanisms of steroid inhibition in the human ρ1 GABAA receptor
http://orcid.org/0000-0003-1211-3427
Fan Chen 12
http://orcid.org/0000-0002-7968-8359
Cowgill John 2
http://orcid.org/0000-0003-2049-3378
Howard Rebecca J. rebecca.howard@dbb.su.se

12
Lindahl Erik erik.lindahl@dbb.su.se

12
1 grid.5037.1 0000000121581746 Dept. of Applied Physics, Science for Life Laboratory, KTH Royal Institute of Technology, Solna, Sweden
2 grid.10548.38 0000 0004 1936 9377 Dept. of Biochemistry and Biophysics, Science for Life Laboratory, Stockholm University, Solna, Sweden
6 9 2024
6 9 2024
2024
15 779512 4 2024
19 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
ρ-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.

Neurological processes from vision to cognition rely on precise control of ion channels, particularly GABAA receptors. Here, the authors report structures of a ρ1 GABAA receptor with naturally occurring steroids, revealing their specific interactions and suggesting approaches to understand and develop drugs.

Subject terms

Cryoelectron microscopy
Ion transport
Ligand-gated ion channels
https://doi.org/10.13039/100004410 European Molecular Biology Organization (EMBO) https://doi.org/10.13039/501100004359 Vetenskapsrådet (Swedish Research Council) 2019-02433, 2021-05806 Lindahl Erik https://doi.org/10.13039/501100009244 Stockholms Universitet (Stockholm University) FV-5.1.2-0523-19 Lindahl Erik issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The neurotransmitter-gated γ-aminobutyric acid-A (GABAA) receptors are anion-permeable pentameric ligand-gated ion channels expressed throughout the nervous system and other tissues. In response to binding the neurotransmitter GABA at an orthosteric site in the extracellular domain (ECD), a series of allosteric conformational changes open a pore over 50 Å away in the transmembrane domain (TMD), allowing anions (typically chloride) to transit the lipid bilayer1. In the continued presence of GABA, this activated open state typically transitions to a more thermodynamically stable desensitized state, with ion permeation occluded at the inner mouth of the TMD pore. Each subunit of the ECD contains 10 strands (β1–β10) interspersed by loops, some of which contribute to agonist binding; each subunit of the TMD contains 4 helices (M1–M4), with the M2 helices surrounding the central pore. In humans, GABAA receptors are homo- or hetero-pentamers formed from a selection of 19 different subunits (α1-6, β1-3, γ1-3, ρ1-3, δ, ε, π and θ).

Although the ρ subtype is similar in sequence and structure to other GABAA receptors, it was previously named GABAC due to its distinct physiological and pharmacological properties2. These include insensitivity to bicuculline and sensitivity to the ρ-type specific inhibitor (1,2,5,6-tetrahydropyridin-4-yl)methylphosphinic acid (TPMPA). Of the three ρ GABAA-receptor subtypes found in mammals, ρ1 is located predominantly in the retina; ρ2 is more widely distributed in the brain, including the cerebellum, thalamus, and frontal cortices; and ρ3 is found in the hippocampus and cerebellum3–5. These channels play roles during earlier postnatal neurodevelopment5 and as potential therapeutic targets for post-stroke motor recovery6. There is increasing interest in developing drugs specific to ρ-type GABAA receptors7. To better understand this system, we recently reported electron cryomicroscopy (cryo-EM) structures of the human ρ1 GABAA receptor (henceforth termed ρ1) to 2.3 Å resolution in the absence and presence of classic agonists and inhibitors8. These structures were facilitated by deleting the flexible N-terminal region and intracellular M3-M4 loop from the wild-type sequence, generating the modified construct ρ1-EM. These modifications improve experimental accessibility while preserving wild-type functional features, enabling opportunities to characterize binding and modulation by pharmacologically relevant compounds.

Interestingly, a number of endogenous neuroactive steroids have been found to modulate GABAA receptors, including ρ19. A site for steroid potentiation at the transmembrane subunit interface, facing the inner membrane leaflet, has been described in some detail; notably, allopregnanolone, which is synthesized from progesterone locally in the brain, was recently resolved by cryo-EM at this site between β and α subunits in classical synaptic α1β2ɣ2 GABAA receptors10,11. The therapeutic relevance of such agents has received increasing attention with the effectiveness of allopregnanolone and its synthetic derivative, zuranolone, in treating post-partum depression12. In addition to positive modulators like allopregnanolone, several neuroactive steroids have been shown to inhibit the ρ1 subtype, although the mechanistic basis of negative modulation remains controversial9. Compounds that negatively modulate ρ1 include sulfated neurosteroids and β-estradiol (E2)13.

Pregnenolone sulfate (PS) was one of the first identified neurosteroids, that is, steroids synthesized locally in the central or peripheral nervous system14. It is thought to exert excitatory effects, in part by suppressing neuro-inhibitory signaling via GABAA receptors15. The specific site(s) and mechanism of PS inhibition are unclear, though physiological, biochemical, and recent structural evidence support a role for pore-facing residues in classical synaptic GABAA receptors4,10,16–18. The estrogen steroid E2 is the major female sex hormone, involved in the development of the reproductive system and secondary sex characteristics, and in regulation of the menstrual cycle19. This hormone is mainly produced in ovaries, but also in other tissues including the brain, and is correlated with mood disorders20. It primarily binds and activates two nuclear receptors21,22, but also mediates rapid and non-genomic effects via membrane proteins such as the G-protein coupled estrogen receptor23. Estrogens have also been shown to mediate rapid actions on ligand-gated ion channels, for example, potentiating human α4β2 neuronal nicotinic24 and NMDA receptors25. In contrast, E2 effects on ρ1 are inhibitory, suggesting a notably distinct mechanism of modulation. Studies employing mutational analysis13, voltage-clamp fluorometry26, and functional modeling27 have shown that E2 and sulfated steroids bind to distinct sites and act through different mechanisms, though their respective details remain to be characterized.

Here, by combining five cryo-EM structures with electrophysiology and molecular dynamics simulations, we characterize the binding sites and negative modulation mechanisms of E2 and PS at ρ1. We find that E2 binds in a pocket at the ECD-TMD interface, apparently specific to the ρ subtypes, and disrupts allosteric conformational changes linking GABA binding to pore opening. In contrast, PS 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.

Results

E2 binds at the ECD-TMD interface of ρ1-EM

To explore distinctive steroid pharmacology in ρ1, we first characterized the functional effects of E2 (Fig. 1a) in our ρ1-EM construct. In agreement with previous reports13, 30 μM E2 reduced ρ1-EM currents in Xenopus oocytes by roughly half in the presence of 1 μM GABA (~ EC50) (Fig. 1b). Hypothesizing that E2 preferentially stabilizes a resting-like state of ρ1, we then determined a cryo-EM structure of ρ1-EM with E2. Like all structures in this and our previous work8, the receptor was reconstituted in saposin nanodiscs with polar brain lipids. We identified a single predominant conformation to an overall resolution of 2.5 Å with C5 symmetry (Supplementary Figs. 1–3 and Table 1). Although E2 has a similar backbone to neurosteroids like allopregnanolone, the intersubunit transmembrane site previously shown to bind allopregnanolone in α1β2ɣ210 only contained tubular densities in the ρ1-EM/E2 complex, similar to those observed in apo ρ1-EM8 and likely corresponding to phospholipid tails (Supplementary Fig. 3c). Instead, we observed a density corresponding in size and shape to E2 at the ECD-TMD interface of each pair of adjacent subunits (Fig. 1c, d).Fig. 1 E2 binds at the ECD-TMD interface of ρ1-EM.

a Chemical structure of β-estradiol (E2). b Sample traces from TEVC recordings of wild-type (black) and F283Y/F284Y (red) ρ1-EM constructs in response to GABA, with and without E2. c Cryo-EM structure of ρ1-EM with E2, viewed from the membrane plane. One subunit of the pentamer is colored dark red for definition. E2 (yellow) and resolved lipids (gray) are shown as thick and thin sticks, respectively. d TMD of ρ1-EM with E2, depicted as in panel c but viewed from the extracellular side. e Zoom view of a single E2 binding site, viewed from the membrane plane relative to the complementary (−) face of a single ρ1-EM subunit. The density assigned to E2 is shown in transparency. E2 and surrounding residues are shown as sticks and labeled. f Zoom view of a single E2 binding site, depicted as in panel e but rotated 90° to show the interface between two subunits from the channel pore. g Superimposed structures of apo (gray, PDB ID: 8OQ6) and E2-bound (red) ρ1-EM, showing no major change upon E2 binding. h Fractional current remaining after treatment with 30 μM E2 with ~ EC50 GABA (1 μM for wild-type, 4 μM for F283Y/F284Y ρ1-EM constructs). Error bars represent SEM from 7 individual oocytes, and stars represent p < 0.0001 (p = 6.016e-9) from a two-way t test.

Table 1 Cryo-EM data collection, refinement, and validation statistics

	E2 (EMD-19167) (PDB 8RH4)	GABA + E2 Primed (EMD-19171) (PDB 8RH7)	GABA + E2 Desensitized (EMD-19172) (PDB 8RH8)	PS (EMD-19173) (PDB 8RH9)	GABA + PS (EMD-19175) (PDB 8RHG)	
Data collection and processing	
Magnification	130,000	130,000	130,000	130,000	130,000	
Voltage (kV)	300	300	300	300	300	
Electron exposure (e–/Å2)	41.58	41.58	41.58	44.61	44.38	
Defocus range (μm)	− 0.8 to − 1.8	− 0.8 to − 1.8	− 0.8 to − 1.8	− 0.8 to − 1.8	− 0.8 to − 1.8	
Pixel size (Å)	0.6725	0.6725	0.6725	0.6725	0.6725	
Symmetry imposed	C5	C5	C5	C1	C1	
Final particles	134,816	100,833	140,148	93,154	88,968	
Map resolution (Å)

FSC threshold

	2.52

0.143

	2.78

0.143

	2.66

0.143

	3.21

0.143

	3.01

0.143

	
Map resolution range (Å)	2.3–2.7	2.6–3.0	2.5–2.9	3.0–3.8	2.9–3.7	
Refinement	
Initial model (PDB code)	8OQ6	8OQ6	8OP9	8OQ6	8OP9	
Model resolution (Å)

FSC threshold

	2.7

0.5

	2.9

0.5

	2.8

0.5

	3.3

0.5

	3.3

0.5

	
Map sharpening B factor (Å2)	− 97.2	− 112.2	− 109.6	− 101.0	− 97.9	
Model composition	
 Non-hydrogen atoms	14292	13866	14067	14247	13897	
 Protein residues	1660	1610	1645	1665	1650	
 Ligands	62	61	42	51	26	
B factors (Å2)	
 Protein	25.67	34.71	41.36	88.77	94.92	
 Ligand	46.63	44.83	72.40	85.43	75.75	
R.m.s. deviations	
 Bond lengths (Å)	0.010	0.005	0.007	0.003	0.003	
 Bond angles (°)	1.261	1.105	1.265	0.628	0.501	
Validation	
 MolProbity score	1.10	1.47	0.94	1.38	1.28	
 Clashscore	2.05	3.48	1.58	4.71	5.10	
 Poor rotamers (%)	0.32	0.34	0.33	0.06	0.07	
Ramachandran plot	
 Favored (%)	97.38	95.28	97.85	97.26	97.98	
 Allowed (%)	2.62	4.72	2.15	2.74	2.02	
 Disallowed (%)	0.00	0.00	0.00	0.00	0.00	

As verified by its protruding C16 methyl group (Fig. 1e), E2 fit unambiguously into this inter-domain density, with its C3 hydroxyl pointing up (toward the extracellular side), and C17 hydroxyl down (toward the intracellular side) (Fig. 1e, f). Overall, the E2 pocket is amphiphilic with local regions of positive charge (Supplementary Fig. 4a, b). The site is capped from the extracellular side by the β1-β2 loop and loop F, particularly the polar side chains of E113 and Q247 proximal to the C3 hydroxyl of E2 (Fig. 1e). From the transmembrane side, each E2 molecule is partially buried in a pocket enclosed by the upper M2-M3 region of the principal subunit, and by the upper M1 and M2 helices of the complementary subunit. On one face, the side chains of S334 and R337 are positioned to make hydrophobic and π-orbital interactions with E2 rings A and D, respectively (Fig. 1e, f). The opposite face approaches the hydrophobic surface of aromatic residues F283 and F284 at the amino-terminus of M1. Notably, substituting tyrosine for phenylalanine at these two positions largely ablated E2 inhibition while preserving GABA activation, indicating the precise geometry of this site critically determines E2 action (Fig. 1b, hand Supplementary Fig. 5b, c).

To our knowledge, small-molecule binding has not been previously shown for this pocket in any other GABAA-receptor structure. In the presence of E2, ρ1-EM is nearly identical to the previously reported apo structure (Fig. 2f), indicating that the steroid does not induce substantial conformational change. Even the local configuration of the binding pocket is preserved, with side chain rotamers of the surrounding residues maintained relative to the apo structure (Fig. 1g). In contrast, GABA binding rearranges residues including R337 in this region8, resulting in a pocket incompatible with E2 binding (Supplementary Fig. 4c). Interestingly, several residues proximal to E2, including F283, S334 and R347, were conserved among ρ1/2 but not α, β or ɣ subfamilies of GABAA receptors (Supplementary Fig. 6a). Moreover, no pocket capable of accommodating E2 was evident at any equivalent interface in the α1β2γ2 type (Supplementary Fig. 6b). In line with previous reports that classical synaptic GABAA receptors are insensitive to direct E2 modulation28, these comparisons suggest a ρ-specific binding site and inhibitory mechanism, which could inform future pharmacological design.Fig. 2 E2 suppresses activating transitions of the ECD upon GABA binding.

a Two different cryo-EM structures were obtained from a single sample of ρ1-EM with GABA and E2, viewed from the membrane plane. Presumed functional state (left, blue: desensitized; right, green: primed) and ensemble contribution as a percent of resolved particles are indicated above each structure. One subunit of each pentamer is colored darker for definition. E2 (yellow) and resolved lipids (gray) are shown as thick and thin sticks, respectively. b Zoom view of a single GABA binding site in the primed state, depicted as in panel a. Density assigned to GABA is shown in transparency. GABA and surrounding residues are shown as sticks and labeled. c Zoom view of a single E2 binding site in the primed state, depicted as in Fig. 1f from the channel pore. The density assigned to E2 is shown in transparency. E2 and surrounding residues are shown as sticks and labeled. d Superimposed structures of ρ1-EM in the apparent desensitized (blue) and primed (green) states, viewed from the membrane plane (left) and extracellular side (right). All but one subunit of each pentamer is rendered transparent for clarity. e GABA concentration-response curves for ρ1-EM in the absence (black) and presence of 30 μM E2 (pink). Error bars represent SEM from 5 individual oocytes. Solid lines represent fits Boltzmann curves with an EC50 of 1.3 μM (GABA alone, 95% confidence interval 1.16–1.41 μM) or 1.8 μM (GABA + E2, 95% confidence interval 1.56–2.21 μM). f Pore-radius profiles of ρ1-EM apo (gray) and E2 (red) structures. Both structures are assigned to a resting-like state. g Pore-radius profiles of ρ1-EM in the presence of GABA alone (black, PDB ID: 8OP9) and GABA + E2 in the primed (green) and desensitized (blue) states.

To explore the specificity of this E2 site, we aligned the rings of several related steroids into the ρ1-EM/E2 complex. The inter-domain site appears to accommodate 17α-estradiol, while the enantiomer ent-17β-estradiol clashes with the side chain of M2 residue S334 (Supplementary Fig. 4d). Consistent with these models, 17α-estradiol was previously shown to inhibit ρ1 similar to E2, while ent-17β-estradiol lacks modulatory effect13. The 5α neurosteroids allopregnanolone and tetrahydrodeoxycorticosterone (THDOC) are among the most structurally similar to E2, yet they have both been shown to potentiate rather than inhibit ρ14; the C3 hydroxyls of both these steroids are predicted to clash with R337 in our structures, suggesting they bind to an alternative site and/or state of the channel.

E2 suppresses activating transitions of the ECD upon GABA binding

To further investigate the structural basis for E2 modulation, we also solved the structures of ρ1-EM in the presence of both E2 and GABA. Under these conditions, we identified two well-resolved classes in the same dataset (Supplementary Fig. 1 and Table 1). One class, comprising 58% of assigned particles, was largely superimposable with the previously reported GABA-bound structure, activated by five molecules of GABA and assigned to a desensitized state8 (Fig. 2a, left). Notably, no E2 could be resolved in this structure. A second class, comprising 42% of assigned particles, also contained GABA in the orthosteric ligand sites, but with a global conformation markedly different from the desensitized state (Fig. 2a, right; Supplementary Fig. 3d, e). GABA binding in this second structure was associated with only a minor ECD twist of 1.2° compared with apo or E2-only conditions (Supplementary Fig. 4e), 7.3° less than in the desensitized state (Fig. 2d), and the pore is closed (Fig. 2g). Accordingly, we assigned this structure to a liganded pre-open state, possibly corresponding to one of the so-called primed states described in other pentameric ligand-gated ion channels29,30. We observed E2 in a site comparable to the E2-only complex (Fig. 2c), suggesting E2 disturbs allosteric GABA activation by wedging into the ECD-TMD interface between each pair of subunits.

The state dependence of E2 binding is reminiscent of the selective stabilization of picrotoxin (PTX) in the closed pore of ρ1-EM8, in line with previous reports that these inhibitors act through related mechanisms13. Indeed, apparent GABA affinity was reduced in the presence of E2, consistent with stabilization of a resting-like state (Fig. 2e). Moreover, fractional E2 inhibition decreased with increasing concentrations of GABA (Supplementary Fig. 5a, c), precluding a purely noncompetitive mechanism (e.g., pore block). A modest apparent reduction in maximal GABA efficacy (Fig. 2e) may represent an artifact of slow desensitization contributing to the steady-state inhibited current; indeed, this effect persisted at all E2 concentrations in both wild-type and F283Y/F284Y constructs (Supplementary Fig. 5a, c). On the other hand, given the limited effect of E2 with saturating GABA in electrophysiology experiments, it may seem surprising that E2 promotes such a substantial population in the primed state by cryo-EM, even under high-GABA conditions that would produce a single desensitized state without E28. However, concentrations of E2 applied in the cryo-EM samples are roughly an order of magnitude higher than concentrations used in electrophysiology, due to the improved solubility of E2 in the presence of lipids and detergents used in grid preparation. Thus, it is difficult to directly assess the functional effect of E2 at cryo-EM concentrations. Interestingly, an overlay of the E2- and PTX-bound structures in the presence of GABA shows that domain twist is even more limited by E2 than by PTX (Supplementary Fig. 4f), in line with previous voltage-clamp fluorometry data showing that the steroid suppresses upper-ECD rearrangements more than the toxin26.

PS occludes the GABA-bound state

To explore alternative mechanisms of ρ1 inhibition, we then characterized the 3β-sulfated neurosteroid PS (Fig. 3a), a negative modulator of this and several other GABAA-receptor subtypes15,31. Similar to previous reports13, 100 μM PS reduced ρ1-EM currents by roughly 35%, though there was a small recovery of the inhibited current in the continued presence of PS (Fig. 3b). This partial reversal of inhibition during PS exposure has also been observed in classical synaptic GABAA receptors, attributed to a redistribution amongst functional states upon PS binding31. Unlike E2, PS wash-out was associated with a transient recovery current ~ 30% larger than steady-state GABA activation prior to treatment (Fig. 3b). This behavior is expected for a pore blocker that preferentially binds to the activated (open or desensitized) state of the pore. Consistent with this model, current recovery upon PS washout was not observed at high concentrations of GABA, where channels are fully activated prior to PS treatment (Supplementary Fig. 7f).Fig. 3 PS occludes the GABA-bound state.

a Chemical structure of pregnenolone sulfate (PS). b Sample trace from TEVC recording of ρ1-EM in response to GABA and PS. The dotted line indicates maximum GABA response prior to PS application to highlight increased current upon washout of PS. c Cryo-EM structure of ρ1 with GABA and PS, viewed from the membrane plane. One subunit of the pentamer is colored dark blue for definition. PS (yellow) and resolved lipids (gray) are shown as thick and thin sticks, respectively. d TMD of ρ1-EM with GABA and PS, depicted as in panel c but viewed from the extracellular side. e Pore-radius profiles of ρ1-EM with GABA alone (black, PDB ID: 8OP9) and GABA + PS (blue). f Zoom view of the PS binding site. The density assigned to PS is shown in transparency. PS and inner pore-lining residues are shown as sticks and labeled. g Superimposition of structures with GABA and PS of ρ1-EM (blue) and α1β2γ2 (gray, PDB ID: 8SGO) GABAA receptors. h Background-subtracted and normalized current-voltage curves for voltage ramps applied to ρ1-EM in the presence of GABA alone (black) or in combination with 100 μM PS (blue). Shaded regions represent SEM from 4 individual oocytes.

To test this pore-block hypothesis, we determined a cryo-EM structure of ρ1-EM in the presence of PS and GABA. To avoid artifacts in a potential pore site, we processed these data without imposing symmetry, resolving a single predominant conformation to 3.2 Å (Supplementary Fig. 1 and Table 1). In addition to five GABA molecules at the orthosteric ECD sites, the resulting map contained a single density capable of accommodating PS, spanning residues P311 (− 2’) to L322 (9’) in the inner half of the channel pore (Fig. 3c, d and Supplementary Fig. 3d, e). At the inward-facing end of PS proximal to the − 2’ side chains, an additional spherical density was modeled as a chloride ion, also observed in our previous structures of ρ1-EM with GABA8 (Fig. 3f).

The PS density could accommodate modeling in two possible poses, with the sulfate group either facing up (toward the 9’ activation gate) or down (toward the −2’ desensitization gate) (Supplementary Fig. 7a). We tested the orientation of PS by running four replicate > 300 ns all-atom molecular dynamics (MD) simulations in each of the two poses (Table 2). Whereas the sulfate-up pose was relatively stable, the sulfate-down pose varied widely, displacing over a 14-Å range up or down the pore axis and > 5 Å median root-mean-square deviation (RMSD) (Supplementary Fig. 7b, c). Accordingly, we modeled PS with the sulfate up for all further analyses. We observed no other steroidal densities in the PS dataset, including in the E2 site or synaptic-subtype allopregnanolone site10.Table 2 System setup of MD simulations

	PS only sulfate up	PS only sulfate down	GABA + PS sulfate up	GABA + PS sulfate down	
Simulation box	118 Å x 118 Å x 188 Å	118 Å x 118 Å x 188 Å	120 Å x 120 Å x 183 Å	120 Å x 120 Å x 183 Å	
Number of atoms	269496	269496	268076	268076	
Number of waters	62349	62349	61958	61958	
Salt concentration	150 mM NaCl	150 mM NaCl	150 mM NaCl	150 mM NaCl	
Number of lipids	186 POPC	186 POPC	186 POPC	186 POPC	
124 POPE	124 POPE	124 POPE	124 POPE	
26 POPS	26 POPS	26 POPS	26 POPS	
70 Cholesterol	70 Cholesterol	70 Cholesterol	70 Cholesterol	
20 Sphingomyelins	20 Sphingomyelins	20 Sphingomyelins	20 Sphingomyelins	
18 PIP2	18 PIP2	18 PIP2	18 PIP2	

The PS pose in ρ1-EM overlapped that in a recently reported complex with the synaptic α1β2γ2 subtype, including the orientation of the sulfate group10 (Fig. 3g). A pore-block mechanism has similarly been proposed in this synaptic subtype, supported by mutations in the inner pore that suppress inhibition16 and disrupt PS stability in MD simulations10. The lower reported sensitivity of ρ1 versus synaptic subtypes to PS inhibition17 may be attributable to sequence differences in the channel pore, particularly at 2’, which is occupied by proline in ρ1, valine in α1, alanine in β2, and serine in ɣ2 (Supplementary Fig. 6a). Indeed, substituting the equivalent α1 residue at 2’ in ρ1 (P315V) has been shown to increase PS sensitivity17. As previously reported8, the ρ1 pore is also expanded relative to classical synaptic GABAA-receptor structures in the desensitized state (Fig. 3g), which could weaken contacts with a pore-bound ligand like PS.

To further validate this blocking mechanism, we compared ρ1-EM functional inhibition by PS to other inhibitors. PS inhibition was more efficacious at more positive potentials (Fig. 3h and Supplementary Fig. 7h), as expected for a negatively charged blocker. In contrast, ρ1-EM inhibition by the neutral blocker PTX was largely independent of voltage (Supplementary Fig. 7h–j). We estimated the fraction of the electric field traversed by the charged moiety of PS upon block by fitting the voltage-dependence of inhibition to the Woodhull model32. This model relies on two free parameters, the fraction of the electric field traversed (δ) and the affinity of the blocker to the pore at 0 mV (K D0mV). Our data for the PS block in ρ1-EM fit well to the Woodhull model with δ = 0.705 and KD0mV = 36.3 μM, indicating the charge traverses ~ 71% of the electric field of the pore. Some increase in the apparent block at positive potentials may actually reflect increased channel activity with increasing voltage that also favors PS block, resulting in a potential overestimate in δ33. Still, this error is likely to be small given that PTX, which shows an inverted activity-dependent block relative to PS, is nearly voltage-independent (Supplementary Fig. 7h–j). Comparison to PTX profiles suggests the activity dependence of the block contributes 0.05–0.09 of the apparent δ value.

The binding of PS in the inner ρ1 pore was reminiscent of our previously reported complex with PTX and GABA, including contacts at 2’. However, PTX selectively stabilizes an intermediate state in which the TMD is locked closed8. In contrast, the complex with PS and GABA was largely superimposable with our previous GABA-only structure, with an all-atom RMSD < 0.8 Å (Fig. 3e). Accordingly, the structure was presumed to be activated, occupying a desensitized state8. Modest changes were observed at either end of the PS site, subtly shifting the − 2’ and 9’ side chains towards the intracellular side and pore axis respectively (Supplementary Fig. 7d). Like PTX, PS decreased maximal GABA efficacy (Supplementary Fig. 7f, g); however, other electrophysiological properties distinguished the mechanisms of these two pore blockers. Whereas apparent GABA affinity decreases with PTX8, it increases with PS (Supplementary Fig. 7g), consistent with the steroid stabilizing an activated- rather than resting-like state. In contrast to E2, fractional inhibition by PS increased with increasing concentrations of GABA (Supplementary Fig. 5c), again consistent with preferential binding upon channel activation. Along with the recovery current observed after PS washout (Fig. 3b), these functional properties support a distinctive mechanism of PS inhibition by entering and binding to stabilize the activated pore.

PS has limited access to the resting-like pore

Finally, we determined a cryo-EM structure of ρ1-EM with PS alone, resolving a single conformation to 3 Å without imposing symmetry (Supplementary Fig. 1 and Table 1). As in the presence of GABA, we observed a PS-like density inside the pore, between the − 2’ and 9’ positions (Fig. 4a, b and Supplementary Fig. 3d, e). PS in this structure was more stable with its sulfate group down rather than up in MD simulations, likely due to the 9’ constriction precluding sulfate occupancy (Fig. 4c, d and Table 2). The steroid was also displaced 1.6 Å down toward the − 2’ gate, compared to its center of mass in the presence of GABA (Fig. 4e). No other densities in this structure were consistent with PS binding.Fig. 4 PS has limited access to the resting-like pore.

a Cryo-EM structure of ρ1-EM with PS, viewed from the membrane plane. One subunit of the pentamer is colored dark red for definition. PS (yellow) and resolved lipids (gray) are shown as thick and thin sticks, respectively. b Zoom views of the inner pore of ρ1-EM with PS, with experimental density assigned to PS and chloride shown in transparency. Two possible poses are shown for PS, either with the sulfate group oriented down towards the cytosol (left) or up towards the 9’ hydrophobic gate (right). PS (yellow), chloride (green), and surrounding residues are shown as sticks and labeled. c Translocation of PS along the pore z-axis in MD simulations launched from the two poses shown in B. Simulation frames are aligned on the Cα atoms of the M2 pore-lining helices, and translocation (ΔZ) calculated for the center of mass of PS non-hydrogen atoms along a linear axis passing through the channel pore. Violin plots represent probability densities from 4 independent simulation replicates of > 300 ns each, sampled every 0.4 ns (n = 3019 and 3271 for simulations with the sulfate oriented down and up, respectively). In each plot, the middle marker indicates the median, and error bars indicate minimum and maximum values. d Mobility of PS in MD simulations, calculated from RMSD of PS non-hydrogen atoms, with sample sizes and plot parameters as in panel (c). e Superimposition of the structures of ρ1-EM with PS alone (red) and with GABA + PS (blue). PS and pore-lining residues are shown as sticks and labeled. f Pore-radius profiles of apo (black, PDB ID: 8OQ6) and PS-bound (red) structures of ρ1-EM. g Potential of mean force-free energy for PS movement along the pore axis in PS-bound structures of ρ1-EM in the presence (red) and absence (blue) of GABA (9′ gate at 0 nm).

The ρ1-EM PS complex was assigned to a resting-like state, with no ligand in the orthosteric ECD sites and a radius < 2 Å at both the − 2’ and 9’ gates (Fig. 4f). It was largely similar to the previously reported apo structure of ρ1-EM8, with an all-atom RMSD of 0.3 Å. The most prominent difference was a modest expansion at the 2’ and 6’ positions, presumably to accommodate the proximal steroid rings of PS (Fig. 4e, f). PS occupancy in a closed pore was surprising, as the steroid radius is at least 5 Å, too large to transit the constrictions at either − 2’ or 9’. Indeed, the bulky steroid rings never fully exited either the − 2’ or 9’ gates in our MD simulations. Using enhanced sampling simulations, we calculated a free-energy barrier > 50 kJ/mol for PS to pass the − 2’ gate, and approaching 100 kJ/mol to pass 9’ (Fig. 4g and Supplementary Fig. 8). The barrier at 9’ was absent in the complex with PS and GABA (Fig. 4g), indicating the steroid can freely enter the pore from the extracellular side upon channel activation. Binding in the pore site was favorable relative to bulk, both in the PS structure and to an even greater extent in the structure with GABA + PS; this profile is consistent with occupancy in both cryo-EM structures and with preferential binding following activation of the 9’ gate.

The structure with PS suggests that transient rearrangements in the course of cryo-EM sample preparation, on the timescale of more than 30 min ligand incubation, allow the steroid to bind in the resting-like state, with apparent alterations in orientation and pose to accommodate the constricted pore. However, structures in the presence of GABA show that a closed pore is not preferentially stabilized by PS as it is by PTX or E2, and the complex with PS and GABA appears to be the more relevant model for ρ1 functional inhibition.

Discussion

Our structural, functional, and computational results reveal distinct sites of action and divergent inhibitory mechanisms for the neuroactive steroids E2 and PS on a ρ1 GABAA receptor. E2 binds at the ECD-TMD interface and appears to act as a wedge, blocking allosteric domain rearrangements that link ECD GABA binding to TMD pore opening (Fig. 5a). The absence of clear E2 density in the desensitized state of the GABA + E2 dataset suggests that full activation of the receptor precludes E2 binding. In contrast, the opening of the 9’ activation gate enables PS to bind inside the pore, blocking ion permeation (Fig. 5b). Preferential stabilization of the activated receptor is clear from the increased apparent GABA affinity and transient increase in current amplitude upon PS washout at low GABA concentrations. No binding was evident for either agent at the inner-leaflet site classically associated with allopregnanolone potentiation of α1β2γ2 GABAA receptors; indeed, this site diverges in sequence especially at a key M1 position (α1-Q242/ρ1-W300), likely accounting for the limited allopregnanolone sensitivity of ρ19. Although our results cannot entirely exclude transient occupancy of additional sites, they highlight the capacity of different steroids to modulate GABAA receptors via structurally distinct, largely exclusive mechanisms. Whereas the mechanism of PS inhibition we report here for ρ1 largely resembles that proposed for the classical synaptic α1β2γ2 GABAA receptor10, the site of E2 inhibition appears specific to ρ subtypes.Fig. 5 Negative modulation mechanisms of ρ-type GABAA receptors by steroids evidenced in this work.

Cartoons are derived from structures of ρ1-EM determined in this and previous work. Structures without GABA (red) in the absence or presence of E2 correspond to resting-like states. One structure with GABA and E2 (green) is assigned to a primed state. Otherwise, activation by GABA (blue) induces agonist-induced transitions in the ECD and 9’ hydrophobic gate, which are retained in the context of the PS block; corresponding experimental structures are parsimoniously assigned to desensitized states. a E2 wedges into the ECD-TMD interface to disrupt allosteric conformational transitions linking GABA binding to ECD rotation and pore opening. b PS binds inside the pore to block ion permeation, with an apparent preference for activated structures.

The buried domain-interface site observed here for E2 is relatively unexplored as a direct mediator of allosteric modulation and, to our knowledge has yet to be visualized in any known pentameric ligand-gated ion channel structure. It is notably distant from steroid sites in previous GABAA-receptor structures, potentially accessible from the extracellular medium rather than upon partitioning into the membrane core. Particularly in the absence of a definitive open state, the limited timescales of classical MD preclude definitive modeling of mechanistic effects of an allosteric ligand like E2; nonetheless, conformational changes in the desensitized state due to expansion of the upper pore are likely to also reflect the open state, and could similarly block E2 binding. The discovery of this evident modulatory site in ρ1 is particularly notable, given that this subfamily is thought to lack classical allosteric sites for benzodiazepines, barbiturates, and general anesthetics. Development of ρ-specific modulators has focused largely on the orthosteric GABA site, where TPMPA and related compounds bind; the E2 site could constitute a novel development target. E2 itself is known to promote excitability in the hippocampus by suppressing GABA signaling, although this effect has been primarily attributed to the alteration of classical synaptic GABAA-receptor expression via nuclear receptors34. Given its IC50 (6.5 μM13) is nearly one hundred times higher than circulating levels (≤ 150 nM in non-pregnant women35), E2 inhibition of ρ1 may play a limited physiological role. Furthermore, therapeutic applications of this site would require selecting against other E2 targets. Nonetheless, this steroid appears to constitute a promising lead compound for the design of ρ-type specific inhibitors, potentially useful in the treatment of visual, sleep, or cognitive disorders7.

Although long suspected, the pore block of GABAA receptors by sulfated neurosteroids has also been controversial, due in part to inconsistent evidence for voltage dependence in synaptic subtypes. Here, we demonstrate that the PS block of ρ1 is indeed mildly voltage-dependent, as expected from the negatively charged sulfate group interacting with the electric field across the pore. The expanded pore of ρ1 versus classical synaptic GABAA receptors in the presence of GABA8 could contribute to the relative robustness of PS block to pore mutations13. Interestingly, inhibition by the related compound pregnanolone sulfate was previously shown to be voltage-dependent at α1β2γ2 GABAA receptors but voltage-independent at ρ113,33. It is possible that PS and pregnanolone sulfate act at different sites; indeed, at least three distinct mechanisms of ρ1 inhibition have been proposed for different steroids13. Alternatively, subtle differences in the position of the sulfate group or local pore structure may position the charged moiety outside the electric field gradient. Although the structural details conferring voltage sensitivity (or lack thereof) on steroid inhibition remain unclear, our structures of PS-bound ρ1 combined with its electrophysiological profile, as well as molecular simulation data coherently support a pore-blocking mechanism for this agent.

Another interesting feature of PS inhibition is its enhancement by GABA activation. Unlike most steroids that primarily modulate ρ1 receptor function at low GABA concentrations13, PS inhibits maximally at saturating GABA. This profile indicates that PS binds preferentially in the context of GABA activation, giving rise in cryo-EM to a partially open or desensitized state. Although PS binding in the resting-like state is evidently possible in the context of prolonged incubation for cryo-EM, it would be disfavored by the permeation barrier at the 9’ gate. This barrier was estimated at 100 kJ/mol by our PMF calculations: although these measurements do not directly report conductance, the central hydrophobic gate is expected to preclude inward transit of PS as well as chloride ions in a physiological setting. Furthermore, the recovery current apparent upon PS washout supports preferential binding to an open rather than desensitized state.

On the other hand, minor rearrangements are apparent in our structures with GABA versus GABA + PS, including a modest constriction of the inner desensitization gate relative to the structure with GABA alone. In contrast, our past and current work with inhibitors favoring the resting-like state of ρ1 (TPMPA, PTX, E2) demonstrates these compounds require little or no rearrangement around the respective binding sites, even on a local scale. These findings may indicate that cryo-EM sample conditions favor a desensitized-like structure that does not perfectly represent the predominant physiological ensemble36, which should include one or more open structures. Alternatively, these results may support a recent hypothesis based on detailed kinetic modeling of the α1β3γ2 GABAA receptor, where PS binding stabilizes a nonconducting state distinct from both open and desensitized31. Such a model would recapitulate several functional features we observe here, including increased apparent GABA affinity in the presence of PS and a transient increase in current amplitude upon PS washout.

Taken together, our findings expand on a growing body of literature demonstrating that despite similar structural backbones, neuroactive steroids can have diverse binding sites and mechanisms of action on GABAA receptors. The importance of neuroactive steroids as building blocks for new therapies is clear, given the recent success of the endogenous modulator allopregnanolone and synthetic derivative zuranolone in the treatment of postpartum depression12. The structures we report here can aid future structure-based drug design to better target ρ-type receptors, which are insensitive to nearly all classical GABAA receptor-targeting therapies.

Methods

Protein expression and purification

The expression-optimized human ρ1 construct (ρ1-EM) was expressed and purified according to previous methods8. Briefly, baculovirus encoding ρ1-EM was amplified by infecting a 300-mL suspension Sf9 cell culture (Novagen). Expi293F cells (Gibco) were infected by baculovirus at a density of 2 × 106 cells/mL. After 6 h incubation at 37 °C, 5 mM sodium butyrate was added and the cells were further cultured at 30 °C for 48 h. Cells were harvested, washed with phosphate-buffered saline, and then flash-frozen until further usage.

For sample preparation of the PS datasets, cell pellets from 2 L culture were resuspended in resuspension buffer (40 mM HEPES pH 7.5, 300 mM NaCl, with cOmplete protease inhibitor tablets (Roche)) and sonicated to break cell membranes. The membranes were pelleted by ultracentrifugation, then resuspended and solubilized in resuspension buffer with 2% lauryl maltose neopentyl glycol (LMNG) and 0.2% cholesteryl hemisuccinate (CHS) for 3 h in a cold room (4–10 °C). The solubilization mixture was ultracentrifuged and the supernatant was applied to 4-mL Strep-Tactin Superflow resin (IBA) and incubated for 90 min. The resin was washed with wash buffer (20 mM HEPES pH 7.5, 300 mM NaCl, 0.005% LMNG, 0.0005% CHS), then the protein was eluted with elution buffer (wash buffer with 10 mM d-Desthiobiotin (Sigma)). The product was further purified by size exclusion chromatography on a Superose 6 column (Cytiva) in flow buffer (20 mM HEPES pH 7.5, 100 mM NaCl, 0.005% LMNG, 0.0005% CHS). Peak fractions were pooled for nanodisc reconstitution. The sample for E2 datasets was purified the same way, except CHS was not included in the purification.

Nanodisc reconstitution

The plasmid for SapA expression was a gift from Salipro Biotech AB. Purification of SapA followed previously published protocols37. For the reconstitution of saposin nanodiscs for the PS datasets, ρ1-EM, SapA, and porcine polar brain lipid (Avanti) were mixed at a molar ratio of 1:15:150, then incubated on ice for 1 h. Bio-Beads SM-2 resin (Bio-Rad) was added into the mixture, then gently rotated overnight at 4 °C. On the next day, the supernatant was collected and further purified by gel-filtration chromatography on a Superose 6 column (Cytiva) with a buffer containing 20 mM HEPES pH 7.5 and 100 mM NaCl. Peak fractions were pooled and concentrated to ~ 5 mg/mL. For the E2 datasets, E2 was mixed with polar brain lipids at a 1:10 molar ratio to form the lipid mixture. The following process was the same as for the PS sample, except the E2 lipid mixture was used.

Cryo-EM grid preparation and data collection

Nanodisc samples were mixed with additive stock solutions in a 9:1 volume ratio, and incubated ≥ 30 min on ice. Stock solutions were prepared for data collection with E2 (2 mM E2, 20 mM fluorinated fos-choline 8 (FFC-8), 0.5% DMSO), GABA + E2 (6 mM GABA, 2 mM E2, 20 mM FFC-8, 0.5% DMSO), PS (10 mM PS, 20 mM FFC-8, 0.5% DMSO) and GABA + PS (6 mM GABA, 10 mM PS, 20 mM FFC-8, 0.5% DMSO).

Right before application to the grid, each mixture was centrifuged to remove potential precipitation. 3 μL of the supernatant was then applied to a glow-discharged grid (R1.2/1.3 300 mesh Au grid, Quantifoil), blotted for 2 s with force 0, and plunged into liquid ethane using a Vitrobot Mark IV (Thermo Fisher Scientific).

Cryo-EM data were collected on a 300 kV Titan Krios (Thermo Fisher Scientific) electron microscope with a K3 Summit detector (Gatan) with magnification 105 k corresponding to 0.8464 Å/px using the software EPU 3.5.0 (Thermo Fisher Scientific). The total dose was ~ 42e-/Å2, and the defocus range was − 0.8 to − 1.8 μm.

Cryo-EM data processing

Dose-fractionated images in super-resolution mode were internally gain-normalized and binned by 2 in EPU during data collection. Cryo-EM data processing was first done in RELION 3.1.438, including motion correction, contrast transfer function (CTF) estimation with CTFFIND 4.139, automatic particle picking with Topaz 0.2.540, particle extraction, 2D classification, 3D classification, 3D refinement, CTF refinement, and polishing. Briefly, two rounds of 2D classification were done to remove junk particles, and 3D classification (classes = 4) was used to assess structural heterogeneity. Particles from classes with protein features were centered and re-extracted and were used for 3D refinement with C5 (E2 datasets) or C1 (PS datasets) symmetry. Multiple rounds of CtfRefine and one or two rounds of Bayesian polishing were executed to improve resolution. Shiny particles were imported into CryoSPARC v4.2.1 for further processing41, including 3D classification in PCA mode and non-uniform refinement42.

Model building and refinement

Model building was started with rigid body fitting of the previously published apo (PDB ID 8OQ6) or GABA-bound (PDB ID 8OP9) structure into the density. The models were manually checked and adjusted in Coot 0.9.543, and ligands, water, and lipids were manually added. The resulting models were further optimized using real-space refinement in PHENIX 1.18.244 and validated by MolProbity45. Pore radius profiles were calculated using CHAP 0.9.146. Structure figures were prepared using UCSF ChimeraX 1.347.

Expression in oocytes and electrophysiology

mRNA encoding the ρ1-EM GABAA receptor was produced by in-vitro transcription using the mMessage mMachine T7 Ultra transcription kit (Ambion) according to the manufacturer protocol. Xenopus laevis oocytes (Ecocyte Bioscience) were injected with 30–50 ng mRNA and incubated 4–8 days at 13 °C in post-injection solution (10 mM HEPES pH 8.5, 88 mM NaCl, 2.4 mM NaHCO3, 1 mM KCl, 0.91 mM CaCl2, 0.82 mM MgSO4, 0.33 mM Ca(NO3)2, 2 mM sodium pyruvate, 0.5 mM theophylline, 0.1 mM gentamicin, 17 mM streptomycin, 10,000 u/L penicillin) prior to two-electrode voltage clamp (TEVC) measurements. Mutagenesis was performed by methods analogous to QuikChange cloning, and the sequence was verified across the entire coding length of the gene.

For TEVC recordings, glass electrodes were pulled and filled with 3 M KCl to give a resistance of 0.5–1.5 MΩ and used to clamp the membrane potential of injected oocytes at − 60 mV with an OC-725C voltage clamp (Warner Instruments). Oocytes were maintained under continuous perfusion with Ringer’s solution (123 mM NaCl, 10 mM HEPES, 2 mM KCl, 2 mM MgSO4, 2 mM CaCl2, pH 7.5) at a flow rate of around 1.5 mL/min. Buffer exchange was accomplished by manually switching the inlet of the perfusion system to the appropriate buffer. Currents were digitized at a sampling rate of 2 kHz and lowpass filtered at 10 Hz with an Axon CNS 1440 A Digidata system controlled by pCLAMP 10 (Molecular Devices).

GABA dose-response curves in the presence and absence of steroids were measured using a 90 s co-application of 30 μM E2 or 100 μM PS during a 3.5–5.5 min pulse of GABA as shown in Fig. 1b. Each oocyte was recorded across the full GABA concentration range of 0.5–16 μM for ρ1-EM or 2–32 μM for the F283Y/F284Y mutant. Currents at the end of the initial GABA only (prior to steroid application) and GABA + steroid pulses were normalized to the maximum current measured from each oocyte. Concentration-response curves were fit using least squares regression considering each Y-value from separate oocytes as individual points using Prism with uncertainty in fitting parameters output as a 95% confidence interval.

Voltage-dependent block experiments were performed similarly to GABA and PS block experiments, with a few modifications. The holding potential for the voltage-dependent block was − 80 mV, and automated voltage ramps from − 80 mV to 30 mV were performed over 4 s in Ringer’s solution only, upon saturation of the 1 μM GABA response, and upon saturation of the 1 μM GABA and 100 μM PS (or 0.5 μM PTX) response. Current elicited in the absence of GABA and PS was subtracted from other responses to remove contributions of leak, capacitive, and endogenous currents for each oocyte. The Woodhull model32 was used to evaluate the fraction of the electric field traversed by the charge of PS using the following equation and fit in Prism:1 p=KD0mV/(KD0mV+PS*e(−n*F*δ*V/(R*T)))

where p is the fraction of receptors not blocked, KD0mV is the dissociation constant at 0 mV, [PS] is the concentration of PS used in the experiments (100 μM), n is the charge, F is Faraday’s constant, R is the gas constant, T is temperature, V is the membrane potential, and δ is the fraction of the field traversed. Voltage dependence of PTX was also fit to a Woodhull model to assess the contribution of the state-dependent block to the measured δ value. The charge was assigned as + 1 due to preference of block in the resting state which is favored at negative potentials according to the rectification behavior in the control conditions. The δ value for PTX was best fit at 0.071 with a 95% confidence interval of 0.053–0.089.

Molecular dynamics simulations

Detailed dataset-specific information can be found in Table 2. Atomic coordinates for the ρ1-EM determined by cryo-EM with different neurosteroid poses were used as starting models for MD simulations. Each subunit was split into two chains for simulation, due to the disconnection between the M3-M4 loop in the structure. The simulation systems were set up in CHARMM-GUI48. The protein was embedded into a lipid mixture mimicking brain lipid compositions49, with the outer leaflet containing 152 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 14 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), 38 cholesterol and 15 sphingomyelin molecules, and the inner leaflet containing 34 POPC, 110 POPE, 26 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine (POPS), 32 cholesterol, 5 sphingomyelin and 18 phosphatidylinositol 4,5-bisphosphate (PIP2) molecules. The protein-lipid complex was subsequently solvated with TIP3P water and 150 mM NaCl. The CHARMM36m forcefield50 was used to describe the proteins. Parameters for the neurosteroid were reproduced from previous work10, in which they were generated by CGenFF51 in CHARMM-GUI48.

Simulations were performed using GROMACS 2021.552 with temperature coupled to 300 K using the velocity-rescaling thermostat53 and pressure of 1 atm using a Parrinello-Rahman barostat54. The LINCS algorithm was used to constrain hydrogen-bond lengths55, and the particle mesh Ewald method56 was used to calculate long-range electrostatic interactions. The systems were energy minimized and then equilibrated for 20 ns, with the position restraints on the protein, and neurosteroids were gradually released. Four replicates each of 300–400 ns were simulated as final unrestrained production runs.

Before analysis, MD simulation trajectories were aligned on the Cα atoms of M2 helices by MDAnalysis57. Root mean square deviations (RMSD) and pore axis movement of ligands were calculated in VMD 1.9.358 and visualized with Matplotlib59.

The potential of mean force (PMF) for PS permeating the pore was calculated using the accelerated weight histogram (AWH) method60, following the protocol outlined in the GROMACS manual (https://tutorials.gromacs.org/docs/awh-tutorial.html). Relevant parameters added to the.mdp file are listed below; for access to the complete inputs, see Data Availability:

Awh = yes; AWH on.

Awh-nbias = 1; One bias, could have multiple.

Awh1-ndim = 1; Dimensionality of the RC

Awh1-dim1-coord-index = 1; Map RC dimension to pull the coordinate index

Awh1-dim1-start  = − 9.1275; Sampling interval min value (nm)

Awh1-dim1-end = 9.1275; Sampling interval max value (nm)

Awh1-dim1-force-constant = 128000; Force const of the harmonic pot. (kJ/(mol*nm^2))

Awh1-dim1-diffusion = 5e-5; Estimate of the diffusion (nm^2/ps)

The simulation was run for 650 ns with 4 walkers sharing biased data and contributing to the same target distribution. The Cα atoms of the protein were restrained to preserve the channel in a state corresponding to experimental conditions. To prevent the neurosteroid from flipping during simulation, a flat-bottom potential of radius 8 Å was added for its upper- and lower-most atoms. Convergence was checked by monitoring the exit time of AWH from initial to final stages (580 and 560 ns for closed and desensitized simulations, respectively), and the stability of PMFs at timepoints in the final stage (Supplementary Fig. 8).

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Reporting Summary

Source data

Source data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51904-7.

Acknowledgements

We thank professors Ryan E Hibbs, Alex Evers and the members of Molecular Biophysics Stockholm for feedback on the project and manuscript, and staff at the Cryo-EM Swedish National Facility for data collection support. Cryo-EM data were collected at the Facility funded by the Knut and Alice Wallenberg, Family Erling Persson and Kempe Foundations, SciLifeLab, and Stockholm University. MD simulations were performed using the computing facilities of the Swedish National Infrastructure for Computing (SNIC 2022/3-40, E.L.), and supported by BioExcel (EuroHPC grant no. 101093290, E.L.). J.C. was supported by an EMBO Postdoctoral Fellowship and C.F. by grant FV-5.1.2-0523-19 from Stockholm University; R.J.H. and E.L. acknowledge grants from the Swedish Research Council (2019-02433, 2021-05806) and Swedish e-Science Research Center.

Author contributions

C.F. and J.C. performed the biochemistry, cryo-EM sample preparation, and data processing. C.F. performed model building, refinement, structural analysis, and MD simulations. J.C. performed electrophysiology. R.J.H. and E.L. supervised the project. All authors contributed to the manuscript writing and revision.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Funding

Open access funding provided by Stockholm University.

Data availability

The cryo-EM maps and the corresponding atomic coordinates have been deposited in the Electron Microscopy Data Bank (EMDB) and the Protein Data Bank (PDB) for E2 (EMD-19167 [https://www.ebi.ac.uk/pdbe/entry/emdb/EMD-19167], PDB-8RH4 [10.2210/pdb8RH4/pdb]), GABA + E2 primed state (EMD-19171 [https://www.ebi.ac.uk/pdbe/entry/emdb/EMD-19171], PDB-8RH7 [10.2210/pdb8RH7/pdb]), GABA + E2 desensitized state (EMD-19172 [https://www.ebi.ac.uk/pdbe/entry/emdb/EMD-19172], PDB-8RH8 [10.2210/pdb8RH8/pdb]), PS (EMD-19173 [https://www.ebi.ac.uk/pdbe/entry/emdb/EMD-19173], PDB-8RH9 [10.2210/pdb8RH9/pdb]), GABA + PS (EMD-19175 [https://www.ebi.ac.uk/pdbe/entry/emdb/EMD-19175], PDB-8RHG [10.2210/pdb8RHG/pdb]). MD simulation trajectories and parameter files are available in Zenodo (10406748 [10.5281/zenodo.10406748]). Source data are provided in this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Chen Fan, John Cowgill.
==== Refs
References

1. Howard RJ Elephants in the dark: Insights and incongruities in pentameric ligand-gated ion channel models J. Mol. Biol. 2021 433 167128 10.1016/j.jmb.2021.167128 34224751
Howard, R. J. Elephants in the dark: Insights and incongruities in pentameric ligand-gated ion channel models. J. Mol. Biol. 433, 167128 (2021).34224751 10.1016/j.jmb.2021.167128
2. Drew CA Johnston GAR Weatherby RP Bicuculline-insensitive GABA receptors: Studies on the binding of (−)-baclofen to rat cerebellar membranes Neurosci. Lett. 1984 52 317 321 10.1016/0304-3940(84)90181-2 6097844
Drew, C. A., Johnston, G. A. R. & Weatherby, R. P. Bicuculline-insensitive GABA receptors: Studies on the binding of (−)-baclofen to rat cerebellar membranes. Neurosci. Lett. 52, 317–321 (1984).6097844 10.1016/0304-3940(84)90181-2
3. Enz R Cutting GR GABAC receptor ρ subunits are heterogeneously expressed in the human CNS and form homo‐ and heterooligomers with distinct physical properties Eur. J. Neurosci. 1999 11 41 50 10.1046/j.1460-9568.1999.00423.x 9987010
Enz, R. & Cutting, G. R. GABAC receptor ρ subunits are heterogeneously expressed in the human CNS and form homo‐ and heterooligomers with distinct physical properties. Eur. J. Neurosci. 11, 41–50 (1999).9987010 10.1046/j.1460-9568.1999.00423.x
4. Johnston G Medicinal chemistry and molecular pharmacology of GABA-C receptors Curr. Top. Med. Chem. 2002 2 903 913 10.2174/1568026023393453 12171579
Johnston, G. Medicinal chemistry and molecular pharmacology of GABA-C receptors. Curr. Top. Med. Chem. 2, 903–913 (2002).12171579 10.2174/1568026023393453
5. Varman DR Soria-Ortíz MB Martínez-Torres A Reyes-Haro D GABAρ3 expression in lobule X of the cerebellum is reduced in the valproate model of autism Neurosci. Lett. 2018 687 158 163 10.1016/j.neulet.2018.09.042 30261230
Varman, D. R., Soria-Ortíz, M. B., Martínez-Torres, A. & Reyes-Haro, D. GABAρ3 expression in lobule X of the cerebellum is reduced in the valproate model of autism. Neurosci. Lett. 687, 158–163 (2018).30261230 10.1016/j.neulet.2018.09.042
6. Van Nieuwenhuijzen PS Targeting GABAC receptors improves post-stroke motor recovery Brain Sci. 2021 11 315 10.3390/brainsci11030315 33801560
Van Nieuwenhuijzen, P. S. et al. Targeting GABAC receptors improves post-stroke motor recovery. Brain Sci. 11, 315 (2021).33801560 10.3390/brainsci11030315
7. Johnston G Chebib M Hanrahan J Mewett K GABAC Receptors as drug targets Curr. Drug Target CNS Neurol. Disord. 2003 2 260 268 10.2174/1568007033482805
Johnston, G., Chebib, M., Hanrahan, J. & Mewett, K. GABAC Receptors as drug targets. Curr. Drug Target CNS Neurol. Disord. 2, 260–268 (2003).10.2174/1568007033482805
8. Cowgill J Structure and dynamics of differential ligand binding in the human ρ-type GABAA receptor Neuron 2023 111 3450 3464 10.1016/j.neuron.2023.08.006 37659407
Cowgill, J. et al. Structure and dynamics of differential ligand binding in the human ρ-type GABAA receptor. Neuron 111, 3450–3464 (2023).37659407 10.1016/j.neuron.2023.08.006
9. Morris KD Moorefield CN Amin J Differential modulation of the gamma-aminobutyric acid type C receptor by neuroactive steroids Mol. Pharmacol. 1999 56 752 759 10496958
Morris, K. D., Moorefield, C. N. & Amin, J. Differential modulation of the gamma-aminobutyric acid type C receptor by neuroactive steroids. Mol. Pharmacol. 56, 752–759 (1999).10496958
10. Legesse DH Structural insights into opposing actions of neurosteroids on GABAA receptors Nat. Commun. 2023 14 5091 10.1038/s41467-023-40800-1 37607940
Legesse, D. H. et al. Structural insights into opposing actions of neurosteroids on GABAA receptors. Nat. Commun. 14, 5091 (2023).37607940 10.1038/s41467-023-40800-1
11. Sun C Zhu H Clark S Gouaux E Cryo-EM structures reveal native GABAA receptor assemblies and pharmacology Nature 2023 622 195 201 10.1038/s41586-023-06556-w 37730991
Sun, C., Zhu, H., Clark, S. & Gouaux, E. Cryo-EM structures reveal native GABAA receptor assemblies and pharmacology. Nature 622, 195–201 (2023).37730991 10.1038/s41586-023-06556-w
12. Marecki R Kałuska J Kolanek A Hakało D Waszkiewicz N Zuranolone - synthetic neurosteroid in treatment of mental disorders: narrative review Front. Psychiatry 2023 14 1298359 10.3389/fpsyt.2023.1298359 38116383
Marecki, R., Kałuska, J., Kolanek, A., Hakało, D. & Waszkiewicz, N. Zuranolone - synthetic neurosteroid in treatment of mental disorders: narrative review. Front. Psychiatry 14, 1298359 (2023).38116383 10.3389/fpsyt.2023.1298359
13. Li W Jin X Covey DF Steinbach JH Neuroactive steroids and human recombinant rho1 GABAC receptors J. Pharmacol. Exp. Ther. 2007 323 236 247 10.1124/jpet.107.127365 17636008
Li, W., Jin, X., Covey, D. F. & Steinbach, J. H. Neuroactive steroids and human recombinant rho1 GABAC receptors. J. Pharmacol. Exp. Ther. 323, 236–247 (2007).17636008 10.1124/jpet.107.127365
14. Schumacher M Pregnenolone sulfate in the brain: a controversial neurosteroid Neurochem. Int. 2008 52 522 540 10.1016/j.neuint.2007.08.022 18068870
Schumacher, M. et al. Pregnenolone sulfate in the brain: a controversial neurosteroid. Neurochem. Int. 52, 522–540 (2008).18068870 10.1016/j.neuint.2007.08.022
15. Majewska MD Mienville JM Vicini S Neurosteroid pregnenolone sulfate antagonizes electrophysiological responses to GABA in neurons Neurosci. Lett. 1988 90 279 284 10.1016/0304-3940(88)90202-9 3138576
Majewska, M. D., Mienville, J. M. & Vicini, S. Neurosteroid pregnenolone sulfate antagonizes electrophysiological responses to GABA in neurons. Neurosci. Lett. 90, 279–284 (1988).3138576 10.1016/0304-3940(88)90202-9
16. Wang M-D Rahman M Zhu D Johansson I-M Bäckström T 3Beta-hydroxysteroids and pregnenolone sulfate inhibit recombinant rat GABA(A) receptor through different channel property Eur. J. Pharmacol. 2007 557 124 131 10.1016/j.ejphar.2006.11.071 17239367
Wang, M.-D., Rahman, M., Zhu, D., Johansson, I.-M. & Bäckström, T. 3Beta-hydroxysteroids and pregnenolone sulfate inhibit recombinant rat GABA(A) receptor through different channel property. Eur. J. Pharmacol. 557, 124–131 (2007).17239367 10.1016/j.ejphar.2006.11.071
17. Seljeset S Bright DP Thomas P Smart TG Probing GABAA receptors with inhibitory neurosteroids Neuropharmacology 2018 136 23 36 10.1016/j.neuropharm.2018.02.008 29447845
Seljeset, S., Bright, D. P., Thomas, P. & Smart, T. G. Probing GABAA receptors with inhibitory neurosteroids. Neuropharmacology 136, 23–36 (2018).29447845 10.1016/j.neuropharm.2018.02.008
18. Akk G Bracamontes J Steinbach JH Pregnenolone sulfate block of GABA(A) receptors: mechanism and involvement of a residue in the M2 region of the alpha subunit J. Physiol. 2001 532 673 684 10.1111/j.1469-7793.2001.0673e.x 11313438
Akk, G., Bracamontes, J. & Steinbach, J. H. Pregnenolone sulfate block of GABA(A) receptors: mechanism and involvement of a residue in the M2 region of the alpha subunit. J. Physiol. 532, 673–684 (2001).11313438 10.1111/j.1469-7793.2001.0673e.x
19. Hall, J. E., Hall, M. E. & Guyton, A. C. Guyton and Hall Textbook of Medical Physiology. (Elsevier, Philadelphia, PA, 2021).
20. Douma SL Husband C O’Donnell ME Barwin BN Woodend AK Estrogen-related mood disorders: reproductive life cycle factors Adv. Nurs. Sci. 2005 28 364 375 10.1097/00012272-200510000-00008
Douma, S. L., Husband, C., O’Donnell, M. E., Barwin, B. N. & Woodend, A. K. Estrogen-related mood disorders: reproductive life cycle factors. Adv. Nurs. Sci. 28, 364–375 (2005).10.1097/00012272-200510000-00008
21. Greene GL Sequence and expression of human estrogen receptor complementary DNA Science 1986 231 1150 1154 10.1126/science.3753802 3753802
Greene, G. L. et al. Sequence and expression of human estrogen receptor complementary DNA. Science 231, 1150–1154 (1986).3753802 10.1126/science.3753802
22. Kuiper GG Enmark E Pelto-Huikko M Nilsson S Gustafsson JA Cloning of a novel receptor expressed in rat prostate and ovary Proc. Natl. Acad. Sci. USA 1996 93 5925 5930 10.1073/pnas.93.12.5925 8650195
Kuiper, G. G., Enmark, E., Pelto-Huikko, M., Nilsson, S. & Gustafsson, J. A. Cloning of a novel receptor expressed in rat prostate and ovary. Proc. Natl. Acad. Sci. USA 93, 5925–5930 (1996).8650195 10.1073/pnas.93.12.5925
23. Prossnitz ER Estrogen signaling through the transmembrane G protein-coupled receptor GPR30 Annu. Rev. Physiol. 2008 70 165 190 10.1146/annurev.physiol.70.113006.100518 18271749
Prossnitz, E. R. et al. Estrogen signaling through the transmembrane G protein-coupled receptor GPR30. Annu. Rev. Physiol. 70, 165–190 (2008).18271749 10.1146/annurev.physiol.70.113006.100518
24. Paradiso K Zhang J Steinbach JH The C terminus of the human nicotinic alpha4beta2 receptor forms a binding site required for potentiation by an estrogenic steroid J. Neurosci. 2001 21 6561 6568 10.1523/JNEUROSCI.21-17-06561.2001 11517245
Paradiso, K., Zhang, J. & Steinbach, J. H. The C terminus of the human nicotinic alpha4beta2 receptor forms a binding site required for potentiation by an estrogenic steroid. J. Neurosci. 21, 6561–6568 (2001).11517245 10.1523/JNEUROSCI.21-17-06561.2001
25. Foy MR 17beta-estradiol enhances NMDA receptor-mediated EPSPs and long-term potentiation J. Neurophysiol. 1999 81 925 929 10.1152/jn.1999.81.2.925 10036289
Foy, M. R. et al. 17beta-estradiol enhances NMDA receptor-mediated EPSPs and long-term potentiation. J. Neurophysiol. 81, 925–929 (1999).10036289 10.1152/jn.1999.81.2.925
26. Eaton MM Lim YB Covey DF Akk G Modulation of the human ρ1 GABAA receptor by inhibitory steroids Psychopharmacology 2014 231 3467 3478 10.1007/s00213-013-3379-z 24317445
Eaton, M. M., Lim, Y. B., Covey, D. F. & Akk, G. Modulation of the human ρ1 GABAA receptor by inhibitory steroids. Psychopharmacology 231, 3467–3478 (2014).24317445 10.1007/s00213-013-3379-z
27. Germann AL Analysis of modulation of the ρ1 GABAA receptor by combinations of inhibitory and potentiating neurosteroids reveals shared and distinct binding sites Mol. Pharmacol. 2020 98 280 291 10.1124/mol.120.119842 32675382
Germann, A. L. et al. Analysis of modulation of the ρ1 GABAA receptor by combinations of inhibitory and potentiating neurosteroids reveals shared and distinct binding sites. Mol. Pharmacol. 98, 280–291 (2020).32675382 10.1124/mol.120.119842
28. Wójtowicz T Mozrzymas JW Estradiol and GABAergic transmission in the hippocampus Vitam. Horm. 2010 82 279 300 10.1016/S0083-6729(10)82015-1 20472144
Wójtowicz, T. & Mozrzymas, J. W. Estradiol and GABAergic transmission in the hippocampus. Vitam. Horm. 82, 279–300 (2010).20472144 10.1016/S0083-6729(10)82015-1
29. Sivilotti LG What single-channel analysis tells us of the activation mechanism of ligand-gated channels: the case of the glycine receptor J. Physiol. 2010 588 45 58 10.1113/jphysiol.2009.178525 19770192
Sivilotti, L. G. What single-channel analysis tells us of the activation mechanism of ligand-gated channels: the case of the glycine receptor. J. Physiol. 588, 45–58 (2010).19770192 10.1113/jphysiol.2009.178525
30. Mukhtasimova N Lee WY Wang H-L Sine SM Detection and trapping of intermediate states priming nicotinic receptor channel opening Nature 2009 459 451 454 10.1038/nature07923 19339970
Mukhtasimova, N., Lee, W. Y., Wang, H.-L. & Sine, S. M. Detection and trapping of intermediate states priming nicotinic receptor channel opening. Nature 459, 451–454 (2009).19339970 10.1038/nature07923
31. Pierce SR Germann AL Steinbach JH Akk G The sulfated steroids pregnenolone sulfate and dehydroepiandrosterone sulfate inhibit the α1β3γ2L GABAA receptor by stabilizing a novel nonconducting state Mol. Pharmacol. 2022 101 68 77 10.1124/molpharm.121.000385 34853153
Pierce, S. R., Germann, A. L., Steinbach, J. H. & Akk, G. The sulfated steroids pregnenolone sulfate and dehydroepiandrosterone sulfate inhibit the α1β3γ2L GABAA receptor by stabilizing a novel nonconducting state. Mol. Pharmacol. 101, 68–77 (2022).34853153 10.1124/molpharm.121.000385
32. Woodhull AM Ionic blockage of sodium channels in nerve J. Gen. Physiol. 1973 61 687 708 10.1085/jgp.61.6.687 4541078
Woodhull, A. M. Ionic blockage of sodium channels in nerve. J. Gen. Physiol. 61, 687–708 (1973).4541078 10.1085/jgp.61.6.687
33. Eisenman LN He Y Fields C Zorumski CF Mennerick S Activation-dependent properties of pregnenolone sulfate inhibition of GABAA receptor-mediated current J. Physiol. 2003 550 679 691 10.1113/jphysiol.2003.043810 12938673
Eisenman, L. N., He, Y., Fields, C., Zorumski, C. F. & Mennerick, S. Activation-dependent properties of pregnenolone sulfate inhibition of GABAA receptor-mediated current. J. Physiol. 550, 679–691 (2003).12938673 10.1113/jphysiol.2003.043810
34. Mukherjee J Estradiol modulates the efficacy of synaptic inhibition by decreasing the dwell time of GABAA receptors at inhibitory synapses Proc. Natl. Acad. Sci. USA 2017 114 11763 11768 10.1073/pnas.1705075114 29078280
Mukherjee, J. et al. Estradiol modulates the efficacy of synaptic inhibition by decreasing the dwell time of GABAA receptors at inhibitory synapses. Proc. Natl. Acad. Sci. USA 114, 11763–11768 (2017).29078280 10.1073/pnas.1705075114
35. Townsend EA Miller VM Prakash YS Sex differences and sex steroids in lung health and disease Endocr. Rev. 2012 33 1 47 10.1210/er.2010-0031 22240244
Townsend, E. A., Miller, V. M. & Prakash, Y. S. Sex differences and sex steroids in lung health and disease. Endocr. Rev. 33, 1–47 (2012).22240244 10.1210/er.2010-0031
36. Dalal V Lipid nanodisc scaffold and size alter the structure of a pentameric ligand-gated ion channel Nat. Commun. 2024 15 25 10.1038/s41467-023-44366-w 38167383
Dalal, V. et al. Lipid nanodisc scaffold and size alter the structure of a pentameric ligand-gated ion channel. Nat. Commun. 15, 25 (2024).38167383 10.1038/s41467-023-44366-w
37. Lyons JA Bøggild A Nissen P Frauenfeld J Saposin-lipoprotein scaffolds for structure determination of membrane transporters Methods Enzymol. 2017 594 85 99 10.1016/bs.mie.2017.06.035 28779844
Lyons, J. A., Bøggild, A., Nissen, P. & Frauenfeld, J. Saposin-lipoprotein scaffolds for structure determination of membrane transporters. Methods Enzymol. 594, 85–99 (2017).28779844 10.1016/bs.mie.2017.06.035
38. Zivanov J New tools for automated high-resolution cryo-EM structure determination in RELION-3 ELife 2018 7 e42166 10.7554/eLife.42166 30412051
Zivanov, J. et al. New tools for automated high-resolution cryo-EM structure determination in RELION-3. ELife 7, e42166 (2018).30412051 10.7554/eLife.42166
39. Grant, T., Rohou, A. & Grigorieff, N. cisTEM, user-friendly software for single-particle image processing. ELife 7, 10.7554/eLife.35383 (2018).
40. Bepler T Positive-unlabeled convolutional neural networks for particle picking in cryo-electron micrographs Nat. Methods 2019 16 1153 1160 10.1038/s41592-019-0575-8 31591578
Bepler, T. et al. Positive-unlabeled convolutional neural networks for particle picking in cryo-electron micrographs. Nat. Methods 16, 1153–1160 (2019).31591578 10.1038/s41592-019-0575-8
41. Punjani A Rubinstein JL Fleet DJ Brubaker M A. cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination Nat. Methods 2017 14 290 296 10.1038/nmeth.4169 28165473
Punjani, A., Rubinstein, J. L., Fleet, D. J. & Brubaker, M. A. cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination. Nat. Methods 14, 290–296 (2017).28165473 10.1038/nmeth.4169
42. Punjani A Zhang H Fleet DJ Non-uniform refinement: adaptive regularization improves single-particle cryo-EM reconstruction Nat. Methods 2020 17 1214 1221 10.1038/s41592-020-00990-8 33257830
Punjani, A., Zhang, H. & Fleet, D. J. Non-uniform refinement: adaptive regularization improves single-particle cryo-EM reconstruction. Nat. Methods 17, 1214–1221 (2020).33257830 10.1038/s41592-020-00990-8
43. Casañal A Lohkamp B Emsley P Current developments in Coot for macromolecular model building of electron cryo-microscopy and crystallographic data Protein Sci. Publ. Protein Soc. 2020 29 1069 1078 10.1002/pro.3791
Casañal, A., Lohkamp, B. & Emsley, P. Current developments in Coot for macromolecular model building of electron cryo-microscopy and crystallographic data. Protein Sci. Publ. Protein Soc. 29, 1069–1078 (2020).10.1002/pro.3791
44. Adams PD PHENIX: a comprehensive Python-based system for macromolecular structure solution Acta Crystallogr. D Biol. Crystallogr. 2010 66 213 221 10.1107/S0907444909052925 20124702
Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. D Biol. Crystallogr. 66, 213–221 (2010).20124702 10.1107/S0907444909052925
45. Williams CJ MolProbity: More and better reference data for improved all-atom structure validation Protein Sci. Publ. Protein Soc. 2018 27 293 315 10.1002/pro.3330
Williams, C. J. et al. MolProbity: More and better reference data for improved all-atom structure validation. Protein Sci. Publ. Protein Soc. 27, 293–315 (2018).10.1002/pro.3330
46. Klesse G Rao S Sansom MSP Tucker SJ CHAP: A versatile tool for the structural and functional annotation of ion channel pores J. Mol. Biol. 2019 431 3353 3365 10.1016/j.jmb.2019.06.003 31220459
Klesse, G., Rao, S., Sansom, M. S. P. & Tucker, S. J. CHAP: A versatile tool for the structural and functional annotation of ion channel pores. J. Mol. Biol. 431, 3353–3365 (2019).31220459 10.1016/j.jmb.2019.06.003
47. Pettersen EF UCSF ChimeraX: Structure visualization for researchers, educators, and developers Protein Sci. Publ. Protein Soc. 2021 30 70 82 10.1002/pro.3943
Pettersen, E. F. et al. UCSF ChimeraX: Structure visualization for researchers, educators, and developers. Protein Sci. Publ. Protein Soc. 30, 70–82 (2021).10.1002/pro.3943
48. Jo S Kim T Iyer VG Im W CHARMM-GUI: a web-based graphical user interface for CHARMM J. Comput. Chem. 2008 29 1859 1865 10.1002/jcc.20945 18351591
Jo, S., Kim, T., Iyer, V. G. & Im, W. CHARMM-GUI: a web-based graphical user interface for CHARMM. J. Comput. Chem. 29, 1859–1865 (2008).18351591 10.1002/jcc.20945
49. Ingólfsson HI Computational lipidomics of the neuronal plasma membrane Biophys. J. 2017 113 2271 2280 10.1016/j.bpj.2017.10.017 29113676
Ingólfsson, H. I. et al. Computational lipidomics of the neuronal plasma membrane. Biophys. J. 113, 2271–2280 (2017).29113676 10.1016/j.bpj.2017.10.017
50. Huang J CHARMM36m: an improved force field for folded and intrinsically disordered proteins Nat. Methods 2017 14 71 73 10.1038/nmeth.4067 27819658
Huang, J. et al. CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nat. Methods 14, 71–73 (2017).27819658 10.1038/nmeth.4067
51. Vanommeslaeghe K CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields J. Comput. Chem. 2010 31 671 690 10.1002/jcc.21367 19575467
Vanommeslaeghe, K. et al. CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. J. Comput. Chem. 31, 671–690 (2010).19575467 10.1002/jcc.21367
52. Páll S Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS J. Chem. Phys. 2020 153 134110 10.1063/5.0018516 33032406
Páll, S. et al. Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS. J. Chem. Phys. 153, 134110 (2020).33032406 10.1063/5.0018516
53. Bussi G Donadio D Parrinello M Canonical sampling through velocity rescaling J. Chem. Phys. 2007 126 014101 10.1063/1.2408420 17212484
Bussi, G., Donadio, D. & Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 126, 014101 (2007).17212484 10.1063/1.2408420
54. Parrinello M Rahman A Crystal structure and pair potentials: A molecular-dynamics study Phys. Rev. Lett. 1980 45 1196 1199 10.1103/PhysRevLett.45.1196
Parrinello, M. & Rahman, A. Crystal structure and pair potentials: A molecular-dynamics study. Phys. Rev. Lett. 45, 1196–1199 (1980).10.1103/PhysRevLett.45.1196
55. Hess B P-LINCS: A parallel linear constraint solver for molecular simulation J. Chem. Theory Comput. 2008 4 116 122 10.1021/ct700200b 26619985
Hess, B. P-LINCS: A parallel linear constraint solver for molecular simulation. J. Chem. Theory Comput. 4, 116–122 (2008).26619985 10.1021/ct700200b
56. Essmann U A smooth particle mesh Ewald method J. Chem. Phys. 1995 103 8577 8593 10.1063/1.470117
Essmann, U. et al. A smooth particle mesh Ewald method. J. Chem. Phys. 103, 8577–8593 (1995).10.1063/1.470117
57. Gowers, R. et al. MDAnalysis: A Python package for the rapid analysis of molecular dynamics simulations. In 98–105 10.25080/Majora-629e541a-00e (Austin, Texas, 2016).
58. Humphrey W Dalke A Schulten K VMD: visual molecular dynamics J. Mol. Graph. 1996 14 27 28 10.1016/0263-7855(96)00018-5
Humphrey, W., Dalke, A. & Schulten, K. VMD: visual molecular dynamics. J. Mol. Graph. 14, 27–28 (1996).10.1016/0263-7855(96)00018-5
59. Hunter JD Matplotlib: A 2D graphics environment Comput. Sci. Eng. 2007 9 90 95 10.1109/MCSE.2007.55
Hunter, J. D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 9, 90–95 (2007).10.1109/MCSE.2007.55
60. Lindahl V Gourdon P Andersson M Hess B Permeability and ammonia selectivity in aquaporin TIP2;1: linking structure to function Sci. Rep. 2018 8 2995 10.1038/s41598-018-21357-2 29445244
Lindahl, V., Gourdon, P., Andersson, M. & Hess, B. Permeability and ammonia selectivity in aquaporin TIP2;1: linking structure to function. Sci. Rep. 8, 2995 (2018).29445244 10.1038/s41598-018-21357-2
