Publications

Janus kinases (JAKs) bind to class I and II cytokine receptors, activating signaling and regulating gene transcription through signal transducer and activator of transcription (STAT) proteins. Type I interferons (IFNs) require the JAK members TYK2 and JAK1, which bind to the receptor subunits IFNAR1 and IFNAR2, respectively. We investigated the role of JAKs in regulating IFNAR signaling activity. Synthetic IFNARs in which the extracellular domains of IFNAR1 and IFNAR2 are replaced with nanobodies had near-native type I IFN signaling, whereas the homomeric variant of IFNAR2 initiated much weaker signaling, despite harboring docking sites for JAKs and STATs. Cells with JAK1 and TYK2 knockout (KO) showed residual signaling, suggesting partial complementation by the remaining JAKs, particularly when they were overexpressed. Live-cell micropatterning experiments confirmed the promiscuous binding of JAK1, JAK2, and TYK2 to IFNAR1 and IFNAR2, and their recruitment correlated with their relative cellular abundances. However, each JAK had a different efficacy in inducing cross-phosphorylation and downstream signaling. JAK binding was also promiscuous for other cytokine receptors, including IFN-L1, IL-10Rβ, TPOR, and GHR, but not for EPOR, which activated different downstream signaling pathways. These findings suggest that competitive binding of JAKs to cytokine receptors together with the varying absolute and relative abundances of the JAKs in different cell types can account for the cell typedependent signaling pleiotropy of cytokine receptors.

SARS-CoV-2 has the capacity to evolve mutations that escape vaccine- and infection-acquired immunity and antiviral drugs. A variant-agnostic therapeutic agent that protects against severe disease without putting selective pressure on the virus would thus be a valuable biomedical tool that would maintain its efficacy despite the ongoing emergence of new variants. Here, we challenge male rhesus macaques with SARS-CoV-2 Deltathe most pathogenic variant in a highly susceptible animal model. At the time of challenge, we also treat the macaques with aerosolized RBD-62, a protein developed through multiple rounds of in vitro evolution of SARS-CoV-2 RBD to acquire 1000-fold enhanced ACE2 binding affinity. RBD-62 treatment equivalently suppresses virus replication in both upper and lower airways, a phenomenon not previously observed with clinically approved vaccines. Importantly, RBD-62 does not block the development of virus-specific T- and B-cell responses and does not elicit anti-drug immunity. These data provide proof-of-concept that RBD-62 can prevent severe disease from a highly virulent variant.

Interferons (IFNs) are essential for defense against viral infections but also drive recruitment of inflammatory cells to sites of infection, a key feature of severe COVID-19. Here, we explore the complexity of the IFN response in COVID-19, examine the effects of manipulating IFN on SARS-CoV-2 viral replication and pathogenesis, and highlight pre-clinical and clinical studies evaluating the therapeutic efficacy of IFN in limiting COVID-19 severity.

Type I and type III interferons (IFNs) are major components in activating the innate immune response. Common to both are two distinct receptor chains (IFNAR1/IFNAR2 and IFNLR1/IL10R2), which form ternary complexes upon binding their respective ligands. This results in close proximity of the intracellularly associated kinases JAK1 and TYK2, which cross phosphorylate each other, the associated receptor chains, and signal transducer and activator of transcriptions, with the latter activating IFN-stimulated genes. While there are clear similarities in the biological responses toward type I and type III IFNs, differences have been found in their tropism, tuning of activity, and induction of the immune response. Here, we focus on how these differences are embedded in the structure/function relations of these two systems in light of the recent progress that provides in-depth information on the structural assembly of these receptors and their functional implications and how these differ between the mouse and human systems.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has substantially diversified during the pandemic, resulting in the successive emergence of variants characterized by various mutations. It has been observed that several epidemic variants, including those classified as variants of concern, share mutations at four key residues (L452R, T478K, E484K, and N501Y) within the receptor binding motif (RBM) region of the spike protein. However, the processes through which these four specific RBM mutations were acquired during the evolution of SARS-CoV-2, as well as the degree to which they enhance viral fitness, remain unclear. Moreover, the effect of these mutations on the properties of the spike protein is not yet fully understood. In this study, we performed a comprehensive phylogenetic analysis and showed that the four RBM mutations have been convergently acquired across various lineages throughout the evolutionary history of SARS-CoV-2. We also found a specific pattern in the order of acquisition for some of these mutations. Additionally, our epidemic dynamic modeling demonstrated that acquiring these mutations leads to an increase in the effective reproduction number of the virus. Furthermore, we engineered mutant spike proteins with all feasible combinations of the four mutations, and examined their properties to uncover the influence that these mutations have on viral characteristics. Our results provide insights into the roles these four mutations play in shaping the viral characteristics, epidemic proliferation, and evolutionary pathway of SARS-CoV-2.

In late 2022, various Omicron subvariants emerged and cocirculated worldwide. These variants convergently acquired amino acid substitutions at critical residues in the spike protein, including residues R346, K444, L452, N460, and F486. Here, we characterize the convergent evolution of Omicron subvariants and the properties of one recent lineage of concern, BQ.1.1. Our phylogenetic analysis suggests that these five substitutions are recurrently acquired, particularly in younger Omicron lineages. Epidemic dynamics modelling suggests that the five substitutions increase viral fitness, and a large proportion of the fitness variation within Omicron lineages can be explained by these substitutions. Compared to BA.5, BQ.1.1 evades breakthrough BA.2 and BA.5 infection sera more efficiently, as demonstrated by neutralization assays. The pathogenicity of BQ.1.1 in hamsters is lower than that of BA.5. Our multiscale investigations illuminate the evolutionary rules governing the convergent evolution for known Omicron lineages as of 2022.

In late 2022, SARS-CoV-2 Omicron subvariants have become highly diversified, and XBB is spreading rapidly around the world. Our phylogenetic analyses suggested that XBB emerged through the recombination of two cocirculating BA.2 lineages, BJ.1 and BM.1.1.1 (a progeny of BA.2.75), during the summer of 2022. XBB.1 is the variant most profoundly resistant to BA.2/5 breakthrough infection sera to date and is more fusogenic than BA.2.75. The recombination breakpoint is located in the receptor-binding domain of spike, and each region of the recombinant spike confers immune evasion and increases fusogenicity. We further provide the structural basis for the interaction between XBB.1 spike and human ACE2. Finally, the intrinsic pathogenicity of XBB.1 in male hamsters is comparable to or even lower than that of BA.2.75. Our multiscale investigation provides evidence suggesting that XBB is the first observed SARS-CoV-2 variant to increase its fitness through recombination rather than substitutions.

Previous studies on the Omicron BA.2 variant suggested that the virological characteristics of BA.2 are determined by the mutations in at least two different regions of the viral genome: in the BA.2 spike gene (enhancing viral fusogenicity and intrinsic pathogenicity) and the non-spike region of the BA.2 genome (leading to intrinsic pathogenicity attenuation). However, the mutations modulating the BA.2 virological properties remain elusive. In this study, we demonstrated that the L371F substitution in the BA.2 spike protein confers greater fusogenicity and intrinsic pathogenicity. Furthermore, we revealed that multiple mutations downstream of the spike gene in the BA.2 genome are responsible for attenuating intrinsic viral pathogenicity and replication capacity. As mutations in the SARS-CoV-2 variant spike proteins could modulate certain virological properties, such as immune evasion and infectivity, most studies have previously focused on spike protein mutations. Our results underpin the importance of non-spike protein-related mutations in SARS-CoV-2 variants. IMPORTANCE Most studies investigating the characteristics of emerging SARS-CoV-2 variants have been focusing on mutations in the spike proteins that affect viral infectivity, fusogenicity, and pathogenicity. However, few studies have addressed how naturally occurring mutations in the non-spike regions of the SARS-CoV-2 genome impact virological properties. In this study, we proved that multiple SARS-CoV-2 Omicron BA.2 mutations, one in the spike protein and another downstream of the spike gene, orchestrally characterize this variant, shedding light on the importance of Omicron BA.2 mutations out of the spike protein.

Type I interferons (IFN-I) are critical mediators of innate control of viral infections but also drive the recruitment of inflammatory cells to sites of infection, a key feature of severe coronavirus disease 2019. Here, IFN-I signaling was modulated in rhesus macaques (RMs) before and during acute SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) infection using a mutated IFN-α2 (IFN-modulator; IFNmod), which has previously been shown to reduce the binding and signaling of endogenous IFN-I. IFNmod treatment in uninfected RMs was observed to induce a modest up-regulation of only antiviral IFN-stimulated genes (ISGs); however, in SARS-CoV-2- infected RMs, IFNmod reduced both antiviral and inflammatory ISGs. IFNmod treatment resulted in a potent reduction in SARS-CoV-2 viral loads both in vitro in Calu-3 cells and in vivo in bronchoalveolar lavage (BAL), upper airways, lung, and hilar lymph nodes of RMs. Furthermore, in SARS-CoV-2-infected RMs, IFNmod treatment potently reduced inflammatory cytokines, chemokines, and CD163+ MRC1- inflammatory macrophages in BAL and expression of Siglec-1 on circulating monocytes. In the lung, IFNmod also reduced pathogenesis and attenuated pathways of inflammasome activation and stress response during acute SARS-CoV-2 infection. Using an intervention targeting both IFN-α and IFN-β pathways, this study shows that, whereas early IFN-I restrains SARS-CoV-2 replication, uncontrolled IFN-I signaling critically contributes to SARS-CoV-2 inflammation and pathogenesis in the moderate disease model of RMs.

Introduction: Chronic lymphocytic leukemia (CLL) is characterized by an aberrant cytokine network that can support tumor growth by triggering janus kinase (JAK)/STAT pathways. Targeting cytokine-signaling should then be a rational therapeutic strategy but the JAK inhibitor ruxolitinib failed to control and seemingly accelerated the disease in clinical trials. Methods: The effect of ruxolitinib on primary human CLL cells was studied in vitro and in vivo. Results: Ruxolitinib increased phosphorylation of IRAK4, an important toll-like receptor (TLR)- signaling intermediate, in circulating CLL cells in vitro. It also enhanced p38 and NFKB1 phosphorylation while lowering STAT3 phosphorylation in CLL cells activated with TLR-7/8 agonists and IL-2. Among the cytokines made by activated CLL cells, high levels of IL-10 contributed strongly to STAT3 phosphorylation and inhibited TLR7 activity. Ruxolitinib limited TLR-mediated IL10 transcription and markedly reduced IL-10 production in vitro. It also decreased blood levels of IL-10 while increasing TNFα along with phospho-p38 expression and gene sets associated with TLR-activation in CLL cells in vivo. The bruton's tyrosine kinase inhibitor ibrutinib decreased IL-10 production in vitro but, in contrast to ruxolitinib, blocked initial IL10 transcription induced by TLR-signaling in vitro, decreased TNFα production, and deactivates CLL cells in vivo. Discussion: These findings suggest the possible benefits of inhibiting growth factors with JAK inhibitors in CLL are outweighed by negative effects on potential tumor suppressors such as IL-10 that allow unrestrained activation of NFκB by drivers such as TLRs. Specific inhibition of growth-promoting cytokines with blocking antibodies or infusing suppressive cytokines like IL-10 might be better strategies to manipulate cytokines in CLL.

Recent studies have revealed the unique virological characteristics of Omicron, particularly those of its spike protein, such as less cleavage efficacy in cells, reduced ACE2 binding affinity, and poor fusogenicity. However, it remains unclear which mutation(s) determine these three virological characteristics of Omicron spike. Here, we show that these characteristics of the Omicron spike protein are determined by its receptor-binding domain. Of interest, molecular phylogenetic analysis revealed that acquisition of the spike S375F mutation was closely associated with the explosive spread of Omicron in the human population. We further elucidated that the F375 residue forms an interprotomer pi-pi interaction with the H505 residue of another protomer in the spike trimer, conferring the attenuated cleavage efficiency and fusogenicity of Omicron spike. Our data shed light on the evolutionary events underlying the emergence of Omicron at the molecular level.

The SARS-CoV-2 Omicron BA.2.75 variant emerged in May 2022. BA.2.75 is a BA.2 descendant but is phylogenetically distinct from BA.5, the currently predominant BA.2 descendant. Here, we show that BA.2.75 has a greater effective reproduction number and different immunogenicity profile than BA.5. We determined the sensitivity of BA.2.75 to vaccinee and convalescent sera as well as a panel of clinically available antiviral drugs and antibodies. Antiviral drugs largely retained potency, but antibody sensitivity varied depending on several key BA.2.75-specific substitutions. The BA.2.75 spike exhibited a profoundly higher affinity for its human receptor, ACE2. Additionally, the fusogenicity, growth efficiency in human alveolar epithelial cells, and intrinsic pathogenicity in hamsters of BA.2.75 were greater than those of BA.2. Our multilevel investigations suggest that BA.2.75 acquired virological properties independent of BA.5, and the potential risk of BA.2.75 to global health is greater than that of BA.5.

Crowded environments are known to affect the diffusion of macromolecules, but their effects on the diffusion of small molecules are largely uncharacterized. We investigate how three protein crowders, bovine serum albumin (BSA), hen egg -white lysozyme, and myoglobin, influence the diffusion rates and interactions of four small molecules: fluorescein, and three drugs, doxorubicin, glycogen syn-thase kinase-3 inhibitor SB216763, and quinacrine. Using Line-FRAP measure-ments, Brownian dynamics simulations, and molecular docking, we find that the diffusion rates of the small molecules are highly affected by self-aggregation, in-teractions with the proteins, and surface adsorption. The diffusion of fluorescein is decreased because of its interactions with the protein crowders and their sur-face adsorption. Protein crowders increase the diffusion rates of doxorubicin and SB216763 by reducing surface interactions and self-aggregation, respectively. Quinacrine diffusion was not affected by protein crowders. The mechanistic in-sights gained here may assist in optimization of compounds for higher mobility in complex macromolecular environments.

After the global spread of the SARS-CoV-2 Omicron BA.2, some BA.2 subvariants, including BA.2.9.1, BA.2.11, BA.2.12.1, BA.4, and BA.5, emerged in multiple countries. Our statistical analysis showed that the effective reproduction numbers of these BA.2 subvariants are greater than that of the original BA.2. Neutralization experiments revealed that the immunity induced by BA.1/2 infections is less effective against BA.4/5. Cell culture experiments showed that BA.2.12.1 and BA.4/5 replicate more efficiently in human alveolar epithelial cells than BA.2, and particularly, BA.4/5 is more fusogenic than BA.2. We further provided the structure of the BA.4/5 spike receptor-binding domain that binds to human ACE2 and considered how the substitutions in the BA.4/5 spike play roles in ACE2 binding and immune evasion. Moreover, experiments using hamsters suggested that BA.4/5 is more pathogenic than BA.2. Our multiscale investigations suggest that the risk of BA.2 subvariants, particularly BA.4/5, to global health is greater than that of original BA.2.

Over half the proteins in the E. coli cytoplasm form homo or hetero-oligomeric structures. Experimentally determined structures are often considered in determining a protein's oligomeric state, but static structures miss the dynamic equilibrium between different quaternary forms. The problem is exacerbated in homo-oligomers, where the oligomeric states are challenging to characterize. Here, we re-evaluated the oligomeric state of 17 different bacterial proteins across a broad range of protein concentrations and solutions by native mass spectrometry (MS), mass photometry (MP), size exclusion chromatography (SEC), and small-angle X-ray scattering (SAXS), finding that most exhibit several oligomeric states. Surprisingly, some proteins did not show mass-action driven equilibrium between the oligomeric states. For approximately half the proteins, the predicted oligomeric forms described in publicly available databases underestimated the complexity of protein quaternary structures in solution. Conversely, AlphaFold multimer provided an accurate description of the potential multimeric states for most proteins, suggesting that it could help resolve uncertainties on the solution state of many proteins.

The elucidation of viral-receptor interactions and an understanding of virus-spreading mechanisms are of great importance, particularly in the era of a pandemic. Indeed, advances in computational chemistry, synthetic biology, and protein engineering have allowed precise prediction and characterization of such interactions. Nevertheless, the hazards of the infectiousness of viruses, their rapid mutagenesis, and the need to study viral-receptor interactions in a complex in vivo setup call for further developments. Here, we show the development of biocompatible genetically engineered extracellular vesicles (EVs) that display the receptor binding domain (RBD) of SARS-CoV-2 on their surface as coronavirus mimetics (EVsRBD). Loading EVsRBD with iron oxide nanoparticles makes them MRI-visible and, thus, allows mapping of the binding of RBD to ACE2 receptors noninvasively in live subjects. Moreover, we show that EVsRBD can be modified to display mutants of the RBD of SARS-CoV-2, allowing rapid screening of currently raised or predicted variants of the virus. The proposed platform thus shows relevance and cruciality in the examination of quickly evolving pathogenic viruses in an adjustable, fast, and safe manner. Relying on MRI for visualization, the presented approach could be considered in the future to map ligand-receptor binding events in deep tissues, which are not accessible to luminescence-based imaging.

Natural evolution is driven by random mutations that improve fitness. In vitro evolution mimics this process, however, on a short time-scale and is driven by the given bait. Here, we used directed in vitro evolution of a random mutant library of Uracil glycosylase (eUNG) displayed on yeast surface to select for binding to chaperones GroEL, DnaK + DnaJ + ATP (DnaKJ) or E. coli cell extract (CE), using binding to the eUNG inhibitor Ugi as probe for native fold. The CE selected population was further divided to Ugi binders (+U) or non-binders (−U). The aim here was to evaluate the sequence space and physical state of the evolved protein binding the different baits. We found that GroEL, DnaKJ and CE-U select and enrich for mutations causing eUNG to misfold, with the three being enriched in mutations in buried and conserved positions, with a tendency to increase positive charge. Still, each selection had its own trajectory, with GroEL and CE-U selecting mutants highly sensitive to protease cleavage while DnaKJ selected partially structured misfolded species with a tendency to refold, making them less sensitive to proteases. More general, our results show that GroEL has a higher tendency to purge promiscuous misfolded protein mutants from the system, while DnaKJ binds misfolding-prone mutant species that are, upon chaperone release, more likely to natively refold. CE-U shares some of the properties of GroEL- and DnaKJ-selected populations, while harboring also unique properties that can be explained by the presence of additional chaperones in CE, such as Trigger factor, HtpG and ClpB.

