Localization Of Binding Sites Research Articles

Research Progress on Spike-Dependent SARS-CoV-2 Fusion Inhibitors and Small Molecules Targeting the S2 Subunit of Spike

Since the beginning of the COVID-19 pandemic, extensive drug repurposing efforts have sought to identify small-molecule antivirals with various mechanisms of action. Here, we aim to review research progress on small-molecule viral entry and fusion inhibitors that directly bind to the SARS-CoV-2 Spike protein. Early in the pandemic, numerous small molecules were identified in drug repurposing screens and reported to be effective in in vitro SARS-CoV-2 viral entry or fusion inhibitors. However, given minimal experimental information regarding the exact location of small-molecule binding sites on Spike, it was unclear what the specific mechanism of action was or where the exact binding sites were on Spike for some inhibitor candidates. The work of countless researchers has yielded great progress, with the identification of many viral entry inhibitors that target elements on the S1 receptor-binding domain (RBD) or N-terminal domain (NTD) and disrupt the S1 receptor-binding function. In this review, we will also focus on highlighting fusion inhibitors that target inhibition of the S2 fusion function, either by disrupting the formation of the postfusion S2 conformation or alternatively by stabilizing structural elements of the prefusion S2 conformation to prevent conformational changes associated with S2 function. We highlight experimentally validated binding sites on the S1/S2 interface and on the S2 subunit. While most substitutions to the Spike protein to date in variants of concern (VOCs) have been localized to the S1 subunit, the S2 subunit sequence is more conserved, with only a few observed substitutions in proximity to S2 binding sites. Several recent small molecules targeting S2 have been shown to have robust activity over recent VOC mutant strains and/or greater broad-spectrum antiviral activity for other more distantly related coronaviruses.

LogoMotif: A Comprehensive Database of Transcription Factor Binding Site Profiles in Actinobacteria

