Protein-protein interaction specificity of Im9 for the endonuclease toxin colicin E9 defined by homologue-scanning mutagenesis

Affinity of disordered protein complexes is modulated by entropy-energy reinforcement

Intrinsically disordered proteins (IDPs), which are very common and essential to many biological activities, sometimes function via interaction with another IDP and form a fuzzy complex, which can be highly stable. It is unclear what the biophysical forces are that govern their thermodynamics and specificity, which are essential for de novo fuzzy complex design. Here, we explored the fuzzy complex formed between ProTα and H1, which are oppositely charged IDPs, by swapping the charges between them, generating variants that have either greater polyampholytic or polyelectrolytic nature as well as different charge patterns. Charge swapping and shuffling dramatically change the affinity of the fuzzy complex, which is contributed to by both enthalpy and entropy, where the latter is dominated by counterion release.

The association between two intrinsically disordered proteins (IDPs) may produce a fuzzy complex characterized by a high binding affinity, similar to that found in the ultrastable complexes formed between two well-structured proteins. Here, using coarse-grained simulations, we quantified the biophysical forces driving the formation of such fuzzy complexes. We found that the high-affinity complex formed between the highly and oppositely charged H1 and ProTα proteins is sensitive to electrostatic interactions. We investigated 52 variants of the complex by swapping charges between the two oppositely charged proteins to produce sequences whose negatively or positively charged residue content was more homogeneous or heterogenous (i.e., polyelectrolytic or polyampholytic, having higher or lower absolute net charges, respectively) than the wild type. We also changed the distributions of oppositely charged residues within each participating sequence to produce variants in which the charges were segregated or well mixed. Both types of changes significantly affect binding affinity in fuzzy complexes, which is governed by both enthalpy and entropy. The formation of H1–ProTa is supported by an increase in configurational entropy and by entropy due to counterion release. The latter can be twice as large as the former, illustrating the dominance of counterion entropy in modulating the binding thermodynamics. Complexes formed between proteins with greater absolute net charges are more stable, both enthalpically and entropically, indicating that enthalpy and entropy have a mutually reinforcing effect. The sensitivity of the thermodynamics of the complex to net charge and the charge pattern within each of the binding constituents may provide a means to achieve binding specificity between IDPs.

Structural design principles for specific ultra-high affinity interactions between colicins/pyocins and immunity proteins

The interactions of the antibiotic proteins colicins/pyocins with immunity proteins is a seminal model system for studying protein–protein interactions and specificity. Yet, a precise and quantitative determination of which structural elements and residues determine their binding affinity and specificity is still lacking. Here, we used comparative structure-based energy calculations to map residues that substantially contribute to interactions across native and engineered complexes of colicins/pyocins and immunity proteins. We show that the immunity protein α1–α2 motif is a unique structurally-dissimilar element that restricts interaction specificity towards all colicins/pyocins, in both engineered and native complexes. This motif combines with a diverse and extensive array of electrostatic/polar interactions that enable the exquisite specificity that characterizes these interactions while achieving ultra-high affinity. Surprisingly, the divergence of these contributing colicin residues is reciprocal to residue conservation in immunity proteins. The structurally-dissimilar immunity protein α1–α2 motif is recognized by divergent colicins similarly, while the conserved immunity protein α3 helix interacts with diverse colicin residues. Electrostatics thus plays a key role in setting interaction specificity across all colicins and immunity proteins. Our analysis and resulting residue-level maps illuminate the molecular basis for these protein–protein interactions, with implications for drug development and rational engineering of these interfaces.

SKEMPI 2.0: an updated benchmark of changes in protein-protein binding energy, kinetics and thermodynamics upon mutation

Motivation

Understanding the relationship between the sequence, structure, binding energy, binding kinetics and binding thermodynamics of protein–protein interactions is crucial to understanding cellular signaling, the assembly and regulation of molecular complexes, the mechanisms through which mutations lead to disease, and protein engineering.

Results

We present SKEMPI 2.0, a major update to our database of binding free energy changes upon mutation for structurally resolved protein–protein interactions. This version now contains manually curated binding data for 7085 mutations, an increase of 133%, including changes in kinetics for 1844 mutations, enthalpy and entropy changes for 443 mutations, and 440 mutations, which abolish detectable binding.

Availability and implementation

The database is available as supplementary data and at https://life.bsc.es/pid/skempi2/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Ultrahigh specificity in a network of computationally designed protein-interaction pairs

Protein networks in all organisms comprise homologous interacting pairs. In these networks, some proteins are specific, interacting with one or a few binding partners, whereas others are multispecific and bind a range of targets. We describe an algorithm that starts from an interacting pair and designs dozens of new pairs with diverse backbone conformations at the binding site as well as new binding orientations and sequences. Applied to a high-affinity bacterial pair, the algorithm results in 18 new ones, with cognate affinities from pico- to micromolar. Three pairs exhibit 3-5 orders of magnitude switch in specificity relative to the wild type, whereas others are multispecific, collectively forming a protein-interaction network. Crystallographic analysis confirms design accuracy, including in new backbones and polar interactions. Preorganized polar interaction networks are responsible for high specificity, thus defining design principles that can be applied to program synthetic cellular interaction networks of desired affinity and specificity.

