Cooperative DNA-binding by Bicoid provides a mechanism for threshold-dependent gene activation in the Drosophila embryo

The Bicoid morphogen directs pattern formation along the anterior-posterior (A-P) axis of the Drosophila embryo. Bicoid is distributed in a concentration gradient that decreases exponentially from the anterior pole, however, it transcribes target genes such as hunchback in a step-function-like pattern; the expression domain is uniform and has a sharply defined posterior boundary. A 'gradient-affinity' model proposed to explain Bicoid action states that (i) cooperative gene activation by Bicoid generates the sharp on/off switch for target gene transcription and (ii) target genes with different affinities for Bicoid are expressed at different positions along the A-P axis. Using an in vivo yeast assay and in vitro methods, we show that Bicoid binds DNA with pairwise cooperativity; Bicoid bound to a strong site helps Bicoid bind to a weak site. These results support the first aspect of the model, providing a mechanism by which Bicoid generates sharp boundaries of gene expression. However, contrary to the second aspect of the model, we find no significant difference between the affinity of Bicoid for the anterior gene hunchback and the posterior gene knirps, We propose, instead, that the arrangement of Bicoids bound to the target gene presents a unique signature to the transcription machinery that, in combination with overall affinity, regulates the extent of gene transcription along the A-P axis

The Drosophila Gap Gene Network Is Composed of Two Parallel Toggle Switches

Drosophila “gap” genes provide the first response to maternal gradients in the early fly embryo. Gap genes are expressed in a series of broad bands across the embryo during first hours of development. The gene network controlling the gap gene expression patterns includes inputs from maternal gradients and mutual repression between the gap genes themselves. In this study we propose a modular design for the gap gene network, involving two relatively independent network domains. The core of each network domain includes a toggle switch corresponding to a pair of mutually repressive gap genes, operated in space by maternal inputs. The toggle switches present in the gap network are evocative of the phage lambda switch, but they are operated positionally (in space) by the maternal gradients, so the synthesis rates for the competing components change along the embryo anterior-posterior axis. Dynamic model, constructed based on the proposed principle, with elements of fractional site occupancy, required 5–7 parameters to fit quantitative spatial expression data for gap gradients. The identified model solutions (parameter combinations) reproduced major dynamic features of the gap gradient system and explained gap expression in a variety of segmentation mutants

Effect of promoter architecture on the cell-to-cell variability in gene expression

According to recent experimental evidence, the architecture of a promoter, defined as the number, strength and regulatory role of the operators that control the promoter, plays a major role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effort that addresses the question of how changes in promoter architecture affect noise in gene expression in a systematic rather than case-by-case fashion. In this article, we make such a systematic investigation, based on a simple microscopic model of gene regulation that incorporates stochastic effects. In particular, we show how operator strength and operator multiplicity affect this variability. We examine different modes of transcription factor binding to complex promoters (cooperative, independent, simultaneous) and how each of these affects the level of variability in transcription product from cell-to-cell. We propose that direct comparison between in vivo single-cell experiments and theoretical predictions for the moments of the probability distribution of mRNA number per cell can discriminate between different kinetic models of gene regulation.Comment: 35 pages, 6 figures, Submitte

Novel derivative of aminobenzenesulfonamide (3c) induces apoptosis in colorectal cancer cells through ROS generation and inhibits cell migration

Background: Colorectal cancer (CRC) is the 3rd most common type of cancer worldwide. New anti-cancer agents are needed for treating late stage colorectal cancer as most of the deaths occur due to cancer metastasis. A recently developed compound, 3c has shown to have potent antitumor effect; however the mechanism underlying the antitumor effect remains unknown. Methods: 3c-induced inhibition of proliferation was measured in the absence and presence NAC using MTT in HT-29 and SW620 cells and xCELLigence RTCA DP instrument. 3c-induced apoptotic studies were performed using flow cytometry. 3c-induced redox alterations were measured by ROS production using fluorescence plate reader and flow cytometry and mitochondrial membrane potential by flow cytometry; NADPH and GSH levels were determined by colorimetric assays. Bcl2 family protein expression and cytochrome c release and PARP activation was done by western blotting. Caspase activation was measured by ELISA. Cell migration assay was done using the real time xCELLigence RTCA DP system in SW620 cells and wound healing assay in HT-29. Results: Many anticancer therapeutics exert their effects by inducing reactive oxygen species (ROS). In this study, we demonstrate that 3c-induced inhibition of cell proliferation is reversed by the antioxidant, N-acetylcysteine, suggesting that 3c acts via increased production of ROS in HT-29 cells. This was confirmed by the direct measurement of ROS in 3c-treated colorectal cancer cells. Additionally, treatment with 3c resulted in decreased NADPH and glutathione levels in HT-29 cells. Further, investigation of the apoptotic pathway showed increased release of cytochrome c resulting in the activation of caspase-9, which in turn activated caspase-3 and −6. 3c also (i) increased p53 and Bax expression, (ii) decreased Bcl2 and BclxL expression and (iii) induced PARP cleavage in human colorectal cancer cells. Confirming our observations, NAC significantly inhibited induction of apoptosis, ROS production, cytochrome c release and PARP cleavage. The results further demonstrate that 3c inhibits cell migration by modulating EMT markers and inhibiting TGFβ-induced phosphorylation of Smad2 and Samd3. Conclusions: Our findings thus demonstrate that 3c disrupts redox balance in colorectal cancer cells and support the notion that this agent may be effective for the treatment of colorectal cancer

