Plasticity of the selectivity filter is essential for permeation in lysosomal TPC2 channels.

Two-pore channels are pathophysiologically important Na+- and Ca2+-permeable channels expressed in lysosomes and other acidic organelles. Unlike most other ion channels, their permeability is malleable and ligand-tuned such that when gated by the signaling lipid PI(3,5)P2, they are more Na+-selective than when gated by the Ca2+ mobilizing messenger nicotinic acid adenine dinucleotide phosphate. However, the structural basis that underlies such plasticity and single-channel behavior more generally remains poorly understood. A recent Cryo-electron microscopy (cryo-EM) structure of TPC2 bound to PI(3,5)P2 in a proposed open-channel conformation provided an opportunity to address this via molecular dynamics (MD) simulation. To our surprise, simulations designed to compute conductance through this structure revealed almost no Na+ permeation events even at very high transmembrane voltages. However further MD simulations identified a spontaneous transition to a dramatically different conformation of the selectivity filter that involved expansion and a flip in the orientation of two core asparagine residues. This alternative filter conformation was remarkably stable and allowed Na+ to flow through the channel leading to a conductance estimate that was in very good agreement with direct single-channel measurements. Furthermore, this conformation was more permeable for Na+ over Ca2+. Our results have important ramifications not just for understanding the control of ion selectivity in TPC2 channels but also more broadly in terms of how ion channels discriminate ions.

Recovery from desensitization in GluA2 AMPA receptors is affected by a single mutation in the N-terminal domain interface.

AMPA-type ionotropic glutamate receptors (AMPARs) are central to various neurological processes, including memory and learning. They assemble as homo- or heterotetramers of GluA1, GluA2, GluA3, and GluA4 subunits, each consisting of an N-terminal domain (NTD), a ligand-binding domain, a transmembrane domain, and a C-terminal domain. While AMPAR gating is primarily controlled by reconfiguration in the ligand-binding domain layer, our study focuses on the NTDs, which also influence gating, yet the underlying mechanism remains enigmatic. In this investigation, we employ molecular dynamics simulations to evaluate the NTD interface strength in GluA1, GluA2, and NTD mutants GluA2-H229N and GluA1-N222H. Our findings reveal that GluA1 has a significantly weaker NTD interface than GluA2. The NTD interface of GluA2 can be weakened by a single point mutation in the NTD dimer-of-dimer interface, namely H229N, which renders GluA2 more GluA1-like. Electrophysiology recordings demonstrate that this mutation also leads to slower recovery from desensitization. Moreover, we observe that lowering the pH induces more splayed NTD states and enhances desensitization in GluA2. We hypothesized that H229 was responsible for this pH sensitivity; however, GluA2-H229N was also affected by pH, meaning that H229 is not solely responsible and that protons exert their effect across multiple domains of the AMPAR. In summary, our work unveils an allosteric connection between the NTD interface strength and AMPAR desensitization.

Exploring the influence of pore shape on conductance and permeation.

There are increasing numbers of ion channel structures featuring heteromeric subunit assembly, exemplified by synaptic α1ßB glycine and α4ß2 nicotinic receptors. These structures exhibit inherent pore asymmetry, but the relevance of this to function is unknown. Furthermore, molecular dynamics simulations performed on symmetrical homomeric channels often lead to thermal distortion whereby conformations of the resulting ensemble are also asymmetrical. When functionally annotating ion channels, researchers often rely on minimal constrictions determined via radius-profile calculations performed with computer programs, such as HOLE or CHAP, coupled with an assessment of pore hydrophobicity. However, such tools typically employ spherical probe particles, limiting their ability to accurately capture pore asymmetry. Here, we introduce an algorithm that employs ellipsoidal probe particles, enabling a more comprehensive representation of the pore geometry. Our analysis reveals that the use of nonspherical ellipsoids for pore characterization provides a more accurate and easily interpretable depiction of conductance. To quantify the implications of pore asymmetry on conductance, we systematically investigated carbon nanotubes with varying degrees of pore asymmetry as model systems. The conductance through these channels shows surprising effects that would otherwise not be predicted with spherical probes. The results have broad implications not only for the functional annotation of biological ion channels but also for the design of synthetic channel systems for use in areas such as water filtration. Furthermore, we make use of the more accurate characterization of channel pores to refine a physical conductance model to obtain a heuristic estimate for single-channel conductance. The code is freely available, obtainable as pip-installable python package and provided as a web service.

Automated Absolute Binding Free Energy Calculation Workflow for Drug Discovery.

