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

Evolutionary paths that link orthogonal pairs of binding proteins

Some protein binding pairs exhibit extreme specificities that functionally insulate them from homologs. Such pairs evolve mostly by accumulating single-point mutations, and mutants are selected if their affinity exceeds the threshold required for function1-4. Thus, homologous and high-specificity binding pairs bring to light an evolutionary conundrum: how does a new specificity evolve while maintaining the required affinity in each intermediate5,6? Until now, a fully functional single-mutation path that connects two orthogonal pairs has only been described where the pairs were mutationally close thus enabling experimental enumeration of all intermediates2. We present an atomistic and graph-theoretical framework for discovering low molecular strain single-mutation paths that connect two extant pairs, enabling enumeration beyond experimental capability. We apply it to two orthogonal bacterial colicin endonuclease-immunity pairs separated by 17 interface mutations7. We were not able to find a strain-free and functional path in the sequence space defined by the two extant pairs. But including mutations that bridge amino acids that cannot be exchanged through single-nucleotide mutations led us to a strain-free 19-mutation trajectory that is completely viable in vivo. Our experiments show that the specificity switch is remarkably abrupt, resulting from only one radical mutation on each partner. Furthermore, each of the critical specificity-switch mutations increases fitness, demonstrating that functional divergence could be driven by positive Darwinian selection. These results reveal how even radical functional changes in an epistatic fitness landscape may evolve.

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.

Polymorphic Toxins and Their Immunity Proteins: Diversity, Evolution, and Mechanisms of Delivery

All bacteria must compete for growth niches and other limited environmental resources. These existential battles are waged at several levels, but one common strategy entails the transfer of growth-inhibitory protein toxins between competing cells. These antibacterial effectors are invariably encoded with immunity proteins that protect cells from intoxication by neighboring siblings. Several effector classes have been described, each designed to breach the cell envelope of target bacteria. Although effector architectures and export pathways tend to be clade specific, phylogenetically distant species often deploy closely related toxin domains. Thus, diverse competition systems are linked through a common reservoir of toxin-immunity pairs that is shared via horizontal gene transfer. These toxin-immunity protein pairs are extraordinarily diverse in sequence, and this polymorphism underpins an important mechanism of self/nonself discrimination in bacteria. This review focuses on the structures, functions, and delivery mechanisms of polymorphic toxin effectors that mediate bacterial competition.

Performance of human and server prediction in CAPRI rounds 38-45

CAPRI challenges offer a variety of blind tests for protein-protein interaction prediction. In CAPRI Rounds 38-45, we generated a set of putative binding modes for each target with an FFT-based docking algorithm, and then scored and ranked these binding modes with a proprietary scoring function, ITScorePP. We have also developed a novel web server, Rebipp. The algorithm utilizes information retrieval to identify relevant biological information to significantly reduce the search space for a particular protein. In parallel, we have also constructed a GPU-based docking server, MDockPP, for protein-protein complex structure prediction. Here, the performance of our protocol in CAPRI rounds 38-45 is reported, which include 16 docking and scoring targets. Among them, three targets contain multiple interfaces: Targets 124, 125, and 136 have 2, 4, and 3 interfaces, respectively. In the predictor experiments, we predicted correct binding modes for nine targets, including one high-accuracy interface, six medium-accuracy binding modes, and six acceptable-accuracy binding modes. For the docking server prediction experiments, we predicted correct binding modes for eight targets, including one high-accuracy, three medium-accuracy, and five acceptable-accuracy binding modes.

Molecular basis of ALK1-mediated signalling by BMP9/BMP10 and their prodomain-bound forms

Activin receptor-like kinase 1 (ALK1)-mediated endothelial cell signalling in response to bone morphogenetic protein 9 (BMP9) and BMP10 is of significant importance in cardiovascular disease and cancer. However, detailed molecular mechanisms of ALK1-mediated signalling remain unclear. Here, we report crystal structures of the BMP10:ALK1 complex at 2.3 Å and the prodomain-bound BMP9:ALK1 complex at 3.3 Å. Structural analyses reveal a tripartite recognition mechanism that defines BMP9 and BMP10 specificity for ALK1, and predict that crossveinless 2 is not an inhibitor of BMP9, which is confirmed by experimental evidence. Introduction of BMP10-specific residues into BMP9 yields BMP10-like ligands with diminished signalling activity in C2C12 cells, validating the tripartite mechanism. The loss of osteogenic signalling in C2C12 does not translate into non-osteogenic activity in vivo and BMP10 also induces bone-formation. Collectively, these data provide insight into ALK1-mediated BMP9 and BMP10 signalling, facilitating therapeutic targeting of this important pathway.

Computational studies of protein-protein dissociation by statistical potential and coarse-grained simulations: a case study on interactions between colicin E9 endonuclease and immunity proteins

Proteins carry out their diverse functions in cells by forming interactions with each other. The dynamics of these interactions are quantified by the measurement of association and dissociation rate constants. Relative to the efforts made to model the association of biomolecules, little has been studied to understand the principles of protein complex dissociation. Using the interaction between colicin E9 endonucleases and immunity proteins as a test system, here we develop a coarse-grained simulation method to explore the dissociation mechanisms of protein complexes. The interactions between proteins in the complex are described by the knowledge-based potential that was constructed by the statistics from available protein complexes in the structural database. Our study provides the supportive evidences to the dual recognition mechanism for the specificity of binding between E9 DNase and immunity proteins, in which the conserved residues of helix III of Im2 and Im9 proteins act as the anchor for binding, while the sequence variations in helix II make positive or negative contributions to specificity. Beyond that, we further suggest that this binding specificity is rooted in the process of complex dissociation instead of association. While we increased the flexibility of protein complexes, we further found that they are less prone to dissociation, suggesting that conformational fluctuations of protein complexes play important functional roles in regulating their binding and dissociation. Our studies therefore bring new insights to the molecule mechanisms of protein-protein interactions, while the method can serve as a new addition to a suite of existing computational tools for the simulations of protein complexes.

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.

The Inclusion of Water Molecules in Residue Interaction Networks Identifies Additional Central Residues

The relevance of water molecules for the recognition and the interaction of biomolecules is widely appreciated. In this paper we address the role that water molecules associated to protein complexes play for the functional relevance of residues by considering their residue interaction networks (RINs). These are commonly defined on the basis of the amino acid composition of the proteins themselves, disregarding the solvation state of the protein. We determine properties of the RINs of two protein complexes, colicin E2/Im2 and barnase/barstar, with and without associated water molecules, using a previously developed methodology and its associated application RINspector. We find that the inclusion of water molecules in RINs leads to an increase in the number of central residues which adds a novel mechanism to the relevance of water molecules for protein function.

