Codon Usage by Transposable Elements and Their Host Genes in Five Species

We compared the codon usage of sequences of transposable elements (TEs) with that of host genes from the species Drosophila melanogaster, Arabidopsis thaliana, Caenorhabditis elegans, Saccharomyces cerevisiae, and Homo sapiens. Factorial correspondence analysis showed that, regardless of the base composition of the genome, the TEs differed from the genes of their host species by their AT-richness. In all species, the percentage of A + T on the third codon position of the TEs was higher than that on the first codon position and lower than that in the noncoding DNA of the genomes. This indicates that the codon choice is not simply the outcome of mutational bias but is also subject to selection constraints. A tendency toward higher A + T on the third position than on the first position was also found in the host genes of A. thaliana, C. elegans, and S. cerevisiae but not in those of D. melanogaster and H. sapiens. This strongly suggests that the AT choice is a host-independent characteristic common to all TEs. The codon usage of TEs generally appeared to be different from the mean of the host genes. In the AT-rich genomes of Arabidopsis thaliana, Caenorhabditis elegans, and Saccharomyces cerevisiae, the codon usage bias of TEs was similar to that of weakly expressed genes. In the GC-rich genome of D. melanogaster, however, the bias in codon usage of the TEs clearly differed from that of weakly expressed genes. These findings suggest that selection acts on TEs and that TEs may display specific behavior within the host genomes. Received: 2 May 2001 / Accepted: 29 October 2001     

Polymerase chain reactions using Saccharomyces, Drosophila and human DNA predict a large family of protein serine/threonine phosphatases.

Using the polymerase chain reaction (PCR) to examine the protein serine/threonine phosphatase (PP) family which includes PP1, PP2A and PP2B, we have identified two, seven, and four novel protein phosphatase genes in Saccharomyces cerevisiae, Drosophila melanogaster and Homo sapiens, respectively. Consequently, the genes in the PP1/PP2A/PP2B family now number 11, 15 and 12 in these species respectively, and the data predicts still more unidentified phosphatases in higher eukaryotes. The PCR analyses also point to the presence in Drosophila and mammals of three or more different genes encoding PP2B, the enzyme recently identified as the target of certain immunosuppressant drugs.     

Protein phosphatase 6 subunit with conserved Sit4-associated protein domain targets IkappaBepsilon

Protein Ser/Thr phosphatases compose a PPP family that includes type-2 PP2A, PP4, and PP6, each with essential functions. The human PP6 gene rescues sit4(ts) mutants of Saccharomyces cerevisiae, and Sit4 phosphatase function depends on multiple Sit4-associated protein (SAP) subunits. We report here finding a SAPS sequence domain encoded in only a single gene each in Schizosaccharomyces pombe, Caenorhabditis elegans, and Drosophila but in three distinct open reading frames in Xenopus, Mus musculus, and Homo sapiens. The SAPS proteins are more divergent in sequence than PP6. Northern hybridization showed differential distribution of the human SAPS-related mRNA in multiple human tissues, named as PP6R1, PP6R2, and PP6R3. Antibodies were generated, distribution of endogenous PP6, PP6R1, PP6R2, and PP6R3 proteins was examined by immunoblotting, and the abundance of mRNA and protein in various tissues did not match. FLAG-tagged PP6R1 and PP6R2 expressed in HEK293 cells co-precipitated endogenous PP6, but not PP2A or PP4, showing specificity for recognition of phosphatases. The SAPS domain of PP6R1 alone was sufficient for association with PP6, and this predicts that conserved sequence motifs in the SAPS domain accounts for the specificity. FLAG-PP6R1 and FLAG-PP6R2 co-precipitated HA-IkappaBepsilon. Knockdown of PP6 or PP6R1 but not PP6R3 with siRNA significantly enhanced degradation of endogenous IkappaBepsilon in response to tumor necrosis factor-alpha. The results show SAPS domain subunits recruit substrates such as IkappaBepsilon as one way to determine specific functions for PP6.     