The molecular basis of ultrahigh specificity in protein-protein interactions remains obscure. The authors present a computational method to design atomically accurate new pairs exhibiting >100,000-fold specificity switches, generating a large and complex interaction network.

Molecular recognition and packing frustration in a helical protein

Biomolecular recognition entails attractive forces for the functional native states and discrimination against potential nonnative interactions that favor alternate stable configurations. The challenge posed by the competition of nonnative stabilization against native-centric forces is conceptualized as frustration. Experiment indicates that frustration is often minimal in evolved biological systems although nonnative possibilities are intuitively abundant. Much of the physical basis of minimal frustration in protein folding thus remains to be elucidated. Here we make progress by studying the colicin immunity protein Im9. To assess the energetic favorability of nonnative versus native interactions, we compute free energies of association of various combinations of the four helices in Im9 (referred to as H1, H2, H3, and H4) by extensive explicit-water molecular dynamics simulations (total simulated time > 300 μs), focusing primarily on the pairs with the largest native contact surfaces, H1-H2 and H1-H4. Frustration is detected in H1-H2 packing in that a nonnative packing orientation is significantly stabilized relative to native, whereas such a prominent nonnative effect is not observed for H1-H4 packing. However, in contrast to the favored nonnative H1-H2 packing in isolation, the native H1-H2 packing orientation is stabilized by H3 and loop residues surrounding H4. Taken together, these results showcase the contextual nature of molecular recognition, and suggest further that nonnative effects in H1-H2 packing may be largely avoided by the experimentally inferred Im9 folding transition state with native packing most developed at the H1-H4 rather than the H1-H2 interface.

Biomolecules need to recognize one another with high specificity: promoting “native” functional intermolecular binding events while avoiding detrimental “nonnative” bound configurations; i.e., “frustration”—the tendency for nonnative interactions—has to be minimized. Folding of globular proteins entails a similar discrimination. To gain physical insight, we computed the binding affinities of helical structures of the protein Im9 in various native or nonnative configurations by atomic simulations, discovering that partial packing of the Im9 core is frustrated. This frustration is overcome when the entire core of the protein is assembled, consistent with experiment indicating no significant kinetic trapping in Im9 folding. Our systematic analysis thus reveals a subtle, contextual aspect of biomolecular recognition and provides a general approach to characterize folding frustration.

The Prodomain-bound Form of Bone Morphogenetic Protein 10 Is Biologically Active on Endothelial Cells

BMP10 is highly expressed in the developing heart and plays essential roles in cardiogenesis. BMP10 deletion in mice results in embryonic lethality because of impaired cardiac development. In adults, BMP10 expression is restricted to the right atrium, though ventricular hypertrophy is accompanied by increased BMP10 expression in a rat hypertension model. However, reports of BMP10 activity in the circulation are inconclusive. In particular, it is not known whether in vivo secreted BMP10 is active or whether additional factors are required to achieve its bioactivity. It has been shown that high-affinity binding of the BMP10 prodomain to the mature ligand inhibits BMP10 signaling activity in C2C12 cells, and it was proposed that prodomain-bound BMP10 (pBMP10) complex is latent. In this study, we demonstrated that the BMP10 prodomain did not inhibit BMP10 signaling activity in multiple endothelial cells, and that recombinant human pBMP10 complex, expressed in mammalian cells and purified under native conditions, was fully active. In addition, both BMP10 in human plasma and BMP10 secreted from the mouse right atrium were fully active. Finally, we confirmed that active BMP10 secreted from mouse right atrium was in the prodomain-bound form. Our data suggest that circulating BMP10 in adults is fully active and that the reported vascular quiescence function of BMP10 in vivo is due to the direct activity of pBMP10 and does not require an additional activation step. Moreover, being an active ligand, recombinant pBMP10 may have therapeutic potential as an endothelial-selective BMP ligand, in conditions characterized by loss of BMP9/10 signaling.

SKEMPI: a Structural Kinetic and Energetic database of Mutant Protein Interactions and its use in empirical models

MOTIVATION

Empirical models for the prediction of how changes in sequence alter protein-protein binding kinetics and thermodynamics can garner insights into many aspects of molecular biology. However, such models require empirical training data and proper validation before they can be widely applied. Previous databases contained few stabilizing mutations and no discussion of their inherent biases or how this impacts model construction or validation.

RESULTS

We present SKEMPI, a database of 3047 binding free energy changes upon mutation assembled from the scientific literature, for protein-protein heterodimeric complexes with experimentally determined structures. This represents over four times more data than previously collected. Changes in 713 association and dissociation rates and 127 enthalpies and entropies were also recorded. The existence of biases towards specific mutations, residues, interfaces, proteins and protein families is discussed in the context of how the data can be used to construct predictive models. Finally, a cross-validation scheme is presented which is capable of estimating the efficacy of derived models on future data in which these biases are not present.

