Normalization of nomenclature for peptide motifs as ligands of modular protein domains

Unveiling multiple copies of MlaC highlights its multifaceted nature

The maintenance of the lipid asymmetry (Mla) system plays a critical role in facilitating the transport of phospholipids between the inner and outer membranes of the Gram-negative bacteria. In E. coli, the system consists of six proteins: MlaA-OmpF/C complex (outer membrane), MlaC (periplasm), and MlaFEDB complex (inner membrane). Despite extensive research on the core proteins (MlaFED) of the Mla system, the occurrence of Mla components like MlaA, MlaB, and MlaC in diderm remains uncertain. Therefore, this gap presents a significant opportunity for further investigation, particularly regarding MlaC, which serves as the sole mobile component of the Mla system. This has led to the identification of multiple copies of MlaC in 63 distinct genera of Proteobacteria and related phyla. Interestingly, amongst these genera, the genetic arrangements of the mla operon were observed to be varying and, thus, were further categorized into four distinct groups. The variations among the genetic organization of the mla operons suggest their evolution through various processes, such as duplications, losses, rearrangements, and fusions. Further, the results of this study highlight the MlaC's substrate promiscuity, illuminating new avenues for the Mla system.

Exploring the druggability of the UEV domain of human TSG101 in search for broad-spectrum antivirals

The ubiquitin E2 variant domain of TSG101 (TSG101-UEV) plays a pivotal role in protein sorting and virus budding by recognizing PTAP motifs within ubiquitinated proteins. Disruption of TSG101-UEV/PTAP interactions has emerged as a promising strategy for the development of host-oriented broad-spectrum antivirals with low susceptibility to resistance. TSG101 is a challenging target characterized by an extended and flat binding interface, low affinity for PTAP ligands, and complex binding energetics. Here, we assess the druggability of the TSG101-UEV/PTAP binding interface by searching for drug-like inhibitors and evaluating their ability to block PTAP recognition, impair budding, and inhibit viral proliferation. A discovery workflow was established by combining in vitro miniaturized HTS assays and a set of cell-based activity assays including high-content bimolecular complementation, virus-like particle release measurement, and antiviral testing in live virus infection. This approach has allowed us to identify a set of chemically diverse molecules that block TSG101-UEV/PTAP binding with IC50s in the low μM range and are able to disrupt the interaction between full-length TSG101 and viral proteins in human cells and inhibit viral replication. State-of-the-art molecular docking studies reveal that the active compounds exploit binding hotspots at the PTAP binding site, unlocking the full binding potential of the TSG101-UEV binding pockets. These inhibitors represent promising hits for the development of novel broad-spectrum antivirals through targeted optimization and are also valuable tools for investigating the involvement of ESCRT in the proliferation of different virus families and study the secondary effects induced by the disruption of ESCRT/virus interactions.

CompariPSSM: a PSSM-PSSM comparison tool for motif-binding determinant analysis

MOTIVATION

Short linear motifs (SLiMs) are compact functional modules that mediate low-affinity protein-protein interactions. SLiMs direct the function of many dynamic signalling and regulatory complexes playing a central role in most biological processes of the cell. Motif-binding determinants describe the contribution of each residue in a motif-containing peptide to the affinity and specificity of binding to the motif-binding partner. Motif-binding determinants are generally defined as a motif consensus pattern or a position-specific scoring matrix (PSSM) encoding quantitative preferences. Motif-binding determinant comparison is an important motif analysis task and can be applied to motif annotation, classification, clustering, discovery and benchmarking. Currently, binding determinant comparison is generally performed by analysing consensus similarity; however, this ignores important quantitative information in both the consensus and non-consensus positions.

RESULTS

We have created a new tool, CompariPSSM, that quantifies the similarity between motif-binding determinants using sliding window PSSM-PSSM comparison and scores PSSM similarity using a randomisation-based probabilistic framework. The tool has been benchmarked on curated data from the eukaryotic linear motif database and experimental data from proteomic peptidephage display. CompariPSSM can be used for peptide classification to validate motif classes, peptide clustering to group functionally related conserved disordered regions, and benchmarking experimental motif discovery methods.

AVAILABILITY AND IMPLEMENTATION

CompariPSSM is available at https://slim.icr.ac.uk/projects/comparipssm.

Reversal of tyrosine-linked ADP-ribosylation by ARH3 and PARG

ADP-ribosylation is an ancient posttranslational modification with exceptional versatility in terms of breadth of modification targets including at least seven different amino acid side chains, various moieties on nucleic acids, and a variety of small chemical compounds. The spatiotemporal signaling dynamic of the different modification variations is tightly regulated and depends on the writers, erases, and readers of each type. Among these, tyrosine ADP-ribosylation (Tyr-ADPr) has been consistently detected as a novel modification type, but systematic analysis of its potential physiological role, modification establishment, and reversal are still lacking. Here we present a re-analysis of recent ADP-ribosylome data and show that Tyr-ADPr sites are conserved and enriched among ribosome biogenesis and mRNA processing proteins and that these sites are affected by the status of the (ADP-ribosyl)hydrolase ARH3. To facilitate the study of Tyr-ADPr, we establish methodologies for the synthesis of well-defined Tyr-ADPr peptides and with these could show that Tyr-ADPr is reversed both by ARH3 and PARG enzymes. Together, our work lays the foundation for the future exploration of the Tyr-ADPr.

Evolutionary and molecular basis of ADP-ribosylation reversal by zinc-dependent macrodomains

Dynamic ADP-ribosylation signaling is a crucial pathway that controls fundamental cellular processes, in particular, the response to cellular stresses such as DNA damage, reactive oxygen species, and infection. In some pathogenic microbes, the response to oxidative stress is controlled by a SirTM/zinc-containing macrodomain (Zn-Macro) pair responsible for establishment and removal of the modification, respectively. Targeting this defence mechanism against the host's innate immune response may lead to novel approaches to support the fight against emerging antimicrobial resistance. Earlier studies suggested that Zn-Macros play a key role in the activation of this defence. Therefore, we used phylogenetic, biochemical, and structural approaches to elucidate the functional properties of these enzymes. Using the substrate mimetic asparagine-ADP-ribose as well as the ADP-ribose product, we characterize the catalytic role of the zinc ion in the removal of the ADP-ribosyl modification. Furthermore, we determined structural properties that contribute to substrate selectivity within the different Zn-Macro branches. Together, our data not only give new insights into the Zn-Macro family but also highlight their distinct features that may be exploited for the development of future therapies.

Dishevelled2 activates WGEF via its interaction with a unique internal peptide motif of the GEF

