Ligand Binding Site Research Articles

Insights Into Molecular Interactions and Biological Effect of Natural Stilbenoids at The TRPA1 Ion Channel.

Natural stilbenoids, polyphenolic compounds notably found in Scots pine and Norway spruce, have been shown to exhibit analgesic and anti-inflammatory effects through the TRPA1 channel, making them promising hits for the development of novel agents to treat inflammatory diseases and pain. In this study, we computationally investigated the putative binding sites of natural stilbenoids at the TRPA1 channel. Specifically, we employed molecular docking and MD simulation approaches to explore three known ligand binding sites at TRPA1. Furthermore, the biological effect of the studied compounds on TRPA1 was assessed in vitro using a fluorescent imaging plate reader (FLIPR™) calcium assay. Our modeling results suggest the stilbenoids exhibit higher affinity to the two agonist binding sites than the antagonistic site. Consistent with this, the in vitro results showed that the stilbenoids act as moderate TRPA1 channel agonists and likely inhibit the channel through a desensitization mechanism rather than act as pure TRPA1 antagonists. Additionally, our bias-force pulling simulations proposed an additional binding pocket for the natural stilbenoids that is distinct from the known ligand binding sites at TRPA1. The results of the study give useful insights into structure-based design and development of novel therapeutic TRPA1 modulators.

Exploration of Cryptic Pockets Using Enhanced Sampling Along Normal Modes: A Case Study of KRAS G12D.

Identification of cryptic pockets has the potential to open new therapeutic opportunities by discovering ligand binding sites that remain hidden in static apo structures of a target protein. Moreover, allosteric cryptic pockets can become valuable for designing target-selective ligands when the natural ligand binding sites are conserved in variants of a protein. For example, before an allosteric cryptic pocket was discovered, KRAS was considered undruggable due to its smooth surface and conservation of the GDP/GTP binding pocket across the wild type and oncogenic isoforms. Recent identification of the Switch-II cryptic pocket in the KRASG12C mutant and FDA approval of anticancer drugs targeting this site underscores the importance of cryptic pockets in solving pharmaceutical challenges. Here, we present a newly developed approach for the exploration of cryptic pockets using weighted ensemble molecular dynamics simulations with inherent normal modes as progress coordinates applied to the wild type KRAS and the G12D mutant. We performed extensive all-atomic simulations (>400 μs) with and without several cosolvents (xenon, ethanol, benzene), and analyzed trajectories using three distinct methods to search for potential binding pockets. These methods have been applied as a proof-of-concept to KRAS and have shown they can predict known cryptic binding sites. Furthermore, we performed ligand-binding simulations of a known inhibitor (MRTX1133) to shed light on the nature of cryptic pockets in KRASG12D and the role of conformational selection vs induced-fit mechanism in the formation of these cryptic pockets.

Molecular recognition of peptides and proteins by cucurbit[n]urils: systems and applications.

The development of methodology for attaching ligand binding sites to proteins of interest has accelerated biomedical science. Such protein tags have widespread applications as well as properties that significantly limit their utility. This review describes the mechanisms and applications of supramolecular systems comprising the synthetic receptors cucurbit[7]uril (Q7) or cucurbit[8]uril (Q8) and their polypeptide ligands. Molecular recognition of peptides and proteins occurs at sites of 1-3 amino acids with high selectivity and affinity via several distinct mechanisms, which are supported by extensive thermodynamic and structural studies in aqueous media. The commercial availability, low cost, high stability, and biocompatibility of these synthetic receptors has led to the development of myriad applications. This comprehensive review compiles the molecular recognition studies and the resulting applications with the goals of providing a valuable resource to the community and inspiring the next generation of innovation.

Cryo-EM investigation of ryanodine receptor type 3.

Ryanodine Receptor isoform 3 (RyR3) is a large ion channel found in the endoplasmic reticulum membrane of many different cell types. Within the hippocampal region of the brain, it is found in dendritic spines and regulates synaptic plasticity. It controls myogenic tone in arteries and is upregulated in skeletal muscle in early development. RyR3 has a unique functional profile with a very high sensitivity to activating ligands, enabling high gain in Ca2+-induced Ca2+ release. Here we solve high-resolution cryo-EM structures of RyR3 in non-activating and activating conditions, revealing structural transitions that occur during channel opening. Addition of activating ligands yields only open channels, indicating an intrinsically high open probability under these conditions. RyR3 has reduced binding affinity to the auxiliary protein FKBP12.6 due to several sequence variations in the binding interface. We map disease-associated sequence variants and binding sites for known pharmacological agents. The N-terminal region contains ligand binding sites for a putative chloride anion and ATP, both of which are targeted by sequence variants linked to epileptic encephalopathy.