AVAILABILITY

The database is available online at http://life.bsc.es/pid/mutation_database/.

Designing specific protein-protein interactions using computation, experimental library screening, or integrated methods

Given the importance of protein-protein interactions for nearly all biological processes, the design of protein affinity reagents for use in research, diagnosis or therapy is an important endeavor. Engineered proteins would ideally have high specificities for their intended targets, but achieving interaction specificity by design can be challenging. There are two major approaches to protein design or redesign. Most commonly, proteins and peptides are engineered using experimental library screening and/or in vitro evolution. An alternative approach involves using protein structure and computational modeling to rationally choose sequences predicted to have desirable properties. Computational design has successfully produced novel proteins with enhanced stability, desired interactions and enzymatic function. Here we review the strengths and limitations of experimental library screening and computational structure-based design, giving examples where these methods have been applied to designing protein interaction specificity. We highlight recent studies that demonstrate strategies for combining computational modeling with library screening. The computational methods provide focused libraries predicted to be enriched in sequences with the properties of interest. Such integrated approaches represent a promising way to increase the efficiency of protein design and to engineer complex functionality such as interaction specificity.

Structure of the ultra-high-affinity colicin E2 DNase--Im2 complex

How proteins achieve high-affinity binding to a specific protein partner while simultaneously excluding all others is a major biological problem that has important implications for protein design. We report the crystal structure of the ultra-high-affinity protein-protein complex between the endonuclease domain of colicin E2 and its cognate immunity (Im) protein, Im2 (K(d)∼10(-)(15) M), which, by comparison to previous structural and biophysical data, provides unprecedented insight into how high affinity and selectivity are achieved in this model family of protein complexes. Our study pinpoints the role of structured water molecules in conjoining hotspot residues that govern stability with residues that control selectivity. A key finding is that a single residue, which in a noncognate context massively destabilizes the complex through frustration, does not participate in specificity directly but rather acts as an organizing center for a multitude of specificity interactions across the interface, many of which are water mediated.

Selection of near-native poses in CAPRI rounds 13-19

In CAPRI rounds 13-19, we submitted models that are of acceptable or higher quality for 6 of the total of 13 targets. This success builds on our record in previous CAPRI rounds. The docking problem can be divided into two steps. In the first, translational/rotational and conformational space is searched to generate a pool of docked poses; the success of this search step is measured by whether near-native poses are included in the pool. In the second step, the pool is selected for near-native poses. In our previous assessment of CAPRI results, we suggested that the search problem is largely solved; a remaining problem is to select near-native poses. Our work in these new rounds of CAPRI was guided by this assessment. To solve the selection problem, we used an assortment of criteria on the interfaces of candidate poses. In one extreme, represented by T29, with very little known interface information, our criterion for top models was based on interface prediction. Poses in which the predicted interface residues occurred in interfaces were selected. Our model 1 for T29 was of medium quality. In the other extreme, represented by T40, with reliably known interface information, our selection was solely based on such information. Nine of the ten models submitted for T40 were of high (3 models), medium (4 models), and acceptable (2 models) quality. Our strategy of mixing predicted and known interface information appears to be widely applicable for the selection of near-native poses.

The structural and energetic basis for high selectivity in a high-affinity protein-protein interaction

High-affinity, high-selectivity protein-protein interactions that are critical for cell survival present an evolutionary paradox: How does selectivity evolve when acquired mutations risk a lethal loss of high-affinity binding? A detailed understanding of selectivity in such complexes requires structural information on weak, noncognate complexes which can be difficult to obtain due to their transient and dynamic nature. Using NMR-based docking as a guide, we deployed a disulfide-trapping strategy on a noncognate complex between the colicin E9 endonuclease (E9 DNase) and immunity protein 2 (Im2), which is seven orders of magnitude weaker binding than the cognate femtomolar E9 DNase-Im9 interaction. The 1.77 A crystal structure of the E9 DNase-Im2 complex reveals an entirely noncovalent interface where the intersubunit disulfide merely supports the crystal lattice. In combination with computational alanine scanning of interfacial residues, the structure reveals that the driving force for binding is so strong that a severely unfavorable specificity contact is tolerated at the interface and as a result the complex becomes weakened through "frustration." As well as rationalizing past mutational and thermodynamic data, comparing our noncognate structure with previous cognate complexes highlights the importance of loop regions in developing selectivity and accentuates the multiple roles of buried water molecules that stabilize, ameliorate, or aggravate interfacial contacts. The study provides direct support for dual-recognition in colicin DNase-Im protein complexes and shows that weakened noncognate complexes are primed for high-affinity binding, which can be achieved by economical mutation of a limited number of residues at the interface.

Following evolutionary paths to protein-protein interactions with high affinity and selectivity

