An ultra‐high affinity protein–protein interface displaying sequence‐robustness

Protein–protein interactions are crucial in biology and play roles in for example, the immune system, signaling pathways, and enzyme regulation. Ultra‐high affinity interactions (K _d <0.1 nM) occur in these systems, however, structures and energetics behind stability of ultra‐high affinity protein–protein complexes are not well understood. Regulation of the starch debranching barley limit dextrinase (LD) and its endogenous cereal type inhibitor (LDI) exemplifies an ultra‐high affinity complex (K _d of 42 pM). In this study the LD–LDI complex is investigated to unveil how robust the ultra‐high affinity is to LDI sequence variation at the protein–protein interface and whether alternative sequences can retain the ultra‐high binding affinity. The interface of LD–LDI was engineered using computational protein redesign aiming at identifying LDI variants predicted to retain ultra‐high binding affinity. These variants present a very diverse set of mutations going beyond conservative and alanine substitutions typically used to probe interfaces. Surface plasmon resonance analysis of the LDI variants revealed that high affinity of LD–LDI requires interactions of several residues at the rim of the protein interface, unlike the classical hotspot arrangement where key residues are found at the center of the interface. Notably, substitution of interface residues in LDI, including amino acids with functional groups different from the wild‐type, could occur without loss of affinity. This demonstrates that ultra‐high binding affinity can be conferred without hotspot residues, thus making complexes more robust to mutational drift in evolution. The present mutational analysis also demonstrates how energetic coupling can emerge between residues at large distances at the interface.     

The free energy folding penalty accompanying binding of intrinsically disordered α‐helical motifs

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes and preform critical roles in many cellular processes, most often through the association with globular proteins. Despite lacking a stable three‐dimensional structure by themselves, they may acquire a defined conformation upon binding globular targets. The most common type of secondary structure acquired by these binding motifs entails formation of an α‐helix. It has been hypothesized that such disorder‐to‐order transitions are associated with a significant free energy penalty due to IDP folding, which reduces the overall IDP‐target affinity. However, the exact magnitude of IDP folding penalty in α‐helical binding motifs has not been systematically estimated. Here, we report the folding penalty contributions for 30 IDPs undergoing folding‐upon‐binding and find that the average IDP folding penalty is +2.0 kcal/mol and ranges from 0.7 to 3.5 kcal/mol. We observe that the folding penalty scales approximately linearly with the change in IDP helicity upon binding, which provides a simple empirical way to estimate folding penalty. We analyze to what extent do pre‐structuring and target‐bound IDP dynamics (fuzziness) reduce the folding penalty and find that these effects combined, on average, reduce the folding cost by around half. Taken together, the presented analysis provides a quantitative basis for understanding the role of folding penalty in IDP‐target interactions and introduces a method estimate this quantity. Estimation and reduction of IDP folding penalty may prove useful in the rational design of helix‐stabilized inhibitors of IDP‐target interactions.StatementThe α‐helical binding motifs are ubiquitous among the intrinsically disordered proteins (IDPs). Upon binding their targets, they undergo a disorder‐to‐order transition, which is accompanied by a significant folding penalty whose magnitude is generally not known. Here, we use recently developed statistical‐thermodynamic model to estimate the folding penalties for 30 IDPs and clarify the roles of IDP pre‐folding and bound‐state dynamics in reducing the folding penalty.     

Single‐molecule imaging reveals Tau trapping at nanometer‐sized dynamic hot spots near the plasma membrane that persists after microtubule perturbation and cholesterol depletion

Accumulation of aggregates of the microtubule‐binding protein Tau is a pathological hallmark of Alzheimer''s disease. While Tau is thought to primarily associate with microtubules, it also interacts with and localizes to the plasma membrane. However, little is known about how Tau behaves and organizes at the plasma membrane of live cells. Using quantitative, single‐molecule imaging, we show that Tau exhibits spatial and kinetic heterogeneity near the plasma membrane of live cells, resulting in the formation of nanometer‐sized hot spots. The hot spots lasted tens of seconds, much longer than the short dwell time (∼ 40 ms) of Tau on microtubules. Pharmacological and biochemical disruption of Tau/microtubule interactions did not prevent hot spot formation, suggesting that these are different from the reported Tau condensation on microtubules. Although cholesterol removal has been shown to reduce Tau pathology, its acute depletion did not affect Tau hot spot dynamics. Our study identifies an intrinsic dynamic property of Tau near the plasma membrane that may facilitate the formation of assembly sites for Tau to assume its physiological and pathological functions.     

Binding of small molecules at interface of protein–protein complex – A newer approach to rational drug design

