Protein-protein Interaction Interface Research Articles

Inhibiting Type VI Secretion System Activity with a Biomimetic Peptide Designed To Target the Baseplate Wedge Complex.

ABSTRACTHuman health is threatened by bacterial infections that are increasingly resistant to multiple drugs. A recently emerged strategy consists of disarming pathogenic bacteria by targeting and blocking their virulence factors. The type VI secretion system (T6SS) is a widespread secretion nanomachine encoded and employed by pathogenic strains to establish their virulence process during host invasion. Given the conservation of T6SS in several human bacterial pathogens, the discovery of an effective broad-spectrum T6SS virulence blocker represents an attractive target for development of antivirulence therapies. Here, we identified and validated a protein-protein interaction interface, TssK-TssG, as a key factor in the assembly of the T6SS baseplate (BP) complex in the pathogen enteroaggregative Escherichia coli (EAEC). In silico and biochemical studies revealed that the determinants of the interface are broadly conserved among pathogenic species, suggesting a role for this interface as a target for T6SS inhibition. Based on the high-resolution structure of the TssKFGE wedge complex, we rationally designed a biomimetic cyclic peptide (BCP) that blocks the assembly of the EAEC BP complex and inhibits the function of T6SS in bacterial cultures. Our BCP is the first compound completely designed from prior structural knowledge with anti-T6SS activity that can be used as a model to target human pathogens.

Reversible Covalent Imine-Tethering for Selective Stabilization of 14-3-3 Hub Protein Interactions.

The stabilization of protein complexes has emerged as a promising modality, expanding the number of entry points for novel therapeutic intervention. Targeting proteins that mediate protein–protein interactions (PPIs), such as hub proteins, is equally challenging and rewarding as they offer an intervention platform for a variety of diseases, due to their large interactome. 14-3-3 hub proteins bind phosphorylated motifs of their interaction partners in a conserved binding channel. The 14-3-3 PPI interface is consequently only diversified by its different interaction partners. Therefore, it is essential to consider, additionally to the potency, also the selectivity of stabilizer molecules. Targeting a lysine residue at the interface of the composite 14-3-3 complex, which can be targeted explicitly via aldimine-forming fragments, we studied the de novo design of PPI stabilizers under consideration of potential selectivity. By applying cooperativity analysis of ternary complex formation, we developed a reversible covalent molecular glue for the 14-3-3/Pin1 interaction. This small fragment led to a more than 250-fold stabilization of the 14-3-3/Pin1 interaction by selective interfacing with a unique tryptophan in Pin1. This study illustrates how cooperative complex formation drives selective PPI stabilization. Further, it highlights how specific interactions within a hub proteins interactome can be stabilized over other interactions with a common binding motif.

Elucidating the Structural Constituents of Fibroblast Growth Factor 14 that Confer its Modulatory Effects on the Voltage-Gated Sodium Channel 1.6

Disruption of protein:protein interactions (PPIs) between voltage-gated Na+ (Nav) channels and their auxiliary proteins is implicated in the etiology of channelopathies. Despite restoration of these perturbed PPIs serving as a novel therapeutic approach, such efforts to develop PPI-based pharmacotherapies are hindered by lack of insights into the structural constituents of regulatory proteins that confer their modulatory effects on Nav channel kinetics. Focusing on the PPI between Nav1.6 and its auxiliary protein fibroblast growth factor 14 (FGF14), we previously identified clusters of amino acids on the β12 sheet and β8-β9 loop of FGF14 that are at its PPI interface with the C-terminal domain (CTD) of Nav1.6 and are predicted to disproportionately contribute to the Gibbs energy of their protein:protein binding. To develop chemical probes, we derived peptidomimetics from the PLEV and EYYV motifs of the β12 sheet and β8-β9 loop, respectively. Crucially, a FGF14PLEV derivative caused a depolarizing shift in the V1/2 of Nav1.6 channel activation and steady-state inactivation, but did not affect long-term inactivation. Conversely, a FGF14EYYV derivate caused an increase in the long-term inactivation of Nav1.6 channels, but did not affect the V1/2 of activation or steady-state inactivation, suggesting that the β12 sheet and β8-β9 loop constituents of FGF14 regulate divergent kinetic properties of Nav1.6 channels. Based upon molecular docking of these two derivatives, it is hypothesized that their divergent regulation of Nav1.6 channel kinetics is conferred by differential interactions with the EF hand like (EFL) and IQ domains of the CTD. Exploration of this hypothesis will enable an increasingly refined understanding of Nav channel regulation and provide insights into developing therapeutics targeted against PPI interfaces between Nav channels and their auxiliary proteins.

