Knowledge-based analysis of microarray gene expression data by using support vector machines

We introduce a method of functionally classifying genes by using gene expression data from DNA microarray hybridization experiments. The method is based on the theory of support vector machines (SVMs). SVMs are considered a supervised computer learning method because they exploit prior knowledge of gene function to identify unknown genes of similar function from expression data. SVMs avoid several problems associated with unsupervised clustering methods, such as hierarchical clustering and self-organizing maps. SVMs have many mathematical features that make them attractive for gene expression analysis, including their flexibility in choosing a similarity function, sparseness of solution when dealing with large data sets, the ability to handle large feature spaces, and the ability to identify outliers. We test several SVMs that use different similarity metrics, as well as some other supervised learning methods, and find that the SVMs best identify sets of genes with a common function using expression data. Finally, we use SVMs to predict functional roles for uncharacterized yeast ORFs based on their expression data.

Identification of thermophilic species by the amino acid compositions deduced from their genomes

The global amino acid compositions as deduced from the complete genomic sequences of six thermophilic archaea, two thermophilic bacteria, 17 mesophilic bacteria and two eukaryotic species were analysed by hierarchical clustering and principal components analysis. Both methods showed an influence of several factors on amino acid composition. Although GC content has a dominant effect, thermophilic species can be identified by their global amino acid compositions alone. This study presents a careful statistical analysis of factors that affect amino acid composition and also yielded specific features of the average amino acid composition of thermophilic species. Moreover, we introduce the first example of a ‘compositional tree’ of species that takes into account not only homologous proteins, but also proteins unique to particular species. We expect this simple yet novel approach to be a useful additional tool for the study of phylogeny at the genome level.

HLA class I and II antigens are partially co-clustered in the plasma membrane of human lymphoblastoid cells

Major histocompatibility complex (MHC) class II molecules displayed clustered patterns at the surfaces of T (HUT-102B2) and B (JY) lymphoma cells characterized by interreceptor distances in the micrometer range as detected by scanning force microscopy of immunogold-labeled antigens. Electron microscopy revealed that a fraction of the MHC class II molecules was also heteroclustered with MHC class I antigens at the same hierarchical level as described by the scanning force microscopy data, after specifically and sequentially labeling the antigens with 30- and 15-nm immunogold beads. On JY cells the estimated fraction of co-clustered HLA II was 0.61, whereas that of the HLA I was 0.24. Clusterization of the antigens was detected by the deviation of their spatial distribution from the Poissonian distribution representing the random case. Fluorescence resonance energy transfer measurements also confirmed partial co-clustering of the HLA class I and II molecules at another hierarchical level characterized by the 2- to 10-nm Förster distance range and providing fine details of the molecular organization of receptors. The larger-scale topological organization of the MHC class I and II antigens may reflect underlying membrane lipid domains and may fulfill significant functions in cell-to-cell contacts and signal transduction.

The use of site-directed fluorophore labeling and donor–donor energy migration to investigate solution structure and dynamics in proteins

The use of molecular genetics for introducing fluorescent molecules enables the use of donor–donor energy migration to determine intramolecular distances in a variety of proteins. This approach can be applied to examine the overall molecular dimensions of proteins and to investigate structural changes upon interactions with specific target molecules. In this report, the donor–donor energy migration method is demonstrated by experiments with the latent form of plasminogen activator inhibitor type 1. Based on the known x-ray structure of plasminogen activator inhibitor type 1, three positions forming the corners of a triangle were chosen. Double Cys substitution mutants (V106C-H185C, H185C-M266C, and M266C-V106C) and corresponding single substitution mutants (V106C, H185C, and M266C) were created and labeled with a sulfhydryl specific derivative of BODIPY (=the D molecule). The side lengths of this triangle were obtained from analyses of the experimental data. The analyses account for the local anisotropic order and rotational motions of the D molecules, as well as for the influence of a partial DD-labeling. The distances, as determined from x-ray diffraction, between the Cα-atoms of the positions V106C–H185C, H185C–M266C, and M266C–V106C were 60.9, 30.8, and 55.1 Å, respectively. These are in good agreement with the distances of 54 ± 4, 38 ± 3, and 55 ± 3 Å, as determined between the BODIPY groups attached via linkers to the same residues. Although the positions of the D-molecules and the Cα-atoms physically cannot coincide, there is a reasonable agreement between the methods.