E9-Im9 Colicin DNase−Immunity Protein Biomolecular Association in Water: A Multiple-Copy and Accelerated Molecular Dynamics Simulation Study

Protein−protein transient and dynamic interactions underlie all biological processes. The molecular dynamics (MD) of the E9 colicin DNase protein, its Im9 inhibitor protein, and their E9-Im9 recognition complex are investigated by combining multiple-copy (MC) MD and accelerated MD (aMD) explicit-solvent simulation approaches, after validation with crystalline-phase and solution experiments. Im9 shows higher flexibility than its E9 counterpart. Im9 displays a significant reduction of backbone flexibility and a remarkable increase in motional correlation upon E9 association. Im9 loops 23−31 and 54−64 open with respect to the E9-Im9 X-ray structure and show high conformational diversity. Upon association a large fraction (∼20 nm2) of E9 and Im9 protein surfaces become inaccessible to water. Numerous salt bridges transiently occurring throughout our six 50 ns long MC-MD simulations are not present in the X-ray model. Among these Im9 Glu31−E9 Arg96 and Im9 Glu41−Lys89 involve interface interactions. Through the use of 10 ns of Im9 aMD simulation, we reconcile the largest thermodynamic impact measured for Asp51Ala mutation with Im9 structure and dynamics. Lys57 acts as an essential molecular switch to shift Im9 surface loop towards an ideal configuration for E9 inhibition. This is achieved by switching Asp60−Lys57 and Asp62−Lys57 hydrogen bonds to Asp51−Lys57 salt bridge. E9-Im9 recognition involves shifts of conformational distributions, reorganization of intramolecular hydrogen bond patterns, and formation of new inter- and intramolecular interactions. The description of key transient biological interactions can be significantly enriched by the dynamic and atomic-level information provided by computer simulations

Thermodynamics-Based Models of Transcriptional Regulation by Enhancers: The Roles of Synergistic Activation, Cooperative Binding and Short-Range Repression

Quantitative models of cis-regulatory activity have the potential to improve our mechanistic understanding of transcriptional regulation. However, the few models available today have been based on simplistic assumptions about the sequences being modeled, or heuristic approximations of the underlying regulatory mechanisms. We have developed a thermodynamics-based model to predict gene expression driven by any DNA sequence, as a function of transcription factor concentrations and their DNA-binding specificities. It uses statistical thermodynamics theory to model not only protein-DNA interaction, but also the effect of DNA-bound activators and repressors on gene expression. In addition, the model incorporates mechanistic features such as synergistic effect of multiple activators, short range repression, and cooperativity in transcription factor-DNA binding, allowing us to systematically evaluate the significance of these features in the context of available expression data. Using this model on segmentation-related enhancers in Drosophila, we find that transcriptional synergy due to simultaneous action of multiple activators helps explain the data beyond what can be explained by cooperative DNA-binding alone. We find clear support for the phenomenon of short-range repression, where repressors do not directly interact with the basal transcriptional machinery. We also find that the binding sites contributing to an enhancer's function may not be conserved during evolution, and a noticeable fraction of these undergo lineage-specific changes. Our implementation of the model, called GEMSTAT, is the first publicly available program for simultaneously modeling the regulatory activities of a given set of sequences

Stable, Precise, and Reproducible Patterning of Bicoid and Hunchback Molecules in the Early Drosophila Embryo

Precise patterning of morphogen molecules and their accurate reading out are of key importance in embryonic development. Recent experiments have visualized distributions of proteins in developing embryos and shown that the gradient of concentration of Bicoid morphogen in Drosophila embryos is established rapidly after fertilization and remains stable through syncytial mitoses. This stable Bicoid gradient is read out in a precise way to distribute Hunchback with small fluctuations in each embryo and in a reproducible way, with small embryo-to-embryo fluctuation. The mechanisms of such stable, precise, and reproducible patterning through noisy cellular processes, however, still remain mysterious. To address these issues, here we develop the one- and three-dimensional stochastic models of the early Drosophila embryo. The simulated results show that the fluctuation in expression of the hunchback gene is dominated by the random arrival of Bicoid at the hunchback enhancer. Slow diffusion of Hunchback protein, however, averages out this intense fluctuation, leading to the precise patterning of distribution of Hunchback without loss of sharpness of the boundary of its distribution. The coordinated rates of diffusion and transport of input Bicoid and output Hunchback play decisive roles in suppressing fluctuations arising from the dynamical structure change in embryos and those arising from the random diffusion of molecules, and give rise to the stable, precise, and reproducible patterning of Bicoid and Hunchback distributions

Effects of macromolecular crowding on intracellular diffusion from a single particle perspective

Compared to biochemical reactions taking place in relatively well-defined aqueous solutions in vitro, the corresponding reactions happening in vivo occur in extremely complex environments containing only 60–70% water by volume, with the remainder consisting of an undefined array of bio-molecules. In a biological setting, such extremely complex and volume-occupied solution environments are termed ‘crowded’. Through a range of intermolecular forces and pseudo-forces, this complex background environment may cause biochemical reactions to behave differently to their in vitro counterparts. In this review, we seek to highlight how the complex background environment of the cell can affect the diffusion of substances within it. Engaging the subject from the perspective of a single particle’s motion, we place the focus of our review on two areas: (1) experimental procedures for conducting single particle tracking experiments within cells along with methods for extracting information from these experiments; (2) theoretical factors affecting the translational diffusion of single molecules within crowded two-dimensional membrane and three-dimensional solution environments. We conclude by discussing a number of recent publications relating to intracellular diffusion in light of the reviewed material