Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Lysophosphatidic Acid Enhances Pulmonary Epithelial Barrier Integrity and Protects Endotoxin-induced Epithelial Barrier Disruption and Lung Injury

Lysophosphatidic acid (LPA), a bioactive phospholipid, induces a wide range of cellular effects, including gene expression, cytoskeletal rearrangement, and cell survival. We have previously shown that LPA stimulates secretion of pro- and anti-inflammatory cytokines in bronchial epithelial cells. This study provides evidence that LPA enhances pulmonary epithelial barrier integrity through protein kinase C (PKC) δ- and ζ-mediated E-cadherin accumulation at cell-cell junctions. Treatment of human bronchial epithelial cells (HBEpCs) with LPA increased transepithelial electrical resistance (TER) by ∼2.0-fold and enhanced accumulation of E-cadherin to the cell-cell junctions through Gα_i-coupled LPA receptors. Knockdown of E-cadherin with E-cadherin small interfering RNA or pretreatment with EGTA (0.1 mm) prior to LPA (1 μm) treatment attenuated LPA-induced increases in TER in HBEpCs. Furthermore, LPA induced tyrosine phosphorylation of focal adhesion kinase (FAK) and overexpression of the FAK inhibitor, and FAK-related non-kinase-attenuated LPA induced increases in TER and E-cadherin accumulation at cell-cell junctions. Overexpression of dominant negative protein kinase δ and ζ attenuated LPA-induced phosphorylation of FAK, accumulation of E-cadherin at cell-cell junctions, and an increase in TER. Additionally, lipopolysaccharide decreased TER and induced E-cadherin relocalization from cell-cell junctions to cytoplasm in a dose-dependent fashion, which was restored by LPA post-treatment in HBEpCs. Intratracheal post-treatment with LPA (5 μm) reduced LPS-induced neutrophil influx, protein leak, and E-cadherin shedding in bronchoalveolar lavage fluids in a murine model of acute lung injury. These data suggest a protective role of LPA in airway inflammation and remodeling.The airway epithelium is the site of first contact for inhaled environmental stimuli, functions as a physical barrier to environmental insult, and is an essential part of innate immunity. Epithelial barrier disruption is caused by inhaled allergens, dust, and irritants, resulting in inflammation, bronchoconstriction, and edema as seen in asthma and other respiratory diseases (1–4). Furthermore, increased epithelial permeability also results in para-cellular leakage of large proteins, such as albumin, immunoglobulin G, and polymeric immunoglobulin A, into the airway lumen (5, 6). The epithelial cell-cell junctional complex is composed of tight junctions, adherens junctions, and desmosomes. These adherens junctions play a pivotal role in regulating the activity of the entire junctional complex because the formation of adherens junctions subsequently leads to the formation of other cell-cell junctions (7–9). The major adhesion molecules in the adherens junctions are the cadherins. E-cadherin is a member of the cadherin family that mediates calcium-dependent cell-cell adhesion. The N-terminal ectodomain of E-cadherin contains homophilic interaction specificity, and the cytoplasmic domain binds to catenins, which interact with actin (10–13). Plasma membrane localization of E-cadherin is critical for the maintenance of epithelial cell-cell junctions and airway epithelium integrity (7, 10, 14). A decrease of adhesive properties of E-cadherin is related to the loss of differentiation and the subsequent acquisition of a higher motility and invasiveness of epithelial cells (10, 14, 15). Dislocation or shedding of E-cadherin in the airway epithelium induces epithelial shedding and increases airway permeability in lung airway diseases (10, 14, 16). In an ovalbumin-challenged guinea pig model of asthma, it has been demonstrated that E-cadherin is dislocated from the lateral margins of epithelial cells (10). Histamine increases airway para-cellular permeability and results in an increased susceptibility of airway epithelial cells to adenovirus infection by interrupting E-cadherin adhesion (14). Serine phosphorylation of E-cadherin by casein kinase II, GSK-3β, and PKD1/PKC² μ enhanced E-cadherin-mediated cell-cell adhesion in NIH3T3 fibroblasts and LNCaP prostate cancer cells (11, 17). However, the regulation and mechanism by which E-cadherin is localized within the pulmonary epithelium is not fully known, particularly during airway remodeling.LPA, a naturally occurring bioactive lipid, is present in body fluids, such as plasma, saliva, follicular fluid, malignant effusions, and bronchoalveolar lavage (BAL) fluids (18–20). Six distinct high affinity cell-surface LPA receptors, LPA-R_1–6, have been cloned and described in mammals (21–26). Extracellular activities of LPA include cell proliferation, motility, and cell survival (27–30). LPA exhibits a wide range of effects on differing cell types, including pulmonary epithelial, smooth muscle, fibroblasts, and T cells (31–35). LPA augments migration and cytokine synthesis in lymphocytes and induces chemotaxis of Jurkat T cells through Matrigel membranes (34). LPA induces airway smooth muscle cell contractility, proliferation, and airway repair and remodeling (35, 36). LPA also potently stimulates IL-8 (31, 37–39), IL-13 receptor α2 (IL-13Rα2) (40), and COX-2 gene expression and prostaglandin E₂ release (41) in HBEpCs. Prostaglandin E₂ and IL-13Rα2 have anti-inflammatory properties in pulmonary inflammation (42, 43). These results suggest that LPA may play a protective role in lung disease by stimulating an innate immune response while simultaneously attenuating the adaptive immune response. Furthermore, intravenous injection with LPA attenuated bacterial endotoxin-induced plasma tumor necrosis factor-α production and myeloperoxidase activity in the lungs of mice (44), suggesting an anti-inflammatory role for LPA in a murine model of sepsis.We have reported that LPA induces E-cadherin/c-Met accumulation in cell-cell contacts and increases TER in HBEpCs (45). Here, for the first time, we report that LPA-induced increases in TER are dependent on PKCδ, PKCζ, and FAK-mediated E-cadherin accumulation at cell-cell junctions. Furthermore, we demonstrate that post-treatment of LPA rescues LPS-induced airway epithelial disruption in vitro and reduces E-cadherin shedding in a murine model of ALI. This study identifies the molecular mechanisms linking the LPA and LPA receptors to maintaining normal pulmonary epithelium barrier function, which is critical in developing novel therapies directed at ameliorating pulmonary diseases.     

