Motifs, modules and games in bacteria

Competition and evolutionary selection among core regulatory motifs in gene expression control

Gene products that are beneficial in one environment may become burdensome in another, prompting the emergence of diverse regulatory schemes that carry their own bioenergetic cost. By ensuring that regulators are only expressed when needed, we demonstrate that autoregulation generally offers an advantage in an environment combining mutation and time-varying selection. Whether positive or negative feedback emerges as dominant depends primarily on the demand for the target gene product, typically to ensure that the detrimental impact of inevitable mutations is minimized. While self-repression of the regulator curbs the spread of these loss-of-function mutations, self-activation instead facilitates their propagation. By analyzing the transcription network of multiple model organisms, we reveal that reduced bioenergetic cost may contribute to the preferential selection of autoregulation among transcription factors. Our results not only uncover how seemingly equivalent regulatory motifs have fundamentally different impact on population structure, growth dynamics, and evolutionary outcomes, but they can also be leveraged to promote the design of evolutionarily robust synthetic gene circuits.

The pathophysiology of major depressive disorder through the lens of systems biology: Network analysis of the psycho-immune-neuroendocrine physiome

BACKGROUND/AIMS

The psycho-immune-neuroendocrine (PINE) network is a predominantly physiological (metabolomic) model constructed from the literature, inter-linking multiple biological processes associated with major depressive disorder (MDD), thereby integrating putative mechanistic pathways for MDD into a single network.

MATERIAL AND METHODS

Previously published metabolomic pathways for the PINE network based on literature searches conducted in 1991-2021 were used to construct an edge table summarizing all physiological pathways in pairs of origin nodes and target nodes. The Gephi software program was used to calculate network metrics from the edge table, including total degree and centrality measures, to ascertain key network nodes and construct a directed network graph.

RESULTS

An edge table and directional network graph of physiological relationships in the PINE network is presented. The network has properties consistent with complex biological systems, with analysis yielding key network nodes comprising pro-inflammatory cytokines (TNF- α, IL6 and IL1), glucocorticoids and corticotropin releasing hormone (CRH). These may represent central structural and regulatory elements in the context of MDD.

CONCLUSION

The identified hubs have a high degree of connection and are known to play roles in the progression from health to MDD. These nodes represent strategic targets for therapeutic intervention or prevention. Future work is required to build a weighted and dynamic simulation of the network PINE.

Hardware, Software, and Wetware Codesign Environment for Synthetic Biology

Synthetic biology is the process of forward engineering living systems. These systems can be used to produce biobased materials, agriculture, medicine, and energy. One approach to designing these systems is to employ techniques from the design of embedded electronics. These techniques include abstraction, standards, modularity, automated design, and formal semantic models of computation. Together, these elements form the foundation of "biodesign automation," where software, robotics, and microfluidic devices combine to create exciting biological systems of the future. This paper describes a "hardware, software, wetware" codesign vision where software tools can be made to act as "genetic compilers" that transform high-level specifications into engineered "genetic circuits" (wetware). This is followed by a process where automation equipment, well-defined experimental workflows, and microfluidic devices are explicitly designed to house, execute, and test these circuits (hardware). These systems can be used as either massively parallel experimental platforms or distributed bioremediation and biosensing devices. Next, scheduling and control algorithms (software) manage these systems' actual execution and data analysis tasks. A distinguishing feature of this approach is how all three of these aspects (hardware, software, and wetware) may be derived from the same basic specification in parallel and generated to fulfill specific cost, performance, and structural requirements.

Data-driven network models for genetic circuits from time-series data with incomplete measurements

Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements of state dynamics. We develop a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a previously reported in vitro genetic circuit and a new E. coli-based transcriptional event detector.

Divide-and-conquer: machine-learning integrates mammalian and viral traits with network features to predict virus-mammal associations

Our knowledge of viral host ranges remains limited. Completing this picture by identifying unknown hosts of known viruses is an important research aim that can help identify and mitigate zoonotic and animal-disease risks, such as spill-over from animal reservoirs into human populations. To address this knowledge-gap we apply a divide-and-conquer approach which separates viral, mammalian and network features into three unique perspectives, each predicting associations independently to enhance predictive power. Our approach predicts over 20,000 unknown associations between known viruses and susceptible mammalian species, suggesting that current knowledge underestimates the number of associations in wild and semi-domesticated mammals by a factor of 4.3, and the average potential mammalian host-range of viruses by a factor of 3.2. In particular, our results highlight a significant knowledge gap in the wild reservoirs of important zoonotic and domesticated mammals' viruses: specifically, lyssaviruses, bornaviruses and rotaviruses.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

Phytoplankton Temporal Strategies Increase Entropy Production in a Marine Food Web Model

We develop a trait-based model founded on the hypothesis that biological systems evolve and organize to maximize entropy production by dissipating chemical and electromagnetic free energy over longer time scales than abiotic processes by implementing temporal strategies. A marine food web consisting of phytoplankton, bacteria, and consumer functional groups is used to explore how temporal strategies, or the lack thereof, change entropy production in a shallow pond that receives a continuous flow of reduced organic carbon plus inorganic nitrogen and illumination from solar radiation with diel and seasonal dynamics. Results show that a temporal strategy that employs an explicit circadian clock produces more entropy than a passive strategy that uses internal carbon storage or a balanced growth strategy that requires phytoplankton to grow with fixed stoichiometry. When the community is forced to operate at high specific growth rates near 2 d⁻¹, the optimization-guided model selects for phytoplankton ecotypes that exhibit complementary for winter versus summer environmental conditions to increase entropy production. We also present a new type of trait-based modeling where trait values are determined by maximizing entropy production rather than by random selection.