How do intricate multi-residue features such as protein-protein interfaces evolve? To address this question, we evolved a new colicin-immunity binding interaction. We started with Im9, which inhibits its cognate DNase ColE9 at 10(-14) M affinity, and evolved it toward ColE7, which it inhibits promiscuously (Kd > 10(-8) M). Iterative rounds of random mutagenesis and selection toward higher affinity for ColE7, and selectivity (against ColE9 inhibition), led to an approximately 10(5)-fold increase in affinity and a 10(8)-fold increase in selectivity. Analysis of intermediates along the evolved variants revealed that changes in the binding configuration of the Im protein uncovered a latent set of interactions, thus providing the key to the rapid divergence of new Im7 variants. Overall, protein-protein interfaces seem to share the evolvability features of enzymes, that is, the exploitation of promiscuous interactions and alternative binding configurations via 'generalist' intermediates, and the key role of compensatory stabilizing mutations in facilitating the divergence of new functions.

Analysis of nitroxide spin label motion in a protein-protein complex using multiple frequency EPR spectroscopy

X- and W-band EPR spectra, at room and low temperatures, are reported for nitroxide spin labels attached to cysteine residues selectively introduced into two proteins, the DNase domain of colicin-E9 and its immunity protein, Im9. The dynamics of each site of attachment on the individual proteins and in the tight DNase-Im9 complex have been analysed by computer simulations of the spectra using a model of Brownian dynamics trajectories for the spin label and protein. Ordering potentials have been introduced to describe mobility of labels restricted by the protein domain. Label mobility varies with position from completely immobilised, to motionally restricted and to freely rotating. Bi-modal dynamics of the spin label have been observed for several sites. We show that W-band spectra are particularly useful for detection of anisotropy of spin label motion. On complex formation significant changes are observed in the dynamics of labels at the binding interface region. This work reveals multi-frequency EPR as a sensitive and valuable tool for detecting conformational changes in protein structure and dynamics especially in protein-protein complexes.

Colicin biology

Colicins are proteins produced by and toxic for some strains of Escherichia coli. They are produced by strains of E. coli carrying a colicinogenic plasmid that bears the genetic determinants for colicin synthesis, immunity, and release. Insights gained into each fundamental aspect of their biology are presented: their synthesis, which is under SOS regulation; their release into the extracellular medium, which involves the colicin lysis protein; and their uptake mechanisms and modes of action. Colicins are organized into three domains, each one involved in a different step of the process of killing sensitive bacteria. The structures of some colicins are known at the atomic level and are discussed. Colicins exert their lethal action by first binding to specific receptors, which are outer membrane proteins used for the entry of specific nutrients. They are then translocated through the outer membrane and transit through the periplasm by either the Tol or the TonB system. The components of each system are known, and their implication in the functioning of the system is described. Colicins then reach their lethal target and act either by forming a voltage-dependent channel into the inner membrane or by using their endonuclease activity on DNA, rRNA, or tRNA. The mechanisms of inhibition by specific and cognate immunity proteins are presented. Finally, the use of colicins as laboratory or biotechnological tools and their mode of evolution are discussed.

Engineering of the Escherichia coli Im7 immunity protein as a loop display scaffold

Protein scaffolds derived from non-immunoglobulin sources are increasingly being adapted and engineered to provide unique binding molecules with a diverse range of targeting specificities. The ColE7 immunity protein (Im7) from Escherichia coli is potentially one such molecule, as it combines the advantages of (i) small size, (ii) stability conferred by a conserved four anti-parallel alpha-helical framework and (iii) availability of variable surface loops evolved to inactivate members of the DNase family of bacterial toxins, forming one of the tightest known protein-protein interactions. Here we describe initial cloning and protein expression of Im7 and its cognate partner the 15 kDa DNase domain of the colicin E7. Both proteins were produced efficiently in E.coli, and their in vitro binding interactions were validated using ELISA and biosensor. In order to assess the capacity of the Im7 protein to accommodate extensive loop region modifications, we performed extensive molecular modelling and constructed a series of loop graft variants, based on transfer of the extended CDR3 loop from the IgG1b12 antibody, which targets the gp120 antigen from HIV-1. Loop grafting in various configurations resulted in chimeric proteins exhibiting retention of the underlying framework conformation, as measured using far-UV circular dichroism spectroscopy. Importantly, there was low but measurable transfer of antigen-specific affinity. Finally, to validate Im7 as a selectable scaffold for the generation of molecular libraries, we displayed Im7 as a gene 3 fusion protein on the surface of fd bacteriophages, the most common library display format. The fusion was successfully detected using an anti-Im7 rabbit polyclonal antibody, and the recombinant phage specifically recognized the immobilized DNase. Thus, Im7 scaffold is an ideal protein display scaffold for the future generation and for the selection of libraries of novel binding proteins.

Trapping the on-pathway folding intermediate of Im7 at equilibrium