Soon after the emergence and global spread of the SARS-CoV-2 Omicron lineage BA.1, another Omicron lineage, BA.2, began outcompeting BA.1. The results of statistical analysis showed that the effective reproduction number of BA.2 is 1.4-fold higher than that of BA.1. Neutralization experiments revealed that immunity induced by COVID vaccines widely administered to human populations is not effective against BA.2, similar to BA.1, and that the antigenicity of BA.2 is notably different from that of BA.1. Cell culture experiments showed that the BA.2 spike confers higher replication efficacy in human nasal epithelial cells and is more efficient in mediating syncytia formation than the BA.1 spike. Furthermore, infection experiments using hamsters indicated that the BA.2 spike-bearing virus is more pathogenic than the BA.1 spike-bearing virus. Altogether, the results of our multiscale investigations suggest that the risk of BA.2 to global health is potentially higher than that of BA.1.

Delivering medication to the lungs via nebulization of pharmaceuticals is a noninvasive and efficient therapy route, particularly for respiratory diseases. The recent worldwide severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) pandemic urges the development of such therapies as an effective alternative to vaccines. The main difficulties in using inhalation therapy are the development of effective medicine and methods to stabilize the biological molecules and transfer them to the lungs efficiently following nebulization. We have developed a high-affinity angiotensin-converting enzyme 2 (ACE2) receptor-binding domain (RBD-62) that can be used as a medication to inhibit infection with SARS-CoV-2 and its variants. In this study, we established a nebulization protocol for drug delivery by inhalation using two commercial vibrating mesh (VM) nebulizers (Aerogen Solo and PARI eFlow) that generate similar mist size distribution in a size range that allows efficient deposition in the small respiratory airway. In a series of experiments, we show the high activity of RBD-62, interferon-α2 (IFN-α2), and other proteins following nebulization. The addition of gelatin significantly stabilizes the proteins and enhances the fractions of active proteins after nebulization, minimizing the medication dosage. Furthermore, hamster inhalation experiments verified the feasibility of the protocol in pulmonary drug delivery. In short, the gelatin-modified RBD-62 formulation in coordination with VM nebulizer can be used as a therapy to cure SARS-CoV-2.

Viruses rapidly co-evolve with their hosts. The 9 million sequenced SARS-CoV-2 genomes by March 2022 provide a detailed account of viral evolution, showing that all amino acids have been mutated many times. However, only a few became prominent in the viral population. Here, we investigated the emergence of the same mutations in unrelated parallel lineages and the extent of such convergent evolution on the molecular level in the spike (S) protein. We found that during the first phase of the pandemic (until mid 2021, before mass vaccination) 31 mutations evolved independently ≥3-times within separated lineages. These included all the key mutations in SARS-CoV-2 variants of concern (VOC) at that time, indicating their fundamental adaptive advantage. The omicron added many more mutations not frequently seen before, which can be attributed to the synergistic nature of these mutations, which is more difficult to evolve. The great majority (24/31) of S-protein mutations under convergent evolution tightly cluster in three functional domains; N-terminal domain, receptor-binding domain, and Furin cleavage site. Furthermore, among the S-protein receptor-binding motif mutations, ACE2 affinity-improving substitutions are favoured. Next, we determined the mutation space in the S protein that has been covered by SARS-CoV-2. We found that all amino acids that are reachable by single nucleotide changes have been probed multiple times in early 2021. The substitutions requiring two nucleotide changes have recently (late 2021) gained momentum and their numbers are increasing rapidly. These provide a large mutation landscape for SARS-CoV-2 future evolution, on which research should focus now.

We hereby describe the process of design and selection of nonantibody protein binders mimicking cytokine signaling. We chose to mimic signaling of IFN-λ1, type 3 interferon (also known as IL-29) for its novelty and the importance of its biological functions. All four known interferons λ signal through binding to the extracellular domains of IL-28 receptor 1 (IL-28R1) and IL-10 receptor 2 (IL-10R2). Our binders were therefore trained to bind both receptors simultaneously. The bifunctional binder molecules were developed by yeast display, a method of directed evolution. The signaling capacity of the bivalent binders was tested by measuring phosphorylation of the JAK/STAT signaling pathway and production of mRNA of six selected genes naturally induced by IFN- λ1 in human cell lines. The newly developed bivalent binders offer opportunities to study cytokine-related biological functions and modulation of the cell behavior by receptor activation on the cell surfaces alternative to the use of natural IFN-λ.

On 24^th November 2021, the sequence of a new SARS-CoV-2 viral isolate Omicron-B.1.1.529 was announced, containing far more mutations in Spike (S) than previously reported variants. Neutralization titers of Omicron by sera from vaccinees and convalescent subjects infected with early pandemic Alpha, Beta, Gamma, or Delta are substantially reduced, or the sera failed to neutralize. Titers against Omicron are boosted by third vaccine doses and are high in both vaccinated individuals and those infected by Delta. Mutations in Omicron knock out or substantially reduce neutralization by most of the large panel of potent monoclonal antibodies and antibodies under commercial development. Omicron S has structural changes from earlier viruses and uses mutations that confer tight binding to ACE2 to unleash evolution driven by immune escape. This leads to a large number of mutations in the ACE2 binding site and rebalances receptor affinity to that of earlier pandemic viruses.

SARS-CoV-2 Lambda, a variant of interest, has spread in some South American countries; however, its virological features and evolutionary traits remain unclear. In this study, we use pseudoviruses and reveal that the spike protein of the Lambda variant is more infectious than that of other variants due to the T76I and L452Q mutations. The RSYLTPGD246-253N mutation, a unique 7-amino acid deletion in the N-terminal domain of the Lambda spike protein, is responsible for evasion from neutralizing antibodies and further augments antibody-mediated enhancement of infection. Although this mutation generates a nascent N-linked glycosylation site, the additional N-linked glycan is dispensable for the virological property conferred by this mutation. Since the Lambda variant has dominantly spread according to the increasing frequency of the isolates harboring the RSYLTPGD246-253N mutation, our data suggest that the RSYLTPGD246-253N mutation is closely associated with the substantial spread of the Lambda variant in South America.

Here, we enhanced the popular yeast display method by multiple rounds of DNA and protein engineering. We introduced surface exposure-tailored reporters, eUnaG2 and DnbALFA, creating a new platform of C and N terminal fusion vectors. The optimization of eUnaG2 resulted in five times brighter fluorescence and 10 °C increased thermostability than UnaG. The optimized DnbALFA has 10-fold the level of expression of the starting protein. Following this, different plasmids were developed to create a complex platform allowing a broad range of protein expression organizations and labeling strategies. Our platform showed up to five times better separation between nonexpressing and expressing cells compared with traditional pCTcon2 and c-myc labeling, allowing for fewer rounds of selection and achieving higher binding affinities. Testing 16 different proteins, the enhanced system showed consistently stronger expression signals over c-myc labeling. In addition to gains in simplicity, speed, and cost-effectiveness, new applications were introduced to monitor protein surface exposure and protein retention in the secretion pathway that enabled successful protein engineering of hard-to-express proteins. As an example, we show how we optimized the WD40 domain of the ATG16L1 protein for yeast surface and soluble bacterial expression, starting from a nonexpressing protein. As a second example, we show how using the here-presented enhanced yeast display method we rapidly selected high-affinity binders toward two protein targets, demonstrating the simplicity of generating new proteinprotein interactions. While the methodological changes are incremental, it results in a qualitative enhancement in the applicability of yeast display for many applications.

The formation of specific protein-protein interactions (PPIs) drive most biological processes. Malfunction of such interactions is the molecular driver of many diseases. Our ability to engineer existing PPIs or create new ones has become a vital research tool. In addition, engineered proteins with new or altered interactions are among the most critical drugs that have been developed in recent years. These include antibodies, cytokines, inhibitors, and others. Here, we provide a perspective on the current status of the methods used to engineer new or altered PPIs. The emergence of the COVID-19 pandemic, which resulted in a worldwide quest to develop specific PPI inhibitors as drugs, provided an up-to-date and state-of-the-art status report on the methodologies for engineering PPIs targeting the interaction of the viral spike protein with its cellular target, ACE2. Multiple, very high affinity binders were generated within a few months using in vitro evolution by itself, or in combination with computational design. The different experimental and computational methods used to block this interaction provide a road map for the future of PPI engineering.

Type I interferons bind to cell surface receptors composed of the subunits IFNAR1 and IFNAR2, the intracellular domains (ICDs) of which are associated with the kinases TYK2 and JAK1, respectively. Ligand binding results in the cross-phosphorylation of TYK2 and JAK1, which then phosphorylate tyrosine residues in the ICDs of the receptor subunits and members of the STAT family of transcription factors. The phosphorylated STATs migrate to the nucleus and drive transcription. We analyzed receptor mutants in knockout cells to study the functional importance of various regions of the receptor ICDs. For IFNAR1, only the TYK2 binding site in the ICD was required for signaling. In contrast, successive truncations of the ICD of IFNAR2 proportionally decreased constitutive STAT binding, STAT phosphorylation, and target gene activation. These findings fit a model in which nonsuccessive stretches along the ICD interact with STATs. Tyrosine residues in the IFNAR1 ICD were not required for signaling, and single tyrosine mutations in the IFNAR2 ICD did not affect signal activation. However, simultaneous mutation of all the tyrosine residues in IFNAR2-ICD reduced STAT phosphorylation, STAT-mediated transcriptional activation, and antiviral activity but not constitutive STAT2 binding. We suggest that tyrosine phosphorylation on IFNAR2-ICD drives the dissociation of phosphorylated STATs, thus maintaining high signaling flux.

SARS-CoV-2 variants of interest and concern will continue to emerge for the duration of the COVID-19 pandemic. To map mutations in the receptor-binding domain (RBD) of the spike protein that affect binding to angiotensin-converting enzyme 2 (ACE2), the receptor for SARS-CoV-2, we applied in vitro evolution to affinity-mature the RBD. Multiple rounds of random mutagenic libraries of the RBD were sorted against decreasing concentrations of ACE2, resulting in the selection of higher affinity RBD binders. We found that mutations present in more transmissible viruses (S477N, E484K and N501Y) were preferentially selected in our high-throughput screen. Evolved RBD mutants include prominently the amino acid substitutions found in the RBDs of B.1.620, B.1.1.7 (Alpha), B1.351 (Beta) and P.1 (Gamma) variants. Moreover, the incidence of RBD mutations in the population as presented in the GISAID database (April 2021) is positively correlated with increased binding affinity to ACE2. Further in vitro evolution increased binding by 1,000-fold and identified mutations that may be more infectious if they evolve in the circulating viral population, for example, Q498R is epistatic to N501Y. We show that our high-affinity variant RBD-62 can be used as a drug to inhibit infection with SARS-CoV-2 and variants Alpha, Beta and Gamma in vitro. In a model of SARS-CoV-2 challenge in hamster, RBD-62 significantly reduced clinical disease when administered before or after infection. A 2.9 Å cryo-electron microscopy structure of the high-affinity complex of RBD-62 and ACE2, including all rapidly spreading mutations, provides a structural basis for future drug and vaccine development and for in silico evaluation of known antibodies.

Type I Interferons (IFN-Is) are a family of cytokines which play a major role in inhibiting viral infection. Resultantly, many viruses have evolved mechanisms in which to evade the IFN-I response. Here we tested the impact of expression of 27 different SARS-CoV-2 genes in relation to their effect on IFN production and activity using three independent experimental methods. We identified six gene products; NSP6, ORF6, ORF7b, NSP1, NSP5 and NSP15, which strongly (>10-fold) blocked MAVS-induced (but not TRIF-induced) IFNβ production. Expression of the first three of these SARS-CoV-2 genes specifically blocked MAVSinduced IFNβ-promoter activity, whereas all six genes induced a collapse in IFNβ mRNA levels, corresponding with suppressed IFNβ protein secretion. Five of these six genes furthermore suppressed MAVS-induced activation of IFNλs, however with no effect on IFNα or IFNγ production. In sharp contrast, SARS-CoV-2 infected cells remained extremely sensitive to anti-viral activity exerted by added IFN-Is. None of the SARS-CoV-2 genes were able to block IFN-I signaling, as demonstrated by robust activation of Interferon Stimulated Genes (ISGs) by added interferon. This, despite the reduced levels of STAT1 and phospho- STAT1, was likely caused by broad translation inhibition mediated by NSP1. Finally, we found that a truncated ORF7b variant that has arisen from a mutant SARS-CoV-2 strain harboring a 382-nucleotide deletion associating with mild disease (Δ382 strain identified in Singapore & Taiwan in 2020) lost its ability to suppress type I and type III IFN production. In summary, our findings support a multi-gene process in which SARS-CoV-2 blocks IFN-production, with ORF7b as a major player, presumably facilitating evasion of host detection during early infection. However, SARS-CoV-2 fails to suppress IFN-I signaling thus providing an opportunity to exploit IFN-Is as potential therapeutic antiviral drugs.

Many SARS-CoV-2 variants with naturally acquired mutations have emerged. These mutations can affect viral properties such as infectivity and immune resistance. Although the sensitivity of naturally occurring SARS-CoV-2 variants to humoral immunity has been investigated, sensitivity to human leukocyte antigen (HLA)-restricted cellular immunity remains largely unexplored. Here, we demonstrate that two recently emerging mutations in the receptor-binding domain of the SARS-CoV-2 spike protein, L452R (in B.1.427/429 and B.1.617) and Y453F (in B.1.1.298), confer escape from HLA-A24-restricted cellular immunity. These mutations reinforce affinity toward the host entry receptor ACE2. Notably, the L452R mutation increases spike stability, viral infectivity, viral fusogenicity, and thereby promotes viral replication. These data suggest that HLA-restricted cellular immunity potentially affects the evolution of viral phenotypes and that a further threat of the SARS-CoV-2 pandemic is escape from cellular immunity.

The crowded cellular milieu affect molecular diffusion through hard (occluded space) and soft (weak, non-specific) interactions. Multiple methods have been developed to measure diffusion coefficients at physiological protein concentrations within cells, each with its limitations. Here, we show that Line-FRAP, combined with rigours data analysis, is able to determine diffusion coefficients in a variety of environments, from in vitro to in vivo. The use of Line mode greatly improves time resolution of FRAP data acquisition, from 20-100 Hz in the classical mode to 800 Hz in the line mode. This improves data analysis, as intensity and radius of the bleach at the first post-bleach frame is critical. We evaluated the method on different proteins labelled chemically or fused to YFP in a wide range of environments. The diffusion coefficients measured in HeLa and in E. coli were ~2.5-fold and 15-fold slower than in buffer, and were comparable to previously published data. Increasing the osmotic pressure on E. coli further decreases diffusion, to the point at which proteins virtually stop moving. The method presented here, which requires a confocal microscope equipped with dual scanners, can be applied to study a large range of molecules with different sizes, and provides robust results in a wide range of environments and protein concentrations for fast diffusing molecules.

Proteins must work together with other proteins to carry out most of their functions in the cell. In the complex biological environment, any correct interaction competes with a huge number of non-relevant macromolecular surfaces. In this chapter, I review the current knowledge on how the binding process occurs, what distinguishes correct binding from the endless numbers of available protein-surfaces and the contributing factors towards uniqueness of the binding interface, and put these into perspective with the different types of protein-protein interactions found in nature. While the gained knowledge provides only averages, it is sufficiently good for the design of new protein-protein interactions, which has become more successful in recent years.

Type I interferons (IFNs) have a central role in innate and adaptive immunities, proliferation, and cancer surveillance. How IFN binding to its specific receptor, the IFN α and β receptor (IFNAR), can drive such variety of processes is an open question. Here, to systematically and thoroughly investigate the molecular mechanism of IFN signaling, we used a CRISPR/Cas9-based approach in a human cell line (HeLa) to generate knockouts (KOs) of the genes participating in the type 1 IFN signaling cascade. We show that both IFNAR chains (IFNAR1 and IFNAR2) are absolutely required for any IFN-induced signaling. Deletion of either signal transducer and activator of transcription 1 (STAT1) or STAT2 had only a partial effect on IFN-induced antiviral activity or gene induction. However, the deletion of both genes completely abrogated any IFN-induced activity. So did a double STAT2IFN regulatory factor 1 (IRF1) KO and, to a large extent, a STAT1 KO together with IRF9 knockdown. KO of any of the STATs had no effect on the phosphorylation of other STATs, indicating that they bound IFNAR independently. STAT3 and STAT6 phosphorylations were fully induced by type 1 IFN in the STAT1STAT2 KO, but did not promote gene induction. Moreover, STAT3 KO did not affect type 1 IFN-induced gene or protein expression. Type 1 IFN also did not activate p38, AKT, or ERK kinase. We conclude that type 1 IFN-induced activities in HeLa cells are mediated by STAT1/STAT2/IRF9, STAT1/STAT1, or STAT2/IRF9 complexes and do not require alternative pathways.

Type I interferons (IFNs) are key players in the antiviral immune response. Interferon alpha (IFN-α) belongs to this class of IFNs and comprises 12 subtypes that differ from each other in their binding affinities for a common receptor and, thus, in their signaling potencies. Recent data suggest that IFN-α6 and -α14 are the most potent IFN-α subtypes in restricting HIV replication when applied exogenously. However, in the context of antiviral therapy, IFNs are administered at high doses, which may compensate for differences in potency seen between IFN-α subtypes. In this study, we reexamined whether IFN-α subtypes induce different biological activities, with a focus on how IFN-α treatment dose affects cellular responses to HIV in primary CD4 ⁺ T cells, peripheral blood mononuclear cells (PBMCs), and macrophages. We found that the subtypes' antiviral activities were dose dependent, with > 0% inhibition of HIV replication at a high dose of all IFN-αs except the weak IFN- α/β receptor (IFNAR) binder, IFN-α1. The quality of the responses engendered by IFN-α1, -α2, -α6, and -α14 was highly comparable, with essentially the same set of genes induced by all four subtypes. Hierarchal cluster analysis revealed that the individual donors were stronger determinants for the IFN-stimulated-gene (ISG) responses than the specific IFN-α subtype used for stimulation. Notably, IFN-α2- derived mutants with substantially reduced IFNAR2 binding still inhibited HIV replication efficiently, whereas mutants with increased IFNAR1 binding potentiated antiviral activity. Overall, our results support the idea that IFN-α subtypes do not induce different biological responses, given that each subtype is exogenously applied at bioequivalent doses.