Fluorescence Anisotropy-Based Tethering for Discovery of Protein-Protein Interaction Stabilizers.

Protein–protein interaction (PPI) networks are fundamental for cellular processes. Small-molecule PPI enhancers have been shown to be powerful tools to fundamentally study PPIs and as starting points for potential new therapeutics. Yet, systematic approaches for their discovery are not widely available, and the design prerequisites of “molecular glues” are poorly understood. Covalent fragment-based screening can identify chemical starting points for these enhancers at specific sites in PPI interfaces. We recently reported a mass spectrometry-based disulfide-trapping (tethering) approach for a cysteine residue in the hub protein 14–3–3, an important regulator of phosphorylated client proteins. Here, we invert the strategy and report the development of a functional read-out for systematic identification of PPI enhancers based on fluorescence anisotropy (FA-tethering) with the reactive handle now on a client-derived peptide. Using the DNA-binding domain of the nuclear receptor Estrogen Related Receptor gamma (ERRγ), we target a native cysteine positioned at the 14–3–3 PPI interface and identify several fragments that form a disulfide bond to ERRγ and stabilize the complex up to 5-fold. Crystallography indicates that fragments bind in a pocket comprised of 14–3–3 and the ERRγ phosphopeptide. FA-tethering presents a streamlined methodology to discover molecular glues for protein complexes.

Experimental and Computational Druggability Exploration of the 14-3-3ζ/SOS1pS1161 PPI Interface.

The exploration of the druggability of certain protein-protein interactions (PPIs) still remains a challenging task in drug discovery. Here, we present a case study using the 14-3-3-PPI, showing how small molecules can be located that are able to modulate this key oncogenic pathway. A workflow embracing biophysical techniques and MD simulations was developed to evaluate the potential of a 14-3-3ζ PPI system to bind new tool compounds. The significance of the use of computational approaches to compensate for the limitations of experimental techniques is demonstrated.

In silico identification of small molecule modulators for disruption of Hsp90–Cdc37 protein–protein interaction interface for cancer therapeutic application

The protein–protein interactions (PPIs) in the biological systems are important to maintain a number of cellular processes. Several disorders including cancer may be developed due to dysfunction in the assembly of PPI networks. Hence, targeting intracellular PPIs can be considered as a crucial drug target for cancer therapy. Among the enormous and diverse group of cancer-enabling PPIs, the Hsp90–Cdc37 is prominent for cancer therapeutic development. The successful inhibition of Hsp90–Cdc37 PPI interface can be an important therapeutic option for cancer management. In the current study, a set of more than sixty thousand compounds belong to four databases were screened through a multi-steps molecular docking study in Glide against the Hsp90–Cdc37 interaction interface. The Glide-score and Prime-MM-GBSA based binding free energy of DCZ3112, standard Hsp90–Cdc37 inhibitor were found to be −6.96 and −40.46 kcal/mol, respectively. The above two parameters were used as cut-off score to reduce the chemical space from all successfully docked molecules. Furthermore, the in-silico pharmacokinetics parameters, common-feature pharmacophore analyses and the molecular binding interactions were used to wipe out the inactive molecules. Finally, four molecules were found to be important to modulate the Hsp90–Cdc37 interface. The potentiality of the final four molecules was checked through several drug-likeness characteristics. The molecular dynamics (MD) simulation study explained that all four molecules retained inside the interface of Hsp90–Cdc37. The binding free energy of each molecule obtained from the MD simulation trajectory was clearly explained the strong affection towards the Hsp90–Cdc37. Hence, the proposed molecule might be crucial for successful inhibition of the Hsp90–Cdc37 interface. Communicated by Ramaswamy H. Sarma

