Interactome-wide prediction of short, disordered protein interaction motifs in humans

Many of the specific functions of intrinsically disordered protein segments are mediated by Short Linear Motifs (SLiMs) interacting with other proteins. Well known examples include SLiMs that interact with 14-3-3, PDZ, SH2, SH3, and WW domains but the true extent and diversity of SLiM-mediated interactions is largely unknown. Here, we attempt to expand our knowledge of human SLiMs by applying in silico SLiM prediction to the human interactome. Combining data from seven different interaction databases, we analysed approximately 6000 protein-centred and 1600 domain-centred human interaction datasets of 3+ unrelated proteins that interact with a common partner. Results were placed in context through comparison to randomised datasets of similar size and composition. The search returned thousands of evolutionarily conserved, intrinsically disordered occurrences of hundreds of significantly enriched recurring motifs, including many that have never been previously identified (). In addition to True Positive results for at least 25 different known SLiMs, a striking number of "off-target" proteins/domains also returned significantly enriched known motifs. Often, this was due to the non-independence of the datasets, with many proteins sharing interaction partners or contributing interactions to multiple domain datasets. The majority of these motif classes, however, were also found to be significantly enriched in one or more randomised datasets. This highlights the need for care when interpreting motif predictions of this nature but also raises the possibility that SLiM occurrences may be successfully identified independently of interaction data. Although not as compositionally biased as previous studies, patterns matching known SLiMs tended to cluster into a few large groups of similar sequence, while novel predictions tended to be more distinctive and less abundant. Whether this is due to ascertainment bias or a true functional composition bias of SLiMs is not clear and warrants further investigation.     

D-SLIMMER: domain-SLiM interaction motifs miner for sequence based protein-protein interaction data

Many biologically important protein-protein interactions (PPIs) have been found to be mediated by short linear motifs (SLiMs). These interactions are mediated by the binding of a protein domain, often with a nonlinear interaction interface, to a SLiM. We propose a method called D-SLIMMER to mine for SLiMs in PPI data on the basis of the interaction density between a nonlinear motif (i.e., a protein domain) in one protein and a SLiM in the other protein. Our results on a benchmark of 113 experimentally verified reference SLiMs showed that D-SLIMMER outperformed existing methods notably for discovering domain-SLiMs interaction motifs. To illustrate the significance of the SLiMs detected, we highlighted two SLiMs discovered from the PPI data by D-SLIMMER that are variants of the known ELM SLiM, as well as a literature-backed SLiM that is yet to be listed in the reference databases. We also presented a novel SLiM predicted by D-SLIMMER that was strongly supported by existing biological literatures. These examples showed that D-SLIMMER is able to find SLiMs that are biologically relevant.     

A comparative study of human TLR 7/8 stimulatory trimer compositions in influenza A viral genomes

Background

Variation in the genomes of single-stranded RNA viruses affects their infectivity and pathogenicity in two ways. First, viral genome sequence variations lead to changes in viral protein sequences and activities. Second, viral genome sequence variation produces diversity at the level of nucleotide composition and diversity in the interactions between viral RNAs and host toll-like receptors (TLRs). A viral genome-typing method based on this type of diversity has not yet been established.

Methodology/Principal Findings

In this study, we propose a novel genomic trait called the “TLR stimulatory trimer composition” (TSTC) and two quantitative indicators, Score S and Score N, named “TLR stimulatory scores” (TSS). Using the complete genome sequences of 10,994 influenza A viruses (IAV) and 251 influenza B viruses, we show that TSTC analysis reveals the diversity of Score S and Score N among the IAVs isolated from various hosts. In addition, we show that low values of Score S are correlated with high pathogenicity and pandemic potential in IAVs. Finally, we use Score S and Score N to construct a logistic regression model to recognize IAV strains that are highly pathogenic or have high pandemic potential.

Conclusions/Significance

Results from the TSTC analysis indicate that there are large differences between human and avian IAV genomes (except for segment 3), as illustrated by Score S. Moreover, segments 1, 2, 3 and 4 may be major determinants of the stimulatory activity exerted on human TLRs 7 and 8. We also find that a low Score S value is associated with high pathogenicity and pandemic potential in IAV. The π value from the TSS-derived logistic regression model is useful for recognizing emerging IAVs that have high pathogenicity and pandemic potential.     

