ST proteins, a new family of plant tandem repeat proteins with a DUF2775 domain mainly found in Fabaceae and Asteraceae

Background

Many proteins with tandem repeats in their sequence have been described and classified according to the length of the repeats: I) Repeats of short oligopeptides (from 2 to 20 amino acids), including structural cell wall proteins and arabinogalactan proteins. II) Repeats that range in length from 20 to 40 residues, including proteins with a well-established three-dimensional structure often involved in mediating protein-protein interactions. (III) Longer repeats in the order of 100 amino acids that constitute structurally and functionally independent units. Here we analyse ShooT specific (ST) proteins, a family of proteins with tandem repeats of unknown function that were first found in Leguminosae, and their possible similarities to other proteins with tandem repeats.

Results

ST protein sequences were only found in dicotyledonous plants, limited to several plant families, mainly the Fabaceae and the Asteraceae. ST mRNAs accumulate mainly in the roots and under biotic interactions. Most ST proteins have one or several Domain(s) of Unknown Function 2775 (DUF2775). All deduced ST proteins have a signal peptide, indicating that these proteins enter the secretory pathway, and the mature proteins have tandem repeat oligopeptides that share a hexapeptide (E/D)FEPRP followed by 4 partially conserved amino acids, which could determine a putative N-glycosylation signal, and a fully conserved tyrosine. In a phylogenetic tree, the sequences clade according to taxonomic group. A possible involvement in symbiosis and abiotic stress as well as in plant cell elongation is suggested, although different STs could play different roles in plant development.

Conclusions

We describe a new family of proteins called ST whose presence is limited to the plant kingdom, specifically to a few families of dicotyledonous plants. They present 20 to 40 amino acid tandem repeat sequences with different characteristics (signal peptide, DUF2775 domain, conservative repeat regions) from the described group of 20 to 40 amino acid tandem repeat proteins and also from known cell wall proteins with repeat sequences. Several putative roles in plant physiology can be inferred from the characteristics found.     

Mass spectrometric and kinetic analysis of ASF/SF2 phosphorylation by SRPK1 and Clk/Sty

Assembly of the spliceosome requires the participation of SR proteins, a family of splicing factors rich in arginine-serine dipeptide repeats. The repeat regions (RS domains) are polyphosphorylated by the SRPK and Clk/Sty families of kinases. The two families of kinases have distinct enzymatic properties, raising the question of how they may work to regulate the function of SR proteins in RNA metabolism in mammalian cells. Here we report the first mass spectral analysis of the RS domain of ASF/SF2, a prototypical SR protein. We found that SRPK1 was responsible for efficient phosphorylation of a short stretch of amino acids in the N-terminal portion of the RS domain of ASF/SF2 while Clk/Sty was able to transfer phosphate to all available serine residues in the RS domain, indicating that SR proteins may be phosphorylated by different kinases in a stepwise manner. Both kinases bind with high affinity and use fully processive catalytic mechanisms to achieve either restrictive or complete RS domain phosphorylation. These findings have important implications on the regulation of SR proteins in vivo by the SRPK and Clk/Sty families of kinases.     

Serine/arginine-rich splicing factors belong to a class of intrinsically disordered proteins

Serine/arginine-rich (SR) splicing factors play an important role in constitutive and alternative splicing as well as during several steps of RNA metabolism. Despite the wealth of functional information about SR proteins accumulated to-date, structural knowledge about the members of this family is very limited. To gain a better insight into structure-function relationships of SR proteins, we performed extensive sequence analysis of SR protein family members and combined it with ordered/disordered structure predictions. We found that SR proteins have properties characteristic of intrinsically disordered (ID) proteins. The amino acid composition and sequence complexity of SR proteins were very similar to those of the disordered protein regions. More detailed analysis showed that the SR proteins, and their RS domains in particular, are enriched in the disorder-promoting residues and are depleted in the order-promoting residues as compared to the entire human proteome. Moreover, disorder predictions indicated that RS domains of SR proteins were completely unstructured. Two different classification methods, the charge-hydropathy measure and the cumulative distribution function (CDF) of the disorder scores, were in agreement with each other, and they both strongly predicted members of the SR protein family to be disordered. This study emphasizes the importance of the disordered structure for several functions of SR proteins, such as for spliceosome assembly and for interaction with multiple partners. In addition, it demonstrates the usefulness of order/disorder predictions for inferring protein structure from sequence.     

