Search PubMed⌕ Search

Biomedical subjects

E D Siggia

Publications and source records attributed to E D Siggia.

12 recordsLinked to original sources

Regulatory element detection using correlation with expression.

We present here a new computational method for discovering cis-regulatory elements that circumvents the need to cluster genes based on their expression profiles. Based on a model in which upstream motifs contribute additively to the log-expression level of a gene, this method requires a single genome-wide set of expression ratios and the upstream sequence for each gene, and outputs statistically significant motifs. Analysis of publicly available expression data for Saccharomyces cerevisiae reveals several new putative regulatory elements, some of which plausibly control the early, transient induction of genes during sporulation. Known motifs generally have high statistical significance.

Algorithms↗

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.

Algorithms↗

Modeling RNA folding paths with pseudoknots: application to hepatitis delta virus ribozyme.

A quantitative understanding of nucleic acid hybridization is essential to many aspects of biotechnology, such as DNA microarrays, as well as to the structure and folding kinetics of RNA. However, predictions of nucleic acid secondary structures have long been impeded by the presence of helices interior to loops, so-called pseudoknots, which impose complex three-dimensional conformational constraints. In this paper we compute the pseudoknot free energies analytically in terms of known standard parameters, and we show how the results can be included in a kinetic Monte Carlo code to follow the succession of secondary structures during quenched or sequential folding. For the hepatitis delta virus ribozyme, we predict several nonnative stems on the folding path, characterize a kinetically trapped state, interpret several experimentally characterized mutations in terms of the folding path, and suggest how hybridization with other parts of the genome inactivates the newly formed ribozyme.

Base Sequence↗

Diffusion in inhomogeneous media: theory and simulations applied to whole cell photobleach recovery.

A continuum description for diffusion in a simple model for an inhomogeneous but isotropic media is derived and implemented numerically. The locally averaged density of diffusible marker is input from experiment to define the sample. Then a single additional parameter, the effective diffusion constant, permits the quantitative simulation of diffusive relaxation from any initial condition. Using this simulation, it is possible to model the recovery of a fluorescently tagged protein in the endoplasmic reticulum (ER) after photobleaching a substantial region of a live cell, and fit an effective diffusion constant which is a property both of the geometry of the ER and the marker. Such quantitative measurements permit inferences about the topology and internal organization of this organelle.

Algorithms↗

Regulatory element detection using a probabilistic segmentation model.

The availability of genome-wide mRNA expression data for organisms whose genome is fully sequenced provides a unique data set from which to decipher how transcription is regulated by the upstream control region of a gene. A new algorithm is presented which decomposes DNA sequence into the most probable "dictionary" of motifs or words. Identification of words is based on a probabilistic segmentation model in which the significance of longer words is deduced from the frequency of shorter words of various length. This eliminates the need for a separate set of reference data to define probabilities, and genome-wide applications are therefore possible. For the 6,000 upstream regulatory regions in the yeast genome, the 500 strongest motifs from a dictionary of size 1,200 match at a significance level of 15 standard deviations to a database of cis-regulatory elements. Analysis of sets of genes such as those up-regulated during sporulation reveals many new putative regulatory sites in addition to identifying previously known sites.

Algorithms↗

Kinetic analysis of secretory protein traffic and characterization of golgi to plasma membrane transport intermediates in living cells.

Quantitative time-lapse imaging data of single cells expressing the transmembrane protein, vesicular stomatitis virus ts045 G protein fused to green fluorescent protein (VSVG-GFP), were used for kinetic modeling of protein traffic through the various compartments of the secretory pathway. A series of first order rate laws was sufficient to accurately describe VSVG-GFP transport, and provided compartment residence times and rate constants for transport into and out of the Golgi complex and delivery to the plasma membrane. For ER to Golgi transport the mean rate constant (i.e., the fraction of VSVG-GFP moved per unit of time) was 2.8% per min, for Golgi to plasma membrane transport it was 3.0% per min, and for transport from the plasma membrane to a degradative site it was 0.25% per min. Because these rate constants did not change as the concentration of VSVG-GFP in different compartments went from high (early in the experiment) to low (late in the experiment), secretory transport machinery was never saturated during the experiments. The processes of budding, translocation, and fusion of post-Golgi transport intermediates carrying VSVG- GFP to the plasma membrane were also analyzed using quantitative imaging techniques. Large pleiomorphic tubular structures, rather than small vesicles, were found to be the primary vehicles for Golgi to plasma membrane transport of VSVG-GFP. These structures budded as entire domains from the Golgi complex and underwent dynamic shape changes as they moved along microtubule tracks to the cell periphery. They carried up to 10,000 VSVG-GFP molecules and had a mean life time in COS cells of 3.8 min. In addition, they fused with the plasma membrane without intersecting other membrane transport pathways in the cell. These properties suggest that the post-Golgi intermediates represent a unique transport organelle for conveying large quantities of protein cargo from the Golgi complex directly to the plasma membrane.

Aluminum Compounds↗

