Base editing enables duplex point mutagenesis in Clostridium autoethanogenum at the price of numerous off-target mutations.

Base editors are recent multiplex gene editing tools derived from the Cas9 nuclease of Streptomyces pyogenes. They can target and modify a single nucleotide in the genome without inducing double-strand breaks (DSB) of the DNA helix. As such, they hold great potential for the engineering of microbes that lack effective DSB repair pathways such as homologous recombination (HR) or non-homologous end-joining (NHEJ). However, few applications of base editors have been reported in prokaryotes to date, and their advantages and drawbacks have not been systematically reported. Here, we used the base editors Target-AID and Target-AID-NG to introduce nonsense mutations into four different coding sequences of the industrially relevant Gram-positive bacterium Clostridium autoethanogenum. While up to two loci could be edited simultaneously using a variety of multiplexing strategies, most colonies exhibited mixed genotypes and most available protospacers led to undesired mutations within the targeted editing window. Additionally, fifteen off-target mutations were detected by sequencing the genome of the resulting strain, among them seven single-nucleotide polymorphisms (SNP) in or near loci bearing some similarity with the targeted protospacers, one 15 nt duplication, and one 12 kb deletion which removed uracil DNA glycosylase (UDG), a key DNA repair enzyme thought to be an obstacle to base editing mutagenesis. A strategy to process prokaryotic single-guide RNA arrays by exploiting tRNA maturation mechanisms is also illustrated.

A genome-scale metabolic model of Cupriavidus necator H16 integrated with TraDIS and transcriptomic data reveals metabolic insights for biotechnological applications.

Exploiting biological processes to recycle renewable carbon into high value platform chemicals provides a sustainable and greener alternative to current reliance on petrochemicals. In this regard Cupriavidus necator H16 represents a particularly promising microbial chassis due to its ability to grow on a wide range of low-cost feedstocks, including the waste gas carbon dioxide, whilst also naturally producing large quantities of polyhydroxybutyrate (PHB) during nutrient-limited conditions. Understanding the complex metabolic behaviour of this bacterium is a prerequisite for the design of successful engineering strategies for optimising product yields. We present a genome-scale metabolic model (GSM) of C. necator H16 (denoted iCN1361), which is directly constructed from the BioCyc database to improve the readability and reusability of the model. After the initial automated construction, we have performed extensive curation and both theoretical and experimental validation. By carrying out a genome-wide essentiality screening using a Transposon-directed Insertion site Sequencing (TraDIS) approach, we showed that the model could predict gene knockout phenotypes with a high level of accuracy. Importantly, we indicate how experimental and computational predictions can be used to improve model structure and, thus, model accuracy as well as to evaluate potential false positives identified in the experiments. Finally, by integrating transcriptomics data with iCN1361 we create a condition-specific model, which, importantly, better reflects PHB production in C. necator H16. Observed changes in the omics data and in-silico-estimated alterations in fluxes were then used to predict the regulatory control of key cellular processes. The results presented demonstrate that iCN1361 is a valuable tool for unravelling the system-level metabolic behaviour of C. necator H16 and can provide useful insights for designing metabolic engineering strategies.

Engineering improved ethylene production: Leveraging systems biology and adaptive laboratory evolution.

Ethylene is a small hydrocarbon gas widely used in the chemical industry. Annual worldwide production currently exceeds 150 million tons, producing considerable amounts of CO2 contributing to climate change. The need for a sustainable alternative is therefore imperative. Ethylene is natively produced by several different microorganisms, including Pseudomonas syringae pv. phaseolicola via a process catalyzed by the ethylene-forming enzyme (EFE), subsequent heterologous expression of EFE has led to ethylene production in non-native bacterial hosts including Escherichia coli and cyanobacteria. However, solubility of EFE and substrate availability remain rate-limiting steps in biological ethylene production. We employed a combination of genome-scale metabolic modelling, continuous fermentation, and protein evolution to enable the accelerated development of a high efficiency ethylene producing E. coli strain, yielding a 49-fold increase in production, the most significant improvement reported to date. Furthermore, we have clearly demonstrated that this increased yield resulted from metabolic adaptations that were uniquely linked to EFE (wild type versus mutant). Our findings provide a novel solution to deregulate metabolic bottlenecks in key pathways, which can be readily applied to address other engineering challenges.

