Role of nucleosome positioning in 3D chromatin organization and loop formation.

We present a physics-based polymer model that can investigate 3D organization of chromatin accounting for DNA elasticity, DNA-bending due to nucleosomes, and 1D organization of nucleosomes along DNA. We find that the packing density of chromatin oscillates between densities corresponding to highly folded and extended configurations as we change the nucleosome organization (length of linker DNA). We compute the looping probability of chromatin and show that the presence of nucleosomes increases the looping probability of the chain compared to that of a bare DNA. We also show that looping probability has a large variability depending on the nature of nucleosome organization and density of linker histones.

Role of transcription factor-mediated nucleosome disassembly in PHO5 gene expression.

Studying nucleosome dynamics in promoter regions is crucial for understanding gene regulation. Nucleosomes regulate gene expression by sterically occluding transcription factors (TFs) and other non-histone proteins accessing genomic DNA. How the binding competition between nucleosomes and TFs leads to transcriptionally compatible promoter states is an open question. Here, we present a computational study of the nucleosome dynamics and organization in the promoter region of PHO5 gene in Saccharomyces cerevisiae. Introducing a model for nucleosome kinetics that takes into account ATP-dependent remodeling activity, DNA sequence effects, and kinetics of TFs (Pho4p), we compute the probability of obtaining different "promoter states" having different nucleosome configurations. Comparing our results with experimental data, we argue that the presence of local remodeling activity (LRA) as opposed to basal remodeling activity (BRA) is crucial in determining transcriptionally active promoter states. By modulating the LRA and Pho4p binding rate, we obtain different mRNA distributions-Poisson, bimodal, and long-tail. Through this work we explain many features of the PHO5 promoter such as sequence-dependent TF accessibility and the role of correlated dynamics between nucleosomes and TFs in opening/coverage of the TATA box. We also obtain possible ranges for TF binding rates and the magnitude of LRA.

Steady-state analysis of glucose repression reveals hierarchical expression of proteins under Mig1p control in Saccharomyces cerevisiae.

Glucose repression is a global transcriptional regulatory mechanism commonly observed in micro-organisms for the repression of enzymes that are not essential for glucose metabolism. In Saccharomyces cerevisiae, Mig1p, a homologue of Wilms' tumour protein, is a global repressor protein dedicated to glucose repression. Mig1p represses genes either by binding directly to the upstream repression sequence of structural genes or by indirectly repressing a transcriptional activator, such as Gal4p. In addition, some genes are repressed by both of the above mechanisms. This raises a fundamental question regarding the physiological relevance of the varied mechanisms of repression that exist involving Mig1p. We address this issue by comparing two well-known glucose-repression systems, that is, SUC2 and GAL gene expression systems, which encompass all the above three mechanisms. We demonstrate using steady-state analysis that these mechanisms lead to a hierarchical glucose repression profile of different family of genes. This switch over from one carbon source to another is well-calibrated as a function of glucose concentration through this hierarchical transcriptional response. The mechanisms prevailing in this repression system can achieve amplification and sensitivity, as observed in the well-characterized MAPK (mitogen-activated protein kinase) cascade system, albeit through a different structure. A critical feature of repression predicted by our steady-state model for the mutant strain of S. cerevisiae lacking Gal80p agrees well with the data reported here as well as that available in the literature.

A steady-state modeling approach to validate an in vivo mechanism of the GAL regulatory network in Saccharomyces cerevisiae.

Cellular regulation is a result of complex interactions arising from DNA-protein and protein-protein binding, autoregulation, and compartmentalization and shuttling of regulatory proteins. Experiments in molecular biology have identified these mechanisms recruited by a regulatory network. Mathematical models may be used to complement the knowledge-base provided by in vitro experimental methods. Interactions identified by in vitro experiments can lead to the hypothesis of multiple candidate models explaining the in vivo mechanism. The equilibrium dissociation constants for the various interactions and the total component concentration constitute constraints on the candidate models. In this work, we identify the most plausible in vivo network by comparing the output response to the experimental data. We demonstrate the methodology using the GAL system of Saccharomyces cerevisiae for which the steady-state analysis reveals that Gal3p neither dimerizes nor shuttles between the cytoplasm and the nucleus.