A comparison of successful and failed protein interface designs highlights the challenges of designing buried hydrogen bonds

The accurate design of new protein–protein interactions is a longstanding goal of computational protein design. However, most computationally designed interfaces fail to form experimentally. This investigation compares five previously described successful de novo interface designs with 158 failures. Both sets of proteins were designed with the molecular modeling program Rosetta. Designs were considered a success if a high‐resolution crystal structure of the complex closely matched the design model and the equilibrium dissociation constant for binding was less than 10 μM. The successes and failures represent a wide variety of interface types and design goals including heterodimers, homodimers, peptide‐protein interactions, one‐sided designs (i.e., where only one of the proteins was mutated) and two‐sided designs. The most striking feature of the successful designs is that they have fewer polar atoms at their interfaces than many of the failed designs. Designs that attempted to create extensive sets of interface‐spanning hydrogen bonds resulted in no detectable binding. In contrast, polar atoms make up more than 40% of the interface area of many natural dimers, and native interfaces often contain extensive hydrogen bonding networks. These results suggest that Rosetta may not be accurately balancing hydrogen bonding and electrostatic energies against desolvation penalties and that design processes may not include sufficient sampling to identify side chains in preordered conformations that can fully satisfy the hydrogen bonding potential of the interface.     

Thermodynamic analysis of subunit interactions in Escherichia coli molybdopterin synthase

The molybdopterin (MPT) synthase complex in Escherichia coli consists of two MoaE subunits and two MoaD subunits in a heterotetrameric structure with the two MoaE subunits forming a central dimer. Each MoaD subunit binds to a single MoaE molecule to form two identical MoaE/MoaD interfaces. Here we define the thermodynamic properties of the interaction between MoaE and MoaD in MPT synthase using a H/D exchange and matrix-assisted laser desorption/ionization (MALDI) mass spectroscopy based method termed SUPREX (stability of unpurified proteins from rates of H/D exchange). SUPREX-derived protein folding free energies and m values are reported for MoaE in the presence and absence of MoaD and MoaD-SH, the thiocarboxylated form of MoaD that is essential for the catalytic activity of MPT synthase. The protein folding free energy measurements were used to calculate a dissociation constant of 17 +/- 7 microM for the binding of MoaD to MoaE in inactive MPT synthase and a dissociation constant of 2.6 +/- 0.9 microM for the binding of MoaD-SH to MoaE in active MPT synthase. The increased binding affinity of MoaD-SH for MoaE is consistent with a previously proposed mechanism for the MPT synthase reaction. Using the increased m values exhibited by MoaE in the presence of either MoaD subunit, the solvent accessible surface area buried upon formation of the subunit interface in MPT synthase was estimated to be 2378 A(2) for inactive MPT synthase and 4117 A(2) for active MPT synthase.     

Dynamic features of homodimer interfaces calculated by normal‐mode analysis

Knowledge of the dynamic features of protein interfaces is necessary for a deeper understanding of protein–protein interactions. We performed normal‐mode analysis (NMA) of 517 nonredundant homodimers and their protomers to characterize dimer interfaces from a dynamic perspective. The motion vector calculated by NMA for each atom of a dimer was decomposed into internal and external motion vectors in individual component subunits, followed by the averaging of time‐averaged correlations between these vectors over atom pairs in the interface. This averaged correlation coefficient (ACC) was defined for various combinations of vectors and investigated in detail. ACCs decrease exponentially with an increasing interface area and r‐value, that is, interface area divided by the entire subunit surface area. As the r‐value reflects the nature of dimer formation, the result suggests that both the interface area and the nature of dimer formation are responsible for the dynamic properties of dimer interfaces. For interfaces with small or medium r‐values and without intersubunit entanglements, ACCs are found to increase on dimer formation when compared with those in the protomer state. In contrast, ACCs do not increase on dimer formation for interfaces with large r‐values and intersubunit entanglements such as in interwinding dimers. Furthermore, relationships between ACCs for intrasubunit atom pairs and for intersubunit atom pairs are found to significantly differ between interwinding and noninterwinding dimers for external motions. External motions are considered as an important factor for characterizing dimer interfaces.     

Exploring functional roles of multibinding protein interfaces

