Protein-protein Binding Interface Research Articles

Conserved biophysical compatibility among the highly variable germline-encoded regions shapes TCR-MHC interactions.

T cells are critically important components of the adaptive immune system primarily responsible for identifying and responding to pathogenic challenges. This recognition of pathogens is driven by the interaction between membrane-bound T cell receptors (TCRs) and antigenic peptides presented on major histocompatibility complex (MHC) molecules. The formation of the TCR-peptide-MHC complex (TCR-pMHC) involves interactions among germline-encoded and hypervariable amino acids. Germline-encoded and hypervariable regions can form contacts critical for complex formation, but only interactions between germline-encoded contacts are likely to be shared across many of all the possible productive TCR-pMHC complexes. Despite this, experimental investigation of these interactions have focused on only a small fraction of the possible interaction space. To address this, we analyzed every possible germline-encoded TCR-MHC contact in humans, thereby generating the first comprehensive characterization of these largely antigen-independent interactions. Our computational analysis suggests that germline-encoded TCR-MHC interactions that are conserved at the sequence level are rare due to the high amino acid diversity of the TCR CDR1 and CDR2 loops, and that such conservation is unlikely to dominate the dynamic protein-protein binding interface. Instead, we propose that binding properties such as the docking orientation are defined by regions of biophysical compatibility between these loops and the MHC surface.

Binding blockade between TLN1 and integrin β1 represses triple-negative breast cancer.

Integrin family are known as key gears in focal adhesion for triple-negative breast cancer (TNBC) metastasis. However, the integrin independent factor TLN1 remains vague in TNBC. Bioinformatics analysis was performed based on TCGA database and Shengjing Hospital cohort. Western blot and RT-PCR were used to detect the expression of TLN1 and integrin pathway in cells. A small-molecule C67399 was screened for blocking TLN1 and integrin β1 through a novel computational screening approach by targeting the protein-protein binding interface. Drug pharmacodynamics were determined through xenograft assay. Upregulation of TLN1 in TNBC samples correlates with metastasis and worse prognosis. Silencing TLN1 in TNBC cells significantly attenuated the migration of tumour cells through interfering the dynamic formation of focal adhesion with integrin β1, thus regulating FAK-AKT signal pathway and epithelial-mesenchymal transformation. Targeting the binding between TLN1 and integrin β1 by C67399 could repress metastasis of TNBC. TLN1 overexpression contributes to TNBC metastasis and C67399 targeting TLN1 may hold promise for TNBC treatment. This study was supported by grants from the National Natural Science Foundation of China (No. 81872159, 81902607, 81874301), Liaoning Colleges Innovative Talent Support Program (Name: Cancer Stem Cell Origin and Biological Behaviour), Outstanding Scientific Fund of Shengjing Hospital (201803), and Outstanding Young Scholars of Liaoning Province (2019-YQ-10).

Designing efficient genetic code expansion in Bacillus subtilis to gain biological insights

Bacillus subtilis is a model gram-positive bacterium, commonly used to explore questions across bacterial cell biology and for industrial uses. To enable greater understanding and control of proteins in B. subtilis, here we report broad and efficient genetic code expansion in B. subtilis by incorporating 20 distinct non-standard amino acids within proteins using 3 different families of genetic code expansion systems and two choices of codons. We use these systems to achieve click-labelling, photo-crosslinking, and translational titration. These tools allow us to demonstrate differences between E. coli and B. subtilis stop codon suppression, validate a predicted protein-protein binding interface, and begin to interrogate properties underlying bacterial cytokinesis by precisely modulating cell division dynamics in vivo. We expect that the establishment of this simple and easily accessible chemical biology system in B. subtilis will help uncover an abundance of biological insights and aid genetic code expansion in other organisms.

Rif1 regulates telomere length through conserved HEAT repeats.

