Specificity and stability of transient protein-protein interactions

OGA mutant aberrantly hydrolyzes O-GlcNAc modification from PDLIM7 to modulate p53 and cytoskeleton in promoting cancer cell malignancy

O-GlcNAcase (OGA) is the only human enzyme that catalyzes the hydrolysis (deglycosylation) of O-linked beta-N-acetylglucosaminylation (O-GlcNAcylation) from numerous protein substrates. OGA has broad implications in many challenging diseases including cancer. However, its role in cell malignancy remains mostly unclear. Here, we report that a cancer-derived point mutation on the OGA's noncatalytic stalk domain aberrantly modulates OGA interactome and substrate deglycosylation toward a specific set of proteins. Interestingly, our quantitative proteomic studies uncovered that the OGA stalk domain mutant preferentially deglycosylated protein substrates with +2 proline in the sequence relative to the O-GlcNAcylation site. One of the most dysregulated substrates is PDZ and LIM domain protein 7 (PDLIM7), which is associated with the tumor suppressor p53. We found that the aberrantly deglycosylated PDLIM7 suppressed p53 gene expression and accelerated p53 protein degradation by promoting the complex formation with E3 ubiquitin ligase MDM2. Moreover, deglycosylated PDLIM7 significantly up-regulated the actin-rich membrane protrusions on the cell surface, augmenting the cancer cell motility and aggressiveness. These findings revealed an important but previously unappreciated role of OGA's stalk domain in protein substrate recognition and functional modulation during malignant cell progression.

Interfacial residues in protein-protein complexes are in the eyes of the beholder

Interactions between proteins are vital in almost all biological processes. The characterization of protein-protein interactions helps us understand the mechanistic basis of biological processes, thereby enabling the manipulation of proteins for biotechnological and clinical purposes. The interface residues of a protein-protein complex are assumed to have the following two properties: (a) they always interact with a residue of a partner protein, which forms the basis for distance-based interface residue identification methods, and (b) they are solvent-exposed in the isolated form of the protein and become buried in the complex form, which forms the basis for Accessible Surface Area (ASA)-based methods. The study interrogates this popular assumption by recognizing interface residues in protein-protein complexes through these two methods. The results show that a few residues are identified uniquely by each method, and the extent of conservation, propensities, and their contribution to the stability of protein-protein interaction varies substantially between these residues. The case study analyses showed that interface residues, unique to distance, participate in crucial interactions that hold the proteins together, whereas the interface residues unique to the ASA method have a potential role in the recognition, dynamics, and specificity of the complex and can also be a hotspot. Overall, the study recommends applying both distance and ASA methods so that some interface residues missed by either method but crucial to the stability, recognition, dynamics, and function of protein-protein complexes are identified in a complementary manner.

Iterative In Silico Screening for Optimizing Stable Conformation of Anti-SARS-CoV-2 Nanobodies