The Wnt-planar cell polarity (Wnt-PCP) pathway is crucial in establishing cell polarity during development and tissue homoeostasis. This pathway is found to be dysregulated in many pathological conditions, including cancer and autoimmune disorders. The central event in Wnt-PCP pathway is the activation of Weak-similarity guanine nucleotide exchange factor (WGEF) by the adapter protein Dishevelled (Dvl). The PDZ domain of Dishevelled2 (Dvl2PDZ) binds and activates WGEF by releasing it from its autoinhibitory state. However, the actual Dvl2PDZ binding site of WGEF and the consequent activation mechanism of the GEF have remained elusive. Using biochemical and molecular dynamics studies, we show that a unique "internal-PDZ binding motif" (IPM) of WGEF mediates the WGEF-Dvl2PDZ interaction to activate the GEF. The residues at P2, P0, P-2 and P-3 positions of IPM play an important role in stabilizing the WGEFpep-Dvl2PDZ interaction. Furthermore, MD simulations of modelled Dvl2PDZ-WGEFIPM peptide complexes suggest that WGEF-Dvl2PDZ interaction may differ from the reported Dvl2PDZ-IPM interactions. Additionally, the apo structure of human Dvl2PDZ shows conformational dynamics different from its IPM peptide bound state, suggesting an induced fit mechanism for the Dvl2PDZ-peptide interaction. The current study provides a model for Dvl2 induced activation of WGEF.

PTEN-regulated PI3K-p110 and AKT isoform plasticity controls metastatic prostate cancer progression

PTEN loss, one of the most frequent mutations in prostate cancer (PC), is presumed to drive disease progression through AKT activation. However, two transgenic PC models with Akt activation plus Rb loss exhibited different metastatic development: Pten/RbPE:-/- mice produced systemic metastatic adenocarcinomas with high AKT2 activation, whereas RbPE:-/- mice deficient for the Src-scaffolding protein, Akap12, induced high-grade prostatic intraepithelial neoplasias and indolent lymph node dissemination, correlating with upregulated phosphotyrosyl PI3K-p85α. Using PC cells isogenic for PTEN, we show that PTEN-deficiency correlated with dependence on both p110β and AKT2 for in vitro and in vivo parameters of metastatic growth or motility, and with downregulation of SMAD4, a known PC metastasis suppressor. In contrast, PTEN expression, which dampened these oncogenic behaviors, correlated with greater dependence on p110α plus AKT1. Our data suggest that metastatic PC aggressiveness is controlled by specific PI3K/AKT isoform combinations influenced by divergent Src activation or PTEN-loss pathways.

Bicyclic Peptide Library Screening for the Identification of Gαi Protein Modulators

Noncanonical G protein activation and inactivation, particularly for the Gαi/s protein subfamilies, have long been a focus of chemical research. Combinatorial libraries were already effectively applied to identify modulators of the guanine-nucleotide exchange, as can be exemplified with peptides such as KB-752 and GPM-1c/d, the so-called guanine-nucleotide exchange modulators. In this study, we identified novel bicyclic peptides from a combinatorial library screening that show prominent properties as molecular switch-on/off modulators of Gαi signaling. Among the series of hits, the exceptional paradigm of GPM-3, a protein and state-specific bicyclic peptide, is the first chemically identified GAP (GTPase-activating protein) modulator with a high binding affinity for Gαi protein. Computational analyses identified and assessed the structure of the bicyclic peptides, novel ligand-protein interaction sites, and their subsequent impact on the nucleotide binding site. This approach can therefore lead the way for the development of efficient chemical biological probes targeting Gαi protein modulation within a cellular context.

Canonical structural-binding modes in the calmodulin-target protein complexes

Intracellular calcium sensor protein calmodulin (CaM) belongs to the large EF-hand protein superfamily. CaM shows a unique and not fully understood ability to bind to multiple targets, allows them to participate in a variety of regulatory processes. The protein has two approximately symmetrical globular domains (the N- and C-lobes). Analysis of the CaM-binding sites of target proteins showed that they have two hydrophobic 'anchor' amino acids separated by 10 to 17 residues. Consequently, several CaM-binding motifs: {1-10}, {1-11}, {1-13}, {1-14}, {1-16}, {1-17}, differing by the distance between the two anchor residues along the amino acid sequence, have been identified. Despite extensive structural information on the role of target-protein amino acid residues in the formation of complexes with CaM, much less is known about the role of amino acids from CaM contributing to these interactions. In this work, a quantitative analysis of the contact surfaces of CaM and target proteins has been carried out for 35 representative three-dimensional structures. It has been shown that, in addition to the two hydrophobic terminal residues of the target fragment, the interaction also involves residues that are 4 residues earlier in the sequence (binding mode {1-5}). It has also been found that the N- and C-lobes of CaM bind the {1-5} motif located at the ends of the target in a structurally identical manner. Methionine residues at positions 51 (corresponding to 124 in the C-lobe), 71 (144), and 72 (145) of the CaM amino acid sequence are key hydrophobic residues for this interaction. They are located at the N- and C-boundaries of the even EF-hand motifs. The hydrophobic core of CaM ('Ф-quatrefoil') consists of 10 amino acids in the N-lobe (and in the C-lobe): Phe16 (Phe89), Phe19 (Phe92), Ile27 (Ile100), Thr29 (Ala102), Leu32 (Leu105), Ile52 (Ile125), Val55 (Ala128), Ile63 (Val136), Phe65 (Tyr138), and Phe68 (Phe141) and do not intersect with the target-binding methionine residues. CaM belongs to the 'dynamic' group of EF-hand proteins, in which calcium and protein ligand binding causes only global conformational changes but does not alter the conservative 'black' and 'grey' clusters described in our earlier works (PLoS One. 2014; 9(10):e109287). The membership of CaM in the 'dynamic' group is determined by the triggering and protective methionine layer: Met51 (Met124), Met71 (Met144) and Met72 (Met145). HIGHLIGHTSInterchain interactions in the unique 35 CaM complex structures were analyzed.Methionine amino acids of the N- and C-lobes of CaM form triggering and protective layers.Interactions of the target terminal residues with these methionine layers are structurally identical.CaM belonging to the 'dynamic' group is determined by the triggering and protective methionine layer.Communicated by Ramaswamy H. Sarma.

PTEN regulated PI3K-p110 and AKT isoform plasticity controls metastatic prostate cancer progression

PTEN loss, one of the most frequent mutations in prostate cancer (PC), is presumed to drive disease progression through AKT activation. However, two transgenic PC models with Akt activation plus Rb loss exhibited different metastasis development: Pten/RbPE:-/- mice produced systemic metastatic adenocarcinomas with high AKT2 activation, whereas RbPE:-/- mice deficient for the Src-scaffolding protein, Akap12, induced high-grade prostatic intraepithelial neoplasias and indolent lymph node disseminations, correlating with upregulated phosphotyrosyl PI3K-p85α. Using PC cells isogenic for PTEN, we show that PTEN-deficiency correlated with dependence on both p110β and AKT2 for in vitro and in vivo parameters of metastatic growth or motility, and with downregulation of SMAD4, a known PC metastasis suppressor. In contrast, PTEN expression, which dampened these oncogenic behaviors, correlated with greater dependence on p110α plus AKT1. Our data suggest that metastatic PC aggressiveness is controlled by specific PI3K/AKT isoform combinations influenced by divergent Src activation or PTEN-loss pathways.

Substrate-binding glycine residues are major determinants for hydrolase and ligase activity of plant legumains