A Survey of Evolutionary Games in Biology

The applications of game theory in biology are numerous and include the comparison and modeling situations between bacteria, viruses, etc. This work provides insights about the connection between biology and evolving populations with classical and quantum evolutionary game theory and explains the benefits of unconventional computing methods in the study of such phenomena. In particular, the introduction of automata brings new possibilities into the decision-making process.

Stochastic models coupling gene expression and partitioning in cell division in Escherichia coli

Regulation of future RNA and protein numbers is a key process by which cells continuously best fit the environment. In bacteria, RNA and proteins exist in small numbers and their regulatory processes are stochastic. Consequently, there is cell-to-cell variability in these numbers, even between sister cells. Traditionally, the two most studied sources of this variability are gene expression and RNA and protein degradation, with evidence suggesting that the latter is subject to little regulation, when compared to the former. However, time-lapse microscopy and single molecule fluorescent tagging have produced evidence that cell division can also be a significant source of variability due to asymmetries in the partitioning of RNA and proteins. Relevantly, the impact of this noise differs from noise in production and degradation since, unlike these, it is not continuous. Rather, it occurs at specific time points, at which moment it can introduce major fluctuations. Several models have now been proposed that integrate noise from cell division, in addition to noise in gene expression, to mimic the dynamics of RNA and protein numbers of cell lineages. This is expected to be particularly relevant in genetic circuits, where significant fluctuations in one component protein, at specific time moments, are expected to perturb near-equilibrium states of the circuits, which can have long-lasting consequences. Here we review stochastic models coupling these processes in Escherichia coli, from single genes to small circuits.

When is a transcription factor a NAP?

Proteins that regulate transcription often also play an architectural role in the genome. Thus, it has been difficult to define with precision the distinctions between transcription factors and nucleoid-associated proteins (NAPs). Anachronistic descriptions of NAPs as 'histone-like' implied an organizational function in a bacterial chromatin-like complex. Definitions based on protein abundance, regulatory mechanisms, target gene number, or the features of their DNA-binding sites are insufficient as marks of distinction, and trying to distinguish transcription factors and NAPs based on their ranking within regulatory hierarchies or positions in gene-control networks is also unsatisfactory. The terms 'transcription factor' and 'NAP' are ad hoc operational definitions with each protein lying along a spectrum of structural and functional features extending from highly specific actors with few gene targets to those with a pervasive influence on the transcriptome. The Streptomyces BldC protein is used to illustrate these issues.

Added value of autoregulation and multi-step kinetics of transcription initiation

Bacterial gene expression regulation occurs mostly during transcription, which has two main rate-limiting steps: the close complex formation, when the RNA polymerase binds to an active promoter, and the subsequent open complex formation, after which it follows elongation. Tuning these steps' kinetics by the action of e.g. transcription factors, allows for a wide diversity of dynamics. For example, adding autoregulation generates single-gene circuits able to perform more complex tasks. Using stochastic models of transcription kinetics with empirically validated parameter values, we investigate how autoregulation and the multi-step transcription initiation kinetics of single-gene autoregulated circuits can be combined to fine-tune steady state mean and cell-to-cell variability in protein expression levels, as well as response times. Next, we investigate how they can be jointly tuned to control complex behaviours, namely, time counting, switching dynamics and memory storage. Overall, our finding suggests that, in bacteria, jointly regulating a single-gene circuit's topology and the transcription initiation multi-step dynamics allows enhancing complex task performance.

Drug treatment efficiency depends on the initial state of activation in nonlinear pathways

An accurate prediction of the outcome of a given drug treatment requires quantitative values for all parameters and concentrations involved as well as a detailed characterization of the network of interactions where the target molecule is embedded. Here, we present a high-throughput in silico screening of all potential networks of three interacting nodes to study the effect of the initial conditions of the network in the efficiency of drug inhibition. Our study shows that most network topologies can induce multiple dose-response curves, where the treatment has an enhanced, reduced or even no effect depending on the initial conditions. The type of dual response observed depends on how the potential bistable regimes interplay with the inhibition of one of the nodes inside a nonlinear pathway architecture. We propose that this dependence of the strength of the drug on the initial state of activation of the pathway may be affecting the outcome and the reproducibility of drug studies and clinical trials.

Systems Level Modeling of the Cell Cycle Using Budding Yeast