Quantitative Bioreactor Monitoring of Intracellular Bacterial Metabolites in Clostridium autoethanogenum Using Liquid Chromatography-Isotope Dilution Mass Spectrometry.

We report a liquid chromatography-isotope dilution mass spectrometry method for the simultaneous quantification of 131 intracellular bacterial metabolites of Clostridium autoethanogenum. A comprehensive mixture of uniformly 13C-labeled internal standards (U-13C IS) was biosynthesized from the closely related bacterium Clostridium pasteurianum using 4% 13C-glucose as a carbon source. The U-13C IS mixture combined with 12C authentic standards was used to validate the linearity, precision, accuracy, repeatability, limits of detection, and quantification for each metabolite. A robust-fitting algorithm was employed to reduce the weight of the outliers on the quantification data. The metabolite calibration curves were linear with R 2 ≥ 0.99, limits of detection were ≤1.0 µM, limits of quantification were ≤10 µM, and precision/accuracy was within RSDs of 15% for all metabolites. The method was subsequently applied for the daily monitoring of the intracellular metabolites of C. autoethanogenum during a CO gas fermentation over 40 days as part of a study to optimize biofuel production. The concentrations of the metabolites were estimated at steady states of different pH levels using the robust-fitting mathematical approach, and we demonstrate improved accuracy of results compared to conventional regression. Metabolic pathway analysis showed that reactions of the incomplete (branched) tricarboxylic acid "cycle" were the most affected pathways associated with the pH shift in the bioreactor fermentation of C. autoethanogenum and the concomitant changes in ethanol production.

Physicochemical and metabolic constraints for thermodynamics-based stoichiometric modelling under mesophilic growth conditions.

Metabolic engineering in the post-genomic era is characterised by the development of new methods for metabolomics and fluxomics, supported by the integration of genetic engineering tools and mathematical modelling. Particularly, constraint-based stoichiometric models have been widely studied: (i) flux balance analysis (FBA) (in silico), and (ii) metabolic flux analysis (MFA) (in vivo). Recent studies have enabled the incorporation of thermodynamics and metabolomics data to improve the predictive capabilities of these approaches. However, an in-depth comparison and evaluation of these methods is lacking. This study presents a thorough analysis of two different in silico methods tested against experimental data (metabolomics and 13C-MFA) for the mesophile Escherichia coli. In particular, a modified version of the recently published matTFA toolbox was created, providing a broader range of physicochemical parameters. Validating against experimental data allowed the determination of the best physicochemical parameters to perform the TFA (Thermodynamics-based Flux Analysis). An analysis of flux pattern changes in the central carbon metabolism between 13C-MFA and TFA highlighted the limited capabilities of both approaches for elucidating the anaplerotic fluxes. In addition, a method based on centrality measures was suggested to identify important metabolites that (if quantified) would allow to further constrain the TFA. Finally, this study emphasised the need for standardisation in the fluxomics community: novel approaches are frequently released but a thorough comparison with currently accepted methods is not always performed.

Gsmodutils: a python based framework for test-driven genome scale metabolic model development.

