Interaction of the myogenic determination factor myogenin with E12 and a DNA target: mechanism and kinetics

Simplifications and approximations in a single-gene circuit modeling

The absence of detailed knowledge about regulatory interactions makes the use of phenomenological assumptions mandatory in cell biology modeling. Furthermore, the challenges associated with the analysis of these models compel the implementation of mathematical approximations. However, the constraints these methods introduce to biological interpretation are sometimes neglected. Consequently, understanding these restrictions is a very important task for systems biology modeling. In this article, we examine the impact of such simplifications, taking the case of a single-gene autoinhibitory circuit; however, our conclusions are not limited solely to this instance. We demonstrate that models grounded in the same biological assumptions but described at varying levels of detail can lead to different outcomes, that is, different and contradictory phenotypes or behaviors. Indeed, incorporating specific molecular processes like translation and elongation into the model can introduce instabilities and oscillations not seen when these processes are assumed to be instantaneous. Furthermore, incorporating a detailed description of promoter dynamics, usually described by a phenomenological regulatory function, can lead to instability, depending on the cooperative binding mechanism that is acting. Consequently, although the use of a regulating function facilitates model analysis, it may mask relevant aspects of the system's behavior. In particular, we observe that the two cooperative binding mechanisms, both compatible with the same sigmoidal function, can lead to different phenotypes, such as transcriptional oscillations with different oscillation frequencies.

Homodimeric and Heterodimeric Interactions among Vertebrate Basic Helix-Loop-Helix Transcription Factors

The basic helix–loop–helix transcription factor (bHLH TF) family is involved in tissue development, cell differentiation, and disease. These factors have transcriptionally positive, negative, and inactive functions by combining dimeric interactions among family members. The best known bHLH TFs are the E-protein homodimers and heterodimers with the tissue-specific TFs or ID proteins. These cooperative and dynamic interactions result in a complex transcriptional network that helps define the cell’s fate. Here, the reported dimeric interactions of 67 vertebrate bHLH TFs with other family members are summarized in tables, including specifications of the experimental techniques that defined the dimers. The compilation of these extensive data underscores homodimers of tissue-specific bHLH TFs as a central part of the bHLH regulatory network, with relevant positive and negative transcriptional regulatory roles. Furthermore, some sequence-specific TFs can also form transcriptionally inactive heterodimers with each other. The function, classification, and developmental role for all vertebrate bHLH TFs in four major classes are detailed.

How p53 Molecules Solve the Target DNA Search Problem: A Review

Interactions between DNA and DNA-binding proteins play an important role in many essential cellular processes. A key function of the DNA-binding protein p53 is to search for and bind to target sites incorporated in genomic DNA, which triggers transcriptional regulation. How do p53 molecules achieve "rapid" and "accurate" target search in living cells? The search dynamics of p53 were expected to include 3D diffusion in solution, 1D diffusion along DNA, and intersegmental transfer between two different DNA strands. Single-molecule fluorescence microscopy enabled the tracking of p53 molecules on DNA and the characterization of these dynamics quantitatively. Recent intensive single-molecule studies of p53 succeeded in revealing each of these search dynamics. Here, we review these studies and discuss the target search mechanisms of p53.

BET-seq: Binding energy topographies revealed by microfluidics and high-throughput sequencing

Biophysical models of transcriptional regulation rely on energetic measurements of the binding affinities between transcription factors (TFs) and target DNA binding sites. Historically, assays capable of measuring TF-DNA binding affinities have been relatively low-throughput (measuring ~10³ sequences in parallel) and have required significant specialized equipment, limiting their use to a handful of laboratories. Recently, we developed an experimental assay and analysis pipeline that allows measurement of binding energies between a single TF and up to 10⁶ DNA species in a single experiment (Binding Energy Topography by sequencing, or BET-seq) (Le et al., 2018). BET-seq employs the Mechanically Induced Trapping of Molecular Interactions (MITOMI) platform to purify DNA bound to a TF at equilibrium followed by high coverage sequencing to reveal relative differences in binding energy for each sequence. While we have previously used BET-seq to refine the binding affinity landscapes surrounding high-affinity DNA consensus target sites, we anticipate this technique will be applied in future work toward measuring a wide variety of TF-DNA landscapes. Here, we provide detailed instructions and general considerations for DNA library design, performing BET-seq assays, and analyzing the resulting data.