Cellular processes are highly interconnected and many proteins are shared in different pathways. Some of these shared proteins or protein families may interact with diverse partners using the same interface regions; such multibinding proteins are the subject of our study. The main goal of our study is to attempt to decipher the mechanisms of specific molecular recognition of multiple diverse partners by promiscuous protein regions. To address this, we attempt to analyze the physicochemical properties of multibinding interfaces and highlight the major mechanisms of functional switches realized through multibinding. We find that only 5% of protein families in the structure database have multibinding interfaces, and multibinding interfaces do not show any higher sequence conservation compared with the background interface sites. We highlight several important functional mechanisms utilized by multibinding families. (a) Overlap between different functional pathways can be prevented by the switches involving nearby residues of the same interfacial region. (b) Interfaces can be reused in pathways where the substrate should be passed from one protein to another sequentially. (c) The same protein family can develop different specificities toward different binding partners reusing the same interface; and finally, (d) inhibitors can attach to substrate binding sites as substrate mimicry and thereby prevent substrate binding.     

Protein-protein interaction and quaternary structure

Protein-protein recognition plays an essential role in structure and function. Specific non-covalent interactions stabilize the structure of macromolecular assemblies, exemplified in this review by oligomeric proteins and the capsids of icosahedral viruses. They also allow proteins to form complexes that have a very wide range of stability and lifetimes and are involved in all cellular processes. We present some of the structure-based computational methods that have been developed to characterize the quaternary structure of oligomeric proteins and other molecular assemblies and analyze the properties of the interfaces between the subunits. We compare the size, the chemical and amino acid compositions and the atomic packing of the subunit interfaces of protein-protein complexes, oligomeric proteins, viral capsids and protein-nucleic acid complexes. These biologically significant interfaces are generally close-packed, whereas the non-specific interfaces between molecules in protein crystals are loosely packed, an observation that gives a structural basis to specific recognition. A distinction is made within each interface between a core that contains buried atoms and a solvent accessible rim. The core and the rim differ in their amino acid composition and their conservation in evolution, and the distinction helps correlating the structural data with the results of site-directed mutagenesis and in vitro studies of self-assembly.     

Formation and carbon monoxide‐dependent dissociation of Allochromatium vinosum cytochrome c′ oligomers using domain‐swapped dimers

The number of artificial protein supramolecules has been increasing; however, control of protein oligomer formation remains challenging. Cytochrome c′ from Allochromatium vinosum (AVCP) is a homodimeric protein in its native form, where its protomer exhibits a four‐helix bundle structure containing a covalently bound five‐coordinate heme as a gas binding site. AVCP exhibits a unique reversible dimer–monomer transition according to the absence and presence of CO. Herein, domain‐swapped dimeric AVCP was constructed and utilized to form a tetramer and high‐order oligomers. The X‐ray crystal structure of oxidized tetrameric AVCP consisted of two monomer subunits and one domain‐swapped dimer subunit, which exchanged the region containing helices αA and αB between protomers. The active site structures of the domain‐swapped dimer subunit and monomer subunits in the tetramer were similar to those of the monomer subunits in the native dimer. The subunit–subunit interactions at the interfaces of the domain‐swapped dimer and monomer subunits in the tetramer were also similar to the subunit–subunit interaction in the native dimer. Reduced tetrameric AVCP dissociated to a domain‐swapped dimer and two monomers upon CO binding. Without monomers, the domain‐swapped dimers formed tetramers, hexamers, and higher‐order oligomers in the absence of CO, whereas the oligomers dissociated to domain‐swapped dimers in the presence of CO, demonstrating that the domain‐swapped dimer maintains the CO‐induced subunit dissociation behavior of native ACVP. These results suggest that protein oligomer formation may be controlled by utilizing domain swapping for a dimer–monomer transition protein.     

Binding linkage in a telomere DNA-protein complex at the ends of Oxytricha nova chromosomes