Postsynaptic clustering of γ-aminobutyric acid type A receptors by the γ3 subunit in vivo

Synaptic localization of γ-aminobutyric acid type A (GABAA) receptors is a prerequisite for synaptic inhibitory function, but the mechanism by which different receptor subtypes are localized to postsynaptic sites is poorly understood. The γ2 subunit and the postsynaptic clustering protein gephyrin are required for synaptic localization and function of major GABAA receptor subtypes. We now show that transgenic overexpression of the γ3 subunit in γ2 subunit-deficient mice restores benzodiazepine binding sites, benzodiazepine-modulated whole cell currents, and postsynaptic miniature currents, suggesting the formation of functional, postsynaptic receptors. Moreover, the γ3 subunit can substitute for γ2 in the formation of GABAA receptors that are synaptically clustered and colocalized with gephyrin in vivo. These clusters were formed even in brain regions devoid of endogenous γ3 subunit, indicating that the factors present for clustering of γ2 subunit-containing receptors are sufficient to cluster γ3 subunit-containing receptors. The GABAA receptor and gephyrin-clustering properties of the ectopic γ3 subunit were also observed for the endogenous γ3 subunit, but only in the absence of the γ2 subunit, suggesting that the γ3 subunit is at a competitive disadvantage with the γ2 subunit for clustering of postsynaptic GABAA receptors in wild-type mice.

Activity-dependent regulation of synaptic clustering in a hippocampal culture system

Currently, there is a limited understanding of the factors that influence the localization and density of individual synapses in the central nervous system. Here we have studied the effects of activity on synapse formation between hippocampal dentate granule cells and CA3 pyramidal neurons in culture, taking advantage of FM1–43 as a fluorescent marker of synaptic boutons. We observed an early tendency for synapses to group together, quickly followed by the appearance of synaptic clusters on dendritic processes. These events were strongly influenced by N-methyl-d-aspartic acid receptor- and cyclic AMP-dependent signaling. The microstructure and localization of the synaptic clusters resembled that found in hippocampus, at mossy fiber synapses of stratum lucidum. Activity-dependent clustering of synapses represents a means for synaptic targeting that might contribute to synaptic organization in the brain.

Kinase domain of the muscle-specific receptor tyrosine kinase (MuSK) is sufficient for phosphorylation but not clustering of acetylcholine receptors: Required role for the MuSK ectodomain?

Formation of the neuromuscular junction (NMJ) depends upon a nerve-derived protein, agrin, acting by means of a muscle-specific receptor tyrosine kinase, MuSK, as well as a required accessory receptor protein known as MASC. We report that MuSK does not merely play a structural role by demonstrating that MuSK kinase activity is required for inducing acetylcholine receptor (AChR) clustering. We also show that MuSK is necessary, and that MuSK kinase domain activation is sufficient, to mediate a key early event in NMJ formation—phosphorylation of the AChR. However, MuSK kinase domain activation and the resulting AChR phosphorylation are not sufficient for AChR clustering; thus we show that the MuSK ectodomain is also required. These results indicate that AChR phosphorylation is not the sole trigger of the clustering process. Moreover, our results suggest that, unlike the ectodomain of all other receptor tyrosine kinases, the MuSK ectodomain plays a required role in addition to simply mediating ligand binding and receptor dimerization, perhaps by helping to recruit NMJ components to a MuSK-based scaffold.