The four-helical protein Im7 folds via a rapidly formed on-pathway intermediate (k(UI)=3000 s(-1) at pH 7.0, 10 degrees C) that contains three (helices I, II and IV) of the four native alpha-helices. The relatively slow (k(IN)=300 s(-1)) conversion of this intermediate into the native structure is driven by the folding and docking of the six residue helix III onto the developing hydrophobic core. Here, we describe the structural properties of four Im7* variants designed to trap the protein in the intermediate state by disrupting the stabilising interactions formed between helix III and the rest of the protein structure. In two of these variants (I54A and L53AI54A), hydrophobic residues within helix III have been mutated to alanine, whilst in the other two mutants the sequence encompassing the native helix III was replaced by a glycine linker, three (H3G3) or six (H3G6) residues in length. All four variants were shown to be monomeric, as judged by analytical ultracentrifugation, and highly helical as measured by far-UV CD. In addition, all the variants denature co-operatively and have a stability (DeltaG(UF)) and buried hydrophobic surface area (M(UF)) similar to those of the on-pathway kinetic intermediate. Structural characterisation of these variants using 1-anilino-8-napthalene sulphonic acid (ANS) binding, near-UV CD and 1D (1)H NMR demonstrate further that the trapped intermediate ensemble is highly structured with little exposed hydrophobic surface area. Interestingly, however, the structural properties of the variants I54A and L53AI54A differ in detail from those of H3G3 and H3G6. In particular, the single tryptophan residue, located near the end of helix IV, and distant from helix III, is in a distinct environment in the two sets of mutants as judged by fluorescence, near-UV CD and the sensitivity of tryptophan fluorescence to iodide quenching. Overall, the results confirm previous kinetic analysis that demonstrated the hierarchical folding of Im7 via an on-pathway intermediate, and show that this species is a highly helical ensemble with a well-formed hydrophobic core. By contrast with the native state, however, the intermediate ensemble is flexible enough to change in response to mutation, its structural properties being tailored by residues in the sequence encompassing the native helix III.

Highly Discriminating Protein–Protein Interaction Specificities in the Context of a Conserved Binding Energy Hotspot

We explore the thermodynamic basis for high affinity binding and specificity in conserved protein complexes using colicin endonuclease-immunity protein complexes as our model system. We investigated the ability of each colicin-specific immunity protein (Im2, Im7, Im8 and Im9) to bind the endonuclease (DNase) domains of colicins E2, E7 and E8 in vitro and compared these to the previously studied colicin E9. We find that high affinity binding (Kd < or = 10(-14) M) is a common feature of cognate colicin DNase-Im protein complexes as are non-cognate protein-protein associations, which are generally 10(6)-10(8)-fold weaker. Comparative alanine scanning of Im2 and Im9 residues involved in binding the E2 DNase revealed similar behaviour to that of the two proteins binding the E9 DNase; helix III forms a conserved binding energy hotspot with specificity residues from helix II only contributing favourably in a cognate interaction, a combination we have termed as "dual recognition". Significant differences are seen, however, in the number and side-chain chemistries of specificity sites that contribute to cognate binding. In Im2, Asp33 from helix II dominates colicin E2 specificity, whereas in Im9 several hydrophobic residues, including position 33 (leucine), help define its colicin specificity. A similar distribution of specificity sites was seen using phage display where, with Im2 as the template, a library of randomised sequences was generated in helix II and the library panned against either the E2 or E9 DNase. Position 33 was the dominant specificity site recovered in all E2 DNase-selected clones, whereas a number of Im9 specificity sites were recovered in E9 DNase-selected clones, including position 33. In order to probe the relationship between biological specificity and in vitro binding affinity we compared the degree of protection afforded to bacteria against colicin E9 toxicity by a set of immunity proteins whose affinities for the E9 DNase differed by up to ten orders of magnitude. This analysis indicated that the Kd required for complete biological protection is <10(-10)M and that the "affinity window" over which the selection of novel immunity protein specificities likely evolves is 10(-6)-10(-10)M. This comprehensive survey of colicin DNase-immunity protein complexes illustrates how high affinity protein-protein interactions can be very discriminating even though binding is dominated by a conserved hotspot, with single or multiple specificity sites modulating the overall binding free energy. We discuss these results in the context of other conserved protein complexes and suggest that they point to a generic specificity mechanism in divergently evolved protein-protein interactions.

Structural analysis of the rate-limiting transition states in the folding of Im7 and Im9: similarities and differences in the folding of homologous proteins