Nanobodies (Nbs or VHHs) are single-domain antibodies (sdAbs) derived from camelid heavy-chain antibodies. Nbs have special and unique characteristics, such as small size, good tissue penetration, and cost-effective production, making Nbs a good candidate for the diagnosis and treatment of viruses and other pathologies. Identifying effective Nbs against COVID-19 would help us control this dangerous virus or other unknown variants in the future. Herein, we introduce an in silico screening strategy for optimizing stable conformation of anti-SARS-CoV-2 Nbs. Firstly, various complexes containing nanobodies were downloaded from the RCSB database, which were identified from immunized llamas. The primary docking between Nbs and the SARS-CoV-2 spike protein receptor-binding domain was performed through the ClusPro program, with the manual screening leaving the reasonable conformation to the next step. Then, the binding distances of atoms between the antigen-antibody interfaces were measured through the NeighborSearch algorithm. Finally, filtered nanobodies were acquired according to HADDOCK scores through HADDOCK docking the COVID-19 spike protein with nanobodies under restrictions of calculated molecular distance between active residues and antigenic epitopes less than 4.5 Å. In this way, those nanobodies with more reasonable conformation and stronger neutralizing efficacy were acquired. To validate the efficacy ranking of the nanobodies we obtained, we calculated the binding affinities (∆G) and dissociation constants (Kd) of all screened nanobodies using the PRODIGY web tool and predicted the stability changes induced by all possible point mutations in nanobodies using the MAESTROWeb server. Furthermore, we examined the performance of the relationship between nanobodies' ranking and their number of mutation-sensitive sites (Spearman correlation > 0.68); the results revealed a robust correlation, indicating that the superior nanobodies identified through our screening process exhibited fewer mutation hotspots and higher stability. This correlation analysis demonstrates the validity of our screening criteria, underscoring the suitability of these nanobodies for future development and practical implementation. In conclusion, this three-step screening strategy iteratively in silico greatly improved the accuracy of screening desired nanobodies compared to using only ClusPro docking or default HADDOCK docking settings. It provides new ideas for the screening of novel antibodies and computer-aided screening methods.

Distinct or Overlapping Areas of Mitochondrial Thioredoxin 2 May Be Used for Its Covalent and Strong Non-Covalent Interactions with Protein Ligands

In silico approaches were employed to examine the characteristics of interactions between human mitochondrial thioredoxin 2 (HsTrx2) and its 38 previously identified mitochondrial protein ligands. All interactions appeared driven mainly by electrostatic forces. The statistically significant residues of HsTrx2 for interactions were characterized as "contact hot spots". Since these were identical/adjacent to putative thermodynamic hot spots, an energy network approach identified their neighbors to highlight possible contact interfaces. Three distinct areas for binding emerged: (i) one around the active site for covalent interactions, (ii) another antipodal to the active site for strong non-covalent interactions, and (iii) a third area involved in both kinds of interactions. The contact interfaces of HsTrx2 were projected as respective interfaces for Escherichia coli Trx1 (EcoTrx1), 2, and HsTrx1. Comparison of the interfaces and contact hot spots of HsTrx2 to the contact residues of EcoTx1 and HsTrx1 from existing crystal complexes with protein ligands supported the hypothesis, except for a part of the cleft/groove adjacent to Trp30 preceding the active site. The outcomes of this study raise the possibility for the rational design of selective inhibitors for the interactions of HsTrx2 with specific protein ligands without affecting the entirety of the functions of the Trx system.

Computational analysis of protein synthesis, diffusion, and binding in compartmental biochips

Protein complex assembly facilitates the combination of individual protein subunits into functional entities, and thus plays a crucial role in biology and biotechnology. Recently, we developed quasi-twodimensional, silicon-based compartmental biochips that are designed to study and administer the synthesis and assembly of protein complexes. At these biochips, individual protein subunits are synthesized from locally confined high-density DNA brushes and are captured on the chip surface by molecular traps. Here, we investigate single-gene versions of our quasi-twodimensional synthesis systems and introduce the trap-binding efficiency to characterize their performance. We show by mathematical and computational modeling how a finite trap density determines the dynamics of protein-trap binding and identify three distinct regimes of the trap-binding efficiency. We systematically study how protein-trap binding is governed by the system's three key parameters, which are the synthesis rate, the diffusion constant and the trap-binding affinity of the expressed protein. In addition, we describe how spatially differential patterns of traps modulate the protein-trap binding dynamics. In this way, we extend the theoretical knowledge base for synthesis, diffusion, and binding in compartmental systems, which helps to achieve better control of directed molecular self-assembly required for the fabrication of nanomachines for synthetic biology applications or nanotechnological purposes.

Cancer-derived mutation in the OGA stalk domain promotes cell malignancy through dysregulating PDLIM7 and p53