MOTIVATION: Genome scale metabolic models (GSMMs) are increasingly important for systems biology and metabolic engineering research as they are capable of simulating complex steady-state behaviour. Constraints based models of this form can include thousands of reactions and metabolites, with many crucial pathways that only become activated in specific simulation settings. However, despite their widespread use, power and the availability of tools to aid with the construction and analysis of large scale models, little methodology is suggested for their continued management. For example, when genome annotations are updated or new understanding regarding behaviour is discovered, models often need to be altered to reflect this. This is quickly becoming an issue for industrial systems and synthetic biotechnology applications, which require good quality reusable models integral to the design, build, test and learn cycle. RESULTS: As part of an ongoing effort to improve genome scale metabolic analysis, we have developed a test-driven development methodology for the continuous integration of validation data from different sources. Contributing to the open source technology based around COBRApy, we have developed the gsmodutils modelling framework placing an emphasis on test-driven design of models through defined test cases. Crucially, different conditions are configurable allowing users to examine how different designs or curation impact a wide range of system behaviours, minimizing error between model versions. AVAILABILITY AND IMPLEMENTATION: The software framework described within this paper is open source and freely available from http://github.com/SBRCNottingham/gsmodutils. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Progress towards platform chemical production using Clostridium autoethanogenum.

Since 2013, there has been an explosion in the number of research articles published on Clostridium autoethanogenum, an acetogen capable of producing platform chemicals such as ethanol and 2,3-butanediol from greenhouse gases. However, no review focusing solely on C. autoethanogenum has appeared in the literature. This review outlines the research conducted into this organism in three broad categories (Enzymology, Genetics, and Systems Biology) and suggestions for future research are offered.

Quantitative Isotope-Dilution High-Resolution-Mass-Spectrometry Analysis of Multiple Intracellular Metabolites in Clostridium autoethanogenum with Uniformly 13C-Labeled Standards Derived from Spirulina.

We have investigated the applicability of commercially available lyophilized spirulina ( Arthrospira platensis), a microorganism uniformly labeled with 13C, as a readily accessible source of multiple 13C-labeled metabolites suitable as internal standards for the quantitative determination of intracellular bacterial metabolites. Metabolites of interest were analyzed by hydrophilic-interaction liquid chromatography coupled with high-resolution mass spectrometry. Multiple internal standards obtained from uniformly (U)-13C-labeled extracts from spirulina were used to enable isotope-dilution mass spectrometry (IDMS) in the identification and quantification of intracellular metabolites. Extraction of the intracellular metabolites of Clostridium autoethanogenum using 2:1:1 chloroform/methanol/water was found to be the optimal method in comparison with freeze-thaw, homogenization, and sonication methods. The limits of quantification were ≤1 µM with excellent linearity for all of the calibration curves ( R2 ≥ 0.99) for 74 metabolites. The precision and accuracy were found to be within relative standard deviations (RSDs) of 15% for 49 of the metabolites and within RSDs of 20% for all of the metabolites. The method was applied to study the effects of feeding different levels of carbon monoxide (as a carbon source) on the central metabolism and Wood-Ljungdahl pathway of C. autoethanogenum grown in continuous culture over 35 days. Using LC-IDMS with U-13C spirulina allowed the successful quantification of 52 metabolites in the samples, including amino acids, carboxylic acids, sugar phosphates, purines, and pyrimidines. The method provided absolute quantitative data on intracellular metabolites that was suitable for computational modeling to understand and optimize the C. autoethanogenum metabolic pathways active in gas fermentation.

Mathematical modelling of clostridial acetone-butanol-ethanol fermentation.

Clostridial acetone-butanol-ethanol (ABE) fermentation features a remarkable shift in the cellular metabolic activity from acid formation, acidogenesis, to the production of industrial-relevant solvents, solventogensis. In recent decades, mathematical models have been employed to elucidate the complex interlinked regulation and conditions that determine these two distinct metabolic states and govern the transition between them. In this review, we discuss these models with a focus on the mechanisms controlling intra- and extracellular changes between acidogenesis and solventogenesis. In particular, we critically evaluate underlying model assumptions and predictions in the light of current experimental knowledge. Towards this end, we briefly introduce key ideas and assumptions applied in the discussed modelling approaches, but waive a comprehensive mathematical presentation. We distinguish between structural and dynamical models, which will be discussed in their chronological order to illustrate how new biological information facilitates the 'evolution' of mathematical models. Mathematical models and their analysis have significantly contributed to our knowledge of ABE fermentation and the underlying regulatory network which spans all levels of biological organization. However, the ties between the different levels of cellular regulation are not well understood. Furthermore, contradictory experimental and theoretical results challenge our current notion of ABE metabolic network structure. Thus, clostridial ABE fermentation still poses theoretical as well as experimental challenges which are best approached in close collaboration between modellers and experimentalists.