In budding yeast, Rif1 negatively regulates telomere length, but the mechanism of this regulation has remained elusive. Previous work identified several functional domains of Rif1, but none of these has been shown to mediate telomere length. To define Rif1 domains responsible for telomere regulation, we localized truncations of Rif1 to a single specific telomere and measured telomere length of that telomere compared to bulk telomeres. We found that a domain in the N-terminus containing HEAT repeats, Rif1177–996, was sufficient for length regulation when tethered to the telomere. Charged residues in this region were previously proposed to mediate DNA binding. We found that mutation of these residues disrupted telomere length regulation even when Rif1 was tethered to the telomere. Mutation of other conserved residues in this region, which were not predicted to interact with DNA, also disrupted telomere length maintenance, while mutation of conserved residues distal to this region did not. Our data suggest that conserved amino acids in the region from 436 to 577 play a functional role in telomere length regulation, which is separate from their proposed DNA binding function. We propose that the Rif1 HEAT repeats region represents a protein-protein binding interface that mediates telomere length regulation.

Structure and Role of BCOR PUFD in Noncanonical PRC1 Assembly and Disease.

Polycomb repression complex 1 (PRC1) is a multiprotein assembly that regulates transcription. The Polycomb group ring finger 1 protein (PCGF1) is central in the assembly of the noncanonical PRC1 variant called PRC1.1 through its direct interaction with BCOR (BCL-6-interacting corepressor) or its paralog, BCOR-like 1 (BCORL1). Previous structural studies revealed that the C-terminal PUFD domain of BCORL1 is necessary and sufficient to heterodimerize with the RAWUL domain of PCGF1 and, together, form a new protein-protein binding interface that associates with the histone demethylase KDM2B. Here, we show that the PUFD of BCOR and BCORL1 differ in their abilities to assemble with KDM2B. Unlike BCORL1, the PUFD of BCOR alone does not stably assemble with KDM2B. Rather, additional residues N-terminal to the BCOR PUFD are necessary for stable association. Nuclear magnetic resonance (NMR) structure determination and 15N T2 relaxation time measurements of the BCOR PUFD alone indicate that the termini of the BCOR PUFD, which are critical for binding PCGF1 and KDM2B, are disordered. This suggests a hierarchical mode of assembly whereby BCOR PUFD termini become structurally ordered upon binding PCGF1, which then allows stable association with KDM2B. Notably, BCOR internal tandem duplications (ITDs) leading to pediatric kidney and brain tumors map to the PUFD termini. Binding studies with the BCOR ITD indicate the ITD would disrupt PRC1.1 assembly, suggesting loss of the ability to assemble PRC1.1 is a critical molecular event driving tumorigenesis.

Probing Protein-Protein Interactions Using Asymmetric Labeling and Carbonyl-Carbon Selective Heteronuclear NMR Spectroscopy.

Protein-protein interactions (PPIs) regulate a plethora of cellular processes and NMR spectroscopy has been a leading technique for characterizing them at the atomic resolution. Technically, however, PPIs characterization has been challenging due to multiple samples required to characterize the hot spots at the protein interface. In this paper, we review our recently developed methods that greatly simplify PPI studies, which minimize the number of samples required to fully characterize residues involved in the protein-protein binding interface. This original strategy combines asymmetric labeling of two binding partners and the carbonyl-carbon label selective (CCLS) pulse sequence element implemented into the heteronuclear single quantum correlation (1H-15N HSQC) spectra. The CCLS scheme removes signals of the J-coupled 15N–13C resonances and records simultaneously two individual amide fingerprints for each binding partner. We show the application to the measurements of chemical shift correlations, residual dipolar couplings (RDCs), and paramagnetic relaxation enhancements (PRE). These experiments open an avenue for further modifications of existing experiments facilitating the NMR analysis of PPIs.

Structural Characterization and Function Prediction of Immunoglobulin-like Fold in Cell Adhesion and Cell Signaling.