Probabilistic prediction of unknown metabolic and signal-transduction networks.

Regulatory networks provide control over complex cell behavior in all kingdoms of life. Here we describe a statistical model, based on representing proteins as collections of domains or motifs, which predicts unknown molecular interactions within these biological networks. Using known protein-protein interactions of Saccharomyces cerevisiae as training data, we were able to predict the links within this network with only 7% false-negative and 10% false-positive error rates. We also use Markov chain Monte Carlo simulation for the prediction of networks with maximum probability under our model. This model can be applied across species, where interaction data from one (or several) species can be used to infer interactions in another. In addition, the model is extensible and can be analogously applied to other molecular data (e.g., DNA sequences).     

Efficient sampling algorithm for estimating subgraph concentrations and detecting network motifs

SUMMARY: Biological and engineered networks have recently been shown to display network motifs: a small set of characteristic patterns that occur much more frequently than in randomized networks with the same degree sequence. Network motifs were demonstrated to play key information processing roles in biological regulation networks. Existing algorithms for detecting network motifs act by exhaustively enumerating all subgraphs with a given number of nodes in the network. The runtime of such algorithms increases strongly with network size. Here, we present a novel algorithm that allows estimation of subgraph concentrations and detection of network motifs at a runtime that is asymptotically independent of the network size. This algorithm is based on random sampling of subgraphs. Network motifs are detected with a surprisingly small number of samples in a wide variety of networks. Our method can be applied to estimate the concentrations of larger subgraphs in larger networks than was previously possible with exhaustive enumeration algorithms. We present results for high-order motifs in several biological networks and discuss their possible functions. AVAILABILITY: A software tool for estimating subgraph concentrations and detecting network motifs (mfinder 1.1) and further information is available at http://www.weizmann.ac.il/mcb/UriAlon/     

Distinctive topologies of partner-switching signaling networks correlate with their physiological roles

Regulatory networks controlling bacterial gene expression often evolve from common origins and share homologous proteins and similar network motifs. However, when functioning in different physiological contexts, these motifs may be re-arranged with different topologies that significantly affect network performance. Here we analyze two related signaling networks in the bacterium Bacillus subtilis in order to assess the consequences of their different topologies, with the aim of formulating design principles applicable to other systems. These two networks control the activities of the general stress response factor sigma(B) and the first sporulation-specific factor sigma(F). Both networks have at their core a "partner-switching" mechanism, in which an anti-sigma factor forms alternate complexes either with the sigma factor, holding it inactive, or with an anti-anti-sigma factor, thereby freeing sigma. However, clear differences in network structure are apparent: the anti-sigma factor for sigma(F) forms a long-lived, "dead-end" complex with its anti-anti-sigma factor and ADP, whereas the genes encoding sigma(B) and its network partners lie in a sigma(B)-controlled operon, resulting in positive and negative feedback loops. We constructed mathematical models of both networks and examined which features were critical for the performance of each design. The sigma(F) model predicts that the self-enhancing formation of the dead-end complex transforms the network into a largely irreversible hysteretic switch; the simulations reported here also demonstrate that hysteresis and slow turn off kinetics are the only two system properties associated with this complex formation. By contrast, the sigma(B) model predicts that the positive and negative feedback loops produce graded, reversible behavior with high regulatory capacity and fast response time. Our models demonstrate how alterations in network design result in different system properties that correlate with regulatory demands. These design principles agree with the known or suspected roles of similar networks in diverse bacteria.     

Dimerization specificity of all 67 B-ZIP motifs in Arabidopsis thaliana: a comparison to Homo sapiens B-ZIP motifs