Triple reassortment increases compatibility among viral ribonucleoprotein genes of contemporary avian and human influenza A viruses

Compatibility among the influenza A virus (IAV) ribonucleoprotein (RNP) genes affects viral replication efficiency and can limit the emergence of novel reassortants, including those with potential pandemic risks. In this study, we determined the polymerase activities of 2,451 RNP reassortants among three seasonal and eight enzootic IAVs by using a minigenome assay. Results showed that the 2009 H1N1 RNP are more compatible with the tested enzootic RNP than seasonal H3N2 RNP and that triple reassortment increased such compatibility. The RNP reassortants among 2009 H1N1, canine H3N8, and avian H4N6 IAVs had the highest polymerase activities. Residues in the RNA binding motifs and the contact regions among RNP proteins affected polymerase activities. Our data indicates that compatibility among seasonal and enzootic RNPs are selective, and enzoosis of multiple strains in the animal-human interface can facilitate emergence of an RNP with increased replication efficiency in mammals, including humans.     

The PDZ-binding motif of the avian NS1 protein affects transmission of the 2009 influenza A(H1N1) virus

By nature of their segmented RNA genome, influenza A viruses (IAVs) have the potential to generate variants through a reassortment process. The influenza nonstructural (NS) gene is critical for a virus to counteract the antiviral responses of the host. Therefore, a newly acquired NS segment potentially determines the replication efficiency of the reassortant virus in a range of different hosts. In addition, the C-terminal PDZ-binding motif (PBM) has been suggested as a pathogenic determinant of IAVs. To gauge the pandemic potential from human and avian IAV reassortment, we assessed the replication properties of NS-reassorted viruses in cultured cells and in the lungs of mice and determined their transmissibility in guinea pigs. Compared with the recombinant A/Korea/01/2009 virus (rK09; 2009 pandemic H1N1 strain), the rK09/VN:NS virus, in which the NS gene was adopted from the A/Vietnam/1203/2004 virus (a human isolate of the highly pathogenic avian influenza H5N1 virus strains), exhibited attenuated virulence and reduced transmissibility. However, the rK09/VN:NS-PBM virus, harboring the PBM in the C-terminus of the NS1 protein, recovered the attenuated virulence of the rK09/VN:NS virus. In a guinea pig model, the rK09/VN:NS-PBM virus showed even greater transmission efficiency than the rK/09 virus. These results suggest that the PBM in the NS1 protein may determine viral persistence in the human and avian IAV interface.     

Attributes of short linear motifs

Traditionally, protein-protein interactions were thought to be mediated by large, structured domains. However, it has become clear that the interactome comprises a wide range of binding interfaces with varying degrees of flexibility, ranging from rigid globular domains to disordered regions that natively lack structure. Enrichment for disorder in highly connected hub proteins and its correlation with organism complexity hint at the functional importance of disordered regions. Nevertheless, they have not yet been extensively characterised. Shifting the attention from globular domains to disordered regions of the proteome might bring us closer to elucidating the dense and complex connectivity of the interactome. An important class of disordered interfaces are the compact mono-partite, short linear motifs (SLiMs, or eukaryotic linear motifs (ELMs)). They are evolutionarily plastic and interact with relatively low affinity due to the limited number of residues that make direct contact with the binding partner. These features confer to SLiMs the ability to evolve convergently and mediate transient interactions, which is imperative to network evolution and to maintain robust cell signalling, respectively. The ability to discriminate biologically relevant SLiMs by means of different attributes will improve our understanding of the complexity of the interactome and aid development of bioinformatics tools for motif discovery. In this paper, the curated instances currently available in the Eukaryotic Linear Motif (ELM) database are analysed to provide a clear overview of the defining attributes of SLiMs. These analyses suggest that functional SLiMs have higher levels of conservation than their surrounding residues, frequently evolve convergently, preferentially occur in disordered regions and often form a secondary structure when bound to their interaction partner. These results advocate searching for small groupings of residues in disordered regions with higher relative conservation and a propensity to form the secondary structure. Finally, the most interesting conclusions are examined in regard to their functional consequences.     