Differential Regulation of Actin Polymerization and Structure by Yeast Formin Isoforms

The budding yeast formins, Bnr1 and Bni1, behave very differently with respect to their interactions with muscle actin. However, the mechanisms underlying these differences are unclear, and these formins do not interact with muscle actin in vivo. We use yeast wild type and mutant actins to further assess these differences between Bnr1 and Bni1. Low ionic strength G-buffer does not promote actin polymerization. However, Bnr1, but not Bni1, causes the polymerization of pyrene-labeled Mg-G-actin in G-buffer into single filaments based on fluorometric and EM observations. Polymerization by Bnr1 does not occur with Ca-G-actin. By cosedimentation, maximum filament formation occurs at a Bnr1:actin ratio of 1:2. The interaction of Bnr1 with pyrene-labeled S265C Mg-actin yields a pyrene excimer peak, from the cross-strand interaction of pyrene probes, which only occurs in the context of F-actin. In F-buffer, Bnr1 promotes much faster yeast actin polymerization than Bni1. It also bundles the F-actin in contrast to the low ionic strength situation where only single filaments form. Thus, the differences previously observed with muscle actin are not actin isoform-specific. The binding of both formins to F-actin saturate at an equimolar ratio, but only about 30% of each formin cosediments with F-actin. Finally, addition of Bnr1 but not Bni1 to pyrene-labeled wild type and S265C Mg-F actins enhanced the pyrene- and pyrene-excimer fluorescence, respectively, suggesting Bnr1 also alters F-actin structure. These differences may facilitate the ability of Bnr1 to form the actin cables needed for polarized delivery of nutrients and organelles to the growing yeast bud.Bni1 and Bnr1 are the two formin isoforms expressed in Saccharomyces cerevisiae (1, 2). These proteins, as other isoforms in the formin family, are large multidomain proteins (3, 4). Several regulatory domains, including one for binding the G-protein rho, are located at the N-terminal half of the protein (4–7). FH1, FH2, and Bud6 binding domains are located in the C-terminal half of the protein (8). The formin homology 1 (FH1)² domain contains several sequential poly-l-proline motifs, and it interacts with the profilin/actin complex to recruit actin monomers and regulate the insertion of actin monomers at the barbed end of actin (9–11). The fomin homology domain 2 (FH2) forms a donut-shaped homodimer, which wraps around actin dimers at the barbed end of actin filaments (12, 13). One important function of formin is to facilitate actin polymerization by stabilizing actin dimers or trimers under polymerization conditions and then to processively associate with the barbed end of the elongating filament to control actin filament elongation kinetics (13–18).A major unsolved protein in the study of formins is the elucidation of the individual functions of different isoforms and their regulation. In vivo, these two budding yeast formins have distinct cellular locations and dynamics (1, 2, 19, 20). Bni1 concentrates at the budding site before the daughter cell buds from the mother cell, moves along with the tip of the daughter cell, and then travels back to the neck between daughter and mother cells at the end of segregation. Bnr1 localizes only at the neck of the budding cell in a very short period of time after bud emergence. Although a key cellular function of these two formins in yeast is to promote actin cable formation (8, 18), the roles of the individual formins in different cellular process is unclear because deleting either individual formin gene has limited impact on cell growth and deleting both genes together is lethal (21).Although each of the two formins can nucleate actin filament formation in vitro, the manner in which they affect polymerization is distinctly isoform-specific. Most of this mechanistic work in vitro has used formin fragments containing the FH1 and FH2 domains. Bni1 alone processively caps the barbed end of actin filaments partially inhibiting polymerization at this end (14, 16, 18). The profilin-actin complex, recruited to the actin barbed end through its binding to Bni1 FH1 domain, possibly raises the local actin concentration and appears to allow this inhibition to be overcome, thereby, accelerating barbed end polymerization. It has also been shown that this complex modifies the kinetics of actin dynamics at the barbed end (9, 11, 18, 22). Moreover, Bni1 participation leads only to the formation of single filaments (8). In comparison, the Bnr1 FH1-FH2 domain facilitates actin polymerization much more efficiently than does Bni1. Moseley and Goode (8) showed Bnr1 accelerates actin polymerization up to 10 times better than does Bni and produces actin filament bundles when the Bnr1/actin molar ratio is above 1:2. Finally, the regulation of Bni1 and Bnr1 by formin binding is different. For example, Bud 6/Aip3, a yeast cell polarity factor, binds to Bni1, but not Bnr1, and also stimulates its activity in vitro.For their studies, Moseley and Goode (8) utilized mammalian skeletal muscle actin instead of the S. cerevisiae actin with which the yeast formins are designed to function. It is entirely possible that the differences observed with the two formins are influenced quantitatively or qualitatively by the nature of the actin used in the study. This possibility must be seriously considered because although yeast and muscle actins are 87% identical in sequence, they display marked differences in their polymerization behavior (23). Yeast actin nucleates filaments better than muscle actin (24, 25). It appears to form shorter and more flexible filaments than does muscle actin (26, 27). Finally, the disposition of the P_i released during the hydrolysis of ATP that occurs during polymerization is different. Yeast actin releases its P_i concomitant with hydrolysis of the bound ATP whereas muscle actin retains the P_i for a significant amount of time following nucleotide hydrolysis (28, 29). This difference is significant because ADP-P_i F-actin has been shown to be more stable than ADP F-actin (30). Another example of this isoform dependence is the interaction of yeast Arp2/3 with yeast versus muscle actins (31). Yeast Arp2/3 complex accelerates polymerization of muscle actin only in the presence of a nucleation protein factor such as WASP. However, with yeast actin, no such auxiliary protein is required. In light of these actin behavioral differences, to better understand the functional differences of these two formins in vivo, we have studied the behavior of Bni 1 and Bnr 1 with WT and mutant yeast actins, and we have also explored the molecular basis underlying the Bnr 1-induced formation of actin nuclei from G-actin.     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]     

