Modulation of the GABA type A receptor (GABAAR) function by cholesterol and other steroids is documented at the functional level, yet its structural basis is largely unknown. Current data on structurally related modulators suggest that cholesterol binds to subunit interfaces between transmembrane domains of the GABAAR. We construct homology models of a human GABAAR based on the structure of the glutamate-gated chloride channel GluCl of Caenorhabditis elegans. The models show the possibility of previously unreported disulfide bridges linking the M1 and M3 transmembrane helices in the α and γ subunits. We discuss the biological relevance of such disulfide bridges. Using our models, we investigate cholesterol binding to intersubunit cavities of the GABAAR transmembrane domain. We find that very similar binding modes are predicted independently by three approaches: analogy with ivermectin in the GluCl crystal structure, automated docking by AutoDock, and spontaneous rebinding events in unbiased molecular dynamics simulations. Taken together, the models and atomistic simulations suggest a somewhat flexible binding mode, with several possible orientations. Finally, we explore the possibility that cholesterol promotes pore opening through a wedge mechanism.

Highlights d b-NADH binding to hVDAC-1 does not change conformation of N-terminal helices d b-NADH binding leads to low VDAC conductance by sterically blocking anion flux d Elevated dynamics of the N-terminal helices are present in the ground state d Movements of helix a2 are important for VDAC gating

Summary The gating of pentameric ligand gated ion channels is sensitive to a variety of allosteric modulators that act on structures that are peripheral to those involved in the allosteric pathway leading from the agonist site to the channel gate. One such structure, the lipid-exposed transmembrane α-helix, M4, is the target of lipids, neurosteroids, and disease-causing mutations. Here we show that M4 interactions with the adjacent transmembrane α-helices, M1 and M3, modulate pLGIC function. Enhanced M4 interactions promote, while ineffective interactions reduce channel function. The interface chemistry governs the intrinsic strength of M4-M1/M3 inter-helical interactions, both influencing channel gating and imparting distinct susceptibilities to the potentiating effects of a lipid-facing M4 congenital myasthenic syndrome mutation. Through aromatic substitutions, functional studies, and molecular dynamics simulations, we elucidate a mechanism by which M4 modulates channel function.

Propofol, an intravenous anesthetic, is a positive modulator of the GABAA receptor, but the mechanistic details, including the relevant binding sites and alternative targets, remain disputed. Here we undertook an in-depth study of alkylphenol-based anesthetic binding to synaptic membranes. We designed, synthesized, and characterized a chemically active alkylphenol anesthetic (2-((prop-2-yn-1-yloxy)methyl)-5-(3-(trifluoromethyl)-3H-diazirin-3-yl)phenol, AziPm-click (1)), for affinity-based protein profiling (ABPP) of propofol-binding proteins in their native state within mouse synaptosomes. The ABPP strategy captured ∼4% of the synaptosomal proteome, including the unbiased capture of five α or β GABAA receptor subunits. Lack of γ2 subunit capture was not due to low abundance. Consistent with this, independent molecular dynamics simulations with alchemical free energy perturbation calculations predicted selective propofol binding to interfacial sites, with higher affinities for α/β than γ-containing interfaces. The simulations indicated hydrogen bonding is a key component leading to propofol-selective binding within GABAA receptor subunit interfaces, with stable hydrogen bonds observed between propofol and α/β cavity residues but not γ cavity residues. We confirmed this by introducing a hydrogen bond-null propofol analogue as a protecting ligand for targeted-ABPP and observed a lack of GABAA receptor subunit protection. This investigation demonstrates striking interfacial GABAA receptor subunit selectivity in the native milieu, suggesting that asymmetric occupancy of heteropentameric ion channels by alkylphenol-based anesthetics is sufficient to induce modulation of activity.

The COVID-19 pandemic has led to unprecedented efforts to identify drugs that can reduce its associated morbidity/mortality rate. Computational chemistry approaches hold the potential for triaging potential candidates far more quickly than their experimental counterparts. These methods have been widely used to search for small molecules that can inhibit critical proteins involved in the SARS-CoV-2 replication cycle. An important target is the SARS-CoV-2 main protease Mpro, an enzyme that cleaves the viral polyproteins into individual proteins required for viral replication and transcription. Unfortunately, standard computational screening methods face difficulties in ranking diverse ligands to a receptor due to disparate ligand scaffolds and varying charge states. Here, we describe full density functional quantum mechanical (DFT) simulations of Mpro in complex with various ligands to obtain absolute ligand binding energies. Our calculations are enabled by a new cloud-native parallel DFT implementation running on computational resources from Amazon Web Services (AWS). The results we obtain are promising: the approach is quite capable of scoring a very diverse set of existing drug compounds for their affinities to M pro and suggest the DFT approach is potentially more broadly applicable to repurpose screening against this target. In addition, each DFT simulation required only ~ 1 h (wall clock time) per ligand. The fast turnaround time raises the practical possibility of a broad application of large-scale quantum mechanics in the drug discovery pipeline at stages where ligand diversity is essential.