The bacterial immunity proteins Im7 and Im9 fold with mechanisms of different kinetic complexity. Whilst Im9 folds in a two-state transition at pH 7.0 and 10 degrees C, Im7 populates an on-pathway intermediate under these conditions. In order to assess the role of sequence versus topology in the folding of these proteins, and to analyse the effect of populating an intermediate on the landscape for folding, we have determined the conformational properties of the rate-limiting transition state for Im9 folding/unfolding using Phi(F)-value analysis and have compared the results with similar data obtained previously for Im7. The data show that the rate-limiting transition states for Im9 and Im7 folding/unfolding are similar: both are compact (beta(T)=0.94 and 0.89, respectively) and contain three of the four native helices docked around a specific hydrophobic core. Significant differences are observed, however, in the magnitude of the Phi(F)-values obtained for the two proteins. Of the 20 residues studied in both proteins, ten have Phi(F)-values in Im7 that exceed those in Im9 by more than 0.2, and of these five differ by more than 0.4. The data suggest that the population of an intermediate in Im7 results in folding via a transition state ensemble that is conformationally restricted relative to that of Im9. The data are consistent with the view that topology is an important determinant of folding. Importantly, however, they also demonstrate that while the folding transition state may be conserved in homologous proteins that fold with two and three-state kinetics, the population of an intermediate can have a significant effect on the breadth of the transition state ensemble.

Structural aspects of the inhibition of DNase and rRNase colicins by their immunity proteins

Nuclease E colicins that exert their cytotoxic activity through either non-specific DNase or specific rRNase action are inhibited by immunity proteins in a high affinity interaction that gives complete protection to the producing host cell from the deleterious effects of the toxin. Previous X-ray crystallographic analysis of these systems has revealed that in both cases, the immunity protein inhibitor forms its highly stable complex with the enzyme by binding as an exosite inhibitor-adjacent to, but not obscuring, the enzyme active site. The structures of the free E9 DNase domain and its complex with an ssDNA substrate now show that inhibition is achieved without deformation of the enzyme and by occlusion of a limited number of residues of the enzyme critical in recognition and binding of the substrate that are 3' to the cleaved scissile phosphodiester. No sequence or structural similarity is evident between the two classes of cytotoxic domain, and the heterodimer interfaces are also dissimilar. Thus, whilst these structures suggest the basis for specificity in each case, they give few indications as to the basis for the remarkably strong binding that is observed. Structural analyses of complexes bearing single site mutations in the immunity protein at the heterodimer interface reveal further differences. For the DNases, a largely plastic interface is suggested, where optimal binding may be achieved in part by rigid body adjustment in the relative positions of inhibitor and enzyme. For the rRNases, a large solvent-filled cavity is found at the immunity-enzyme interface, suggesting that other considerations, such as that arising from the entropy contribution from bound water molecules, may have greater significance in the determination of rRNase complex affinity than for the DNases.

Translocation of a functional protein by a voltage-dependent ion channel

The voltage-dependent gating of the colicin channel involves a substantial structural rearrangement that results in the transfer of about 35% of the 200 residues in its pore-forming domain across the membrane. This transfer appears to represent an unusual type of protein translocation that does not depend on a large, multimeric, protein pore. To investigate the ability of this system to transport arbitrary proteins, we made use of a pair of strongly interacting proteins, either of which could serve as a translocated cargo or as a probe to detect the other. Here we show that both an 86-residue and a 134-residue hydrophilic protein inserted into the translocated segment of colicin A are themselves translocated and are functional on the trans side of the bilayer. The disparate features of these proteins suggest that the colicin channel has a general protein translocation mechanism.

Using chimeric immunity proteins to explore the energy landscape for alpha-helical protein folding

To address the role of sequence in the folding of homologous proteins, the folding and unfolding kinetics of the all-helical bacterial immunity proteins Im2 and Im9 were characterised, together with six chimeric derivatives of these proteins. We show that both Im2 and Im9 fold rapidly (k(UN)(H(2)O)) approximately 2000 s(-1) at pH 7.0, 25 degrees C) in apparent two-state transitions, through rate-limiting transition states that are highly compact (beta(TS)0.93 and 0.96, respectively). Whilst the folding and unfolding properties of three of the chimeras (Im2 (1-44)(Im9), Im2 (1-64)(Im9 )and Im2 (25-44)(Im9)) are similar to their parental counterparts, in other chimeric proteins the introduced sequence variation results in altered kinetic behaviour. At low urea concentrations, Im2 (1-29)(Im9) and Im2 (56-64)(Im9) fold in two-state transitions via transition states that are significantly less compact (beta(TS) approximately 0.7) than those characterised for the other immunity proteins presented here. At higher urea concentrations, however, the rate-limiting transition state for these two chimeras switches or moves to a more compact species (beta(TS) approximately 0.9). Surprisingly, Im2 (30-64)(Im9) populates a highly collapsed species (beta(I)=0.87) in the dead-time (2.5 ms) of stopped flow measurements. These data indicate that whilst topology may place significant constraints on the folding process, specific inter-residue interactions, revealed here through multiple sequence changes, can modulate the ruggedness of the folding energy landscape.

NMR investigation of the interaction of the inhibitor protein Im9 with its partner DNase