Viral RNA-binding ability conferred by SUMOylation at PB1 K612 of influenza A virus is essential for viral pathogenesis and transmission

Posttranslational modifications, such as SUMOylation, play specific roles in the life cycle of invading pathogens. However, the effect of SUMOylation on the adaptation, pathogenesis, and transmission of influenza A virus (IAV) remains largely unknown. Here, we found that a conserved lysine residue at position 612 (K612) of the polymerase basic protein 1 (PB1) of IAV is a bona fide SUMOylation site. SUMOylation of PB1 at K612 had no effect on the stability or cellular localization of PB1, but was critical for viral ribonucleoprotein (vRNP) complex activity and virus replication in vitro. When tested in vivo, we found that the virulence of SUMOylation-defective PB1/K612R mutant IAVs was highly attenuated in mice. Moreover, the airborne transmission of a 2009 pandemic H1N1 PB1/K612R mutant virus was impaired in ferrets, resulting in reversion to wild-type PB1 K612. Mechanistically, SUMOylation at K612 was essential for PB1 to act as the enzymatic core of the viral polymerase by preserving its ability to bind viral RNA. Our study reveals an essential role for PB1 K612 SUMOylation in the pathogenesis and transmission of IAVs, which can be targeted for the design of anti-influenza therapies.     

SLiM-binding pockets: an attractive target for broad-spectrum antivirals

Short linear motif (SLiM)-mediated interactions offer a unique strategy for viral intervention due to their compact interfaces, ease of convergent evolution, and key functional roles. Consequently, many viruses extensively mimic host SLiMs to hijack or deregulate cellular pathways and the same motif-binding pocket is often targeted by numerous unrelated viruses. A toolkit of therapeutics targeting commonly mimicked SLiMs could provide prophylactic and therapeutic broad-spectrum antivirals and vastly improve our ability to treat ongoing and future viral outbreaks. In this opinion article, we discuss the therapeutic relevance of SLiMs, advocating their suitability as targets for broad-spectrum antiviral inhibitors.     

Oseltamivir-Resistant Influenza A (H1N1) Virus Strain with an H274Y Mutation in Neuraminidase Persists without Drug Pressure in Infected Mallards

Influenza A virus (IAV) has its natural reservoir in wild waterfowl, and emerging human IAVs often contain gene segments from avian viruses. The active drug metabolite of oseltamivir (oseltamivir carboxylate [OC]), stockpiled as Tamiflu for influenza pandemic preparedness, is not removed by conventional sewage treatment and has been detected in river water. There, it may exert evolutionary pressure on avian IAV in waterfowl, resulting in the development of resistant viral variants. A resistant avian IAV can circulate among wild birds only if resistance does not restrict viral fitness and if the resistant virus can persist without continuous drug pressure. In this in vivo mallard (Anas platyrhynchos) study, we tested whether an OC-resistant avian IAV (H1N1) strain with an H274Y mutation in the neuraminidase (NA-H274Y) could retain resistance while drug pressure was gradually removed. Successively infected mallards were exposed to decreasing levels of OC, and fecal samples were analyzed for the neuraminidase sequence and phenotypic resistance. No reversion to wild-type virus was observed during the experiment, which included 17 days of viral transmission among 10 ducks exposed to OC concentrations below resistance induction levels. We conclude that resistance in avian IAV that is induced by exposure of the natural host to OC can persist in the absence of the drug. Thus, there is a risk that human-pathogenic IAVs that evolve from IAVs circulating among wild birds may contain resistance mutations. An oseltamivir-resistant pandemic IAV would pose a substantial public health threat. Therefore, our observations underscore the need for prudent oseltamivir use, upgraded sewage treatment, and surveillance for resistant IAVs in wild birds.     

Hidden Multivalency in Phosphatase Recruitment by a Disordered AKAP Scaffold