Basic region-leucine zipper (B-ZIP) proteins are a class of dimeric sequence-specific DNA-binding proteins unique to eukaryotes. We have identified 67 B-ZIP proteins in the Arabidopsis thaliana genome. No A.thaliana B-ZIP domains are homologous with any Homo sapiens B-ZIP domains. Here, we predict the dimerization specificity properties of the 67 B-ZIP proteins in the A.thaliana genome based on three structural properties of the dimeric alpha-helical leucine zipper coiled coil structure: (i) length of the leucine zipper, (ii) placement of asparagine or a charged amino acid in the hydrophobic interface and (iii) presence of interhelical electrostatic interactions. Many A.thaliana B-ZIP leucine zippers are predicted to be eight or more heptads in length, in contrast to the four or five heptads typically found in H.sapiens, a prediction experimentally verified by circular dichroism analysis. Asparagine in the a position of the coiled coil is typically observed in the second heptad in H.sapiens. In A.thaliana, asparagine is abundant in the a position of both the second and fifth heptads. The particular placement of asparagine in the a position helps define 14 families of homodimerizing B-ZIP proteins in A.thaliana, in contrast to the six families found in H.sapiens. The repulsive interhelical electrostatic interactions that are used to specify heterodimerizing B-ZIP proteins in H.sapiens are not present in A.thaliana. Instead, we predict that plant leucine zippers rely on charged amino acids in the a position to drive heterodimerization. It appears that A.thaliana define many families of homodimerizing B-ZIP proteins by having long leucine zippers with asparagine judiciously placed in the a position of different heptads.     

Dynamic network rewiring determines temporal regulatory functions in Drosophila melanogaster development processes

The identification of network motifs has been widely considered as a significant step towards uncovering the design principles of biomolecular regulatory networks. To date, time‐invariant networks have been considered. However, such approaches cannot be used to reveal time‐specific biological traits due to the dynamic nature of biological systems, and hence may not be applicable to development, where temporal regulation of gene expression is an indispensable characteristic. We propose a concept of a “temporal sequence of network motifs”, a sequence of network motifs in active sub‐networks constructed over time, and investigate significant network motifs in the active temporal sub‐networks of Drosophila melanogaster . Based on this concept, we find a temporal sequence of network motifs which changes according to developmental stages and thereby cannot be identified from the whole static network. Moreover, we show that the temporal sequence of network motifs corresponding to each developmental stage can be used to describe pivotal developmental events.     

Phosphoinositide-dependent kinase-1 orthologues from five eukaryotes are activated by the hydrophobic motif in AGC kinases

Phosphoinositide-dependent kinase-1 (PDK1) mediates activation of many AGC kinases by docking onto a phosphorylated hydrophobic motif located C-terminal of the catalytic domain in the AGC kinase. The interaction shifts PDK1 into a conformation with increased catalytic activity and leads to autophosphorylation of PDK1. We demonstrate here that addition of a hydrophobic motif peptide increases the catalytic activity of PDK1 orthologues from Homo sapiens, Aplysia californica, Arabidopsis thaliana, Schizosaccharomyces pombe (ksg1), and Saccharomyces cerevisiae (Pkh1 and Pkh2) 2- to 12-fold. Furthermore, the hydrophobic motif peptide increases autophosphorylation of PDK1 from Homo sapiens, S. pombe, and S. cerevisiae (Phk2). Our results suggest that PDK1 interaction and activation by the hydrophobic motif of AGC kinases is a central mechanism in PDK1 function, which is conserved during eukaryotic evolution.     

NLSdb: database of nuclear localization signals

NLSdb is a database of nuclear localization signals (NLSs) and of nuclear proteins. NLSs are short stretches of residues mediating transport of nuclear proteins into the nucleus. The database contains 114 experimentally determined NLSs that were obtained through an extensive literature search. Using 'in silico mutagenesis' this set was extended to 308 experimental and potential NLSs. This final set matched over 43% of all known nuclear proteins and matches no currently known non-nuclear protein. NLSdb contains over 6000 predicted nuclear proteins and their targeting signals from the PDB and SWISS-PROT/TrEMBL databases. The database also contains over 12 500 predicted nuclear proteins from six entirely sequenced eukaryotic proteomes (Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Arabidopsis thaliana and Saccharomyces cerevisiae). NLS motifs often co-localize with DNA-binding regions. This observation was used to also annotate over 1500 DNA-binding proteins. NLSdb can be accessed via the web site: http://cubic.bioc.columbia.edu/db/NLSdb/.     