Protein-protein interactions mediate the vast majority of cellular processes. Though protein interactions obey basic chemical principles also within the cell, the in vivo physiological environment may not allow for equilibrium to be reached. Thus, in vitro measured thermodynamic affinity may not provide a complete picture of protein interactions in the biological context. Binding kinetics composed of the association and dissociation rate constants are relevant and important in the cell. Therefore, changes in protein-protein interaction kinetics have a significant impact on the in vivo activity of the proteins. The common protocol for the selection of tighter binders from a mutant library selects for protein complexes with slower dissociation rate constants. Here we describe a method to specifically select for variants with faster association rate constants by using pre-equilibrium selection, starting from a large random library. Toward this end, we refine the selection conditions of a TEM1-beta-lactamase library against its natural nanomolar affinity binder beta-lactamase inhibitor protein (BLIP). The optimal selection conditions depend on the ligand concentration and on the incubation time. In addition, we show that a second sort of the library helps to separate signal from noise, resulting in a higher percent of faster binders in the selected library. Fast associating protein variants are of particular interest for drug development and other biotechnological applications.

Chronic activation of the immune system in HIV infection is one of the strongest predictors of morbidity and mortality. As such, approaches that reduce immune activation have received considerable interest. Previously, we demonstrated that administration of a type I interferon receptor antagonist (IFN-1ant) during acute SIV infection of rhesus macaques results in increased virus replication and accelerated disease progression. Here, we administered a long half-life PASylated IFN-1ant to ART-treated and ART-naïve macaques during chronic SIV infection and measured expression of interferon stimulated genes (ISG) by RNA sequencing, plasma viremia, plasma cytokines, T cell activation and exhaustion as well as cell-associated virus in CD4 T cell subsets sorted from peripheral blood and lymph nodes. Our study shows that IFN-1ant administration in both ART-suppressed and ART-untreated chronically SIV-infected animals successfully results in reduction of IFN-I-mediated inflammation as defined by reduced expression of ISGs but had no effect on plasma levels of IL-1β, IL-1ra, IL-6 and IL-8. Unlike in acute SIV infection, we observed no significant increase in plasma viremia up to 25 weeks after IFN-1ant administration or up to 15 weeks after ART interruption. Likewise, cell-associated virus measured by SIV gag DNA copies was similar between IFN-1ant and placebo groups. In addition, evaluation of T cell activation and exhaustion by surface expression of CD38, HLA-DR, Ki67, LAG-3, PD-1 and TIGIT, as well as transcriptome analysis showed no effect of IFN-I blockade. Thus, our data show that blocking IFN-I signaling during chronic SIV infection suppresses IFN-I-related inflammatory pathways without increasing virus replication, and thus may constitute a safe therapeutic intervention in chronic HIV infection.

The HIV-1-encoded accessory protein Vpu exerts several immunomodulatory functions, including counteraction of the host restriction factor tetherin, downmodulation of CD4, and inhibition of NF-κB activity to facilitate HIV-1 infection. However, the relative contribution of individual Vpu functions to HIV-1 infection in vivo remained unclear. Here, we used a humanized mouse model and HIV-1 strains with selective mutations in vpu to demonstrate that the anti-tetherin activity of Vpu is a prerequisite for efficient viral spread during the early phase of infection. Mathematical modeling and gain-of-function mutations in SIVcpz, the simian precursor of pandemic HIV-1, corroborate this finding. Blockage of interferon signaling combined with transcriptome analyses revealed that basal tetherin levels are sufficient to control viral replication. These results establish tetherin as a key effector of the intrinsic immune defense against HIV-1, and they demonstrate that Vpu-mediated tetherin antagonism is critical for efficient viral spread during the initial phase of HIV-1 replication. The HIV-1-encoded accessory protein Vpu exerts several functions. Using a humanized mouse model and HIV-1 Vpu mutant viruses, Yamada et al. demonstrate that Vpu-mediated antagonism of the interferon-induced antiviral protein tetherin is critical for efficient viral spread during the initial phase of HIV-1 replication in vivo.

Proteins have evolved to balance efficient binding of desired partners with rejection of unwanted interactions. To investigate the evolution of protein-protein interactions, we selected a random library of pre-stabilized TEM1 beta-lactamase against wild-type TEM1 using yeast surface display. Three mutations were sufficient to achieve micromolar affinity binding between the two. The X-ray structure emphasized that the main contribution of the selected mutations was to modify the protein fold, specifically removing the N'-terminal helix, which consequently allowed protein coupling via a beta-sheet-mediated interaction resembling amyloid interaction mode. The only selected mutation located at the interaction interface (E58V) is reminiscent of the single mutation commonly causing sickle-cell anemia. Interestingly, the evolved mutations cannot be inserted into the wild-type protein due to reduced thermal stability of the resulting mutant protein. These results reveal a simple mechanism by which undesirable binding is purged by loss of thermal stability.

For over a century, enzymatic activity has been studied in vitro, assuming similar activity in the crowded cellular milieu. Here, we determined in real time the catalytic activity of TEM1-β-lactamase inside living cells and compared the values to those obtained in vitro. We found the apparent in vivo catalytic efficiency, k_cat/K_m, to be lower than in vitro, with significant cell-to-cell variability. Surprisingly, the results show that inside the cell the apparent catalytic efficiency decreases, and K_m increases with increasing enzyme concentration. To rationalize these findings, we measured enzyme and substrate diffusion rates in the cell and found the latter to be slower than expected. Simulations showed that for attenuated diffusion the substrate flux becomes rate-limiting, explaining why reaction rates in vivo can be independent on enzyme concentrations. The octanol/water partition of the substrate is 4.5, which is in the range of Food and Drug Administrationapproved drugs. This suggests substrate-limited reaction rates to be common. These findings indicate that in vitro data cannot be simply extrapolated to the crowded in vivo environment.

How structural dynamics affects cytokine signaling is under debate. Here, we investigated the dynamics of the type I interferon (IFN) receptor, IFNAR1, and its effect on signaling upon binding IFN and IFNAR2 using a combination of structure-based mechanistic studies, in situ binding, and gene induction assays. Our study reveals that IFNAR1 flexibility modulates ligand-binding affinity, which, in turn, regulates biological signaling. We identified the hinge sites and key interactions implicated in IFNAR1 inter-subdomain (SD1-SD4) movements. We showed that the predicted cooperative movements are essential to accommodate intermolecular interactions. Engineered disulfide bridges, computationally predicted to interfere with IFNAR1 dynamics, were experimentally confirmed. Notably, introducing disulfide bonds between subdomains SD2 and SD3 modulated IFN binding and activity in accordance with the relative attenuation of cooperative movements with varying distance from the hinge center, whereas locking the SD3 SD4 interface flexibility in favor of an extended conformer increased activity.

Type I interferons (IFN-1) are cytokines that affect the expression of thousands of genes, resulting in profound cellular changes. IFN-1 activates the cell by dimerizing its two-receptor chains, IFNAR1 and IFNAR2, which are expressed on all nucleated cells. Despite a similar mode of binding, the different IFN-1s activate a spectrum of activities. The causes for differential activation may stem from differences in IFN-1-binding affinity, duration of binding, number of surface receptors, induction of feedbacks, and cell type-specific variations. All together these will alter the signal that is transmitted from the extracellular domain inward. The intracellular domain binds, directly or indirectly, different effector proteins that transmit signals. The composition of effector molecules deviates between different cell types and tissues, inserting an additional level of complexity to the system. Moreover, IFN-1s do not act on their own, and clearly there is much cross-talk between the activated effector molecules by IFN-1 and other cytokines. The outcome generated by all of these factors (processing step) is an observed phenotype, which can be the transformation of the cell to an antiviral state, differentiation of the cell to a specific immune cell, senescence, apoptosis, and many more. IFN-1 activities can be divided into robust and tunable. Antiviral activity, which is stimulated by minute amounts of IFN-1 and is common to all cells, is termed robust. The other activities, which we term tunable, are cell type-specific and often require more stringent modes of activation. In this review, I summarize the current knowledge on the mode of activation and processing that is initiated by IFN-1, in perspective of the resulting phenotypes.

With computer-based data-fitting methods becoming a standard tool in biochemistry, progress curve analysis of enzyme kinetics is a feasible, yet seldom used tool. Here we present a versatile Matlab-based tool (PCAT) to analyze catalysis progress curves with three complementary model approaches. The first two models are based on the known closed-form solution for this problem: the first describes the required Lambert W function with an analytical approximation and the second provides a numerical solution of the Lambert W function. The third model is a direct simulation of the enzyme kinetics. Depending on the chosen model, the tools excel in speed, accuracy or initial value requirements. Using simulated and experimental data, we show the strengths and pitfalls of the different fitting models. Direct simulation proves to have the highest level of accuracy, but it also requires reasonable initial values to converge. Finally, we propose a standard procedure to obtain optimized enzyme kinetic parameters from single progress curves.

Type-I interferons (IFNs) play a key role in the immune defences against viral and bacterial infections, and in cancer immunosurveillance. We have established that clathrin-dependent endocytosis of the type-I interferon (IFN-α/β) receptor (IFNAR) is required for JAK/STAT signalling. Here we show that the internalized IFNAR1 and IFNAR2 subunits of the IFNAR complex are differentially sorted by the retromer at the early endosome. Binding of the retromer VPS35 subunit to IFNAR2 results in IFNAR2 recycling to the plasma membrane, whereas IFNAR1 is sorted to the lysosome for degradation. Depletion of VPS35 leads to abnormally prolonged residency and association of the IFNAR subunits at the early endosome, resulting in increased activation of STAT1-and IFN-dependent gene transcription. These experimental data establish the retromer complex as a key spatiotemporal regulator of IFNAR endosomal sorting and a new factor in type-I IFN-induced JAK/STAT signalling and gene transcription.

Protein-protein interactions occur via well-defined interfaces on the protein surface. Whereas the location of homologous interfaces is conserved, their composition varies, suggesting that multiple solutions may support high-affinity binding. In this study, we examined the plasticity of the interface of TEM1 β-lactamase with its protein inhibitor BLIP by low-stringency selection of a random TEM1 library using yeast surface display. Our results show that most interfacial residues could be mutated without a loss in binding affinity, protein stability, or enzymatic activity, suggesting plasticity in the interface composition supporting high-affinity binding. Interestingly, many of the selected mutations promoted faster association. Further selection for faster binders was achieved by drastically decreasing the library-ligand incubation time to 30 s. Preequilibrium selection as suggested here is a novel methodology for specifically selecting faster-associating protein complexes.

Type I interferons serve as the first line of defense against pathogen invasion. Binding of IFNs to its receptors, IFNAR1 and IFNAR2, is leading to activation of the IFN response. To determine whether structural perturbations observed during binding are propagated to the cytoplasmic domain, multiple mutations were introduced into the transmembrane helix and its surroundings. Insertion of one to five alanine residues near either the N or C terminus of the transmembrane domain (TMD) likely promotes a rotation of 100 degrees and a translation of 1.5 per added residue. Surprisingly, the added alanines had little effect on the binding affinity of IFN to the cell surface receptors, STAT phosphorylation, or gene induction. Similarly, substitution of the juxtamembrane residues of the TMD with alanines, or replacement of the TMD of IFNAR1 with that of IFNAR2, did not affect IFN binding or activity. Finally, only the addition of 10 serine residues (but not 2 or 4) between the extracellular domain of IFNAR1 and the TMD had some effect on signaling. Bioinformatic analysis shows a correlation between high sequence conservation of TMDs of cytokine receptors and the ability to transmit structural signals. Sequence conservation near the TMD of IFNAR1 is low, suggesting limited functional importance for this region. Our results suggest that IFN binding to the extracellular domains of IFNAR1 and IFNAR2 promotes proximity between the intracellular domains and that differential signaling is a function of duration of activation and affinity of binding rather than specific conformational changes transmitted from the outside to the inside of the cell.

Type I interferons are multi-potent cytokines that serve as first line of defense against viruses and other pathogens, posses immunomudolatory functions and elicit a growth inhibitory response. In recent years it has been shown that interferons are also detrimental, for example in lupus, AIDS, tuberculosis and cognitive decline, highlighted the need to develop interferon antagonists. We have previously developed the antagonist IFN-1ant, with much reduced binding to the IFNAR1 receptor and enhanced binding to IFNAR2. Here, we further tune the IFN-1ant by producing three additional antagonists based on IFN-1ant but with altered activity profiles. We show that in all three cases the antiproliferative activity of interferons is blocked and the induction of gene transcription of immunomudolatory and antiproliferative associated genes are substantially decreased. Conversely, each of the new antagonists elicits a different degree of antiviral response, STAT phosphorylation and related gene induction. Two of the new antagonists promote decreased activity in relation to the original IFN-1ant, while one of them promotes increased activity. As we do not know the exact causes of the detrimental effects of IFNs, the four antagonists that were produced and analyzed provide the opportunity to investigate the extent of antagonistic and agonistic activity optimal for a given condition.

We analysed gene expression microarray data from whole blood samples from 228 multiple sclerosis (MS) patients either untreated or treated with one of three alternative commonly used interferon beta (IFNβ) disease modifying drugs: Avonex® (×1 weekly), Betaseron® (every second day) or Rebif® (×3weekly). Patient injectionswere not timed to coordinate sample collections, thus providing a global transcriptomic profile for each population of patients studied. Three hundred and fifty one genes were significantly differentially expressed by at least one of the IFNβ drugs. Despite the different drug sources with distinct injection and dosage protocols, a striking similaritywas found in the identity and functional classes of the differentially expressed genes induced. Using the 25 most-upregulated genes, we defined a robust IFNβ gene expression signature that quantifies the IFN activation state per blood sample collected irrespective of the type of IFNβ therapy. This 25-gene signature also defined basal IFN activation states among untreated MS patients, which differed among individuals but remained relatively constant per patient with time. The maximum drug-induced IFN-activation state was similar for all three drugs despite a 1.7-2.0-fold diminished average effect for Avonex. This and a more erratic effect of Avonex per patient across longitudinal measurements is likely a result of its reduced injection frequency. In summary,we have defined a robust blood-derived type I IFN gene signature fromMS patients. This signature could potentially serve to generically quantify the systemic Type I IFN activation status for any other clinical manifestation, inclusive of other autoimmune diseases.

Type I interferons (IFNs) are best known for their role in innate immunity, but they are also involved in other functions including immunomodulation, restricting proliferation, cancer surveillance, and the regulation of the adaptive immune response. All these responses are mediated through the interaction with a single cell surface receptor, albeit at different ligand and receptor concentrations, ligand subtypes, and time of activation. Here we review the functional plasticity of IFN signaling from a quantitative perspective, showing how variations in different ingredients of the system lead to differential IFN responses and how cells tune the system to maximize efficiency while minimizing detrimental effects. We present a basic model wherein the integrated action of different feedback mechanisms can provide sufficient temporal control to differentially drive cellular decisions.

Agrobacterium is a pathogen that genetically transforms plants. The bacterial VirE2 protein envelopes the T-DNA of Agrobacterium and protects it from degradation. Within the transfected cells, VirE2 interacts with the plant VIP1 leading to nuclear transport of the T-DNA complex. Active VirE2 is an oligomer with a tendency to aggregate, hampering its studies at the molecular level. In addition, no structural or quantitative information is available regarding VIP1 or its interactions. The lack of information is mainly because both VIP1 and VirE2 are difficult to express and purify. Here, we present the development of efficient protocols that resulted in pure and stable Histagged VIP1 and VirE2. Circular dichroism spectroscopy and computational predictions indicated that VIP1 is mostly intrinsically disordered. This may explain the variety of protein-protein interactions it participates in. Size exclusion chromatography revealed that VirE2 exists in a two-state equilibrium between a monomer and an oligomeric form. Using the purified proteins, we performed peptide array screening and revealed the binding sites on both proteins. VirE2 binds the disordered regions of VIP1, while the site in VirE2 that binds VIP1 is different from the VirE2 DNA-binding site. Peptides derived from these sites may be used as lead compounds that block Agrobacterium infection of plants.

IFNβ is a common therapeutic option to treat multiple sclerosis. It is unique among the family of type I IFNs in that it binds to the interferon receptors with high affinity, conferring exceptional biological properties. We have previously reported the generation of an interferon superagonist (dubbed YNSα8) that is built on the backbone of a low affinity IFNα but modified to exhibit higher receptor affinity than even for IFNβ. Here, YNSα8 was fused with a 600-residue hydrophilic, unstructured N-terminal polypeptide chain comprising proline, alanine, and serine (PAS) to prolong its plasma half-life via "PASylation." PAS-YNSα8 exhibited a 10-fold increased half-life in both pharmacodynamic and pharmacokinetic assays in a transgenic mouse model harboring the human receptors, notably without any detectable loss in biological potency or bioavailability. This long-lived superagonist conferred significantly improved protection from MOG_35-55-induced experimental autoimmune encephalomyelitis compared with IFNβ, despite being injected with a 4-fold less frequency and at an overall 16-fold lower dosage. These data were corroborated by FACS measurements showing a decrease of CD11b⁺/CD45^hi myeloid lineage cells detectable in the CNS, as well as a decrease in IBA⁺ cells in spinal cord sections determined by immunohistochemistry for PAS-YNSα8-treated animals. Importantly, PAS-YNSα8 did not induce antibodies upon repeated administration, and its biological efficacy remained unchanged after 21 days of treatment. A striking correlation between increased levels of CD274 (PD-L1) transcripts from spleen-derived CD4⁺ cells and improved clinical response to autoimmune encephalomyelitis was observed, indicating that, at least in this mouse model of multiple sclerosis, CD274 may serve as a biomarker to predict the effectiveness of IFN therapy to treat this complex disease.

Type I interferons (IFNs), including various IFN-α isoforms and IFN-β, are a family of homologous, multifunctional cytokines. IFNs activate different cellular responses by binding to a common receptor that consists of two subunits, IFNAR1 and IFNAR2. In addition to stimulating antiviral responses, they also inhibit cell proliferation and modulate other immune responses. We characterized various IFNs, including a mutant IFN-α2 (IFN-1ant) that bound tightly to IFNAR2 but had markedly reduced binding to IFNAR1. Whereas IFN-1ant stimulated antiviral activity in a range of cell lines, it failed to elicit immunomodulatory and antiproliferative activities. The antiviral activities of the various IFNs tested depended on a set of IFN-sensitive genes (the "robust" genes) that were controlled by canonical IFN response elements and responded at low concentrations of IFNs. Conversely, these elements were not found in the promoters of genes required for the antiproliferative responses of IFNs (the " tunable" genes). The extent of expression of tunable genes was cell type-specific and correlated with the magnitude of the antiproliferative effects of the various IFNs. Although IFN-1ant induced the expression of robust genes similarly in five different cell lines, its antiviral activity was virus- and cell type-specific. Our findings suggest that IFN-1ant may be a therapeutic candidate for the treatment of specific viral infections without inducing the immunomodulatory and antiproliferative functions of wild-type IFN.