O-GlcNAcase (OGA) is the sole enzyme that hydrolyzes O-GlcNAcylation from thousands of proteins and is dysregulated in many diseases including cancer. However, the substrate recognition and pathogenic mechanisms of OGA remain largely unknown. Here we report the first discovery of a cancer-derived point mutation on the OGA's non-catalytic stalk domain that aberrantly regulated a small set of OGA-protein interactions and O-GlcNAc hydrolysis in critical cellular processes. We uncovered a novel cancer-promoting mechanism in which the OGA mutant preferentially hydrolyzed the O-GlcNAcylation from modified PDLIM7 and promoted cell malignancy by down-regulating p53 tumor suppressor in different types of cells through transcription inhibition and MDM2-mediated ubiquitination. Our study revealed the OGA deglycosylated PDLIM7 as a novel regulator of p53-MDM2 pathway, offered the first set of direct evidence on OGA substrate recognition beyond its catalytic site, and illuminated new directions to interrogate OGA's precise role without perturbing global O-GlcNAc homeostasis for biomedical applications.

Tight Complex Formation of the Fumarate Sensing DcuS-DcuR Two-Component System at the Membrane and Target Promoter Search by Free DcuR Diffusion

Signaling of two-component systems by phosphoryl transfer requires interaction of the sensor kinase with the response regulator. Interaction of the C4-dicarboxylate-responsive and membrane-integral sensor kinase DcuS with the response regulator DcuR was studied. In vitro, the cytoplasmic part of DcuS (PAS_C-Kin) was employed. Stable complexes were formed, when either DcuS or DcuR were phosphorylated (K_d 22 ± 11 and 28 ± 7 nM, respectively). The unphosphorylated proteins produced a more labile complex (K_d 1380 ± 395 nM). Bacterial two-hybrid studies confirm interaction of DcuR with DcuS (and PAS_C-Kin) in vivo. The absolute contents of DcuR (197-979 pmol mg⁻¹ protein) in the bacteria exceeded those of DcuS by more than 1 order of magnitude. According to the K_d values, DcuS exists in complex, with phosphorylated but also unphosphorylated DcuR. In live cell imaging, the predominantly freely diffusing DcuR becomes markedly less mobile after phosphorylation and activation of DcuS by fumarate. Portions of the low mobility fraction accumulated at the cell poles, the preferred location of DcuS, and other portions within the cell, representing phosphorylated DcuR bound to promoters. In the model, acitvation of DcuS increases the affinity toward DcuR, leading to DcuS-P × DcuR formation and phosphorylation of DcuR. The complex is stable enough for phosphate-transfer, but labile enough to allow exchange between DcuR from the cytosol and DcuR-P of the complex. Released DcuR-P diffuses to target promoters and binds. Uncomplexed DcuR-P in the cytosol binds to nonactivated DcuS and becomes dephosphorylated. The lower affinity between DcuR and DcuS avoids blocking of DcuS and allows rapid exchange of DcuR.

IMPORTANCE Complex formation of membrane-bound sensor kinases with the response regulators represents an inherent step of signaling from the membrane to the promoters on the DNA. In the C4-dicarboxylate-sensing DcuS-DcuR two-component system, complex formation is strengthened by activation (phosphorylation) in vitro and in vivo, with trapping of the response regulator DcuR at the membrane. Single-molecule tracking of DcuR in the bacterial cell demonstrates two populations of DcuR with decreased mobility in the bacteria after activation: one at the membrane, but a second in the cytosol, likely representing DNA-bound DcuR. The data suggest a model with binding of DcuR to DcuS-P for phosphorylation, and of DcuR-P to DcuS for dephosphorylation, allowing rapid adaptation of the DcuR phosphorylation state. DcuR-P is released and transferred to DNA by 3D diffusion.

Are transient protein-protein interactions more dispensable?