Domains that belong to an immunoglobulin (Ig) fold are extremely abundant in cell surface receptors, which play significant roles in cell-cell adhesion and signaling. Although the structures of domains in an Ig fold share common topology of β-barrels, functions of receptors in adhesion and signaling are regulated by the very heterogeneous binding between these domains. Additionally, only a small number of domains are directly involved in the binding between two multidomain receptors. It is challenging and time consuming to experimentally detect the binding partners of a given receptor and further determine which specific domains in this receptor are responsible for binding. Therefore, current knowledge in the binding mechanism of Ig-fold domains and their impacts on cell adhesion and signaling is very limited. A bioinformatics study can shed light on this topic from a systematic point of view. However, there is so far no computational analysis on the structural and functional characteristics of the entire Ig fold. We constructed nonredundant structural data sets for all domains in Ig fold, depending on their functions in cell adhesion and signaling. We found that data sets of domains in adhesion receptors show different binding preference from domains in signaling receptors. Using structural alignment, we further built a common structural template for each group of a domain data set. By mapping the protein-protein binding interface of each domain in a group onto the surface of its structural template, we found binding interfaces are highly overlapped within each specific group. These overlapped interfaces, we called consensus binding interfaces, are distinguishable among different data sets of domains. Finally, the residue compositions on the consensus interfaces were used as indicators for multiple machine learning algorithms to predict if they can form homotypic interactions with each other. The overall performance of the cross-validation tests shows that our prediction accuracies ranged between 0.6 and 0.8.

Native Hydrophobic Binding Interactions at the Transition State for Association between the TAZ1 Domain of CBP and the Disordered TAD-STAT2 Are Not a Requirement.

A significant fraction of the eukaryotic proteome consists of proteins that are either partially or completely disordered under native-like conditions. Intrinsically disordered proteins (IDPs) are common in protein-protein interactions and are involved in numerous cellular processes. Although many proteins have been identified as disordered, much less is known about the binding mechanisms of the coupled binding and folding reactions involving IDPs. Here we have analyzed the rate-limiting transition state for binding between the TAZ1 domain of CREB binding protein and the intrinsically disordered transactivation domain of STAT2 (TAD-STAT2) by site-directed mutagenesis and kinetic experiments (Φ-value analysis) and found that the native protein-protein binding interface is not formed at the transition state for binding. Instead, native hydrophobic binding interactions form late, after the rate-limiting barrier has been crossed. The association rate constant in the absence of electrostatic enhancement was determined to be rather high. This is consistent with the Φ-value analysis, which showed that there are few or no obligatory native contacts. Also, linear free energy relationships clearly demonstrate that native interactions are cooperatively formed, a scenario that has usually been observed for proteins that fold according to the so-called nucleation-condensation mechanism. Thus, native hydrophobic binding interactions at the rate-limiting transition state for association between TAD-STAT2 and TAZ1 are not a requirement, which is generally in agreement with previous findings on other IDP systems and might be a common mechanism for IDPs.

Insights from the First Principles Based Large Scale Protein Thermostability Calculations

Proteins under evolutionary pressure have acquired a subtle balance of thermodynamic properties that optimizes functioning at a physiological range of the environmental conditions. The fields of biotechnology and bioengineering, however, aim to recalibrate, and sometimes design, biological molecules by tailoring certain physico-chemical properties. Enhancing protein thermal stability is one of the frequently desired features both by the industry and scientific community. Increased thermostability may be achieved by introducing mutations into a protein, and identification of such mutations is thus of key interest. In the current study we utilize a recently developed automated approach [1] to perform a large scale protein mutation scan assessing changes in thermostability from free energy molecular dynamics simulations. A remarkable agreement with experimental data was achieved and further enhanced by adopting a consensus force field approach. Subsequently, we demonstrate the strength of the first principles based methods in capturing mutation induced free energy changes in a membrane protein (neurotensin receptor 1). Finally, gaining access to accurate protein thermostabilities enables us to estimate free energy changes in protein-protein binding upon an amino acid mutation. We report on a number of cases examplifying the use of first principles approaches for mutations at a protein-protein binding interface.[1] V. Gapsys, S. Michielssens, D. Seeliger, and B. L. de Groot. pmx: Automated protein structure and topology generation for alchemical perturbations. J. Comput. Chem., 36(5):348-354, 2015.

