A unified statistical framework for sequence comparison and structure comparison

We present an approach for assessing the significance of sequence and structure comparisons by using nearly identical statistical formalisms for both sequence and structure. Doing so involves an all-vs.-all comparison of protein domains [taken here from the Structural Classification of Proteins (scop) database] and then fitting a simple distribution function to the observed scores. By using this distribution, we can attach a statistical significance to each comparison score in the form of a P value, the probability that a better score would occur by chance. As expected, we find that the scores for sequence matching follow an extreme-value distribution. The agreement, moreover, between the P values that we derive from this distribution and those reported by standard programs (e.g., blast and fasta validates our approach. Structure comparison scores also follow an extreme-value distribution when the statistics are expressed in terms of a structural alignment score (essentially the sum of reciprocated distances between aligned atoms minus gap penalties). We find that the traditional metric of structural similarity, the rms deviation in atom positions after fitting aligned atoms, follows a different distribution of scores and does not perform as well as the structural alignment score. Comparison of the sequence and structure statistics for pairs of proteins known to be related distantly shows that structural comparison is able to detect approximately twice as many distant relationships as sequence comparison at the same error rate. The comparison also indicates that there are very few pairs with significant similarity in terms of sequence but not structure whereas many pairs have significant similarity in terms of structure but not sequence.

Quantitative insight into proliferation and differentiation of oligodendrocyte type 2 astrocyte progenitor cells in vitro

As part of our attempts at understanding fundamental principles that underlie the generation of nondividing terminally differentiated progeny from dividing precursor cells, we have developed approaches to a quantitative analysis of proliferation and differentiation of oligodendrocyte type 2 astrocyte (O-2A) progenitor cells at the clonal level. Owing to extensive previous studies of clonal differentiation in this lineage, O-2A progenitor cells represent an excellent system for such an analysis. Previous studies have resulted in two competing hypotheses; one of them suggests that progenitor cell differentiation is symmetric, the other hypothesis introduces an asymmetric process of differentiation. We propose a general model that incorporates both such extreme hypotheses as special cases. Our analysis of experimental data has shown, however, that neither of these extreme cases completely explains the observed kinetics of O-2A progenitor cell proliferation and oligodendrocyte generation in vitro. Instead, our results indicate that O-2A progenitor cells become competent for differentiation after they complete a certain number of critical mitotic cycles that represent a period of symmetric development. This number varies from clone to clone and may be thought of as a random variable; its probability distribution was estimated from experimental data. Those O-2A cells that have undergone the critical divisions then may differentiate into an oligodendrocyte in each of the subsequent mitotic cycles with a certain probability, thereby exhibiting the asymmetric type of differentiation.

The nature of the excited state of the reaction center of photosystem II of green plants: A high-resolution fluorescence spectroscopy study

We studied the electronically excited state of the isolated reaction center of photosystem II with high-resolution fluorescence spectroscopy at 5 K and compared the obtained spectral features with those obtained earlier for the primary electron donor. The results show that there is a striking resemblance between the emitting and charge-separating states in the photosystem II reaction center, such as a very similar shape of the phonon wing with characteristic features at 19 and 80 cm−1, almost identical frequencies of a number of vibrational modes, a very similar double-Gaussian shape of the inhomogeneous distribution function, and relatively strong electron-phonon coupling for both states. We suggest that the emission at 5 K originates either from an exciton state delocalized over the inactive branch of the photosystem or from a fraction of the primary electron donor that is long-lived at 5 K. The latter possibility can be explained by a distribution of the free energy difference of the primary charge separation reaction around zero. Both possibilities are in line with the idea that the state that drives primary charge separation in the reaction center of photosystem II is a collective state, with contributions from all chlorophyll molecules in the central part of the complex.

A wholly empirical explanation of perceived motion

Because the retinal activity generated by a moving object cannot specify which of an infinite number of possible physical displacements underlies the stimulus, its real-world cause is necessarily uncertain. How, then, do observers respond successfully to sequences of images whose provenance is ambiguous? Here we explore the hypothesis that the visual system solves this problem by a probabilistic strategy in which perceived motion is generated entirely according to the relative frequency of occurrence of the physical sources of the stimulus. The merits of this concept were tested by comparing the directions and speeds of moving lines reported by subjects to the values determined by the probability distribution of all the possible physical displacements underlying the stimulus. The velocities reported by observers in a variety of stimulus contexts can be accounted for in this way.

An information theory model of hydrophobic interactions.

A molecular model of poorly understood hydrophobic effects is heuristically developed using the methods of information theory. Because primitive hydrophobic effects can be tied to the probability of observing a molecular-sized cavity in the solvent, the probability distribution of the number of solvent centers in a cavity volume is modeled on the basis of the two moments available from the density and radial distribution of oxygen atoms in liquid water. The modeled distribution then yields the probability that no solvent centers are found in the cavity volume. This model is shown to account quantitatively for the central hydrophobic phenomena of cavity formation and association of inert gas solutes. The connection of information theory to statistical thermodynamics provides a basis for clarification of hydrophobic effects. The simplicity and flexibility of the approach suggest that it should permit applications to conformational equilibria of nonpolar solutes and hydrophobic residues in biopolymers.

Brush border myosin-I truncated in the motor domain impairs the distribution and the function of endocytic compartments in an hepatoma cell line.

Myosins I, a ubiquitous monomeric class of myosins that exhibits actin-based motor properties, are associated with plasma and/or vesicular membranes and have been suggested as players for trafficking events between cell surface and intracellular membranous structures. To investigate the function of myosins 1, we have transfected a mouse hepatoma cell line (BWTG3) with cDNAs encoding the chicken brush border myosin-I (BBMI) and two variants truncated in the motor domain. One variant is deleted of the first 446 amino acids and thereby lacks the ATP binding site, whereas the other is deleted of the entire motor domain and lacks the ATP and actin binding sites. We have observed (i) that significant amounts of the truncated variants are recovered with membrane fractions after cell fractionation, (ii) that they codistribute with a compartment containing alpha2-macroglobulin internalized for 30 min as determined by fluorescent microscopy, (iii) that the production of BBMI-truncated variants impairs the distribution of the acidic compartment and ligands internalized for 30 min, and (iv) that the production of the truncated variant containing the actin binding site decreases the rate of alpha2-macroglobulin degradation whereas the production of the variant lacking the ATP binding site and the actin binding site increases the rate of a2-macroglobulin degradation. These observations indicate that the two truncated variants have a dominant negative effect on the distribution and the function of the endocytic compartments. We propose that an unidentified myosin-I might contribute to the distribution of endocytic compartments in a juxtanuclear position and/or to the regulation of the delivery of ligands to the degradative compartment in BWTG3 cells.