Implications of a possible clustering of highest-energy cosmic rays

Recently, a possible clustering of a subset of observed ultra-high energy cosmic rays above ≃40 EeV (4 × 1019 eV) in pairs near the supergalactic plane was reported. We show that a confirmation of this effect would provide information on the origin and nature of these events and, in case of charged primaries, imply interesting constraints on the extragalactic magnetic field. Possible implications for the most common models of ultra-high energy cosmic ray production in the literature are discussed.

Involvement of unique leucine-zipper motif of PSD-Zip45 (Homer 1c/vesl-1L) in group 1 metabotropic glutamate receptor clustering

Several scaffold proteins for neurotransmitter receptors have been identified as candidates for receptor targeting. However, the molecular mechanism underlying such receptor clustering and targeting to postsynaptic specializations remains unknown. PSD-Zip45 (also named Homer 1c/vesl-1L) consists of the NH2 terminus containing the enabled/VASP homology 1 domain and the COOH terminus containing the leucine zipper. Here, we demonstrate immunohistochemically that metabotropic glutamate receptor 1α (mGluR1α) and PSD-Zip45/Homer 1c are colocalized to synapses in the cerebellar molecular layer but not in the hippocampus. In cultured hippocampal neurons, PSD-Zip45/Homer1c and N-methyl-d-aspartate receptors are preferentially colocalized to dendritic spines. Cotransfection of mGluR1α or mGluR5 and PSD-Zip45/Homer 1c into COS-7 cells results in mGluR clustering induced by PSD-Zip45/Homer 1c. An in vitro multimerization assay shows that the extreme COOH-terminal leucine zipper is involved in self-multimerization of PSD-Zip45/Homer 1c. A clustering assay of mGluRs in COS-7 cells also reveals a critical role of this leucine-zipper motif of PSD-Zip45/Homer 1c in mGluR clustering. These results suggest that the leucine zipper of subsynaptic scaffold protein is a candidate motif involved in neurotransmitter receptor clustering at the central synapse.

A hierarchical neuronal network for planning behavior

Planning a goal-directed sequence of behavior is a higher function of the human brain that relies on the integrity of prefrontal cortical areas. In the Tower of London test, a puzzle in which beads sliding on pegs must be moved to match a designated goal configuration, patients with lesioned prefrontal cortex show deficits in planning a goal-directed sequence of moves. We propose a neuronal network model of sequence planning that passes this test and, when lesioned, fails in a way that mimics prefrontal patients’ behavior. Our model comprises a descending planning system with hierarchically organized plan, operation, and gesture levels, and an ascending evaluative system that analyzes the problem and computes internal reward signals that index the correct/erroneous status of the plan. Multiple parallel pathways connecting the evaluative and planning systems amend the plan and adapt it to the current problem. The model illustrates how specialized hierarchically organized neuronal assemblies may collectively emulate central executive or supervisory functions of the human brain.

A transmembrane form of annexin XII detected by site-directed spin labeling

Previous studies of the annexin family of Ca2+ binding proteins identified a soluble monomer in the absence of Ca2+ and a trimer adsorbed on the membrane surface in the presence of Ca2+. On the basis of site-directed spin-labeling studies of annexin XII at low pH, we now report a membrane-inserted form of the protein with a dramatically different structure. The data suggest that upon insertion a continuous transmembrane α-helix is reversibly formed from a helix–loop–helix motif in the solution structure. Other regions with similar membrane-insertion potential were identified in the amino acid sequence, and we propose that the corresponding helices come together to form an aqueous pore that mediates the ion channel activity reported for several annexins.

Clustering of meiotic double-strand breaks on yeast chromosome III