Alpha and beta protein subunits of the telomere end binding protein from Oxytricha nova (OnTEBP) combine with telomere single strand DNA to form a protective cap at the ends of chromosomes. We tested how protein-protein interactions seen in the co-crystal structure relate to DNA binding through use of fusion proteins engineered as different combinations of domains and subunits derived from OnTEBP. Joining alpha and beta resulted in a protein that bound single strand telomere DNA with high affinity (K(D-DNA)=1.4 nM). Another fusion protein, constructed without the C-terminal protein-protein interaction domain of alpha, bound DNA with 200-fold diminished affinity (K(D-DNA)=290 nM) even though the DNA-binding domains of alpha and beta were joined through a peptide linker. Adding back the alpha C-terminal domain as a separate protein restored high-affinity DNA binding. The binding behaviors of these fusion proteins and the native protein subunits are consistent with cooperative linkage between protein-association and DNA-binding equilibria. Linking DNA-protein stability to protein-protein contacts at a remote site may provide a trigger point for DNA-protein disassembly during telomere replication when the single strand telomere DNA must exchange between a very stable OnTEBP complex and telomerase.     

Plant Kinesin-Like Calmodulin Binding Protein Employs Its Regulatory Domain for Dimerization

Kinesin-like calmodulin binding protein (KCBP), a Kinesin-14 family motor protein, is involved in the structural organization of microtubules during mitosis and trichome morphogenesis in plants. The molecular mechanism of microtubule bundling by KCBP remains unknown. KCBP binding to microtubules is regulated by Ca²⁺-binding proteins that recognize its C-terminal regulatory domain. In this work, we have discovered a new function of the regulatory domain. We present a crystal structure of an Arabidopsis KCBP fragment showing that the C-terminal regulatory domain forms a dimerization interface for KCBP. This dimerization site is distinct from the dimerization interface within the N-terminal domain. Side chains of hydrophobic residues of the calmodulin binding helix of the regulatory domain form the C-terminal dimerization interface. Biochemical experiments show that another segment of the regulatory domain located beyond the dimerization interface, its negatively charged coil, is unexpectedly and absolutely required to stabilize the dimers. The strong microtubule bundling properties of KCBP are unaffected by deletion of the C-terminal regulatory domain. The slow minus-end directed motility of KCBP is also unchanged in vitro. Although the C-terminal domain is not essential for microtubule bundling, we suggest that KCBP may use its two independent dimerization interfaces to support different types of bundled microtubule structures in cells. Two distinct dimerization sites may provide a mechanism for microtubule rearrangement in response to Ca²⁺ signaling since Ca²⁺- binding proteins can disengage KCBP dimers dependent on its C-terminal dimerization interface.     

Both ends of Escherichia coli ribosomal protein S13 are immunochemically accessible in situ.

To investigate the structure ofEscherichia coli ribosomal protein S13 in 30S ribosomal subunits, we have previously generated 22 S13 specific monoclonal antibodies and mapped their specific epitopes to the S13 sequence. The availability of these S13 epitopesin situ has been further examined by incubating these monoclonal antibodies with 30S ribosomal subunits and analyzing formation of monoclonal antibody-linked ribosome dimers by sucrose gradients centrifugation. We have found that none of the 22 monoclonal antibodies makes ribosome dimers individually as do typical antisera. However, one monoclonal antibody, designated AS13-MAb 2, reacts with 30S ribosomal subunits to form immunocomplexes sedimenting faster than subunit monomers. When AS13-MAb 2 is paired with any one of three monoclonal antibodies directed to the S13 C-terminal epitopes, dimer formation is observed. Other pairs of monoclonal antibodies directed to distinct S13 epitopes have been tested similarly for dimer formation. Monoclonal antibody AS13-MAb 22, directed to the N-terminal region of 22 residues, also causes subunits to form typical dimers, but only if paired with one of the three monoclonal antibodies directed to the S13 C-terminal region. The close proximity of the epitopes recognized by AS13-MAbs 2 and 22 has been established by the mutual competition between the antibodies binding to intact 30S subunits. These results corroborate our previous observation, using polyclonal antibodies, that S13 has more than one epitope exposed on 30S subunits. Our finding that sequences on both ends of the S13 molecule are immunochemically accessible provides information about the molecular organization of S13in situ.     

Protein-protein interaction regulates proteins' mechanical stability