Protein-protein interactions (PPIs) are key drivers of cell function and evolution. While it is widely assumed that most permanent PPIs are important for cellular function, it remains unclear whether transient PPIs are equally important. Here, we estimate and compare dispensable content among transient PPIs and permanent PPIs in human. Starting with a human reference interactome mapped by experiments, we construct a human structural interactome by building three-dimensional structural models for PPIs, and then distinguish transient PPIs from permanent PPIs using several structural and biophysical properties. We map common mutations from healthy individuals and disease-causing mutations onto the structural interactome, and perform structure-based calculations of the probabilities for common mutations (assumed to be neutral) and disease mutations (assumed to be mildly deleterious) to disrupt transient PPIs and permanent PPIs. Using Bayes' theorem we estimate that a similarly small fraction (<~20%) of both transient and permanent PPIs are completely dispensable, i.e., effectively neutral upon disruption. Hence, transient and permanent interactions are subject to similarly strong selective constraints in the human interactome.

Modeling the Dynamics of Protein–Protein Interfaces, How and Why?

Protein–protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein–protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein–protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein–protein interface.

Structure and thermodynamics of transient protein-protein complexes by chemometric decomposition of SAXS datasets

Transient biomolecular interactions play crucial roles in many cellular signaling and regulation processes. However, deciphering the structure of these assemblies is challenging owing to the difficulties in isolating complexes from the individual partners. The additive nature of small-angle X-ray scattering (SAXS) data allows for probing the species present in these mixtures, but decomposition into structural and thermodynamic information is difficult. We present a chemometric approach enabling the decomposition of titration SAXS data into species-specific information. Using extensive synthetic SAXS data, we demonstrate that robust decomposition can be achieved for titrations with a maximum fraction of complex of 0.5 that can be extended to 0.3 when two orthogonal titrations are simultaneously analyzed. The effect of the structural features, titration points, relative concentrations, and noise are thoroughly analyzed. The validation of the strategy with experimental data highlights the power of the approach to provide unique insights into this family of biomolecular assemblies.

Molecular and Structural Basis of Cross-Reactivity in M. tuberculosis Toxin-Antitoxin Systems

Mycobacterium tuberculosis genome encodes over 80 toxin-antitoxin (TA) systems. While each toxin interacts with its cognate antitoxin, the abundance of TA systems presents an opportunity for potential non-cognate interactions. TA systems mediate manifold interactions to manage pathogenicity and stress response network of the cell and non-cognate interactions may play vital roles as well. To address if non-cognate and heterologous interactions are feasible and to understand the structural basis of their interactions, we have performed comprehensive computational analyses on the available 3D structures and generated structural models of paralogous M. tuberculosis VapBC and MazEF TA systems. For a majority of the TA systems, we show that non-cognate toxin-antitoxin interactions are structurally incompatible except for complexes like VapBC15 and VapBC11, which show similar interfaces and potential for cross-reactivity. For TA systems which have been experimentally shown earlier to disfavor non-cognate interactions, we demonstrate that they are structurally and stereo-chemically incompatible. For selected TA systems, our detailed structural analysis identifies specificity conferring residues. Thus, our work improves the current understanding of TA interfaces and generates a hypothesis based on congenial binding site, geometric complementarity, and chemical nature of interfaces. Overall, our work offers a structure-based explanation for non-cognate toxin-antitoxin interactions in M. tuberculosis.

Soft regions of protein surface are potent for stable dimer formation