In the yeast Saccharomyces cerevisiae, meiotic recombination is initiated by transient DNA double-strand breaks (DSBs) that are repaired by interaction of the broken chromosome with its homologue. To identify a large number of DSB sites and gain insight into the control of DSB formation at both the local and the whole chromosomal levels, we have determined at high resolution the distribution of meiotic DSBs along the 340 kb of chromosome III. We have found 76 DSB regions, mostly located in intergenic promoter-containing intervals. The frequency of DSBs varies at least 50-fold from one region to another. The global distribution of DSB regions along chromosome III is nonrandom, defining large (39–105 kb) chromosomal domains, both hot and cold. The distribution of these localized DSBs indicates that they are likely to initiate most crossovers along chromosome III, but some discrepancies remain to be explained.

Calcium-dependent Clustering of Inositol 1,4,5-Trisphosphate Receptors

Rat basophilic leukemia (RBL-2H3) cells predominantly express the type II receptor for inositol 1,4,5-trisphosphate (InsP3), which operates as an InsP3-gated calcium channel. In these cells, cross-linking the high-affinity immunoglobulin E receptor (FcεR1) leads to activation of phospholipase C γ isoforms via tyrosine kinase- and phosphatidylinositol 3-kinase-dependent pathways, release of InsP3-sensitive intracellular Ca2+ stores, and a sustained phase of Ca2+ influx. These events are accompanied by a redistribution of type II InsP3 receptors within the endoplasmic reticulum and nuclear envelope, from a diffuse pattern with a few small aggregates in resting cells to large isolated clusters after antigen stimulation. Redistribution of type II InsP3 receptors is also seen after treatment of RBL-2H3 cells with ionomycin or thapsigargin. InsP3 receptor clustering occurs within 5–10 min of stimulus and persists for up to 1 h in the presence of antigen. Receptor clustering is independent of endoplasmic reticulum vesiculation, which occurs only at ionomycin concentrations >1 μM, and maximal clustering responses are dependent on the presence of extracellular calcium. InsP3 receptor aggregation may be a characteristic cellular response to Ca2+-mobilizing ligands, because similar results are seen after activation of phospholipase C-linked G-protein-coupled receptors; cholecystokinin causes type II receptor redistribution in rat pancreatoma AR4–2J cells, and carbachol causes type III receptor redistribution in muscarinic receptor-expressing hamster lung fibroblast E36M3R cells. Stimulation of these three cell types leads to a reduction in InsP3 receptor levels only in AR4–2J cells, indicating that receptor clustering does not correlate with receptor down-regulation. The calcium-dependent aggregation of InsP3 receptors may contribute to the previously observed changes in affinity for InsP3 in the presence of elevated Ca2+ and/or may establish discrete regions within refilled stores with varying capacity to release Ca2+ when a subsequent stimulus results in production of InsP3.

Direct Involvement of N-Cadherin–mediated Signaling in Muscle Differentiation

Cell–cell interactions, mediated by members of the cadherin family of Ca2+-dependent adhesion molecules, play key roles in morphogenetic processes as well as in the transduction of long-range growth and differentiation signals. In muscle differentiation cell adhesion is involved in both early stages of myogenic induction and in later stages of myoblast interaction and fusion. In this study we have explored the involvement of a specific cadherin, namely N-cadherin, in myogenic differentiation. For that purpose we have treated different established lines of cultured myoblasts with beads coated with N-cadherin–specific ligands, including a recombinant N-cadherin extracellular domain, and anti-N-cadherin antibodies. Immunofluorescent labeling for cadherins and catenins indicated that treatment with the cadherin-reactive beads for several hours enhances the assembly of cell–cell adherens-type junctions. Moreover, immunofluorescence and immunoblotting analyses indicated that treatment with the beads for 12–24 h induces myogenin expression and growth arrest, which are largely independent of cell plating density. Upon longer incubation with the beads (2–3 d) a major facilitation in the expression of several muscle-specific sarcomeric proteins and in cell fusion into myotubes was observed. These results suggest that surface clustering or immobilization of N-cadherin can directly trigger signaling events, which promote the activation of a myogenic differentiation program.