Proteins involved in the regulation of the cell cycle are highly conserved across all eukaryotes, and so a relatively simple eukaryote such as yeast can provide insight into a variety of cell cycle perturbations including those that occur in human cancer. To date, the budding yeast Saccharomyces cerevisiae has provided the largest amount of experimental and modeling data on the progression of the cell cycle, making it a logical choice for in-depth studies of this process. Moreover, the advent of methods for collection of high-throughput genome, transcriptome, and proteome data has provided a means to collect and precisely quantify simultaneous cell cycle gene transcript and protein levels, permitting modeling of the cell cycle on the systems level. With the appropriate mathematical framework and sufficient and accurate data on cell cycle components, it should be possible to create a model of the cell cycle that not only effectively describes its operation, but can also predict responses to perturbations such as variation in protein levels and responses to external stimuli including targeted inhibition by drugs. In this review, we summarize existing data on the yeast cell cycle, proteomics technologies for quantifying cell cycle proteins, and the mathematical frameworks that can integrate this data into representative and effective models. Systems level modeling of the cell cycle will require the integration of high-quality data with the appropriate mathematical framework, which can currently be attained through the combination of dynamic modeling based on proteomics data and using yeast as a model organism.

Fate-Regulating Circuits in Viruses: From Discovery to New Therapy Targets

Current antivirals effectively target diverse viruses at various stages of their life cycles. Nevertheless, curative therapy has remained elusive for important pathogens, such as human immunodeficiency virus type 1 (HIV-1) and herpesviruses, in large part due to viral latency and the evolution of resistance to existing therapies. Here, we review the discovery of viral master circuits: virus-encoded autoregulatory gene networks that autonomously control viral expression programs (i.e., between active, latent, and abortive fates). These circuits offer the opportunity for a new class of antivirals that could lead to intrinsic combination-antiviral therapies within a single molecule-evolutionary escape from such circuit-disrupting antivirals would require simultaneous evolution of both the viral cis regulatory element (e.g., the DNA-binding site) and the trans element (e.g., the transcription factor) in order for the virus to recapitulate a circuit that would not be disrupted. We review the architectures of these fate-regulating master circuits in HIV-1 and the human herpesvirus cytomegalovirus along with potential circuit-disruption strategies that may ultimately enable escape-resistant antiviral therapies.

Positive Autoregulation of an Acyl-Homoserine Lactone Quorum-Sensing Circuit Synchronizes the Population Response

Many proteobacteria utilize acyl-homoserine lactone quorum-sensing signals. At low population densities, cells produce a basal level of signal, and when sufficient signal has accumulated in the surrounding environment, it binds to its receptor, and quorum-sensing-dependent genes can be activated. A common characteristic of acyl-homoserine lactone quorum sensing is that signal production is positively autoregulated. We have examined the role of positive signal autoregulation in Pseudomonas aeruginosa. We compared population responses and individual cell responses in populations of wild-type P. aeruginosa to responses in a strain with the signal synthase gene controlled by an arabinose-inducible promoter so that signal was produced at a constant rate per cell regardless of cell population density. At a population level, responses of the wild type and the engineered strain were indistinguishable, but the responses of individual cells in a population of the wild type showed greater synchrony than the responses of the engineered strain. Although sufficient signal is required to activate expression of quorum-sensing-regulated genes, it is not sufficient for activation of certain genes, the late genes, and their expression is delayed until other conditions are met. We found that late gene responses were reduced in the engineered strain. We conclude that positive signal autoregulation is not a required element in acyl-homoserine lactone quorum sensing, but it functions to enhance synchrony of the responses of individuals in a population. Synchrony might be advantageous in some situations, whereas a less coordinated quorum-sensing response might allow bet hedging and be advantageous in other situations.

There are many quorum-sensing systems that involve a transcriptional activator, which responds to an acyl-homoserine lactone signal. In all of the examples studied, the gene coding for signal production is positively autoregulated by the signal, and it has even been described as essential for a quorum-sensing response. We have used the opportunistic pathogen Pseudomonas aeruginosa as a model to show that positive autoregulation is not required for a robust quorum-sensing response. We also show that positive autoregulation of signal production enhances the synchrony of the response. This information enhances our general understanding of the biological significance of how acyl-homoserine lactone quorum-sensing circuits are arranged.

Microscopic agents programmed by DNA circuits

Information stored in synthetic nucleic acids sequences can be used in vitro to create complex reaction networks with precisely programmed chemical dynamics. Here, we scale up this approach to program networks of microscopic particles (agents) dispersed in an enzymatic solution. Agents may possess multiple stable states, thus maintaining a memory and communicate by emitting various orthogonal chemical signals, while also sensing the behaviour of neighbouring agents. Using this approach, we can produce collective behaviours involving thousands of agents, for example retrieving information over long distances or creating spatial patterns. Our systems recapitulate some fundamental mechanisms of distributed decision making and morphogenesis among living organisms and could find applications in cases where many individual clues need to be combined to reach a decision, for example in molecular diagnostics.

Boosting functionality of synthetic DNA circuits with tailored deactivation

Molecular programming takes advantage of synthetic nucleic acid biochemistry to assemble networks of reactions, in vitro, with the double goal of better understanding cellular regulation and providing information-processing capabilities to man-made chemical systems. The function of molecular circuits is deeply related to their topological structure, but dynamical features (rate laws) also play a critical role. Here we introduce a mechanism to tune the nonlinearities associated with individual nodes of a synthetic network. This mechanism is based on programming deactivation laws using dedicated saturable pathways. We demonstrate this approach through the conversion of a single-node homoeostatic network into a bistable and reversible switch. Furthermore, we prove its generality by adding new functions to the library of reported man-made molecular devices: a system with three addressable bits of memory, and the first DNA-encoded excitable circuit. Specific saturable deactivation pathways thus greatly enrich the functional capability of a given circuit topology.