Design, Synthesis, and Pharmacological Evaluation of Analogues Derived from the PLEV Tetrapeptide as Protein-Protein Interaction Modulators of Voltage-Gated Sodium Channel 1.6.

The voltage-gated Na+ (Nav) channel is the molecular determinant of excitability. Disruption of protein-protein interactions (PPIs) between Nav1.6 and fibroblast growth factor 14 (FGF14) leads to impaired excitability of neurons in clinically relevant brain areas associated with channelopathies. Here, we designed, synthesized, and pharmacologically characterized new peptidomimetics based on a PLEV tetrapeptide scaffold derived from the FGF14:Nav1.6 PPI interface. Addition of an N-terminal 1-adamantanecarbonyl pharmacophore significantly improved peptidomimetic inhibitory potency. Surface plasmon resonance studies revealed that while this moiety was sufficient to confer binding to FGF14, altering the C-terminal moiety from methoxy (21a) to π bond-containing (23a and 23b) or cycloalkane substituents (23e) abrogated the binding to Nav1.6. Whole-cell patch-clamp electrophysiology subsequently revealed that 21a had functionally relevant interactions with both the C-terminal tail of Nav1.6 and FGF14. Collectively, these findings support that 21a (PW0564) may serve as a promising lead to develop target-selective neurotherapeutics by modulating protein-channel interactions.

Fragment-Based Stabilizers of Protein-Protein Interactions through Imine-Based Tethering.

Small‐molecule stabilization of protein–protein interactions (PPIs) is a promising concept in drug discovery, however the question how to identify or design chemical starting points in a “bottom‐up” approach is largely unanswered. We report a novel concept for identifying initial chemical matter for PPI stabilization based on imine‐forming fragments. The imine bond offers a covalent anchor for site‐directed fragment targeting, whereas its transient nature enables efficient analysis of structure–activity relationships. This bond enables fragment identification and optimisation using protein crystallography. We report novel fragments that bind specifically to a lysine at the PPI interface of the p65‐subunit‐derived peptide of NF‐κB with the adapter protein 14‐3‐3. Those fragments that subsequently establish contacts with the p65‐derived peptide, rather than with 14‐3‐3, efficiently stabilize the 14‐3‐3/p65 complex and offer novel starting points for molecular glues.

Fragment‐Based Stabilizers of Protein–Protein Interactions through Imine‐Based Tethering

AbstractSmall‐molecule stabilization of protein–protein interactions (PPIs) is a promising concept in drug discovery, however the question how to identify or design chemical starting points in a “bottom‐up” approach is largely unanswered. We report a novel concept for identifying initial chemical matter for PPI stabilization based on imine‐forming fragments. The imine bond offers a covalent anchor for site‐directed fragment targeting, whereas its transient nature enables efficient analysis of structure–activity relationships. This bond enables fragment identification and optimisation using protein crystallography. We report novel fragments that bind specifically to a lysine at the PPI interface of the p65‐subunit‐derived peptide of NF‐κB with the adapter protein 14‐3‐3. Those fragments that subsequently establish contacts with the p65‐derived peptide, rather than with 14‐3‐3, efficiently stabilize the 14‐3‐3/p65 complex and offer novel starting points for molecular glues.

Exploring CIP2A modulators using multiple molecular modeling approaches

