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.

Congruent Docking of Dimeric Kinesin and ncd into Three-dimensional Electron Cryomicroscopy Maps of Microtubule–Motor ADP Complexes

We present a new map showing dimeric kinesin bound to microtubules in the presence of ADP that was obtained by electron cryomicroscopy and image reconstruction. The directly bound monomer (first head) shows a different conformation from one in the more tightly bound empty state. This change in the first head is amplified as a movement of the second (tethered) head, which tilts upward. The atomic coordinates of kinesin·ADP dock into our map so that the tethered head associates with the bound head as in the kinesin dimer structure seen by x-ray crystallography. The new docking orientation avoids problems associated with previous predictions; it puts residues implicated by proteolysis-protection and mutagenesis studies near the microtubule but does not lead to steric interference between the coiled-coil tail and the microtubule surface. The observed conformational changes in the tightly bound states would probably bring some important residues closer to tubulin. As expected from the homology with kinesin, the atomic coordinates of nonclaret disjunctional protein (ncd)·ADP dock in the same orientation into the attached head in a map of microtubules decorated with dimeric ncd·ADP. Our results support the idea that the observed direct interaction between the two heads is important at some stages of the mechanism by which kinesin moves processively along microtubules.

Controlling cell attachment on contoured surfaces with self-assembled monolayers of alkanethiolates on gold.

This paper describes a method based on experimentally simple techniques--microcontact printing and micromolding in capillaries--to prepare tissue culture substrates in which both the topology and molecular structure of the interface can be controlled. The method combines optically transparent contoured surfaces with self-assembled monolayers (SAMs) of alkanethiolates on gold to control interfacial characteristics; these tailored interfaces, in turn, control the adsorption of proteins and the attachment of cells. The technique uses replica molding in poly(dimethylsiloxane) molds having micrometer-scale relief patterns on their surfaces to form a contoured film of polyurethane supported on a glass slide. Evaporation of a thin (< 12 nm) film of gold on this surface-contoured polyurethane provides an optically transparent substrate, on which SAMs of terminally functionalized alkanethiolates can be formed. In one procedure, a flat poly(dimethylsiloxane) stamp was used to form a SAM of hexadecanethiolate on the raised plateaus of the contoured surface by contact printing hexadecanethiol [HS(CH2)15CH3]; a SAM terminated in tri(ethylene glycol) groups was subsequently formed on the bare gold remaining in the grooves by immersing the substrate in a solution of a second alkanethiol [HS(CH2)11(OCH2CH2)3OH]. Then this patterned substrate was immersed in a solution of fibronectin, the protein adsorbed only on the methyl-terminated plateau regions of the substrate [the tri(ethylene glycol)-terminated regions resisted the adsorption of protein]; bovine capillary endothelial cells attached only on the regions that adsorbed fibronectin. A complementary procedure confined protein adsorption and cell attachment to the grooves in this substrate.

Ordered restriction endonuclease maps of yeast artificial chromosomes created by optical mapping on surfaces.

We have developed a surface mounting technology for the rapid construction of ordered restriction maps from individual DNA molecules. Optical restriction maps constructed from yeast artificial chromosome DNA molecules mounted on specially derivatized glass surfaces are accurate and reproducible, and the technology is amenable to automation. The mounting procedures described here should also be useful for fluorescence in situ hybridization studies. We believe these improvements to optical mapping will further stimulate the development of nonelectrophoretic approaches to genome analysis.

A secondary-structure model for the self-cleaving region of Neurospora VS RNA.

Neurospora VS RNA performs an RNA-mediated self-cleavage reaction whose products contain 2',3'-cyclic phosphate and 5'-hydroxyl termini. This reaction is similar to those of hammerhead, hairpin, and hepatitis delta virus ribozymes; however, VS RNA is not similar in sequence to these other self-cleaving motifs. Here we propose a model for the secondary structure of the self-cleaving region of VS RNA, supported by site-directed mutagenesis and chemical modification structure probing data. The secondary structure of VS RNA is distinct from those of the other naturally occurring RNA self-cleaving domains. In addition to a unique secondary structure, several Mg-dependent interactions occur during the folding of VS RNA into its active tertiary conformation.

Self-organized defensive behavior in honeybees

We investigated the defensive behavior of honeybees under controlled experimental conditions. During an attack on two identical targets, the spatial distribution of stings varied as a function of the total number of stings, evincing the classic “pitchfork bifurcation” phenomenon of nonlinear dynamics. The experimental results support a model of defensive behavior based on a self-organizing mechanism. The model helps to explain several of the characteristic features of the honeybee defensive response: (i) the ability of the colony to localize and focus its attack, (ii) the strong variability between different hives in the intensity of attack, as well as (iii) the variability observed within the same hive, and (iv) the ability of the colony to amplify small differences between the targets.

Residual delay maps unveil global patterns of atmospheric nonlinearity and produce improved local forecasts

We use residual-delay maps of observational field data for barometric pressure to demonstrate the structure of latitudinal gradients in nonlinearity in the atmosphere. Nonlinearity is weak and largely lacking in tropical and subtropical sites and increases rapidly into the temperate regions where the time series also appear to be much noisier. The degree of nonlinearity closely follows the meridional variation of midlatitude storm track frequency. We extract the specific functional form of this nonlinearity, a V shape in the lagged residuals that appears to be a basic feature of midlatitude synoptic weather systems associated with frontal passages. We present evidence that this form arises from the relative time scales of high-pressure versus low-pressure events. Finally, we show that this nonlinear feature is weaker in a well regarded numerical forecast model (European Centre for Medium-Range Forecasts) because small-scale temporal and spatial variation is smoothed out in the grided inputs. This is significant, in that it allows us to demonstrate how application of statistical corrections based on the residual-delay map may provide marked increases in local forecast accuracy, especially for severe weather systems.