Peptide asparaginyl ligases (PALs) are useful tools for precision modifications of proteins and live-cell surfaces by ligating peptides after Asn/Asp (Asx). They share high sequence and structural similarity to plant legumains that are generally known as asparaginyl endopeptidases (AEPs), thus making it challenging to identify PALs from AEPs. In this study, we investigate 875 plant species from algae to seed plants with available sequence data in public databases to identify new PALs. We conducted evolutionary trace analysis on 1500 plant legumains, including eight known PALs, to identify key residues that could differentiate ligases and proteases, followed by recombinant expression and functional validation of 16 novel legumains. Previously, we showed that the substrate-binding sequences flanking the catalytic site can strongly influence the enzymatic direction of a legumain and which we named as ligase-activity determinants (LADs). Here, we show that two conserved substrate-binding Gly residues of LADs are critical, but negative determinants for ligase activity. Our results suggest that specific glycine residues are molecular determinants to identify PALs and AEPs as two different legumain subfamilies, accounting for c. 1% and 88%, respectively.

Reconfigurable Peptide Analogs of Apolipoprotein A-I Reveal Tunable Features of Nanodisc Assembly

Nanodisc (ND)-forming membrane scaffold proteins or peptides developed from apolipoprotein A-I (apoA-I) have led to considerable promise in structural biology and therapeutic applications. However, the rationale and regularity characteristics in peptide sequence design remain inconclusive. Here, we proposed a consensus-based normalization approach through the reversed engineering of apoA-IΔ1-45 to design reconfigurable apoA-I peptide analogs (APAs) for tunable ND assembly. We present extensive morphological validations and computational simulation analyses on divergent APA-NDs that are generated by our method. Fifteen divergent APAs were generated accordingly to study the assembly machinery of NDs. We show that APA designs exhibit multifactorial influence in terms of varying APA tandem repeats, sequence composition, and lipid-to-APA ratio to form tunable diameters of NDs. There is a strong positive correlation between DMPC-to-APA ratios and ND diameters. Longer APA with more tandem repeats tends to yield higher particle size homogeneity. Our results also suggest proline is a dispensable residue for the APA-ND formation. Interestingly, proline-rich substitution not only provides an inward-bending effect in forming smaller NDs but also induces the cumulative chain flexibility that enables larger ND formation at higher lipid ratios. Additionally, proline-tryptophan residues in APAs play a dominant role in forming larger NDs. Molecular simulation shows that enriched basic and acidic residues in APAs evoke abundant hydrogen-bond and salt bridge networks to reinforce the structural stability of APA-NDs. Together, our findings provide a rational basis for understanding APA design. The proposed model could be extended to other apolipoproteins for desired ND engineering.

Design and synthesis of a tetracyclic tripeptide mimetic frozen in a polyproline type II (PP2) helix conformation

A synthesis of the new tetracyclic scaffold ProM-19, which represents a XPP tripeptide unit frozen in a PPII helix conformation, was developed. As a key building block, N-Boc-protected ethyl (1S,3S,4R)-2-azabicyclo[2.2.1]hept-5-ene-2-carboxylate was prepared through a diastereoselective aza-Diels-Alder reaction and subsequent hydrogenolytic removal of the chiral N-1-phenylethyl substituent under temporary protection of the double bond through dihydroxylation and reconstitution by Corey-Winter olefination. The target compound Boc-[ProM-19]-OMe was then prepared via subsequent peptide coupling and Ru-catalyzed ring-closing metathesis steps employing (S)-N-Boc-allylgylcine and cis-5-vinyl-proline methyl ester as additional building blocks. In addition, Ac-[2-Cl-Phe]-[Pro]-[ProM-19]-OMe was prepared by solution phase peptide synthesis as a potential ligand for the ena-VASP EVH1 domain.

Kiaa1024L/Minar2 is essential for hearing by regulating cholesterol distribution in hair bundles

Unbiased genetic screens implicated a number of uncharacterized genes in hearing loss, suggesting some biological processes required for auditory function remain unexplored. Loss of Kiaa1024L/Minar2, a previously understudied gene, caused deafness in mice, but how it functioned in the hearing was unclear. Here, we show that disruption of kiaa1024L/minar2 causes hearing loss in the zebrafish. Defects in mechanotransduction, longer and thinner hair bundles, and enlarged apical lysosomes in hair cells are observed in the kiaa1024L/minar2 mutant. In cultured cells, Kiaa1024L/Minar2 is mainly localized to lysosomes, and its overexpression recruits cholesterol and increases cholesterol labeling. Strikingly, cholesterol is highly enriched in the hair bundle membrane, and loss of kiaa1024L/minar2 reduces cholesterol localization to the hair bundles. Lowering cholesterol levels aggravates, while increasing cholesterol levels rescues the hair cell defects in the kiaa1024L/minar2 mutant. Therefore, cholesterol plays an essential role in hair bundles, and Kiaa1024L/Minar2 regulates cholesterol distribution and homeostasis to ensure normal hearing.

Phage display identification of nanomolar ligands for human NEDD4-WW3: Energetic and dynamic implications for the development of broad-spectrum antivirals

The recognition of PPxY viral Late domains by the third WW domain of the human HECT-E3 ubiquitin ligase NEDD4 (NEDD4-WW3) is essential for the budding of many viruses. Blocking these interactions is a promising strategy to develop broad-spectrum antivirals. As all WW domains, NEDD4-WW3 is a challenging therapeutic target due to the low binding affinity of its natural interactions, its high conformational plasticity, and its complex thermodynamic behavior. In this work, we set out to investigate whether high affinity can be achieved for monovalent ligands binding to the isolated NEDD4-WW3 domain. We show that a competitive phage-display set-up allows for the identification of high-affinity peptides showing inhibitory activity of viral budding. A detailed biophysical study combining calorimetry, nuclear magnetic resonance, and molecular dynamic simulations reveals that the improvement in binding affinity does not arise from the establishment of new interactions with the domain, but is associated to conformational restrictions imposed by a novel C-terminal -LFP motif in the ligand, unprecedented in the PPxY interactome. These results, which highlight the complexity of WW domain interactions, provide valuable insight into the key elements for high binding affinity, of interest to guide virtual screening campaigns for the identification of novel therapeutics targeting NEDD4-WW3 interactions.

Targeting Gαi/s Proteins with Peptidyl Nucleotide Exchange Modulators

Chemical probes that specifically modulate the activity of heterotrimeric G proteins provide excellent tools for investigating G protein-mediated cell signaling. Herein, we report a family of selective peptidyl Gαi/s modulators derived from peptide library screening and optimization. Conjugation to a cell-penetrating peptide rendered the peptides cell-permeable and biologically active in cell-based assays. The peptides exhibit potent guanine-nucleotide exchange modulator-like activity toward Gαi and Gαs. Molecular docking and dynamic simulations revealed the molecular basis of the protein-ligand interactions and their effects on GDP binding. This study demonstrates the feasibility of developing direct Gαi/s modulators and provides a novel chemical probe for investigating cell signaling through GPCRs/G proteins.

Structural basis of cyclic oligoadenylate binding to the transcription factor Csa3 outlines cross talk between type III and type I CRISPR systems