The bacterial toxin colicin E9 is secreted by producing Escherichia coli cells with its 9.5 kDa inhibitor protein Im9 bound tightly to its 14.5 kDa C-terminal DNase domain. Double- and triple-resonance NMR spectra of the 24 kDa complex of uniformly 13C and 15N labeled Im9 bound to the unlabeled DNase domain have provided sufficient constraints for the solution structure of the bound Im9 to be determined. For the final ensemble of 20 structures, pairwise RMSDs for residues 3-84 were 0.76 +/- 0.14 A for the backbone atoms and 1.36 +/- 0.15 A for the heavy atoms. Representative solution structures of the free and bound Im9 are highly similar, with backbone and heavy atom RMSDs of 1.63 and 2.44 A, respectively, for residues 4-83, suggesting that binding does not cause a major conformational change in Im9. The NMR studies have also allowed the DNase contact surface on Im9 to be investigated through changes in backbone chemical shifts and NOEs between the two proteins determined from comparisons of 1H-1H-13C NOESY-HSQC spectra with and without 13C decoupling. The NMR-defined interface agrees well with that determined in a recent X-ray structure analysis with the major difference being that a surface loop of Im9, which is at the interface, has a different conformation in the solution and crystal structures. Tyr54, a key residue on the interface, is shown to exhibit NMR characteristics indicative of slow rotational flipping. A mechanistic description of the influence binding of Im9 has on the dynamic behavior of E9 DNase, which is known to exist in two slowly interchanging conformers in solution, is proposed.

Specificity in protein-protein interactions: the structural basis for dual recognition in endonuclease colicin-immunity protein complexes

Bacteria producing endonuclease colicins are protected against their cytotoxic activity by virtue of a small immunity protein that binds with high affinity and specificity to inactivate the endonuclease. DNase binding by the immunity protein occurs through a "dual recognition" mechanism in which conserved residues from helix III act as the binding-site anchor, while variable residues from helix II define specificity. We now report the 1.7 A crystal structure of the 24.5 kDa complex formed between the endonuclease domain of colicin E9 and its cognate immunity protein Im9, which provides a molecular rationale for this mechanism. Conserved residues of Im9 form a binding-energy hotspot through a combination of backbone hydrogen bonds to the endonuclease, many via buried solvent molecules, and hydrophobic interactions at the core of the interface, while the specificity-determining residues interact with corresponding specificity side-chains on the enzyme. Comparison between the present structure and that reported recently for the colicin E7 endonuclease domain in complex with Im7 highlights how specificity is achieved by very different interactions in the two complexes, predominantly hydrophobic in nature in the E9-Im9 complex but charged in the E7-Im7 complex. A key feature of both complexes is the contact between a conserved tyrosine residue from the immunity proteins (Im9 Tyr54) with a specificity residue on the endonuclease directing it toward the specificity sites of the immunity protein. Remarkably, this tyrosine residue and its neighbour (Im9 Tyr55) are the pivots of a 19 degrees rigid-body rotation that relates the positions of Im7 and Im9 in the two complexes. This rotation does not affect conserved immunity protein interactions with the endonuclease but results in different regions of the specificity helix being presented to the enzyme.

Hierarchical order of critical residues on the immunity-determining region of the Im7 protein which confer specific immunity to its cognate colicin

The directed mutagenesis study of the Im7 protein of colicin E7 revealed that three residues, D31, D35, and E39, located in the loop 1 and helix 2 regions of the protein were critical for initiating the complex formation with its cognate colicin E7. Interestingly, the importance of these three critical residues in conferring specific immunity to its own colicin was exhibited in a hierarchical order, respectively. Moreover, we found that existence of the three critical residues was common among the DNase-type Im proteins. Most likely the three residues of the DNase-type immunity proteins are critical for initiating the unique protein-protein interactions with their cognate colicin. In addition, replacement of the helix 2 of Im7 by the corresponding region of Im8 produced a phenotype of the mutant protein very similar to that of Im8. This result suggests that the DNase-type Im proteins indeed share a "homologous-structural framework" and evolution of the Im proteins may be engendered by minor amino acid changes in this specific immunity-determining region without causing structural alteration of the proteins.

The crystal structure of the DNase domain of colicin E7 in complex with its inhibitor Im7 protein

BACKGROUND

Colicin E7 (ColE7) is one of the bacterial toxins classified as a DNase-type E-group colicin. The cytotoxic activity of a colicin in a colicin-producing cell can be counteracted by binding of the colicin to a highly specific immunity protein. This biological event is a good model system for the investigation of protein recognition.

RESULTS

The crystal structure of a one-to-one complex between the DNase domain of colicin E7 and its cognate immunity protein Im7 has been determined at 2.3 A resolution. Im7 in the complex is a varied four-helix bundle that is identical to the structure previously determined for uncomplexed Im7. The structure of the DNase domain of ColE7 displays a novel alpha/beta fold and contains a Zn2+ ion bound to three histidine residues and one water molecule in a distorted tetrahedron geometry. Im7 has a V-shaped structure, extending two arms to clamp the DNase domain of ColE7. One arm (alpha1(*)-loop12-alpha2(*); where * represents helices in Im7) is located in the region that displays the greatest sequence variation among members of the immunity proteins in the same subfamily. This arm mainly uses acidic sidechains to interact with the basic sidechains in the DNase domain of ColE7. The other arm (loop 23-alpha3(*)-loop 34) is more conserved and it interacts not only with the sidechain but also with the mainchain atoms of the DNase domain of ColE7.