Protein–protein interaction is a vital process which drives many important physiological processes in the cell and has also been implicated in several diseases. Though the protein–protein interaction network is quite complex but understanding its interacting partners using both in silico as well as molecular biology techniques can provide better insights for targeting such interactions. Targeting protein–protein interaction with small molecules is a challenging task because of druggability issues. Nevertheless, several studies on the kinetics as well as thermodynamic properties of protein–protein interactions have immensely contributed toward better understanding of the affinity of these complexes. But, more recent studies on hot spots and interface residues have opened up new avenues in the drug discovery process. This approach has been used in the design of hot spot based modulators targeting protein–protein interaction with the objective of normalizing such interactions.     

ATP regulates RNA‐driven cold inducible RNA binding protein phase separation

Intrinsically disordered proteins and proteins containing intrinsically disordered regions are highly abundant in the proteome of eukaryotes and are extensively involved in essential biological functions. More recently, their role in the organization of biomolecular condensates has become evident and along with their misregulation in several neurologic disorders. Currently, most studies involving these proteins are carried out in vitro and using purified proteins. Given that in cells, condensate‐forming proteins are exposed to high, millimolar concentrations of cellular metabolites, we aimed to reveal the interactions of cellular metabolites and a representative condensate‐forming protein. Here, using the arginine–glycine/arginine–glycine–glycine (RG/RGG)‐rich cold inducible RNA binding protein (CIRBP) as paradigm, we studied binding of the cellular metabolome to CIRBP. We found that most of the highly abundant cellular metabolites, except nucleotides, do not directly bind to CIRBP. ATP, ADP, and AMP as well as NAD⁺, NADH, NADP⁺, and NADPH directly interact with CIRBP, involving both the folded RNA‐recognition motif and the disordered RG/RGG region. ATP binding inhibited RNA‐driven phase separation of CIRBP. Thus, it might be beneficial to include cellular metabolites in in vitro liquid–liquid phase separation studies of RG/RGG and other condensate‐forming proteins in order to better mimic the cellular environment in the future.     

Small-world network approach to identify key residues in protein-protein interaction

We show that protein complexes can be represented as small-world networks, exhibiting a relatively small number of highly central amino-acid residues occurring frequently at protein-protein interfaces. We further base our analysis on a set of different biological examples of protein-protein interactions with experimentally validated hot spots, and show that 83% of these predicted highly central residues, which are conserved in sequence alignments and nonexposed to the solvent in the protein complex, correspond to or are in direct contact with an experimentally annotated hot spot. The remaining 17% show a general tendency to be close to an annotated hot spot. On the other hand, although there is no available experimental information on their contribution to the binding free energy, detailed analysis of their properties shows that they are good candidates for being hot spots. Thus, highly central residues have a clear tendency to be located in regions that include hot spots. We also show that some of the central residues in the protein complex interfaces are central in the monomeric structures before dimerization and that possible information relating to hot spots of binding free energy could be obtained from the unbound structures.     

Disease‐linked TDP‐43 hyperphosphorylation suppresses TDP‐43 condensation and aggregation

Post‐translational modifications (PTMs) have emerged as key modulators of protein phase separation and have been linked to protein aggregation in neurodegenerative disorders. The major aggregating protein in amyotrophic lateral sclerosis and frontotemporal dementia, the RNA‐binding protein TAR DNA‐binding protein (TDP‐43), is hyperphosphorylated in disease on several C‐terminal serine residues, a process generally believed to promote TDP‐43 aggregation. Here, we however find that Casein kinase 1δ‐mediated TDP‐43 hyperphosphorylation or C‐terminal phosphomimetic mutations reduce TDP‐43 phase separation and aggregation, and instead render TDP‐43 condensates more liquid‐like and dynamic. Multi‐scale molecular dynamics simulations reveal reduced homotypic interactions of TDP‐43 low‐complexity domains through enhanced solvation of phosphomimetic residues. Cellular experiments show that phosphomimetic substitutions do not affect nuclear import or RNA regulatory functions of TDP‐43, but suppress accumulation of TDP‐43 in membrane‐less organelles and promote its solubility in neurons. We speculate that TDP‐43 hyperphosphorylation may be a protective cellular response to counteract TDP‐43 aggregation.     

SARS‐CoV‐2 nucleocapsid protein phase‐separates with RNA and with human hnRNPs