Nuclear membrane dynamics and reassembly in living cells: targeting of an inner nuclear membrane protein in interphase and mitosis.

The mechanisms of localization and retention of membrane proteins in the inner nuclear membrane and the fate of this membrane system during mitosis were studied in living cells using the inner nuclear membrane protein, lamin B receptor, fused to green fluorescent protein (LBR-GFP). Photobleaching techniques revealed the majority of LBR-GFP to be completely immobilized in the nuclear envelope (NE) of interphase cells, suggesting a tight binding to heterochromatin and/or lamins. A subpopulation of LBR-GFP within ER membranes, by contrast, was entirely mobile and diffused rapidly and freely (D = 0. 41 +/- 0.1 microm2/s). High resolution confocal time-lapse imaging in mitotic cells revealed LBR-GFP redistributing into the interconnected ER membrane system in prometaphase, exhibiting the same high mobility and diffusion constant as observed in interphase ER membranes. LBR-GFP rapidly diffused across the cell within the membrane network defined by the ER, suggesting the integrity of the ER was maintained in mitosis, with little or no fragmentation and vesiculation. At the end of mitosis, nuclear membrane reformation coincided with immobilization of LBR-GFP in ER elements at contact sites with chromatin. LBR-GFP-containing ER membranes then wrapped around chromatin over the course of 2-3 min, quickly and efficiently compartmentalizing nuclear material. Expansion of the NE followed over the course of 30-80 min. Thus, selective changes in lateral mobility of LBR-GFP within the ER/NE membrane system form the basis for its localization to the inner nuclear membrane during interphase. Such changes, rather than vesiculation mechanisms, also underlie the redistribution of this molecule during NE disassembly and reformation in mitosis.

Animals↗

Driving proteins off DNA using applied tension.

Proteins that bind DNA so as to reduce its end-to-end length can be dissociated by application of force. The thermodynamics of this process are discussed, with special attention to the case of histones bound to DNA (i.e., a string of nucleosomes, or chromatin fiber). The histone octamer is predicted to be driven off chromatin fiber for tensions >2 piconewtons.

Biophysical Phenomena↗

Electrostatics of lipid bilayer bending.

The electrostatic contribution to spontaneous membrane curvature is calculated within Poisson-Boltzmann theory under a variety of assumptions and emphasizing parameters in the physiological range. Asymmetrical surface charges can be fixed with respect to bilayer midplane area or with respect to the lipid-water area, but induce curvatures of opposite signs. Unequal screening layers on the two sides of a vesicle (e.g., multivalent cationic proteins on one side and monovalent salt on the other) also induce bending. For reasonable parameters, tubules formed by electrostatically induced bending can have radii in the 50-100-nm range, often seen in many intracellular organelles. Thus membrane associated proteins may induce curvature and subsequent budding, without themselves being intrinsically curved. Furthermore, we derive the previously unexplored effects of respecting the strict conservation of charge within the interior of a vesicle. The electrostatic component of the bending modulus is small under most of our conditions and is left as an experimental parameter. The large parameter space of conditions is surveyed in an array of graphs.

Biophysical Phenomena↗

Polymer models of meiotic and mitotic chromosomes.

Polymers tied together by constraints exhibit an internal pressure; this idea is used to analyze physical properties of the bottle-brush-like chromosomes of meiotic prophase that consist of polymer-like flexible chromatin loops, attached to a central axis. Using a minimal number of experimental parameters, semiquantitative predictions are made for the bending rigidity, radius, and axial tension of such brushes, and the repulsion acting between brushes whose bristles are forced to overlap. The retraction of lampbrush loops when the nascent transcripts are stripped away, the oval shape of diplotene bivalents between chiasmata, and the rigidity of pachytene chromosomes are all manifestations of chromatin pressure. This two-phase (chromatin plus buffer) picture that suffices for meiotic chromosomes has to be supplemented by a third constituent, a chromatin glue to understand mitotic chromosomes, and explain how condensation can drive the resolution of entanglements. This process resembles a thermal annealing in that a parameter (the affinity of the glue for chromatin and/or the affinity of the chromatin for buffer) has to be tuned to achieve optimal results. Mechanical measurements to characterize this protein-chromatin matrix are proposed. Finally, the propensity for even slightly chemically dissimilar polymers to phase separate (cluster like with like) can explain the apparent segregation of the chromatin into A + T- and G + C-rich regions revealed by chromosome banding.

Animals↗

Fluctuations and supercoiling of DNA.

Frequently, DNA in vivo is organized into loops that are partially underwound and consequently form interwound helical supercoils. Methods from polymer statistical mechanics are used to show how the competition between entropy (thermal fluctuations) and elastic energy determines supercoil radius and pitch, in good agreement with recent experiments and simulations. Supercoil reorganization by means of slithering (reptation) of the DNA along the supercoil is argued to be a slow process. Extension of supercoiled DNA by an applied force shows a number of unexpected features, including coexistence of interwound and helical states.

DNA, Superhelical↗