Like other human oncoproteins, Cancerous Inhibitor of Protein Phosphatase 2A (CIP2A) exerts cancer promoting function through interaction with other partner proteins, such as MYC and Protein Phosphatase 2A (PP2A). CIP2A regulates several MYC-independent and/or dependent gene expression programs. Broadly, CIP2A can inhibit PP2A, and especially it has been shown to inhibit MYC-associated PP2A, precisely to increase MYC stability and function. Availability of crystal structure has broached the research focus to develop new therapeutics targeting CIP2A. In the present study, structural information of the protein has been used for identification of modulators for homo-dimer CIP2A using advanced structure-based drug design approaches. The compound library, ‘Maybridge Screening Collection’ database (∼62,000 compounds) has been virtually screened to find out potent modulators for CIP2A. Identification of hotspot region on CIP2A protein-protein interaction interface has been performed using three different tools (HotPoint, SiteMap and icmPocketFinder). Thereafter, molecular docking (Extra Precision and Induced Fit Docking), and long range molecular dynamics simulation analysis, and ADME profile analysis have been carried out for screening purposes. Calculations of MM-PBSA based binding free energy (ΔG), and Density Functional Theory for quantum chemical simulations have been carried out for the hit compounds obtained through multi-step molecular docking based virtual screening technique. The multi-chemometric studies suggested that hit modulators have formed significant numbers of molecular interactions with hotspot residues in the homo-dimer interface region, which enable to hold CIP2A binding stability. Compounds with average ΔG values (ranging –3.4 x 10−2 to −1.1 x 10−2 KJ/mol) signifying promising CIP2A modulators. Communicated by Ramaswamy H. Sarma

Diameter Is a Key 3D Characteristic for Assessments of Efficient Inhibitors of Protein-Protein Interactions.

Three-dimensional (3D) molecular descriptors, including physicochemical and shape properties, for protein-protein interaction (PPI) interface inhibitors have become a topic of discussion. However, the relationships between such properties and binding free energy have not been adequately investigated. In this study, we focused on identifying key 3D molecular descriptors related to the binding free energy and/or the ligand efficiency (LE) of PPI interface inhibitors. A positive correlation was found between the binding free energy and the diameter (D) of cylindrical 3D molecules, in addition to a correlation between LE and D/heavy atom count (HAC). In addition, we showed a correlation between LE and D/HAC for macrocyclic compounds, suggesting that the present findings could be applied during assessments of the potential of macrocyclic PPI interface inhibitors in drug discovery processes.

Discovery of Protein-Protein Interaction Inhibitors by Integrating Protein Engineering and Chemical Screening Platforms

Protein-protein interactions (PPIs) govern intracellular life, and identification of PPI inhibitors is challenging. Roadblocks in assay development stemming from weak binding affinities of natural PPIs impede progress in this field. We postulated that enhancing binding affinity of natural PPIs via protein engineering will aid assay development and hit discovery. This proof-of-principle study targets PPI between linear ubiquitin chains and NEMO UBAN domain, which activates NF-κB signaling. Using phage display, we generated ubiquitin variants that bind to the functional UBAN epitope with high affinity, act as competitive inhibitors, and structurally maintain the existing PPI interface. When utilized in assay development, variants enable generation of robust cell-based assays for chemical screening. Top compounds identified using this approach directly bind to UBAN and dampen NF-κB signaling. This study illustrates advantages of integrating protein engineering and chemical screening in hit identification, a development that we anticipate will have wide application in drug discovery.

Mapping Protein-Protein Interaction Interface Peptides with Jun-Fos Assisted Phage Display and Deep Sequencing.