Evolution of In Silico Strategies for Protein-Protein Interaction Drug Discovery

The advent of advanced molecular modeling software, big data analytics, and high-speed processing units has led to the exponential evolution of modern drug discovery and better insights into complex biological processes and disease networks. This has progressively steered current research interests to understanding protein-protein interaction (PPI) systems that are related to a number of relevant diseases, such as cancer, neurological illnesses, metabolic disorders, etc. However, targeting PPIs are challenging due to their "undruggable" binding interfaces. In this review, we focus on the current obstacles that impede PPI drug discovery, and how recent discoveries and advances in in silico approaches can alleviate these barriers to expedite the search for potential leads, as shown in several exemplary studies. We will also discuss about currently available information on PPI compounds and systems, along with their usefulness in molecular modeling. Finally, we conclude by presenting the limits of in silico application in drug discovery and offer a perspective in the field of computer-aided PPI drug discovery.

Novel Immunity Proteins Associated with Colicin M-like Bacteriocins Exhibit Promiscuous Protection in Pseudomonas

Bacteriocins related to colicin M, acting via cleavage of the cell wall precursor lipid II, have been characterized in γ- and β-proteobacteria. Depending on the species, immunity is provided by either an inner membrane-anchored periplasmic protein or by an integral membrane protein. In Pseudomonas however, the immunity partner of colicin M-like bacteriocins remains unknown. Based on an in silico analysis in pseudomonad genomes, we here identify a gene encoding a putative immunity partner that represents a novel type of integral membrane protein (PmiA, Pseudomonas colicin M-like immunity type A). By heterologous expression of pmiA genes in susceptible strains, we show that immunity to colicin M-like bacteriocins is indeed provided by the cognate PmiA. Sequence homology among PmiA proteins is essentially absent, except for a short motif with a conserved periplasm-exposed aspartate residue. However, PmiA's protective function is not abolished by changing this acidic residue to the uncharged alanine. Immunity by PmiAs appears promiscuous to the extent that PmiA homologs from a clade sharing <40% pairwise amino acid identity, equally provide protection against the bacteriocin linked to the original PmiA. This study shows that multiple immunity factors have evolved independently to silence lipid II-targeting enzymatic bacteriocins. Their relaxed bacteriocin immunization capacity contrasts to the strict specificity of immunity proteins shielding the enzymatic domain of nuclease bacteriocins. The nature of associated immune functions needs consideration when using such natural protein antibiotics or designing novel variants.

Diversification of β-Augmentation Interactions between CDI Toxin/Immunity Proteins

Contact-dependent growth inhibition (CDI) is a widespread mechanism of inter-bacterial competition mediated by the CdiB/CdiA family of two-partner secretion proteins. CdiA effectors carry diverse C-terminal toxin domains (CdiA-CT), which are delivered into neighboring target cells to inhibit growth. CDI(+) bacteria also produce CdiI immunity proteins that bind specifically to cognate CdiA-CT toxins and protect the cell from auto-inhibition. Here, we compare the structures of homologous CdiA-CT/CdiI complexes from Escherichia coli EC869 and Yersinia pseudotuberculosis YPIII to explore the evolution of CDI toxin/immunity protein interactions. Both complexes share an unusual β-augmentation interaction, in which the toxin domain extends a β-hairpin into the immunity protein to complete a six-stranded anti-parallel sheet. However, the specific contacts differ substantially between the two complexes. The EC869 β-hairpin interacts mainly through direct H-bond and ion-pair interactions, whereas the YPIII β-hairpin pocket contains more hydrophobic contacts and a network of bridging water molecules. In accord with these differences, we find that each CdiI protein only protects target bacteria from its cognate CdiA-CT toxin. The compact β-hairpin binding pocket within the immunity protein represents a tractable system for the rationale design of small molecules to block CdiA-CT/CdiI complex formation. We synthesized a macrocyclic peptide mimic of the β-hairpin from EC869 toxin and solved its structure in complex with cognate immunity protein. These latter studies suggest that small molecules could potentially be used to disrupt CDI toxin/immunity complexes.

Structures of the Ultra-High-Affinity Protein-Protein Complexes of Pyocins S2 and AP41 and Their Cognate Immunity Proteins from Pseudomonas aeruginosa

How ultra-high-affinity protein–protein interactions retain high specificity is still poorly understood. The interaction between colicin DNase domains and their inhibitory immunity (Im) proteins is an ultra-high-affinity interaction that is essential for the neutralisation of endogenous DNase catalytic activity and for protection against exogenous DNase bacteriocins. The colicin DNase–Im interaction is a model system for the study of high-affinity protein–protein interactions. However, despite the fact that closely related colicin-like bacteriocins are widely produced by Gram-negative bacteria, this interaction has only been studied using colicins from Escherichia coli. In this work, we present the first crystal structures of two pyocin DNase–Im complexes from Pseudomonas aeruginosa, pyocin S2 DNase–ImS2 and pyocin AP41 DNase–ImAP41. These structures represent divergent DNase–Im subfamilies and are important in extending our understanding of protein–protein interactions for this important class of high-affinity protein complex. A key finding of this work is that mutations within the immunity protein binding energy hotspot, helix III, are tolerated by complementary substitutions at the DNase–Immunity protein binding interface. Im helix III is strictly conserved in colicins where an Asp forms polar interactions with the DNase backbone. ImAP41 contains an Asp-to-Gly substitution in helix III and our structures show the role of a co-evolved substitution where Pro in DNase loop 4 occupies the volume vacated and removes the unfulfilled hydrogen bond. We observe the co-evolved mutations in other DNase–Immunity pairs that appear to underpin the split of this family into two distinct groups.

•

We have identified two different bacteriocin DNase–Im subfamilies.

•

First structures of pyocin DNase domains in complex with neutralising Im proteins.

•

The subfamilies are characterised by distinct Im helix III motifs.

•

ImAP41 lacks the key Asp in Im helix III and one of the conserved interfacial waters.

•

New DNase–Im family expands the region that governs bacteriocin selectivity.

Flexibility and small pockets at protein-protein interfaces: New insights into druggability