Effects of Membrane Mimetics on Cytochrome P450-Cytochrome b5 Interactions Characterized by NMR Spectroscopy

Mammalian cytochrome P450 (P450) is a membrane-bound monooxygenase whose catalytic activities require two electrons to be sequentially delivered from its redox partners: cytochrome b5 (cytb5) and cytochrome P450 reductase, both of which are membrane proteins. Although P450 functional activities are known to be affected by lipids, experimental evidence to reveal the effect of membrane on P450-cytb5 interactions is still lacking. Here, we present evidence for the influence of phospholipid bilayers on complex formation between rabbit P450 2B4 (CYP2B4) and rabbit cytb5 at the atomic level, utilizing NMR techniques. General line broadening and modest chemical shift perturbations of cytb5 resonances characterize CYP2B4-cytb5 interactions on the intermediate time scale. More significant intensity attenuation and a more specific protein-protein binding interface are observed in bicelles as compared with lipid-free solution, highlighting the importance of the lipid bilayer in stabilizing stronger and more specific interactions between CYP2B4 and cytb5, which may lead to a more efficient electron transfer. Similar results observed for the interactions between CYP2B4 lacking the transmembrane domain (tr-CYP2B4) and cytb5 imply interactions between tr-CYP2B4 and the membrane surface, which might assist in CYP2B4-cytb5 complex formation by orienting tr-CYP2B4 for efficient contact with cytb5. Furthermore, the observation of weak and nonspecific interactions between CYP2B4 and cytb5 in micelles suggests that lipid bilayer structures and low curvature membrane surface are preferable for CYP2B4-cytb5 complex formation. Results presented in this study provide structural insights into the mechanism behind the important role that the lipid bilayer plays in the interactions between P450s and their redox partners.

Combining Water Percolation Analysis and Molecular Dynamics Simulations for Protein-Protein Binding Interface Prediction

We consider ab initio prediction of protein-protein binding interfaces of weakly-associating proteins in aqueous environment. Systems of interest are important in cellular signalling pathways and manifest Kd values in the micromolar range (weak binding). We use molecular dynamics simulations of individual proteins at varying hydration levels to determine percolation thresholds for water networks. By mapping average local water density around all residues at the percolation threshold, we are able to determine clusters of residues with low water density which we consider as potential binding regions. We have performed analysis on a series of 20 proteins and compare to existing automated protein-protein binding interface prediction servers.

Protein docking prediction using predicted protein-protein interface

BackgroundMany important cellular processes are carried out by protein complexes. To provide physical pictures of interacting proteins, many computational protein-protein prediction methods have been developed in the past. However, it is still difficult to identify the correct docking complex structure within top ranks among alternative conformations.ResultsWe present a novel protein docking algorithm that utilizes imperfect protein-protein binding interface prediction for guiding protein docking. Since the accuracy of protein binding site prediction varies depending on cases, the challenge is to develop a method which does not deteriorate but improves docking results by using a binding site prediction which may not be 100% accurate. The algorithm, named PI-LZerD (using Predicted Interface with Local 3D Zernike descriptor-based Docking algorithm), is based on a pair wise protein docking prediction algorithm, LZerD, which we have developed earlier. PI-LZerD starts from performing docking prediction using the provided protein-protein binding interface prediction as constraints, which is followed by the second round of docking with updated docking interface information to further improve docking conformation. Benchmark results on bound and unbound cases show that PI-LZerD consistently improves the docking prediction accuracy as compared with docking without using binding site prediction or using the binding site prediction as post-filtering.ConclusionWe have developed PI-LZerD, a pairwise docking algorithm, which uses imperfect protein-protein binding interface prediction to improve docking accuracy. PI-LZerD consistently showed better prediction accuracy over alternative methods in the series of benchmark experiments including docking using actual docking interface site predictions as well as unbound docking cases.