This Podcast features an interview with Gideon Schreiber, senior author of a Research Article that appears in the 27 May 2014 issue of Science Signaling, about an engineered variant of a type I interferon that elicits only antiviral responses. Interferons are major components of the innate immune response and serve as a link between the innate and adaptive immune responses. Type I interferons stimulate the expression of many genes that affect various cellular processes, including limiting the spread of viruses between cells, suppressing cell proliferation, and modulating the immune response. Levin et al. identified a variant of the type I interferon IFN-alpha 2 that elicits antiviral responses without suppressing cell proliferation or activating immunomodulatory genes.

In the last years, a growing interest has been gathering around the ability of Molecular Dynamics (MD) to provide insight into the paths of long-range structural communication in biomolecules. The knowledge of the mechanisms related to structural communication helps in the rationalization in atomistic details of the effects induced by mutations, ligand binding, and the intrinsic dynamics of proteins. We here present PyInteraph, a tool for the analysis of structural ensembles inspired by graph theory. PyInteraph is a software suite designed to analyze MD and structural ensembles with attention to binary interactions between residues, such as hydrogen bonds, salt bridges, and hydrophobic interactions. PyInteraph also allows the different classes of intra- and intermolecular interactions to be represented, combined or alone, in the form of interaction graphs, along with performing network analysis on the resulting interaction graphs. The program also integrates the network description with a knowledge-based force field to estimate the interaction energies between side chains in the protein. It can be used alone or together with the recently developed xPyder PyMOL plugin through an xPyder-compatible format. The software capabilities and associated protocols are here illustrated by biologically relevant cases of study. The program is available free of charge as Open Source software via the GPL v3 license at http://linux.btbs.unimib.it/ pyinteraph/.

We have generated transgenic mice that harbor humanized type I interferon receptors (IFNARs) enabling the study of type I human interferons (Hu-IFN-Is) in mice. These "HyBNAR" (Hybrid IFNAR) mice encode transgenic variants of IFNAR1 and IFNAR2 with the human extracellular domains being fused to transmembrane and cytoplasmic segments of mouse sequence. B16F1 mouse melanoma cells harboring the HyBNAR construct specifically bound Hu-IFN-Is and were rendered sensitive to Hu-IFN-I stimulated anti-proliferation, STAT1 activation and activation of a prototypical IFN-I response gene (MX2). HyBNAR mice were crossed with a transgenic strain expressing the luciferase reporter gene under the control of the IFN-responsive MX2 promoter (MX2-Luciferase). Both the HyBNAR and HyBNAR/MX2-Luciferase mice were responsive to all Hu-IFN-Is tested, inclusive of IFNα2A, IFNβ, and a human superagonist termed YNSα8. The mice displayed dose-dependent pharmacodynamic responses to Hu-IFN-I injection, as assessed by measuring the expression of IFN-responsive genes. Our studies also demonstrated a weak activation of endogenous mouse interferon response, especially after high dose administration of Hu-IFNs. In sharp contrast to data published for humans, our pharmacodynamic readouts demonstrate a very short-lived IFN-I response in mice, which is not enhanced by sub-cutaneous (SC) injections in comparison to other administration routes. With algometric differences between humans and mice taken into account, the HyBNAR mice provides a convenient non-primate pre-clinical model to advance the study of human IFN-Is.

Inflammation in HIV infection is predictive of non-AIDS morbidity and death, higher set point plasma virus load² and virus acquisition; thus, therapeutic agents are in development to reduce its causes and consequences. However, inflammation may simultaneously confer both detrimental and beneficial effects. This dichotomy is particularly applicable to type I interferons (IFN-I) which, while contributing to innate control of infection^4-10, also provide target cells for the virus during acute infection, impair CD4 T-cell recovery, and are associated with disease progression^6,7,11-19. Here we manipulated IFN-I signalling in rhesus macaques (Macaca mulatta) during simian immunodeficiency virus (SIV) transmission and acute infection with two complementary in vivo interventions. We show that blockade of the IFN-I receptor caused reduced antiviral gene expression, increased SIV reservoir size and accelerated CD4 T-cell depletion with progression to AIDS despite decreased T-cell activation. In contrast, IFN-α2a administration initially upregulated expression of antiviral genes and prevented systemic infection. However, continued IFN-α2a treatment induced IFN-I desensitization and decreased antiviral gene expression, enabling infection with increased SIV reservoir size and accelerated CD4 T-cell loss. Thus, the timing of IFN-induced innate responses in acute SIV infection profoundly affects overall disease course and outweighs the detrimental consequences of increased immune activation. Yet, the clinical consequences of manipulation of IFN signalling are difficult to predict in vivo and therapeutic interventions in human studies should be approached with caution.

A long-term aim of computational design is to generate specific protein-protein interactions at desired affinity, specificity, and kinetics. The past three years have seen the first reports on atomically accurate de novo interactions. These were based on advances in design algorithms and the ability to harness high-throughput experimental characterization of design variants to optimize binding. Current state-of-the-art in computational design lacks precision, and therefore requires intensive experimental optimization to achieve parity with natural binders. Recent successes (and failures) point the way to future progress in design methodology that would enable routine and robust design of binders and inhibitors, while also shedding light on the essential features of biomolecular recognition.

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

Traditionally, biochemical studies are performed in dilute homogenous solutions, which are very different from the dense mixture of molecules found in cells. Thus, the physiological relevance of these studies is in question. This recognition motivated scientists to formulate the effect of crowded solutions in general, and excluded volume in particular, on biochemical processes. Using polymers or proteins as crowders, it was shown that while crowding tends to significantly enhance the formation of complexes containing many subunits, dimerizations are only mildly affected. Computer simulations, together with experimental evidence, indicate soft interactions and diffusion as critical factors that operate in a concerted manner with excluded volume to modulate protein binding. Yet, these approaches do not truly mimic the cellular environment. In vivo studies may overcome this shortfall. The few studies conducted thus far suggest that in cells, binding and folding occur at rates close to those determined in dilute solutions. Obtaining quantitative biochemical information on reactions inside living cells is currently a main challenge of the field, as the complexity of the intracellular milieu was what motivated crowding research to begin with.

Microscale thermophoresis (MST) allows for quantitative analysis of protein interactions in free solution and with low sample consumption. The technique is based on thermophoresis, the directed motion of molecules in temperature gradients. Thermophoresis is highly sensitive to all types of binding-induced changes of molecular properties, be it in size, charge, hydration shell or conformation. In an all-optical approach, an infrared laser is used for local heating, and molecule mobility in the temperature gradient is analyzed via fluorescence. In standard MST one binding partner is fluorescently labeled. However, MST can also be performed label-free by exploiting intrinsic protein UV-fluorescence. Despite the high molecular weight ratio, the interaction of small molecules and peptides with proteins is readily accessible by MST. Furthermore, MST assays are highly adaptable to fit to the diverse requirements of different biomolecules, such as membrane proteins to be stabilized in solution. The type of buffer and additives can be chosen freely. Measuring is even possible in complex bioliquids like cell lysate allowing close to in vivo conditions without sample purification. Binding modes that are quantifiable via MST include dimerization, cooperativity and competition. Thus, its flexibility in assay design qualifies MST for analysis of biomolecular interactions in complex experimental settings, which we herein demonstrate by addressing typically challenging types of binding events from various fields of life science.

Interferons induce a pleiotropy of responses through binding the same cell surface receptor. Here we investigated the molecular mechanism driving interferon-induced apoptosis. Using a nonbiased small interferingRNA(siRNA) screen, we show that silencing genes whose products are directly engaged in the initiation of interferon signaling completely abrogate the interferon antiproliferative response. Apoptosis-related genes such as the caspase-8, cFLIP, and DR5 genes specifically interfere with interferon-induced apoptosis, which we found to be independent of the activity of death ligands. The one gene for which silencing resulted in the strongest proapoptotic effect upon interferon signaling is the cFLIP gene, where silencing shortened the time of initiation of apoptosis from days to hours and increased dramatically the population of apoptotic cells. Thus, cFLIP serves as a regulator for interferon-induced apoptosis. A shift over time in the balance between cFLIP and caspase-8 results in downstream caspase activation and apoptosis. While gamma interferon (IFN- γ) also causes caspase-8 upregulation, we suggest that it follows a different path to apoptosis.

Protein-protein interactions (PPIs) are essential for cellular viability and activity. Here, we present a rapid, semi-quantitative method (termed FRETex) to analyze PPIs, taking advantage of the strong and specific FRET signal between fused CyPET donor and YPET acceptor molecules. To demonstrate the robustness of this approach, we analyzed the interactions between three protein pairs and their muteins: TEM1-β-lactamase binding its inhibitor BLIP, barnase binding barstar and ornithine decarboxylase binding its inhibitor antizyme. The CyPET/YPET fused proteins were produced in small quantities, and the measurements were conducted directly in the proteins crude Escherichia coli lysates without any purification step. Protein concentrations were determined from the fluorescence intensities of the lysates. While binding titration curves were produced, the resulting affinities were not always precise. Therefore, we also conducted time-resolved chase experiments using non-labeled binding partners as chasers. The acquired dissociation rate constants were in a good agreement with those measured by surface plasmon resonance. Due to the simplicity of FRETex, and the ability to obtain semi-quantitative binding data, FRETex is a suitable method for tasks such as mutant scans, protein-engineering, scanning for inhibitors and more.

Type I interferons (IFNs) form a network of homologous cytokines that bind to a shared, heterodimeric cell surface receptor and engage signaling pathways that activate innate and adaptive immune responses. The ability of IFNs to mediate differential responses through the same cell surface receptor has been subject of a controversial debate and has important medical implications. During the past decade, a comprehensive insight into the structure, energetics, and dynamics of IFN recognition by its two-receptor subunits, as well as detailed correlations with their functional properties on the level of signal activation, gene expression, and biological responses were obtained. All type I IFNs bind the two-receptor subunits at the same sites and form structurally very similar ternary complexes. Differential IFN activities were found to be determined by different lifetimes and ligand affinities toward the receptor subunits, which dictate assembly and dynamics of the signaling complex in the plasma membrane. We present a simple model, which explains differential IFN activities based on rapid endocytosis of signaling complexes and negative feedback mechanisms interfering with ternary complex assembly. More insight into signaling pathways as well as endosomal signaling and trafficking will be required for a comprehensive understanding, which will eventually lead to therapeutic applications of IFNs with increased efficacy.

The crowded environment of cells poses a challenge for rapid protein-protein association. Yet, it has been established that the rates of association are similar in crowded and in dilute solutions. Here we probe the pathway leading to fast association between TEM1 β-lactamase and its inhibitor protein BLIP in crowded solutions. We show that the affinity of the encounter complex, the rate of final complex formation, and the structure of the transition state are similar in crowded solutions and in buffer. The experimental results were reproduced by calculations based on the transient-complex theory for protein association. Both experiments and calculations suggest that while crowding agents decrease the diffusion constant of the associating proteins, they also induce an effective excluded-volume attraction between them. The combination of the two opposing effects thus results in nearly identical overall association rates in diluted and crowded solutions.

Historically, rate constants were determined in vitro and it was unknown whether they were valid for in vivo biological processes. Here, we bridge this gap by measuring binding dynamics between a pair of proteins in living HeLa cells. Binding of a β-lactamase to its protein inhibitor was initiated by microinjection and monitored by Förster resonance energy transfer. Association rate constants for the wild-type and an electrostatically optimized mutant were only 25% and 50% lower than in vitro values, whereas no change in the rate constant was observed for a slower binding mutant. These changes are much smaller than might be anticipated considering the high macromolecular crowding within the cell. Single-cell analyses of association rate constants and fluorescence recovery after photobleaching reveals a naturally occurring variation in cell density, which is translated to an up to a twofold effect on binding rate constants. The data show that for this model protein interaction the intracellular environment had only a small effect on the association kinetics, justifying the extrapolation of in vitro data to processes in the cell.

The CAPRI (Critical Assessment of Predicted Interactions) and CASP (Critical Assessment of protein Structure Prediction) experiments have demonstrated the power of community-wide tests of methodology in assessing the current state of the art and spurring progress in the very challenging areas of protein docking and structure prediction. We sought to bring the power of community-wide experiments to bear on a very challenging protein design problem that provides a complementary but equally fundamental test of current understanding of protein-binding thermodynamics. We have generated a number of designed protein-protein interfaces with very favorable computed binding energies but which do not appear to be formed in experiments, suggesting that there may be important physical chemistry missing in the energy calculations. A total of 28 research groups took up the challenge of determining what is missing: we provided structures of 87 designed complexes and 120 naturally occurring complexes and asked participants to identify energetic contributions and/or structural features that distinguish between the two sets. The community found that electrostatics and solvation terms partially distinguish the designs from the natural complexes, largely due to the nonpolar character of the designed interactions. Beyond this polarity difference, the community found that the designed binding surfaces were, on average, structurally less embedded in the designed monomers, suggesting that backbone conformational rigidity at the designed surface is important for realization of the designed function. These results can be used to improve computational design strategies, but there is still much to be learned; for example, one designed complex, which does form in experiments, was classified by all metrics as a nonbinder.

Vertebrates have multiple genes encoding Type I interferons (IFN), for reasons that are not fully understood. The Type I IFN appear to bind to the same heterodimeric receptor and the subtypes have been shown to have different potencies in various experimental systems. To put this concept on a quantitative basis, we have determined the binding affinities and rate constants of 12 human Alpha-IFN subtypes to isolated interferon receptor chains 1 and 2. Alpha-IFNs bind IFNAR1 and IFNAR2 at affinities of 0.5-5μM and 0.4-5nM respectively (except for IFN-alpha1 - 220nM). Additionally we have examined the biological activity of these molecules in several antiviral and antiproliferative models. Particularly for antiproliferative potency, the binding affinity and activity correlate. However, the EC ₅₀ values differ significantly (1.5nM versus 0.1nM for IFN-alpha2 in WISH versus OVCAR cells). For antiviral potency, there are several instances where the relationship appears to be more complicated than simple binding. These results will serve as a point of reference for further understanding of this multiple ligand/receptor system.

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

Type I Interferons (IFNs) are important cytokines for innate immunity against viruses and cancer. Sixteen human type I IFN variants signal through the same cell-surface receptors, IFNAR1 and IFNAR2, yet they can evoke markedly different physiological effects. The crystal structures of two human type I IFN ternary signaling complexes containing IFNα2 and IFNω reveal recognition modes and heterotrimeric architectures that are unique among the cytokine receptor superfamily but conserved between different type I IFNs. Receptor-ligand cross-reactivity is enabled by conserved receptor-ligand "anchor points" interspersed among ligand-specific interactions that "tune" the relative IFN-binding affinities, in an apparent extracellular "ligand proofreading" mechanism that modulates biological activity. Functional differences between IFNs are linked to their respective receptor recognition chemistries, in concert with a ligand-induced conformational change in IFNAR1, that collectively control signal initiation and complex stability, ultimately regulating differential STAT phosphorylation profiles, receptor internalization rates, and downstream gene expression patterns. PaperFlick:

Type I interferons trigger diverse biological effects by binding a common receptor, composed of IFNAR1 and IFNAR2. Intriguingly, while the activation of an antiviral state is common to all cells, antiproliferative activity and apoptosis affect only part of the population, even when cells are stimulated with saturating interferon concentrations. Manipulating receptor expression by different small interfering RNA (siRNA) concentrations reduced the fraction of responsive cells independent of the interferon used, including a newly generated, extremely tight-binding variant. Reduced receptor numbers increased 50% effective concentrations (EC ₅₀s) for alpha interferon 2 (IFN-α2) but not for the tight-binding variant. A correlation between receptor numbers, STAT activation, and gene induction is observed. Our data suggest that for a given cell, the response is binary (+/-) and dependent on the stochastic expression levels of the receptors on an individual cell. A low number of receptors suffices for antiviral response and is thus a robust feature common to all cells. Conversely, a high number of receptors is required for antiproliferative activity, which allows for fine-tuning on a single-cell level.

The "4D Biology Workshop for Health and Disease" , held on 16-17th of March 2010 in Brussels, aimed at finding the best organising principles for large-scale proteomics, interactomics and structural genomics/biology initiatives, and setting the vision for future high-throughput research and large-scale data gathering in biological and medical science.Major conclusions of the workshop include the following. (i) Development of new technologies and approaches to data analysis is crucial. Biophysical methods should be developed that span a broad range of time/spatial resolution and characterise structures and kinetics of interactions. Mathematics, physics, computational and engineering tools need to be used more in biology and new tools need to be developed. (ii) Database efforts need to focus on improved definitions of ontologies and standards so that system-scale data and associated metadata can be understood and shared efficiently. (iii) Research infrastructures should play a key role in fostering multidisciplinary research, maximising knowledge exchange between disciplines and facilitating access to diverse technologies. (iv) Understanding disease on a molecular level is crucial. System approaches may represent a new paradigm in the search for biomarkers and new targets in human disease. (v) Appropriate education and training should be provided to help efficient exchange of knowledge between theoreticians, experimental biologists and clinicians. These conclusions provide a strong basis for creating major possibilities in advancing research and clinical applications towards personalised medicine.

Interactions between macromolecules in general, and between proteins in particular, are essential for any life process. Examples include transfer of information, inhibition or activation of function, molecular recognition as in the immune system, assembly of macromolecular structures and molecular machines, and more. Proteins interact with affinities ranging from millimolar to femtomolar and, because affinity determines the concentration required to obtain 50% binding, the amount of different complexes formed is very much related to local concentrations. Although the concentration of a specific binding partner is usually quite low in the cell (nanomolar to micromolar), the total concentration of other macromolecules is very high, allowing weak and non-specific interactions to play important roles. In this review we address the question of binding specificity, that is, how do some proteins maintain monogamous relations while others are clearly polygamous. We examine recent work that addresses the molecular and structural basis for specificity versus promiscuity. We show through examples how multiple solutions exist to achieve binding via similar interfaces and how protein specificity can be tuned using both positive and negative selection (specificity by demand). Binding of a protein to numerous partners can be promoted through variation in which residues are used for binding, conformational plasticity and/or post-translational modification. Natively unstructured regions represent the extreme case in which structure is obtained only upon binding. Many natively unstructured proteins serve as hubs in protein-protein interaction networks and such promiscuity can be of functional importance in biology.