The transient assembly of multiprotein complexes mediates many aspects of cell regulation and signalling in living organisms. Modulation of the formation of these complexes through targeting protein-protein interfaces can offer greater selectivity than the inhibition of protein kinases, proteases or other post-translational regulatory enzymes using substrate, co-factor or transition state mimetics. However, capitalising on protein-protein interaction interfaces as drug targets has been hindered by the nature of interfaces that tend to offer binding sites lacking the well-defined large cavities of classical drug targets. In this review we posit that interfaces formed by concerted folding and binding (disorder-to-order transitions on binding) of one partner and other examples of interfaces where a protein partner is bound through a continuous epitope from a surface-exposed helix, flexible loop or chain extension may be more tractable for the development of "orthosteric", competitive chemical modulators; these interfaces tend to offer small-volume but deep pockets and/or larger grooves that may be bound tightly by small chemical entities. We discuss examples of such protein-protein interaction interfaces for which successful chemical modulators are being developed.

Mathematical model accurately predicts protein release from an affinity-based delivery system

Affinity-based controlled release modulates the delivery of protein or small molecule therapeutics through transient dissociation/association. To understand which parameters can be used to tune release, we used a mathematical model based on simple binding kinetics. A comprehensive asymptotic analysis revealed three characteristic regimes for therapeutic release from affinity-based systems. These regimes can be controlled by diffusion or unbinding kinetics, and can exhibit release over either a single stage or two stages. This analysis fundamentally changes the way we think of controlling release from affinity-based systems and thereby explains some of the discrepancies in the literature on which parameters influence affinity-based release. The rate of protein release from affinity-based systems is determined by the balance of diffusion of the therapeutic agent through the hydrogel and the dissociation kinetics of the affinity pair. Equations for tuning protein release rate by altering the strength (KD) of the affinity interaction, the concentration of binding ligand in the system, the rate of dissociation (koff) of the complex, and the hydrogel size and geometry, are provided. We validated our model by collapsing the model simulations and the experimental data from a recently described affinity release system, to a single master curve. Importantly, this mathematical analysis can be applied to any single species affinity-based system to determine the parameters required for a desired release profile.

A "fuzzy"-logic language for encoding multiple physical traits in biomolecules

To carry out their activities, biological macromolecules balance different physical traits, such as stability, interaction affinity, and selectivity. How such often opposing traits are encoded in a macromolecular system is critical to our understanding of evolutionary processes and ability to design new molecules with desired functions. We present a framework for constraining design simulations to balance different physical characteristics. Each trait is represented by the equilibrium fractional occupancy of the desired state relative to its alternatives, ranging from none to full occupancy, and the different traits are combined using Boolean operators to effect a "fuzzy"-logic language for encoding any combination of traits. In another paper, we presented a new combinatorial backbone design algorithm AbDesign where the fuzzy-logic framework was used to optimize protein backbones and sequences for both stability and binding affinity in antibody-design simulation. We now extend this framework and find that fuzzy-logic design simulations reproduce sequence and structure design principles seen in nature to underlie exquisite specificity on the one hand and multispecificity on the other hand. The fuzzy-logic language is broadly applicable and could help define the space of tolerated and beneficial mutations in natural biomolecular systems and design artificial molecules that encode complex characteristics.

Colicin import into E. coli cells: a model system for insights into the import mechanisms of bacteriocins

Bacteriocins are a diverse group of ribosomally synthesized protein antibiotics produced by most bacteria. They range from small lanthipeptides produced by lactic acid bacteria to much larger multi domain proteins of Gram negative bacteria such as the colicins from Escherichia coli. For activity bacteriocins must be released from the producing cell and then bind to the surface of a sensitive cell to instigate the import process leading to cell death. For over 50years, colicins have provided a working platform for elucidating the structure/function studies of bacteriocin import and modes of action. An understanding of the processes that contribute to the delivery of a colicin molecule across two lipid membranes of the cell envelope has advanced our knowledge of protein-protein interactions (PPI), protein-lipid interactions and the role of order-disorder transitions of protein domains pertinent to protein transport. In this review, we provide an overview of the arrangement of genes that controls the synthesis and release of the mature protein. We examine the uptake processes of colicins from initial binding and sequestration of binding partners to crossing of the outer membrane, and then discuss the translocation of colicins through the cell periplasm and across the inner membrane to their cytotoxic site of action. This article is part of a Special Issue entitled: Protein trafficking and secretion in bacteria. Guest Editors: Anastassios Economou and Ross Dalbey.

Blind prediction of interfacial water positions in CAPRI

We report the first assessment of blind predictions of water positions at protein-protein interfaces, performed as part of the critical assessment of predicted interactions (CAPRI) community-wide experiment. Groups submitting docking predictions for the complex of the DNase domain of colicin E2 and Im2 immunity protein (CAPRI Target 47), were invited to predict the positions of interfacial water molecules using the method of their choice. The predictions-20 groups submitted a total of 195 models-were assessed by measuring the recall fraction of water-mediated protein contacts. Of the 176 high- or medium-quality docking models-a very good docking performance per se-only 44% had a recall fraction above 0.3, and a mere 6% above 0.5. The actual water positions were in general predicted to an accuracy level no better than 1.5 Å, and even in good models about half of the contacts represented false positives. This notwithstanding, three hotspot interface water positions were quite well predicted, and so was one of the water positions that is believed to stabilize the loop that confers specificity in these complexes. Overall the best interface water predictions was achieved by groups that also produced high-quality docking models, indicating that accurate modelling of the protein portion is a determinant factor. The use of established molecular mechanics force fields, coupled to sampling and optimization procedures also seemed to confer an advantage. Insights gained from this analysis should help improve the prediction of protein-water interactions and their role in stabilizing protein complexes.

Using the concept of transient complex for affinity predictions in CAPRI rounds 20-27 and beyond

Predictions of protein-protein binders and binding affinities have traditionally focused on features pertaining to the native complexes. In developing a computational method for predicting protein-protein association rate constants, we introduced the concept of transient complex after mapping the interaction energy surface. The transient complex is located at the outer boundary of the bound-state energy well, having near-native separation and relative orientation between the subunits but not yet formed most of the short-range native interactions. We found that the width of the binding funnel and the electrostatic interaction energy of the transient complex are among the features predictive of binders and binding affinities. These ideas were very promising for the five affinity-related targets (T43-45, 55, and 56) of CAPRI rounds 20-27. For T43, we ranked the single crystallographic complex as number 1 and were one of only two groups that clearly identified that complex as a true binder; for T44, we ranked the only design with measurable binding affinity as number 4. For the nine docking targets, continuing on our success in previous CAPRI rounds, we produced 10 medium-quality models for T47 and acceptable models for T48 and T49. We conclude that the interaction energy landscape and the transient complex in particular will complement existing features in leading to better prediction of binding affinities.