Absolute binding free energies play a crucial role in drug development, particularly as part of the lead discovery process. In recent work, we showed how in silico predictions directly could support drug development by ranking and recommending favorable ideas over unfavorable ones. Here, we demonstrate a Python workflow that enables the calculation of ABFEs with minimal manual input effort, such as the receptor PDB and ligand SDF files, and outputs a .tsv file containing the ranked ligands and their corresponding binding free energies. The implementation uses Snakemake to structure and control the execution of tasks, allowing for dynamic control of parameters and execution patterns. We provide an example of a benchmark system that demonstrates the effectiveness of the automated workflow.

The mechanism of mammalian proton-coupled peptide transporters.

Proton-coupled oligopeptide transporters (POTs) are of great pharmaceutical interest owing to their promiscuous substrate binding site that has been linked to improved oral bioavailability of several classes of drugs. Members of the POT family are conserved across all phylogenetic kingdoms and function by coupling peptide uptake to the proton electrochemical gradient. Cryo-EM structures and alphafold models have recently provided new insights into different conformational states of two mammalian POTs, SLC15A1, and SLC15A2. Nevertheless, these studies leave open important questions regarding the mechanism of proton and substrate coupling, while simultaneously providing a unique opportunity to investigate these processes using molecular dynamics (MD) simulations. Here, we employ extensive unbiased and enhanced-sampling MD to map out the full SLC15A2 conformational cycle and its thermodynamic driving forces. By computing conformational free energy landscapes in different protonation states and in the absence or presence of peptide substrate, we identify a likely sequence of intermediate protonation steps that drive inward-directed alternating access. These simulations identify key differences in the extracellular gate between mammalian and bacterial POTs, which we validate experimentally in cell-based transport assays. Our results from constant-PH MD and absolute binding free energy (ABFE) calculations also establish a mechanistic link between proton binding and peptide recognition, revealing key details underpining secondary active transport in POTs. This study provides a vital step forward in understanding proton-coupled peptide and drug transport in mammals and pave the way to integrate knowledge of solute carrier structural biology with enhanced drug design to target tissue and organ bioavailability.

The cells in our body are sealed by a surrounding membrane that allows them to control which molecules can enter or leave. Desired molecules are often imported via transport proteins that require a source of energy. One way that transporter proteins achieve this is by simultaneously moving positively charged particles called protons across the membrane. Proteins called POTs (short for proton-coupled oligopeptide transporters) use this mechanism to import small peptides and drugsin to the cells of the kidney and small intestine. Sitting in the centre of these transporters is a pocket that binds to the imported peptide which has a gate on either side: an outer gate that opens towards the outside of the cell, and an inner gate that opens towards the cell's interior. The movement of protons from the outer to the inner gate is thought to shift the shape of the transporter from an outwards to an inwards-facing state. However, the molecular details of this energetic coupling are not well understood. To explore this, Lichtinger et al. used computer simulations to pinpoint where protons bind on POTs to trigger the gates to open. The simulations proposed that two sites together make up the outward-facing gate, which opens upon proton binding. Lichtinger et al. then validated these sites experimentally in cultured human cells that produce mutant POTs. After the desired peptide/drug has attached to the binding pocket, the protons then move to two more sites further down the transporter. This triggers the inner gate to open, which ultimately allows the small molecule to move into the cell. These findings represent a significant step towards understanding how POTs transport their cargo. Since POTs can transport a range of drugs from the digestive tract into the body, these results could help researchers design molecules that are better absorbed. This could lead to more orally available medications, making it easier for patients to adhere to their treatment regimen.

Cryo-EM structures of prokaryotic ligand-gated ion channel GLIC provide insights into gating in a lipid environment.

GLIC, a proton-activated prokaryotic ligand-gated ion channel, served as a model system for understanding the eukaryotic counterparts due to their structural and functional similarities. Despite extensive studies conducted on GLIC, the molecular mechanism of channel gating in the lipid environment requires further investigation. Here, we present the cryo-EM structures of nanodisc-reconstituted GLIC at neutral and acidic pH in the resolution range of 2.6 - 3.4 Å. In our apo state at pH 7.5, the extracellular domain (ECD) displays conformational variations compared to the existing apo structures. At pH 4.0, three distinct conformational states (C1, C2 and O states) are identified. The protonated structures exhibit a compacted and counter-clockwise rotated ECD compared with our apo state. A gradual widening of the pore in the TMD is observed upon reducing the pH, with the widest pore in O state, accompanied by several layers of water pentagons. The pore radius and molecular dynamics (MD) simulations suggest that the O state represents an open conductive state. We also observe state-dependent interactions between several lipids and proteins that may be involved in the regulation of channel gating. Our results provide comprehensive insights into the importance of lipids impact on gating.

Selective Potentiation of the (α4)3(ß2)2 Nicotinic Acetylcholine Receptor Response by NS9283 Analogues.