Nonlinearity in synthetic molecular circuits is usually achieved by manipulation of network topology or of production kinetics. Here, the authors achieve bistability and other nonlinear behaviours by manipulating the individual degradation rate laws of circuit components using saturable pathways.

Control of MarRAB Operon in Escherichia coli via Autoactivation and Autorepression

Choice of network topology for gene regulation has been a question of interest for a long time. How do simple and more complex topologies arise? In this work, we analyze the topology of the marRAB operon in Escherichia coli, which is associated with control of expression of genes associated with conferring resistance to low-level antibiotics to the bacterium. Among the 2102 promoters in E. coli, the marRAB promoter is the only one that encodes for an autoactivator and an autorepressor. What advantages does this topology confer to the bacterium? In this work, we demonstrate that, compared to control by a single regulator, the marRAB regulatory arrangement has the least control cost associated with modulating gene expression in response to environmental stimuli. In addition, the presence of dual regulators allows the regulon to exhibit a diverse range of dynamics, a feature that is not observed in genes controlled by a single regulator.

Biophysical and biological perspective in biosemiotics

The cell and its basic constituents are introduced here through a biophysical and information communication theoretic approach in biology and biosemiotics. With this purpose, the requirements of primordial cellular structures, single binding events, and signalling cascades are first mentioned stepwise, in relation to the model of the cellular sensing mechanism. This is followed by the concepts of cross reactions in sensing and pattern recognitions, wherein an information theoretic approach is addressed and the features of multicellularity are discussed along. Multicellularity is introduced as the path that leads to the loss of the direct causal relations. The loss of true causal relation is considered as a form of translation that enables meaning-encoded communication over the informative processes. In this sense, semiosis may not be exclusive. Synthetic biology is exemplified as a form of artificial selection mechanisms for the generation of 'self-reproducing' systems with information coding and processing machineries. These discussions are summarised at the end.

Sparse cliques trump scale-free networks in coordination and competition

Cooperative behavior, a natural, pervasive and yet puzzling phenomenon, can be significantly enhanced by networks. Many studies have shown how global network characteristics affect cooperation; however, it is difficult to understand how this occurs based on global factors alone, low-level network building blocks, or motifs are necessary. In this work, we systematically alter the structure of scale-free and clique networks and show, through a stochastic evolutionary game theory model, that cooperation on cliques increases linearly with community motif count. We further show that, for reactive stochastic strategies, network modularity improves cooperation in the anti-coordination Snowdrift game and the Prisoner's Dilemma game but not in the Stag Hunt coordination game. We also confirm the negative effect of the scale-free graph on cooperation when effective payoffs are used. On the flip side, clique graphs are highly cooperative across social environments. Adding cycles to the acyclic scale-free graph increases cooperation when multiple games are considered; however, cycles have the opposite effect on how forgiving agents are when playing the Prisoner's Dilemma game.

Single-cell kinetics of a repressilator when implemented in a single-copy plasmid

Synthetic genetic clocks, such as the Elowitz-Leibler repressilator, will be key regulatory components of future synthetic circuits. We constructed a single-copy repressilator (SCR) by implementing the original repressilator circuit on a single-copy F-plasmid. After verifying its functionality, we studied its behaviour as a function of temperature and compared it with that of the original low-copy-number repressilator (LCR). Namely, we compared the period of oscillations, functionality (the fraction of cells exhibiting oscillations) and robustness to internal fluctuations (the fraction of expected oscillations that would occur). We found that, under optimal temperature conditions, the dynamics of the two systems differs significantly, although qualitatively they respond similarly to temperature changes. Exception to this is in the functionality, in which the SCR is higher at lower temperatures but lower at higher temperatures. Next, by adding IPTG to the medium at low and high concentrations during microscopy sessions, we showed that the functionality of the SCR is more robust to external perturbations, which indicates that the oscillatory behaviour of the LCR can be disrupted by affecting only a few of the copies in a cell. We conclude that the SCR, the first functional, synthetic, single-copy, ring-type genetic clock, is more robust to lower temperatures and to external perturbations than the original LCR. The SCR will be of use in future synthetic circuits, since it complements the array of tasks that the LCR can perform.

Slow and steady wins the race: a bacterial exploitative competition strategy in fluctuating environments

One promising frontier for synthetic biology is the development of synthetic ecologies, whereby interacting species form an additional layer of connectivity for engineered gene circuits. Toward this goal, an important step is to understand different types of bacterial interactions in natural settings, among which competition is the most prevalent. By constructing a two-species population dynamics model, here, we mimicked bacterial growth in nature with resource-limited fluctuating environments and searched for optimal strategies for bacterial exploitative competition. In a simple game with two strategy options (constant or susceptible growth), we found that the species playing the constant growth strategy always outplays or is evenly matched with its competitor, suggesting that constant growth is a "no-loss" good bet. We also showed that adoption of sophisticated strategies enables a species to maximize its fitness when its competitor grows susceptibly. The pursuit of fitness maximization is, however, associated with potential loss if both species are capable of strategy adjustment, indicating an intrinsic risk-return trade-off. These findings offer new insights into bacterial competition and may also facilitate the engineering of microbial consortia for synthetic biology applications.