Phosphorylation and functions of inhibitor-2 family of proteins

Protein phosphatase-1 (PP1) is an essential protein Ser/Thr phosphatase that is extraordinarily conserved from yeast to human, and Inhibitor-2 (I-2) is the most ancient of the heat-stable proteins specific for PP1. We identified novel I-2 homologues in Caenorhabditis elegans (Ce) and Xenopus laevis (Xe) and compared them to the I-2 proteins from Homo sapiens (Hs), Saccharomyces cerevisiae (GLC8), and Drosophila melanogaster (Dm). The Ce I-2 and Dm I-2 showed the highest potency inhibition of rabbit PP1 with IC50 near 5 nM compared to Hs I-2 and Xe I-2 with IC50 between 10 and 50 nM and GLC8 with >100-fold lower activity. Inhibition of PP1 bound to Nek2 kinase activated the kinase to phosphorylate a C-Nap1 domain substrate. All the species of I-2 except GLC8 activated the Nek2::PP1 to the same extent as microcystin-LR. Only Hs I-2 and Xe I-2, not the I-2 proteins more divergent in sequence, directly activated human Aurora-A kinase. Various species of I-2 have a common PxTP phosphorylation site that showed a wide range of reactivity with GSK3, ERK, or CDC2/cyclinB1 kinases. The Suc1 subunit of CDC2/cyclinB1 enhanced reactivity with I-2, consistent with this being a site of mitotic phosphorylation. The results show species specificity among the I-2 family within the context of conserved PP1 inhibitory activity and variable phosphorylation by Pro-directed kinases.     

Propagation of Fluctuations in Biochemical Systems,I: Linear SSC Networks

We investigate the propagation of random fluctuations through biochemical networks in which the number of molecules of each species is large enough so that the concentrations are well modeled by differential equations. We study the effect of network topology on the emergent properties of the reaction system by characterizing the behavior of variance as fluctuations propagate down chains and studying the effect of side chains and feedback loops. We also investigate the asymptotic behavior of the system as one reaction becomes fast relative to the others.     

New insight into the evolution of aquaporins from flowering plants and vertebrates: orthologous identification and functional transfer is possible

Aquaporins (AQPs) represent a family of channel proteins that transport water and/or small solutes across cell membranes in the three domains of life. In all previous phylogenetic analysis of aquaporin, trees constructed using proteins with very low amino acid identity (<15%) were incongruent with rRNA data. In this work, restricting the evolutionary study of aquaporins to proteins with high amino acid identity (>25%), we showed congruence between AQPs and organismal trees. On the basis of this analysis, we defined 19 orthologous gene clusters in flowering plant species (3 PIP-like, 7 TIP-like, 6 NIP-like and 3 SIP-like). We described specific conserved motifs for each subfamily and each cluster, which were used to develop a method for automatic classification. Analysis of amino acid identity between orthologous monocotyledon and dicotyledon AQPs from each cluster, suggested that PIPs are under high evolutionary constraint. The phylogenetic analysis allowed us the assignment of orthologous aquaporins for very distant animal lineages (tetrapods-fishes). We also demonstrated that the location of all vertebrate AQPs in the ortholog clusters could be predicted by comparing their amino acid identity with human AQPs. We defined four AQP subfamilies in animals: AQP1-like, AQP8-like, AQP3-like and AQP11-like. Phylogenetic analysis showed that the four animal AQPs subfamilies are related with PIP-like, TIP-like, NIP-like and SIP-like subfamilies, respectively. Thus, this analysis would allow the prediction of individual AQPs function on the basis of orthologous genes from Arabidopsis thaliana and Homo sapiens.