NS9283, 3-(3-pyridyl)-5-(3-cyanophenyl)-1,2,4-oxadiazole, is a selective positive allosteric modulator of (α4)3(ß2)2 nicotinic acetylcholine receptors (nAChRs). It has good subtype selective therapeutic potential afforded by its specific binding to the unique α4-α4 subunit interface present in the (α4)3(ß2)2 nAChR. However, there is currently a lack of structure activity relationship (SAR) studies aimed at developing a class of congeners endowed with the same profile of activity that can help consolidate the druggability of the α4-α4 subunit interface. In this study, new NS9283 analogues were designed, synthesized, and characterized for their ability to selectively potentiate the ACh activity at heterologous (α4)3(ß2)2 nAChRs vs nAChR subtypes (α4)2(ß2)3, α5α4ß2, and α7. With few exceptions, all the NS9283 analogues exerted positive modulation of the (α4)3(ß2)2 nAChR ACh-evoked responses. Above all, those modified at the 3-cyanophenyl moiety by replacement with 3-nitrophenyl (4), 4-cyanophenyl (10), and N-formyl-4-piperidinyl (20) showed the same efficacy as NS9283, although with lower potency. Molecular dynamics simulations of NS9283 and some selected analogues highlighted consistency between potentiation activity and pose of the ligand inside the α4-α4 site with the main interaction being with the complementary (-) side and induction of a significant conformational change of the Trp156 residue in the principal (+) side.

Structural basis for antibiotic transport and inhibition in PepT2, the mammalian proton-coupled peptide transporter.

The uptake and elimination of beta-lactam antibiotics in the human body are facilitated by the proton-coupled peptide transporters PepT1 (SLC15A1) and PepT2 (SLC15A2). The mechanism by which SLC15 family transporters recognize and discriminate between different drug classes and dietary peptides remains unclear, hampering efforts to improve antibiotic pharmacokinetics through targeted drug design and delivery. Here, we present cryo-EM structures of the mammalian proton-coupled peptide transporter, PepT2, in complex with the widely used beta-lactam antibiotics cefadroxil, amoxicillin and cloxacillin. Our structures, combined with pharmacophore mapping, molecular dynamics simulations and biochemical assays, establish the mechanism of antibiotic recognition and the important role of protonation in drug binding and transport.

Molecular basis for pH sensing in the KDEL trafficking receptor.

Trafficking receptors control protein localization through the recognition of specific signal sequences that specify unique cellular locations. Differences in luminal pH are important for the vectorial trafficking of cargo receptors. The KDEL receptor is responsible for maintaining the integrity of the ER by retrieving luminally localized folding chaperones in a pH-dependent mechanism. Structural studies have revealed the end states of KDEL receptor activation and the mechanism of selective cargo binding. However, precisely how the KDEL receptor responds to changes in luminal pH remains unclear. To explain the mechanism of pH sensing, we combine analysis of X-ray crystal structures of the KDEL receptor at neutral and acidic pH with advanced computational methods and cell-based assays. We show a critical role for ordered water molecules that allows us to infer a direct connection between protonation in different cellular compartments and the consequent changes in the affinity of the receptor for cargo.

Small-molecule positive allosteric modulation of homomeric kainate receptors GluK1-3: development of screening assays and insight into GluK3 structure.

The kainate receptors GluK1-3 (glutamate receptor ionotropic, kainate receptors 1-3) belong to the family of ionotropic glutamate receptors and are essential for fast excitatory neurotransmission in the brain, and are associated with neurological and psychiatric diseases. How these receptors can be modulated by small-molecule agents is not well understood, especially for GluK3. We show that the positive allosteric modulator BPAM344 can be used to establish robust calcium-sensitive fluorescence-based assays to test agonists, antagonists, and positive allosteric modulators of GluK1-3. The half-maximal effective concentration (EC50) of BPAM344 for potentiating the response of 100 µm kainate was determined to be 26.3 µm for GluK1, 75.4 µm for GluK2, and 639 µm for GluK3. Domoate was found to be a potent agonist for GluK1 and GluK2, with an EC50 of 0.77 and 1.33 µm, respectively, upon co-application of 150 µm BPAM344. At GluK3, domoate acts as a very weak agonist or antagonist with a half-maximal inhibitory concentration (IC50) of 14.5 µm, in presence of 500 µm BPAM344 and 100 µm kainate for competition binding. Using H523A-mutated GluK3, we determined the first dimeric structure of the ligand-binding domain by X-ray crystallography, allowing location of BPAM344, as well as zinc-, sodium-, and chloride-ion binding sites at the dimer interface. Molecular dynamics simulations support the stability of the ion sites as well as the involvement of Asp761, Asp790, and Glu797 in the binding of zinc ions. Using electron microscopy, we show that, in presence of glutamate and BPAM344, full-length GluK3 adopts a dimer-of-dimers arrangement.