Deletion of the Chloride Transporter Slc26a7 Causes Distal Renal Tubular Acidosis and Impairs Gastric Acid Secretion

SLC26A7 (human)/Slc26a7 (mouse) is a recently identified chloride-base exchanger and/or chloride transporter that is expressed on the basolateral membrane of acid-secreting cells in the renal outer medullary collecting duct (OMCD) and in gastric parietal cells. Here, we show that mice with genetic deletion of Slc26a7 expression develop distal renal tubular acidosis, as manifested by metabolic acidosis and alkaline urine pH. In the kidney, basolateral Cl⁻/HCO₃⁻ exchange activity in acid-secreting intercalated cells in the OMCD was significantly decreased in hypertonic medium (a normal milieu for the medulla) but was reduced only mildly in isotonic medium. Changing from a hypertonic to isotonic medium (relative hypotonicity) decreased the membrane abundance of Slc26a7 in kidney cells in vivo and in vitro. In the stomach, stimulated acid secretion was significantly impaired in isolated gastric mucosa and in the intact organ. We propose that SLC26A7 dysfunction should be investigated as a potential cause of unexplained distal renal tubular acidosis or decreased gastric acid secretion in humans.The collecting duct segment of the distal kidney nephron plays a major role in systemic acid base homeostasis by acid secretion and bicarbonate absorption. The acid secretion occurs via H⁺-ATPase and H-K-ATPase into the lumen and bicarbonate is absorbed via basolateral Cl⁻/HCO₃⁻ exchangers (1–4). The tubules, which are located within the outer medullary region of the kidney collecting duct (OMCD),² have the highest rate of acid secretion among the distal tubule segments and are therefore essential to the maintenance of acid base balance (2).The gastric parietal cell is the site of generation of acid and bicarbonate through the action of cytosolic carbonic anhydrase II (5, 6). The intracellular acid is secreted into the lumen via gastric H-K-ATPase, which works in conjunction with a chloride channel and a K⁺ recycling pathway (7–10). The intracellular bicarbonate is transported to the blood via basolateral Cl⁻/HCO₃⁻ exchangers (11–14).SLC26 (human)/Slc26 (mouse) isoforms are members of a conserved family of anion transporters that display tissue-specific patterns of expression in epithelial cells (15–24). Several SLC26 members can function as chloride/bicarbonate exchangers. These include SLC26A3 (DRA), SLC26A4 (pendrin), SLC26A6 (PAT1 or CFEX), SLC26A7, and SLC26A9 (25–31). SLC26A7 and SLC26A9 can also function as chloride channels (32–34).SLC26A7/Slc26a7 is predominantly expressed in the kidney and stomach (28, 29). In the kidney, Slc26a7 co-localizes with AE1, a well-known Cl⁻/HCO₃⁻ exchanger, on the basolateral membrane of (acid-secreting) A-intercalated cells in OMCD cells (29, 35, 36) (supplemental Fig. 1). In the stomach, Slc26a7 co-localizes with AE2, a major Cl⁻/HCO₃⁻ exchanger, on the basolateral membrane of acid secreting parietal cells (28). To address the physiological function of Slc26a7 in the intact mouse, we have generated Slc26a7 ko mice. We report here that Slc26a7 ko mice exhibit distal renal tubular acidosis and impaired gastric acidification in the absence of morphological abnormalities in kidney or stomach.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Phosphoproteome Analysis of Functional Mitochondria Isolated from Resting Human Muscle Reveals Extensive Phosphorylation of Inner Membrane Protein Complexes and Enzymes