Whole genome sequence and manual annotation of Clostridium autoethanogenum, an industrially relevant bacterium.

BACKGROUND: Clostridium autoethanogenum is an acetogenic bacterium capable of producing high value commodity chemicals and biofuels from the C1 gases present in synthesis gas. This common industrial waste gas can act as the sole energy and carbon source for the bacterium that converts the low value gaseous components into cellular building blocks and industrially relevant products via the action of the reductive acetyl-CoA (Wood-Ljungdahl) pathway. Current research efforts are focused on the enhancement and extension of product formation in this organism via synthetic biology approaches. However, crucial to metabolic modelling and directed pathway engineering is a reliable and comprehensively annotated genome sequence. RESULTS: We performed next generation sequencing using Illumina MiSeq technology on the DSM10061 strain of Clostridium autoethanogenum and observed 243 single nucleotide discrepancies when compared to the published finished sequence (NCBI: GCA_000484505.1), with 59.1 % present in coding regions. These variations were confirmed by Sanger sequencing and subsequent analysis suggested that the discrepancies were sequencing errors in the published genome not true single nucleotide polymorphisms. This was corroborated by the observation that over 90 % occurred within homopolymer regions of greater than 4 nucleotides in length. It was also observed that many genes containing these sequencing errors were annotated in the published closed genome as encoding proteins containing frameshift mutations (18 instances) or were annotated despite the coding frame containing stop codons, which if genuine, would severely hinder the organism's ability to survive. Furthermore, we have completed a comprehensive manual curation to reduce errors in the annotation that occur through serial use of automated annotation pipelines in related species. As a result, different functions were assigned to gene products or previous functional annotations rejected because of missing evidence in various occasions. CONCLUSIONS: We present a revised manually curated full genome sequence for Clostridium autoethanogenum DSM10061, which provides reliable information for genome-scale models that rely heavily on the accuracy of annotation, and represents an important step towards the manipulation and metabolic modelling of this industrially relevant acetogen.

Coenzyme A-transferase-independent butyrate re-assimilation in Clostridium acetobutylicum-evidence from a mathematical model.

The hetero-dimeric CoA-transferase CtfA/B is believed to be crucial for the metabolic transition from acidogenesis to solventogenesis in Clostridium acetobutylicum as part of the industrial-relevant acetone-butanol-ethanol (ABE) fermentation. Here, the enzyme is assumed to mediate re-assimilation of acetate and butyrate during a pH-induced metabolic shift and to faciliate the first step of acetone formation from acetoacetyl-CoA. However, recent investigations using phosphate-limited continuous cultures have questioned this common dogma. To address the emerging experimental discrepancies, we investigated the mutant strain Cac-ctfA398s::CT using chemostat cultures. As a consequence of this mutation, the cells are unable to express functional ctfA and are thus lacking CoA-transferase activity. A mathematical model of the pH-induced metabolic shift, which was recently developed for the wild type, is used to analyse the observed behaviour of the mutant strain with a focus on re-assimilation activities for the two produced acids. Our theoretical analysis reveals that the ctfA mutant still re-assimilates butyrate, but not acetate. Based upon this finding, we conclude that C. acetobutylicum possesses a CoA-tranferase-independent butyrate uptake mechanism that is activated by decreasing pH levels. Furthermore, we observe that butanol formation is not inhibited under our experimental conditions, as suggested by previous batch culture experiments. In concordance with recent batch experiments, acetone formation is abolished in chemostat cultures using the ctfa mutant.