Formation and maintenance of nitrogen-fixing cell patterns in filamentous cyanobacteria

Cyanobacteria forming one-dimensional filaments are paradigmatic model organisms of the transition between unicellular and multicellular living forms. Under nitrogen-limiting conditions, in filaments of the genus Anabaena, some cells differentiate into heterocysts, which lose the possibility to divide but are able to fix environmental nitrogen for the colony. These heterocysts form a quasiregular pattern in the filament, representing a prototype of patterning and morphogenesis in prokaryotes. Recent years have seen advances in the identification of the molecular mechanism regulating this pattern. We use these data to build a theory on heterocyst pattern formation, for which both genetic regulation and the effects of cell division and filament growth are key components. The theory is based on the interplay of three generic mechanisms: local autoactivation, early long-range inhibition, and late long-range inhibition. These mechanisms can be identified with the dynamics of hetR, patS, and hetN expression. Our theory reproduces quantitatively the experimental dynamics of pattern formation and maintenance for wild type and mutants. We find that hetN alone is not enough to play the role as the late inhibitory mechanism: a second mechanism, hypothetically the products of nitrogen fixation supplied by heterocysts, must also play a role in late long-range inhibition. The preponderance of even intervals between heterocysts arises naturally as a result of the interplay between the timescales of genetic regulation and cell division. We also find that a purely stochastic initiation of the pattern, without a two-stage process, is enough to reproduce experimental observations.

Models of buffering of dosage imbalances in protein complexes

Background

Stoichiometric imbalances in macromolecular complexes can lead to altered function. Such imbalances stem from under- or over-expression of a subunit of a complex consequent to a deletion, duplication or regulatory mutation of an allele encoding the relevant protein. In some cases, the phenotypic perturbations induced by such alterations can be subtle or be lacking because nonlinearities in the process of protein complex assembly can provide some degree of buffering.

Results

We explore with biochemical models of increasing plausibility how buffering can be elicited. Specifically, we analyze the formation of a dimer AB and show that there are particular sets of parameters so that decreasing/increasing the input amount of either A or B translates into a non proportional (buffered) change of AB. The buffer effect also appears in higher-order structures provided that there are intermediate subcomplexes in the assembly process.

Conclusions

We highlight the importance of protein degradation and/or conformational inactivation for buffering to appear. The models sketched here have experimental support but can be further tested with existing biological resources.

Reviewers

This article was reviewed by Eugene Koonin, Berend Snel and Csaba Pal.

Comparative expression profiling identifies differential roles for Myogenin and p38α MAPK signaling in myogenesis

Skeletal muscle differentiation is mediated by a complex gene expression program requiring both the muscle-specific transcription factor Myogenin (Myog) and p38α MAPK (p38α) signaling. However, the relative contribution of Myog and p38α to the formation of mature myotubes remains unknown. Here, we have uncoupled the activity of Myog from that of p38α to gain insight into the individual roles of these proteins in myogenesis. Comparative expression profiling confirmed that Myog activates the expression of genes involved in muscle function. Furthermore, we found that in the absence of p38α signaling, Myog expression leads to the down-regulation of genes involved in cell cycle progression. Consistent with this, the expression of Myog is sufficient to induce cell cycle exit. Interestingly, p38α-defective, Myog-expressing myoblasts fail to form multinucleated myotubes, suggesting an important role for p38α in cell fusion. Through the analysis of p38α up-regulated genes, the tetraspanin CD53 was identified as a candidate fusion protein, a role confirmed both ex vivo in primary myoblasts, and in vivo during myofiber regeneration in mice. Thus, our study has revealed an unexpected role for Myog in mediating cell cycle exit and has identified an essential role for p38α in cell fusion through the up-regulation of CD53.

Topology and dynamics of the zebrafish segmentation clock core circuit

By combining biochemical, embryological, and mathematical approaches, this work uncovers an important role for protein-protein interactions in determining the dynamics of the somite-forming segmentation clock in vertebrates.