PIRLs: a novel class of plant intracellular leucine-rich repeat proteins

Leucine-rich repeat (LRR) proteins feature tandem leucine-rich motifs that form a protein-protein interaction domain. Plants contain diverse classes of LRR proteins, many of which take part in signal transduction. We have identified a novel family of nine Arabidopsis LRR proteins that, based on predicted intracellular location and LRR motif consensus sequence, are related to Ras-binding LRR proteins found in signaling complexes in animals and yeast. This new class has been named plant intracellular Ras group-related LRR proteins (PIRLs). We have characterized PIRL cDNAs, rigorously defined gene and protein annotations, investigated gene family evolution and surveyed mRNA expression. While LRR regions suggested a relationship to Ras group LRR proteins, outside of their LRR domains PIRLs differed from Ras group proteins, exhibiting N- and C-terminal regions containing low complexity stretches and clusters of charged amino acids. PIRL genes grouped into three subfamilies based on sequence relationships and gene structures. Related gene pairs and dispersed chromosomal locations suggested family expansion by ancestral genomic or segmental duplications. Expression surveys revealed that all PIRL mRNAs are actively transcribed, with three expressed differentially in leaves, roots or flowers. These results define PIRLs as a distinct, plant-specific class of intracellular LRR proteins that probably mediate protein interactions, possibly in the context of signal transduction. T-DNA knock-out mutants have been isolated as a starting point for systematic functional analysis of this intriguing family.     

Axo-glial interactions regulate the localization of axonal paranodal proteins

The SR proteins, a group of abundant arginine/serine (RS)-rich proteins, are essential pre-mRNA splicing factors that are localized in the nucleus. The RS domain of these proteins serves as a nuclear localization signal. We found that RS domain-bearing proteins do not utilize any of the known nuclear import receptors and identified a novel nuclear import receptor specific for SR proteins. The SR protein import receptor, termed transportin-SR (TRN-SR), binds specifically and directly to the RS domains of ASF/SF2 and SC35 as well as several other SR proteins. The nuclear transport regulator RanGTP abolishes this interaction. Recombinant TRN-SR mediates nuclear import of RS domain- bearing proteins in vitro. TRN-SR has amino acid sequence similarity to several members of the importin beta/transportin family. These findings strongly suggest that TRN-SR is a nuclear import receptor for the SR protein family.     

Molecular modeling of pathogenesis-related proteins of family 5

The family of pathogenesis-related (PR) 5 proteins have diverse functions, and some of them are classified as thaumatins, osmotins, and inhibitors of α-amylase or trypsin. Although the specific function of many PR5 in plants is unknown, they are involved in the acquired systemic resistance and response to biotic stress, causing the inhibition of hyphal growth and reduction of spore germination, probably by a membrane permeabilization mechanism or by interaction with pathogen receptors. We have constructed three-dimensional models of four proteins belonging to the Rosaceae and Fagaceae botanical families by using the technique of comparative molecular modelling by homology. There are four main structural differences between all the PR5, corresponding to regions with replacements of amino acids. Folding and the secondary structures are very similar for all of them. However, the isoelectric point and charge distributions differ for earch protein.     

The novel flightless-I gene brings together two gene families, actin- binding proteins related to gelsolin and leucine-rich-repeat proteins involved in Ras signal transduction