RNA interference by type III CRISPR systems results in the synthesis of cyclic oligoadenylate (cOA) second messengers, which are known to bind and regulate various CARF domain-containing nuclease receptors. The CARF domain-containing Csa3 family of transcriptional factors associated with the DNA-targeting type I CRISPR systems regulate expression of various CRISPR and DNA repair genes in many prokaryotes. In this study, we extend the known receptor repertoire of cOA messengers to include transcriptional factors by demonstrating specific binding of cyclic tetra-adenylate (cA4) to Saccharolobus solfataricus Csa3 (Csa3Sso). Our 2.0-Å resolution X-ray crystal structure of cA4-bound full-length Csa3Sso reveals the binding of its CARF domain to an elongated conformation of cA4. Using cA4 binding affinity analyses of Csa3Sso mutants targeting the observed Csa3Sso•cA4 structural interface, we identified a Csa3-specific cA4 binding motif distinct from a more widely conserved cOA-binding CARF motif. Using a rational surface engineering approach, we increased the cA4 binding affinity of Csa3Sso up to ∼145-fold over the wildtype, which has potential applications for future second messenger-driven CRISPR gene expression and editing systems. Our in-solution Csa3Sso structural analysis identified cA4-induced allosteric and asymmetric conformational rearrangement of its C-terminal winged helix-turn-helix effector domains, which could potentially be incompatible to DNA binding. However, specific in vitro binding of the purified Csa3Sso to its putative promoter (PCas4a) was found to be cA4 independent, suggesting a complex mode of Csa3Sso regulation. Overall, our results support cA4-and Csa3-mediated cross talk between type III and type I CRISPR systems.

Conserved C-terminal motifs in odorant receptors instruct their cell surface expression and cAMP signaling

The highly individual plasma membrane expression and cAMP signaling of odorant receptors have hampered their ligand assignment and functional characterization in test cell systems. Chaperones have been identified to support the cell surface expression of only a portion of odorant receptors, with mechanisms remaining unclear. The presence of amino acid motifs that might be responsible for odorant receptors' individual intracellular retention or cell surface expression, and thus, for cAMP signaling, is under debate: so far, no such protein motifs have been suggested. Here, we demonstrate the existence of highly conserved C-terminal amino acid motifs, which discriminate at least between class-I and class-II odorant receptors, with their numbers of motifs increasing during evolution, by comparing C-terminal protein sequences from 4808 receptors across eight species. Truncation experiments and mutation analysis of C-terminal motifs, largely overlapping with helix 8, revealed single amino acids and their combinations to have differential impact on the cell surface expression and on stimulus-dependent cAMP signaling of odorant receptors in NxG 108CC15 cells. Our results demonstrate class-specific and individual C-terminal motif equipment of odorant receptors, which instruct their functional expression in a test cell system, and in situ may regulate their individual cell surface expression and intracellular cAMP signaling.

Time series computational prediction of vaccines for influenza A H3N2 with recurrent neural networks

Influenza viruses are persistently threatening public health, causing annual epidemics and sporadic pandemics due to rapid viral evolution. Vaccines are used to prevent influenza infections but the composition of the influenza vaccines have to be updated regularly to ensure its efficacy. Computational tools and analyses have become increasingly important in guiding the process of vaccine selection. By constructing time-series training samples with splittings and embeddings, we develop a computational method for predicting suitable strains as the recommendation of the influenza vaccines using recurrent neural networks (RNNs). The Encoder-decoder architecture of RNN model enables us to perform sequence-to-sequence prediction. We employ this model to predict the prevalent sequence of the H3N2 viruses sampled from 2006 to 2017. The identity between our predicted sequence and recommended vaccines is greater than 98% and the Pepitope<0.2 indicates their antigenic similarity. The multi-step vaccine prediction further demonstrates the robustness of our method which achieves comparable results in contrast to single step prediction. The results show significant matches of the recommended vaccine strains to the circulating strains. We believe it would facilitate the process of vaccine selection and surveillance of seasonal influenza epidemics.

Analysis of the SNARE Stx8 recycling reveals that the retromer-sorting motif has undergone evolutionary divergence

Fsv1/Stx8 is a Schizosaccharomyces pombe protein similar to mammalian syntaxin 8. stx8Δ cells are sensitive to salts, and the prevacuolar endosome (PVE) is altered in stx8Δ cells. These defects depend on the SNARE domain, data that confirm the conserved function of syntaxin8 and Stx8 in vesicle fusion at the PVE. Stx8 localizes at the trans-Golgi network (TGN) and the prevacuolar endosome (PVE), and its recycling depends on the retromer component Vps35, and on the sorting nexins Vps5, Vps17, and Snx3. Several experimental approaches demonstrate that Stx8 is a cargo of the Snx3-retromer. Using extensive truncation and alanine scanning mutagenesis, we identified the Stx8 sorting signal. This signal is an IEMeaM sequence that is located in an unstructured protein region, must be distant from the transmembrane (TM) helix, and where the ¹³³I, ¹³⁴E, ¹³⁵M, and ¹³⁸M residues are all essential for recycling. This sorting motif is different from those described for most retromer cargoes, which include aromatic residues, and resembles the sorting motif of mammalian polycystin-2 (PC2). Comparison of Stx8 and PC2 motifs leads to an IEMxx(I/M) consensus. Computer-assisted screening for this and for a loose Ψ(E/D)ΨXXΨ motif (where Ψ is a hydrophobic residue with large aliphatic chain) shows that syntaxin 8 and PC2 homologues from other organisms bear variation of this motif. The phylogeny of the Stx8 sorting motifs from the Schizosaccharomyces species shows that their divergence is similar to that of the genus, showing that they have undergone evolutionary divergence. A preliminary analysis of the motifs in syntaxin 8 and PC2 sequences from various organisms suggests that they might have also undergone evolutionary divergence, what suggests that the presence of almost-identical motifs in Stx8 and PC2 might be a case of convergent evolution.

Eukaryotes possess membranous intracellular compartments, whose communication is essential for cellular homeostasis. Protein complexes that facilitate the generation, transport, and fusion of coated vesicles mediate this communication. Since alterations in these processes lead to human disease, their characterization is of biological and medical interest. Retromer is a protein complex that facilitates retrograde trafficking from the prevacuolar endosome to the Golgi, being essential for the functionality of the endolysosomal system. SNAREs are required for vesicle fusion and, after facilitating membrane merging, are supposed to return to their donor organelle for new rounds of fusion. However, little is known about this recycling. We have found that Stx8, a fungal SNARE similar to human syntaxin 8, is a retromer cargo, and have identified its retromer binding motif. Sequence screening and comparison has determined that this sorting motif is conserved mainly in fungal Stx8 sequences. Notably, this motif is similar to the retromer sorting motif that is present in a family of vertebrate ion transporters. Our initial phylogenetic analyses suggest that, although retromer and some of its cargoes are conserved, the sorting motif in the cargoes might have undergone evolutionary divergence.

Recognising the signals for endosomal trafficking