During vertebrate embryogenesis, the rhythmic and sequential segmentation of the body axis is regulated by an oscillating genetic network termed the segmentation clock. We describe a new dynamic model for the core pace-making circuit of the zebrafish segmentation clock based on a systematic biochemical investigation of the network's topology and precise measurements of somitogenesis dynamics in novel genetic mutants. We show that the core pace-making circuit consists of two distinct negative feedback loops, one with Her1 homodimers and the other with Her7:Hes6 heterodimers, operating in parallel. To explain the observed single and double mutant phenotypes of her1, her7, and hes6 mutant embryos in our dynamic model, we postulate that the availability and effective stability of the dimers with DNA binding activity is controlled in a “dimer cloud” that contains all possible dimeric combinations between the three factors. This feature of our model predicts that Hes6 protein levels should oscillate despite constant hes6 mRNA production, which we confirm experimentally using novel Hes6 antibodies. The control of the circuit's dynamics by a population of dimers with and without DNA binding activity is a new principle for the segmentation clock and may be relevant to other biological clocks and transcriptional regulatory networks.

The segmented pattern of the vertebral column, one of the defining features of the vertebrate body, is established during embryogenesis. The embryo's segments, called somites, form sequentially and rhythmically from head to tail. The periodicity of somite formation is regulated by the segmentation clock, a genetic oscillator that ticks in the posterior-most embryonic tissue: for each tick of the clock, one new bilateral pair of segments is made. The period of the clock appears to determine the number and the length of segments, but what controls this periodicity? In this article, we have investigated the interactions of three transcription factors that form the core of the clock's regulatory circuit, and have measured how the period of segmentation changes when these factors are mutated alone or in combination. We find that these three factors contribute to a “dimer cloud” that contains all possible dimeric combinations; however, only two dimers in this cloud can bind DNA, which allows them to directly regulate the oscillatory gene expression that underpins the periodicity of segment formation. Nevertheless, a mathematical model of the clock's dynamics based on our experimental findings indicates that the non-DNA-binding dimers also influence the stability, and hence the function, of the two DNA-binding dimers controlling the segmentation clock's period. Such involvement of non-DNA-binding dimers is a novel regulatory principle for the segmentation clock, which might also be a general mechanism that operates in other biological clocks.

Snail regulates MyoD binding-site occupancy to direct enhancer switching and differentiation-specific transcription in myogenesis

In skeletal myogenesis, the transcription factor MyoD activates distinct transcriptional programs in progenitors compared to terminally differentiated cells. Using ChIP-Seq and gene expression analyses, we show that in primary myoblasts, Snail-HDAC1/2 repressive complex binds and excludes MyoD from its targets. Notably, Snail binds E box motifs that are G/C rich in their central dinucleotides, and such sites are almost exclusively associated with genes expressed during differentiation. By contrast, Snail does not bind the A/T-rich E boxes associated with MyoD targets in myoblasts. Thus, Snai1-HDAC1/2 prevent MyoD occupancy on differentiation-specific regulatory elements, and the change from Snail to MyoD binding often results in enhancer switching during differentiation. Furthermore, we show that a regulatory network involving myogenic regulatory factors (MRFs), Snai1/2, miR-30a, and miR-206 acts as a molecular switch that controls entry into myogenic differentiation. Together, these results reveal a regulatory paradigm that directs distinct gene expression programs in progenitors versus terminally differentiated cells.

Next generation microfluidic platforms for high-throughput protein biochemistry

DNA technologies such as cloning, DNA microarrays, and next generation sequencing have transformed the life sciences. Protein technologies on the other hand have not seen such explosive progress. This is mainly due to the inherent difficulty of working with proteins because of their manifold physical characteristics as opposed to the well behaved and well understood DNA polymer. Recent technological advancements have increased the throughput of protein biochemistry to levels where it is becoming of interest to systems biology. Here I review methods for high-throughput in situ synthesis and characterization of proteins and their integration with microfluidic devices. In the near future, the use of gene synthesis, microfluidic based protein synthesis and characterization will give rise to a resurgence of protein biochemistry in the current world of high-throughput genomics.

Stochastic simulation of delay-induced circadian rhythms in Drosophila