Elastomeric proteins are molecular springs found not only in a variety of biological machines and tissues, but also in biomaterials of superb mechanical properties. Regulating the mechanical stability of elastomeric proteins is not only important for a range of biological processes, but also critical for the use of engineered elastomeric proteins as building blocks to construct nanomechanical devices and novel materials of well-defined mechanical properties. Here we demonstrate that protein-protein interactions can potentially serve as an effective means to regulate the mechanical properties of elastomeric proteins. We show that the binding of fragments of IgG antibody to a small protein, GB1, can significantly enhance the mechanical stability of GB1. The regulation of the mechanical stability of GB1 by IgG fragments is not through direct modification of the interactions in the mechanically key region of GB1; instead, it is accomplished via the long-range coupling between the IgG binding site and the mechanically key region of GB1. Although Fc and Fab bind GB1 at different regions of GB1, their binding to GB1 can increase the mechanical stability of GB1 significantly. Using alanine point mutants of GB1, we show that the amplitude of mechanical stability enhancement of GB1 by Fc does not correlate with the binding affinity, suggesting that binding affinity only affects the population of GB1/human Fc (hFc) complex at a given concentration of hFc, but does not affect the intrinsic mechanical stability of the GB1/hFc complex. Furthermore, our results indicate that the mechanical stability enhancement by IgG fragments is robust and can tolerate sequence/structural perturbation to GB1. Our results demonstrate that the protein-protein interaction is an efficient approach to regulate the mechanical stability of GB1-like proteins and we anticipate that this new methodology will help to develop novel elastomeric proteins with tunable mechanical stability and compliance.     

An atypical heterotrimeric Gα protein has substantially reduced nucleotide binding but retains nucleotide-independent interactions with its cognate RGS protein and Gβγ dimer

Abstract

Plants uniquely have a family of proteins called extra-large G proteins (XLG) that share homology in their C-terminal half with the canonical Gα subunits; we carefully detail here that Arabidopsis XLG2 lacks critical residues requisite for nucleotide binding and hydrolysis which is consistent with our quantitative analyses. Based on microscale thermophoresis, Arabidopsis XLG2 binds GTPγS with an affinity 100 times lower than that to canonical Gα subunits. This means that given the concentration range of guanine nucleotide in plant cells, XLG2 is not likely bound by GTP in vivo. Homology modeling and molecular dynamics simulations provide a plausible mechanism for the poor nucleotide binding affinity of XLG2. Simulations indicate substantially stronger salt bridge networks formed by several key amino-acid residues of AtGPA1 which are either misplaced or missing in XLG2. These residues in AtGPA1 not only maintain the overall shape and integrity of the apoprotein cavity but also increase the frequency of favorable nucleotide-protein interactions in the nucleotide-bound state. Despite this loss of nucleotide dependency, XLG2 binds the RGS domain of AtRGS1 with an affinity similar to the Arabidopsis AtGPA1 in its apo-state and about 2 times lower than AtGPA1 in its transition state. In addition, XLG2 binds the Gβγ dimer with an affinity similar to that of AtGPA1. XLG2 likely acts as a dominant negative Gα protein to block G protein signaling. We propose that XLG2, independent of guanine nucleotide binding, regulates the active state of the canonical G protein pathway directly by sequestering Gβγ and indirectly by promoting heterodimer formation.

Communicated by Ramaswamy H. Sarma     

The Subunit Interfaces of Weakly Associated Homodimeric Proteins

We analyzed subunit interfaces in 315 homodimers with an X-ray structure in the Protein Data Bank, validated by checking the literature for data that indicate that the proteins are dimeric in solution and that, in the case of the “weak” dimers, the homodimer is in equilibrium with the monomer. The interfaces of the 42 weak dimers, which are smaller by a factor of 2.4 on average than in the remainder of the set, are comparable in size with antibody-antigen or protease-inhibitor interfaces. Nevertheless, they are more hydrophobic than in the average transient protein-protein complex and similar in amino acid composition to the other homodimer interfaces. The mean numbers of interface hydrogen bonds and hydration water molecules per unit area are also similar in homodimers and transient complexes. Parameters related to the atomic packing suggest that many of the weak dimer interfaces are loosely packed, and we suggest that this contributes to their low stability. To evaluate the evolutionary selection pressure on interface residues, we calculated the Shannon entropy of homologous amino acid sequences at 60% sequence identity. In 93% of the homodimers, the interface residues are better conserved than the residues on the protein surface. The weak dimers display the same high degree of interface conservation as other homodimers, but their homologs may be heterodimers as well as homodimers. Their interfaces may be good models in terms of their size, composition, and evolutionary conservation for the labile subunit contacts that allow protein assemblies to share and exchange components, allosteric proteins to undergo quaternary structure transitions, and molecular machines to operate in the cell.     