The endosomal compartment is a major sorting station controlling the balance between endocytic recycling and lysosomal degradation, and its homeostasis is emerging as a central factor in various neurodegenerative diseases such as Alzheimer's and Parkinson's. Membrane trafficking is generally coordinated by the recognition of specific signals in transmembrane protein cargos by different transport machineries. A number of different protein trafficking complexes are essential for sequence-specific recognition and retrieval of endosomal cargos, recycling them to other compartments and acting to counter-balance the default endosomal sorting complex required for transport-mediated degradation pathway. In this review, we provide a summary of the key endosomal transport machineries, and the molecular mechanisms by which different cargo sequences are specifically recognised.

NBCZone: Universal three-dimensional construction of eleven amino acids near the catalytic nucleophile and base in the superfamily of (chymo)trypsin-like serine fold proteases

(Chymo)trypsin-like serine fold proteases belong to the serine/cysteine proteases found in eukaryotes, prokaryotes, and viruses. Their catalytic activity is carried out using a triad of amino acids, a nucleophile, a base, and an acid. For this superfamily of proteases, we propose the existence of a universal 3D structure comprising 11 amino acids near the catalytic nucleophile and base - Nucleophile-Base Catalytic Zone (NBCZone). The comparison of NBCZones among 169 eukaryotic, prokaryotic, and viral (chymo)trypsin-like proteases suggested the existence of 15 distinct groups determined by the combination of amino acids located at two "key" structure-functional positions 54_T and 55_T near the catalytic base His57_T. Most eukaryotic and prokaryotic proteases fell into two major groups, [ST]A and TN. Usually, proteases of [ST]A group contain a disulfide bond between cysteines Cys42_T and Cys58_T of the NBCZone. In contrast, viral proteases were distributed among seven groups, and lack this disulfide bond. Furthermore, only the [ST]A group of eukaryotic proteases contains glycine at position 43_T, which is instrumental for activation of these enzymes. In contrast, due to the side chains of residues at position 43_T prokaryotic and viral proteases do not have the ability to carry out the structural transition of the eukaryotic zymogen-zyme type.

(ADP-ribosyl)hydrolases: structure, function, and biology

ADP-ribosylation is an intricate and versatile posttranslational modification involved in the regulation of a vast variety of cellular processes in all kingdoms of life. Its complexity derives from the varied range of different chemical linkages, including to several amino acid side chains as well as nucleic acids termini and bases, it can adopt. In this review, we provide an overview of the different families of (ADP-ribosyl)hydrolases. We discuss their molecular functions, physiological roles, and influence on human health and disease. Together, the accumulated data support the increasingly compelling view that (ADP-ribosyl)hydrolases are a vital element within ADP-ribosyl signaling pathways and they hold the potential for novel therapeutic approaches as well as a deeper understanding of ADP-ribosylation as a whole.

Tempel: time-series mutation prediction of influenza A viruses via attention-based recurrent neural networks

Motivation

Influenza viruses are persistently threatening public health, causing annual epidemics and sporadic pandemics. The evolution of influenza viruses remains to be the main obstacle in the effectiveness of antiviral treatments due to rapid mutations. The goal of this work is to predict whether mutations are likely to occur in the next flu season using historical glycoprotein hemagglutinin sequence data. One of the major challenges is to model the temporality and dimensionality of sequential influenza strains and to interpret the prediction results.

Results

In this article, we propose an efficient and robust time-series mutation prediction model (Tempel) for the mutation prediction of influenza A viruses. We first construct the sequential training samples with splittings and embeddings. By employing recurrent neural networks with attention mechanisms, Tempel is capable of considering the historical residue information. Attention mechanisms are being increasingly used to improve the performance of mutation prediction by selectively focusing on the parts of the residues. A framework is established based on Tempel that enables us to predict the mutations at any specific residue site. Experimental results on three influenza datasets show that Tempel can significantly enhance the predictive performance compared with widely used approaches and provide novel insights into the dynamics of viral mutation and evolution.

Availability and implementation

The datasets, source code and supplementary documents are available at: https://drive.google.com/drive/folders/15WULR5__6k47iRotRPl3H7ghi3RpeNXH.

Supplementary information

Supplementary data are available at Bioinformatics online.

Using the Object-Oriented PowerShell for Simple Proteomics Data Analysis

Scripting languages such as Python and Bash are appreciated for solving simple, everyday tasks in bioinformatics. A more recent, object-oriented command shell and scripting language, PowerShell, has many attractive features: an object-oriented interactive command line, fluent navigation and manipulation of XML files, ability to explore and consume Web services from the command line, consistent syntax and grammar, rich regular expressions, and advanced output formatting. The key difference between classical command shells and scripting languages, such as bash, and object-oriented ones, such as PowerShell, is that in the latter the result of a command is a structured object with inherited properties and methods rather than a simple stream of characters. Conveniently, PowerShell is included in all new releases of Microsoft Windows and is available for Linux and macOS, making any data processing script portable. In this chapter we demonstrate how PowerShell in particular allows easy interaction with mass spectrometry data in XML formats, connection to Web services for tools such as BLAST, and presentation of results as formatted text or graphics. These features make PowerShell much more than "yet another scripting language."

Comprehensive analysis of all evolutionary paths between two divergent PDZ domain specificities

To understand the molecular evolution of functional diversity in protein families, we comprehensively investigated the consequences of all possible mutation combinations separating two peptide-binding domains with highly divergent specificities. We analyzed the Erbin PDZ domain (Erbin-PDZ), which exhibits canonical type I specificity, and a synthetic Erbin-PDZ variant (E-14) that differs at six positions and exhibits an atypical specificity that closely resembles that of the natural Pdlim4 PDZ domain (Pdlim4-PDZ). We constructed a panel of 64 PDZ domains covering all possible transitions between Erbin-PDZ and E-14 (i.e., the panel contained variants with all possible combinations of either the Erbin-PDZ or E-14 sequence at the six differing positions). We assessed the specificity profiles of the 64 PDZ domains using a C-terminal phage-displayed peptide library containing all possible genetically encoded heptapeptides. The specificity profiles clustered into six distinct groups, showing that intermediate domains can be nodes for the evolution of divergent functions. Remarkably, three substitutions were sufficient to convert the specificity of Erbin-PDZ to that of Pdlim4-PDZ, whereas Pdlim4-PDZ contains 71 differences relative to Erbin-PDZ. X-ray crystallography revealed the structural basis for specificity transition: a single substitution in the center of the binding site, supported by contributions from auxiliary substitutions, altered the main chain conformation of the peptide ligand to resemble that of ligands bound to Pdlim4-PDZ. Our results show that a very small set of mutations can dramatically alter protein specificity, and these findings support the hypothesis whereby complex protein functions evolve by gene duplication followed by cumulative mutations.

A PLPPV sequence in the p8 region of Gag provides late domain function for mouse mammary tumor virus