Tightly packed complexes of nucleocapsid protein and genomic RNA form the core of viruses and assemble within viral factories, dynamic compartments formed within the host cells associated with human stress granules. Here, we test the possibility that the multivalent RNA‐binding nucleocapsid protein (N) from severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) condenses with RNA via liquid–liquid phase separation (LLPS) and that N protein can be recruited in phase‐separated forms of human RNA‐binding proteins associated with SG formation. Robust LLPS with RNA requires two intrinsically disordered regions (IDRs), the N‐terminal IDR and central‐linker IDR, as well as the folded C‐terminal oligomerization domain, while the folded N‐terminal domain and the C‐terminal IDR are not required. N protein phase separation is induced by addition of non‐specific RNA. In addition, N partitions in vitro into phase‐separated forms of full‐length human hnRNPs (TDP‐43, FUS, hnRNPA2) and their low‐complexity domains (LCs). These results provide a potential mechanism for the role of N in SARS‐CoV‐2 viral genome packing and in host‐protein co‐opting necessary for viral replication and infectivity.     

Using collections of structural models to predict changes of binding affinity caused by mutations in protein–protein interactions

Protein–protein interactions (PPIs) in all the molecular aspects that take place both inside and outside cells. However, determining experimentally the structure and affinity of PPIs is expensive and time consuming. Therefore, the development of computational tools, as a complement to experimental methods, is fundamental. Here, we present a computational suite: MODPIN, to model and predict the changes of binding affinity of PPIs. In this approach we use homology modeling to derive the structures of PPIs and score them using state‐of‐the‐art scoring functions. We explore the conformational space of PPIs by generating not a single structural model but a collection of structural models with different conformations based on several templates. We apply the approach to predict the changes in free energy upon mutations and splicing variants of large datasets of PPIs to statistically quantify the quality and accuracy of the predictions. As an example, we use MODPIN to study the effect of mutations in the interaction between colicin endonuclease 9 and colicin endonuclease 2 immune protein from Escherichia coli. Finally, we have compared our results with other state‐of‐art methods.     

Systematic discovery of linear binding motifs targeting an ancient protein interaction surface on MAP kinases

Mitogen‐activated protein kinases (MAPK) are broadly used regulators of cellular signaling. However, how these enzymes can be involved in such a broad spectrum of physiological functions is not understood. Systematic discovery of MAPK networks both experimentally and in silico has been hindered because MAPKs bind to other proteins with low affinity and mostly in less‐characterized disordered regions. We used a structurally consistent model on kinase‐docking motif interactions to facilitate the discovery of short functional sites in the structurally flexible and functionally under‐explored part of the human proteome and applied experimental tools specifically tailored to detect low‐affinity protein–protein interactions for their validation in vitro and in cell‐based assays. The combined computational and experimental approach enabled the identification of many novel MAPK‐docking motifs that were elusive for other large‐scale protein–protein interaction screens. The analysis produced an extensive list of independently evolved linear binding motifs from a functionally diverse set of proteins. These all target, with characteristic binding specificity, an ancient protein interaction surface on evolutionarily related but physiologically clearly distinct three MAPKs (JNK, ERK, and p38). This inventory of human protein kinase binding sites was compared with that of other organisms to examine how kinase‐mediated partnerships evolved over time. The analysis suggests that most human MAPK‐binding motifs are surprisingly new evolutionarily inventions and newly found links highlight (previously hidden) roles of MAPKs. We propose that short MAPK‐binding stretches are created in disordered protein segments through a variety of ways and they represent a major resource for ancient signaling enzymes to acquire new regulatory roles.     

The mechanism of action of the SSB interactome reveals it is the first OB‐fold family of genome guardians in prokaryotes

The single‐stranded DNA binding protein (SSB) is essential to all aspects of DNA metabolism in bacteria. This protein performs two distinct, but closely intertwined and indispensable functions in the cell. SSB binds to single‐stranded DNA (ssDNA) and at least 20 partner proteins resulting in their regulation. These partners comprise a family of genome guardians known as the SSB interactome. Essential to interactome regulation is the linker/OB‐fold network of interactions. This network of interactions forms when one or more PXXP motifs in the linker of SSB bind to an OB‐fold in a partner, with interactome members involved in competitive binding between the linker and ssDNA to their OB‐fold. Consequently, when linker‐binding occurs to an OB‐fold in an interactome partner, proteins are loaded onto the DNA. When linker/OB‐fold interactions occur between SSB tetramers, cooperative ssDNA‐binding results, producing a multi‐tetrameric complex that rapidly protects the ssDNA. Within this SSB‐ssDNA complex, there is an extensive and dynamic network of linker/OB‐fold interactions that involves multiple tetramers bound contiguously along the ssDNA lattice. The dynamic behavior of these tetramers which includes binding mode changes, sliding as well as DNA wrapping/unwrapping events, are likely coupled to the formation and disruption of linker/OB‐fold interactions. This behavior is essential to facilitating downstream DNA processing events. As OB‐folds are critical to the essence of the linker/OB‐fold network of interactions, and they are found in multiple interactome partners, the SSB interactome is classified as the first family of prokaryotic, oligosaccharide/oligonucleotide binding fold (OB‐fold) genome guardians.     