Computational Assessment of Protein–protein Binding Affinity by Reversely Engineering the Energetics in Protein Complexes

The cellular functions of proteins are maintained by forming diverse complexes. The stability of these complexes is quantified by the measurement of binding affinity, and mutations that alter the binding affinity can cause various diseases such as cancer and diabetes. As a result, accurate estimation of the binding stability and the effects of mutations on changes of binding affinity is a crucial step to understanding the biological functions of proteins and their dysfunctional consequences. It has been hypothesized that the stability of a protein complex is dependent not only on the residues at its binding interface by pairwise interactions but also on all other remaining residues that do not appear at the binding interface. Here, we computationally reconstruct the binding affinity by decomposing it into the contributions of interfacial residues and other non-interfacial residues in a protein complex. We further assume that the contributions of both interfacial and non-interfacial residues to the binding affinity depend on their local structural environments such as solvent-accessible surfaces and secondary structural types. The weights of all corresponding parameters are optimized by Monte-Carlo simulations. After cross-validation against a large-scale dataset, we show that the model not only shows a strong correlation between the absolute values of the experimental and calculated binding affinities, but can also be an effective approach to predict the relative changes of binding affinity from mutations. Moreover, we have found that the optimized weights of many parameters can capture the first-principle chemical and physical features of molecular recognition, therefore reversely engineering the energetics of protein complexes. These results suggest that our method can serve as a useful addition to current computational approaches for predicting binding affinity and understanding the molecular mechanism of protein–protein interactions.     

A miniaturized technique for assessing protein thermodynamics and function using fast determination of quantitative cysteine reactivity

Protein thermodynamic stability is a fundamental physical characteristic that determines biological function. Furthermore, alteration of thermodynamic stability by macromolecular interactions or biochemical modifications is a powerful tool for assessing the relationship between protein structure, stability, and biological function. High-throughput approaches for quantifying protein stability are beginning to emerge that enable thermodynamic measurements on small amounts of material, in short periods of time, and using readily accessible instrumentation. Here we present such a method, fast quantitative cysteine reactivity, which exploits the linkage between protein stability, sidechain protection by protein structure, and structural dynamics to characterize the thermodynamic and kinetic properties of proteins. In this approach, the reaction of a protected cysteine and thiol-reactive fluorogenic indicator is monitored over a gradient of temperatures after a short incubation time. These labeling data can be used to determine the midpoint of thermal unfolding, measure the temperature dependence of protein stability, quantify ligand-binding affinity, and, under certain conditions, estimate folding rate constants. Here, we demonstrate the fQCR method by characterizing these thermodynamic and kinetic properties for variants of Staphylococcal nuclease and E. coli ribose-binding protein engineered to contain single, protected cysteines. These straightforward, information-rich experiments are likely to find applications in protein engineering and functional genomics.     

Comparing interfacial dynamics in protein-protein complexes: an elastic network approach

Background  

The transient, or permanent, association of proteins to form organized complexes is one of the most common mechanisms of regulation of biological processes. Systematic physico-chemical studies of the binding interfaces have previously shown that a key mechanism for the formation/stabilization of dimers is the steric and chemical complementarity of the two semi-interfaces. The role of the fluctuation dynamics at the interface of the interacting subunits, although expectedly important, proved more elusive to characterize. The aim of the present computational study is to gain insight into salient dynamics-based aspects of protein-protein interfaces.     

A Combinatorial Histidine Scanning Library Approach to Engineer Highly pH-Dependent Protein Switches

There is growing interest in the development of protein switches, which are proteins whose function, such as binding a target molecule, can be modulated through environmental triggers. Efforts to engineer highly pH sensitive protein–protein interactions typically rely on the rational introduction of ionizable groups in the protein interface. Such experiments are typically time intensive and often sacrifice the protein's affinity at the permissive pH. The underlying thermodynamics of proton‐linkage dictate that the presence of multiple ionizable groups, which undergo a pK_a change on protein binding, are necessary to result in highly pH‐dependent binding. To test this hypothesis, a novel combinatorial histidine library was developed where every possible combination of histidine and wild‐type residue is sampled throughout the interface of a model anti‐RNase A single domain VHH antibody. Antibodies were coselected for high‐affinity binding and pH‐sensitivity using an in vitro, dual‐function selection strategy. The resulting antibodies retained near wild‐type affinity yet became highly sensitive to small decreases in pH, drastically decreasing their binding affinity, due to the incorporation of multiple histidine groups. Several trends were observed, such as histidine “hot‐spots,” which will help enhance the development of pH switch proteins as well as increase our understanding of the role of ionizable residues in protein interfaces. Overall, the combinatorial approach is rapid, general, and robust and should be capable of producing highly pH‐sensitive protein affinity reagents for a number of different applications.     