A force-activated trip switch triggers rapid dissociation of a colicin from its immunity protein

A single-molecule force study shows that rapid dissociation of a high-affinity protein interaction can be triggered by site-specific remodelling of one protein partner, and that prevention of remodelling maintains avidity.

Colicins are protein antibiotics synthesised by Escherichia coli strains to target and kill related bacteria. To prevent host suicide, colicins are inactivated by binding to immunity proteins. Despite their high avidity (K_d≈fM, lifetime ≈4 days), immunity protein release is a pre-requisite of colicin intoxication, which occurs on a timescale of minutes. Here, by measuring the dynamic force spectrum of the dissociation of the DNase domain of colicin E9 (E9) and immunity protein 9 (Im9) complex using an atomic force microscope we show that application of low forces (<20 pN) increases the rate of complex dissociation 10⁶-fold, to a timescale (lifetime ≈10 ms) compatible with intoxication. We term this catastrophic force-triggered increase in off-rate a trip bond. Using mutational analysis, we elucidate the mechanism of this switch in affinity. We show that the N-terminal region of E9, which has sparse contacts with the hydrophobic core, is linked to an allosteric activator region in E9 (residues 21–30) whose remodelling triggers immunity protein release. Diversion of the force transduction pathway by the introduction of appropriately positioned disulfide bridges yields a force resistant complex with a lifetime identical to that measured by ensemble techniques. A trip switch within E9 is ideal for its function as it allows bipartite complex affinity, whereby the stable colicin:immunity protein complex required for host protection can be readily converted to a kinetically unstable complex whose dissociation is necessary for cellular invasion and competitor death. More generally, the observation of two force phenotypes for the E9:Im9 complex demonstrates that force can re-sculpt the underlying energy landscape, providing new opportunities to modulate biological reactions in vivo; this rationalises the commonly observed discrepancy between off-rates measured by dynamic force spectroscopy and ensemble methods.

Many proteins interact with other proteins as part of their function. One method of modulating the activity of protein complexes is to break them apart. Some complexes, however, are extremely kinetically stable and it is unclear how these can dissociate on a biologically relevant timescale. In this study we address this question using protein complexes between colicin E9 (a bacterial toxin) and its immunity protein Im9. These highly avid complexes (with a lifetime of days) must be broken apart for colicin to be activated. By using single-molecule force methods we show that pulling on one end of colicin E9 drastically destabilises the complex so that it dissociates a million-fold faster than its intrinsic rate. We then show that preventing this destabilisation (by the insertion of cross-links that pin the N-terminus of E9 in place) yields a kinetically stable complex. It has previously been postulated that force can destabilise a protein complex by partially unfolding one or more binding partners. Our work provides new experimental evidence that shows this is the case and provides a mechanism for this phenomenon, which we term a trip bond. For the E9:Im9 complex, trip bond behaviour allows a stable complex to be rapidly dissociated by application of a surprisingly small force.

Role of tyrosine hot-spot residues at the interface of colicin E9 and immunity protein 9: a comparative free energy simulation study

The endonuclease activity of the bacterial colicin 9 enzyme is controlled by the specific and high-affinity binding of immunity protein 9 (Im9). Molecular dynamics simulation studies in explicit solvent were used to investigate the free energy change associated with the mutation of two hot-spot interface residues [tyrosine (Tyr): Tyr54 and Tyr55] of Im9 to Ala. In addition, the effect of several other mutations (Leu33Ala, Leu52Ala, Val34Ala, Val37Ala, Ser48Ala, and Ile53Ala) with smaller influence on binding affinity was also studied. Good qualitative agreement of calculated free energy changes and experimental data on binding affinity of the mutations was observed. The simulation studies can help to elucidate the molecular details on how the mutations influence protein-protein binding affinity. The role of solvent and conformational flexibility of the partner proteins was studied by comparing the results in the presence or absence of solvent and with or without positional restraints. Restriction of the conformational mobility of protein partners resulted in significant changes of the calculated free energies but of similar magnitude for isolated Im9 and for the complex and therefore in only modest changes of binding free energy differences. Although the overall binding free energy change was similar for the two Tyr-Ala mutations, the physical origin appeared to be different with solvation changes contributing significantly to the Tyr55Ala mutation and to a loss of direct protein-protein interactions dominating the free energy change due to the Tyr54Ala mutation.

Conformational properties of the unfolded state of Im7 in nondenaturing conditions

The unfolded ensemble in aqueous solution represents the starting point of protein folding. Characterisation of this species is often difficult since the native state is usually predominantly populated at equilibrium. Previous work has shown that the four-helix protein, Im7 (immunity protein 7), folds via an on-pathway intermediate. While the transition states and folding intermediate have been characterised in atomistic detail, knowledge of the unfolded ensemble under the same ambient conditions remained sparse. Here, we introduce destabilising amino acid substitutions into the sequence of Im7, such that the unfolded state becomes predominantly populated at equilibrium in the absence of denaturant. Using far- and near-UV CD, fluorescence, urea titration and heteronuclear NMR experiments, we show that three amino acid substitutions (L18A-L19A-L37A) are sufficient to prevent Im7 folding, such that the unfolded state is predominantly populated at equilibrium. Using measurement of chemical shifts, (15)N transverse relaxation rates and sedimentation coefficients, we show that the unfolded species of L18A-L19A-L37A deviates significantly from random-coil behaviour. Specifically, we demonstrate that this unfolded species is compact (R(h)=25 Å) relative to the urea-denatured state (R(h)≥30 Å) and contains local clusters of hydrophobic residues in regions that correspond to the four helices in the native state. Despite these interactions, there is no evidence for long-range stabilising tertiary interactions or persistent helical structure. The results reveal an unfolded ensemble that is conformationally restricted in regions of the polypeptide chain that ultimately form helices I, II and IV in the native state.

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.

Nuclease colicins and their immunity proteins