Self-Assembled Antibody-Oligonucleotide Conjugates for Targeted Delivery of Complementary Antisense Oligonucleotides.

Antibody-oligonucleotide conjugate (AOC) affords preferential cell targeting and enhanced cellular uptake of antisense oligonucleotide (ASO). Here, we have developed a modular AOC (MAOC) approach based on accurate self-assembly of separately prepared antibody and ASO modules. Homogeneous multimeric AOC with defined ASO-to-antibody ratio were generated by L-DNA scaffold mediated precise self-assembly of antibodies and ASOs. The MAOC approach has been implemented to deliver exon skipping ASOs via transferrin receptor (TfR1) mediated internalization. We discovered an anti-TfR1 sdAb that can greatly enhance nuclear delivery of ASOs. Cryo-EM structure of the sdAb-TfR1 complex showed a new epitope that does not overlap with the binding sites of endogenous TfR1 ligands. In vivo functional analyses of MAOCs with one ASO for single exon skipping and two ASOs for double exon skipping showed that both ASO concentration and exon skipping efficacy of MAOC in cardiac and skeletal muscles are dramatically higher than conventional ASOs in the transgenic Duchenne muscular dystrophy (DMD) mouse model. MAOC treatment was well tolerated in vivo and not associated with any toxicity-related morbidity or mortality. Collectively, our data suggest that the self-assembled MAOC is a viable option for broadening the therapeutic application of ASO via multi-specific targeting and delivery.

Prediction of protein-ligand binding sites modulating activity of MAST protein kinases

Aim. Identification of the protein-ligand binding sites, that may be the target of compounds, affecting individual human protein kinases of the MAST family (MAST1, 2, 3, 4 and MASTL / GWL). Methods. Literature and database search. Comparison of protein and ligand structures. Protein structure modeling, structural superimposition, etc. Results. The structural alignment demonstrates significant similarity of catalytic domains in MAST1, 2, 3, 4 and MASTL (GWL). It justifies transferring of reference ligands from PDB structures to human MASTs, discovering potential sites of ligand binding. 13 sites of ligand-binding were specified based on refrence ligands, transferred from RCSB Protein Data Bank structures, and differences in sites amino acid composition of MAST family members were discovered. Сonclusions. Based on the differences in the amino acid composition of the studied pockets in MAST1, 2, 3, 4 and MASTL (GWL), the sites B, C, D, E, F, were selected for further study and virtual screening for new selective inhibitors of individual members of MAST protein kinase family.

Structural features of carrot α-tubulin predetermining the natural resistance to dinitroaniline herbicides

Aim. To explain the natural resistance of Daucus carota L. to dinitroaniline herbicides. To clarify features of the carrot α-tubulin that may affect formation of the ligand-protein complex based on the structural and electrostatic analysis of the ligand-binding site. Methods. Reconstruction of the spatial structure of α-tubulin from D. carota and Toxoplasma gondii using profile (Swiss-Model) and de-novo (AlphaFold2) modeling. Molecular dynamics (MD) simulations of the built 3D-models in Gromacs. Analysis of the molecular electrostatics with PDB2PQR/APBS tools. Visualization and analysis of molecular structures in PyMOL. Results. It has been shown that along with the typical positive charge of the dinitroaniline-binding pocket, all isotypes of carrot α-tubulin demonstrate negatively charged regions that may cause conflicts with the nitro-groups of the ligands. Also, the MD-stable negatively charged "bridge" between Cys316 and the aryl-fragment of Phe255 was observed in all α-tubulin isotypes. In our opinion, it not only competes with the cyclic fragment of dinitroanilines, but overall prevent the opening of the site pocket in carrot α-tubulin. Conclusions. It was clarified that natural resistance of D. carota to dinitroaniline herbicides may be associated with steric and electrostatic conflicts observed in in all α-tubulin isotypes. In our opinion, it prevents interaction with dinitroaniline compounds at the stage of primary site recognition on the early stages of protein-ligand complex formation.

Mut-Map: Comprehensive Computational Pipeline for Structural Mapping and Analysis of Cancer-Associated Mutations.

