Pleiotropic functions of a conserved insect-specific Hox peptide motif

The proteins that regulate developmental processes in animals have generally been well conserved during evolution. A few cases are known where protein activities have functionally evolved. These rare examples raise the issue of how highly conserved regulatory proteins with many roles evolve new functions while maintaining old functions. We have investigated this by analyzing the function of the ;QA' peptide motif of the Hox protein Ultrabithorax (Ubx), a motif that has been conserved throughout insect evolution since its establishment early in the lineage. We precisely deleted the QA motif at the endogenous locus via allelic replacement in Drosophila melanogaster. Although the QA motif was originally characterized as involved in the repression of limb formation, we have found that it is highly pleiotropic. Curiously, deleting the QA motif had strong effects in some tissues while barely affecting others, suggesting that QA function is preferentially required for a subset of Ubx target genes. QA deletion homozygotes had a normal complement of limbs, but, at reduced doses of Ubx and the abdominal-A (abd-A) Hox gene, ectopic limb primordia and adult abdominal limbs formed when the QA motif was absent. These results show that redundancy and the additive contributions of activity-regulating peptide motifs play important roles in moderating the phenotypic consequences of Hox protein evolution, and that pleiotropic peptide motifs that contribute quantitatively to several functions are subject to intense purifying selection.     

The conservation profile of a protein bears the imprint of the molecule that is evolutionarily coupled to the protein

The conservation profile of a protein is a curve of the conservation levels of amino acids along the sequence. Biologists are usually more interested in individual points on the curve (namely, the conserved amino acids) than the overall shape of the curve. Here, we show that the conservation curves of proteins bear the imprints of molecules that are evolutionarily coupled to the proteins. Our method is based on recent studies that a sequence conservation profile is quantitatively linked to its structural packing profile. We find that the conservation profiles of nucleic acid (NA) binding proteins are better correlated with the packing profiles of the protein–NA complexes than those of the proteins alone. This indicates that a nucleic acid binding protein evolves to accommodate the nucleic acid in such a way that the residues involved in binding have their conservation levels closely coupled with the specific nucleotides. Proteins 2015; 83:1407–1413. © 2015 Wiley Periodicals, Inc.     

Proteome-wide analysis of single-nucleotide variations in the N-glycosylation sequon of human genes

N-linked glycosylation is one of the most frequent post-translational modifications of proteins with a profound impact on their biological function. Besides other functions, N-linked glycosylation assists in protein folding, determines protein orientation at the cell surface, or protects proteins from proteases. The N-linked glycans attach to asparagines in the sequence context Asn-X-Ser/Thr, where X is any amino acid except proline. Any variation (e.g. non-synonymous single nucleotide polymorphism or mutation) that abolishes the N-glycosylation sequence motif will lead to the loss of a glycosylation site. On the other hand, variations causing a substitution that creates a new N-glycosylation sequence motif can result in the gain of glycosylation. Although the general importance of glycosylation is well known and acknowledged, the effect of variation on the actual glycoproteome of an organism is still mostly unknown. In this study, we focus on a comprehensive analysis of non-synonymous single nucleotide variations (nsSNV) that lead to either loss or gain of the N-glycosylation motif. We find that 1091 proteins have modified N-glycosylation sequons due to nsSNVs in the genome. Based on analysis of proteins that have a solved 3D structure at the site of variation, we find that 48% of the variations that lead to changes in glycosylation sites occur at the loop and bend regions of the proteins. Pathway and function enrichment analysis show that a significant number of proteins that gained or lost the glycosylation motif are involved in kinase activity, immune response, and blood coagulation. A structure-function analysis of a blood coagulation protein, antithrombin III and a protease, cathepsin D, showcases how a comprehensive study followed by structural analysis can help better understand the functional impact of the nsSNVs.     

Helicase motif Ia is involved in single-strand DNA-binding and helicase activities of the herpes simplex virus type 1 origin-binding protein,UL9

UL9 is a multifunctional protein essential for herpes simplex virus type 1 (HSV-1) replication in vivo. UL9 is a member of the superfamily II helicases and exhibits helicase and origin-binding activities. It is thought that UL9 binds the origin of replication and unwinds it in the presence of ATP and the HSV-1 single-stranded DNA (ssDNA)-binding protein. We have previously characterized the biochemical properties of mutants in all helicase motifs except for motif Ia (B. Marintcheva and S. Weller, J. Biol. Chem. 276:6605-6615, 2001). Structural information for other superfamily I and II helicases indicates that motif Ia is involved in ssDNA binding. By analogy, we hypothesized that UL9 motif Ia is important for the ssDNA-binding function of the protein. On the basis of sequence conservation between several UL9 homologs within the Herpesviridae family and distant homology with helicases whose structures have been solved, we designed specific mutations in motif Ia and analyzed them genetically and biochemically. Mutant proteins with residues predicted to be involved in ssDNA binding (R112A and R113A/F115A) exhibited wild-type levels of intrinsic ATPase activity and moderate to severe defects in ssDNA-stimulated ATPase activity and ssDNA binding. The S110T mutation targets a residue not predicted to contact ssDNA directly. The mutant protein with this mutation exhibited wild-type levels of intrinsic ATPase activity and near wild-type levels of ssDNA-stimulated ATPase activity and ssDNA binding. All mutant proteins lack helicase activity but were able to dimerize and bind the HSV-1 origin of replication as well as wild-type UL9. Our results indicate that residues from motif Ia contribute to the ssDNA-binding and helicase activities of UL9 and are essential for viral growth. This work represents the successful application of an approach based on a combination of bioinformatics and structural information from related proteins to deduce valuable information about a protein of interest.     