Interactions of peptides and proteins with inorganic surfaces are important to both natural and artificial systems; however, a detailed understanding of such interactions is lacking. In this study, we applied new approaches to quantitatively measure the binding of amino acids and proteins to gold surfaces. Real-time surface plasmon resonance (SPR) measurements showed that TEM1-β-lactamase inhibitor protein (BLIP) interacts only weakly with Au nanoparticles (NPs). However, fusion of three histidine residues to BLIP (3H-BLIP) resulted in a significant increase in the binding to the Au NPs, which further increased when the histidine tail was extended to six histidines (6H-BLIP). Further increasing the number of His residues had no effect on the binding. A parallel study using continuous (111)-textured Au surfaces and single-crystalline, (111)-oriented, Au islands by ellipsometry, FTIR, and localized surface plasmon resonance (LSPR) spectroscopy further confirmed the results, validating the broad applicability of Au NPs as model surfaces. Evaluating the binding of all other natural amino acid homotripeptides fused to BLIP (except Cys and Pro) showed that aromatic and positively-charged residues bind preferentially to Au with respect to small aliphatic and negatively charged residues, and that the rate of association is related to the potency of binding. The binding of all fusions was irreversible. These findings were substantiated by SPR measurements of synthesized, free, soluble tripeptides using Au-NP-modified SPR chips. Here, however, the binding was reversible allowing for determination of binding affinities that correlate with the binding potencies of the related BLIP fusions. Competition assays performed between 3H-BLIP and the histidine tripeptide (3 His) suggest that Au binding residues promote the adsorption of proteins on the surface, and by this facilitate the irreversible interaction of the polypeptide chain with Au. The binding of amino acids to Au was simulated by using a continuum solvent model, showing agreement with the experimental values. These results, together with the observed binding potencies and kinetics of the BLIP fusions and free peptides, suggest a binding mechanism that is markedly different from biological protein-protein interactions.

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

Cytokines are essential proteins that exert potent control over entire cell populations to fight infections and other pathologies, but can by themselves cause disease. Therefore, cytokine-related drugs act either by stimulating or blocking their activities. Our knowledge of the structures of cytokine-receptor complexes, the biophysical basis of their binding, and their mode of biological activation has substantially increased in recent years. This knowledge has been translated into new drugs and drug candidates. This review summarizes our current understanding of the receptor-mediated activity of cytokines, their relation to health and disease, and the agents in use to activate and block their actions.

Background: Accurate evaluation and modelling of residue-residue interactions within and between proteins is a key aspect of computational structure prediction including homology modelling, protein-protein docking, refinement of low-resolution structures, and computational protein design.Results: Here we introduce a method for accurate protein structure modelling and evaluation based on a novel 4-distance description of residue-residue interaction geometry. Statistical 4-distance preferences were extracted from high-resolution protein structures and were used as a basis for a knowledge-based potential, called Hunter. We demonstrate that 4-distance description of side chain interactions can be used reliably to discriminate the native structure from a set of decoys. Hunter ranked the native structure as the top one in 217 out of 220 high-resolution decoy sets, in 25 out of 28 "Decoys 'R' Us" decoy sets and in 24 out of 27 high-resolution CASP7/8 decoy sets. The same concept was applied to side chain modelling in protein structures. On a set of very high-resolution protein structures the average RMSD was 1.47 Å for all residues and 0.73 Å for buried residues, which is in the range of attainable accuracy for a model. Finally, we show that Hunter performs as good or better than other top methods in homology modelling based on results from the CASP7 experiment. The supporting web site http://bioinfo.weizmann.ac.il/hunter/ was developed to enable the use of Hunter and for visualization and interactive exploration of 4-distance distributions.Conclusions: Our results suggest that Hunter can be used as a tool for evaluation and for accurate modelling of residue-residue interactions in protein structures. The same methodology is applicable to other areas involving high-resolution modelling of biomolecules.

Protein-surface interactions are fundamental in natural processes, and have great potential for applications ranging from nanotechnology to medicine. A recent workshop highlighted the current achievements and the main challenges in the field.

This chapter discusses the energy landscape separating the unbound from the bound state of proteinprotein interactions. Specific proteinprotein interactions provide a major part of the basic organization of living cells. The structure of a protein complex embeds the information of the relative mutual organization of two proteins in a frozen state while the affinity allows analysis of the equilibrium thermodynamics of the interaction. Analysis of the kinetics of association and dissociation in solution or within membranes allows for characterizing the landscape in more detail. Single-point mutation analysis has often resulted in perplexing results, which are difficult to simulate using computation. This can be attributed to the cooperative effect of the binding of a group of residues, as is the case in proteinprotein interfaces, and to potential structural changes resulting from mutations. This cooperativity is explained in terms of the modular architecture of binding sites. Compared to interaction in solution, anchoring of proteins into the membrane severely changes the energetic landscapes of proteinprotein interaction. While the long-range interactions formed during association are a global feature of the protein, the architecture of the interface is modular, with the individual modules being cooperative within themselves but additive between modules.

Summary: Protein structures can be viewed as networks of contacts (edges) between amino-acid residues (nodes). Here we dissect proteins into sub-graphs consisting of six nodes and their corresponding edges, with an edge being either a backbone hydrogen bond (H-bond) or a covalent interaction. Six thousand three hundred and twenty-two such sub-graphs were found in a large non-redundant dataset of high-resolution structures, from which 35 occur much more frequently than in a random model. Many of these significant sub-graphs (also called network motifs) correspond to sub-structures of α helices and β-sheets, as expected. However, others correspond to more exotic sub-structures such as 3₁₀ helix, Schellman motif and motifs that were not defined previously. This topological characterization of patterns is very useful for producing a detailed differences map to compare protein structures. Here we analyzed in details the differences between NMR, molecular dynamics (MD) simulations and X-ray structures for Lysozyme, SH3 and the lambda repressor. In these cases, the same structures solved by NMR and simulated by MD showed small but consistent differences in their motif composition from the crystal structures, despite a very small root mean square deviation (RMSD) between them. This may be due to differences in the pair-wise energy functions used and the dynamic nature of these proteins.

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

Methods for protein modeling and design advanced rapidly in recent years. At the heart of these computational methods is an energy function that calculates the free energy of the system. Many of these functions were also developed to estimate the consequence of mutation on protein stability or binding affinity. In the current study, we chose six different methods that were previously reported as being able to predict the change in protein stability (AAG) upon mutation: CC/PBSA, EGAD, FoldX, I-Mutant2.0, Rosetta and Hunter. We evaluated their performance on a large set of 2156 single mutations, avoiding for each program the mutations used for training. The correlation coefficients between experimental and predicted AAG values were in the range of 0.59 for the best and 0.26 for the worst performing method. All the tested computational methods showed a correct trend in their predictions, but failed in providing the precise values. This is not due to lack in precision of the experimental data, which showed a correlation coefficient of 0.86 between different measurements. Combining the methods did not significantly improve prediction accuracy compared to a single method. These results suggest that there is still room for improvement, which is crucial if we want forcefields to perform better in their various tasks.

Multiple type I interferons (IFN-α/β) elicit Jak/Stat activation, rapid gene induction, and pleiotropic effects, such as differentiation, antiviral protection, and blocks in proliferation, which are dependent on the IFN subtype and the cellular context. To date, ligand- and receptor-specific molecular determinants underlying IFN-α/β differential activities or potencies have been well characterized. To analyze cellular determinants that impact subtype-specific potency, human fibrosarcoma U5A-derived clones, exhibiting a gradient of IFN sensitivity by virtue of increasing receptor levels, were monitored for Jak/Stat signaling, gene induction, cell cycle lengthening, and apoptosis. In cells with scarce receptors, IFN-β was more potent than IFN-α2 in antiproliferative activity, while the two subtypes were equipotent in all other readouts. Conversely, in cells with abundant receptors, IFN-α2 matched or even surpassed IFN-β in all readouts tested. Our results suggest that the differential activities of the IFN subtypes are dictated not only by the intrinsic ligand/receptor binding kinetics but also by the density of cell surface receptor components.

The three-dimensional structures of proteins are stabilized by the interactions between amino acid residues. Here we report a method where four distances are calculated between any two side chains to provide an exact spatial definition of their bonds. The data were binned into a four-dimensional grid and compared to a random model, from which the preference for specific four-distances was calculated. A clear relation between the quality of the experimental data and the tightness of the distance distribution was observed, with crystal structure data providing far tighter distance distributions than NMR data. Since the four-distance data have higher information content than classical bond descriptions, we were able to identify many unique inter-residue features not found previously in proteins. For example, we found that the side chains of Arg, Glu, Val and Leu are not symmetrical in respect to the interactions of their head groups. The described method may be developed into a function, which computationally models accurately protein structures.

Antizyme (Az) is a highly conserved key regulatory protein bearing a major role in regulating polyamine levels in the cell. It has the ability to bind and inhibit ornithine decarboxylase (ODC), targeting it for degradation. Az inhibitor (AzI) impairs the activity of Az. In this study, we mapped the binding sites of ODC and AzI on Az using Ala scan mutagenesis and generated models of the two complexes by constrained computational docking. In order to scan a large number of mutants in a short time, we developed a workflow combining high-throughput mutagenesis, small-scale parallel partial purification of His-tagged proteins and their immobilization on a tris-nitrilotriacetic-acid-coated surface plasmon resonance chip. This combination of techniques resulted in a significant reduction in time for production and measurement of large numbers of mutant proteins. The data-driven docking results suggest that both proteins occupy the same binding site on Az, with Az binding within a large groove in AzI and ODC. However, single-mutant data provide information concerning the location of the binding sites only, not on their relative orientations. Therefore, we generated a large number of double-mutant cycles between residues on Az and ODC and used the resulting interaction energies to restrict docking. The model of the complex is well defined and accounts for the mutant data generated here, and previously determined biochemical data for this system. Insights on the structure and function of the complexes, as well as general aspects of the method, are discussed.

A review of some of the advancements made in dealing with kinetic pathway of protein complex formation are presented. The review also focused on the nature of the precomplex formed through diffusion, the transition state, other intermediates, and the association pathway for the protein. It was revealed that the association rate played a key role in such processes, along with other processes that needed speed. It was also revealed that the reorganization of the Actin cytoskeleton demonstrated the significance of rapid protein association. Actin cytoskeleton reorganization was attained through Actin polymerization, which was nucleated by the Arp2/3 complex. It was observed that the speed of upstream signaling had a critical impact on the rate of polymer formation, as Actin polymerization was initiated with a nucleation process.

Studies of protein-protein interactions, carried out in polymer solutions, are designed to mimic the crowded environment inside living cells. It was shown that crowding enhances oligomerization and polymerization of macromolecules. Conversely, we have shown that crowding has only a small effect on the rate of association of protein complexes. Here, we investigated the equilibrium effects of crowding on protein heterodimerization of TEM1-β-lactamase with β-lactamase inhibitor protein (BLIP) and barnase with barstar. We also contrasted these with the effect of crowding on the weak binding pair CyPet-YPet. We measured the association and dissociation rates as well as the affinities and thermodynamic parameters of these interactions in polyethylene glycol and dextran solutions. For TEM1-BLIP and for barnase-barstar, only a minor reduction in association rate constants compared to that expected based on solution viscosity was found. Dissociation rate constants showed similar levels of reduction. Overall, this resulted in a binding affinity that is quite similar to that in aqueous solutions. On the other hand, for the CyPet-YPet pair, aggregation, and not enhanced dimerization, was detected in polyethylene glycol solutions. The results suggest that typical crowding agents have only a small effect on specific protein-protein dimerization reactions. Although crowding in the cell results from proteins and other macromolecules, one may still speculate that binding in vivo is not very different from that measured in dilute solutions.

The association reaction between pairs of proteins proceeds through an encounter complex that develops into the final complex. Here, we combined Brownian dynamics simulations with experimental studies to analyze the structures of the encounter complexes along the association reaction between TEM1-β-lactamase and its inhibitor, β-lactamase-inhibitor protein. The encounter complex can be considered as an ensemble of short-lived low free-energy states that are stabilized primarily by electrostatic forces and desolvation. For the wild-type, the simulation showed two main encounter regions located outside the physical binding site. One of these regions was located near the experimentally determined transition state. To validate whether these encounters are fruitful or futile, we examined three groups of mutations that altered the encounter. The first group consisted of mutations that increased the experimental rate of association through electrostatic optimization. This resulted in an increase in the size of the encounter region located near the experimentally determined transition state, as well as a decrease in the energy of this region and an increase in the number of successful trajectories (i.e., encounters that develop into complex). A second group of mutations was specifically designed to either increase or decrease the size and energy of the second encounter complex, but either way it did not affect k_on. A third group of mutations consisted of residues that increased k_on without significantly affecting the encounter complexes. These results indicate that the size and energy of the encounter regions are only two of several parameters that lead to fruitful association, and that electrostatic optimization is a major driving force in fast association.

The GTP-binding protein Ras plays a central role in the regulation of various cellular processes, acting as a molecular switch that triggers signaling cascades. Only Ras bound to GTP is able to interact strongly with effector proteins like Raf kinase, phosphatidylinositol 3-kinase, and RalGDS, whereas in the GDP-bound state, the stability of the complex is strongly decreased, and signaling is interrupted. To determine whether this process is only controlled by the stability of the complex, we used computer-aided protein design to improve the interaction between Ras and effector. We challenged the Raf·Raf complex in this study because Raf among all effectors shows the highest Ras affinity and the fastest association kinetics. The proposed mutations were characterized as to their changes in dynamics and binding strength. We demonstrate that Ras-Raf interaction can only be improved at the cost of a loss in specificity of Ras·GTP versus Ras·GDP. As shown by NMR spectroscopy, the Raf mutation A85K leads to a shift of Ras switch I in the GTP-bound as well as in the GDP-bound state, thereby increasing the complex stability. In a luciferase-based reporter gene assay, Raf A85K is associated with higher signaling activity, which appears to be a mere matter of Ras-Raf affinity.

A new method is presented for the redesign of protein-protein interfaces, resulting in specificity of the designed pair while maintaining high affinity. The design is based on modular interface architecture and was carried out on the interaction between TEM1 β-lactamase and its inhibitor protein, β-lactamase inhibitor protein. The interface between these two proteins is composed of several mostly independent modules. We previously showed that it is possible to delete a complete module without affecting the overall structure of the interface. Here, we replace a complete module with structure fragments taken from nonrelated proteins. Nature-optimized fragments were chosen from 10⁷ starting templates found in the Protein Data Bank. A procedure was then developed to identify sets of interacting template residues with a backbone arrangement mimicking the original module. This generated a final list of 361 putative replacement modules that were ranked using a novel scoring function based on grouped atom-atom contact surface areas. The top-ranked designed complex exhibited an affinity of at least the wild-type level and a mode of binding that was remarkably specific despite the absence of negative design in the procedure. In retrospect, the combined application of three factors led to the success of the design approach: utilizing the modular construction of the interface, capitalizing on native rather than artificial templates, and ranking with an accurate atom-atom contact surface scoring function.

Surface plasmon resonance (SPR) technology is a central, widely used tool for kinetic studies of interactions between unlabeled biomolecules in real time. The ProteOn XPR36 protein interaction array system possesses all the qualities of high-level SPR biosensing technology combined with high-throughput and multiplexing capabilities. This system is based on the built-in orientation controlled multi-channel module of the ProteOn XPR36 system, which allows parallel measurement of multiple binding interactions between as many as six protein pairs. This kind of multiplexing is done efficiently on the ProteOn XPR36 system using the innovative technique of One-shot Kinetics (TM), which allows simultaneous monitoring of multiple protein pair interactions. Here we demonstrate the capabilities of the ProteOn XPR36 for the measurements of a wide range of biomolecules interactions including proteins, DNA oligos, small molecules and antibodies and discuss its optimal mode of use.

Type I interferons (IFNs) signal for their diverse biological effects by binding a common receptor on target cells, composed of the two transmembrane IFNAR1 and IFNAR2 proteins. We have previously differentially enhanced the antiproliferative activity of IFN by increasing the weak binding affinity of IFN to IFNAR1. In this study, we further explored the affinity interdependencies between the two receptor subunits and the role of IFNAR1 in differential IFN activity. For this purpose, we generated a panel of mutations targeting the IFNAR2 binding site on the background of the IFNα2 YNS mutant, which increases the affinity to IFNAR1 by 60-fold, resulting in IFNAR2-to-IFNAR1 binding affinity ratios ranging from 1000:1 to 1:1000. Both the antiproliferative and antiviral potencies of the interferon mutants clearly correlated to the in situ binding IC₅₀ values, independently of the relative contributions of the individual receptors, thus relating to the integral lifetime of the complex. However, the antiproliferative potency correlated throughout the entire range of affinities, as well as with prolonged IFNAR1 receptor down-regulation, whereas the antiviral potency reached a maximum at binding affinities equivalent to that of wild-type IFNα2. Our data suggest that (i) the specific activity of interferon is related to the ternary complex binding affinity and not to affinity toward individual receptor components and (ii) although the antiviral pathway is strongly dependent on pSTAT1 activity, the cytostatic effect requires additional mechanisms that may involve IFNAR1 down-regulation. This differential interferon response is ultimately mediated through distinct gene expression profiling.

Type I interferons (IFNs) are multifunctional cytokines that activate cellular responses by binding a common receptor consisting of two subunits, IFNAR-1 and IFNAR-2. Although the binding of IFNs to IFNAR-2 is well characterized, the binding to the lower affinity IFNAR-1 remains less well understood. Previous reports identified a region of human IFN-α2 on the B and C helices ("site IA": N65, L80, Y85, Y89) that plays a key role in binding IFNAR-1 and contributes strongly to differential activation by various type I IFNs. The current studies demonstrate that residues on the D helix are also involved in IFNAR-1 binding. In particular, residue 120 (Arg in IFN-α2; Lys in IFN-α2/α1) appears to be a "hot-spot" residue: substitution by alanine significantly decreased biological activity, and the charge-reversal mutation of residue 120 to Glu caused drastic loss of antiviral and antiproliferative activity for both IFN-α2 and IFN-α2/α1. Mutations in residues of helix D maintained their affinity for IFNAR-2 but had decreased affinity for IFNAR-1. Single-site or multiple-site mutants in the IFNAR-1 binding site that had little or no detectable in vitro biological activity were capable of blocking in vitro antiviral and antiproliferative activity of native IFN-α2; i.e., they are type I IFN antagonists. These prototype IFN antagonists can be developed further for possible therapeutic use in systemic lupus erythematosus, and analogous molecules can be designed for use in animal models.