A shift in the dominant phenotype governs the pH-induced metabolic switch of Clostridium acetobutylicumin phosphate-limited continuous cultures.

In response to changing extracellular pH levels, phosphate-limited continuous cultures of Clostridium acetobutylicum reversibly switches its metabolism from the dominant formation of acids to the prevalent production of solvents. Previous experimental and theoretical studies have revealed that this pH-induced metabolic switch involves a rearrangement of the intracellular transcriptomic, proteomic and metabolomic composition of the clostridial cells. However, the influence of the population dynamics on the observations reported has so far been neglected. Here, we present a method for linking the pH shift, clostridial growth and the acetone-butanol-ethanol fermentation metabolic network systematically into a model which combines the dynamics of the external pH and optical density with a metabolic model. Furthermore, the recently found antagonistic expression pattern of the aldehyde/alcohol dehydrogenases AdhE1/2 and pH-dependent enzyme activities have been included into this combined model. Our model predictions reveal that the pH-induced metabolic shift under these experimental conditions is governed by a phenotypic switch of predominantly acidogenic subpopulation towards a predominantly solventogenic subpopulation. This model-driven explanation of the pH-induced shift from acidogenesis to solventogenesis by population dynamics casts an entirely new light on the clostridial response to changing pH levels. Moreover, the results presented here underline that pH-dependent growth and pH-dependent specific enzymatic activity play a crucial role in this adaptation. In particular, the behaviour of AdhE1 and AdhE2 seems to be the key factor for the product formation of the two phenotypes, their pH-dependent growth, and thus, the pH-induced metabolic switch in C. acetobutylicum.

Simulations of stressosome activation emphasize allosteric interactions between RsbR and RsbT.

BACKGROUND: The stressosome is a bacterial signalling complex that responds to environmental changes by initiating a protein partner switching cascade, which leads to the release of the alternative sigma factor, σB. Stress perception increases the phosphorylation of the stressosome sensor protein, RsbR, and the scaffold protein, RsbS, by the protein kinase, RsbT. Subsequent dissociation of RsbT from the stressosome activates the σB cascade. However, the sequence of physical events that occur in the stressosome during signal transduction is insufficiently understood. RESULTS: Here, we use computational modelling to correlate the structure of the stressosome with the efficiency of the phosphorylation reactions that occur upon activation by stress. In our model, the phosphorylation of any stressosome protein is dependent upon its nearest neighbours and their phosphorylation status. We compare different hypotheses about stressosome activation and find that only the model representing the allosteric activation of the kinase RsbT, by phosphorylated RsbR, qualitatively reproduces the experimental data. CONCLUSIONS: Our simulations and the associated analysis of published data support the following hypotheses: (i) a simple Boolean model is capable of reproducing stressosome dynamics, (ii) different stressors induce identical stressosome activation patterns, and we also confirm that (i) phosphorylated RsbR activates RsbT, and (ii) the main purpose of RsbX is to dephosphorylate RsbS-P.

Integrative modelling of pH-dependent enzyme activity and transcriptomic regulation of the acetone-butanol-ethanol fermentation of Clostridium acetobutylicum in continuous culture.

In a continuous culture under phosphate limitation the metabolism of Clostridium acetobutylicum depends on the external pH level. By comparing seven steady-state conditions between pH 5.7 and pH 4.5 we show that the switch from acidogenesis to solventogenesis occurs between pH 5.3 and pH 5.0 with an intermediate state at pH 5.1. Here, an integrative study is presented investigating how a changing external pH level affects the clostridial acetone-butanol-ethanol (ABE) fermentation pathway. This is of particular interest as the biotechnological production of n-butanol as biofuel has recently returned into the focus of industrial applications. One prerequisite is the furthering of the knowledge of the factors determining the solvent production and their integrative regulations. We have mathematically analysed the influence of pH-dependent specific enzyme activities of branch points of the metabolism on the product formation. This kinetic regulation was compared with transcriptomic regulation regarding gene transcription and the proteomic profile. Furthermore, both regulatory mechanisms were combined yielding a detailed projection of their individual and joint effects on the product formation. The resulting model represents an important platform for future developments of industrial butanol production based on C. acetobutylicum.