Synthetic biology at the interface of functional genomics

Functional genomics is considered a powerful tool that helps understand the relation between an organism's genotype and possible phenotypes. Volumes of data generated on several 'omics' platforms have revealed the network complexities underlying biological processes. Systems and synthetic biology have garnered much attention because of the ability to infer and comprehend the uncertainties associated with such complexities. Also, part-wise characterization of the network components (e.g. DNA, RNA, protein) has rendered an engineering perspective in life sciences to build modular and functional devices. This approach can be used to combat one of the many concerns of the world, i.e. in the area of biomedical translational research by designing and constructing novel therapeutic devices to intervene network perturbation in a diseased state to transform to a healthy state.

An evolutionarily conserved RNase-based mechanism for repression of transcriptional positive autoregulation

It is known that environmental context influences the degree of regulation at the transcriptional and post-transcriptional levels. However, the principles governing the differential usage and interplay of regulation at these two levels are not clear. Here, we show that the integration of transcriptional and post-transcriptional regulatory mechanisms in a characteristic network motif drives efficient environment-dependent state transitions. Through phenotypic screening, systems analysis, and rigorous experimental validation, we discovered an RNase (VNG2099C) in Halobacterium salinarum that is transcriptionally co-regulated with genes of the aerobic physiologic state but acts on transcripts of the anaerobic state. Through modelling and experimentation we show that this arrangement generates an efficient state-transition switch, within which RNase-repression of a transcriptional positive autoregulation (RPAR) loop is critical for shutting down ATP-consuming active potassium uptake to conserve energy required for salinity adaptation under aerobic, high potassium, or dark conditions. Subsequently, we discovered that many Escherichia coli operons with energy-associated functions are also putatively controlled by RPAR indicating that this network motif may have evolved independently in phylogenetically distant organisms. Thus, our data suggest that interplay of transcriptional and post-transcriptional regulation in the RPAR motif is a generalized principle for efficient environment-dependent state transitions across prokaryotes.

Gene co-expression modules as clinically relevant hallmarks of breast cancer diversity

Co-expression modules are groups of genes with highly correlated expression patterns. In cancer, differences in module activity potentially represent the heterogeneity of phenotypes important in carcinogenesis, progression, or treatment response. To find gene expression modules active in breast cancer subpopulations, we assembled 72 breast cancer-related gene expression datasets containing ∼5,700 samples altogether. Per dataset, we identified genes with bimodal expression and used mixture-model clustering to ultimately define 11 modules of genes that are consistently co-regulated across multiple datasets. Functionally, these modules reflected estrogen signaling, development/differentiation, immune signaling, histone modification, ERBB2 signaling, the extracellular matrix (ECM) and stroma, and cell proliferation. The Tcell/Bcell immune modules appeared tumor-extrinsic, with coherent expression in tumors but not cell lines; whereas most other modules, interferon and ECM included, appeared intrinsic. Only four of the eleven modules were represented in the PAM50 intrinsic subtype classifier and other well-established prognostic signatures; although the immune modules were highly correlated to previously published immune signatures. As expected, the proliferation module was highly associated with decreased recurrence-free survival (RFS). Interestingly, the immune modules appeared associated with RFS even after adjustment for receptor subtype and proliferation; and in a multivariate analysis, the combination of Tcell/Bcell immune module down-regulation and proliferation module upregulation strongly associated with decreased RFS. Immune modules are unusual in that their upregulation is associated with a good prognosis without chemotherapy and a good response to chemotherapy, suggesting the paradox of high immune patients who respond to chemotherapy but would do well without it. Other findings concern the ECM/stromal modules, which despite common themes were associated with different sites of metastasis, possibly relating to the “seed and soil” hypothesis of cancer dissemination. Overall, co-expression modules provide a high-level functional view of breast cancer that complements the “cancer hallmarks” and may form the basis for improved predictors and treatments.

Dynamic DNA-toolbox reaction circuits: a walkthrough

In living organisms, the integration of signals from the environment and the molecular computing leading to a cellular response are orchestrated by Gene Regulatory Networks (GRN). However, the molecular complexity of in vivo genetic regulation makes it next to impossible to describe in a quantitative manner. Reproducing, in vitro, reaction networks that could mimic the architecture and behavior of in vivo networks, yet lend themselves to mathematical modeling, represents a useful strategy to understand, and even predict, the function of GRN. In this paper, we define a set of in vitro, DNA-based molecular transformations that can be linked to each other in such a way that the product of one transformation can activate or inhibit the production of one or several other DNA compounds. Therefore, these reactions can be wired in arbitrary networks. This approach provides an experimental way to reproduce the dynamic features of genetic regulation in a test tube. We introduce the rules to design the necessary DNA species, a guide to implement the chemical reactions and ways to optimize the experimental conditions. We finally show how this framework, or "DNA toolbox", can be used to generate an inversion module, though many other behaviors, including oscillators and bistable switches, can be implemented.