Disordered scaffold proteins provide multivalent landing pads that, via a series of embedded Short Linear Motifs (SLiMs), bring together the components of a complex to orchestrate precise spatial and temporal regulation of cellular processes. One such protein is AKAP5 (previously AKAP79), which contains SLiMs that anchor PKA and Calcineurin, and recruit substrate (the TRPV1 receptor). Calcineurin is anchored to AKAP5 by a well-characterised PxIxIT SLiM. Here we show, using a combination of biochemical and biophysical approaches, that the Calcineurin PxIxIT-binding groove also recognises several hitherto unknown lower-affinity SLiMs in addition to the PxIxIT motif. We demonstrate that the assembly is in reality a complex system with conserved SLiMs spanning a wide affinity range. The capture is analogous to that seen for many DNA-binding proteins that have a weak non-specific affinity for DNA outside the canonical binding site, but different in that it involves (i) two proteins, and (ii) hydrophobic rather than electrostatic interactions. It is also compatible with the requirement for both stable anchoring of the enzyme and responsive downstream signalling. We conclude that the AKAP5 C-terminus is enriched in lower-affinity/mini-SLiMs that, together with the canonical SLiM, maintain a structurally disordered but tightly regulated signalosome.     

Non-Avian Animal Reservoirs Present a Source of Influenza A PB1-F2 Proteins with Novel Virulence-Enhancing Markers

PB1-F2 protein, expressed from an alternative reading frame of most influenza A virus (IAV) PB1 segments, may possess specific residues associated with enhanced inflammation (L62, R75, R79, and L82) and cytotoxicity (I68, L69, and V70). These residues were shown to increase the pathogenicity of primary viral and secondary bacterial infections in a mouse model. In contrast to human seasonal influenza strains, virulence-associated residues are present in PB1-F2 proteins from pandemic H1N1 1918, H2N2 1957, and H3N2 1968, and highly pathogenic H5N1 strains, suggesting their contribution to viruses'' pathogenic phenotypes. Non-human influenza strains may act as donors of virulent PB1-F2 proteins. Previously, avian influenza strains were identified as a potential source of inflammatory, but not cytotoxic, PB1-F2 residues. Here, we analyze the frequency of virulence-associated residues in PB1-F2 sequences from IAVs circulating in mammalian species in close contact with humans: pigs, horses, and dogs. All four inflammatory residues were found in PB1-F2 proteins from these viruses. Among cytotoxic residues, I68 was the most common and was especially prevalent in equine and canine IAVs. Historically, PB1-F2 from equine (about 75%) and canine (about 20%) IAVs were most likely to have combinations of the highest numbers of residues associated with inflammation and cytotoxicity, compared to about 7% of swine IAVs. Our analyses show that, in addition to birds, pigs, horses, and dogs are potentially important sources of pathogenic PB1-F2 variants. There is a need for surveillance of IAVs with genetic markers of virulence that may be emerging from these reservoirs in order to improve pandemic preparedness and response.     

Identification of a highly conserved H1 subtype-specific epitope with diagnostic potential in the hemagglutinin protein of influenza A virus

Subtype specificity of influenza A virus (IAV) is determined by its two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). For HA, 16 distinct subtypes (H1–H16) exist, while nine exist for NA. The epidemic strains of H1N1 IAV change frequently and cause annual seasonal epidemics as well as occasional pandemics, such as the notorious 1918 influenza pandemic. The recent introduction of pandemic A/H1N1 IAV (H1N1pdm virus) into humans re-emphasizes the public health concern about H1N1 IAV. Several studies have identified conserved epitopes within specific HA subtypes that can be used for diagnostics. However, immune specific epitopes in H1N1 IAV have not been completely assessed. In this study, linear epitopes on the H1N1pdm viral HA protein were identified by peptide scanning using libraries of overlapping peptides against convalescent sera from H1N1pdm patients. One epitope, P5 (aa 58–72) was found to be immunodominant in patients and to evoke high titer antibodies in mice. Multiple sequence alignments and in silico coverage analysis showed that this epitope is highly conserved in influenza H1 HA [with a coverage of 91.6% (9,860/10,767)] and almost completely absent in other subtypes [with a coverage of 3.3% (792/23,895)]. This previously unidentified linear epitope is located outside the five well-recognized antigenic sites in HA. A peptide ELISA method based on this epitope was developed and showed high correlation (χ² = 51.81, P<0.01, Pearson correlation coefficient R = 0.741) with a hemagglutination inhibition test. The highly conserved H1 subtype-specific immunodominant epitope may form the basis for developing novel assays for sero-diagnosis and active surveillance against H1N1 IAVs.     