It is more than 80 years since Gratia first described 'a remarkable antagonism between two strains of Escherichia coli'. Shown subsequently to be due to the action of proteins (or peptides) produced by one bacterium to kill closely related species with which it might be cohabiting, such bacteriocins have since been shown to be commonplace in the internecine warfare between bacteria. Bacteriocins have been studied primarily from the twin perspectives of how they shape microbial communities and how they penetrate bacteria to kill them. Here, we review the modes of action of a family of bacteriocins that cleave nucleic acid substrates in E. coli, known collectively as nuclease colicins, and the specific immunity (inhibitor) proteins that colicin-producing organisms make in order to avoid committing suicide. In a process akin to targeting in mitochondria, nuclease colicins engage in a variety of cellular associations in order to translocate their cytotoxic domains through the cell envelope to the cytoplasm. As well as informing on the process itself, the study of nuclease colicin import has also illuminated functional aspects of the host proteins they parasitize. We also review recent studies where nuclease colicins and their immunity proteins have been used as model systems for addressing fundamental problems in protein folding and protein-protein interactions, areas of biophysics that are intimately linked to the role of colicins in bacterial competition and to the import process itself.

Sequence and structural features of binding site residues in protein-protein complexes: comparison with protein-nucleic acid complexes

Background

Protein-protein interactions are important for several cellular processes. Understanding the mechanism of protein-protein recognition and predicting the binding sites in protein-protein complexes are long standing goals in molecular and computational biology.

Methods

We have developed an energy based approach for identifying the binding site residues in protein–protein complexes. The binding site residues have been analyzed with sequence and structure based parameters such as binding propensity, neighboring residues in the vicinity of binding sites, conservation score and conformational switching.

Results

We observed that the binding propensities of amino acid residues are specific for protein-protein complexes. Further, typical dipeptides and tripeptides showed high preference for binding, which is unique to protein-protein complexes. Most of the binding site residues are highly conserved among homologous sequences. Our analysis showed that 7% of residues changed their conformations upon protein-protein complex formation and it is 9.2% and 6.6% in the binding and non-binding sites, respectively. Specifically, the residues Glu, Lys, Leu and Ser changed their conformation from coil to helix/strand and from helix to coil/strand. Leu, Ser, Thr and Val prefer to change their conformation from strand to coil/helix.

Conclusions

The results obtained in this study will be helpful for understanding and predicting the binding sites in protein-protein complexes.

Integrating energy calculations with functional assays to decipher the specificity of G protein-RGS protein interactions

The diverse Regulator of G protein Signaling (RGS) family sets the timing of G protein signaling. To understand how the structure of RGS proteins determines their common ability to inactivate G proteins and their selective G protein recognition, we combined structure-based energy calculations with biochemical measurements of RGS activity. We found a previously unidentified group of variable 'Modulatory' residues that reside at the periphery of the RGS domain-G protein interface and fine-tune G protein recognition. Mutations of Modulatory residues in high-activity RGS proteins impaired RGS function, whereas redesign of low-activity RGS proteins in critical Modulatory positions yielded complete gain of function. Therefore, RGS proteins combine a conserved core interface with peripheral Modulatory residues to selectively optimize G protein recognition and inactivation. Finally, we show that our approach can be extended to analyze interaction specificity across other large protein families.

ReplicOpter: a replicate optimizer for flexible docking

We present a computationally efficient method for flexible refinement of docking predictions that reflects observed motions within a protein's structural class. Using structural homologs, we derive deformation models that capture likely motions. The models or "replicates" typically align along a rigid core, with a handful of flexible loops, linkers and tails. A few replicates can generate a much larger number of conformers, by exchanging each flexible region independently of the others. In this way, 10 replicates of a protein having 6 flexible regions can be used to generate a million conformations of a molecule. While this has obvious advantages in terms of sampling, the cost of assessing energies at every conformer is prohibitive, particularly when both molecules are flexible. Our approach addresses this combinatorial explosion, using key assumptions to compress the sampling by many orders of magnitude. ReplicOpter can perform hierarchical clustering from a list of rigid docking predictions and find nearby structures to any promising cluster representatives. These predicted complexes can then be refined and rescored. ReplicOpter's scoring function includes a Lennard-Jones potential softened using the Anderson-Chandler-Weeks decomposition, a desolvation term derived from the Atomic Contact Energy function, Coulombic electrostatics, hydrogen bonding, and terms to model pi-pi and pi-cation interactions. ReplicOpter has performed well on several recent CAPRI systems. We are presently benchmarking ReplicOpter on the complete docking benchmark set to fully establish its utility in refining rigid docking predictions and identifying near-native solutions.

Systematic mutation and thermodynamic analysis of central tyrosine pairs in polyspecific NKG2D receptor interactions

The homodimeric, activating natural killer cell receptor NKG2D interacts with multiple monomeric ligands polyspecifically, yet without central conformational flexibility. Crystal structures of multiple NKG2D-ligand interactions have identified the NKG2D tyrosine pair Tyr 152 and Tyr 199 as forming multiple specific but diverse interactions with MICA and related proteins. Here we systematically altered each tyrosine to tryptophan, phenylalanine, isoleucine, leucine, valine, serine, and alanine to measure the effect of mutation on affinity and thermodynamics for binding a range of similar ligands: MICA, the higher-affinity ligand MICB, and MICdesign, a high-affinity version of MICA that shares all NKG2D contact residues with MICA. Affinity and residue size were related: tryptophan could often substitute for tyrosine without loss of affinity; loss of the tyrosine hydroxyl through mutation to phenylalanine was tolerated more at position 152 than 199; and the smallest residues coincide with lowest affinities in general. NKG2D mutant van't Hoff binding thermodynamics generally show that substitution of other residues for tyrosine causes a moderate positive or flat van't Hoff slope consistent with moderate loss of binding enthalpy. One set of NKG2D mutations caused MICA to adopt a positive van't Hoff slope corresponding to absorption of heat, and another set caused MICB to adopt a negative slope of greater heat release than wild-type. MICdesign shared one example of the first set with MICA and one of the second set with MICB. When the NKG2D mutation affinities were arranged according to change in nonpolar surface area and compared to results from specific antibody-antigen and protein-peptide interactions, it was found that hydrophobic surface loss in NKG2D reduced binding affinity less than reported in the other contexts. The hydrophobic effect at the center of the NKG2D binding appears more similar to that at the periphery of an antibody-antigen binding site than at its center. Therefore the polyspecific NKG2D binding site is more tolerant of structural alteration in general than either an antibody-antigen or protein-peptide binding site, and this tolerance may adapt NKG2D to a broad range of protein surfaces with micromolar affinity.

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.

Energy based approach for understanding the recognition mechanism in protein-protein complexes