The Drosophila melanogaster gene flightless-I, involved in gastrulation and muscle degeneration, has Caenorhabditis elegans and human homologues. In these highly conserved genes, two previously known gene families have been brought together, families encoding the actin- binding proteins related to gelsolin and the leucine-rich-repeat (LRR) group of proteins involved in protein-protein interactions. Both these gene families exhibit characteristics of molecular changes involving replication slippage and exon shuffling. Phylogenetic analyses of 19 amino acid sequences of 6 related protein types indicate that actin- associated proteins related to gelsolin are monophyletic to a common ancestor and include flightless proteins. Conversely, comparison of 24 amino acid sequences of LRR proteins including the flightless proteins indicates that flightless proteins are members of a structurally related subgroup. Included in the flightless cluster are human and mouse rsp-1 proteins involved in suppressing v-Ras transformation of cells and the membrane-associated yeast (Saccharomyces cerevisae) adenylate cyclase whose analogous LRRs are required for interaction with Ras proteins. There is a strong possibility that ligands for this group could be related and that flightless may have a similar role in Ras signal transduction. It is hypothesized that an ancestral monomeric gelsolin precursor protein has undergone at least four independent gene reorganization events to account for the structural diversity of the extant family of gelsolin-related proteins and that gene duplication and exon shuffling events occurred prior to or at the beginning of multicellular life, resulting in the evolution of some members of the family soon after the appearance of actin-type proteins.      

Co-evolution of the branch site and SR proteins in eukaryotes

Serine-arginine-rich (SR) proteins are essential for splicing in metazoans but are absent in yeast. By contrast, many fungi have SR protein homologs with variable arginine-rich regions analogous to the arginine-serine-rich (RS) domain in metazoans. The density of RS repeats in these regions correlates with the conservation of the branch site signal, providing evidence for an ancestral origin of SR proteins and indicating that the SR proteins and the branch site co-evolved.     

A new method for identification of protein (sub)families in a set of proteins based on hydropathy distribution in proteins

Structural similarity among proteins is reflected in the distribution of hydropathicity along the amino acids in the protein sequence. Similarities in the hydropathy distributions are obvious for homologous proteins within a protein family. They also were observed for proteins with related structures, even when sequence similarities were undetectable. Here we present a novel method that employs the hydropathy distribution in proteins for identification of (sub)families in a set of (homologous) proteins. We represent proteins as points in a generalized hydropathy space, represented by vectors of specifically defined features. The features are derived from hydropathy of the individual amino acids. Projection of this space onto principal axes reveals groups of proteins with related hydropathy distributions. The groups identified correspond well to families of structurally and functionally related proteins. We found that this method accurately identifies protein families in a set of proteins, or subfamilies in a set of homologous proteins. Our results show that protein families can be identified by the analysis of hydropathy distribution, without the need for sequence alignment.     

U box proteins as a new family of ubiquitin-protein ligases.

The U box is a domain of approximately 70 amino acids that is present in proteins from yeast to humans. The prototype U box protein, yeast Ufd2, was identified as a ubiquitin chain assembly factor that cooperates with a ubiquitin-activating enzyme (E1), a ubiquitin-conjugating enzyme (E2), and a ubiquitin-protein ligase (E3) to catalyze ubiquitin chain formation on artificial substrates. E3 enzymes are thought to determine the substrate specificity of ubiquitination and have been classified into two families, the HECT and RING finger families. Six mammalian U box proteins have now been shown to mediate polyubiquitination in the presence of E1 and E2 and in the absence of E3. These U box proteins exhibited different specificities for E2 enzymes in this reaction. Deletion of the U box or mutation of conserved amino acids within it abolished ubiquitination activity. Some U box proteins catalyzed polyubiquitination by targeting lysine residues of ubiquitin other than lysine 48, which is utilized by HECT and RING finger E3 enzymes for polyubiquitination that serves as a signal for proteolysis by the 26 S proteasome. These data suggest that U box proteins constitute a third family of E3 enzymes and that E4 activity may reflect a specialized type of E3 activity.     

Circular dichroism studies on synthetic peptides corresponding to the cleavage site region of precursor proteins