Cryo‐EM structures of pentameric autoinducer‐2 exporter from Escherichia coli reveal its transport mechanism

Bacteria utilize small extracellular molecules to communicate in order to collectively coordinate their behaviors in response to the population density. Autoinducer‐2 (AI‐2), a universal molecule for both intra‐ and inter‐species communication, is involved in the regulation of biofilm formation, virulence, motility, chemotaxis, and antibiotic resistance. While many studies have been devoted to understanding the biosynthesis and sensing of AI‐2, very little information is available on its export. The protein TqsA from Escherichia coli, which belongs to the AI‐2 exporter superfamily, has been shown to export AI‐2. Here, we report the cryogenic electron microscopic structures of two AI‐2 exporters (TqsA and YdiK) from E. coli at 3.35 Å and 2.80 Å resolutions, respectively. Our structures suggest that the AI‐2 exporter exists as a homo‐pentameric complex. In silico molecular docking and native mass spectrometry experiments were employed to demonstrate the interaction between AI‐2 and TqsA, and the results highlight the functional importance of two helical hairpins in substrate binding. We propose that each monomer works as an independent functional unit utilizing an elevator‐type transport mechanism.     

Thermodynamic coupling between neighboring binding sites in homo‐oligomeric ligand sensing proteins from mass resolved ligand‐dependent population distributions

Homo‐oligomeric ligand‐activated proteins are ubiquitous in biology. The functions of such molecules are commonly regulated by allosteric coupling between ligand‐binding sites. Understanding the basis for this regulation requires both quantifying the free energy ΔG transduced between sites, and the structural basis by which it is transduced. We consider allostery in three variants of the model ring‐shaped homo‐oligomeric trp RNA‐binding attenuation protein (TRAP). First, we developed a nearest‐neighbor statistical thermodynamic binding model comprising microscopic free energies for ligand binding to isolated sites ΔG ₀, and for coupling between adjacent sites, ΔG _α . Using the resulting partition function (PF) we explored the effects of these parameters on simulated population distributions for the 2 ^N possible liganded states. We then experimentally monitored ligand‐dependent population shifts using conventional spectroscopic and calorimetric methods and using native mass spectrometry (MS). By resolving species with differing numbers of bound ligands by their mass, native MS revealed striking differences in their ligand‐dependent population shifts. Fitting the populations to a binding polynomial derived from the PF yielded coupling free energy terms corresponding to orders of magnitude differences in cooperativity. Uniquely, this approach predicts which of the possible 2 ^N liganded states are populated at different ligand concentrations, providing necessary insights into regulation. The combination of statistical thermodynamic modeling with native MS may provide the thermodynamic foundation for a meaningful understanding of the structure–thermodynamic linkage that drives cooperativity.     

A human kinase yeast array for the identification of kinases modulating phosphorylation‐dependent protein–protein interactions

Protein kinases play an important role in cellular signaling pathways and their dysregulation leads to multiple diseases, making kinases prime drug targets. While more than 500 human protein kinases are known to collectively mediate phosphorylation of over 290,000 S/T/Y sites, the activities have been characterized only for a minor, intensively studied subset. To systematically address this discrepancy, we developed a human kinase array in Saccharomyces cerevisiae as a simple readout tool to systematically assess kinase activities. For this array, we expressed 266 human kinases in four different S. cerevisiae strains and profiled ectopic growth as a proxy for kinase activity across 33 conditions. More than half of the kinases showed an activity‐dependent phenotype across many conditions and in more than one strain. We then employed the kinase array to identify the kinase(s) that can modulate protein–protein interactions (PPIs). Two characterized, phosphorylation‐dependent PPIs with unknown kinase–substrate relationships were analyzed in a phospho‐yeast two‐hybrid assay. CK2α1 and SGK2 kinases can abrogate the interaction between the spliceosomal proteins AAR2 and PRPF8, and NEK6 kinase was found to mediate the estrogen receptor (ERα) interaction with 14‐3‐3 proteins. The human kinase yeast array can thus be used for a variety of kinase activity‐dependent readouts.     

HIV‐1 infection of CD4 T cells impairs antigen‐specific B cell function