Proteins fold into a well-defined structure as a result of the collapse of the polypeptide chain, while transient protein-complex formation mainly is a result of binding of two folded individual monomers. Therefore, a protein-protein interface does not resemble the core of monomeric proteins, but has a more polar nature. Here, we address the question of whether the physico-chemical characteristics of intraprotein versus interprotein bonds differ, or whether interfaces are different from folded monomers only in the preference for certain types of interactions. To address this question we assembled a high resolution, nonredundant, protein-protein interaction database consisting of 1374 homodimer and 572 heterodimer complexes, and compared the physico-chemical properties of these interactions between protein interfaces and monomers. We performed extensive statistical analysis of geometrical properties of interatomic interactions of different types: hydrogen bonds, electrostatic interactions, and aromatic interactions. Our study clearly shows that there is no significant difference in the chemistry, geometry, or packing density of individual interactions between interfaces and monomeric structures. However, the distribution of different bonds differs. For example, side-chain-side-chain interactions constitute over 62% of all interprotein interactions, while they make up only 36% of the bonds stabilizing a protein structure. As on average, properties of backbone interactions are different from those of side chains, a quantitative difference is observed. Our findings clearly show that the same knowledge-based potential can be used for protein-binding sites as for protein structures. However, one has to keep in mind the different architecture of the interfaces and their unique bond preference.

Protein-protein interactions networks has come to be a buzzword associated with nets containing edges that represent a pair of interacting proteins (e.g. hormone-receptor, enzyme-inhibitor, antigen-antibody, and a subset of multichain biological machines). Yet, each such interaction composes its own unique network, in which vertices represent amino acid residues, and edges represent atomic contacts. Recent studies have shown that analyses of the data encapsulated in these detailed networks may impact predictions of structure-function correlation. Here, we study homologous families of protein-protein interfaces, which share the same fold but vary in sequence. In this context, we address what properties of the network are shared among relatives with different sequences (and hence different atomic interactions) and which are not. Herein, we develop the general mathematical framework needed to compare the modularity of homologous networks. We then apply this analysis to the structural data of a few interface families, including hemoglobin α-β, growth hormone-receptor, and Serine protease-inhibitor. Our results suggest that interface modularity is an evolutionarily conserved property. Hence, protein-protein interfaces can be clustered down to a few modules, with the boundaries being evolutionarily conserved along homologous complexes. This suggests that protein engineering of protein-protein binding sites may be simplified by varying each module, but retaining the overall modularity of the interface.