By having knowledge about the characteristics of protein interaction interfaces, we will be able to manipulate protein complexes for therapies. Dimer state is considered as the primary alphabet of the most proteins' quaternary structure. The properties of binding interface between subunits and of noninterface region define the specificity and stability of the intended protein complex. Considering some topological properties and amino acids' affinity for binding in interfaces of protein dimers, we construct the interface-specific recurrence plots. The data obtained from recurrence quantitative analysis, and accessibility-related metrics help us to classify the protein dimers into four distinct classes. Some mechanical properties of binding interfaces are computed for each predefined class of the dimers. The computed mechanical characteristics of binding patch region are compared with those of nonbinding region of proteins. Our observations indicate that the mechanical properties of protein binding sites have a decisive impact on determining the dimer stability. We introduce a new concept in analyzing protein structure by considering mechanical properties of protein structure. We conclude that the interface region between subunits of stable dimers is usually mechanically softer than the interface of unstable protein dimers. AbbreviationsAABaverage affinity for bindingANManisotropic network modelAPCaffinity propagation clusteringASAaccessible surface areaCCDinter residues distanceCSCcomplex stability codeDMdistance matrixΔG_dissPISA-computed dissociation free energyGNMGaussian normal mode analysisNMAnormal mode analysisPBPprotein binding patchPISAproteins, interfaces, structures and assembliesrASArelative accessible area in respect to unfolded state of residuesRMrecurrence matrixrPrelative protrusionRPrecurrence plotRQArecurrence quantitative analysisSEMstandard error of meanCommunicated by Ramaswamy H. Sarma.

Interface residues of transient protein-protein complexes have extensive intra-protein interactions apart from inter-protein interactions

Background

Protein-protein interactions are crucial for normal biological processes and to regulate cellular reactions that affect gene expression and function. Several previous studies have emphasized the roles of residues at the interface of protein-protein complexes in conferring stability and specificity to the complex. Interface residues in a protein are well known for their interactions with sidechain and main chain atoms with the interacting protein. However, the extent of intra-protein interactions involving interface residues in a protein-protein complex and their relative contribution in comparison to inter-protein interactions are not clearly understood. This paper probes this feature using a dataset of protein-protein complexes of known 3-D structure.

Results

We have analysed a dataset of 45 transient protein-protein complex structures with at least one of the interacting proteins with a known structure available also in the unbound form. We observe that a large proportion of interface residues (1608 out of 2137 interface residues, 75%) are involved in intra and inter-protein interactions simultaneously. The amino acid propensities of such interfacial residues involved in bifurcated interactions are found to be highly similar to the general propensities to occur at protein-protein interfaces. Finally, we observe that a majority (83%) of intra-protein interactions of interface residues with bifurcated interactions, are also observed in the protein uncomplexed form.

Conclusions

We have shown, to the best of our knowledge for the first time, that a vast majority of the protein-protein interface residues are involved in extensive intra-protein interactions apart from inter-protein interactions. For a majority of such interface residues the microenvironment in the tertiary structure is pre-formed and retained upon complex formation with its cognate partner during transient interactions.

Reviewers

This article was reviewed by Arumay Pal and Mallur Madhusudhan.

Electronic supplementary material

The online version of this article (10.1186/s13062-019-0232-2) contains supplementary material, which is available to authorized users.

Residue-level determinants of angiopoietin-2 interactions with its receptor Tie2

We combined computational and experimental methods to interrogate the binding determinants of angiopoietin-2 (Ang2) to its receptor tyrosine kinase (RTK) Tie2-a central signaling system in angiogenesis, inflammation, and tumorigenesis. We used physics-based electrostatic and surface-area calculations to identify the subset of interfacial Ang2 and Tie2 residues that can affect binding directly. Using random and site-directed mutagenesis and yeast surface display (YSD), we validated these predictions and identified additional Ang2 positions that affected receptor binding. We then used burial-based calculations to classify the larger set of Ang2 residues that are buried in the Ang2 core, whose mutations can perturb the Ang2 structure and thereby affect interactions with Tie2 indirectly. Our analysis showed that the Ang2-Tie2 interface is dominated by nonpolar contributions, with only three Ang2 and two Tie2 residues that contribute electrostatically to intermolecular interactions. Individual interfacial residues contributed only moderately to binding, suggesting that engineering of this interface will require multiple mutations to reach major effects. Conversely, substitutions in substantially buried Ang2 residues were more prevalent in our experimental screen, reduced binding substantially, and are therefore more likely to have a deleterious effect that might contribute to oncogenesis. Computational analysis of additional RTK-ligand complexes, c-Kit-SCF and M-CSF-c-FMS, and comparison to previous YSD results, further show the utility of our combined methodology.