Understanding system dynamics of an adaptive enzyme network from globally profiled kinetic parameters

Background

A major challenge in mathematical modeling of biological systems is to determine how model parameters contribute to systems dynamics. As biological processes are often complex in nature, it is desirable to address this issue using a systematic approach. Here, we propose a simple methodology that first performs an enrichment test to find patterns in the values of globally profiled kinetic parameters with which a model can produce the required system dynamics; this is then followed by a statistical test to elucidate the association between individual parameters and different parts of the system’s dynamics.

Results

We demonstrate our methodology on a prototype biological system of perfect adaptation dynamics, namely the chemotaxis model for Escherichia coli. Our results agreed well with those derived from experimental data and theoretical studies in the literature. Using this model system, we showed that there are motifs in kinetic parameters and that these motifs are governed by constraints of the specified system dynamics.

Conclusions

A systematic approach based on enrichment statistical tests has been developed to elucidate the relationships between model parameters and the roles they play in affecting system dynamics of a prototype biological network. The proposed approach is generally applicable and therefore can find wide use in systems biology modeling research.

A conceptual review on systems biology in health and diseases: from biological networks to modern therapeutics

Human physiology is an ensemble of various biological processes spanning from intracellular molecular interactions to the whole body phenotypic response. Systems biology endures to decipher these multi-scale biological networks and bridge the link between genotype to phenotype. The structure and dynamic properties of these networks are responsible for controlling and deciding the phenotypic state of a cell. Several cells and various tissues coordinate together to generate an organ level response which further regulates the ultimate physiological state. The overall network embeds a hierarchical regulatory structure, which when unusually perturbed can lead to undesirable physiological state termed as disease. Here, we treat a disease diagnosis problem analogous to a fault diagnosis problem in engineering systems. Accordingly we review the application of engineering methodologies to address human diseases from systems biological perspective. The review highlights potential networks and modeling approaches used for analyzing human diseases. The application of such analysis is illustrated in the case of cancer and diabetes. We put forth a concept of cell-to-human framework comprising of five modules (data mining, networking, modeling, experimental and validation) for addressing human physiology and diseases based on a paradigm of system level analysis. The review overtly emphasizes on the importance of multi-scale biological networks and subsequent modeling and analysis for drug target identification and designing efficient therapies.

Effects of rate-limiting steps in transcription initiation on genetic filter motifs

The behavior of genetic motifs is determined not only by the gene-gene interactions, but also by the expression patterns of the constituent genes. Live single-molecule measurements have provided evidence that transcription initiation is a sequential process, whose kinetics plays a key role in the dynamics of mRNA and protein numbers. The extent to which it affects the behavior of cellular motifs is unknown. Here, we examine how the kinetics of transcription initiation affects the behavior of motifs performing filtering in amplitude and frequency domain. We find that the performance of each filter is degraded as transcript levels are lowered. This effect can be reduced by having a transcription process with more steps. In addition, we show that the kinetics of the stepwise transcription initiation process affects features such as filter cutoffs. These results constitute an assessment of the range of behaviors of genetic motifs as a function of the kinetics of transcription initiation, and thus will aid in tuning of synthetic motifs to attain specific characteristics without affecting their protein products.

RMOD: a tool for regulatory motif detection in signaling network

Regulatory motifs are patterns of activation and inhibition that appear repeatedly in various signaling networks and that show specific regulatory properties. However, the network structures of regulatory motifs are highly diverse and complex, rendering their identification difficult. Here, we present a RMOD, a web-based system for the identification of regulatory motifs and their properties in signaling networks. RMOD finds various network structures of regulatory motifs by compressing the signaling network and detecting the compressed forms of regulatory motifs. To apply it into a large-scale signaling network, it adopts a new subgraph search algorithm using a novel data structure called path-tree, which is a tree structure composed of isomorphic graphs of query regulatory motifs. This algorithm was evaluated using various sizes of signaling networks generated from the integration of various human signaling pathways and it showed that the speed and scalability of this algorithm outperforms those of other algorithms. RMOD includes interactive analysis and auxiliary tools that make it possible to manipulate the whole processes from building signaling network and query regulatory motifs to analyzing regulatory motifs with graphical illustration and summarized descriptions. As a result, RMOD provides an integrated view of the regulatory motifs and mechanism underlying their regulatory motif activities within the signaling network. RMOD is freely accessible online at the following URL: http://pks.kaist.ac.kr/rmod.

Network motifs provide signatures that characterize metabolism

Motifs are repeating patterns that determine the local properties of networks. In this work, we characterized all 3-node motifs using enzyme commission numbers of the International Union of Biochemistry and Molecular Biology to show that motif abundance is related to biochemical function. Further, we present a comparative analysis of motif distributions in the metabolic networks of 21 species across six kingdoms of life. We found the distribution of motif abundances to be similar between species, but unique across cellular organelles. Finally, we show that motifs are able to capture inter-species differences in metabolic networks and that molecular differences between some biological species are reflected by the distribution of motif abundances in metabolic networks.

The underappreciated role of allostery in the cellular network