The conformations of synthetic peptides which span the region in which the precursor part of proteins (signal sequences) destined for export are cleaved by signal peptidases, were investigated by circular dichroism spectroscopy. Pentapeptides comprising amino acids only from the carboxy-terminus of signal sequences or the amino terminus of the mature protein do not have any preferred conformation in a variety of solvents. Octa- and nonapeptides containing amino acids from the carboxy-terminal protion of signal sequences and the amino-terminus of the mature portions of precursor proteins tend to adopt beta-turn conformations in trifluoroethanol and micelles of sodium dodecylsulphate. Hence, in addition to the distribution of amino acids with small side chains at the carboxy terminus of signal sequences, it is conceivable that signal peptidases also recognize a beta-turn conformation in the cleavage site region of precursor proteins.     

Characterization of functional domains of human EB1 family proteins

EB1 family proteins are evolutionarily conserved proteins that bind microtubule plus-ends and centrosomes and regulate the dynamics and organization of microtubules. Human EB1 family proteins, which include EB1, EBF3, and RP1, also associate with the tumor suppressor protein adenomatous polyposis coli (APC) and p150glued, a component of the dynactin complex. The structural basis for interaction between human EB1 family proteins and their associated proteins has not been defined in detail. EB1 family proteins have a calponin homology (CH) domain at their N terminus and an EB1-like C-terminal motif at their C terminus; the functional importance of these domains has not been determined. To better understand functions of human EB1 family proteins and to reveal functional similarities and differences among these proteins, we performed detailed characterizations of interactions between human EB1 family proteins and their associated proteins. We show that amino acids 1-133 of EB1 and EBF3 and the corresponding region of RP1, which contain a CH domain, are necessary and sufficient for binding microtubules, thus demonstrating for the first time that a CH domain contributes to binding microtubules. EB1 family proteins use overlapping but different regions that contain the EB1-like C-terminal motif to associate with APC and p150glued. Neither APC nor p150glued binding domain is necessary for EB1 or EBF3 to induce microtubule bundling, which requires amino acids 1-181 and 1-185 of EB1 and EBF3, respectively. We also determined that the EB1 family protein-binding regions are amino acids 2781-2820 and 18-111 of APC and p150glued, respectively.     

Molecular identification and chromosomal localization of genes encoding<Emphasis Type="Italic"> Triticum aestivum</Emphasis> xylanase inhibitor I-like proteins in cereals

TAXI (Triticum aestivum xylanase inhibitor) proteins are present in wheat flour and are known to inhibit glycosyl hydrolase family 11 endoxylanases, enzymes which are commonly applied in grain processing. Here, we describe the PCR-based molecular identification of genes encoding endoxylanase inhibitors HVXI and SCXI, the TAXI-like proteins from barley (Hordeum vulgare) and rye (Secale cereale) respectively. The HVXI coding sequence encodes a mature protein of 384 amino acids preceded by a 19 amino acid long signal sequence. SCXI-II/III has an open reading frame encoding a signal peptide of 21 amino acids and a mature protein of 375 amino acids. As for TAXI-I, no introns were detected in the untranslated regions and coding sequences identified. These newly identified sequences allowed us to perform a multiple sequence alignment with TAXI-I and similar proteins. Rice TAXI-type proteins clustered together with the cereal endoxylanase inhibitors. Dicotyledonous proteins with sequence similarity to TAXI-I, including the tomato xyloglucan-specific endoglucanase inhibiting protein, formed a different clade. The TAXI-type proteins may hence be part of a superfamily of proteins all involved in plant responses to biotic or abiotic stress and for which a function as glycosyl hydrolase inhibitors can be suggested. The chromosomal localization of the TAXI-I gene identified on wheat chromosome 3B, of the SCXI-II/III gene identified on rye chromosome 6R, and the presence of a cluster of TAXI-like genes on rice chromosome 1, allowed us to assign the location of TAXI-like genes to the wheat-rye translocation area 3BL/6RL characterized by RFLP markers XGlb33 and Xpsr454 and isozyme Est-5. In rice, RFLP marker C1310S corresponds to a TAXI-like protein encoding sequence.Communicated by P. Langridge     