Protein-protein interactions govern many cellular processes, and identifying binding interaction sites on proteins can facilitate the discovery of inhibitors to block such interactions. Here we identify peptides from a randomly fragmented plasmid encoding the β-lactamase inhibitory protein (BLIP) and the Lac repressor (LacI) that represent regions of protein-protein interactions. We utilized a Jun-Fos-assisted phage display system that has previously been used to screen cDNA and genomic libraries to identify antibody antigens. Affinity selection with polyclonal antibodies against LacI or BLIP resulted in the rapid enrichment of in-frame peptides from various regions of the proteins. Further, affinity selection with β-lactamase enriched peptides that encompass regions of BLIP previously shown to contribute strongly to the binding energy of the BLIP/β-lactamase interaction, i.e., hotspot residues. Further, one of the regions enriched by affinity selection encompassed a disulfide-constrained region of BLIP that forms part of the BLIP interaction surface in the native complex that we show also binds to β-lactamase as a disulfide-constrained macrocycle peptide with a KD of ∼1 μM. Fragmented open reading frame (ORF) libraries may efficiently identify such naturally constrained peptides at protein-protein interaction interfaces. With sufficiently deep coverage of ORFs by peptide-coding inserts, phage display and deep sequencing can provide detailed information on the domains or peptides that contribute to an interaction. Such information should enable the design of potentially therapeutic macrocycles or peptidomimetics that block the interaction.

A Novel Index of Contact Frequency from Noise Protein-Protein Interaction Data Help for Accurate Interface Residue Pair Prediction.

Protein-protein interactions are important for most biological processes and have been studied for decades. However, the detailed formation mechanism of protein-protein interaction interface is still ambiguous, which makes it difficult to accurately predict the protein-protein interaction interface residue pairs. Here, we extract the interface residue-residue contacts from the decoys in the ZDOCK protein-protein complex decoy set with RMSD mostly larger than 3Å. To accurately compute the interface residue-residue contacts, we define a new constant called interface residue pairs frequency, which counts the atom contact numbers between two interface residues. We normalize interface residue pairs frequency to pick out the top residue-residue pairs from all the possible pairs preferential to be on correct protein-protein interaction interface. When tested on 37 protein dimers from the decoy set where most decoys are incorrect, our method successfully predicts 30 protein dimers with a success rate of up to 81.1%. Higher accuracy than some other state-of-the-art methods confirmed the performance of our method.

Interaction Energetics and Druggability of the Protein-Protein Interaction between Kelch-like ECH-Associated Protein 1 (KEAP1) and Nuclear Factor Erythroid 2 Like 2 (Nrf2).

Development of small molecule inhibitors of protein-protein interactions (PPIs) is hampered by our poor understanding of the druggability of PPI target sites. Here, we describe the combined application of alanine-scanning mutagenesis, fragment screening, and FTMap computational hot spot mapping to evaluate the energetics and druggability of the highly charged PPI interface between Kelch-like ECH-associated protein 1 (KEAP1) and nuclear factor erythroid 2 like 2 (Nrf2), an important drug target. FTMap identifies four binding energy hot spots at the active site. Only two of these are exploited by Nrf2, which alanine scanning of both proteins shows to bind primarily through E79 and E82 interacting with KEAP1 residues S363, R380, R415, R483, and S508. We identify fragment hits and obtain X-ray complex structures for three fragments via crystal soaking using a new crystal form of KEAP1. Combining these results provides a comprehensive and quantitative picture of the origins of binding energy at the interface. Our findings additionally reveal non-native interactions that might be exploited in the design of uncharged synthetic ligands to occupy the same site on KEAP1 that has evolved to bind the highly charged DEETGE binding loop of Nrf2. These include π-stacking with KEAP1 Y525 and interactions at an FTMap-identified hot spot deep in the binding site. Finally, we discuss how the complementary information provided by alanine-scanning mutagenesis, fragment screening, and computational hot spot mapping can be integrated to more comprehensively evaluate PPI druggability.

Structures of the AMPA receptor in complex with its auxiliary subunit cornichon.