Timing Is Everything: Coordinated Control of Host Shutoff by Influenza A Virus NS1 and PA-X Proteins

Like all viruses, influenza viruses (IAVs) use host translation machinery to decode viral mRNAs. IAVs ensure efficient translation of viral mRNAs through host shutoff, a process whereby viral proteins limit the accumulation of host proteins through subversion of their biogenesis. Despite its small genome, the virus deploys multiple host shutoff mechanisms at different stages of infection, thereby ensuring successful replication while limiting the communication of host antiviral responses. In this Gem, we review recent data on IAV host shutoff proteins, frame the outstanding questions in the field, and propose a temporally coordinated model of IAV host shutoff.     

Porcine pulmonary collectins show distinct interactions with influenza A viruses: role of the N-linked oligosaccharides in the carbohydrate recognition domain

Influenza A virus (IAV) infections are a major cause of respiratory disease of humans and animals. Pigs can serve as important intermediate hosts for transmission of avian IAV strains to humans, and for the generation of reassortant strains; this may result in the appearance of new pandemic IAV strains in humans. We have studied the role of the porcine lung collectins surfactant proteins D and A (pSP-D and pSP-A), two important components of the innate immune response against IAV. Hemagglutination inhibition assays revealed that both pSP-D and pSP-A display substantially greater inhibitory activity against IAV strains isolated from human, swine, and horse, than lung collectins from other animal species. The more potent activity of pSP-D results from interactions mediated by the asparagine-linked oligosaccharide located in the carbohydrate recognition domain of pSP-D, which is absent in SP-Ds from other species characterized to date. Presence of this sialylated oligosaccharide moiety enhances the anti-influenza activity of pSP-D, as demonstrated by assays of viral aggregation, inhibition of infectivity, and neutrophil response to IAV. The greater hemagglutination inhibitory activity of pSP-A is due to porcine-specific structural features of the conserved asparagine-linked oligosaccharide in the carbohydrate recognition domain of SP-A. A more efficient lung collectin-mediated immune response against IAV in pigs may play a role in providing conditions by which pigs can act as "mixing vessel" hosts that can lead to the production of reassortant, pandemic strains of IAV.     

Universal and strain specific structure features of segment 8 genomic RNA of influenza A virus—application of 4-thiouridine photocrosslinking

RNA structure in the influenza A virus (IAV) has been the focus of several studies that have shown connections between conserved secondary structure motifs and their biological function in the virus replication cycle. Questions have arisen on how to best recognize and understand the pandemic properties of IAV strains from an RNA perspective, but determination of the RNA secondary structure has been challenging. Herein, we used chemical mapping to determine the secondary structure of segment 8 viral RNA (vRNA) of the pandemic A/California/04/2009 (H1N1) strain of IAV. Additionally, this long, naturally occurring RNA served as a model to evaluate RNA mapping with 4-thiouridine (4sU) crosslinking. We explored 4-thiouridine as a probe of nucleotides in close proximity, through its incorporation into newly transcribed RNA and subsequent photoactivation. RNA secondary structural features both universal to type A strains and unique to the A/California/04/2009 (H1N1) strain were recognized. 4sU mapping confirmed and facilitated RNA structure prediction, according to several rules: 4sU photocross-linking forms efficiently in the double-stranded region of RNA with some flexibility, in the ends of helices, and across bulges and loops when their structural mobility is permitted. This method highlighted three-dimensional properties of segment 8 vRNA secondary structure motifs and allowed to propose several long-range three-dimensional interactions. 4sU mapping combined with chemical mapping and bioinformatic analysis could be used to enhance the RNA structure determination as well as recognition of target regions for antisense strategies or viral RNA detection.     

MotifAnalyzer‐PDZ: A computational program to investigate the evolution of PDZ‐binding target specificity