Failures to produce neutralizing antibodies upon HIV‐1 infection result in part from B‐cell dysfunction due to unspecific B‐cell activation. How HIV‐1 affects antigen‐specific B‐cell functions remains elusive. Using an adoptive transfer mouse model and ex vivo HIV infection of human tonsil tissue, we found that expression of the HIV‐1 pathogenesis factor NEF in CD4 T cells undermines their helper function and impairs cognate B‐cell functions including mounting of efficient specific IgG responses. NEF interfered with T cell help via a specific protein interaction motif that prevents polarized cytokine secretion at the T‐cell–B‐cell immune synapse. This interference reduced B‐cell activation and proliferation and thus disrupted germinal center formation and affinity maturation. These results identify NEF as a key component for HIV‐mediated dysfunction of antigen‐specific B cells. Therapeutic targeting of the identified molecular surface in NEF will facilitate host control of HIV infection.     

Artificial intelligence based methods for hot spot prediction

Proteins interact through their interfaces to fulfill essential functions in the cell. They bind to their partners in a highly specific manner and form complexes that have a profound effect on understanding the biological pathways they are involved in. Any abnormal interactions may cause diseases. Therefore, the identification of small molecules which modulate protein interactions through their interfaces has high therapeutic potential. However, discovering such molecules is challenging. Most protein–protein binding affinity is attributed to a small set of amino acids found in protein interfaces known as hot spots. Recent studies demonstrate that drug-like small molecules specifically may bind to hot spots. Therefore, hot spot prediction is crucial. As experimental data accumulates, artificial intelligence begins to be used for computational hot spot prediction. First, we review machine learning and deep learning for computational hot spot prediction and then explain the significance of hot spots toward drug design.     

The leucine zipper EF‐hand containing transmembrane protein‐1 EF‐hand is a tripartite calcium,temperature, and pH sensor

Leucine Zipper EF‐hand containing transmembrane protein‐1 (LETM1) is an inner mitochondrial membrane protein that mediates mitochondrial calcium (Ca²⁺)/proton exchange. The matrix residing carboxyl (C)‐terminal domain contains a sequence identifiable EF‐hand motif (EF1) that is highly conserved among orthologues. Deletion of EF1 abrogates LETM1 mediated mitochondrial Ca²⁺ flux, highlighting the requirement of EF1 for LETM1 function. To understand the mechanistic role of this EF‐hand in LETM1 function, we characterized the biophysical properties of EF1 in isolation. Our data show that EF1 exhibits α‐helical secondary structure that is augmented in the presence of Ca²⁺. Unexpectedly, EF1 features a weak (~mM), but specific, apparent Ca²⁺‐binding affinity, consistent with the canonical Ca²⁺ coordination geometry, suggested by our solution NMR. The low affinity is, at least in part, due to an Asp at position 12 of the binding loop, where mutation to Glu increases the affinity by ~4‐fold. Further, the binding affinity is sensitive to pH changes within the physiological range experienced by mitochondria. Remarkably, EF1 unfolds at high and low temperatures. Despite these unique EF‐hand properties, Ca²⁺ binding increases the exposure of hydrophobic regions, typical of EF‐hands; however, this Ca²⁺‐induced conformational change shifts EF1 from a monomer to higher order oligomers. Finally, we showed that a second, putative EF‐hand within LETM1 is unreactive to Ca²⁺ either in isolation or tandem with EF1. Collectively, our data reveal that EF1 is structurally and biophysically responsive to pH, Ca²⁺ and temperature, suggesting a role as a multipartite environmental sensor within LETM1.     

A weakened interface in the P182L variant of HSP27 associated with severe Charcot‐Marie‐Tooth neuropathy causes aberrant binding to interacting proteins

HSP27 is a human molecular chaperone that forms large, dynamic oligomers and functions in many aspects of cellular homeostasis. Mutations in HSP27 cause Charcot‐Marie‐Tooth (CMT) disease, the most common inherited disorder of the peripheral nervous system. A particularly severe form of CMT disease is triggered by the P182L mutation in the highly conserved IxI/V motif of the disordered C‐terminal region, which interacts weakly with the structured core domain of HSP27. Here, we observed that the P182L mutation disrupts the chaperone activity and significantly increases the size of HSP27 oligomers formed in vivo, including in motor neurons differentiated from CMT patient‐derived stem cells. Using NMR spectroscopy, we determined that the P182L mutation decreases the affinity of the HSP27 IxI/V motif for its own core domain, leaving this binding site more accessible for other IxI/V‐containing proteins. We identified multiple IxI/V‐bearing proteins that bind with higher affinity to the P182L variant due to the increased availability of the IxI/V‐binding site. Our results provide a mechanistic basis for the impact of the P182L mutation on HSP27 and suggest that the IxI/V motif plays an important, regulatory role in modulating protein–protein interactions.