Mitochondria play a central role in energy metabolism and cellular survival, and consequently mitochondrial dysfunction is associated with a number of human pathologies. Reversible protein phosphorylation emerges as a central mechanism in the regulation of several mitochondrial processes. In skeletal muscle, mitochondrial dysfunction is linked to insulin resistance in humans with obesity and type 2 diabetes. We performed a phosphoproteomics study of functional mitochondria isolated from human muscle biopsies with the aim to obtain a comprehensive overview of mitochondrial phosphoproteins. Combining an efficient mitochondrial isolation protocol with several different phosphopeptide enrichment techniques and LC-MS/MS, we identified 155 distinct phosphorylation sites in 77 mitochondrial phosphoproteins, including 116 phosphoserine, 23 phosphothreonine, and 16 phosphotyrosine residues. The relatively high number of phosphotyrosine residues suggests an important role for tyrosine phosphorylation in mitochondrial signaling. Many of the mitochondrial phosphoproteins are involved in oxidative phosphorylation, tricarboxylic acid cycle, and lipid metabolism, i.e. processes proposed to be involved in insulin resistance. We also assigned phosphorylation sites in mitochondrial proteins involved in amino acid degradation, importers and transporters, calcium homeostasis, and apoptosis. Bioinformatics analysis of kinase motifs revealed that many of these mitochondrial phosphoproteins are substrates for protein kinase A, protein kinase C, casein kinase II, and DNA-dependent protein kinase. Our results demonstrate the feasibility of performing phosphoproteome analysis of organelles isolated from human tissue and provide novel targets for functional studies of reversible phosphorylation in mitochondria. Future comparative phosphoproteome analysis of mitochondria from healthy and diseased individuals will provide insights into the role of abnormal phosphorylation in pathologies, such as type 2 diabetes.Mitochondria are the primary energy-generating systems in eukaryotes. They play a crucial role in oxidative metabolism, including carbohydrate metabolism, fatty acid oxidation, and urea cycle, as well as in calcium signaling and apoptosis (1, 2). Mitochondrial dysfunction is centrally involved in a number of human pathologies, such as type 2 diabetes, Parkinson disease, and cancer (3). The most prevalent form of cellular protein post-translational modifications (PTMs),1 reversible phosphorylation (4–6), is emerging as a central mechanism in the regulation of mitochondrial functions (7, 8). The steadily increasing numbers of reported mitochondrial kinases, phosphatases, and phosphoproteins imply an important role of protein phosphorylation in different mitochondrial processes (9–11).Mass spectrometry (MS)-based proteome analysis is a powerful tool for global profiling of proteins and their PTMs, including protein phosphorylation (12, 13). A variety of proteomics techniques have been developed for specific enrichment of phosphorylated proteins and peptides and for phosphopeptide-specific data acquisition techniques at the MS level (14). Enrichment methods based on affinity chromatography, such as titanium dioxide (TiO₂) (15–17), zwitterionic hydrophilic interaction chromatography (ZIC-HILIC) (18), immobilized metal affinity chromatography (IMAC) (19, 20), and ion exchange chromatography (strong anion exchange and strong cation exchange) (21, 22), have shown high efficiencies for enrichment of phosphopeptides (14). Recently, we demonstrated that calcium phosphate precipitation (CPP) is highly effective for enriching phosphopeptides (23). It is now generally accepted that no single method is comprehensive, but combinations of different enrichment methods produce distinct overlapping phosphopeptide data sets to enhance the overall results in phosphoproteome analysis (24, 25). Phosphopeptide sequencing by mass spectrometry has seen tremendous advances during the last decade (26). For example, MS/MS product ion scanning, multistage activation, and precursor ion scanning are effective methods for identifying serine (Ser), threonine (Thr), and tyrosine (Tyr) phosphorylated peptides (14, 26).A “complete” mammalian mitochondrial proteome was reported by Mootha and co-workers (27) and included 1098 proteins. The mitochondrial phosphoproteome has been characterized in a series of studies, including yeast, mouse and rat liver, porcine heart, and plants (19, 28–31). To date, the largest data set by Deng et al. (30) identified 228 different phosphoproteins and 447 phosphorylation sites in rat liver mitochondria. However, the in vivo phosphoproteome of human mitochondria has not been determined. A comprehensive mitochondrial phosphoproteome is warranted for further elucidation of the largely unknown mechanisms by which protein phosphorylation modulates diverse mitochondrial functions.The percutaneous muscle biopsy technique is an important tool in the diagnosis and management of human muscle disorders and has been widely used to investigate metabolism and various cellular and molecular processes in normal and abnormal human muscle, in particular the molecular mechanism underlying insulin resistance in obesity and type 2 diabetes (32). Skeletal muscle is rich in mitochondria and hence a good source for a comprehensive proteomics and functional analysis of mitochondria (32, 33).The major aim of the present study was to obtain a comprehensive overview of site-specific phosphorylation of mitochondrial proteins in functionally intact mitochondria isolated from human skeletal muscle. Combining an efficient protocol for isolation of skeletal muscle mitochondria with several different state-of-the-art phosphopeptide enrichment methods and high performance LC-MS/MS, we identified 155 distinct phosphorylation sites in 77 mitochondrial phosphoproteins, many of which have not been reported before. We characterized this mitochondrial phosphoproteome by using bioinformatics tools to classify functional groups and functions, including kinase substrate motifs.