The late (L) domain sequence used by mouse mammary tumor virus (MMTV) remains undefined. Similar to other L domain-containing proteins, MMTV p8 and p14NC proteins are monoubiquitinated, suggesting L domain function. Site-directed mutagenesis of p8, PLPPV, and p14NC, PLPPL, sequences in MMTV Gag revealed a requirement only for the PLPPV sequence in virion release in a position-dependent manner. Electron microscopy of a defective Gag mutant confirmed an L domain budding defect morphology. The equine infectious anemia virus (EIAV) YPDL core L domain sequence and PLPPV provided L domain function in reciprocal MMTV and EIAV Gag exchange mutants, respectively. Alanine scanning of the PLPPV sequence revealed a strict requirement for the valine residue but only minor requirements for any one of the other residues. Thus, PLPPV provides MMTV L domain function, representing a fourth type of retroviral L domain that enables MMTV Gag proteins to co-opt cellular budding pathways for release.

Minimotif Miner 4: a million peptide minimotifs and counting

Minimotif Miner (MnM) is a database and web system for analyzing short functional peptide motifs, termed minimotifs. We present an update to MnM growing the database from ∼300 000 to >1 000 000 minimotif consensus sequences and instances. This growth comes largely from updating data from existing databases and annotation of articles with high-throughput approaches analyzing different types of post-translational modifications. Another update is mapping human proteins and their minimotifs to know human variants from the dbSNP, build 150. Now MnM 4 can be used to generate mechanistic hypotheses about how human genetic variation affect minimotifs and outcomes. One example of the utility of the combined minimotif/SNP tool identifies a loss of function missense SNP in a ubiquitylation minimotif encoded in the excision repair cross-complementing 2 (ERCC2) nucleotide excision repair gene. This SNP reaches genome wide significance for many types of cancer and the variant identified with MnM 4 reveals a more detailed mechanistic hypothesis concerning the role of ERCC2 in cancer. Other updates to the web system include a new architecture with migration of the web system and database to Docker containers for better performance and management. Weblinks:minimotifminer.org and mnm.engr.uconn.edu

Peptide Sequence Influence on the Differentiation of Valine and Norvaline by Hot Electron Capture Dissociation

Isomeric amino acid residues such as valine (Val) and norvaline (Nva) are common in recombinant proteins. The mis-incorporation of Nva for leucine (Leu) causes heterogeneity and in some cases even toxicity. Previous studies have shown that hot electron capture dissociation (HECD) is able to differentiate Val from Nva by producing diagnostic w ions on custom designed synthetic model peptides. To broaden the utilization of HECD in proteomic studies and to define the critical structural features, a thorough investigation was performed on representative peptides including specifically designed synthetic peptides as well as biological peptides bearing tryptic digest-like features and peptides with post-translational modifications. Experimental evidence confirmed that the formation of a w ion is directly dependent upon the presence of the corresponding z ion. The results suggested that a charge carrier residue at the C-terminus is promoting the formation of diagnostic w ions for Nva. Thus, peptides resulting from trypsin digestion, with arginine (Arg) or lysine (Lys) at the C-terminus, can be analyzed using the HECD method. Post-translational modification (PTM) such as phosphorylation did not prevent the generation of the requisite side chain fragmentation w ions. These results suggest the general applicability of HECD for unambiguous identification of Val and Nva especially in structure characterization of therapeutic proteins.

Sequential cleavage of the proteins encoded by HNOT/ALG3, the human counterpart of the Drosophila NOT and yeast ALG3 gene, results in products acting in distinct cellular compartments

This study provides first insights into the biosynthesis, structure, biochemistry and complex processing of the proteins encoded by hNOT/ALG3, the human counterpart of the Drosophila Neighbour of TID (NOT) and the yeast asparagine linked glycosylation 3 gene (ALG3), which encodes a mannosyltransferase. Unambiguous evidence that both the fly and human proteins act as mannosyltransferases has not been provided yet. Previously, we showed that hNOT/ALG3 encodes two alternatively spliced main transcripts, hNOT-1/ALG3-1 and hNOT-4/ALG3-4, and their 15 truncated derivatives that lack diverse sets of exons and/or carry point mutations that result in premature termination codons. Here we show that the truncated transcripts are not translated. The two main forms hNOT-1/ALG3-1 and -4, distinguishable by alternative exon 1, encode full-length precursors that undergo a complex posttranslational processing. To specifically detect the two full-length hNOT/ALG3 proteins and their distinct derivatives and to examine their expression profiles and cellular location we generated polyclonal antibodies against diverse parts of the putative full-length proteins. We provide experimental evidence for the N-glycosylation of the two precursors. This modification seems to be a prerequisite for their sequential cleavage resulting in derivatives destined to distinct cellular compartments and links them with the N-glycosylation machinery not as its functional component but as molecules functionally dependent on its action. We present the expression profiles and subcellular location of the two full-length proteins, their N-glycosylated forms and distinct cleavage products. Furthermore, using diverse bioinformatics tools, we characterize the properties and predict the 2D and 3D structure of the two proteins and, for comparative purposes, of their Drosophila counterpart.

PSSMSearch: a server for modeling, visualization, proteome-wide discovery and annotation of protein motif specificity determinants

There is a pressing need for in silico tools that can aid in the identification of the complete repertoire of protein binding (SLiMs, MoRFs, miniMotifs) and modification (moiety attachment/removal, isomerization, cleavage) motifs. We have created PSSMSearch, an interactive web-based tool for rapid statistical modeling, visualization, discovery and annotation of protein motif specificity determinants to discover novel motifs in a proteome-wide manner. PSSMSearch analyses proteomes for regions with significant similarity to a motif specificity determinant model built from a set of aligned motif-containing peptides. Multiple scoring methods are available to build a position-specific scoring matrix (PSSM) describing the motif specificity determinant model. This model can then be modified by a user to add prior knowledge of specificity determinants through an interactive PSSM heatmap. PSSMSearch includes a statistical framework to calculate the significance of specificity determinant model matches against a proteome of interest. PSSMSearch also includes the SLiMSearch framework's annotation, motif functional analysis and filtering tools to highlight relevant discriminatory information. Additional tools to annotate statistically significant shared keywords and GO terms, or experimental evidence of interaction with a motif-recognizing protein have been added. Finally, PSSM-based conservation metrics have been created for taxonomic range analyses. The PSSMSearch web server is available at http://slim.ucd.ie/pssmsearch/.

Analysis of MAPK and MAPKK gene families in wheat and related Triticeae species

Background

The mitogen-activated protein kinase (MAPK) family is involved in signal transduction networks that underpin many different biological processes in plants, ranging from development to biotic and abiotic stress responses. To date this class of enzymes has received little attention in Triticeae species, which include important cereal crops (wheat, barley, rye and triticale) that represent over 20% of the total protein food-source worldwide.

Results

The work presented here focuses on two subfamilies of Triticeae MAPKs, the MAP kinases (MPKs), and the MAPK kinases (MKKs) whose members phosphorylate the MPKs. In silico analysis of multiple Triticeae sequence databases led to the identification of 152 MAPKs belonging to these two sub-families. Some previously identified MAPKs were renamed to reflect the literature consensus on MAPK nomenclature. Two novel MPKs, MPK24 and MPK25, have been identified, including the first example of a plant MPK carrying the TGY activation loop sequence common to mammalian p38 MPKs. An EF-hand calcium-binding domain was found in members of the Triticeae MPK17 clade, a feature that appears to be specific to Triticeae species. New insights into the novel MEY activation loop identified in MPK11s are offered. When the exon-intron patterns for some MPKs and MKKs of wheat, barley and ancestors of wheat were assembled based on transcript data in GenBank, they showed deviations from the same sequence predicted in Ensembl. The functional relevance of MAPKs as derived from patterns of gene expression, MPK activation and MKK-MPK interaction is discussed.