Circadian rhythms are ubiquitous in all eukaryotes and some prokaryotes. Several computational models with or without time delays have been developed for circadian rhythms. Exact stochastic simulations have been carried out for several models without time delays, but no exact stochastic simulation has been done for models with delays. In this paper, we proposed a detailed and a reduced stochastic model with delays for circadian rhythms in Drosophila based on two deterministic models of Smolen et al. and employed exact stochastic simulation to simulate circadian oscillations. Our simulations showed that both models can produce sustained oscillations and that the oscillation is robust to noise in the sense that there is very little variability in oscillation period although there are significant random fluctuations in oscillation peaks. Moreover, although average time delays are essential to simulation of oscillation, random changes in time delays within certain range around fixed average time delay cause little variability in the oscillation period. Our simulation results also showed that both models are robust to parameter variations and that oscillation can be entrained by light/dark circles. Our simulations further demonstrated that within a reasonable range around the experimental result, the rates that dclock and per promoters switch back and forth between activated and repressed sites have little impact on oscillation period.

Nitric Oxide Down-regulates Caveolin-3 Levels through the Interaction with Myogenin, Its Transcription Factor

Certain patients suffering from chronic diseases such as AIDS or cancer experience a constant cellular secretion of tumor necrosis factor alpha and other pro-inflammatory cytokines that results in a continuous release of nitric oxide (*NO) to the bloodstream. One immediate consequence of the deleterious action of *NO is weight loss and the progressive destruction of muscular mass in a process known as cachexia. We have previously reported that caveolin-3, a specific marker of muscle cells, becomes down-regulated by the action of *NO on muscular myotubes. We describe herein that the changes observed in caveolin-3 levels are due to the alteration of the DNA binding activity of the muscular transcription factor myogenin. In the presence of *NO, the binding of transcription factors from cell nuclear extracts of muscular tissues to the E boxes present in the caveolin-3 promoter become substantially reduced. When we purified recombinant myogenin and treated it with *NO donors, we could detect its S-nitrosylation by three independent methods, suggesting that very likely one of the cysteine residues of the molecule is being modified. Given the role of myogenin as a regulatory protein that determines the level of multiple muscle genes expressed during late myogenesis, our results might represent a novel mode of regulation of muscle development under conditions of nitric oxide-mediated toxicity.

A systems approach to measuring the binding energy landscapes of transcription factors

A major goal of systems biology is to predict the function of biological networks. Although network topologies have been successfully determined in many cases, the quantitative parameters governing these networks generally have not. Measuring affinities of molecular interactions in high-throughput format remains problematic, especially for transient and low-affinity interactions. We describe a high-throughput microfluidic platform that measures such properties on the basis of mechanical trapping of molecular interactions. With this platform we characterized DNA binding energy landscapes for four eukaryotic transcription factors; these landscapes were used to test basic assumptions about transcription factor binding and to predict their in vivo function.

Inhibition of MuSK expression by CREB interacting with a CRE-like element and MyoD

The type I receptor-like protein tyrosine kinase MuSK is essential for the neuromuscular junction formation. MuSK expression is tightly regulated during development, but the underlying mechanisms were unclear. Here we identified a novel mechanism by which MuSK expression may be regulated. A cyclic AMP response element (CRE)-like element in the 5'-flanking region of the MuSK gene binds to CREB1 (CRE-binding protein 1). Mutation of this element increases the MuSK promoter activity, suggesting a role for CREB1 in attenuation of MuSK expression. Interestingly, CREB mutants unable to bind to DNA also inhibit MuSK promoter activity, suggesting a CRE-independent inhibitory mechanism. In agreement, CREB1 could inhibit a mutant MuSK transgene reporter whose CRE site was mutated. We provide evidence that CREB interacts directly with MyoD, a myogenic factor essential for MuSK expression in muscle cells. Suppression of CREB expression by small interfering RNA increases MuSK promoter activity. These results demonstrate an important role for CREB1 in the regulation of MuSK expression.

Temperature and the expression of myogenic regulatory factors (MRFs) and myosin heavy chain isoforms during embryogenesis in the common carp Cyprinus carpio L