Coexistence of linear zones and pinwheels within orientation maps in cat visual cortex

Revealing the layout of cortical maps is important both for understanding the processes involved in their development and for uncovering the mechanisms underlying neural computation. The typical organization of orientation maps in the cat visual cortex is radial; complete orientation cycles are mapped around orientation singularities. In contrast, long linear zones of orientation representation have been detected in the primary visual cortex of the tree shrew. In this study, we searched for the existence of long linear sequences and wide linear zones within orientation preference maps of the cat visual cortex. Optical imaging based on intrinsic signals was used. Long linear sequences and wide linear zones of preferred orientation were occasionally detected along the border between areas 17 and 18, as well as within area 18. Adjacent zones of distinct radial and linear organizations were observed across area 18 of a single hemisphere. However, radial and linear organizations were not necessarily segregated; long (7.5 mm) linear sequences of preferred orientation were found embedded within a typical pinwheel-like organization of orientation. We conclude that, although the radial organization is dominant, perfectly linear organization may develop and perform the processing related to orientation in the cat visual cortex.

Direct observation of the enhancement of noncooperative protein self-assembly by macromolecular crowding: Indefinite linear self-association of bacterial cell division protein FtsZ

Recent measurements of sedimentation equilibrium and sedimentation velocity have shown that the bacterial cell division protein FtsZ self-associates to form indefinitely long rod-like linear aggregates in the presence of GDP and Mg2+. In the present study, the newly developed technique of non-ideal tracer sedimentation equilibrium was used to measure the effect of high concentrations—up to 150 g/liter—of each of two inert “crowder” proteins, cyanmethemoglobin or BSA, on the thermodynamic activity and state of association of dilute FtsZ under conditions inhibiting (−Mg2+) and promoting (+Mg2+) FtsZ self-association. Analysis of equilibrium gradients of both FtsZ and crowder proteins indicates that, under the conditions of the present experiment, FtsZ interacts with each of the two crowder proteins essentially entirely via steric repulsion, which may be accounted for quantitatively by a simple model in which hemoglobin, albumin, and monomeric FtsZ are modeled as effective spherical hard particles, and each oligomeric species of FtsZ is modeled as an effective hard spherocylinder. The functional dependence of the sedimentation of FtsZ on the concentrations of FtsZ and either crowder indicates that, in the presence of high concentrations of crowder, both the weight-average degree of FtsZ self-association and the range of FtsZ oligomer sizes present in significant abundance are increased substantially.

Self-assembled monolayers from a designed combinatorial library of de novo β-sheet proteins

A variety of naturally occurring biomaterials owe their unusual structural and mechanical properties to layers of β-sheet proteins laminated between layers of inorganic mineral. To explore the possibility of fabricating novel two-dimensional protein layers, we studied the self-assembly properties of de novo proteins from a designed combinatorial library. Each protein in the library has a distinct 63 amino acid sequence, yet they all share an identical binary pattern of polar and nonpolar residues, which was designed to favor the formation of six-stranded amphiphilic β-sheets. Characterization of proteins isolated from the library demonstrates that (i) they self assemble into monolayers at an air/water interface; (ii) the monolayers are dominated by β-sheet secondary structure, as shown by both circular dichroism and infrared spectroscopies; and (iii) the measured areas (500- 600 Å2) of individual protein molecules in the monolayers match those expected for proteins folded into amphiphilic β-sheets. The finding that similar structures are formed by distinctly different protein sequences suggests that assembly into β-sheet monolayers can be encoded by binary patterning of polar and nonpolar amino acids. Moreover, because the designed binary pattern is compatible with a wide variety of different sequences, it may be possible to fabricate β-sheet monolayers by using combinations of side chains that are explicitly designed to favor particular applications of novel biomaterials.

How amoeboids self-organize into a fruiting body: Multicellular coordination in Dictyostelium discoideum

When individual amoebae of the cellular slime mold Dictyostelium discoideum are starving, they aggregate to form a multicellular migrating slug, which moves toward a region suitable for culmination. The culmination of the morphogenesis involves complex cell movements that transform a mound of cells into a globule of spores on a slender stalk. The movement has been likened to a “reverse fountain,” whereby prestalk cells in the upper part form a stalk that moves downwards and anchors to the substratum, while prespore cells in the lower part move upwards to form the spore head. So far, however, no satisfactory explanation has been produced for this process. Using a computer simulation that we developed, we now demonstrate that the processes that are essential during the earlier stages of the morphogenesis are in fact sufficient to produce the dynamics of the culmination stage. These processes are cAMP signaling, differential adhesion, cell differentiation, and production of extracellular matrix. Our model clarifies the processes that generate the observed cell movements. More specifically, we show that periodic upward movements, caused by chemotactic motion, are essential for successful culmination, because the pressure waves they induce squeeze the stalk downwards through the cell mass. The mechanisms revealed by our model have a number of self-organizing and self-correcting properties and can account for many previously unconnected and unexplained experimental observations.

Invariance and restriction toward a limited set of self-antigens characterize neonatal IgM antibody repertoires and prevail in autoreactive repertoires of healthy adults.

Analysis of the reactivity of IgM with self-antigens in tissues by a quantitative immunoblotting technique showed striking invariance among newborns in the human and in the mouse. The self-reactive repertoire of IgM of adults was also markedly conserved; it comprised most anti-self reactivities that prevailed among neonates. Multivariate analysis confirmed the homogeneity of IgM repertoires of neonates toward self- and non-self-antigens. Multivariate analysis discriminated between newborn and adult repertoires for reactivity with two of five sources of self-proteins and with non-self-antigens. Our observations support the concept that naturally activated B lymphocytes are selected early in development and throughout life for reactivity with a restricted set of self-antigens.