Prediction of short linear protein binding regions

Short linear motifs in proteins (typically 3-12 residues in length) play key roles in protein-protein interactions by frequently binding specifically to peptide binding domains within interacting proteins. Their tendency to be found in disordered segments of proteins has meant that they have often been overlooked. Here we present SLiMPred (short linear motif predictor), the first general de novo method designed to computationally predict such regions in protein primary sequences independent of experimentally defined homologs and interactors. The method applies machine learning techniques to predict new motifs based on annotated instances from the Eukaryotic Linear Motif database, as well as structural, biophysical, and biochemical features derived from the protein primary sequence. We have integrated these data sources and benchmarked the predictive accuracy of the method, and found that it performs equivalently to a predictor of protein binding regions in disordered regions, in addition to having predictive power for other classes of motif sites such as polyproline II helix motifs and short linear motifs lying in ordered regions. It will be useful in predicting peptides involved in potential protein associations and will aid in the functional characterization of proteins, especially of proteins lacking experimental information on structures and interactions. We conclude that, despite the diversity of motif sequences and structures, SLiMPred is a valuable tool for prioritizing potential interaction motifs in proteins.     

Slipknotted and unknotted monovalent cation-proton antiporters evolved from a common ancestor

While the slipknot topology in proteins has been known for over a decade, its evolutionary origin is still a mystery. We have identified a previously overlooked slipknot motif in a family of two-domain membrane transporters. Moreover, we found that these proteins are homologous to several families of unknotted membrane proteins. This allows us to directly investigate the evolution of the slipknot motif. Based on our comprehensive analysis of 17 distantly related protein families, we have found that slipknotted and unknotted proteins share a common structural motif. Furthermore, this motif is conserved on the sequential level as well. Our results suggest that, regardless of topology, the proteins we studied evolved from a common unknotted ancestor single domain protein. Our phylogenetic analysis suggests the presence of at least seven parallel evolutionary scenarios that led to the current diversity of proteins in question. The tools we have developed in the process can now be used to investigate the evolution of other repeated-domain proteins.     

The Arabidopsis TRM1-TON1 interaction reveals a recruitment network common to plant cortical microtubule arrays and eukaryotic centrosomes

Land plant cells assemble microtubule arrays without a conspicuous microtubule organizing center like a centrosome. In Arabidopsis thaliana, the TONNEAU1 (TON1) proteins, which share similarity with FOP, a human centrosomal protein, are essential for microtubule organization at the cortex. We have identified a novel superfamily of 34 proteins conserved in land plants, the TON1 Recruiting Motif (TRM) proteins, which share six short conserved motifs, including a TON1-interacting motif present in all TRMs. An archetypal member of this family, TRM1, is a microtubule-associated protein that localizes to cortical microtubules and binds microtubules in vitro. Not all TRM proteins can bind microtubules, suggesting a diversity of functions for this family. In addition, we show that TRM1 interacts in vivo with TON1 and is able to target TON1 to cortical microtubules via its C-terminal TON1 interaction motif. Interestingly, three motifs of TRMs are found in CAP350, a human centrosomal protein interacting with FOP, and the C-terminal M2 motif of CAP350 is responsible for FOP recruitment at the centrosome. Moreover, we found that TON1 can interact with the human CAP350 M2 motif in yeast. Taken together, our results suggest conservation of eukaryotic centrosomal components in plant cells.     

The THAP domain: a novel protein motif with similarity to the DNA-binding domain of P element transposase

We have identified a novel evolutionarily conserved protein motif - designated the THAP domain - that defines a new family of cellular factors. We have found that the THAP domain presents striking similarities with the site-specific DNA-binding domain (DBD) of Drosophila P element transposase, including a similar size, N-terminal location, and conservation of the residues that define the THAP motif, such as the C2CH signature (Cys-Xaa(2-4)-Cys-Xaa(35-50)-Cys-Xaa(2)-His). Our results suggest that the THAP domain is a novel example of a DBD that is shared between cellular proteins and transposases from mobile genomic parasites.     

PAS domains. Common structure and common flexibility