Understanding the functional impact of genetic mutations on protein structures is essential for advancing cancer research and developing targeted therapies. The main challenge lies in accurately mapping these mutations to protein structures and analysing their effects on protein function. To address this, Mut-Map (https://genemutation.org/) is a comprehensive computational pipeline designed to integrate mutation data from the Catalogue Of Somatic Mutations In Cancer database with protein structural data from the Protein Data Bank and AlphaFold models. The pipeline begins by taking a UniProt ID and proceeds through mapping corresponding Protein Data Bank structures, renumbering residues, and assessing disorder percentages. It then overlays mutation data, categorizes mutations based on structural context, and visualizes them using advanced tools like MolStar. This approach allows for a detailed analysis of how mutations may disrupt protein function by affecting key regions such as DNA interfaces, ligand-binding sites, and dimer interactions. To validate the pipeline, a case study on the TP53 gene, a critical tumour suppressor often mutated in cancers, was conducted. The analysis highlighted the most frequent mutations occurring at the DNA-binding interface, providing insights into their potential role in cancer progression. Mut-Map offers a powerful resource for elucidating the structural implications of cancer-associated mutations, paving the way for more targeted therapeutic strategies and advancing our understanding of protein structure-function relationships.

Combined Physics- and Machine-Learning-Based Method to Identify Druggable Binding Sites Using SILCS-Hotspots.

Identifying druggable binding sites on proteins is an important and challenging problem, particularly for cryptic, allosteric binding sites that may not be obvious from X-ray, cryo-EM, or predicted structures. The Site-Identification by Ligand Competitive Saturation (SILCS) method accounts for the flexibility of the target protein using all-atom molecular simulations that include various small molecule solutes in aqueous solution. During the simulations, the combination of protein flexibility and comprehensive sampling of the water and solute spatial distributions can identify buried binding pockets absent in experimentally determined structures. Previously, we reported a method for leveraging the information in the SILCS sampling to identify binding sites (termed Hotspots) of small mono- or bicyclic compounds, a subset of which coincide with known binding sites of drug-like molecules. Here, we build on that physics-based approach and present a ML model for ranking the Hotspots according to the likelihood they can accommodate drug-like molecules (e.g., molecular weight >200 Da). In the independent validation set, which includes various enzymes and receptors, our model recalls 67% and 89% of experimentally validated ligand binding sites in the top 10 and 20 ranked Hotspots, respectively. Furthermore, we show that the model's output Decision Function is a useful metric to predict binding sites and their potential druggability in new targets. Given the utility the SILCS method for ligand discovery and optimization, the tools presented represent an important advancement in the identification of orthosteric and allosteric binding sites and the discovery of drug-like molecules targeting those sites.

Allosteric modulation of the CXCR4:CXCL12 axis by targeting receptor nanoclustering via the TMV-TMVI domain

CXCR4 is a ubiquitously expressed chemokine receptor that regulates leukocyte trafficking and arrest in both homeostatic and pathological states. It also participates in organogenesis, HIV-1 infection, and tumor development. Despite the potential therapeutic benefit of CXCR4 antagonists, only one, plerixafor (AMD3100), which blocks the ligand-binding site, has reached the clinic. Recent advances in imaging and biophysical techniques have provided a richer understanding of the membrane organization and dynamics of this receptor. Activation of CXCR4 by CXCL12 reduces the number of CXCR4 monomers/dimers at the cell membrane and increases the formation of large nanoclusters, which are largely immobile and are required for correct cell orientation to chemoattractant gradients. Mechanistically, CXCR4 activation involves a structural motif defined by residues in TMV and TMVI. Using this structural motif as a template, we performed in silico molecular modeling followed by in vitro screening of a small compound library to identify negative allosteric modulators of CXCR4 that do not affect CXCL12 binding. We identified AGR1.137, a small molecule that abolishes CXCL12-mediated receptor nanoclustering and dynamics and blocks the ability of cells to sense CXCL12 gradients both in vitro and in vivo while preserving ligand binding and receptor internalization.

Allosteric modulation of the CXCR4:CXCL12 axis by targeting receptor nanoclustering via the TMV-TMVI domain

CXCR4 is a ubiquitously expressed chemokine receptor that regulates leukocyte trafficking and arrest in both homeostatic and pathological states. It also participates in organogenesis, HIV-1 infection, and tumor development. Despite the potential therapeutic benefit of CXCR4 antagonists, only one, plerixafor (AMD3100), which blocks the ligand-binding site, has reached the clinic. Recent advances in imaging and biophysical techniques have provided a richer understanding of the membrane organization and dynamics of this receptor. Activation of CXCR4 by CXCL12 reduces the number of CXCR4 monomers/dimers at the cell membrane and increases the formation of large nanoclusters, which are largely immobile and are required for correct cell orientation to chemoattractant gradients. Mechanistically, CXCR4 activation involves a structural motif defined by residues in TMV and TMVI. Using this structural motif as a template, we performed in silico molecular modeling followed by in vitro screening of a small compound library to identify negative allosteric modulators of CXCR4 that do not affect CXCL12 binding. We identified AGR1.137, a small molecule that abolishes CXCL12-mediated receptor nanoclustering and dynamics and blocks the ability of cells to sense CXCL12 gradients both in vitro and in vivo while preserving ligand binding and receptor internalization.

Aminomethyl Salicylaldehydes Lock onto a Surface Lysine by Forming an Extended Intramolecular Hydrogen Bond Network.

The development of electrophilic ligands that rapidly modify specific lysine residues remains a major challenge. Salicylaldehyde-based inhibitors have been reported to form stable imine adducts with the catalytic lysine of protein kinases. However, the targeted lysine in these examples is buried in a hydrophobic environment. A key unanswered question is whether this strategy can be applied to a lysine on the surface of a protein, where rapid hydrolysis of the resulting salicylaldimine is more likely. Here, we describe a series of aminomethyl-substituted salicylaldehydes that target a fully solvated lysine on the surface of the ATPase domain of Hsp90. By systematically varying the orientation of the salicylaldehyde, we discovered ligands with long residence times, the best of which engages Hsp90 in a quasi-irreversible manner. Crystallographic analysis revealed a daisy-chain network of intramolecular hydrogen bonds in which the salicylaldimine is locked into position by the adjacent piperidine linker. This study highlights the potential of aminomethyl salicylaldehydes to generate conformationally stabilized, hydrolysis-resistant imines, even when the targeted lysine is far from the ligand binding site and is exposed to bulk solvent.

Pulse EPR Methods in the Angstrom to Nanometre Scale Shed Light on the Conformational Flexibility of a Fluoride Riboswitch.

Riboswitches control gene regulation upon external stimuli such as environmental factors or ligand binding. The fluoride sensing riboswitch from Thermotoga petrophila is a complex regulatory RNA proposed to be involved in resistance to F- cytotoxicity. The details of structure and dynamics underpinning the regulatory mechanism are currently debated. Here we demonstrate that a combination of pulsed electron paramagnetic resonance (ESR/EPR) spectroscopies, detecting distances in the angstrom to nanometre range, can probe distinct regions of conformational flexibility in this riboswitch. PELDOR (pulsed electron-electron double resonance) revealed a similar preorganisation of the sensing domain in three forms, i.e. the free aptamer, the Mg2+-bound apo, and the F--bound holo form. 19F ENDOR (electron-nuclear double resonance) was used to investigate the active site structure of the F--bound holo form. Distance distributions without a priori structural information were compared with in silico modelling of spin label conformations based on the crystal structure. While PELDOR, probing the periphery of the RNA fold, revealed conformational flexibility of the RNA backbone, ENDOR indicated low structural heterogeneity at the ligand binding site. Overall, the combination of PELDOR and ENDOR with sub-angstrom precision gave insight into structural organisation and flexibility of a riboswitch, not easily attainable by other biophysical techniques.

Pulse EPR Methods in the Angstrom to Nanometre Scale Shed Light on the Conformational Flexibility of a Fluoride Riboswitch

Riboswitches control gene regulation upon external stimuli such as environmental factors or ligand binding. The fluoride sensing riboswitch from Thermotoga petrophila is a complex regulatory RNA proposed to be involved in resistance to F‐ cytotoxicity. The details of structure and dynamics underpinning the regulatory mechanism are currently debated. Here we demonstrate that a combination of pulsed electron paramagnetic resonance (ESR/EPR) spectroscopies, detecting distances in the angstrom to nanometre range, can probe distinct regions of conformational flexibility in this riboswitch. PELDOR (pulsed electron‐electron double resonance) revealed a similar preorganisation of the sensing domain in three forms, i.e. the free aptamer, the Mg2+‐bound apo, and the F‐‐bound holo form. 19F ENDOR (electron‐nuclear double resonance) was used to investigate the active site structure of the F‐‐bound holo form. Distance distributions without a priori structural information were compared with in silico modelling of spin label conformations based on the crystal structure. While PELDOR, probing the periphery of the RNA fold, revealed conformational flexibility of the RNA backbone, ENDOR indicated low structural heterogeneity at the ligand binding site. Overall, the combination of PELDOR and ENDOR with sub‐angstrom precision gave insight into structural organisation and flexibility of a riboswitch, not easily attainable by other biophysical techniques.

D1-Tyr246 and D2-Tyr244 in photosystem II: Insights into bicarbonate binding and electron transfer from QA•− to QB

In photosystem II (PSII), D1-Tyr246 and D2-Tyr244 are symmetrically located at the binding site of the bicarbonate ligand of the non-heme Fe complex. Here, we investigated the role of the symmetrically arranged tyrosine pair, D1-Tyr246 and D2-Tyr244, in the function of PSII, by generating four chloroplast mutants of PSII from Chlamydomonas reinhardtii: D1-Y246F, D1-Y246T, D2-Y244F, and D2-Y244T. The mutants exhibited altered photoautotrophic growth, reduced PSII protein accumulation, and impaired O2-evolving activity. Flash-induced fluorescence yield decay kinetics indicated a significant slowdown in electron transfer from QA•− to QB in all mutants. Bicarbonate reconstitution resulted in enhanced O2-evolving activity, suggesting destabilization of bicarbonate binding in the mutants. Structural analyses based on a quantum mechanical/molecular mechanical approach identified the existence of a water channel that leads to incorporation of bulk water molecules and destabilization of the bicarbonate binding site. The water intake channels, crucial for bicarbonate stability, exhibited distinct paths in the mutants. These findings shed light on the essential role of the tyrosine pair in maintaining bicarbonate stability and facilitating efficient electron transfer in native PSII.

Hormonal and Allosteric Regulation of the Luteinizing Hormone/Chorionic Gonadotropin Receptor.

Luteinizing hormone (LH) and human chorionic gonadotropin (CG), like follicle-stimulating hormone, are the most important regulators of the reproductive system. They exert their effect on the cell through the LH/CG receptor (LHCGR), which belongs to the family of G protein-coupled receptors. Binding to gonadotropin induces the interaction of LHCGR with various types of heterotrimeric G proteins (Gs, Gq/11, Gi) and β-arrestins, which leads to stimulation (Gs) or inhibition (Gi) of cyclic adenosine monophosphate-dependent cascades, activation of the phospholipase pathway (Gq/11), and also to the formation of signalosomes that mediate the stimulation of mitogen-activated protein kinases (β-arrestins). The efficiency and selectivity of activation of intracellular cascades by different gonadotropins varies, which is due to differences in their interaction with the ligand-binding site of LHCGR. Gonadotropin signaling largely depends on the status of N- and O-glycosylation of LH and CG, on the formation of homo- and heterodimeric receptor complexes, on the cell-specific microenvironment of LHCGR and the presence of autoantibodies to it, and allosteric mechanisms are important in the implementation of these influences, which is due to the multiplicity of allosteric sites in different loci of the LHCGR. The development of low-molecular-weight allosteric regulators of LHCGR with different profiles of pharmacological activity, which can be used in medicine for the correction of reproductive disorders and in assisted reproductive technologies, is promising. These and other issues regarding the hormonal and allosteric regulation of LHCGR are summarized and discussed in this review.

Comprehensive characterization of a novel small-molecule activator for the nuclear receptor Nur77: Chemical, molecular, and biological insights

The nuclear receptor Nur77 is a transcription factor belonging to the NR4A subfamily. Upon activation, it regulates a wide array of biological and pathophysiological processes by modulating the expression of its target genes. Previous findings have classified Nur77 as an orphan receptor because of the discovery of a structurally atypical ligand-binding domain and the lack of identification of an endogenous ligand. Nevertheless, recent studies have uncovered several endogenous, natural, and small synthetic molecules that can bind to and activate Nur77. However, developing selective and potent Nur77 activators remains a significant challenge. The discovery of novel and potential small synthetic molecules that modulate Nur77 activity will facilitate therapeutic interventions of Nur77 against several human diseases. In this study, we have reported the development of a novel and effective Nur77 ligand. Our data show that (1E,4E)-1,5-di(pyrazin-2-yl)penta-1,4-dien-3-one (PB) induces Nur77 transcriptional activity by interacting with a putative Nur77 ligand binding site by forming solid hydrogen bonding. Calculated chemical parameters denote that PB has sophisticated chemical features that will enhance its interaction with the Nur77 ligand-binding domain. Molecular docking simulations showed that PB fits in the Nur77 putative ligand binding site with solid hydrogen bonding, and molecular dynamics simulations indicate that PB binding would stabilize the Nur77 ligand binding domain. Further, in vitro studies revealed that PB increases Nur77 nuclear expression and activity, inhibits cigarette smoke-induced inflammatory phenotype of airway epithelial cells, and protects against apoptosis. These findings provide insights into developing an effective Nur77 small-molecule activator which could be developed into a therapeutic agent against inflammatory diseases.

Promoted electrochemical performance by MOF on MOF composite catalyst of microbial fuel cell: CuCo-MOF@ZIF-8 and the comparison between two-step hydrothermal method and dual-solution method

The microbial fuel cell (MFC) is a device that simultaneously achieves electricity generation and sewage degradation. In this study, a novel cathode catalyst metal-organic frameworks (MOFs) have been fabricated by two-step hydrothermal and dual-solution method (CuCo-MOF@ZIF-8). The synthesized trimetal MOFs exhibited a 3D badminton-like structure morphology and porosity. The results of the characterizations showed that CuCo-MOF@ZIF-8 possesses greater surface area porosity and novel functional groups. The Trimetal MOF-on-MOF mode not only demonstrated the stability of the structure but also enhanced its mechanism. Molecular mechanism analysis revealed changes in the organic ligand and metal binding site due to the transformation of Cu2+ to Cu+, Co2+ to Co3+, and Zn–N to Zn–O organic connection. Furthermore, differences between the two fabrication methods were compared. The solid-state single preparation (CuCo-MOF@ZIF-8-1), was synthesized using the two-step hydrothermal method; the liquid mixed preparation material (CuCo-MOF@ZIF-8-2), was synthesized using the dual-solution method; they exhibited completely different chemical structures and morphologies during material testing and characterization. The maximum output power density of CuCo-MOF@ZIF-8-2-MFC was 246.38 mW/m2, about 2.49 times of ZIF-8 (98.72 mW/m2). The output voltage of CuCo-MOF@ZIF-8-1-MFC was measured at 357 mV over 10 d, while that of CuCo-MOF@ZIF-8-2-MFC reached 365 mV over 12 d.

Osteocalcin binds to a GPRC6A Venus fly trap allosteric site to positively modulate GPRC6A signaling.

GPRC6A, a member of the Family C G-protein coupled receptors, regulates energy metabolism and sex hormone production and is activated by diverse ligands, including cations, L-amino acids, the osteocalcin (Ocn) peptide and the steroid hormone testosterone. We sought a structural framework for the ability of multiple distinct classes of ligands to active GPRC6A. We created a structural model of GPRC6A using Alphafold2. Using this model we explored a putative orthosteric ligand binding site in the bilobed Venus fly trap (VFT) domain of GPRC6A and two positive allosteric modulator (PAM) sites, one in the VFT and the other in the 7 transmembrane (7TM) domain. We provide evidence that Ocn peptides act as a PAM for GPRC6A by binding to a site in the VFT that is distinct from the orthosteric site for calcium and L-amino acids. In agreement with this prediction, alternatively spliced GPRC6A isoforms 2 and 3, which lack regions of the VFT, and mutations in the computationally predicted Ocn binding site, K352E and H355P, prevent Ocn activation of GPRC6A. These observations explain how dissimilar ligands activate GPRC6A and set the stage to develop novel molecules to activate and inhibit this previously poorly understood receptor.

Structural insight into binding site access and ligand recognition by human ABCB1.

ABCB1 is a broad-spectrum efflux pump central to cellular drug handling and multidrug resistance in humans. However, its mechanisms of poly-specific substrate recognition and transport remain poorly resolved. Here we present cryo-EM structures of lipid embedded human ABCB1 in its apo, substrate-bound, inhibitor-bound, and nucleotide-trapped states at 3.4-3.9 Å resolution without using stabilizing antibodies or mutations and each revealing a distinct conformation. The substrate binding site is located within one half of the molecule and, in the apo state, is obstructed by transmembrane helix (TM) 4. Substrate and inhibitor binding are distinguished by major differences in TM arrangement and ligand binding chemistry, with TM4 playing a central role in all conformational transitions. Our data offer fundamental new insights into the role structural asymmetry, secondary structure breaks, and lipid interactions play in ABCB1 function and have far-reaching implications for ABCB1 inhibitor design and predicting its substrate binding profiles.