Proteolysis of beta-galactosidase following SigmaB activation in Bacillus subtilis.

In Bacillus subtilis the σ(B) mediated general stress response provides protection against various environmental and energy related stress conditions. To better understand the general stress response, we need to explore the mechanism by which the components interact. Here, we performed experiments in B. subtilis wild type and mutant strains to test and validate a mathematical model of the dynamics of σ(B) activity. In the mutant strain BSA115, σ(B) transcription is inducible by the addition of IPTG and negative control of σ(B) activity by the anti-sigma factor RsbW is absent. In contrast to our expectations of a continuous ß-galactosidase activity from a ctc::lacZ fusion, we observed a transient activity in the mutant. To explain this experimental finding, we constructed mathematical models reflecting different hypotheses regarding the regulation of σ(B) and ß-galactosidase dynamics. Only the model assuming instability of either ctc::lacZ mRNA or ß-galactosidase protein is able to reproduce the experiments in silico. Subsequent Northern blot experiments revealed stable high-level ctc::lacZ mRNA concentrations after the induction of the σ(B) response. Therefore, we conclude that protein instability following σ(B) activation is the most likely explanation for the experimental observations. Our results thus support the idea that B. subtilis increases the cytoplasmic proteolytic degradation to adapt the proteome in face of environmental challenges following activation of the general stress response. The findings also have practical implications for the analysis of stress response dynamics using lacZ reporter gene fusions, a frequently used strategy for the σ(B) response.

A systems biology approach to investigate the effect of pH-induced gene regulation on solvent production by Clostridium acetobutylicum in continuous culture.

BACKGROUND: Clostridium acetobutylicum is an anaerobic bacterium which is known for its solvent-producing capabilities, namely regarding the bulk chemicals acetone and butanol, the latter being a highly efficient biofuel. For butanol production by C. acetobutylicum to be optimized and exploited on an industrial scale, the effect of pH-induced gene regulation on solvent production by C. acetobutylicum in continuous culture must be understood as fully as possible. RESULTS: We present an ordinary differential equation model combining the metabolic network governing solvent production with regulation at the genetic level of the enzymes required for this process. Parameterizing the model with experimental data from continuous culture, we demonstrate the influence of pH upon fermentation products: at high pH (pH 5.7) acids are the dominant product while at low pH (pH 4.5) this switches to solvents. Through steady-state analyses of the model we focus our investigations on how alteration in gene expression of C. acetobutylicum could be exploited to increase butanol yield in a continuous culture fermentation. CONCLUSIONS: Incorporating gene regulation into the model of solvent production by C. acetobutylicum enables an accurate representation of the pH-induced switch to solvent production to be obtained and theoretical investigations of possible synthetic-biology approaches to be pursued. Steady-state analyses suggest that, to increase butanol yield, alterations in the expression of single solvent-associated genes are insufficient; a more complex approach targeting two or more genes is required.

Modeling of cellular processes: methods, data, and requirements.

Systems biology is a comprehensive quantitative analysis how the components of a biological system interact over time which requires an interdisciplinary team of investigators. System-theoretic methods are applied to investigate the system's behavior. Using known information about the considered system, a conceptual model is defined. It is transferred in a mathematical model that can be simulated (analytically or numerically) and analyzed using system-theoretic tools. Finally, simulation results are compared with experimental data. However, assumptions, approximations, and requirements to available experimental data are crucial ingredients of this systems biology workflow. Consequently, the modeling of cellular processes creates special demands on the design of experiments: the quality, the amount, and the completeness of data. The relation between models and data is discussed in this chapter. Thereby, we focus on the requirements on experimental data from the perspective of systems biology projects.