Conclusions

A comprehensive resource of accurately annotated and curated Triticeae MPK and MKK sequences has been created for wheat, barley, rye, triticale, and two ancestral wheat species, goat grass and red wild einkorn. The work we present here offers a central information resource that will resolve existing confusion in the literature and sustain expansion of MAPK research in the crucial Triticeae grains.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4545-9) contains supplementary material, which is available to authorized users.

Modulation of ABA Signaling by Altering VxGΦL Motif of PP2Cs in Oryza sativa

The abscisic acid (ABA) signaling pathway is regulated by clade A type 2C protein phosphatases (PP2CAs) in plants. In the presence of ABA, PP2Cs release stress/ABA-activated protein kinases by binding to ABA-bound receptors (PYL/RCARs) for activation. Although the wedging tryptophan in PP2Cs is critical in the interaction with PYL/RCARs in Arabidopsis and rice, it remains elusive as to how other interface regions are involved in the interaction. Here, we report the identification of a conserved region on PP2Cs, termed the VxGΦL motif, which modulates the interaction with PYL/RCARs through its second and fourth residues. The effects of the second and fourth residues on the interaction of OsPP2C50 with several OsPYL/RCAR proteins were investigated by systematic mutagenesis. One OsPP2C50 mutant, VFGML ("FM") mutant, lowered the affinity to OsPYL/RCAR3 by ∼15-fold in comparison with the wild-type. Comparison of the crystal structures of wild-type OsPP2C50:ABA:OsPYL/RCAR3 with those composed of FM mutant revealed local conformational changes near the VxGΦL motif, further supported by hydrogen-deuterium exchange mass spectrometry. In rice protoplasts, ABA signaling was altered by mutations in the VxGΦL motif. Transgenic Arabidopsis plants overexpressing OsPP2C50 and OsPP2C50FM showed altered ABA sensitivity. Taken together, the VxGΦL motif of PP2Cs appears to modulate the affinity of PP2Cs with PYL/RCARs and thus likely to alter the ABA signaling, leading to the differential sensitivity to ABA in planta.

Secreted Effectors Encoded within and outside of the Francisella Pathogenicity Island Promote Intramacrophage Growth

The intracellular bacterial pathogen Francisella tularensis causes tularemia, a zoonosis that can be fatal. The type VI secretion system (T6SS) encoded by the Francisella pathogenicity island (FPI) is critical for the virulence of this organism. Existing studies suggest that the complete repertoire of T6SS effectors delivered to host cells is encoded by the FPI. Using a proteome-wide approach, we discovered that the FPI-encoded T6SS exports at least three effectors encoded outside of the island. These proteins share features with virulence determinants of other pathogens, and we provide evidence that they can contribute to intramacrophage growth. The remaining proteins that we identified are encoded within the FPI. Two of these FPI-encoded proteins constitute effectors, whereas the others form a unique complex required for core function of the T6SS apparatus. The discovery of secreted effectors mediating interactions between Francisella and its host significantly advances our understanding of the pathogenesis of this organism.

The AAA+ ATPase p97, a cellular multitool

The AAA+ (ATPases associated with diverse cellular activities) ATPase p97 is essential to a wide range of cellular functions, including endoplasmic reticulum-associated degradation, membrane fusion, NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells) activation and chromatin-associated processes, which are regulated by ubiquitination. p97 acts downstream from ubiquitin signaling events and utilizes the energy from ATP hydrolysis to extract its substrate proteins from cellular structures or multiprotein complexes. A multitude of p97 cofactors have evolved which are essential to p97 function. Ubiquitin-interacting domains and p97-binding domains combine to form bi-functional cofactors, whose complexes with p97 enable the enzyme to interact with a wide range of ubiquitinated substrates. A set of mutations in p97 have been shown to cause the multisystem proteinopathy inclusion body myopathy associated with Paget's disease of bone and frontotemporal dementia. In addition, p97 inhibition has been identified as a promising approach to provoke proteotoxic stress in tumors. In this review, we will describe the cellular processes governed by p97, how the cofactors interact with both p97 and its ubiquitinated substrates, p97 enzymology and the current status in developing p97 inhibitors for cancer therapy.

It All Starts with a Sandwich: Identification of Sialidases with Trans-Glycosylation Activity

Sialidases (3.2.1.18) may exhibit trans-sialidase activity to catalyze sialylation of lactose if the active site topology is congruent with that of the Trypanosoma cruzi trans-sialidase (EC 2.4.1.-). The present work was undertaken to test the hypothesis that a particular aromatic sandwich structure of two amino acids proximal to the active site of the T. cruzi trans-sialidase infers trans-sialidase activity. On this basis, four enzymes with putative trans-sialidase activity were identified through an iterative alignment from 2909 native sialidases available in GenBank, which were cloned and expressed in Escherichia coli. Of these, one enzyme, SialH, derived from Haemophilus parasuis had an aromatic sandwich structure on the protein surface facing the end of the catalytic site (Phe168; Trp366), and was indeed found to exhibit trans-sialidase activity. SialH catalyzed production of the human milk oligosaccharide 3’-sialyllactose as well as the novel trans-sialylation product 3-sialyllactose using casein glycomacropeptide as sialyl donor and lactose as acceptor. The findings corroborated that Tyr119 and Trp312 in the T. cruzi trans-sialidase are part of an aromatic sandwich structure that confers trans-sialylation activity for lactose sialylation. The in silico identification of trans-glycosidase activity by rational active site topology alignment thus proved to be a quick tool for selecting putative trans-sialidases amongst a large group of glycosyl hydrolases. The approach moreover provided data that help understand structure-function relations of trans-sialidases.

Structural Basis of the Recruitment of Ubiquitin-specific Protease USP15 by Spliceosome Recycling Factor SART3

Ubiquitin-specific proteases (USPs) USP15 and USP4 belong to a subset of USPs featuring an N-terminal tandem domain in USP (DUSP) and ubiquitin-like (UBL) domain. Squamous cell carcinoma antigen recognized by T-cell 3 (SART3), a spliceosome recycling factor, binds to the DUSP-UBL domain of USP15 and USP4, recruiting them to the nucleus from the cytosol to control deubiquitination of histone H2B and spliceosomal proteins, respectively. To provide structural insight, we solved crystal structures of SART3 in the apo-form and in complex with the DUSP-UBL domain of USP15 at 2.0 and 3.0 Å, respectively. Structural analysis reveals SART3 contains 12 half-a-tetratricopeptide (HAT) repeats, organized into two subdomains, HAT-N and HAT-C. SART3 dimerizes through the concave surface of HAT-C, whereas the HAT-C convex surface binds USP15 in a novel bipartite mode. Isothermal titration calorimetry measurements and mutagenesis analysis confirmed key residues of USP15 involved in the interaction and indicated USP15 binds 20-fold stronger than USP4.