Allosteric propagation results in communication between distinct sites in the protein structure; it also encodes specific effects on cellular pathways, and in this way it shapes cellular response. One example of long-range effects is binding of morphogens to cell surface receptors, which initiates a cascade of protein interactions that leads to genome activation and specific cellular action. Allosteric propagation results from combinations of multiple factors, takes place through dynamic shifts of conformational ensembles, and affects the equilibria of macromolecular interactions. Here, we (a) emphasize the well-known yet still underappreciated role of allostery in conveying explicit signals across large multimolecular assemblies and distances to specify cellular action; (b) stress the need for quantitation of the allosteric effects; and finally, (c) propose that each specific combination of allosteric effectors along the pathway spells a distinct function. The challenges are colossal; the inspiring reward will be predicting function, misfunction, and outcomes of drug regimes.

Noise propagation in gene regulation networks involving interlinked positive and negative feedback loops

It is well known that noise is inevitable in gene regulatory networks due to the low-copy numbers of molecules and local environmental fluctuations. The prediction of noise effects is a key issue in ensuring reliable transmission of information. Interlinked positive and negative feedback loops are essential signal transduction motifs in biological networks. Positive feedback loops are generally believed to induce a switch-like behavior, whereas negative feedback loops are thought to suppress noise effects. Here, by using the signal sensitivity (susceptibility) and noise amplification to quantify noise propagation, we analyze an abstract model of the Myc/E2F/MiR-17-92 network that is composed of a coupling between the E2F/Myc positive feedback loop and the E2F/Myc/miR-17-92 negative feedback loop. The role of the feedback loop on noise effects is found to depend on the dynamic properties of the system. When the system is in monostability or bistability with high protein concentrations, noise is consistently suppressed. However, the negative feedback loop reduces this suppression ability (or improves the noise propagation) and enhances signal sensitivity. In the case of excitability, bistability, or monostability, noise is enhanced at low protein concentrations. The negative feedback loop reduces this noise enhancement as well as the signal sensitivity. In all cases, the positive feedback loop acts contrary to the negative feedback loop. We also found that increasing the time scale of the protein module or decreasing the noise autocorrelation time can enhance noise suppression; however, the systems sensitivity remains unchanged. Taken together, our results suggest that the negative/positive feedback mechanisms in coupled feedback loop dynamically buffer noise effects rather than only suppressing or amplifying the noise.

Effects of gene length on the dynamics of gene expression

In Escherichia coli, the nucleotide length of a gene is bound to affect its expression dynamics. From simulations of a stochastic model of gene expression at the nucleotide and codon levels we show that, within realistic parameter values, the nucleotide length affects RNA and protein mean levels, as well as the expected transient time for RNA and protein numbers to change, following a signal. Fluctuations in RNA and protein numbers are found to be minimized for a small range of lengths, which matches the means of the distributions of lengths found in E. coli of both essential and non-essential genes. The variance of the length distribution for essential genes is found to be smaller than for non-essential genes, implying that these distributions are far from random. Finally, gene lengths are shown to affect the kinetics of a genetic switch, namely, the correlation between temporal proteins numbers, the stability of the two noisy attractors of the switch, and how biased is the choice of noisy attractor. The stability increases with gene length due to increased 'memory' about the previous states of the switch. We argue that, by affecting the dynamics of gene expression and of genetic circuits, gene lengths are subject to selection.

Synthetic multicellularity

The ability to synthesize biological constructs on the scale of the organisms we observe unaided is probably one of the more outlandish, yet recurring, dreams humans have had since they began to modify genes. This review brings together recent developments in synthetic biology, cell and developmental biology, computation, and technological development to provide context and direction for the engineering of rudimentary, autonomous multicellular ensembles.

A DNA network as an information processing system

Biomolecular systems that can process information are sought for computational applications, because of their potential for parallelism and miniaturization and because their biocompatibility also makes them suitable for future biomedical applications. DNA has been used to design machines, motors, finite automata, logic gates, reaction networks and logic programs, amongst many other structures and dynamic behaviours. Here we design and program a synthetic DNA network to implement computational paradigms abstracted from cellular regulatory networks. These show information processing properties that are desirable in artificial, engineered molecular systems, including robustness of the output in relation to different sources of variation. We show the results of numerical simulations of the dynamic behaviour of the network and preliminary experimental analysis of its main components.

DYNAMICS OF THE TRANSITION TO PATHOGENICITY INERWINIA CHRYSANTHEMI

The enterobacteria Erwinia chrysanthemi and other soft-rot Erwiniae cause soft-rot disease in plants by secreting extracellular enzymes among which the main virulence factors are pectate lyases (Pels). These pectic enzymes are produced by the activation of the pel genes whose transcription is controlled by a complex regulatory network. Using the knowledge acquired in a previous work, a simplified regulatory network is proposed, keeping only the key variables for the transition to pathogenicity. We identify that the core mechanism for the onset of Pel is governed by a small metabolico-genetic network involving the repressor KdgR and the inductor KDG. Next we consider that the triggering of Pel synthesis is relayed by a quorum sensing (QS) phenomenon describing the ability of bacteria to use the size and density of their colonies to regulate the production of pectate lyases. The simplified network is described by only a few differential equations, thereby allowing the use of standard bifurcation analysis in the phase space. From this modeling emerges a qualitative but generic mechanism for the transition to virulence of a pectinolytic bacterium when it infects a plant.