Recognition of short linear motifs (SLiMs) or peptides by proteins is an important component of many cellular processes. However, due to limited and degenerate binding motifs, prediction of cellular targets is challenging. In addition, many of these interactions are transient and of relatively low affinity. Here, we focus on one of the largest families of SLiM‐binding domains in the human proteome, the PDZ domain. These domains bind the extreme C‐terminus of target proteins, and are involved in many signaling and trafficking pathways. To predict endogenous targets of PDZ domains, we developed MotifAnalyzer‐PDZ, a program that filters and compares all motif‐satisfying sequences in any publicly available proteome. This approach enables us to determine possible PDZ binding targets in humans and other organisms. Using this program, we predicted and biochemically tested novel human PDZ targets by looking for strong sequence conservation in evolution. We also identified three C‐terminal sequences in choanoflagellates that bind a choanoflagellate PDZ domain, the Monsiga brevicollis SHANK1 PDZ domain (mbSHANK1), with endogenously‐relevant affinities, despite a lack of conservation with the targets of a homologous human PDZ domain, SHANK1. All three are predicted to be signaling proteins, with strong sequence homology to cytosolic and receptor tyrosine kinases. Finally, we analyzed and compared the positional amino acid enrichments in PDZ motif‐satisfying sequences from over a dozen organisms. Overall, MotifAnalyzer‐PDZ is a versatile program to investigate potential PDZ interactions. This proof‐of‐concept work is poised to enable similar types of analyses for other SLiM‐binding domains (e.g., MotifAnalyzer‐Kinase). MotifAnalyzer‐PDZ is available at http://motifAnalyzerPDZ.cs.wwu.edu .     

Cleavage of Hemagglutinin-Bearing Lentiviral Pseudotypes and Their Use in the Study of Influenza Virus Persistence

Influenza A viruses (IAVs) are a major cause of infectious respiratory human diseases and their transmission is dependent upon the environment. However, the role of environmental factors on IAV survival outside the host still raises many questions. In this study, we used lentiviral pseudotypes to study the influence of the hemagglutinin protein in IAV survival. High-titered and cleaved influenza-based lentiviral pseudoparticles, through the use of a combination of two proteases (HAT and TMPRSS2) were produced. Pseudoparticles bearing hemagglutinin proteins derived from different H1N1, H3N2 and H5N1 IAV strains were subjected to various environmental parameters over time and tested for viability through single-cycle infectivity assays. We showed that pseudotypes with different HAs have different persistence profiles in water as previously shown with IAVs. Our results also showed that pseudotypes derived from H1N1 pandemic virus survived longer than those derived from seasonal H1N1 virus from 1999, at high temperature and salinity, as previously shown with their viral counterparts. Similarly, increasing temperature and salinity had a negative effect on the survival of the H3N2 and H5N1 pseudotypes. These results showed that pseudotypes with the same lentiviral core, but which differ in their surface glycoproteins, survived differently outside the host, suggesting a role for the HA in virus stability.     

Avian influenza A virus susceptibility,infection, transmission,and antibody kinetics in European starlings

Avian influenza A viruses (IAVs) pose risks to public, agricultural, and wildlife health. Bridge hosts are spillover hosts that share habitat with both maintenance hosts (e.g., mallards) and target hosts (e.g., poultry). We conducted a comprehensive assessment of European starlings (Sturnus vulgaris), a common visitor to both urban and agricultural environments, to assess whether this species might act as a potential maintenance or bridge host for IAVs. First, we experimentally inoculated starlings with a wild bird IAV to investigate susceptibility and replication kinetics. Next, we evaluated whether IAV might spill over to starlings from sharing resources with a widespread IAV reservoir host. We accomplished this using a specially designed transmission cage to simulate natural environmental transmission by exposing starlings to water shared with IAV-infected mallards (Anas platyrhynchos). We then conducted a contact study to assess intraspecies transmission between starlings. In the initial experimental infection study, all inoculated starlings shed viral RNA and seroconverted. All starlings in the transmission study became infected and shed RNA at similar levels. All but one of these birds seroconverted, but detectable antibodies were relatively transient, falling to negative levels in a majority of birds by 59 days post contact. None of the contact starlings in the intraspecies transmission experiment became infected. In summary, we demonstrated that starlings may have the potential to act as IAV bridge hosts if they share water with IAV-infected waterfowl. However, starlings are unlikely to act as maintenance hosts due to limited, if any, intraspecies transmission. In addition, starlings have a relatively brief antibody response which should be considered when interpreting serology from field samples. Further study is needed to evaluate the potential for transmission from starlings to poultry, a possibility enhanced by starling’s behavioral trait of forming very large flocks which can descend on poultry facilities when natural resources are scarce.