Actinobacteria undergo a complex multicellular life cycle and produce a wide range of specialized metabolites, including the majority of the antibiotics. These biological processes are controlled by intricate regulatory pathways, and to better understand how they are controlled we need to augment our insights into the transcription factor binding sites. Here, we present LogoMotif (https://logomotif.bioinformatics.nl), an open-source database for characterized and predicted transcription factor binding sites in Actinobacteria, along with their cognate position weight matrices and hidden Markov models. Genome-wide predictions of binding site locations in Streptomyces model organisms are supplied and visualized in interactive regulatory networks. In the web interface, users can freely access, download and investigate the underlying data. With this curated collection of actinobacterial regulatory interactions, LogoMotif serves as a basis for binding site predictions, thus providing users with clues on how to elicit the expression of genes of interest and guide genome mining efforts.

Control of polymer-protein interactions by tuning the composition and length of polymer chains.

Nanomoduling the 3D shape and chemical functionalities in a synthetic polymer may create recognition cavities for biomacromolecule binding, which serves as an attractive alternative to natural antibodies with much less cost. To obtain fundamental understanding and predict molecular design rules of the polymer antibody, we analyze the complex structure between the biomarker protein epithelial cell adhesion molecule (EpCAM) and a series of polymer ligands via molecular dynamics (MD) simulations. For monomeric ligands, strong enrichment of aromatic residues in protein binding sites is revealed, in line with the reported observations for natural antibodies. Yet, for linear polymers with a growing degree of polymerization, for the first time, a drastic change is revealed on the type of enriched protein residues and the location of protein binding sites, driven by the increasing steric hindrance effect that makes the adsorption of the polymer in the protein exterior feasible. Varying the polymer length and monomeric composition also significantly affects the ligand binding affinity. Here, we have captured three distinct dependences of the ligand binding free energy on the degree of polymerization: for NIPAm based hydrophilic polymers, TBAm dominated hydrophobic polymers and AAc dominated charged polymers. These results can be rationalized by the complex structure and the composition of protein residues at the binding interface. The entire analysis demonstrates unique binding features for polymer ligands and the possibility to modulate their binding sites and affinity by engineering the polymer structure.

Elucidation of the structural basis for ligand binding and translocation in conserved insect odorant receptor co-receptors

In numerous insects, the olfactory receptor family forms a unique class of heteromeric cation channels. Recent progress in resolving the odorant receptor structures offers unprecedented opportunities for deciphering their molecular mechanisms of ligand recognition. Unexpectedly, these structures in apo or ligand-bound states did not reveal the pathway taken by the ligands between the extracellular space and the deep internal cavities. By combining molecular modeling with electrophysiological recordings, we identified amino acids involved in the dynamic entry pathway and the binding of VUAA1 to Drosophila melanogaster’s odorant receptor co-receptor (Orco). Our results provide evidence for the exact location of the agonist binding site and a detailed and original mechanism of ligand translocation controlled by a network of conserved residues. These findings would explain the particularly high selectivity of Orcos for their ligands.

Structures of the DarR transcription regulator reveal unique modes of second messenger and DNA binding

The mycobacterial repressor, DarR, a TetR family regulator (TFR), was the first transcription regulator shown to bind c-di-AMP. However, the molecular basis for this interaction and the mechanism involved in DNA binding by DarR remain unknown. Here we describe DarR-c-di-AMP and DarR-DNA structures and complementary biochemical assays. The DarR-c-di-AMP structure reveals a unique effector binding site for a TFR, located between DarR dimer subunits. Strikingly, we show this motif also binds cAMP. The location of the adenine nucleotide binding site between subunits suggests this interaction may facilitate dimerization and hence DNA binding. Indeed, biochemical assays show cAMP enhances DarR DNA binding. Finally, DarR-DNA structures reveal a distinct TFR DNA-binding mechanism involving two interacting dimers on the DNA. Thus, the combined data unveil a newly described second messenger binding motif and DNA binding mode for this important family of regulators.

Insight into the mechanism of H+-coupled nucleobase transport

Members of the nucleobase/ascorbic acid transporter (NAT) gene family are found in all kingdoms of life. In mammals, the concentrative uptake of ascorbic acid (vitamin C) by members of the NAT family is driven by the Na+ gradient, while the uptake of nucleobases in bacteria is powered by the H+ gradient. Here, we report the structure and function of PurTCp, a NAT family member from Colwellia psychrerythraea. The structure of PurTCp was determined to 2.80 Å resolution by X-ray crystallography. PurTCp forms a homodimer, and each protomer has 14 transmembrane segments folded into a transport domain (core domain) and a scaffold domain (gate domain). A purine base is present in the structure and defines the location of the substrate binding site. Functional studies reveal that PurTCp transports purines but not pyrimidines and that purine binding and transport is dependent on the pH. Mutation of a conserved aspartate residue close to the substrate binding site reveals the critical role of this residue in H+-dependent transport of purines. Comparison of the PurTCp structure with transporters of the same structural fold suggests that rigid-body motions of the substrate-binding domain are central for substrate translocation across the membrane.

Development of Riluzole Analogs with Improved Use-Dependent Inhibition of Skeletal Muscle Sodium Channels.

Several commercially available and newly synthesized riluzole analogs were evaluated in vitro as voltage-gated skeletal muscle sodium-channel blockers. Data obtained from the patch-clamp technique demonstrated that potency is well correlated with lipophilicity and the introduction of a protonatable amino function in the benzothiazole 2-position enhances the use-dependent behavior. The most interesting compound, the 2-piperazine analog of riluzole (14), although slightly less potent than the parent compound in the patch-clamp assay as well as in an in vitro model of myotonia, showed greater use-dependent Nav1.4 blocking activity. Docking studies allowed the identification of the key interactions that 14 makes with the amino acids of the local anesthetic binding site within the pore of the channel. The reported results pave the way for the identification of novel compounds useful in the treatment of cell excitability disorders.

Identification of a Cryptic Binding Site in CRISPR-Cas9 for Targeted Inhibition.

The need for precision control of CRISPR-Cas9 genome editing has created a demand for anti-CRISPR molecules. Recently, the first class of small-molecule Cas9 inhibitors has been identified, verifying the feasibility of regulating CRISPR-Cas9 activity using direct-acting small molecules. However, it remains enigmatic as to the location of the ligand binding site(s) on CRISPR-Cas9 and how the ligand binding leads to Cas9 functional inhibition. Here, we established an integrative computational protocol, including massive binding site mapping, molecular docking, molecular dynamics simulations, and free energy calculations. Ultimately, a Cas9 ligand binding site was discovered from the dynamics trajectories that is hidden within its carboxyl-terminal domain (CTD), a domain recognizing the protospacer adjacent motif (PAM). Using the top inhibitor BRD0539 as a probe, we demonstrated that the ligand binding induces significant CTD structural rearrangements toward an incompetent conformation for PAM DNA engagement. The revealed molecular mechanism of BRD0539 inhibiting Cas9 is in well agreement with the experimental data. This study provides a structural and mechanistic basis for the potency improvement of existing ligands and the rational discovery of novel small-molecule brakes for developing safer CRISPR-Cas9 technologies.

AttentionDTA: Drug-Target Binding Affinity Prediction by Sequence-Based Deep Learning With Attention Mechanism.

The identification of drug-target relations (DTRs) is substantial in drug development. A large number of methods treat DTRs as drug-target interactions (DTIs), a binary classification problem. The main drawback of these methods are the lack of reliable negative samples and the absence of many important aspects of DTR, including their dose dependence and quantitative affinities. With increasing number of publications of drug-protein binding affinity data recently, DTRs prediction can be viewed as a regression problem of drug-target affinities (DTAs) which reflects how tightly the drug binds to the target and can present more detailed and specific information than DTIs. The growth of affinity data enables the use of deep learning architectures, which have been shown to be among the state-of-the-art methods in binding affinity prediction. Although relatively effective, due to the black-box nature of deep learning, these models are less biologically interpretable. In this study, we proposed a deep learning-based model, named AttentionDTA, which uses attention mechanism to predict DTAs. Different from the models using 3D structures of drug-target complexes or graph representation of drugs and proteins, the novelty of our work is to use attention mechanism to focus on key subsequences which are important in drug and protein sequences when predicting its affinity. We use two separate one-dimensional Convolution Neural Networks (1D-CNNs) to extract the semantic information of drug's SMILES string and protein's amino acid sequence. Furthermore, a two-side multi-head attention mechanism is developed and embedded to our model to explore the relationship between drug features and protein features. We evaluate our model on three established DTA benchmark datasets, Davis, Metz, and KIBA. AttentionDTA outperforms the state-of-the-art deep learning methods under different evaluation metrics. The results show that the attention-based model can effectively extract protein features related to drug information and drug features related to protein information to better predict drug target affinities. It is worth mentioning that we test our model on IC50 dataset, which provides the binding sites between drugs and proteins, to evaluate the ability of our model to locate binding sites. Finally, we visualize the attention weight to demonstrate the biological significance of the model. The source code of AttentionDTA can be downloaded from https://github.com/zhaoqichang/AttentionDTA_TCBB.

Structural insights for the modulation of KCNQ1 channel by ML277.

KCNQ1 is a voltage-gated potassium ion channel protein which co-assembles with its accessory regulatory KCNE subunits, KCNE1-5, in order to fulfill its diverse physiological roles in various tissues. In the heart, KCNQ1 helps maintain a normal QT interval in the electrocardiogram, and loss of function mutations can lead to long QT syndrome type 1 (LQT1), seizures, and sudden death. Activators of KCNQ1 may be a promising strategy to oppose the effect of LQT mutations. In this work, we determined the cryo-EM structure of Xenopus KCNQ1 bound to Ca2+/calmodulin with or without the KCNQ1 activator, ML277. Our results showed that a binding site for ML277 lies in the transmembrane domain of the channel which is lined by residues of the S4-S5 linker, and the S5 and S6 helices of two adjacent subunits. MD simulations, along with point mutations, support the binding site location and showed that ML277 binds almost identically in open and closed channels. Comparative structural analyses highlighted important residues for the selectivity of ML277 towards KCNQ1, among other KCNQ subfamily members, and provided an explanation for the limited efficacy of ML277 in the presence of KCNE3. This work paves the way for a better understanding of KCNQ1 regulation by small molecules and may lead to the development of novel and potent activators of KCNQ1, providing a framework for therapeutic intervention targeting KCNQ1 loss of function pathologies.

Exploring molecular interaction of cefpirome with human serum albumin: In vitro and in silico approaches

Cefpirome (CFP), a fourth-generation cephalosporin for parenteral administration, is highly influential against gram-negative and gram-positive bacteria. Multiple approaches, such as fluorescence, absorption, voltammetric, and molecular docking were applied to assess the binding of CFP to the primary carrier protein, human serum albumin (HSA). The CFP-triggered quenching of HSA fluorescence reported the CFP−HSA complex formation, which was accomplished by a static quenching mechanism. A moderate binding affinity (Kf = 1.78 × 104 M−1 at 298 K) was found in the CFP−HSA interaction. The complex was anticipated to be stabilized by the presence of van der Waals interactions and hydrogen bonds. Upon CFP binding, microenvironmental alterations surrounding Trp and Tyr residues of HSA were observed. To establish the location of the CFP binding site in HSA, molecular docking studies were conducted, and the results indicated that CFP attaches to site I of HSA mostly through hydrophobic interactions and hydrogen bonding. The results of fluorescence spectroscopy and molecular docking supported one another quite well.

Glycosaminoglycan Analysis: Purification, Structural Profiling, and GAG-Protein Interactions.

Glycosaminoglycans (GAGs) are long, linear polysaccharides that are ubiquitously expressed on the cell surface and in the extracellular matrix of all animal cells. These complex carbohydrates are composed of alternating glucosamine and uronic acids that can be heterogeneously N- and O-sulfated. The arrangement and orientation of the sulfated sugar residues specify the location of distinct ligand binding sites on the cell surface, and their capacity to bind ligands impacts cell growth and development, the ability to form tissues and organs, and normal physiology. The heterogeneous nature of GAGs and their inherent structural diversity across different tissues, cell types, and disease states creates challenges to characterizing their structure and function. Here, we describe detailed methods to investigate GAG-protein interactions in vitro and evaluate the structural composition of two classes of sulfated GAGs, heparan sulfate and chondroitin/dermatan sulfate, using liquid chromatography, mass spectrometry, and radiolabeling techniques. Overall, these methods facilitate the evaluation of GAG structure and function to uncover the unique roles these molecules play in cell biology and human disease.

Local Anesthetic Like Inhibition of the Cardiac Na+ Channel Nav1.5 by Chloroquine and Hydroxychloroquine.

Chloroquine (CQ) and its derivate hydroxychloroquine (HCQ) are successfully deployed for different diseases beyond the prophylaxis and treatment of malaria. Both substances exhibit antiviral properties and have been proposed for prophylaxis and treatment of COVID-19 caused by SARS-CoV-2. CQ and HCQ cause similar adverse events including life-threatening cardiac arrhythmia generally based on QT-prolongation, which is one of the most reported adverse events for both agents associated with the treatment of COVID-19. Various drugs known to induce QT-prolongation have been proven to exert local anesthetic (LA)-like properties regarding their impact on the cardiac Na+ channel Nav1.5. Inhibition of Nav1.5 is considered as the primary mechanism of cardiotoxicity caused by LAs. However, the mechanism of the arrhythmogenic effects of CQ and HCQ related to Nav1.5 has not yet been fully investigated. Therefore, the exact mechanism of how CQ and HCQ affect the sodium currents generated by Nav1.5 need to be further elucidated. This in vitro study aims to investigate the effects of CQ and HCQ on Nav1.5-generated sodium currents to identify possible LA-like mechanisms that might contribute to their arrhythmogenic properties. The effects of CQ and HCQ on Nav1.5-generated sodium currents by HEK-293 cells expressing either wild-type human Nav1.5 or mutant Nav1.5 F1760A are measured using the whole-cell patch-clamp technique. Both agents induce a state-dependent inhibition of Nav1.5. Furthermore, CQ and HCQ produce a use-dependent block of Nav1.5 and a shift of fast and slow inactivation. Results of experiments investigating the effect on the LA-insensitive mutant Nav1.5-F1760A indicate that both agents at least in part employ the proposed LA-binding site of Nav1.5 to induce inhibition. This study demonstrated that CQ and HCQ exert LA-typical effects on Nav1.5 involving the proposed LA binding site, thus contributing to their arrhythmogenic properties.

Homology Modeling Odorant-binding Protein-1 (OBP1) Anopheles Farauti Protein Target for Mosquito Repellent

The odorant binding protein (OBP) has a role as a target protein for potential interaction mechanism activity for the development of repellent compounds. The purpose of this study was to analyze the physico-chemical properties of the protein, the stability of the three-dimensional structure of the OBP1 Anopheles farauti protein, and to predict the binding site pocket as the target of the active protein site against inhibitors. Analysis of physico-chemical properties was carried out by the ProtParam Expasy server. The theoretical calculated isoelectric point (pI) was found to be less than 7 indicating the acidic nature of this protein. The aliphatic index of 78 indicates the thermal stability of the protein. The Grand Average of Hydropathicity (GRAVY) is estimated at -0.355; This lower GRAVY value indicates a possible better interaction of this protein with water. Secondary structure analysis was carried out by SOPMA which revealed that Alpha helix (55.86%) predominated among the secondary structural elements followed by Random coil (32.41%), Extended strand (8.97%), and Beta turn (2.76%). Three-dimensional structure modeling of OBP1 Anopheles farauti was performed with the Swiss-Model server and the protein refine Galaxy server. Homology modeling results obtained PDB ID 2ERB template with sequence identity 94.4%. The model was validated for the three-dimensional structure of the protein using the MolProbity, ProSA, ProQ, ERRAT, Verfy3D, and PROCHECK servers. The prediction results of pocket binding sites using DoSiteScore obtained three pocket binding site locations, namely P_0 (Drug score 0.84); P_1 ((0.75); P_2 (Drug score 0.64). Conclusion Homology modeling of the protein OBP1 Anopheles farauti has been carried out and the three-dimensional structure of the Model_OBP1_04 protein has been obtained that meets the criteria for valid structural parameters so that this structure can be used In-silico molecular docking and molecular dynamic studies for the development of mosquitoes repellent.

Structure-Based Discovery of Negative Allosteric Modulators of the Metabotropic Glutamate Receptor 5.

Recently determined structures of class C G protein-coupled receptors (GPCRs) revealed the location of allosteric binding sites and opened new opportunities for the discovery of novel modulators. In this work, molecular docking screens for allosteric modulators targeting the metabotropic glutamate receptor 5 (mGlu5) were performed. The mGlu5 receptor is activated by the main excitatory neurotransmitter of the nervous central system, L-glutamate, and mGlu5 receptor activity can be allosterically modulated by negative or positive allosteric modulators. The mGlu5 receptor is a promising target for the treatment of psychiatric and neurodegenerative diseases, and several allosteric modulators of this GPCR have been evaluated in clinical trials. Chemical libraries containing fragment- (1.6 million molecules) and lead-like (4.6 million molecules) compounds were docked to an allosteric binding site of mGlu5 identified in X-ray crystal structures. Among the top-ranked compounds, 59 fragments and 59 lead-like compounds were selected for experimental evaluation. Of these, four fragment- and seven lead-like compounds were confirmed to bind to the allosteric site with affinities ranging from 0.43 to 8.6 μM, corresponding to a hit rate of 9%. The four compounds with the highest affinities were demonstrated to be negative allosteric modulators of mGlu5 signaling in functional assays. The results demonstrate that virtual screens of fragment- and lead-like chemical libraries have complementary advantages and illustrate how access to high-resolution structures of GPCRs in complex with allosteric modulators can accelerate lead discovery.

Spectrofluorometric and computational approaches for the interaction studies of aclonifen and bifenox with human serum albumin

Interaction of two broadly used herbicides, aclonifen (ACF) and bifenox (BIF) with the major transporter in human circulation, human serum albumin (HSA) were examined using fluorescence and absorption spectral measurements combined with in silico analyses. Assessment of the fluorescence and absorption spectral results affirmed the complexation between ACF/BIF and HSA. Increase in the KSV value with temperature characterized the ACF/BIF-induced quenching of the protein fluorescence as dynamic quenching. The moderate binding affinities (Kf = 1.74×104 – 1.95×106 M−1 for ACF–HSA complex; Kf = 2.00×103 – 1.02×106 M−1 for BIF–HSA complex) were pointed out between ACF/BIF and HSA, showing a relatively higher binding constant values with increasing temperatures. Quantitative evaluation of thermodynamic data (ΔS = +0.86 kJ mol−1 K−1 and ΔH = +225.43 kJ mol−1 for ACF–HSA complex; ΔS = +1.11 kJ mol−1 K−1 and ΔH = +304.63 kJ mol−1 for BIF–HSA complex) predicted the contribution of hydrophobic interactions in the ACF–HSA and BIF–HSA association processes, which were well supported by our molecular docking results. In silico analyses were made to acquire insight details into the ACF and BIF binding to HSA at the binding sites and suggested the locations of ACF and BIF binding sites as both subdomain IIA (site I) and subdomain IIIA (site II) of HSA, showing more preference toward site I.

CavitySpace: A Database of Potential Ligand Binding Sites in the Human Proteome.

Location and properties of ligand binding sites provide important information to uncover protein functions and to direct structure-based drug design approaches. However, as binding site detection depends on the three-dimensional (3D) structural data of proteins, functional analysis based on protein ligand binding sites is formidable for proteins without structural information. Recent developments in protein structure prediction and the 3D structures built by AlphaFold provide an unprecedented opportunity for analyzing ligand binding sites in human proteins. Here, we constructed the CavitySpace database, the first pocket library for all the proteins in the human proteome, using a widely-applied ligand binding site detection program CAVITY. Our analysis showed that known ligand binding sites could be well recovered. We grouped the predicted binding sites according to their similarity which can be used in protein function prediction and drug repurposing studies. Novel binding sites in highly reliable predicted structure regions provide new opportunities for drug discovery. Our CavitySpace is freely available and provides a valuable tool for drug discovery and protein function studies.

Low potency inhibition of NaV1.7 by externally applied QX-314 via a depolarizing shift in the voltage-dependence of activation

QX-314 is a quaternary permanently charged lidocaine derivative that inhibits voltage-gated sodium channels (NaV). As it is membrane impermeable, it is generally considered that QX-314 applied externally is inactive, unless it can gain access to the internal local anesthetic binding site via another entry pathway. Here, we characterized the electrophysiological effects of QX-314 on NaV1.7 heterologously expressed in HEK293 cells, and found that at high concentrations, external QX-314 inhibited NaV1.7 current (IC50 2.0 ± 0.3 mM) and shifted the voltage-dependence to more depolarized potentials (ΔV50 +10.6 mV). Unlike lidocaine, the activity of external QX-314 was not state- or use-dependent. The effect of externally applied QX-314 on NaV1.7 channel biophysics differed to that of internally applied QX-314, suggesting QX-314 has an additional externally accessible site of action. In line with this hypothesis, disruption of the local anesthetic binding site in a [F1748A]NaV1.7 mutant reduced the potency of lidocaine by 40-fold, but had no effect on the potency or activity of externally applied QX-314. Therefore, we conclude, using an expression system where QX-314 was unable to cross the membrane, that externally applied QX-314 is able to inhibit NaV1.7 peak current at low millimolar concentrations.

Computational Modelling of the Interactions Between Polyoxometalates and Biological Systems.

Polyoxometalates (POMs) structures have raised considerable interest for the last years in their application to biological processes and medicine. Within this area, our mini-review shows that computational modelling is an emerging tool, which can play an important role in understanding the interaction of POMs with biological systems and the mechanisms responsible of their activity, otherwise difficult to achieve experimentally. During recent years, computational studies have mainly focused on the analysis of POM binding to proteins and other systems such as lipid bilayers and nucleic acids, and on the characterization of reaction mechanisms of POMs acting as artificial metalloproteases and phosphoesterases. From early docking studies locating binding sites, molecular dynamics (MD) simulations have allowed to characterize the nature of POM···protein interactions, and to evaluate the effect of the charge, size, and shape of the POM on protein affinity, including also, the atomistic description of chaotropic character of POM anions. Although these studies rely on the interaction with proteins and nucleic acid models, the results could be extrapolated to other biomolecules such as carbohydrates, triglycerides, steroids, terpenes, etc. Combining MD simulations with quantum mechanics/molecular mechanics (QM/MM) methods and DFT calculations on cluster models, computational studies are starting to shed light on the factors governing the activity and selectivity for the hydrolysis of peptide and phosphoester bonds catalysed by POMs.

Structural insights into recognition of chemokine receptors by Staphylococcus aureus leukotoxins.

Staphylococcus aureus (SA) leukocidin ED (LukED) belongs to a family of bicomponent pore forming toxins that play important roles in SA immune evasion and nutrient acquisition. LukED targets specific G protein-coupled chemokine receptors to lyse human erythrocytes (red blood cells) and leukocytes (white blood cells). The first recognition step of receptors is critical for specific cell targeting and lysis. The structural and molecular bases for this mechanism are not well understood but could constitute essential information to guide antibiotic development. Here, we characterized the interaction of LukE with chemokine receptors ACKR1, CCR2, and CCR5 using a combination of structural, pharmacological, and computational approaches. First, crystal structures of LukE in complex with a small molecule mimicking sulfotyrosine side chain (p-cresyl sulfate) and with peptides containing sulfotyrosines issued from receptor sequences revealed the location of receptor sulfotyrosine binding sites in the toxins. Then, by combining previous and novel experimental data with protein docking, classical and accelerated weight histogram (AWH) molecular dynamics we propose models of the ACKR1-LukE and CCR5-LukE complexes. This work provides novel insights into chemokine receptor recognition by leukotoxins and suggests that the conserved sulfotyrosine binding pocket could be a target of choice for future drug development.