In the brain, AMPA-type glutamate receptors (AMPARs) form complexes with their auxiliary subunits and mediate the majority of fast excitatory neurotransmission. Signals transduced by these complexes are critical for synaptic plasticity, learning, and memory. The two major categories of AMPAR auxiliary subunits are transmembrane AMPAR regulatory proteins (TARPs) and cornichon homologs (CNIHs); these subunits share little homology and play distinct roles in controlling ion channel gating and trafficking of AMPAR. Here, I report high-resolution cryo-electron microscopy structures of AMPAR in complex with CNIH3. Contrary to its predicted membrane topology, CNIH3 lacks an extracellular domain and instead contains four membrane-spanning helices. The protein-protein interaction interface that dictates channel modulation and the lipids surrounding the complex are revealed. These structures provide insights into the molecular mechanism for ion channel modulation and assembly of AMPAR/CNIH3 complexes.

DUF3669, a \u201cdomain of unknown function\u201d within ZNF746 and ZNF777, oligomerizes and contributes to transcriptional repression

BackgroundZNF746 and ZNF777 belong to a subset of the large Krüppel-associated box (KRAB) zinc finger (ZNF) transcription factor family. They contain, like four other members in human, an additional conserved domain, the “domain of unknown function 3669” (DUF3669). Previous work on members of this subfamily suggested involvement in transcriptional regulation and aberrant ZNF746 overexpression leads to neuronal cell death in Parkinson’s disease.ResultsHere we demonstrate that N-terminal protein segments of the ZNF746a major isoform and ZNF777 act in concert to exert moderate transcriptional repression activities. Full potency depended on the intact configuration consisting of DUF3669, a variant KRAB domain and adjacent sequences. While DUF3669 contributes an intrinsic weak inhibitory activity, the isolated KRAB-AB domains did not repress. Importantly, DUF3669 provides a novel protein-protein interaction interface and mediates direct physical interaction between the members of the subfamily in oligomers. The ZNF746 protein segment encoded by exons 5 and 6 boosted repressor potency, potentially due to the presence of an acceptor lysine for sumoylation at K189. Repressor activity of the potent canonical ZNF10 KRAB domain was not augmented by heterologous transfer of DUF3669, pointing to the importance of context for DUF3669’s impact on transcription. Neither ZNF746a nor ZNF777 protein segments stably associated with TRIM28 within cells. Isoform ZNF746b that contains, unlike the major isoform, a full-length KRAB-A subdomain, displayed substantially increased repressor potency. This increase is due to canonical mechanisms known for KRAB domains since it did not take place in HAP1 knockout models of TRIM28 and SETDB1. A glycine to glutamic acid replacement that complies with a bona fide conserved “MLE” sequence within KRAB-A led to a further strong gain in repressor potency to levels comparable to those of the canonical ZNF10 KRAB domain. Each gain of repressive activity was accompanied by an enhanced interaction with TRIM28 protein.ConclusionDUF3669 adds a protein-protein interaction surface to a subgroup of KRAB-ZNF proteins within an N-terminal configuration with variant KRAB and adjacent sequences likely regulated by sumoylation. DUF3669 contributes to transcriptional repression strength and its homo- and hetero-oligomerization characteristics probably extended the regulatory repertoire of KRAB-ZNF transcription factors during amniote evolution.

Attention mechanism enhanced LSTM with residual architecture and its application for protein-protein interaction residue pairs prediction

BackgroundRecurrent neural network(RNN) is a good way to process sequential data, but the capability of RNN to compute long sequence data is inefficient. As a variant of RNN, long short term memory(LSTM) solved the problem in some extent. Here we improved LSTM for big data application in protein-protein interaction interface residue pairs prediction based on the following two reasons. On the one hand, there are some deficiencies in LSTM, such as shallow layers, gradient explosion or vanishing, etc. With a dramatic data increasing, the imbalance between algorithm innovation and big data processing has been more serious and urgent. On the other hand, protein-protein interaction interface residue pairs prediction is an important problem in biology, but the low prediction accuracy compels us to propose new computational methods.ResultsIn order to surmount aforementioned problems of LSTM, we adopt the residual architecture and add attention mechanism to LSTM. In detail, we redefine the block, and add a connection from front to back in every two layers and attention mechanism to strengthen the capability of mining information. Then we use it to predict protein-protein interaction interface residue pairs, and acquire a quite good accuracy over 72%. What’s more, we compare our method with random experiments, PPiPP, standard LSTM, and some other machine learning methods. Our method shows better performance than the methods mentioned above.ConclusionWe present an attention mechanism enhanced LSTM with residual architecture, and make deeper network without gradient vanishing or explosion to a certain extent. Then we apply it to a significant problem– protein-protein interaction interface residue pairs prediction and obtain a better accuracy than other methods. Our method provides a new approach for protein-protein interaction computation, which will be helpful for related biomedical researches.