Embryos of the common carp, Cyprinus carpio L., were reared from fertilization of the eggs to inflation of the swim bladder in the larval stage at 18 and 25 degrees C. cRNA probes were used to detect transcripts of the myogenic regulatory factors MyoD, Myf-5 and myogenin, and five myosin heavy chain (MyHC) isoforms during development. The genes encoding Myf-5 and MyoD were switched on first in the unsegmented mesoderm, followed by myogenin as the somites developed. Myf-5 and MyoD transcripts were initially limited to the adaxial cells, but Myf-5 expression spread laterally into the presomitic mesoderm before somite formation. Two distinct bands of staining could be seen corresponding to the cellular fields of the forming somites, but as each furrow delineated, Myf-5 mRNA levels declined. Upon somite formation, MyoD expression spread laterally to encompass the full somite width. Expression of the myogenin gene was also switched on during somite formation, and expression of both transcripts persisted until the somites became chevron-shaped. Expression of MyoD was then downregulated shortly before myogenin. The expression patterns of the carp myogenic regulatory factor (MRF) genes most-closely resembled that seen in the zebrafish rather than the rainbow trout (where expression of MyoD remains restricted to the adaxial domain of the somite for a prolonged period) or the herring (where expression of MyoD persists longer than that of myogenin). Expression of two embryonic forms of MyHC began simultaneously at the 25-30 somite stage and continued until approximately two weeks post-hatch. However, the three adult isoforms of fast muscle MyHC were not detected in any stage examined, emphasizing a developmental gap that must be filled by other, as yet uncharacterised, MyHC isoform(s). No differences in the timing of expression of any mRNA transcripts were seen between temperature groups. A phylogenetic analysis of the MRFs was conducted using all available full-length amino acid sequences. A neighbour-joining tree indicated that all four members evolved from a common ancestral gene, which first duplicated into two lineages, each of which underwent a further duplication to produce Myf-5 and MyoD, and myogenin and MRF4. Parologous copies of MyoD from trout and Xenopus clustered closely together within clades, indicating recent duplications. By contrast, MyoD paralogues from gilthead seabream were more divergent, indicating a more-ancient duplication.

A New Class of Transcription Factors Mediates Brassinosteroid-Regulated Gene Expression in Arabidopsis

Brassinosteroids (BRs) signal through a plasma membrane-localized receptor kinase to regulate plant growth and development. We showed previously that a novel protein, BES1, accumulates in the nucleus in response to BRs, where it plays a role in BR-regulated gene expression; however, the mechanism by which BES1 regulates gene expression is unknown. In this study, we dissect BES1 subdomains and establish that BES1 is a transcription factor that binds to and activates BR target gene promoters both in vitro and in vivo. BES1 interacts with a basic helix-loop-helix protein, BIM1, to synergistically bind to E box (CANNTG) sequences present in many BR-induced promoters. Loss-of-function and gain-of-function mutants of BIM1 and its close family members display BR response phenotypes. Thus, BES1 defines a new class of plant-specific transcription factors that cooperate with transcription factors such as BIM1 to regulate BR-induced genes.

Id proteins negatively regulate basic helix-loop-helix transcription factor function by disrupting subnuclear compartmentalization

Id helix-loop-helix (HLH) proteins act as global regulators of metazoan cell fate, cell growth, and differentiation. They heterodimerize with and inhibit the DNA-binding function of members of the basic helix-loop-helix (bHLH) family of transcription factors. Using real time fluorescence microscopy techniques in single living cells, we show here that nuclear pools of chromatin-associated bHLH transcription factor are freely exchangeable and in constant flux. The existence of a dynamic equilibrium between DNA-bound and free bHLH protein is also directly demonstrable in vitro. By contrast, Id protein is not associated with any subcellular, macromolecular structures and displays a more highly mobile, diffuse nuclear-cytoplasmic distribution. When co-expressed with antagonist Id protein, the chromatin-associated sublocalization of bHLH protein is abolished, and there is an accompanying 100-fold increase in its nuclear mobility to a level expected for freely diffusible Id-bHLH heterodimer. These results suggest that nuclear Id protein acts by sequestering pools of transiently diffusing bHLH protein to prevent reassociation with chromatin domains. Such a mechanism would explain how Id proteins are able to overcome the large DNA-binding free energy of bHLH proteins that is necessary to accomplish their inhibitory effect.