Tethered-Hopping Model for Protein-DNA Binding and Unbinding Based on Sox2-Oct1-Hoxb1 Ternary Complex Simulations

The sliding and hopping models encapsulate the essential protein-DNA binding process for binary complex formation and dissociation. However, the effects of a cofactor protein on the protein-DNA binding process that leads to the formation of a ternary complex remain largely unknown. Here we investigate the effect of the cofactor Sox2 on the binding and unbinding of Oct1 with the Hoxb1 control element. We simulate the association of Oct1 with Sox2-Hoxb1 using molecular dynamics simulations, and the dissociation of Oct1 from Sox2-Hoxb1 using steered molecular dynamics simulations, in analogy to a hopping event of Oct1. We compare the kinetic and thermodynamic properties of three model complexes (the wild-type and two mutants) in which the Oct1-DNA base-specific interactions or the Sox2-Oct1 protein-protein interactions are largely abolished. We find that Oct1-DNA base-specific interactions contribute significantly to the total interaction energy of the ternary complex, and that nonspecific Oct1-DNA interactions are sufficient for driving the formation of the protein-DNA interface. The Sox2-Oct1 protein-protein binding interface is largely hydrophobic, with remarkable shape complementarity. This interface promotes the formation of the ternary complex and slows the dissociation of Oct1 from its DNA-binding site. We propose a simple two-step reaction model of protein-DNA binding, called the tethered-hopping model, that explains the importance of the cofactor Sox2 and may apply to similar ternary protein-DNA complexes.

Conformational Plasticity in IL-2 Is Captured Via in Silico Fragment-Based Drug Design Using SILCS Fragmaps

The IL-2 protein is an experimentally well-characterized example of conformational plasticity in a protein-protein binding interface. Whereas the original crystal structure of the IL-2 protein lacked a particular binding pocket at the interface, subsequent structures of the protein complexed with small-molecule inhibitors have shown the creation of a hydrophobic pocket in order to accommodate inhibitor binding poses. A new in silico fragment-based drug design methodology (SILCS FragMaps) was applied to the original IL-2 crystal structure and has been found to recapitulate this pocket observed in the complexes.

Selection and analysis of an optimized anti-VEGF antibody: crystal structure of an affinity-matured fab in complex with antigen

The Fab portion of a humanized antibody (Fab-12; IgG form known as rhuMAb VEGF) to vascular endothelial growth factor (VEGF) has been affinity-matured through complementarity-determining region (CDR) mutation, followed by affinity selection using monovalent phage display. After stringent binding selections at 37°C, with dissociation (off-rate) selection periods of several days, high affinity variants were isolated from CDR-H1, H2, and H3 libraries. Mutations were combined to obtain cumulatively tighter-binding variants. The final variant identified here, Y0317, contained six mutations from the parental antibody. In vitro cell-based assays show that four mutations yielded an improvement of about 100-fold in potency for inhibition of VEGF-dependent cell proliferation by this variant, consistent with the equilibrium binding constant determined from kinetics experiments at 37°C. Using X-ray crystallography, we determined a high-resolution structure of the complex between VEGF and the affinity-matured Fab fragment. The overall features of the binding interface seen previously with wild-type are preserved, and many contact residues are maintained in precise alignment in the superimposed structures. However, locally, we see evidence for improved contacts between antibody and antigen, and two mutations result in increased van der Waals contact and improved hydrogen bonding. Site-directed mutants confirm that the most favorable improvements as judged by examination of the complex structure, in fact, have the greatest impact on free energy of binding. In general, the final antibody has improved affinity for several VEGF variants as compared with the parental antibody; however, some contact residues on VEGF differ in their contribution to the energetics of Fab binding. The results show that small changes even in a large protein-protein binding interface can have significant effects on the energetics of interaction.