A mathematical analysis of nuclear intensity dynamics for Mig1-GFP under consideration of bleaching effects and background noise in Saccharomyces cerevisiae.

Fluorescence microscopy is an imaging technique that provides insights into signal transduction pathways through the generation of quantitative data, such as the spatiotemporal distribution of GFP-tagged proteins in signaling pathways. The data acquired are, however, usually a composition of both the GFP-tagged proteins of interest and of an autofluorescent background, which both undergo photobleaching during imaging. We here present a mathematical model based on ordinary differential equations that successfully describes the shuttling of intracellular Mig1-GFP under changing environmental conditions regarding glucose concentration. Our analysis separates the different bleaching rates of Mig1-GFP and background, and the background-to-Mig1-GFP ratio. By applying our model to experimental data, we can thus extract the Mig1-GFP signal from the overall acquired signal and investigate the influence of kinase and phosphatase on Mig1. We found a stronger regulation of Mig1 through its kinase than through its phosphatase when controlled by the glucose concentration, with a constant (de)phosphorylation rate independent of the glucose concentration. By replacing the term for decreasing excited Mig1-GFP concentration with a constant, we were able to reconstruct the dynamics of Mig1-GFP, as it would occur without bleaching and background noise. Our model effectively demonstrates how data, acquired with an optical microscope, can be processed and used for a systems biology analysis of signal transduction pathways.

How mathematical modelling elucidates signalling in Bacillus subtilis.

Appropriate stimulus perception, signal processing and transduction ensure optimal adaptation of bacteria to environmental challenges. In the Gram-positive model bacterium Bacillus subtilis signalling networks and molecular interactions therein are well-studied, making this species a suitable candidate for the application of mathematical modelling. Here, we review systems biology approaches, focusing on chemotaxis, sporulation, σ(B) -dependent general stress response and competence. Processes like chemotaxis and Z-ring assembly depend critically on the subcellular localization of proteins. Environmental response strategies, including sporulation and competence, are characterized by phenotypic heterogeneity in isogenic cultures. The examples of mathematical modelling also include investigations that have demonstrated how operon structure and signalling dynamics are intricately interwoven to establish optimal responses. Our review illustrates that these interdisciplinary approaches offer new insights into the response of B. subtilis to environmental challenges. These case studies reveal modelling as a tool to increase the understanding of complex systems, to help formulating hypotheses and to guide the design of more directed experiments that test predictions.

Dynamics of receptor and protein transducer homodimerisation.

BACKGROUND: Signalling pathways are complex systems in which not only simple monomeric molecules interact, but also more complex structures that include constitutive or induced protein assemblies. In particular, the hetero-and homo-dimerisation of proteins is a commonly encountered motif in signalling pathways. Several authors have suggested in recent times that dimerisation relates to a series of physical and biological outcomes used by the cell in the regulation of signal transduction. RESULTS: In this paper we investigate the role of homodimerisation in receptor-protein transducer interactions. Towards this end, mathematical modelling is used to analyse the features of such kind of interactions and to predict the behaviour of the system under different experimental conditions. A kinetic model in which the interaction between homodimers provokes a dual mechanism of activation (single and double protein transducer activation at the same time) is proposed. In addition, we analyse under which conditions the use of a power-law representation for the system is useful. Furthermore, we investigate the dynamical consequences of this dual mechanism and compare the performance of the system in different simulated experimental conditions. CONCLUSION: The analysis of our mathematical model suggests that in receptor-protein interacting systems with dual mechanism there may be a shift between double and single activation in a way that intense double protein transducer activation could initiate and dominate the signal in the short term (getting a fast intense signal), while single protein activation could control the system in the medium and long term (when input signal is weaker and decreases slowly). Our investigation suggests that homodimerisation and oligomerisation are mechanisms used to enhance and regulate the dynamic properties of the initial steps in signalling pathways.