Investigation of protein-protein interactions and hot spot region between PD-1 and PD-L1 by fragment molecular orbital method

Inhibitors to interfere protein-protein interactions (PPI) between programmed cell death 1 (PD-1) and programmed death ligand-1 (PD-L1) block evasion of cancers from immune surveillance. Analyzing hot spot residues in PPI is important for small-molecule drug development. In order to find out hot spots on PPI interface in PD-1/PD-L1 complex, we analyzed PPI in PD-1/PD-L1 with a new analysis method, 3-dimensional scattered pair interactions energies (3D-SPIEs), which assorts significant interactions with fragment molecular orbital (FMO) method. By additionally analyzing PPI in PD-1/antibody and PD-L1/antibody complexes, and small-ligand interactions in PD-L1/peptide and PD-L1/small-molecule complexes, we narrowed down the hot spot region with 3D-SPIEs-based interaction map, which integrates PPI and small-ligand interactions. Based on the map, there are two hot spot regions in PPI of PD-1/PD-L1 and the first hot spot region is important for inhibitors. In particular, LY56, LE58, and LN66 in the first hot spot of PD-L1 are important for PD-L1-antibodies and small-inhibitors in common, while LM115 is important for small-inhibitors. Therefore, the 3D-SPIEs-based map would provide valuable information for designing new small-molecule inhibitors to inhibit PPI of PD-1/PD-L1 and the FMO/3D-SPIEs method provides an effectual tool to understand PPI and integrate PPI and small-ligand interactions at a quantum mechanical level.

Probing the Druggablility on the Interface of the Protein-Protein Interaction and Its Allosteric Regulation Mechanism on the Drug Screening for the CXCR4 Homodimer.

Modulating protein–protein interactions (PPIs) with small drug-like molecules targeting it exhibits great promise in modern drug discovery. G protein-coupled receptors (GPCRs) are the largest family of targeted proteins and could form dimers in living biological cells through PPIs. However, compared to drug development of the orthosteric site, there has been lack of investigations on the druggability of the PPI interface for GPCRs and its functional implication on experiments. Thus, in order to address these issues, we constructed a novel computational strategy, which involved in molecular dynamics simulation, virtual screening and protein structure network (PSN), to study one representative GPCR homodimer (CXCR4). One druggable pocket was identified in the PPI interface and one small molecule targeting it was screened, which could strengthen PPI mainly through hydrophobic interaction between the benzene rings of the PPI molecule and TM4 of the receptor. The PSN results further reveals that the PPI molecule could increase the number of the allosteric regulation pathways between the druggable pocket of the dimer interface to the orthostatic site for the subunit A but only play minor role for the other subunit B, leading to the asymmetric change in the volume of the binding pockets for the two subunits (increase for the subunit A and minor change for the subunit B). Consequently, the screening performance of the subunit A to the antagonists is enhanced while the subunit B is unchanged nearly, implying that the PPI molecule may be beneficial to enhance the drug efficacies of the antagonists. In addition, one main regulation pathway with the highest frequency was identified for the subunit A, which consists of Trp1955.34–Tyr190ECL2–Val1965.35–Gln2005.39–Asp2626.58–Cys28N-term, revealing their importance in the allosteric regulation from the PPI molecule. The observations from the work could provide valuable information for the development of the PPI drug-like molecule for GPCRs.