Mathematical model identifies effective P53 accumulation with target gene binding affinity in DNA damage response for cell fate decision

Functional p53 signaling is essential for appropriate responses to diverse stimuli. P53 dynamics employs the information from the stimulus leading to selective gene expression and cell fate decision. However, the decoding mechanism of p53 dynamics under DNA damage challenge remains poorly understood. Here we mathematically modeled the recently dual-phase p53 dynamics under doxorubicin treatment. We found that p53 could perform sequential pulses followed by a high-amplitude terminal pulse at relatively low doxorubicin treatment, whereas p53 became steadily accumulated when damage level was high. The effective p53 integral above a threshold but not the absolute accumulation of p53 precisely discriminated survival and death. Silencing negative regulators in p53 network might promote the occurrence of terminal pulse. Furthermore, lower binding affinity and degradation rate of p53 target genes could favorably discriminate high and low dose doxorubicin treatment. Grouping by temporal profiles suggested that the p53 dynamics rather than the doxorubicin doses could better discriminate cellular outcomes and confer less variation for effective p53 integral reemphasizing the importance of p53 level regulation. Our model has established a theoretical framework that p53 dynamics can work cooperatively with its binding affinity to target genes leading to cell fate choice, providing new clues of optimized clinical design by manipulating p53 dynamics.

Modeling and Dynamic Behavior of eIF2 Dependent Regulatory System With Disturbances

Eukaryotic initiation factor 2 (eIF2) is a central controller of the eukaryotic translational machinery. To sustain the on-going translation activity, eIF2 cycles between its GTP and GDP bound states. However, in response to cellular stresses, the phosphorylation of eIF2 takes place, which acts as an inhibitor of the guanine nucleotide exchange factor eIF2B and switches the translation activity on physiological timescales. The main objective of this paper is to investigate the stability of the regulatory system under nominal conditions, parametric fluctuations, and structural damages. In this paper, a mathematical model of eIF2-dependent regulatory system is used to identify the stability-conferring features within the system with the help of direct and indirect methods of Lyapunov stability theory. To investigate the impact of intrinsic fluctuations and structural damages on the stability of regulatory system, the mathematical model has been linearized around feasible equilibrium point and the variation of system poles has been observed. The investigations have revealed that the regulatory model is stable and able to tolerate the intrinsic stressors but becomes unstable when particular complex is targeted to override the undesirable interaction. Our analyses indicate that, the stability is a collective property and damage in the structure of the system changes the stability of the system.

Specific DNA sequences allosterically enhance protein-protein interaction in a transcription factor through modulation of protein dynamics: implications for specificity of gene regulation

Most genes are regulated by multiple transcription factors, often assembling into multi-protein complexes in the gene regulatory region. Understanding of the molecular origin of specificity of gene regulatory complex formation in the context of the whole genome is currently inadequate. A phage transcription factor λ-CI forms repressive multi-protein complexes by binding to multiple binding sites in the genome to regulate the lifecycle of the phage. The protein-protein interaction between two DNA-bound λ-CI molecules is stronger when they are bound to the correct pair of binding sites, suggesting allosteric transmission of recognition of correct DNA sequences to the protein-protein interaction interface. Exploration of conformation and dynamics by time-resolved fluorescence anisotropy decay and molecular dynamics suggests a change in protein dynamics to be a crucial factor in mediating allostery. A lattice-based model suggests that DNA-sequence induced allosteric effects could be crucial underlying factors in differentially stabilizing the correct site-specific gene regulatory complexes. We conclude that transcription factors have evolved multiple mechanisms to augment the specificity of DNA-protein interactions in order to achieve an extraordinarily high degree of spatial and temporal specificities of gene regulatory complexes, and DNA-sequence induced allostery plays an important role in the formation of sequence-specific gene regulatory complexes.