Publications

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.

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.

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.

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.

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.

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.

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.

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. [Display omitted] •The effective reproduction number of BA.2 is 1.4-fold higher than that of BA.1•BA.2 is resistant to BA.1-induced humoral immunity•The BA.2 spike is more fusogenic than BA.1 spike•BA.2 spike-bearing virus is more pathogenic than BA.1 spike-bearing virus Yamasoba and G2P-Japan Consortium et al. elucidate the characteristics of SARS-CoV-2 Omicron variant BA.2—transmissibility, immune resistance, virological property, and pathogenicity. The effective population number of BA.2 is higher than that of BA.1, and the antigenicity of BA.2 is different from that of BA.1. The BA.2 spike is more fusogenic than the BA.1 spike, and notably, the BA.2 spike-bearing virus is more pathogenic than the BA.1 spike-bearing virus. This multiscale investigation suggests the potential risk of BA.2 to global health.

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.

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 (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[beta] production. Expression of the first three of these SARS-CoV-2 genes specifically blocked MAVS-induced IFN[beta]-promoter activity, whereas all six genes induced a collapse in IFN[beta] mRNA levels, corresponding with suppressed IFN[beta] protein secretion. Five of these six genes furthermore suppressed MAVS-induced activation of IFN[lambda]s, however with no effect on IFN[alpha] or IFN[gamma] 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 ([DELTA]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.

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 protein–protein interactions. While the methodological changes are incremental, it results in a qualitative enhancement in the applicability of yeast display for many applications.

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-λ.

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.

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.

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.

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 Administration–approved 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.

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 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.

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/.

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.

IFN is a common therapeutic option to treat multiple sclerosis (MS). It is unique amongst 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 YNS8) that is built on the backbone of a low affinity IFN, but modified to exhibit higher receptor affinity than IFN. Here, YNS8 was fused with a 600-residue hydrophilic, unstructured N-terminal polypeptide chain, comprising Proline, Alanine and Serine (PAS), in order to prolong its plasma half-life via “PASylation”. PAS-IFN8 exhibited a 10-fold increased half-life in both pharmacodynamic and pharmacokinetic assays in a transgenic mouse model harboring the human receptors [1], notably without any detectable loss in biological potency or bioavailability. This long-life superagonist conferred significantly improved protection from MOG35–55 peptide-induced experimental autoimmune encephalomyelitis (EAE) compared to IFN, despite being injected with a 4-fold decrease in frequency and at an overall 16-fold lower dosage. These data were corroborated by FACS and immunohistochemistry, showing a decrease of myeloid lineage cells in the central nervous system (CNS) for PAS-YNS8 treated animals. Importantly, PAS-IFN8 did not induce antibodies upon repeated administration, and its biological efficacy remained unchanged after 21 days of treatment. We next extracted RNA from spleen-derived CD4 + cells and performed high throughput qPCR for 200 IFN-response genes enriched for immune function. A striking strong correlation of gene expression and EAE clinical response was observed for two genes - CD274 (PD-L1) and CXCR3. An IFN-responsive role for PD-L1 in EAE clinical response has been recently reported elsewhere [2]. This preclinical study thus supports that we have generated a novel, potent and pharmacologically long-acting IFN-variant with potential to treat multiple sclerosis with greater efficacy than current IFN therapies available for use today. It is also provides us with a means to unravel the central IFN-response genes that are providing clinical benefit in EAE and potentially for MS.

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.

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

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.

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.

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.

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.

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.

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.

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.

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.

Keywords: Biochemistry & Molecular Biology; Cell Biology; Immunology

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.

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.

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 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.

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.

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.

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.

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?

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.

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.

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.

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.

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.

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-β.

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.

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.

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.

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.

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.

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.

[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.

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.

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

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.

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.

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.

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.

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.