Protein-protein interactions play an essential role in the regulation of various cellular processes. Understanding the recognition mechanism of protein-protein complexes is a challenging task in molecular and computational biology. In this work, we have developed an energy based approach for identifying the binding sites and important residues for binding in protein-protein complexes. The new approach is different from the traditional distance based contacts in which the repulsive interactions are treated as binding sites as well as the contacts within a specific cutoff have been treated in the same way. We found that the residues and residue-pairs with charged and aromatic side chains are important for binding. These residues influence to form cation-, electrostatic and aromatic interactions. Our observation has been verified with the experimental binding specificity of protein-protein complexes and found good agreement with experiments. Based on these results we have proposed a novel mechanism for the recognition of protein-protein complexes: the charged and aromatic residues in receptor and ligand initiate recognition by making suitable interactions between them; the neighboring hydrophobic residues assist the stability of complex along with other hydrogen bonding partners by the polar residues. Further, the propensity of residues in the binding sites of receptors and ligands, atomic contributions and the influence on secondary structure will be discussed.

Assembly dynamics and stability of the pneumococcal epsilon zeta antitoxin toxin (PezAT) system from Streptococcus pneumoniae

The pneumococcal epsilon zeta antitoxin toxin (PezAT) system is a chromosomally encoded, class II toxin antitoxin system from the human pathogen Streptococcus pneumnoniae. Neutralization of the bacteriotoxic protein PezT is carried out by complex formation with its cognate antitoxin PezA. Here we study the stability of the inhibitory complex in vivo and in vitro. We found that toxin release is impeded in Escherichia coli and Bacillus subtilis due to the proteolytic resistance of PezA once bound to PezT. These findings are supported by in vitro experiments demonstrating a strong thermodynamic stabilization of both proteins upon binding. A detailed kinetic analysis of PezAT assembly revealed that these particular features of PezAT are based on a strong, electrostatically guided binding mechanism leading to a stable toxin antitoxin complex with femtomolar affinity. Our data show that PezAT complex formation is distinct to all other conventional toxin antitoxin modules and a controlled mode of toxin release is required for activation.

Fundamental link between folding states and functional states of proteins

Folding and function of proteins are two aspects of proteins which are usually considered as basically unrelated phenomena that are optimized by evolution independently. From the funnel model of folding/unfolding and the associated energy landscape, we infer the paradigm that the minimum number of folding intermediates is determined by the number of all functional states of a protein ("essential" folding intermediates). Here, we demonstrate the supposed fundamental link using the Ras protein complexed with the GTP analogue GppNHp that occurs in two structural states coexisting in solution. State 2 was shown earlier to represent the effector interacting state, and the function of state 1 was hitherto unknown. By (31)P NMR spectroscopy, we demonstrate that state 1 represents the conformation interacting with guanine nucleotide exchange factors (GEFs). Denaturation experiments of the protein with a chaotropic reagent show that both functional states coexist during folding and unfolding. Application of high pressure represents another perturbation of the energy landscape, leading to an increased population of the state 1 as observed by NMR spectroscopy. The specific volume difference between the two states DeltaV(12) is 17.2 +/- 0.5 mL mol(-1), indicating that state 1 represents a more open conformation of the protein. The free energies of stabilization for state 1 and state 2 at 278 K can be determined as 8.3 and 9.8 kJ mol(-1), respectively.

Sequence and structural analysis of binding site residues in protein-protein complexes

The binding sites in protein-protein complexes have been identified with different methods including atomic contacts, reduction in solvent accessibility and interaction energy between the interacting partners. In our earlier work, we have developed an energy-based criteria for identifying the binding sites in protein-protein complexes, which showed that the interacting residues are different from that obtained with distance-based methods. In this work, we analyzed the binding site residues based on sequence and structural properties, such as, neighboring residues, secondary structure, solvent accessibility, conservation of residues, medium and long-range contacts and surrounding hydrophobicity. Our results showed that the neighboring residues of binding sites in proteins and ligands are different from each other although the interacting pairs of residues have a common behavior. The analysis on surrounding hydrophobicity reveals that the binding residues are less hydrophobic than non-binding sites, which suggests that the hydrophobic core are important for folding and stability whereas the surface seeking residues play a critical role in binding. This tendency has been verified with the number of contacts in binding sites. In addition, the binding site residues are highly conserved compared with non-binding residues. We suggest that the incorporation of sequence and structure-based features may improve the prediction accuracy of binding sites in protein-protein complexes.

Optimizing protein stability in vivo

Identifying mutations that stabilize proteins is challenging because most substitutions are destabilizing. In addition to being of immense practical utility, the ability to evolve protein stability in vivo may indicate how evolution has formed today's protein sequences. Here we describe a genetic selection that directly links the in vivo stability of proteins to antibiotic resistance. It allows the identification of stabilizing mutations within proteins. The large majority of mutants selected for improved antibiotic resistance are stabilized both thermodynamically and kinetically, indicating that similar principles govern stability in vivo and in vitro. The approach requires no prior structural or functional knowledge and allows selection for stability without a need to maintain function. Mutations that enhance thermodynamic stability of the protein Im7 map overwhelmingly to surface residues involved in binding to colicin E7, showing how the evolutionary pressures that drive Im7-E7 complex formation have compromised the stability of the isolated Im7 protein.

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.

Energetic determinants of protein binding specificity: insights into protein interaction networks

One of the challenges of the postgenomic era is to provide a more realistic representation of cellular processes by combining a systems biology description of functional networks with information on their interacting components. Here we carried out a systematic large-scale computational study on a structural protein-protein interaction network dataset in order to dissect thermodynamic characteristics of binding determining the interplay between protein affinity and specificity. As expected, interactions involving specific binding sites display higher affinities than those of promiscuous binding sites. Next, in order to investigate a possible role of modular distribution of hot spots in binding specificity, we divided binding sites into modules previously shown to be energetically independent. In general, hot spots that interact with different partners are located in different modules. We further observed that common hot spots tend to interact with partners exhibiting common binding motifs, whereas different hot spots tend to interact with partners with different motifs. Thus, energetic properties of binding sites provide insights into the way proteins modulate interactions with different partners. Knowledge of those factors playing a role in protein specificity is important for understanding how proteins acquire additional partners during evolution. It should also be useful in drug design.

Analogous interactions in initiating complexes of the classical and lectin pathways of complement