PAS (PER-ARNT-SIM) domains are a family of sensor protein domains involved in signal transduction in a wide range of organisms. Recent structural studies have revealed that these domains contain a structurally conserved alpha/beta-fold, whereas almost no conservation is observed at the amino acid sequence level. The photoactive yellow protein, a bacterial light sensor, has been proposed as the PAS structural prototype yet contains an N-terminal helix-turn-helix motif not found in other PAS domains. Here we describe the atomic resolution structure of a photoactive yellow protein deletion mutant lacking this motif, revealing that the PAS domain is indeed able to fold independently and is not affected by the removal of these residues. Computer simulations of currently known PAS domain structures reveal that these domains are not only structurally conserved but are also similar in their conformational flexibilities. The observed motions point to a possible common mechanism for communicating ligand binding/activation to downstream transducer proteins.     

KH domain: one motif, two folds

The K homology (KH) module is a widespread RNA-binding motif that has been detected by sequence similarity searches in such proteins as heterogeneous nuclear ribonucleoprotein K (hnRNP K) and ribosomal protein S3. Analysis of spatial structures of KH domains in hnRNP K and S3 reveals that they are topologically dissimilar and thus belong to different protein folds. Thus KH motif proteins provide a rare example of protein domains that share significant sequence similarity in the motif regions but possess globally distinct structures. The two distinct topologies might have arisen from an ancestral KH motif protein by N- and C-terminal extensions, or one of the existing topologies may have evolved from the other by extension, displacement and deletion. C-terminal extension (deletion) requires β-sheet rearrangement through the insertion (removal) of a β-strand in a manner similar to that observed in serine protease inhibitors serpins. Current analysis offers a new look on how proteins can change fold in the course of evolution.     

A conserved motif controls nuclear localization of Drosophila Muscleblind

Human Muscleblind-like proteins are alternative splicing regulators that are functionally altered in the RNA-mediated disease myotonic dystrophy. There are different Muscleblind protein isoforms in Drosophila and we previously determined that these have different subcellular localizations in the COS-M6 cell line. Here, we describe the conservation of the sequence motif KRAEK in isoforms C and E and propose a specific function for this motif. Different Muscleblind isoforms localize to the peri-plasma membrane (MblA), cytoplasm (MblB), or show no preference for the nuclear or cytoplasmic compartment (MblC and MblD) in Drosophila S2 cells transiently transfected with Musclebind expression plasmids. Mutation of the KRAEK motif reduces MblC nuclear localization, whereas fusion of a single KRAEK motif to the heterologous protein β-galactosidase is sufficient to target the reporter protein to the nucleus of S2 cells. This motif is not exclusive to Muscleblind proteins and is detected in several other protein types. Taken together, these results suggest that the KRAEK motif regulates nuclear translocation of Muscleblind and may constitute a new class of nuclear localization signal.     

Protein–protein interactions leave evolutionary footprints: High molecular coevolution at the core of interfaces

Protein–protein interactions are essential to all aspects of life. Specific interactions result from evolutionary pressure at the interacting interfaces of partner proteins. However, evolutionary pressure is not homogeneous within the interface: for instance, each residue does not contribute equally to the binding energy of the complex. To understand functional differences between residues within the interface, we analyzed their properties in the core and rim regions. Here, we characterized protein interfaces with two evolutionary measures, conservation and coevolution, using a comprehensive dataset of 896 protein complexes. These scores can detect different selection pressures at a given position in a multiple sequence alignment. We also analyzed how the number of interactions in which a residue is involved influences those evolutionary signals. We found that the coevolutionary signal is higher in the interface core than in the interface rim region. Additionally, the difference in coevolution between core and rim regions is comparable to the known difference in conservation between those regions. Considering proteins with multiple interactions, we found that conservation and coevolution increase with the number of different interfaces in which a residue is involved, suggesting that more constraints (i.e., a residue that must satisfy a greater number of interactions) allow fewer sequence changes at those positions, resulting in higher conservation and coevolution values. These findings shed light on the evolution of protein interfaces and provide information useful for identifying protein interfaces and predicting protein–protein interactions.     

Functional genomic analysis reveals the utility of the I/LWEQ module as a predictor of protein:actin interaction

The I/LWEQ module is a conserved sequence that we have identified as an actin-binding motif in the metazoan focal adhesion protein talin and the yeast protein Sla2p. Both of these proteins are associated with the actin cytoskeleton in cells. To better establish the value of the I/LWEQ module for prediction of actin-binding function, we have applied a functional genomics approach. Analysis of the 23 available I/LWEQ module sequences supports the division of I/LWEQ protein superfamily into four groups: (1) metazoan talin, (2) Dictyostelium discoideum talin homologs TalA/B, (3) metazoan Hip1p, and (4) yeast Sla2p. We show here that I/LWEQ modules from each major group bind to F-actin in vitro and that GFP-fusion proteins of the I/LWEQ modules of talin and Sla2p bind to F-actin in vivo. Therefore, the presence of an I/LWEQ module is strongly predictive of protein-actin interactions. The structural and functional conservation of the I/LWEQ module across the phylogenetic distance between cellular slime molds and mammals implies that the role of the I/LWEQ module is to connect diverse proteins involved in distinct cellular processes, including cell adhesion, cytoskeletal organization, and cell differentiation, to the actin cytoskeleton.