Disorder predictors also predict backbone dynamics for a family of disordered proteins

Several algorithms have been developed that use amino acid sequences to predict whether or not a protein or a region of a protein is disordered. These algorithms make accurate predictions for disordered regions that are 30 amino acids or longer, but it is unclear whether the predictions can be directly related to the backbone dynamics of individual amino acid residues. The nuclear Overhauser effect between the amide nitrogen and hydrogen (NHNOE) provides an unambiguous measure of backbone dynamics at single residue resolution and is an excellent tool for characterizing the dynamic behavior of disordered proteins. In this report, we show that the NHNOE values for several members of a family of disordered proteins are highly correlated with the output from three popular algorithms used to predict disordered regions from amino acid sequence. This is the first test between an experimental measure of residue specific backbone dynamics and disorder predictions. The results suggest that some disorder predictors can accurately estimate the backbone dynamics of individual amino acids in a long disordered region.     

Poly(ADP-ribose) binds to specific domains in DNA damage checkpoint proteins

Poly(ADP-ribose) is formed in possibly all multicellular organisms by a familiy of poly(ADP-ribose) polymerases (PARPs). PARP-1, the best understood and until recently the only known member of this family, is a DNA damage signal protein catalyzing its automodification with multiple, variably sized ADP-ribose polymers that may contain up to 200 residues and several branching points. Through these polymers, PARP-1 can interact noncovalently with other proteins and alter their functions. Here we report the discovery of a poly(ADP-ribose)-binding sequence motif in several important DNA damage checkpoint proteins. The 20-amino acid motif contains two conserved regions: (i) a cluster rich in basic amino acids and (ii) a pattern of hydrophobic amino acids interspersed with basic residues. Using a combination of alanine scanning, polymer blot analysis, and photoaffinity labeling, we have identified poly(ADP-ribose)-binding sites in the following proteins: p53, p21(CIP1/WAF1), xeroderma pigmentosum group A complementing protein, MSH6, DNA ligase III, XRCC1, DNA polymerase epsilon, DNA-PK(CS), Ku70, NF-kappaB, inducible nitric-oxide synthase, caspase-activated DNase, and telomerase. The poly(ADP-ribose)-binding motif was found to overlap with five important functional domains responsible for (i) protein-protein interactions, (ii) DNA binding, (iii) nuclear localization, (iv) nuclear export, and (v) protein degradation. Thus, PARPs may target specific signal network proteins via poly(ADP-ribose) and regulate their domain functions.     

Molecular characterization of cold adaptation based on ortholog protein sequences from <Emphasis Type="Italic">Vibrionaceae</Emphasis> species

A set of 298 protein families from psychrophilic Vibrio salmonicida was compiled to identify genotypic characteristics that discern it from orthologous sequences from the mesophilic Vibrio/Photobacterium branch of the gamma-Proteobacteria (Vibrionaceae family). In our comparative exploration we employed alignment based bioinformatical and statistical methods. Interesting information was found in the substitution matrices, and the pattern of asymmetries in the amino acid substitution process. Together with the compositional difference, they identified the amino acids Ile, Asn, Ala and Gln as those having the most psycrophilic involvement. Ile and Asn are enhanced whereas Gln and Ala are suppressed. The inflexible Pro residue is also suppressed in loop regions, as expected in a flexible structure. The dataset were also classified and analysed according to the predicted subcellular location, and we made an additional study of 183 intracellular and 65 membrane proteins. Our results revealed that the psychrophilic proteins have similar hydrophobic and charge contributions in the core of the protein as mesophilic proteins, while the solvent-exposed surface area is significantly more hydrophobic. In addition, the psychrophilic intracellular (but not the membrane) proteins are significantly more negatively charged at the surface. Our analysis supports the hypothesis of preference for more flexible amino acids at the molecular surface. Life in cold climate seems to be obtained through many minor structural modifications rather than certain amino acids substitutions.     

Proline-rich proteins and glycoproteins: expression of salivary gland multigene families

