Equilibrium distributions of topological states in circular DNA: Interplay of supercoiling and knotting

Two variables define the topological state of closed double-stranded DNA: the knot type, K, and ΔLk, the linking number difference from relaxed DNA. The equilibrium distribution of probabilities of these states, P(ΔLk, K), is related to two conditional distributions: P(ΔLk|K), the distribution of ΔLk for a particular K, and P(K|ΔLk) and also to two simple distributions: P(ΔLk), the distribution of ΔLk irrespective of K, and P(K). We explored the relationships between these distributions. P(ΔLk, K), P(ΔLk), and P(K|ΔLk) were calculated from the simulated distributions of P(ΔLk|K) and of P(K). The calculated distributions agreed with previous experimental and theoretical results and greatly advanced on them. Our major focus was on P(K|ΔLk), the distribution of knot types for a particular value of ΔLk, which had not been evaluated previously. We found that unknotted circular DNA is not the most probable state beyond small values of ΔLk. Highly chiral knotted DNA has a lower free energy because it has less torsional deformation. Surprisingly, even at |ΔLk| > 12, only one or two knot types dominate the P(K|ΔLk) distribution despite the huge number of knots of comparable complexity. A large fraction of the knots found belong to the small family of torus knots. The relationship between supercoiling and knotting in vivo is discussed.

An empirical explanation of brightness

In this second part of our study on the mechanism of perceived brightness, we explore the effects of manipulating three-dimensional geometry. The additional scenes portrayed here demonstrate that the same luminance profile can elicit different sensations of brightness as a function of how the objects in the scene are arranged in space. This further evidence confirms the implication of the scenes presented in the accompanying paper, namely that sensations of relative brightness—including standard demonstrations of simultaneous brightness contrast—cannot arise by computations of local contrast. The most plausible explanation of the full range of perceptual phenomena we have described is an empirical strategy that links the luminance profile in a visual stimulus with an association (the percept) that represents the profile’s most probable real-world source.

An equilibrium statistical model for the spreading phase of open-ocean convection

A “most probable state” equilibrium statistical theory for random distributions of hetons in a closed basin is developed here in the context of two-layer quasigeostrophic models for the spreading phase of open-ocean convection. The theory depends only on bulk conserved quantities such as energy, circulation, and the range of values of potential vorticity in each layer. The simplest theory is formulated for a uniform cooling event over the entire basin that triggers a homogeneous random distribution of convective towers. For a small Rossby deformation radius typical for open-ocean convection sites, the most probable states that arise from this theory strongly resemble the saturated baroclinic states of the spreading phase of convection, with a stabilizing barotropic rim current and localized temperature anomaly.

Automated de novo sequencing of proteins by tandem high-resolution mass spectrometry

A de novo sequencing program for proteins is described that uses tandem MS data from electron capture dissociation and collisionally activated dissociation of electrosprayed protein ions. Computer automation is used to convert the fragment ion mass values derived from these spectra into the most probable protein sequence, without distinguishing Leu/Ile. Minimum human input is necessary for the data reduction and interpretation. No extra chemistry is necessary to distinguish N- and C-terminal fragments in the mass spectra, as this is determined from the electron capture dissociation data. With parts-per-million mass accuracy (now available by using higher field Fourier transform MS instruments), the complete sequences of ubiquitin (8.6 kDa) and melittin (2.8 kDa) were predicted correctly by the program. The data available also provided 91% of the cytochrome c (12.4 kDa) sequence (essentially complete except for the tandem MS-resistant region K13–V20 that contains the cyclic heme). Uncorrected mass values from a 6-T instrument still gave 86% of the sequence for ubiquitin, except for distinguishing Gln/Lys. Extensive sequencing of larger proteins should be possible by applying the algorithm to pieces of ≈10-kDa size, such as products of limited proteolysis.

Building a dictionary for genomes: Identification of presumptive regulatory sites by statistical analysis

The availability of complete genome sequences and mRNA expression data for all genes creates new opportunities and challenges for identifying DNA sequence motifs that control gene expression. An algorithm, “MobyDick,” is presented that decomposes a set of DNA sequences into the most probable dictionary of motifs or words. This method is applicable to any set of DNA sequences: for example, all upstream regions in a genome or all genes expressed under certain conditions. Identification of words is based on a probabilistic segmentation model in which the significance of longer words is deduced from the frequency of shorter ones of various lengths, eliminating the need for a separate set of reference data to define probabilities. We have built a dictionary with 1,200 words for the 6,000 upstream regulatory regions in the yeast genome; the 500 most significant words (some with as few as 10 copies in all of the upstream regions) match 114 of 443 experimentally determined sites (a significance level of 18 standard deviations). When analyzing all of the genes up-regulated during sporulation as a group, we find many motifs in addition to the few previously identified by analyzing the subclusters individually to the expression subclusters. Applying MobyDick to the genes derepressed when the general repressor Tup1 is deleted, we find known as well as putative binding sites for its regulatory partners.

Training and search methods for speech recognition.

Speech recognition involves three processes: extraction of acoustic indices from the speech signal, estimation of the probability that the observed index string was caused by a hypothesized utterance segment, and determination of the recognized utterance via a search among hypothesized alternatives. This paper is not concerned with the first process. Estimation of the probability of an index string involves a model of index production by any given utterance segment (e.g., a word). Hidden Markov models (HMMs) are used for this purpose [Makhoul, J. & Schwartz, R. (1995) Proc. Natl. Acad. Sci. USA 92, 9956-9963]. Their parameters are state transition probabilities and output probability distributions associated with the transitions. The Baum algorithm that obtains the values of these parameters from speech data via their successive reestimation will be described in this paper. The recognizer wishes to find the most probable utterance that could have caused the observed acoustic index string. That probability is the product of two factors: the probability that the utterance will produce the string and the probability that the speaker will wish to produce the utterance (the language model probability). Even if the vocabulary size is moderate, it is impossible to search for the utterance exhaustively. One practical algorithm is described [Viterbi, A. J. (1967) IEEE Trans. Inf. Theory IT-13, 260-267] that, given the index string, has a high likelihood of finding the most probable utterance.

Nuclear Pore Complex Number and Distribution throughout the Saccharomyces cerevisiae Cell Cycle by Three-Dimensional Reconstruction from Electron Micrographs of Nuclear Envelopes

The number of nuclear pore complexes (NPCs) in individual nuclei of the yeast Saccharomyces cerevisiae was determined by computer-aided reconstruction of entire nuclei from electron micrographs of serially sectioned cells. Nuclei of 32 haploid cells at various points in the cell cycle were modeled and found to contain between 65 and 182 NPCs. Morphological markers, such as cell shape and nuclear shape, were used to determine the cell cycle stage of the cell being examined. NPC number was correlated with cell cycle stage to reveal that the number of NPCs increases steadily, beginning in G1-phase, suggesting that NPC assembly occurs continuously throughout the cell cycle. However, the accumulation of nuclear envelope observed during the cell cycle, indicated by nuclear surface area, is not continuous at the same rate, such that the density of NPCs per unit area of nuclear envelope peaks in apparent S-phase cells. Analysis of the nuclear envelope reconstructions also revealed no preferred NPC-to-NPC distance. However, NPCs were found in large clusters over regions of the nuclear envelope. Interestingly, clusters of NPCs were most pronounced in early mitotic nuclei and were found to be associated with the spindle pole bodies, but the functional significance of this association is unknown.