Recognition of a Key Anchor Residue by a Conserved Hydrophobic Pocket Ensures Subunit Interface Integrity in DNA Clamps

Sliding clamp proteins encircle duplex DNA and are involved in processive DNA replication and the DNA damage response. Clamp proteins are ring-shaped oligomers (dimers or trimers) and are loaded onto DNA by an ATP-dependent clamp loader complex that ruptures the interface between two adjacent subunits. Here we measured the solution dynamics of the human clamp protein, proliferating cell nuclear antigen, by monitoring the change in the fluorescence of a site-specifically labeled. To unravel the origins of clamp subunit interface stability, we carried out comprehensive comparative analysis of the interfaces of seven sliding clamps. We used computational modeling (molecular dynamic simulations and MM/GBSA binding energy decomposition analyses) to identify conserved networks of hydrophobic residues critical for clamp stability and ring-opening dynamics. The hydrophobic network is shared among clamp proteins and exhibits a “key in a keyhole” pattern where a bulky aromatic residue from one clamp subunit is anchored into a hydrophobic pocket of the opposing subunit. Bioinformatics and dynamic network analyses showed that this oligomeric latch is conserved across DNA sliding clamps from all domains of life and dictates the dynamics of clamp opening and closing.     

Action-at-a-distance interactions enhance protein binding affinity

The identification of protein mutations that enhance binding affinity may be achieved by computational or experimental means, or by a combination of the two. Sources of affinity enhancement may include improvements to the net balance of binding interactions of residues forming intermolecular contacts at the binding interface, such as packing and hydrogen-bonding interactions. Here we identify noncontacting residues that make substantial contributions to binding affinity and that also provide opportunities for mutations that increase binding affinity of the TEM1 beta-lactamase (TEM1) to the beta-lactamase inhibitor protein (BLIP). A region of BLIP not on the direct TEM1-binding surface was identified for which changes in net charge result in particularly large increases in computed binding affinity. Some mutations to the region have previously been characterized, and our results are in good correspondence with this results of that study. In addition, we propose novel mutations to BLIP that were computed to improve binding significantly without contacting TEM1 directly. This class of noncontacting electrostatic interactions could have general utility in the design and tuning of binding interactions.     

Design, production and molecular structure of a new family of artificial alpha-helicoidal repeat proteins (αRep) based on thermostable HEAT-like repeats

Repeat proteins have a modular organization and a regular architecture that make them attractive models for design and directed evolution experiments. HEAT repeat proteins, although very common, have not been used as a scaffold for artificial proteins, probably because they are made of long and irregular repeats. Here, we present and validate a consensus sequence for artificial HEAT repeat proteins. The sequence was defined from the structure-based sequence analysis of a thermostable HEAT-like repeat protein. Appropriate sequences were identified for the N- and C-caps. A library of genes coding for artificial proteins based on this sequence design, named αRep, was assembled using new and versatile methodology based on circular amplification. Proteins picked randomly from this library are expressed as soluble proteins. The biophysical properties of proteins with different numbers of repeats and different combinations of side chains in hypervariable positions were characterized. Circular dichroism and differential scanning calorimetry experiments showed that all these proteins are folded cooperatively and are very stable (T_m > 70 °C). Stability of these proteins increases with the number of repeats. Detailed gel filtration and small-angle X-ray scattering studies showed that the purified proteins form either monomers or dimers. The X-ray structure of a stable dimeric variant structure was solved. The protein is folded with a highly regular topology and the repeat structure is organized, as expected, as pairs of alpha helices. In this protein variant, the dimerization interface results directly from the variable surface enriched in aromatic residues located in the randomized positions of the repeats. The dimer was crystallized both in an apo and in a PEG-bound form, revealing a very well defined binding crevice and some structure flexibility at the interface. This fortuitous binding site could later prove to be a useful binding site for other low molecular mass partners.     

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.