Protein-water interactions have long been recognized as a major determinant of chain folding, conformational stability, binding specificity and catalysis. However, the detailed effects of water on stabilizing protein-protein interactions remain elusive. A way to test experimentally the contribution of water-mediated interactions is by applying double mutant cycle analysis on pairs of residues that do not form direct interactions, but are bridged by water. Seven such interactions within the interface between TEM1 and BLIP proteins were evaluated. No significant interaction free energy was found between either of them. Water can bridge interactions, but also stabilize the structure of the monomer. To distinguish between these, we performed a bioinformatic analysis using AQUAPROT (http://bioinfo.weizmann.ac.il/ aquaprot) to determine the degree of water conservation between the bound and unbound states. 29 structures of twelve complexes and 20 related monomers were analyzed. Of the 262 water molecules located within the interfaces, 145 were conserved between the unbound and bound structures. Strikingly, all 50 buried or partially buried waters in the monomer structures were conserved at the same location in the bound structures. Thus, buried waters have an important role in stabilizing the monomer fold rather than contributing to protein-protein binding, and are not replaced by residues from the incoming protein. Taking together the experimental and bioinformatics evidence suggests that exposed waters within the interface may be good sites for protein engineering, while buried or mostly buried waters should be left unchanged.

Site-specific conjugation of proteins to surfaces, spectroscopic probes, or other functional units is a key task for implementing biochemical assays. The streptavidin-biotin interaction has proven a highly versatile tool for detection, quantification, and functional analysis of proteins. We have developed an approach for site-specific reversible biotinylation of recombinant proteins through their histidine tag using biotin conjugated to the multivalent chelator trisnitrilotriacetic acid (^BTtris-NTA). Stable binding of ^BTtris-NTA to His-tagged proteins was demonstrated, which is readily reversed by addition of imidazole, enabling versatile conjugation schemes in solution as well as at interfaces. Gel filtration experiments revealed that His-tagged proteins bind to streptavidin doped with ^BT-tris-NTA in a 2:1 stoichiometry. Furthermore, an increased binding affinity toward His-tagged proteins was observed for ^BTtris-NTA linked to streptavidin compared to tris-NTA in solution and on surfaces. These results indicate an efficient cooperative interaction of two adjacent tris-NTA moieties with a single His-tag, yielding an extremely tight complex with a lifetime of several days. We demonstrate several applications of ^BTtris-NTA including multiplexed capturing of proteins to biosensor surfaces, cell surface labeling, and Western blot detection. The remarkable selectivity of the His-tag-specific biotinylation, as well as the highly stable, yet reversible complex provides the basis for numerous further applications for functional protein analysis.

EmrE is an Escherichia coli H⁺-coupled multidrug transporter that provides a unique experimental paradigm because of its small size and stability, and because its activity can be studied in detergent solution. In this work, we report a study of the transient kinetics of substrate binding and substrate-induced proton release in EmrE. For this purpose, we measured transient changes in the tryptophan fluorescence upon substrate binding and the rates of substrate-induced proton release. The fluorescence of the essential and fully conserved Trp residue at position 63 is sensitive to the occupancy of the binding site with either protons or substrate. The maximal rate of binding to detergent-solubilized EmrE of TPP⁺, a high-affinity substrate, is 2 × 10⁷ M^-1·s^-1, a rate typical of diffusion-limited reactions. Rate measurements with medium- and low-affinity substrates imply that the affinity is determined mainly by the k_off of the substrate. The rates of substrate binding and substrate-induced release of protons are faster at basic pHs and slower at lower pHs. These findings imply that the substrate-binding rates are determined by the generation of the species capable of binding; this is controlled by the high affinity to protons of the glutamate at position 14, because an Asp replacement with a lower pK is faster at the same pHs.

The process of protein-protein association starts with their random collision, which may develop into an encounter complex followed by a transition state and final complex formation. Here we aim to experimentally characterize the nature of the transition state of protein-protein association for three different protein-protein interactions; Barnase-Barstar, TEM1-BLIP and IFNα2-IFNAR2, and use the data to model the transition state structures. To model the transition state, we determined inter-protein distance-constraints of the activated complex by using double mutant cycles (DMC) assuming that interacting residues are spatially close. Significant ΔΔG_int values were obtained only between residues on Barnase and Barstar. However, introducing specific mutations that optimize the charge complementarity between BLIP and TEM1 resulted in the introduction of significant ΔΔG_int values also between residues of these two proteins. While electrostatic interactions make major contributions towards stabilizing the transition state, we show two examples where steric hindrance exerts an effect on the transition state as well. To model the transition-state structures from the experimental ΔΔG_int values, we introduced a method for structure perturbation, searching for those inter-protein orientations that best support the experimental ΔΔG_int values. Two types of transition states were found, specific as observed for Barnase-Barstar and the electrostatically optimized TEM1-BLIP mutants, and diffusive as shown for wild-type TEM1-BLIP and IFNα2-IFNAR2. The specific transition states are characterized by defined inter-protein orientations, which cannot be modeled for the diffusive transition states. Mutations introduced through rational design can change the transition state from diffusive to specific. Together, these data provide a structural view of the mechanism allowing rates of association to differ by five orders of magnitude between different protein complexes.

Interferons (IFNs) stand in the frontline of defense against viral infections. In this study, we aimed at characterizing the gene expression profile specific to the antiviral effect out of the hundreds of genes involved also in other IFN activities. We found that the IFN-induced antiviral state is maintained for a prolonged time even after IFN occlusion. This was achieved through the active expression of a small set of

The development of bioinformatic tools by individual labs results in the abundance of parallel programs for the same task. For example, identification of binding site regions between interacting proteins is done using: ProMate, WHISCY, PPI-Pred, PINUP and others. All servers first identify unique properties of binding sites and then incorporate them into a predictor. Obviously, the resulting prediction would improve if the most suitable parameters from each of those predictors would be incorporated into one server. However, because of the variation in methods and databases, this is currently not feasible. Here, the protein-binding site prediction server is extended into a general protein-binding sites research tool, ProMateus. This web tool, based on ProMate's infrastructure enables the easy exploration and incorporation of new features and databases by the user, providing an evaluation of the benefit of individual features and their combination within a set framework. This transforms the individual research into a community exercise, bringing out the best from all users for optimized predictions. The analysis is demonstrated on a database of protein protein and protein-DNA interactions. This approach is basically different from that used in generating meta-servers. The implications of the open-research approach are discussed. ProMateus is available at http://bip.weizmann.ac.il/promate.

The 3D structures of macromolecules are difficult to grasp and also to communicate. By their nature, movies or animations are particularly useful for highlighting key features by offering a 'guided tour' of structures and conformation changes. However, high-quality movies are rarely seen because they are currently difficult and time consuming to make. By adopting the traditional movie 'storyboard' concept, which gives guidance and direction to filming, eMovie makes the creation of lengthy molecular animations much easier. This tool is a plug-in for the open-source molecular graphics program PyMOL, and enables experts and novices alike to produce informative and high-quality molecular animations.

All α-interferons (IFNα) bind the IFNAR1 receptor subunit with low affinity. Increasing the binding affinity was shown to specifically increase the antiproliferative potency of IFNα2. Here, we constructed a phage display library by randomizing three positions on IFNα2 previously shown to confer weak binding to IFNAR1. The tightest binding variant selected, comprised of mutations H57Y, E58N, and Q61S (YNS), was shown to bind IFNAR1 60-fold tighter compared with wild-type IFNα2, and 3-fold tighter compared with IFNβ. Binding of YNS to IFNAR2 was comparable with wild-type IFNα2. The YNS mutant conferred a 150-fold higher antiproliferative potency in WISH cells compared with wild-type IFNα2, whereas its antiviral activity was increased by only 3.5-fold. The high antiproliferative activity was related to an induction of apoptosis, as demonstrated by annexin V binding assays, and to specific gene induction, particularly TRAIL. To determine the potency of the YNS mutant in a xenograft cancer model, we injected it twice a week to nude mice carrying transplanted MDA231 human breast cancer cells. After 5 weeks, no tumors remained in mice treated with YNS, whereas most mice treated with wild-type IFNα2 showed visible tumors. Histological analysis of these tumors showed a significant antiangiogenic effect of YNS, compared with wild-type IFNα2. This work demonstrates the application of detailed biophysical understanding in the process of protein engineering, yielding an interferon variant with highly increased biological potency.

In a typical cell, proteins function in the crowded cytoplasmic environment where 30% of the space is occupied by macromolecules of varying size and nature. This environment may be simulated in vitro using synthetic polymers. Here, we followed the association and diffusion rates of TEM1-β-lactamase (TEM) and the β-lactamase inhibitor protein (BLIP) in the presence of crowding agents of varying molecular mass, from monomers (ethylene glycol, glycerol, or sucrose) to polymeric agents such as different polyethylene glycols (PEGs, 0.2-8 kDa) and Ficoll. An inverse linear relation was found between translational diffusion of the proteins and viscosity in all solutions tested, in accordance with the Stokes-Einstein (SE) relation. Conversely, no simple relation was found between either rotational diffusion rates or association rates (k_on) and viscosity. To assess the translational diffusion-independent steps along the association pathway, we introduced a new factor, α, which corrects the relative change in k_on by the relative change in solution viscosity, thus measuring the deviations of the association rates from SE behavior. We found that these deviations were related to the three regimes of polymer solutions: dilute, semidilute, and concentrated. In the dilute regime PEGs interfere with TEM-BLIP association by introducing a repulsive force due to solvophobic preferential hydration, which results in slower association than predicted by the SE relation. Crossing over from the dilute to the semidilute regime results in positive deviations from SE behavior, i.e., relatively faster association rates. These can be attributed to the depletion interaction, which results in an effective attraction between the two proteins, winning over the repulsive force. In the concentrated regime, PEGs again dramatically slow down the association between TEM and BLIP, an effect that does not depend on the physical dimensions of PEGs, but rather on their mass concentration. This is probably a manifestation of the monomer-like repulsive depletion effect known to occur in concentrated polymer solutions. As a transition from moderate to high crowding agent concentration can occur in the cellular milieu, this behavior may modulate protein association in vivo, thereby modulating biological function.

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?

Proteins bind one another in aqua's solution to form tight and specific complexes. Previously we have shown that this is achieved through the modular architecture of the interaction network formed by the interface residues, where tight cooperative interactions are found within modules but not between them. Here we extend this study to cover the entire interface of TEM1 β-lactamase and its protein inhibitor BLIP using an improved method for deriving interaction maps based on REDUCE to add hydrogen atoms and then by evaluating the interactions using modifications of the programs PROBE, NCI and PARE. An extensive mutagenesis study of the interface residues indeed showed that each module is energetically independent on other modules, and that cooperativity is found only within a module. By solving the X-ray structure of two interface mutations affecting two different modules, we demonstrated that protein--protein binding occur via the structural reorganization of the binding modules, either by a "lock and key" or an induced fit mechanism. To explain the cooperativity within a module, we performed multiple-mutant cycle analysis of cluster 2 resulting in a high-resolution energy map of this module. Mutant studies are usually done in reference to alanine, which can be regarded as a deletion of a side-chain. However, from a biological perspective, there is a major interest to understand non-Ala substitutions, as they are most common. Using X-ray crystallography and multiple-mutant cycle analysis we demonstrated the added complexity in understanding non-Ala mutations. Here, a double mutation replacing the wild-type Glu,Tyr to Tyr,Asn on TEM1 (res id 104,105) caused a major backbone structural rearrangement of BLIP, changing the composition of two modules but not of other modules within the interface. This shows the robustness of the modular approach, yet demonstrates the complexity of in silico protein design.

Motivation: Secondary structures are key descriptors of a protein fold and its topology. In recent years, they facilitated intensive computational tasks for finding structural homologues, fold prediction and protein design. Their popularity stems from an appealing regularity in patterns of geometry and chemistry. However, the definition of secondary structures is of subjective nature. An unsupervised de-novo discovery of these structures would shed light on their nature, and improve the way we use these structures in algorithms of structural bioinformatics. Methods: We developed a new method for unsupervised partitioning of undirected graphs, based on patterns of small recurring network motifs. Our input was the network of all H-bonds and covalent interactions of protein backbones. This method can be also used for other biological and non-biological networks. Results: In a fully unsupervised manner, and without assuming any explicit prior knowledge, we were able to rediscover the existence of conventional α-helices, parallel β-sheets, anti-parallel sheets and loops, as well as various non-conventional hybrid structures. The relation between connectivity and crystallographic temperature factors establishes the existence of novel secondary structures.

All type I IFNs act through a single cell surface receptor composed of the IFNAR1 and IFNAR2 subunits and two associated cytoplasmic tyrosine kinases of the Janus family, Tyk2 and Jak1. A central issue in type I IFN biology is to understand how 1 a multitude of subtypes can generate similar signaling outputs but also govern specific 1 cellular responses. This review summarizes results from the last decade that contributed to our current state of knowledge of IFN-receptor complex structure and assembly.

The expression, purification, detection, and assay of recombinant proteins have been made more convenient and rapid by the use of small affinity tags. To facilitate the purification of interferon-α2c (IFN-α2c) by metal chelate affinity chromatography, N-terminal 6-histidine tag was introduced via genetic manipulation. Two preparations of IFN material were purified; one contained IFN-α2c with the 6-histidine tag, and the other contained IFN-α2c without the 6-histidine tag. The antigenic properties of the human IFN-α2c subvariant with and without the 6-histidine tag, as well as the effects of the N-terminal 6-histidine tag on IFN-α2c interaction with the extracellular domain of human IFN-α receptor chain 2 (IFNAR2-EC) were examined. For the purposes of this study, IFNs were characterized by Western blots with anti-IFN monoclonal antibodies (mAb) and bioassays. Immunoblot analyses showed differences between IFN-α2c-6-histidine tag and IFN-α2a, b, c in their interaction with IFNAR2-EC. We also observed differences between IFN-α2c-6-histidine tag and IFN-α2a, b, c in bioactivities. This study is the first report that shows that an N-terminal 6-histidine tag on IFN-α2c can affect its interaction with receptor and cause a different bioactivity.

Type I interferons (IFNs) elicit antiviral, antiproliferative and immunomodulatory properties in cells. All of them bind to the same receptor proteins, IFNAR1 and IFNAR2, with different affinities. While the 13 known IFNαs are highly conserved, the C-terminal unstructured tail was found to have large variation in its net charge, from neutral to +4. This led us to speculate that the tail may have a role in modulation of the IFN biological activity, through fine-tuning the binding to IFNAR2. To evaluate this hypothesis, we replaced the tail of IFNα2 with that of IFNα8 and IFNβ tails, or deleted the last five residues of this segment. Mutations to the more positively charged tail of IFNα8 resulted in a 20-fold higher affinity to IFNAR2, which results in a higher antiviral and antiproliferative activity. Double and multiple mutant cycle analysis placed the tail near a negatively charged loop on IFNAR2, comprising of residues Glu 132-134. Deleting the tail resulted in only twofold reduction in binding compared to the wild-type. Next, we modeled the location of the tail using a two-step procedure: first we generated 200 models of the tail docked on IFNAR2 using HADDOCK, second the models were scored according to the fit between experimentally determined rates of association of nine mutant complexes, and their calculated rates using the PARE software. From the results we suggest that the unstructured tail of IFNα is gaining a specific structure in the bound state, binding to a groove below the 132-134 loop in IFNAR2.

Alpha and beta interferons (IFN-α and IFN-β) are multifunctional cytokines that exhibit differential activities through a common receptor composed of the subunits IFNAR1 and IFNAR2. Here we combined biophysical and functional studies to explore the mechanism that allows the alpha and beta IFNs to act differentially. For this purpose, we have engineered an IFN-α2 triple mutant termed the HEQ mutant that mimics the biological properties of IFN-β. Compared to wild-type (wt) IFN-α2, the HEQ mutant confers a 30-fold higher binding affinity towards IFNAR1, comparable to that measured for IFN-β, resulting in a much higher stability of the ternary complex as measured on model membranes. The HEQ mutant, like IFN-β, promotes a differentially higher antiproliferative effect than antiviral activity. Both bring on a down-regulation of the IFNAR2 receptor upon induction, confirming an increased ternary complex stability of the plasma membrane. Oligonucleotide microarray experiments showed similar gene transcription profiles induced by the HEQ mutant and IFN-β and higher levels of gene induction or repression than those for wt IFN-α2. Thus, we show that the differential activities of IFN-β are directly related to the binding affinity for IFNAR1. Conservation of the residues mutated in the HEQ mutant within IFN-α subtypes suggests that IFN-α has evolved to bind IFNAR1 weakly, apparently to sustain differential levels of biological activities compared to those induced by IFN-β.

De novo design and redesign of proteins and protein complexes have made promising progress in recent years. Here, we give an overview of how to use available computer-based tools to design proteins to bind faster and tighter to their protein-complex partner by electrostatic optimization between the two proteins. Electrostatic optimization is possible because of the simple relation between the Debye-Huckel energy of interaction between a pair of proteins and their rate of association. This can be used for rapid, structure-based calculations of the electrostatic attraction between the two proteins in the complex. Using these principles, we developed two computer programs that predict the change in k(on), and as such the affinity, on introducing charged mutations. The two programs have a web interface that is available at www.weizmann.ac.il/home/bcges/PARE.html and http://bip.weizmann.ac.il/hypare. When mutations leading to charge optimization are introduced outside the physical binding site, the rate of dissociation is unchanged and therefore the change in k(on) parallels that of the affinity. This design method was evaluated on a number of different protein complexes resulting in binding rates and affinities of hundreds of fold faster and tighter compared to wild type. In this chapter, we demonstrate the procedure and go step by step over the methodology of using these programs for protein-association design. Finally, the way to easily implement the principle of electrostatic design for any protein complex of choice is shown.

In this issue of Structure, Cho et al. (2005) provide a detailed structural analysis of the molecular mechanism of affinity maturation for the murine T cell receptor Vβ domain. Information from nine X-ray structures of combinations along the maturation pathway and comprehensive thermodynamics reveal the fine details of high-affinity binding.

The association of two proteins is preceded by a mutual diffusional search in solution. The role of translational and rotational diffusion in this process has been studied theoretically for many years. However, systematic experimental verification of theoretical results is still lacking. We report here measurements of association rates of the proteins β-lactamase (TEM) and β-lactamase inhibitor protein (BLIP) in solutions of glycerol and poly(ethylene glycol) of increasing viscosity. We also measured translational and rotational diffusion in the same solutions, using fluorescence correlation spectroscopy and fluorescence anisotropy, respectively. It is found that in glycerol both translational and rotational diffusion rates are inversely dependent on viscosity, as predicted by the classical Stokes-Einstein relations, while the association rate depends nonlinearly on viscosity. In contrast, the association rate depends only weakly on the viscosity of the polymer solutions, which results in a similar weak dependence of k_on on viscosity. The data are modeled using the theory of diffusion-limited association. Deviations from the theory are explained by a short-range solute-induced repulsion between the proteins in glycerol solution and an attractive depletion interaction generated by the polymers. These results open the way to the creation of a unified framework for all nonspecific effects involved in the protein association process, as well as to better theoretical understanding of these effects. Further, they reflect on the complex factors controlling protein association within the crowded environment of cells and suggest that a high concentration of macromolecules does not significantly impede protein association.

Type I interferons activate cellular responses by forming a ternary complex with two receptor components, IFNAR1 and IFNAR2. While the binding of the IFNAR2 receptor to interferon is of high affinity and well characterized, the binding to IFNAR1 is weak, transient, and poorly understood. Here, we mapped the complete binding region of IFNAR1 on IFNα2 by creating a panel of 21 single alanine mutant proteins, and determined their binding affinities. The IFNAR1 binding site on IFNα2 maps to the center of the B and C helices, opposite to the binding site for IFNAR2. No hot spots for binding were found in the interface, with individual mutations having an up to fivefold effect on binding. Of the nine residues that affected binding, three adjacent conserved residues, located on the B helix, conferred an increase in the binding affinity to IFNAR1, as well as an increase in the biological activity of the interferon mutant. This suggests that binding of alpha interferons to the IFNAR1 receptor is sub-optimal. A correlation between binding affinity and biological activity was found, albeit not across the whole range of affinities. In WISH cells, but not DAUDI cells, the anti-proliferative activity was markedly affected by fluctuations in the IFNα2 affinity towards the IFNAR1 receptor. On the other hand, the antiviral activity of interferons on WISH cells seems to change in accordance to the binding affinity towards IFNAR1 only as long as the binding affinity is not beyond twofold of the wild-type. In accordance, the biological roles of the two interferon-receptor subunits are discussed.

The principal goal of the Israel Structural Proteomics Center (ISPC) is to determine the structures of proteins related to human health in their functional context. Emphasis is on the solution of structures of proteins complexed with their natural partner proteins and/or with DNA. To date, the ISPC has solved the structures of 14 proteins, including two protein complexes. It has adopted automated high-throughput (HTP) cloning and expression techniques and is now expressing in Escherichia coli, Pichia pastoris and baculovirus, and in a cell-free E. coli system. Protein expression in E. coli is the primary system of choice in which different parameters are tested in parallel. Much effort is being devoted to development of automated refolding of proteins expressed as inclusion bodies in E. coli. The current procedure utilizes tagged proteins from which the tag can subsequently be removed by TEV protease, thus permitting streamlined purification of a large number of samples. Robotic protein crystallization screens and optimization utilize both the batch method under oil and vapour diffusion. In order to record and organize the data accumulated by the ISPC, a laboratory information-management system (LIMS) has been developed which facilitates data monitoring and analysis. This permits optimization of conditions at all stages of protein production and structure determination. A set of bioinformatics tools, which are implemented in our LIMS, is utilized to analyze each target.

The rate of association of a protein complex is a function of an intrinsic basal rate and of the magnitude of electrostatic steering. In the present study we analyze the contribution of electrostatics towards the association rate of proteins in a database of 68 transient hetero-protein-protein complexes. Our calculations are based on an upgraded version of the computer algorithm PARE, which was shown to successfully predict the impact of mutations on k _on by calculating the difference in Columbic energy of interaction of a pair of proteins. HyPare (http://bip.weizmann.ac.il/HyPare), automatically calculates the impact of mutations on a perresidue basis for all residues of a protein-protein interaction, achieving a precision similar to that of PARE. Our calculations show that electrostatics play a marginal role (

Protein-protein interactions are essential for life. Yet, our understanding of the general principles governing binding is not complete. In the present study, we show that the interface between proteins is built in a modular fashion; each module is comprised of a number of closely interacting residues, with few interactions between the modules. The boundaries between modules are defined by clustering the contact map of the interface. We show that mutations in one module do not affect residues located in a neighboring module. As a result, the structural and energetic consequences of the deletion of entire modules are surprisingly small. To the contrary, within their module, mutations cause complex energetic and structural consequences. Experimentally, this phenomenon is shown on the interaction between TEM1-beta-lactamase and beta-lactamase inhibitor protein (BLIP) by using multiple-mutant analysis and x-ray crystallography. Replacing an entire module of five interface residues with Ala created a large cavity in the interface, with no effect on the detailed structure of the remaining interface. The modular architecture of binding sites, which resembles human engineering design, greatly simplifies the design of new protein interactions and provides a feasible view of how these interactions evolved.

Many peptide and protein drugs have a short circulatory half-life in vivo. The covalent attachment of polyethylene glycol (PEG) chains (PEGylation) can overcome this deficiency, but pegylated peptides and proteins are often inactive. In this study, we present a novel PEG-IFNα2 conjugate, PEG ₄₀-FMS-IFNα2, capable of regenerating native interferon α2 (IFNα2) at a slow rate under physiological conditions. A 2-sulfo-9-fluorenylmethoxycarbonyl (FMS) containing bifunctional reagent, MAL-FMS-NHS, has been synthesized, enabling the linkage of a 40 kDa PEG-SH to IFNα2 through a slowly hydrolyzable bond. By use of a BIAcore binding assay, the in vitro rate of regeneration of native interferon was estimated to have a half-life of 65 h. Following subcutaneous administration to rats and monitoring circulating antiviral activity, active IFNα2 levels peaked at 50 h, with substantial levels still being detected 200 h after administration. This value contrasts with a half-life of about 1 h measured for unmodified interferon. The concentration of active IFNα2 scaled linearly with the quantity injected. Comparing subcutaneous to intravenous administration of PEG₄₀-FMS-IFNα2, we found that the long circulatory lifetime of IFNα2 was affected both by the slow rate of absorption of the PEGylated protein from the subcutaneous volume and by the slow rate of discharge from the PEG in circulation. A numerical simulation of the results was in good agreement with the results observed in vivo. The pharmacokinetic profile of this novel IFNα2 conjugate combines a prolonged maintenance in vivo with the regeneration of active-native IFNα2, ensuring ready access to peripheral tissues and thus an overall advantage over currently used formulations.

Association of two proteins can be described as a two-step process, with the formation of an encounter complex followed by desolvation and establishment of a tight complex. Here, by using the computer algorithm PARE, we designed a set of mutants of the Ras effector protein Ral guanine nucleotide dissociation stimulator (RalGDS) with optimized electrostatic steering. The fastest binding RalGDS mutant, M26K,D47K,E54K, binds Ras 14-fold faster and 25-fold tighter compared with WT. A linear correlation was found between the calculated and experimental data, with a correlation coefficient of 0.97 and a slope of 0.65 for the 24 mutants produced. The data suggest that increased electrostatic steering specifically stabilizes the encounter complex and transition state. This conclusion is backed up by Φ analysis of the encounter complex and transition state of the RalGDS^{M26K,D47K,E54K}/Ras complex, with both values being close to 1. Upon further formation of the final complex, the increased Coulombic interactions are probably counterbalanced by the cost of desolvation of charges, keeping the dissociation rate constant almost unchanged. This mechanism is also reflected by the mutual compensation of enthalpy and entropy changes quantified by isothermal titration calorimetry. The binding constants of the faster binding RalGDS mutants toward Ras are similar to those of Raf, the most prominent Ras effector, suggesting that the design methodology may be used to switch between signal transduction pathways.

Prolonging the circulatory half-life of low mass protein drugs can be achieved either by administration of a pro-drug or through co-injection with a carrier protein that will slowly release the active protein. The rate of release is concentration and affinity dependent. The optimal relationship between these two in prolonging the half-life of a pro-drug is the focus of this work. Interferon (IFN) β is one of the most widely used protein drugs in the clinic. Here, we show that the circulatory half-life of IFNβ can be significantly extended by co-administration with the extracellular domain of the IFN receptor ifnar2 (ifnar2-EC). To investigate the concentration/affinity relation, a range of tighter binding ifnar2-EC mutants was designed that bind IFNβ, but not IFNα2, up to 50-fold tighter compared with the wild-type ifnar2-EC. This increased affinity is related to a slower dissociation rate, whereas the association of IFNβ with ifnar2-EC is already near optimum. Using the wild-type and mutant receptors, we investigated their potential in occluding IFNβ from circulation in a tissue culture assay, as well as in rats. To determine the potential of ifnar2-EC as a carrier protein, we co-administered a mixture of IFNβ and ifnar2-EC to rats both intravenously and subcutaneously, and followed the blood plasma concentrations of IFNβ over time. The tighter binding ifnar2-EC mutant had a clear advantage in prolonging the half-life of IFNβ in circulation, even when lower protein concentrations were administered. A numerical simulation of the in vivo data demonstrates that the optimal binding affinity of a carrier protein is around the concentration needed to obtain optimal activity of the ligand.

Is the whole protein surface available for interaction with other proteins, or are specific sites pre-assigned according to their biophysical and structural character? And if so, is it possible to predict the location of the binding site from the surface properties? These questions are answered quantitatively by probing the surfaces of proteins using spheres of radius of 10Å on a database (DB) of 57 unique, non-homologous proteins involved in heteromeric, transient protein-protein interactions for which the structures of both the unbound and bound states were determined. In structural terms, we found the binding site to have a preference for β-sheets and for relatively long non-structured chains, but not for α-helices. Chemically, aromatic side-chains show a clear preference for binding sites. While the hydrophobic and polar content of the interface is similar to the rest of the surface, hydrophobic and polar residues tend to cluster in interfaces. In the crystal, the binding site has more bound water molecules surrounding it, and a lower B-factor already in the unbound protein. The same biophysical properties were found to hold for the unbound and bound DBs. All the significant interface properties were combined into ProMate, an interface prediction program. This was followed by an optimization step to choose the best combination of properties, as many of them are correlated. During optimization and prediction, the tested proteins were not used for data collection, to avoid over-fitting. The prediction algorithm is fully automated, and is used to predict the location of potential binding sites on unbound proteins with known structures. The algorithm is able to successfully predict the location of the interface for about 70% of the proteins. The success rate of the predictor was equal whether applied on the unbound DB or on the disjoint bound DB. A prediction is assumed correct if over half of the predicted continuous interface patch is indeed interface. The ability to predict the location of protein-protein interfaces has far reaching implications both towards our understanding of specificity and kinetics of binding, as well as in assisting in the analysis of the proteome.

Physiological media constitutes a crowded environment that serves as the field of action for protein-protein interaction in vivo. Measuring protein-protein interaction in crowded solutions can mimic this environment. In this work we follow the process of protein-protein association and its rate constants (k_on) of the β-lactamase (TEM)-β-lactamase inhibitor protein (BLIP) complex in crowded solution using both low and high molecular mass crowding agents. In all crowded solutions (0-40% (w/w) of ethylene glycol (EG), poly(ethylene glycol) (PEG) 200, 1000, 3350, 8000 Da Ficoll-70 and Haemaccel the measured absolute k_on, but not k _off values, were found to be slower as compared to buffer. However, there is a fundamental difference between low and high mass crowding agents. In the presence of low mass crowding agents and Haemaccel k_on depends inversely on the solution viscosity. In high mass polymer solutions k _on changes only slightly, even at viscosities 12-fold higher than water. The border between low and high molecular mass polymers is sharp and is dictated by the ratio between the polymer length (L) and its persistence length (L_p). Polymers that are long enough to form a flexible coil (L/L _p>2) behave as high molecular mass polymers and those who are unable to do so (L/L_p

Docking algorithms produce many possible structures of a protein-protein complex. In most cases some of them resemble the correct structure within an r.m.s.d. of

The potent antiviral and antiproliferative activities of human type I interferons (IFNs) are mediated by a single receptor comprising two subunits, IFNAR1 and IFNAR2. The structure of the IFNAR2 IFN binding ectodomain (IFNAR2-EC), the first helical cytokine receptor structure determined in solution, reveals the molecular basis for IFN binding. The atypical perpendicular orientation of its two fibronectin domains explains the lack of C domain involvement in ligand binding. A model of the IFNAR2-EC/IFNα2 complex based on double mutant cycle-derived constraints uncovers an extensive and predominantly aliphatic hydrophobic patch on the receptor that interacts with a matching hydrophobic surface of IFNα2. An adjacent motif of alternating charged side chains guides the two proteins into a tight complex. The binding interface may account for crossreactivity and ligand specificity of the receptor. This molecular description of IFN binding should be invaluable for study and design of IFN-based biomedical agents.

[No abstract available]

The human interferon receptor (IFNAR) mediates the antiviral and antiproliferative activities of type I interferons (IFNs). This receptor is comprised of subunits IFNAR1 and IFNAR2, the latter exhibiting nanomolar affinity for IFNs. Here the extracellular domain of IFNAR2 (IFNAR2-EC), a soluble 25 kDa IFN-binding polypeptide, and its complex with IFN-α2 were studied using multidimensional NMR. IFNAR2-EC is comprised of two fibronectin-III (FN-III) domains connected by a helical hinge region. The deduced global fold was utilized to improve the alignment of IFNAR2-EC against structurally related receptors and to model its structure. A striking feature of IFNAR2-EC is the limited and localized deviations in chemical shifts exhibited upon ligand binding, observed for only 15% of its backbone ¹H and ¹⁵N nuclei. Analysis of these deviations maps the IFN-α2 binding site upon IFNAR2-EC to a contiguous surface on the N-terminal domain, including the S3-S4 loop (residues 44-53), the S5-S6 loop and S6 β-strand (residues 74-82), and the S7 β-strand and the hinge region (residues 95-105). The C-terminal domain contributes only marginally to ligand binding, and no change in the hypothesized interdomain interface is observed. The proposed binding domain encompasses all residues implicated by mutagenesis studies in IFN binding, and suggests adjacent residues cooperate in forming the binding surface. D₂O-exchange experiments indicate that binding of IFN-α2 induces tightening of the N-terminal domain of IFNAR2-EC. This increase in receptor rigidity may play an important role in initiating the intracellular stage of the IFN signaling cascade.

The structure of a protein-protein interaction, its affinity and thermodynamic characteristics depict a 'frozen' state of a complex. This picture ignores the kinetic nature of complex formation and dissociation, which are of major biological and biophysical interest. This review highlights recent advances in deciphering the kinetic pathway of protein-protein complexation, the nature of the encounter complex, transition state and intermediate along the reaction, and the effects of mutation, viscosity, pH and salt on association.

The pleiotropic activity of type I interferons has been attributed to the specific interaction of IFN with the cell-surface receptor components ifnarl and ifnar2. To date, the structure of IFN has been solved, but not that of the receptor or the complex. In this study, the structure of the IFN-α2-ifnar2 complex was generated with a docking procedure, using nuclear Overhauser effect-like distance constraints obtained from double-mutant cycle experiments. The interaction free energy between 13 residues of the ligand and 11 of the receptor was measured by double-mutant cycles. Of the 100 pairwise interactions probed, five pairs of residues were found to interact. These five interactions were incorporated as distance constraints into the flexible docking program PRODOCK by using fixed and movable energy-gradient grids attached to the receptor and ligand, respectively. Multistart minimization and Monte Carlo minimization docking of IFN-α2 onto ifnar2 converged to a well-defined average structure, with the five distance constraints being satisfied. Furthermore, no structural artifacts or intraloop energy strain were observed. The mutual binding sites on IFN-α2 and ifnar2 predicted from the model showed an almost complete superposition with the ones determined from mutagenesis studies. Based on this structure, differences in IFN-α2 versus IFN-β binding are discussed.

Association of a protein complex follows a two-step mechanism, with the first step being the formation of an encounter complex that evolves into the final complex. Here, we analyze recent experimental data of the association of TEM1-β-lactamase with BLIP using theoretical calculations and simulation. We show that the calculated Debye-Hückel energy of interaction for a pair of proteins during association resembles an energy funnel, with the final complex at the minima. All attraction is lost at inter-protein distances of 20 Å, or rotation angles of >60° from the orientation of the final complex. For faster-associating protein complexes, the energy funnel deepens and its volume increases. Mutations with the largest impact on association (hotspots for association) have the largest effect on the size and depth of the energy funnel. Analyzing existing evidence, we suggest that the transition state along the association pathway is the formation of the final complex from the encounter complex. Consequently, pairs of proteins forming an encounter complex will tend to dissociate more readily than to evolve into the final complex. Increasing directional diffusion by increasing favorable electrostatic attraction results in a faster forming and slower dissociating encounter complex. The possible applicability of electrostatic calculations for protein-protein docking is discussed.

Association of a protein complex follows a two step reaction mechanism, with the first step being the formation of an encounter complex which evolves into the final complex. Here we present new experimental data for the association of the bacterial ribonuclease barnase and its polypeptide inhibitor barstar which shed light on the thermodynamics and structure of the transition state and preceding encounter complex of association at diminishing electrostatic attraction. We show that the activation entropy at the transition state is close to zero, with the activation enthalpy being equal to the free energy of binding. This observation was independent of the magnitude of the mutual electrostatic attraction, which were altered by mutagenesis or by addition of salt. The low activation entropy implies that the transition state is mostly solvated at all ionic strengths. The structure of the transition state was probed by measuring pairwise interaction energies using double-mutant-cycles. While at low ionic strength all proximal charge-pairs form contacts, at high salt only a subset of these interactions are maintained. More specifically, charge-charge interactions between partially buried residues are lost, while exposed charged residues maintain their ability to form specific interactions even at the highest salt concentration. Uncharged residues do not interact at any ionic strength. The results presented here suggest that the barnase-barstar binding sites are correctly aligned during the transition state even at diminishing electrostatic attraction, although specific short range interactions of uncharged residues are not yet formed. Furthermore, most of the interface desolvation (which contributes to the entropy of the system) has not yet occurred. This picture seems to be valid at low and high salt. However, at high salt, interactions of the activated complex are limited to a more restricted set of residues which are easier approached during diffusion, prior to final docking. This suggest that the steering region at high salt is more limited, albeit maintaining its specificity.

Investigating protein-protein interactions by mutational analysis requires practical techniques for quantifying rate constants and equilibrium constants over several orders of magnitude with reasonably high sample throughput. We have employed spectroscopic interferometry for label-free monitoring of the interaction between the cytokine interferon alpha2 (IFN alpha2) and the extracellular domain of its receptor ifnar2 (ifnar2-EC), We implemented a versatile surface chemistry for the glass substrate of this transducer for covalent immobilization of proteins. Affinity capturing with a monoclonal anti-ifnar2-EC antibody (mAb) followed by crosslinking with a second, noncompetitive mAb provided stable, but still reversible, immobilization of ifnar2-EC, We measured kinetics and affinity of numerous of mutants of IFN alpha2 and ifnar2-EC. Dissociation rate constants up to 0.3 s(-1) and association rate constants up to 3 x 10(6) M(-)1 s(-1) were resolved by the system. Dissociation constants down to 200 muM were measured with protein concentrations up to 50 muM without no background signal or nonspecific binding. The instrument detection limit is similar to 10 pm without the need for temperature stabilization or referencing channels. The system proved effective for large-scale mutational analysis involving alanine scanning mutagenesis and double mutant cycles. (C) 2001 Academic Press.

Polypeptide drugs are generally short-lived species in circulation. In this study, we have covalently linked seven moieties of 2-sulfo-9-fluorenylmethoxycarbonyl (FMS) to the amino groups of human interferon-α2. The derivative thus obtained (FMS₇-IFN-α2) has ≈4% the biological potency and 33 ± 4% the receptor binding capacity of the native cytokine. Upon incubation, FMS₇-IFN-α2 undergoes time-dependent spontaneous hydrolysis, generating active interferon with t_1/2 values of 24 ± 2 h at pH 8.5 and 98 ± 10 h at pH 7.4. When native IFN-α2 is intravenously administered to mice, circulating antiviral activity is maintained for a short duration and then declines with t_1/2 = 4 ± 0.5 h, reaching undetectable values at ≈18 h after administration. With intravenously administered FMS₇-IFN-α2, there is a lag period of 2 h, followed by a progressive elevation in circulating antiviral-active protein, which peaked at 20 h and declined with t_1/2 = 35 ± 4 h. FMS₇-IFN-α2 is resistant to α-chymotrypsin digest and to proteolytic inactivation by human serum proteases in vitro. We have thus introduced here an inactive IFN-α2 derivative, which is resistant to in situ inactivation and has the capability of slowly reverting to the native active protein at physiological conditions in vivo and in vitro. Having these attributes, FMS₇-IFN-α2 maintains prolonged circulating antiviral activity in mice, exceeding 7-8 times the activity of intravenously administered native cytokine.

Type I interferons bind to two cell surface receptors, ifnar1 and ifnar2, as the first step in the activation of several signal transduction pathways that elicit an antiviral state and an anti-proliferative response. Here, we quantitatively mapped the complete binding region of ifnar2 on interferon (IFN)α2 by 35 individual mutations to alanine and isosteric residues. Of the six "hot-spot" residues identified (Leu-30, Arg-33, Arg-144, Ala-145, Met-148, and Arg-149), four are located on the E-helix, which is located at the center of the binding site flanked by residues on the A-helix and the AB-loop. The contribution of residues of the D-helix, which have been previously implicated in binding, proved to be marginal for the interaction with the extracellular domain of ifnar2. Interestingly, the ifnar2 binding site overlaps the largest continuous hydrophobic patch on IFNα2. Thus, hydrophobic interactions seem to play a significant role stabilizing this interaction, with the charged residues contributing toward the rapid association of the complex. Relating the anti-viral and anti-proliferative activity of the various interferon mutants with their affinity toward ifnar2 results in linear function over the whole range of affinities investigated, suggesting that ifnar2 binding is the rate-determining step in cellular activation. Dose-time analysis of the anti-viral response revealed that shortening the incubation time of low-level activation cannot be compensated by higher IFN doses. Considering the strict dependence of the cellular response on affinity, these results suggest that for maintaining transcription of IFN-responsive genes over a longer time period, low but continuous signaling through the IFN receptor is essential.

A protein design strategy was developed to specifically enhance the rate of association (k(on)) between a pair of proteins without affecting the rate of dissociation (k(off)). The method is based on increasing the electrostatic attraction between the proteins by incorporating charged residues in the vicinity of the binding interface. The contribution of mutations towards the rate of association was calculated using a newly developed computer algorithm, which predicted accurately the rate of association of mutant protein complexes relative to the wild type. Using this design strategy, the rate of association and the affinity between TEM 1 β-lactamase and its protein inhibitor BLIP was enhanced 250-fold, while the dissociation rate constant was unchanged. The results emphasize that long range electrostatic forces specifically alter k(on), but do not effect k(off). The design strategy presented here is applicable for increasing rates of association and affinities of protein complexes in general.

An experimental approach to evaluate the net binding free energy of buried hydrogen bonds and salt bridges is presented. The approach, which involves a modified multiple-mutant cycle protocol, was applied to selected interactions between TEM-1-β-lactamase and its protein inhibitor, BLIP. The selected interactions (two salt bridges and two hydrogen bonds) all involving BLIP-D49, define a distinct binding unit. The penta mutant, where all side-chains constructing the binding unit were mutated to Ala, was used as a reference state to which combinations of side-chains were introduced. At first, pairs of interacting residues were added allowing the determination of interaction energies in the absence of neighbors, using double mutant cycles. Addition of neighboring residues allowed the evaluation of their cooperative effects on the interaction. The two isolated salt bridges were either neutral or repulsive whereas the two hydrogen bonds contribute 0.3 kcal mol^-1 each. Conversely, a double mutant cycle analysis of these interactions in their native environment showed that they all stabilize the complex by 1-1.5 kcal mol^-1. Examination of the effects of neighboring residues on each of the interactions revealed that the formation of a salt bridge triad, which involves two connected salt bridges, had a strong cooperative effect on stabilizing the complex independent of the presence or absence of additional neighbors. These results demonstrate the importance of forming networks of buried salt bridges. We present theoretical electrostatic calculations which predict the observed mode of cooperativity, and suggest that the cooperative networking effect results from the favorable contribution of the protein to the interaction. Furthermore, a good correlation between calculated and experimentally determined interaction energies for the two salt bridges, and to a lesser extent for the two hydrogen bonds, is shown. The data analysis was performed on values of ΔΔG(Kd)/(+)) which reflect the strength of short range interactions, while ΔΔG((O)/(KD)) values which include the effects of long range electrostatic forces that alter specifically ΔΔG((Ka)/(+)) were treated separately. (C) 2000 Academic Press.

Type I interferons (IFN) exert pleiotropic activities through binding to two cell surface receptors, ifnar1 and ifnar2. We are investigating the biophysical basis of IFN signaling by characterizing the complex of the extra-cellular domain of ifnar2 (ifnar2-EC) with IFNs on the level of purified recombinant proteins in vitro. Here, we present a detailed mutational study on the functional epitopes on both IFN and ifnar2. Kinetic and thermodynamic parameters were determined by label-free heterogeneous phase detection. On IFNα2, a relatively small functional epitope comprising ten amino acid residues was localized, which is nearly entirely formed by residues on the AB loop. Two hot-spot residues, L30 and R33, account for two-thirds of the total interaction energy. Comparing the anti-viral potency of the various mutants to the binding affinity towards ifnar2 revealed a proportional correlation between the true, suggesting a rate-limiting role of ifnar2 binding in IFN signaling. On ifnar2, residues T46, I47 and M48 were identified as hot-spots in the interaction with IFNα2. For another ten residues on ifnar2, significant contribution of interaction energy was determined. Based on these data, the functional epitope on ifnar2 was defined according to a homology model based on other members of the class II hCR family in good agreement with the complementary functional epitope on IFNα2. Although IFNα2 and IFNβ bind competitively to the same functional epitope, mutational analysis revealed distinct centers of binding for these IFNs on ifnar2. This small shift of the binding site may result in different angular orientation, which can be critically coupled to cytoplasmic signaling.

Type I interferons are cytokines which activate an anti-viral response by binding to two specific cell surface receptors, ifnar1 and ifnar2. Here, we report purification and refolding of the extracellular part of human ifnar2 (ifnar2-EC) expressed in Escherichia coli and its characterization with respect to its interaction with interferon α2 (IFNα2). The 25 kDa, non-glycosylated ifnar2-EC is a stable, fully active protein, which inhibits antiviral activity of IFNα2. The stoichiometry of binding IFNα2 is 1:1, as determined by gel filtration, chemical cross-linking and solid-phase detection. The affinity of this interaction is 10 nM, which is similar to the affinity measured for the cell surface-bound ifnar2 receptor. No difference in affinity was found throughout various assays using optical detection as BIAcore or reflectometric interference spectorscopy. However, the binding kinetics as measured in homogeneous phase by fluorescence dequenching was about three times faster than that measured on a sensor surface. The rate of complex formation is relatively high compared to other cytokine-receptor interactions. The salt dependence of the association kinetics suggest a limited but significant contribution of electrostatic forces towards the rate of complex formation. The dissociation constant increases with decreasing pH according to the protonation of a base with a pK(a) of 6.7. The surface properties of the IFNα2 binding surface on ifnar2 were interpreted according to the pH and salt dependence of the interaction.

The rate of association of proteins is dictated by diffusion, but can be enhanced by favorable electrostatic forces. Here the relationship between the electrostatic energy of interaction, and the kinetics of protein-complex formation was analyzed for the protein pairs of: hirudin-thrombin, acetyl-cholinesterase-fasciculin and barnase-barstar, and for a panel of point mutants of these proteins. Electrostatic energies of interaction were calculated as the difference between the electrostatic energy of the complex and the sum of the energies of the two individual proteins, using the computer simulation package DelPhi. Calculated electrostatic energies of interaction were compared to experimentally determined rates of association. One kcal/mol of Coulombic interaction energy increased the rate of association by a factor of 2.8, independent of the protein-complex or mutant analyzed. Electrostatic energies of interaction were also determined from the salt dependence of the association rate constant, using the same basic equation as for the theoretical calculation. A Bronsted analysis of the electrostatic energies of interactions plotted versus experimentally determined ln(rate)s of association shows a linear relation between the two, with a β value close to 1. This is interpreted as the energy of the transition state varies according to the electrostatic interaction energy, fitting a two state model for the association reaction. Calculating electrostatic rate enhancement from the electrostatic interaction energy can be used as a powerful tool to design protein complexes with altered rates of association and affinities.

BLIP is a secreted protein from Streptomyces clavuligerus that inhibits a wide range of β-lactamases. Here we investigate the tight interaction of BLIP, expressed heterologousely in E. coli, with TEM-1. Kinetic and thermodynamic constants were determined using methods with the proteins either in a homogeneous or in a heterogeneous phase. While values of ΔΔG((mut-wt)) are similar whether measured by fluorescence quench, enzyme inhibition, or surface plasmon resonance, absolute values of ΔG and kinetic constants vary. Association and dissociation rate constants of 10⁵ M^-1 s^{- 1} and 10^-4 s^-1, respectively, and a nanomolar affinity were determined for the wild-type proteins. The highest affinity is measured at pH 7.5, with a decreasing association rate constant at higher pH values, and an increasing dissociation rate constant at lower pH values. The marginal effect of salt on the kinetics of binding, as well as the calculated surface potentials, suggests a limited role for electrostatic forces in guiding this reaction. Still, mutations of interfacial residues affect the rate of association significantly, so that an increase in the net negative charge on either protein reduces the association rate constant. We show that simple electrostatic rules can explain this behavior. BLIP inhibits the catalytic activity of TEM-1 by binding its active site. Yet, mutations of active site residues on TEM-1 only have a moderate though cooperative effect on the binding energy. This can be explained in light of the peripheral location of the active site in the interface between the two proteins.

The electrostatic enhancement of the association rate of barnase and barstar is calculated using a transition-state theory like expression and atomic-detail modeling of the protein molecules. This expression predicts that the rate enhancement is simply the average Boltzmann factor in the region of configurational space where association occurs instantaneously in the diffusion-controlled limit. Based on experimental evidence, this 'transition state' is defined by configurations in which, relative to the stereospecifically bound complex, the two proteins are shifted apart by ~8 Å (so a layer of water can be accommodated in the interface) and the two binding surfaces are rotated away by 0°to 3°. The values of the average Boltzmann factor, calculated by solving the Poisson-Boltzmann equation, for the wild-type complex and 16 complexes with single mutations are found to correlate well with experimental results for the electrostatic rate enhancement. The predicted rate enhancement is found to be somewhat insensitive to the precise definition of the transition state, due to the long-range nature of electrostatic interactions. The experimental ionic strength dependence of the rate enhancement is also reasonably reproduced.

Barnase, a small extracellular ribonuclease from BacilIus amyloliquefaciens and its intracellular inhibitor barstar have co-evolved to bind tightly and rapidly. Barnase has also evolved to be catalytically active. The active site of barnase and its binding site for barstar use the same subset of amino acids. The exception is Glu73 (the general base in catalysis), which although located at the centre of the binding site, is separated by three ordered water molecules from barstar. We examined in this work the contribution of Glu73 to both catalysis and barstar binding. Truncation mutants of the general base (Glu73 → Ala or Ser) retain a residual RNase activity of about 0.3% while mutants with larger hydrophobic replacements (Glu 73 → Trp or Phe) have virtually no catalytic activity. This, and binding data of 3'-GMP with the different barnase mutants suggest that the loss in activity results from the elimination of the general base, which can be substituted to some extent by water or other polar side-chains in truncation mutants. All of the Glu73 mutations lead to a weakening of the free energy of complex formation with barstar by 1.4 to 3.0 kcal/mol (including Gln). This is surprising, since Glu73 does not interact directly with barstar and there is an electrostatic repulsion between Glu73 on barnase and the negatively charged binding surface of barstar. A newly developed method of constructing double mutant cycles between multiple mutations at the same site appears to pinpoint a favourable interaction between Glu73 and one of its nearest neighbours in barstar, Asp39. The coupling energy between those residues is presumably indirect: the carboxylate of Glu73 organizes neighbouring positively charged groups in barnase, Lys27, Arg83, and Arg87 to interact with Asp39 in barstar. This emphasizes that an apparent interaction between a pair of residues as measured with double mutant cycles is the sum of their direct and indirect interactions.

We have studied the thermodynamics of the interaction between the ribonuclease barnase and its natural polypeptide inhibitor barstar. The contribution of specific residues and interactions within the barnase-barstar interface to the enthalpy of binding has been examined using isothermal titration calorimetry and protein engineering. The enthalpy of association of the wild-type proteins is -18.9 (± 0.1) kcal/mol at pH 8 and at 25°C. The enthalpy of binding remains favourable for 31 different combinations of mutations in the interface. The effects on the binding enthalpy upon replacing a side-chain involved in the interaction of barnase and barstar are, however, always unfavourable and in most cases larger than the effects on the free energy of binding. Interaction enthalpies calculated by double mutant cycle analysis are in some cases much larger than the interaction free energies. The interaction enthalpies for complexes between different barnase mutants with amino acid substitutions of the general base residue glutamic acid 73 and a barstar variant (D39A) vary by as much as 8.3 kcal/mol while the coupling free energies differ only by 1 kcal/mol. The use of enthalpies for the analysis of structure-activity relationships appears to be complicated by enthalpy-entropy compensation of weak intermolecular interactions. These tend to cancel out in measurements of free energy, which is thus the preferred quantity for simple analysis of interactions.

We have documented the folding pathway of the 10-kDa protein barstar from the first few microseconds at the resolution of individual residues from its well characterized denatured state. The denatured state had been shown from NMR to have flickering native-like structure in the first two of its four α-helices. Φ-value analysis shows that the first helix becomes substantially consolidated as the intermediate is formed in a few hundred microseconds, as does the second to a lesser extent. A native-like structure then is formed in a few hundred milliseconds as the whole structure consolidates. Peptide fragments corresponding to sequences containing the first two helices separately and together as a helix-loop-helix motif have little helical structure under conditions that favor folding. The early stages of folding fit the nucleation-condensation model that was proposed for the smaller chymotrypsin inhibitor 2, which is a single module of structure and folds by two-state kinetics. The early stages of the multistate folding of the larger, multimodular, barnase have proved experimentally inaccessible. The folding pathway of barstar links those of CI2 and barnase to give a unified scheme for folding.

The rapid association of barnase and its intracellular inhibitor barstar has been analysed from the effects of mutagenesis and electrostatic screening. A basal association rate constant of 10⁵ M^-1 s^-1 is increased to over 5 x 10⁹ M^-1 s^-1 by electrostatic forces. The association between the oppositely charged proteins proceeds through the rate-determining formation of an early, weakly specific complex, which is dominated by long- range electrostatic interactions, followed by precise docking to form the high affinity complex. This mode of binding is likely to be used widely in nature to increase association rate constants between molecules and its principles may be used for protein design.