CONCLUSIONS

The protein interfaces between the DNase domain of ColE7 and Im7 are charge-complementary and charge interactions contribute significantly to the tight and specific binding between the two proteins. The more variable arm in Im7 dominates the binding specificity of the immunity protein to its cognate colicin. Biological and structural data suggest that the DNase active site for ColE7 is probably near the metal-binding site.

Enzymological characterization of the nuclease domain from the bacterial toxin colicin E9 from Escherichia coli

The cytotoxicity of the bacterial toxin colicin E9 is due to a non-specific DNase that penetrates the cytoplasm of the infected organism and causes cell death. We report the first enzymological characterization of the overexpressed and purified 15 kDa DNase domain (E9 DNase) from this class of toxin. CD spectroscopy shows the E9 DNase to be structured in solution, and analytical ultracentrifugation data indicate that the enzyme is a monomer. The nuclease activity of the E9 DNase was compared with the well-studied, non-specific DNase I by using a spectrophotometric assay with calf thymus DNA as the substrate. Both enzymes require divalent metal ions for activity but, unlike DNase I, the E9 DNase is not activated by Ca2+ ions. Somewhat surprisingly, the E9 DNase shows optimal activity and linear kinetics in the presence of transition metals such as Ni2+ and Co2+ but displays non-linear kinetics with metals such as Mg2+ and Ca2+. Conversely, Ni2+ and other transition metals showed poor activity in a plasmid-based nicking assay, yielding significant amounts of linearized plasmid, whereas Mg2+ was very active, with the main intermediate being open-circle DNA. The results suggest that, on entry into bacterial cells, the E9 DNase is likely to exhibit primarily Mg2+-dependent nicking activity against chromosomal DNA, although other metals could also be utilized to introduce both single- and double-strand cleavages.

A structural comparison of the colicin immunity proteins Im7 and Im9 gives new insights into the molecular determinants of immunity-protein specificity

We report the first detailed comparison of two immunity proteins which, in conjunction with recent protein engineering data, begins to explain how these structurally similar proteins are able to bind and inhibit the endonuclease domain of colicin E9 (E9 DNase) with affinities that differ by 12 orders of magnitude. In the present work, we have determined the X-ray structure of the Escherichia coli colicin E7 immunity protein Im7 to 2.0 A resolution by molecular replacement, using as a trial model the recently determined NMR solution structure of Im9. Whereas the two proteins adopt similar four-helix structures, subtle structural differences, in particular involving a conserved tyrosine residue critical for E9 DNase binding, and the identity of key residues in the specificity helix, lie at the heart of their markedly different ability to bind the E9 DNase. Two other crystal structures were reported recently for Im7; in one, Im7 was a monomer and was very similar to the structure reported here, whereas in the other it was a dimer to which functional significance was assigned. Since this previous work suggested that Im7 could exist either as a monomer or a dimer, we used analytical ultracentrifugation to investigate this question further. Under a variety of solution conditions, we found that Im7 only ever exists in solution as a monomer, even up to protein concentrations of 15 mg/ml, casting doubt on the functional significance of the crystallographically observed dimer. This work provides a structural framework with which we can understand immunity-protein specificity, and in addition we believe it to be the first successfully refined crystal structure solved by molecular replacement using an NMR trial model with less than 100% sequence identity.

Immunity proteins and their specificity for endonuclease colicins: telling right from wrong in protein-protein recognition

Immunity proteins inhibit colicins, protein toxins released by bacteria during times of environmental stress, by binding and inactivating their cytotoxic domains. This protects the producing organism as it attempts to kill off competing bacteria. The cytotoxic domains of related colicins share a high degree of sequence identity, as do their corresponding immunity proteins, yet specificity and affinity are also high, with little non-cognate biological cross-protection evident under physiological conditions. We review recent work on DNase-specific immunity proteins, which shows that, although both cognate and non-cognate proteins can bind a single toxin, their affinities can differ by as much as 12 orders of magnitude. We have termed this mode of binding dual recognition, because the DNase-binding surface of an immunity protein is made up of two components, one conserved and the other variable. The strength of the binding interaction is dominated by the conserved residues, while neighbouring variable residues control specificity. Similar dual recognition systems may exist in other biological contexts, particularly where a protein must discriminate the right binding partner from numerous, structurally homologous alternatives.

Measuring protein-protein interactions

The binding of one protein to another provokes a variety of biophysical changes that can then be used as a measure of the binding reaction. Optical spectroscopy, particularly fluorescence, is the most flexible technique, but surface plasmon resonance biosensors, microcalorimetry and mass spectroscopy have recently shown significant development.