Our recent research interests have focused on a group of unusual proteins and glycoproteins high in proline content, or the so-called proline-rich proteins (PRPs). The PRPs are tissue-specific expressions of salivary gland multigene families. Normally PRPs are not detected or are present in very low amounts in rat, mouse and hamster salivary glands, but these unusual proteins are dramatically induced by treatment with the catecholamine isoproterenol. The structures and organizations of several PRP mRNAs and PRP genes have been determined. The amino acid sequences of all PRPs show 4 distinct regions, namely, a signal peptide, a transition region, a repeat region and a carboxyl-terminal region. Glycoproteins induced by isoproterenol treatment may be N-glycosylated or O-glycosylated. The N-glycosylated glycoprotein GP-158 from rat submandibular glands has a 12 amino acid glycopeptide which repeats possibly 49 times. Proline-rich proteins of the parotid glands of rats and mice are also greatly induced by dietary tannins. The apparent unique occurrence of PRPs in saliva suggests that one biological role is to neutralize the detrimental effects of dietary tannins and other polyphenols. The upstream regions of the mouse and hamster PRP genes contain cyclic AMP-regulated sequences as demonstrated by deletions and transient transfections. The PRP multigene family members of mouse are all located on chromosome 8.     

An evolutionarily ancient Oatp: insights into conserved functional domains of these proteins

Cellular uptake of organic solutes is mediated in large part by a gene family of membrane transporters called OATPs (SLC21A). To study the structural determinants and evolutionary development of the SLC21A family, we have cloned and functionally characterized a highly expressed evolutionarily primitive Oatp from the liver of the small skate, Raja erinacea. A full-length cDNA (2.3 kb) was obtained that encodes a protein of 689 amino acids. The characteristics of this novel skate Oatp, including tissue expression, subcellular localization, substrate selectivity, Na(+) dependence, and inhibitor selectivity were generally similar to liver-specific human OATP-C and rat Oatp4. However, sequence comparisons with other OATPs indicate that this skate Oatp shares only approximately 40-50% amino acid identity with the liver-specific OATPs/Oatps and with human OATP-F. Further computer analysis revealed that the highest amino acid identities reside in the first external (78%) and internal loops (75%) and transmembrane domains 2 (76%), 3 (62%), 4 (70%), and 11 (64%). We propose that the conserved regions of the SLC21A transporter family may be critical structural determinants of substrate specificity and function.     

Higher-order interhelical spatial interactions in membrane proteins

Higher-order interactions are important for protein folding and assembly. We introduce the concept of interhelical three-body interactions as derived from Delaunay triangulation and alpha shapes of protein structures. In addition to glycophorin A, where triplets are strongly correlated with protein stability, we found that tight interhelical triplet interactions exist extensively in other membrane proteins, where many types of triplets occur far more frequently than in soluble proteins. We developed a probabilistic model for estimating the value of membrane helical interaction triplet (MHIT) propensity. Because the number of known structures of membrane proteins is limited, we developed a bootstrap method for determining the 95% confidence intervals of estimated MHIT values. We identified triplets that have high propensity for interhelical interactions and are unique to membrane proteins, e.g. AGF, AGG, GLL, GFF and others. A significant fraction (32%) of triplet types contains triplets that may be involved in interhelical hydrogen bond interactions, suggesting the prevalent and important roles of H-bond in the assembly of TM helices. There are several well-defined spatial conformations for triplet interactions on helices with similar parallel or antiparallel orientations and with similar right-handed or left-handed crossing angles. Often, they contain small residues and correspond to the regions of the closest contact between helices. Sequence motifs such as GG4 and AG4 can be part of the three-body interactions that have similar conformations, which in turn can be part of a higher-order cooperative four residue spatial motif observed in helical pairs from different proteins. In many cases, spatial motifs such as serine zipper and polar clamp are part of triplet interactions. On the basis of the analysis of the archaeal rhodopsin family of proteins, tightly packed triplet interactions can be achieved with several different choices of amino acid residues.