The classical and lectin pathways of complement activation neutralize pathogens and stimulate key immunological processes. Both pathways are initiated by collagen-containing, soluble pattern recognition molecules associated with specific serine proteases. In the classical pathway, C1q binds to Ab-Ag complexes or bacterial surfaces to activate C1r and C1s. In the lectin pathway, mannan-binding lectin and ficolins bind to carbohydrates on pathogens to activate mannan-binding lectin-associated serine protease 2. To characterize the interactions leading to classical pathway activation, we have analyzed binding between human C1q, C1r, and C1s, which associate to form C1, using full-length and truncated protease components. We show that C1r and C1s bind to C1q independently. The CUB1-epidermal growth factor fragments contribute most toward binding, but CUB2 of C1r, but not of C1s, is also important. Each C1rs tetramer presents a total of six binding sites, one for each of the collagenous domains of C1q. We also demonstrate that subcomponents of the lectin and classical pathways cross-interact. Thus, although the stoichiometries of complexes differ, interactions are analogous, with equivalent contacts between recognition and protease subcomponents. Importantly, these new data are contrary to existing models of C1 and enable us to propose a new model using mannan-binding lectin-mannan-binding lectin-associated serine protease interactions as a template.

The mechanism of folding of Im7 reveals competition between functional and kinetic evolutionary constraints

Many proteins reach their native state through pathways involving the presence of folding intermediates. It is not clear whether this type of folding landscape results from insufficient evolutionary pressure to optimize folding efficiency, or arises from a conflict between functional and folding constraints. Here, using protein-engineering, ultra-rapid mixing and stopped-flow experiments combined with restrained molecular dynamics simulations, we characterize the transition state for the formation of the intermediate populated during the folding of the bacterial immunity protein, Im7, and the subsequent molecular steps leading to the native state. The results provide a comprehensive view of the folding process of this small protein. An analysis of the contributions of native and non-native interactions at different stages of folding reveals how the complexity of the folding landscape arises from concomitant evolutionary pressures for function and folding efficiency.

Chemical and genetic wrappers for improved phage and RNA display

An Achilles heel inherent to all molecular display formats, background binding between target and display system introduces false positives into screens and selections. For example, the negatively charged surfaces of phage, mRNA, and ribosome display systems bind with unacceptably high nonspecificity to positively charged target molecules, which represent an estimated 35% of proteins in the human proteome. Here we report the first systematic attempt to understand why a broad class of molecular display selections fail, and then solve the underlying problem for both phage and RNA display. Firstly, a genetic strategy was used to introduce a short, charge-neutralizing peptide into the solvent-exposed, negatively charged phage coat. The modified phage (KO7(+)) reduced or eliminated nonspecific binding to the problematic high-pI proteins. In the second, chemical approach, nonspecific interactions were blocked by oligolysine wrappers in the cases of phage and total RNA. For phage display applications, the peptides Lys(n) (where n=16 to 24) emerged as optimal for wrapping the phage. Lys(8), however, provided effective wrappers for RNA binding in assays against the RNA binding protein HIV-1 Vif. The oligolysine peptides blocked nonspecific binding to allow successful selections, screens, and assays with five previously unworkable protein targets.

E9-Im9 colicin DNase-immunity protein biomolecular association in water: a multiple-copy and accelerated molecular dynamics simulation study

Protein-protein transient and dynamic interactions underlie all biological processes. The molecular dynamics (MD) of the E9 colicin DNase protein, its Im9 inhibitor protein, and their E9-Im9 recognition complex are investigated by combining multiple-copy (MC) MD and accelerated MD (aMD) explicit-solvent simulation approaches, after validation with crystalline-phase and solution experiments. Im9 shows higher flexibility than its E9 counterpart. Im9 displays a significant reduction of backbone flexibility and a remarkable increase in motional correlation upon E9 association. Im9 loops 23-31 and 54-64 open with respect to the E9-Im9 X-ray structure and show high conformational diversity. Upon association a large fraction (approximately 20 nm2) of E9 and Im9 protein surfaces become inaccessible to water. Numerous salt bridges transiently occurring throughout our six 50 ns long MC-MD simulations are not present in the X-ray model. Among these Im9 Glu31-E9 Arg96 and Im9 Glu41-Lys89 involve interface interactions. Through the use of 10 ns of Im9 aMD simulation, we reconcile the largest thermodynamic impact measured for Asp51Ala mutation with Im9 structure and dynamics. Lys57 acts as an essential molecular switch to shift Im9 surface loop towards an ideal configuration for E9 inhibition. This is achieved by switching Asp60-Lys57 and Asp62-Lys57 hydrogen bonds to Asp51-Lys57 salt bridge. E9-Im9 recognition involves shifts of conformational distributions, reorganization of intramolecular hydrogen bond patterns, and formation of new inter- and intramolecular interactions. The description of key transient biological interactions can be significantly enriched by the dynamic and atomic-level information provided by computer simulations.

Colicins exploit native disorder to gain cell entry: a hitchhiker's guide to translocation

The translocation of protein toxins into a cell relies on a myriad of protein–protein interactions. One such group of toxins are enzymatic E colicins, protein antibiotics produced by Escherichia coli in times of stress. These proteins subvert ordinary nutrient uptake mechanisms to enter the cell and unleash nuclease activity. We, and others, have previously shown that uptake of ColE9 (colicin E9) is dependent on engagement of the OM (outer membrane) receptors BtuB and OmpF as well as recruitment of the periplasmic protein TolB, forming a large supramolecular complex. Intriguingly, colicins bind TolB using a natively disordered region to mimic the interaction of TolB with Pal (peptidoglycan-associated lipoprotein). This is thought to trigger OM instability and prime the system for translocation. Here, we review key interactions in the assembly of this ‘colicin translocon’ and discuss the key role disorder plays in achieving uptake.

Hot-spot residues at the E9/Im9 interface help binding via different mechanisms

Protein-protein association involves many interface interactions, but they do not contribute equally. Ala scanning experiments reveal that only a few mutations significantly lower binding affinity. These key residues, which appear to drive protein-protein association, are called hot-spot residues. Molecular dynamics simulations of the Colicin E9/Im9 complex show Im9 Glu41 and Im9 Ser50, both hot-spots, bind via different mechanisms. The results suggest that Im9 Ser50 restricts Glu41 in a conformation auspicious for salt-bridge formation across the interface. This type of model may be helpful in engineering hot-spot clusters at protein-protein interfaces and, consequently, the design of specificity.