MicroRNA-mediated positive feedback loop and optimized bistable switch in a cancer network Involving miR-17-92

MicroRNAs (miRNAs) are small, noncoding RNAs that play an important role in many key biological processes, including development, cell differentiation, the cell cycle and apoptosis, as central post-transcriptional regulators of gene expression. Recent studies have shown that miRNAs can act as oncogenes and tumor suppressors depending on the context. The present work focuses on the physiological significance of miRNAs and their role in regulating the switching behavior. We illustrate an abstract model of the Myc/E2F/miR-17-92 network presented by Aguda et al. (2008), which is composed of coupling between the E2F/Myc positive feedback loops and the E2F/Myc/miR-17-92 negative feedback loop. By systematically analyzing the network in close association with plausible experimental parameters, we show that, in the presence of miRNAs, the system bistability emerges from the system, with a bistable switch and a one-way switch presented by Aguda et al. instead of a single one-way switch. Moreover, the miRNAs can optimize the switching process. The model produces a diverse array of response-signal behaviors in response to various potential regulating scenarios. The model predicts that this transition exists, one from cell death or the cancerous phenotype directly to cell quiescence, due to the existence of miRNAs. It was also found that the network involving miR-17-92 exhibits high noise sensitivity due to a positive feedback loop and also maintains resistance to noise from a negative feedback loop.

LsrR quorum sensing "switch" is revealed by a bottom-up approach

Quorum sensing (QS) enables bacterial multicellularity and selective advantage for communicating populations. While genetic “switching” phenomena are a common feature, their mechanistic underpinnings have remained elusive. The interplay between circuit components and their regulation are intertwined and embedded. Observable phenotypes are complex and context dependent. We employed a combination of experimental work and mathematical models to decipher network connectivity and signal transduction in the autoinducer-2 (AI-2) quorum sensing system of E. coli. Negative and positive feedback mechanisms were examined by separating the network architecture into sub-networks. A new unreported negative feedback interaction was hypothesized and tested via a simple mathematical model. Also, the importance of the LsrR regulator and its determinant role in the E. coli QS “switch”, normally masked by interfering regulatory loops, were revealed. Our simple model allowed mechanistic understanding of the interplay among regulatory sub-structures and their contributions to the overall native functioning network. This “bottom up” approach in understanding gene regulation will serve to unravel complex QS network architectures and lead to the directed coordination of emergent behaviors.

Quorum sensing is a mechanism by which bacterial cells communicate within a population. One particular form of communication in E. coli is through a universal signaling molecule known as autoinducer 2. Although the importance of this form of cell-cell interaction has been recognized in the formation of biofilms and virulent infections, the mechanisms by which this form of communication is regulated is still not well understood. In this paper, we presented a method of unraveling these mechanisms by using a combination of experimental work and mathematical models. We took apart the network architecture and isolated the different components. The examination of these isolated sub-networks provided us with a better understanding of the underlying mechanisms that control and regulate bacterial quorum sensing. We were also able to predict new network interactions with the help of our mathematical models. This bottom up approach, combined with our modeling efforts, proved effective in unraveling the mechanisms of quorum sensing in E. coli.

Robust network topologies for generating switch-like cellular responses

Signaling networks that convert graded stimuli into binary, all-or-none cellular responses are critical in processes ranging from cell-cycle control to lineage commitment. To exhaustively enumerate topologies that exhibit this switch-like behavior, we simulated all possible two- and three-component networks on random parameter sets, and assessed the resulting response profiles for both steepness (ultrasensitivity) and extent of memory (bistability). Simulations were used to study purely enzymatic networks, purely transcriptional networks, and hybrid enzymatic/transcriptional networks, and the topologies in each class were rank ordered by parametric robustness (i.e., the percentage of applied parameter sets exhibiting ultrasensitivity or bistability). Results reveal that the distribution of network robustness is highly skewed, with the most robust topologies clustering into a small number of motifs. Hybrid networks are the most robust in generating ultrasensitivity (up to 28%) and bistability (up to 18%); strikingly, a purely transcriptional framework is the most fragile in generating either ultrasensitive (up to 3%) or bistable (up to 1%) responses. The disparity in robustness among the network classes is due in part to zero-order ultrasensitivity, an enzyme-specific phenomenon, which repeatedly emerges as a particularly robust mechanism for generating nonlinearity and can act as a building block for switch-like responses. We also highlight experimentally studied examples of topologies enabling switching behavior, in both native and synthetic systems, that rank highly in our simulations. This unbiased approach for identifying topologies capable of a given response may be useful in discovering new natural motifs and in designing robust synthetic gene networks.

Dynamic modelling of oestrogen signalling and cell fate in breast cancer cells

Cancers of the breast and other tissues arise from aberrant decision-making by cells regarding their survival or death, proliferation or quiescence, damage repair or bypass. These decisions are made by molecular signalling networks that process information from outside and from within the breast cancer cell and initiate responses that determine the cell's survival and reproduction. Because the molecular logic of these circuits is difficult to comprehend by intuitive reasoning alone, we present some preliminary mathematical models of the basic decision circuits in breast cancer cells that may aid our understanding of their susceptibility or resistance to endocrine therapy.