The Functional Human C-Terminome

All translated proteins end with a carboxylic acid commonly called the C-terminus. Many short functional sequences (minimotifs) are located on or immediately proximal to the C-terminus. However, information about the function of protein C-termini has not been consolidated into a single source. Here, we built a new "C-terminome" database and web system focused on human proteins. Approximately 3,600 C-termini in the human proteome have a minimotif with an established molecular function. To help evaluate the function of the remaining C-termini in the human proteome, we inferred minimotifs identified by experimentation in rodent cells, predicted minimotifs based upon consensus sequence matches, and predicted novel highly repetitive sequences in C-termini. Predictions can be ranked by enrichment scores or Gene Evolutionary Rate Profiling (GERP) scores, a measurement of evolutionary constraint. By searching for new anchored sequences on the last 10 amino acids of proteins in the human proteome with lengths between 3-10 residues and up to 5 degenerate positions in the consensus sequences, we have identified new consensus sequences that predict instances in the majority of human genes. All of this information is consolidated into a database that can be accessed through a C-terminome web system with search and browse functions for minimotifs and human proteins. A known consensus sequence-based predicted function is assigned to nearly half the proteins in the human proteome. Weblink: http://cterminome.bio-toolkit.com.

KATP Channels in the Cardiovascular System

KATP channels are integral to the functions of many cells and tissues. The use of electrophysiological methods has allowed for a detailed characterization of KATP channels in terms of their biophysical properties, nucleotide sensitivities, and modification by pharmacological compounds. However, even though they were first described almost 25 years ago (Noma 1983, Trube and Hescheler 1984), the physiological and pathophysiological roles of these channels, and their regulation by complex biological systems, are only now emerging for many tissues. Even in tissues where their roles have been best defined, there are still many unanswered questions. This review aims to summarize the properties, molecular composition, and pharmacology of KATP channels in various cardiovascular components (atria, specialized conduction system, ventricles, smooth muscle, endothelium, and mitochondria). We will summarize the lessons learned from available genetic mouse models and address the known roles of KATP channels in cardiovascular pathologies and how genetic variation in KATP channel genes contribute to human disease.

How pathogens use linear motifs to perturb host cell networks

Molecular mimicry is one of the powerful stratagems that pathogens employ to colonise their hosts and take advantage of host cell functions to guarantee their replication and dissemination. In particular, several viruses have evolved the ability to interact with host cell components through protein short linear motifs (SLiMs) that mimic host SLiMs, thus facilitating their internalisation and the manipulation of a wide range of cellular networks. Here we present convincing evidence from the literature that motif mimicry also represents an effective, widespread hijacking strategy in prokaryotic and eukaryotic parasites. Further insights into host motif mimicry would be of great help in the elucidation of the molecular mechanisms behind host cell invasion and the development of anti-infective therapeutic strategies.

Phylogenetic analysis of glycerol 3-phosphate acyltransferases in opisthokonts reveals unexpected ancestral complexity and novel modern biosynthetic components

Glycerolipid synthesis represents a central metabolic process of all forms of life. In the last decade multiple genes coding for enzymes responsible for the first step of the pathway, catalyzed by glycerol 3-phosphate acyltransferase (GPAT), have been described, and characterized primarily in model organisms like Saccharomyces cerevisiae and mice. Notoriously, the fungal enzymes share low sequence identity with their known animal counterparts, and the nature of their homology is unclear. Furthermore, two mitochondrial GPAT isoforms have been described in animal cells, while no such enzymes have been identified in Fungi. In order to determine if the yeast and mammalian GPATs are representative of the set of enzymes present in their respective groups, and to test the hypothesis that metazoan orthologues are indeed absent from the fungal clade, a comparative genomic and phylogenetic analysis was performed including organisms spanning the breadth of the Opisthokonta supergroup. Surprisingly, our study unveiled the presence of ‘fungal’ orthologs in the basal taxa of the holozoa and ‘animal’ orthologues in the basal holomycetes. This includes a novel clade of fungal homologues, with putative peroxisomal targeting signals, of the mitochondrial/peroxisomal acyltransferases in Metazoa, thus potentially representing an undescribed metabolic capacity in the Fungi. The overall distribution of GPAT homologues is suggestive of high relative complexity in the ancestors of the opisthokont clade, followed by loss and sculpting of the complement in the descendent lineages. Divergence from a general versatile metabolic model, present in ancestrally deduced GPAT complements, points to distinctive contributions of each GPAT isoform to lipid metabolism and homeostasis in contemporary organisms like humans and their fungal pathogens.

The development and application of a quantitative peptide microarray based approach to protein interaction domain specificity space

Protein interaction domain (PID) linear peptide motif interactions direct diverse cellular processes in a specific and coordinated fashion. PID specificity, or the interaction selectivity derived from affinity preferences between possible PID-peptide pairs is the basis of this ability. Here, we develop an integrated experimental and computational cellulose peptide conjugate microarray (CPCMA) based approach for the high throughput analysis of PID specificity that provides unprecedented quantitative resolution and reproducibility. As a test system, we quantify the specificity preferences of four Src Homology 2 domains and 124 physiological phosphopeptides to produce a novel quantitative interactome. The quantitative data set covers a broad affinity range, is highly precise, and agrees well with orthogonal biophysical validation, in vivo interactions, and peptide library trained algorithm predictions. In contrast to preceding approaches, the CPCMAs proved capable of confidently assigning interactions into affinity categories, resolving the subtle affinity contributions of residue correlations, and yielded predictive peptide motif affinity matrices. Unique CPCMA enabled modes of systems level analysis reveal a physiological interactome with expected node degree value decreasing as a function of affinity, resulting in minimal high affinity binding overlap between domains; uncover that Src Homology 2 domains bind ligands with a similar average affinity yet strikingly different levels of promiscuity and binding dynamic range; and parse with unprecedented quantitative resolution contextual factors directing specificity. The CPCMA platform promises broad application within the fields of PID specificity, synthetic biology, specificity focused drug design, and network biology.

Mix and match: investigating heteromeric and heterotypic gap junction channels in model systems and native tissues

This review is based in part on a roundtable discussion session: "Physiological roles for heterotypic/heteromeric channels" at the 2013 International Gap Junction Conference (IGJC 2013) in Charleston, South Carolina. It is well recognized that multiple connexins can specifically co-assemble to form mixed gap junction channels with unique properties as a means to regulate intercellular communication. Compatibility determinants for both heteromeric and heterotypic gap junction channel formation have been identified and associated with specific connexin amino acid motifs. Hetero-oligomerization is also a regulated process; differences in connexin quality control and monomer stability are likely to play integral roles to control interactions between compatible connexins. Gap junctions in oligodendrocyte:astrocyte communication and in the cardiovascular system have emerged as key systems where heterotypic and heteromeric channels have unique physiologic roles. There are several methodologies to study heteromeric and heterotypic channels that are best applied to either heterologous expression systems, native tissues or both. There remains a need to use and develop different experimental approaches in order to understand the prevalence and roles for mixed gap junction channels in human physiology.