Inactivation of colicin Y by intramembrane helix-helix interaction with its immunity protein

The construction of hybrids between colicins U and Y and the mutagenesis of the colicin Y gene (cya) have revealed amino acid residues important for interactions between colicin Y and its cognate immunity protein (Cyi). Four such residues (I578, T582, Y586 and V590) were found in helices 8 and 9 of the colicin Y pore-forming domain. To verify the importance of these residues, the corresponding amino acids in the colicin B protein were mutated to the residues present in colicin Y. An Escherichia coli strain with cloned colicin Y immunity gene (cyi) inactivated this mutant, but not the wild-type colicin B. In addition, interacting amino acid pairs in Cya and Cyi were identified using a set of Cyi point mutant strains. These data are consistent with antiparallel helix-helix interactions between Cyi helix T3 and Cya helix 8 of the pore-forming domain as a molecular mechanism of colicin Y inactivation by its immunity protein.

Experimental and computational analyses of the energetic basis for dual recognition of immunity proteins by colicin endonucleases

Colicin endonucleases (DNases) are bound and inactivated by immunity (Im) proteins. Im proteins are broadly cross-reactive yet specific inhibitors binding cognate and non-cognate DNases with K(d) values that vary between 10(-4) and 10(-14) M, characteristics that are explained by a 'dual-recognition' mechanism. In this work, we addressed for the first time the energetics of Im protein recognition by colicin DNases through a combination of E9 DNase alanine scanning and double-mutant cycles (DMCs) coupled with kinetic and calorimetric analyses of cognate Im9 and non-cognate Im2 binding, as well as computational analysis of alanine scanning and DMC data. We show that differential DeltaDeltaGs observed for four E9 DNase residues cumulatively distinguish cognate Im9 association from non-cognate Im2 association. E9 DNase Phe86 is the primary specificity hotspot residue in the centre of the interface, which is coordinated by conserved and variable hotspot residues of the cognate Im protein. Experimental DMC analysis reveals that only modest coupling energies to Im9 residues are observed, in agreement with calculated DMCs using the program ROSETTA and consistent with the largely hydrophobic nature of E9 DNase-Im9 specificity contacts. Computed values for the 12 E9 DNase alanine mutants showed reasonable agreement with experimental DeltaDeltaG data, particularly for interactions not mediated by interfacial water molecules. DeltaDeltaG predictions for residues that contact buried water molecules calculated using solvated rotamer models met with mixed success; however, we were able to predict with a high degree of accuracy the location and energetic contribution of one such contact. Our study highlights how colicin DNases are able to utilise both conserved and variable amino acids to distinguish cognate from non-cognate Im proteins, with the energetic contributions of the conserved residues modulated by neighbouring specificity sites.

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.

The molecular architecture of protein-protein binding sites

The formation of specific protein interactions plays a crucial role in most, if not all, biological processes, including signal transduction, cell regulation, the immune response and others. Recent advances in our understanding of the molecular architecture of protein-protein binding sites, which facilitates such diversity in binding affinity and specificity, are enabling us to address key questions. What is the amino acid composition of binding sites? What are interface hotspots? How are binding sites organized? What are the differences between tight and weak interacting complexes? How does water contribute to binding? Can the knowledge gained be translated into protein design? And does a universal code for binding exist, or is it the architecture and chemistry of the interface that enable diverse but specific binding solutions?

Differential recognition of natural and nonnatural substrate by molecular chaperone α-crystallin—A subunit exchange study

alpha-Crystallin is a molecular chaperone that recognizes proteins substrates in stress. It binds to the unstable conformer of a large variety of related or unrelated substrates and thus prevents them aggregating and holds them in a folding competent state. In this article, we have tried to critically analyze, from experimental point of view, whether alpha-crystallin has any preference for its natural substrates compared to the nonnatural one. Our results clearly show that alpha-crystallin is exceptionally active and sensitive in preventing aggregation of its natural substrates and can fully prevent such an aggregation in a substoichiometric ratio, but nonnatural substrates require a considerably higher amount of alpha-crystallin. Using suitable fluorescent-labeled alpha-crystallins and performing fluorescence resonance energy transfer experiments, we were able to determine the subunit exchange kinetics between the alpha-crystallin oligomers. It was found that while alpha-crystallin was bound to its natural substrate, the rate of subunit exchange was slightly decreased. But, when a nonnatural substrate carbonic anhydrase remained bound to the chaperone, further loss in subunit exchange rate was observed. Nonnatural substrate was found to create higher activation energy barrier for the subunit exchange reaction compared to the native substrates. Similarities in major beta-sheet structure of both alpha-crystallin and its natural substrates may be the reason for the preference in molecular recognition in comparison with the nonnatural substrate.

Synthetic anti-BR3 antibodies that mimic BAFF binding and target both human and murine B cells

BR3, which is expressed on all mature B cells, is a specific receptor for the B-cell survival and maturation factor BAFF (B-cell–activating factor belonging to the tumor necrosis factor [TNF] family). In order to investigate the consequences of targeting BR3 in murine models and to assess the potential of BR3 antibodies as human therapeutics, synthetic antibody phage libraries were employed to identify BAFF-blocking antibodies cross-reactive to murine and human BR3, which share 52% identity in their extracellular domains. We found an antibody, CB1, which exhibits μM affinity for murine BR3 and very weak affinity for the human receptor. CB3s, an affinity-matured variant of CB1, has sub-nM affinity for BR3 from both species. Alanine scanning and crystallographic structural analysis of the CB3s/BR3 complex reveal that CB3s mimics BAFF by interacting with a similar region of the BR3 surface. Despite this similarity in binding epitopes, CB1 variants antagonize BAFF-dependent human B-cell proliferation in vitro and are effective at reducing murine B-cell populations in vivo, showing significant promise as therapeutics for human B-cell–mediated diseases.

Recognition of pore-forming colicin Y by its cognate immunity protein

Construction of hybrid immunity genes between colicin U (cui) and Y (cyi) immunity genes and site-directed mutagenesis of cyi were used to identify amino-acid residues of the colicin Y immunity protein (Cyi) involved in recognition of colicin Y. These amino-acid residues were localized close to the cytoplasmic site of the Cyi transmembrane helices T3 (S104, S107, F110, A112) and T4 (A159). Mutations in cui, which converted Cui sequence to Cyi sequence in positions 104, 107, 110, 112 and 159, resulted in an immunity gene that also conferred (besides immunity to colicin U) a high degree of immunity to colicin Y.