Computational modelling and analysis of the molecular network regulating sporulation initiation in Bacillus subtilis
© Ihekwaba et al.; licensee BioMed Central Ltd. 2014
Received: 8 August 2014
Accepted: 13 October 2014
Published: 24 October 2014
Bacterial spores are important contaminants in food, and the spore forming bacteria are often implicated in food safety and food quality considerations. Spore formation is a complex developmental process involving the expression of more than 500 genes over the course of 6 to 8 hrs. The process culminates in the formation of resting cells capable of resisting environmental extremes and remaining dormant for long periods of time, germinating when conditions promote further vegetative growth. Experimental observations of sporulation and germination are problematic and time consuming so that reliable models are an invaluable asset in terms of prediction and risk assessment. In this report we develop a model which assists in the interpretation of sporulation dynamics.
This paper defines and analyses a mathematical model for the network regulating Bacillus subtilis sporulation initiation, from sensing of sporulation signals down to the activation of the early genes under control of the master regulator Spo0A. Our model summarises and extends other published modelling studies, by allowing the user to execute sporulation initiation in a scenario where Isopropyl ?-D-1-thiogalactopyranoside (IPTG) is used as an artificial sporulation initiator as well as in modelling the induction of sporulation in wild-type cells. The analysis of the model results and the comparison with experimental data indicate that the model is good at predicting inducible responses to sporulation signals. However, the model is unable to reproduce experimentally observed accumulation of phosphorelay sporulation proteins in wild type B. subtilis. This model also highlights that the phosphorelay sub-component, which relays the signals detected by the sensor kinases to the master regulator Spo0A, is crucial in determining the response dynamics of the system.
We show that there is a complex connectivity between the phosphorelay features and the master regulatory Spo0A. Additional we discovered that the experimentally observed regulation of the phosphotransferase Spo0B for wild-type B. subtilis may be playing an important role in the network which suggests that modelling of sporulation initiation may require additional experimental support.
Spore-forming bacteria are a major cause of food spoilage and disease, and food industries actively employ strategies to ensure adequate inactivation of spores and the control of outgrowth, both for species that potentially lead to spoilage and for foodborne pathogens. The contamination of foods with bacterial spores is well documented in many scenarios ,. Spore formers need specific conditions for germination and growth, which can occur in food, such as low or high temperatures and anaerobic or acidophilic conditions ,, sometimes in combination.
Most microbial spore forming bacteria respond to stress (e.g., nutrient deprivation) by inducing the expression of an appropriate suit of adaptive (stress-response) genes to help them cope with adverse environmental circumstances; an extreme example is endospore formation. The initiation of sporulation is one of the decisive moments in the life cycle of sporulating bacteria, whereby an extremely durable cell called a spore is formed, as exemplified by the bacterium Bacillus subtilis.
The decision to abandon vegetative growth and enter sporulation has far reaching consequences for B. subtilis. It involves a switch between two completely different genetic programs requiring energy-expensive changes in gene expression and cell morphology ,. These changes are regulated by a complex network involving more than 120 genes ,. The size and complexity of the sporulation network make it difficult to interpret the switch from vegetative growth to sporulation in terms of the interactions between genes, proteins, and small molecules. This is where computer modelling and simulation is able to help. Mathematical modelling is utilized in understanding how a global response can emerge from a network of local interactions -. In addition, the use of modelling and simulation tools may permit further, the formulation of hypotheses on missing components and interactions which, after a process of reiterative computer simulation, can guide future experimentation.
Various mathematical models of B. subtilis sporulation mechanisms can be found in the literature ,-, however, for the purpose of our study we only mention a small list of them here. We start with the work of Narula and colleagues, who in  defined a computational model which was used to generate a molecular explanation of how the cell integrates and relays its environmental signals to the master regulator of the sporulation initiation process, Spo0A. They employed an inducible promoter system whereby the activation of KinA, the upstream phosphoryl donor in the sporulation network, was regulated through the addition of IPTG to the medium. Moreover, in their more recent study , they further considered a doubly inducible system whereby another kinase, KinC, was expressed from an inducible IPTG promoter, while Spo0A was expressed from a xylose-inducible promoter. Based on the experimental analysis of these systems, they were able to define a mechanistic model supporting a conjecture on how important the proper activation pattern of Spo0A was to the system and how its early accumulation may impair the ability of the cell to complete sporulation. Other groups, Sen and colleagues, as well as Kothamachu and colleagues performed also modelling studies, however, their main focus was on the four component B. subtilis phosphorelay, the key module in the molecular network that transduces the environmental signals captured by the sensing kinases towards the master regulator Spo0A ,. While Sen and colleagues  considered the complex positive and negative feedback loops thought to exist between Spo0A and the upstream proteins of the phosphorelay; as well as how the modulation induced by B. subtilis cell cycle on the synthesis of phosphorelay proteins affects the dynamics of Spo0A; Kothamachu et. al., in their work  modelled the interactions between the four components of the phosphorelay to study the dependence of the Spo0A steady state activation level on the kinetic parameters of the phospho-transfer reactions - without taking into consideration the feedback loops among species. Another group, Kuchina and colleagues , took a different perspective in their modelling, by trying to solve the problem of whether the molecular machinery regulating the sporulation process in B. subtilis proceeds through a sequence of reversible or irreversible decision-making steps. To do this, the authors employed a different modelling approach, whereby the impact of fluctuating environmental conditions was modelled with the aim of identifying the exact point at which cells become committed to sporulation - abstracting the molecular details of the sporulation network.
Since it is known that B subtilis exhibits ultra-sensitivity to sporulation stimuli , characterizing the step at which sporulation decision is made could provide hints on the upstream molecular dynamics. The experimental data provided by  has shown that the best matching scenario is one when the initial progression towards sporulation commitment is reversible (with an irreversible later stage), with commitment arising downstream of Spo0A activation. This finding would therefore suggest that the switch mechanism leading to the ultra-sensitivity is to be placed downstream of Spo0A activation, in the area of the sporulation network which should lead to the sustained expression of the sigma factors that control later sporulation genes.
After reviewing the aforementioned works, the purpose of our study was to understand how the current state-of-the-art in the modelling of B. subtilis sporulation - which so far has been based on the use of inducible systems - can be used to explain the behaviour of the organism, including wild type strains. To this end, we first defined a sporulation model that borrowed generously information from the above-mentioned studies and summarised the current knowledge of the B. subtilis sporulation network at the molecular level. The scope of our model is limited to the chain of molecular interactions that are triggered by the inducers, which are then forwarded through the phosphorelay, ultimately leading to the activation of the master regulator Spo0A, and to the expression of the spoIIA, spoIIE and spoIIG genes. This model, which is implemented in the modelling and simulation package COPASI, is available for download in the COPASI native format (Additional file 1) and in SBML as well (Additional file 2), to facilitate checks of all the results presented in this paper. We analyzed the model output and performed a validation against published experimental data. We use transient, steady state and parameter sensitivity analysis to determine, among model components and parameters, the most important quantities in determining the predicted behaviour.
The results of our analysis suggest that the response to artificial sporulation signals in inducible systems is quantitatively and qualitatively different from the one observed in wild type and cannot be predicted with models of current biological knowledge. Specifically, in wild type a slower and almost linear pattern of Spo0A activation is observed, in contrast with the switch-like behaviour that is observed for the inducible models. This conclusion is supported by a number of remarks that can be found dispersed in the literature, which point to the need for a gradual accumulation of active Spo0A (see for instance  and ) or the influence that the cell-cycle has on slowing-down Spo0A accumulation (see  and ). Moreover, our analysis highlights components and interactions in the phosphorelay that require a further characterization. Indeed, some experimental evidence of regulated synthesis for the phosphorelay component Spo0B was not included in the published models, and the reaction kinetics used in our model, which are largely extracted from the literature, may not provide a faithful representation of the molecular mechanism for the wild type behaviour.
In this section, we provide a detailed description of the mathematical model as well as the methods used for analysis.
B. subtilis sporulation network
The DNA-binding protein Spo0A is the master regulator for entry into sporulation. The capacity of the protein to alter transcription is governed both by accumulation of the protein and its state of phosphorylation. Phosphorylation of Spo0A causes the re-orientation of a phenylalanine residue in the molecule, allowing an effector domain to become active .
Spo0A phosphorylation is controlled by a multicomponent phosphorelay capable of integrating diverse physiological and environmental signals. The B. subtilis phosphorelay is comprised of at least five sensor histidine kinases - termed KinA, KinB, KinC, KinD, and KinE - a secondary messenger, Spo0F, a phosphotransferase, Spo0B, and the response regulator Spo0A. The presence of multiple phosphorelay components in this organism is thought to provide more points of regulation and increased precision for the response as each step reinforces discrimination . Spo0A is both an activator and a repressor of gene transcription. As an activator, Spo0A acts in conjunction with RNA polymerase containing either the housekeeping sigma factors sigma A or the alternative sigma factor, sigmaH, by binding to cognate sequences in the promoter regions of genes under its control . Binding of Spo0A~P to DNA directly affects the expression of approximately 121 genes . Among these are key genes that drive the positive regulation of sporulation - particularly the spoIIA,, spoIIE and spoIIG, genes involved in establishing compartment-specific transcription under the control of σF (spoIIA operon and the spoIIE gene) and σE (spoIIG operon) . Acting as a repressor, Spo0A blocks expression of the abr B gene, a gene which encodes the ‘transition state’ regulator, AbrB, which is itself a repressor of sigmaH ,. Spo0A mediated repression of abr B leads to the depletion of the AbrB protein from the cell and to the accumulation of sigmaH, which results in the activation of genes under its positive control - which includes KinA, Spo0F and the Spo0A gene itself ,. Therefore, Spo0A-mediated repression of the abr B sets up a self-reinforcing cycle that contributes to the accumulation of Spo0A protein at the start of the sporulation ,.
Brief summary of the model
Out of the five kinases identified as capable of initiating sporulation in B. subtilis, we have considered KinA, the major kinase responsible for initiation of sporulation in our model. KinA overexpression during exponential growth is sufficient to induce entry into sporulation . Observations on how specific sporulation-inducing signals affect the activation of this protein are available in literature ,.
We consider two distinct types of induction signals that drive the cell into sporulation. One is a sporulation signal (SS, in the model), which we use to represent the influence of the adverse environment on the auto-phosphorylation of KinA (the molecular details of SS are still unknown). The other signal is the commonly used artificial sporulation initiation system where KinA expression level is controlled from an isopropyl-D-thiogalactopyranoside (IPTG) inducible promoter. IPTG is a popular chemical inducer that can be used to exogenously control the transcriptional activity of wild-type or synthetic lactose repressor (LacI) in order to manipulate the production of KinA proteins in bacteria . Use of IPTG allows for tuneable control of KinA expression ,. Hence, we can feed our model with different stimuli that together can represent either an artificial induction, independent of the environmental conditions, or a natural induction occurring in wild type.
The sporulation signal we introduce is modelled as an abstract species SS, which controls the autophosphorylation of KinA dimers -. The IPTG regulation effect on the tuning of the KinA levels has been modelled by following some of the assumptions available in literature ,,,. We model IPTG regulation of sporulation initiation through the indirect release of KinA synthesis inhibition by LacI, which when successfully incorporated into a promoter, has been used to regulate the expression of many genes. LacI is a DNA-binding transcription factor that represses transcription of the operon involved in transport and catabolism of lactose ,. This binding is altered by the presence of sugar ligands (inducers), which elicit a conformational change in the protein to a state with lower affinity for the operator DNA sequence. When lactose becomes available, it is converted into allolactose, binds to the lac repressor, causing an allosteric change in its shape (LacI_d). In its changed state, the lac repressor is unable to bind tightly to its cognate operator. This effect is referred to as induction, because it induces, rather than represses expression of the metabolic genes. In vitro, IPTG is commonly used as an allolactose mimic to induce transcription of genes being regulated by lac repressor. IPTG binds to Lac-repressor (LacI) and deactivates it by bringing it to the inactive form LacI_d, this way reducing the effective free Lac-repressor concentration that would negatively affect KinA levels. Therefore, in the presence of IPTG, expressed KinA levels can be increased.
In our model, KinA is the phosphoryl donor that activates the B. subtilis phosphorelay. We include the main phosphorelay species Spo0F, Spo0B and Spo0A , abstracting the effects of the transition-state regulators Sin I/R (sporulation inhibitor (SinR) and its antagonist (SinI) -), based on the results reported by Chastanet et al.. Chastanet and colleagues reported in  that these regulators were quickly repressed early on in sporulation and thus have little effect on the dynamics of Spo0A. Similarly, we do not explicitly consider the abrB repression induced by Spo0A~P, but we model the impact of AbrB depletion on the phosphorelay species KinA, Spo0F and Spo0A, by a Spo0A~P induced positive feedback .
Dephosphorylation of Spo0A has been shown to be carefully monitored by cells at all times during the cell cycle through two families of aspartyl-phosphate phosphatases - the Rap and Spo0E families ,, which provide opportunity for negative signals to influence the cell's decision whether to sporulate or continue with vegetative growth. RapA, B and E indirectly inhibit the phosphorylation of Spo0A via dephosphorylation of Spo0F ,,,, whilst Spo0E dephosphorylates Spo0A directly ,,. To simplify our modelling, Rap and Spo0E phosphatases are only taken into consideration in an abstract manner. The dephosphorylation of Spo0F and Spo0A proteins are modelled as first order reactions assuming that Rap and Spo0E concentrations are not rate limiting.
Spo0A~P up-regulates transcription from spoIIA, spoIIE, and spoIIG promoters -. spollA_t mRNAs are translated into AA, AB and AC proteins, spollE_t mRNAs are translated into IIE protein molecules and spollG_t mRNAs are translated into GA and GB protein molecules. Following the modelling approach proposed by Narula et. al. in , we model the dynamics of transcription up-regulation as being the result of a cooperative binding of active Spo0A to the promoters. Moreover, for each species, apart from SS and IPTG, we model the synthesis and degradation processes. Transcription is always modelled as a single chemical reaction synthesizing mRNA where the synthesis rate is controlled by promoter strengths. Translation is modelled as a single chemical reaction with a rule proportional to the amount of corresponding mRNA.
Design and implementation
Parameter values used for the modelling
k=0.1 nMs -1
mass action: k.[lacI_t]
k = 0.2 s -1
LacI inactivation by IPTG
mass action: k .[lacI_t].[IPTG]
k=0.000175 nM -1 ?s -
mass action: k.[LacI_d]
k=0.016 s -1
k1=0.0275 nMs -1
k2=0.24 nMs -1
k3=1.95 nMs -1
mass action: k.[kinA_t]
k=0.0659 s -1
mass action: k.[KinA] 2
k=0.001 nM -1 ?s -1
KinA dimer dissociation
mass action: k.[KinA:KinA]
k=0.25 s -1
KinA dimer auto-phosphorylation
mass action: k.[KinA:KinA].[SS]
k=0.001 nM -1 ?s -1
mass action: k.[Spo0F]
mass action: k.[Spo0F] [KinA:KinA~P]
k=0.00039 nM -1 ?s -1
mass action: k.[Spo0F~P]
k=0.05 s -1
k=0.2384 nMs -1
mass action: k.[Spo0B]
k=0.1076 s -1
mass action: k.[Spo0B] [Spo0F~P]
k=0.00001 nM -1 ?s -1
k1=0.01388 nMs -1
k=0.1388 nMs -1 ,
mass action: k.[Spo0A]
k=0.2143 s -1
mass action: k.[Spo0A] [Spo0B~P]
k=0.0008 nM -1 ?s -1
mass action: k.[Spo0~P]
k=0.05 s -1
k1=0.0277 nMs -1
k2=0.4166 nMs -1
mass action: k [spollA_t]
k=0.1250 s -1
mass action: k [spollA_t]
k=0.0555 s -1
mass action: k [spollA_t]
k=0.0138 s -1
k1=0.0208 nMs -1
k2=0.3125 nMs -1
mass action: k [spollE_t]
k=0.0138 s -1
k1=0.0222 nMs -1
k2=0.7290 nMs -1
mass action: k [spollG_t]
k=0.0034 s -1
mass action: k [spollG_t]
k=0.0138 s -1
As shown in Figure 2, the LacI gene is constitutively expressed. The LacI repression of kinA gene expression is encoded as a negative cooperative effect, whereby the kinA_t mRNA is produced at a rate inversely proportional to LacI n , with n=2. The deactivation effect of IPTG on LacI is modelled through the reversible sequestering of LacI into an inactive form LacI_d that is not capable of repressing the kinA gene expression . The transcription of the kinA gene is also regulated by the feedback loop exerted by Spo0A~P. According to , we assume that KinA forms a homo-complex  via its N- and C-terminal domains and that KinA auto-phosphorylation is activated in response to an unspecified sporulation signal modelled by the SS species.
The phosphorelay circuit of our model includes the main phosphorelay species Spo0F, Spo0B and Spo0A, the phosphorylation, dephosphorylation and phosphotransfer reactions, and also synthesis and degradation of all species. Transcription of spo0F and spo0A genes is regulated by the amount of Spo0A~P. The dephosphorylation of Spo0F and Spo0A proteins are modelled as simple first order reactions, and due to the lack of knowledge about the regulation of phosphotransferase Spo0B, we do not include in the model any transcriptional control nor dephosphorylation reactions for this species. However, it is important to note that, according to experimental evidence , the concentration of Spo0B increases considerably upon sporulation initiation in wild type B. subtilis, suggesting the existence of additional regulation mechanisms. Finally, Spo0A ~ P up-regulates transcription from spoIIA, spoIIE, and spoIIG promoters -. Following the assumptions in Narula et.al, we have reflected these regulations in our model through cooperative activation effects on the transcription rates. spollA_t mRNAs are translated into AA, AB and AC proteins, spollE_t mRNAs are translated into IIE protein molecules and spollG_t mRNAs are translated into GA and GB protein molecules. As for the species and reactions considered here, this model is an extension of the model presented by  for IPTG inducible systems, in that details of the IPTG activity taken from , an additional stimulus (SS) representing the natural induction of sporulation in wild-type, and a positive feedback loop from Spo0A~P to KinA were all added to the model. Concerning the details of the kinetics of reactions, we obtained a number of the parameters from the literature. It is noteworthy that there is no consensus on some important parameters e.g., the details of the kinetics of the regulated transcription of the phosphorelay species and of the downstream spoll targets, which are not exactly known and are inconsistently reported. We report in Table 1 the source of the specific kinetic information for each of the reactions we defined in our model.
Model analysis methods
The overall model includes 13 species and taking into consideration transcripts, dimerization and post-translational modifications, therefore gives rise to 27 distinct forms. The total number of unidirectional reactions included in the model is 55. To analyse such a large model, we encoded the biochemical reactions in the COPASI software package , and used the algorithms implemented within the software for analysis.
The model was analysed using COPASI by adopting a continuous interpretation of the variables, with values that represent concentrations of molecular species, and a deterministic interpretation of the reaction rates, so that the reaction rates provided in the previous tables represent variation terms in the first time derivative of the concentrations. The model is encoded as a set of coupled ordinary differential equations. Therefore, the time-dependent transient behaviour of the model has been computed through numerical integration of the differential equations and the steady-state analysis by the solution of the set of polynomial equations that is obtained by setting to zero the value of the time derivatives in the differential equations.
In modelling and simulation, parameters are distinguished as being entities of a model that are either constant, under our direct control or vary independently - examples of which are kinetic constants and time. Variables are entities whose values are entirely determined by parameters - examples of which are the internal reactant concentrations and fluxes (and of course any other quantities calculated from these) . Sensitivity analysis is an important tool used to study the dependence of a system on external parameters ,, and sensitivity considerations often play an important role in the design of control systems . Parameter sensitivity analysis can also be utilized to validate a model s response, and iteratively to design experiments that support the estimation of parameters .
We conducted parametric sensitivity analyses to determine how our model s behaviour depends on its parameter values. Parametric sensitivity analysis can be divided into global sensitivity analysis (which addresses wide variations in parameter values), and local sensitivity analysis (which addresses small variations around a nominal operating condition). For our study, we exploited the local sensitivity analysis features of COPASI, which allows predicting the effect of small changes to the parameter on the model s behaviour through the computation of relative sensitivity coefficients. The absolute values of coefficients provide insights into the model s robustness, with small sensitivity coefficient indicating robustness with respect to perturbations of a parameter, and large coefficients suggesting high sensitivity. We also conducted global sensitivity analysis for selected parameters, following the indications provided by the local sensitivity results.
In this section we report and comment on the model outcomes, which we obtained by transient, steady-state and sensitivity analyses of the network model. We considered a reference scenario, which is defined by the instantiation of model parameters (initial state and rates) as specified in Table 1. When presenting and discussing our model results, we explored in some cases modelling options that lead to the definition of model parameters different from those in the reference scenarios.
The values of several kinetic parameters (i.e. those that are marked as ‘Fitted’ in Table 1) were determined by systematically exploring the parameter space to find a satisfactory match with the experimental data. The dataset used for the fitting included measurements of one final output effector species (e.g., the spollG transcript), and one intermediate species in the network (KinA). Experimental data for the dynamics of these species under IPTG stimulation were obtained from . The returned predictions of all the other experimental data used for evaluating the quality of the model were not used for the fitting. It is fair to note, that a clear consensus does not exist in the literature for the various parameters for which we determined the nominal value through fitting. For example, for the hill function parameters for the Spo0A regulation of spoll operons transcription, different sets of values are provided, see for instance the two papers published from the same group  and . These discrepancies in the published values of parameters are clearly indicative of areas for further experimental quantification of the kinetics of various biochemical processes.
The model is valid for inducible systems
Since this model is based on previous studies on inducible systems, we first analysed its ability to reproduce the available experimental observations for that class of systems.
Furthermore, our model is able to predict the concentration of the active Spo0A protein in the cell. For the inducible system being modelled, there are no direct experimental measurements of the active form of Spo0A available. Therefore, we used as an indirect reporter the levels of the Spo0A~P directly regulated spollG transcripts. We illustrate in Figure 4 the predictions returned by our model with the experimental data published by Narula and colleagues on the IPTG-dependent levels of expression of spollG. As presented in , the experimental concentration of spollG transcript was measured after 3 hours upon induction with an increasing level of IPTG (on the horizontal axis). The data in Figure 4 are normalized, dividing each value by the concentration when IPTG is zero. The goodness of the agreement is 8.84%.
This steady-state is representative of a situation where the cell is sensing unfavourable conditions of the environment but has not yet committed to sporulation. When the IPTG stimulus is included, the levels of KinA are artificially increased and so the amount of phosphorylated KinA is increased, leading to a stronger activation of the phosphorelay, as shown by the grey and black bars in Figure 5A. Moreover, the positive feedback exerted by Spo0A~P on KinA, Spo0F and Spo0A genes transcription further shifts the cell response to sporulation commitment. As shown in Figure 5B, AA, AB, AC and IIE proteins achieve their maximum levels when IPTG is equal to 10 μM, whereas GA and GB concentration levels show small additional growth when IPTG is increased from 10 μM to 20 μM. This difference in the dynamics is due to the different affinities that Spo0A~P has for the spoll promoters , which lead to the expression of GA and GB only at very high concentrations of Spo0A~P. These results are in line for the observed sporulation rate of IPTG stimulated cells.
The Spo0A positive feedback on KinA is not essential for model output
We included in our model the positive feedback of Spo0A~P on KinA expression. From a point of view of the sporulation network, this is an important structural element, because it establishes a loop between the final effector output of the phosphorelay Spo0A and the originator species, i.e. KinA. However, the authors in  are able to obtain the same output behaviour in terms of the predicted Spo0A activation, without considering this feedback. So, we quantified the relevance of this feedback in the overall architecture of the model.
To check what the consequence of the exclusion of this link would be on the dynamics of the system, we ran a simulation where the contribution of Spo0A~ to KinA expression is set to null (parameter k3 of rate tr2 in Table 1). What we obtained is a response to IPTG of reduced intensity, but qualitatively very similar. Quite interestingly, the net effect of this reduction can be compensated by a very local change of the model. It is indeed sufficient to modify the values of the parameters that regulate KinA transcription through LacI repression (k2 and KL of rate tr2 in Table 1), without any other change in the model, to obtain levels of KinA expression that match the experimental ones.
Since the levels of phosphorylated KinA are practically the same in the two modelling scenarios, it is expected that the level of the phosphorylated forms downstream in the phosphorelay will be the same. Therefore, the two modelling scenarios provide the same response.
The phosphorelay motif transfers the dynamics of kinase activation onto the final effector Spo0A
In both cases, all the species follow the same increasing trend, apart from Spo0B?~?P which is not monotonic. Notice that the vertical scale is logarithmic, and that in both cases (A and B) the stimulus is increased in a linear way. In Figure 7A, we can postulate that the sigmoidal response of the system is determined by the release of the cooperative inhibition on KinA transcription, whereas in Figure 7B the logarithmic shape of the three increasing curves is indicating that the response is linear. In both cases, the shapes of the curves for the first phosphorelay active species KinA?~?P is relayed in the phosphorelay down to effector species Spo0A~P.
The model cannot reproduce both wild type and inducible systems
As shown in Figure 7, our model allows for activated sporulation through the increase of the SS signal. With a proper modulation of the SS signal, this form of induction should allow reproduction of the activation of the master regulatory protein Spo0A in the wild type.
The experimental data set does not show a sigmoidal shape. Although the data is for a time course and not for equilibrium, we would expect the positive feedback acting on KinA, Spo0F and Spo0A, to result in sharp increase of concentrations at the sporulation onset as opposed to linear trends. On the contrary, we find in this experimental data for wild-type a confirmation of the need for a gradual accumulation of the master regulatory protein Spo0A for a proper activation of the downstream targets, which has been suggested by several studies, for instance  and . The amplitude of this linear increase in the total amount of species cannot be reproduced with our model. In particular, the total concentration of Spo0A cannot exceed 4.5 μM, while the experimental data from  that we report in Figure 8 show accumulation of total Spo0A up to 37.5 μM. This important discrepancy is indicative of a profoundly different mechanism of sporulation initiation and highlights important conflicts in the modelling studies our work is based on.
The activation of the final effector Spo0A is not sensitive to Spo0B regulation
We note furthermore, that the experimental data reported by Eswaramoorthy and colleagues in  show a variation of Spo0B concentration as the sporulation commitment of the cell is enforced. More precisely, the experimental evidence shows the total amount of Spo0B has a linear increase as the sporulation process progresses.
This experimental evidence is in contrast to the consensus in literature that says that Spo0B is not regulated. The only work that reports a regulation of Spo0B is a recent paper by Carabetta et al., which states that the phosphotransfer activity of Spo0B may be modulated by a ternary complex formed by YlbF, YmcA and YaaT proteins. The role of this complex appears to be the one of accelerating the phosphorylation of Spo0A, probably by interacting with Spo0F and Spo0B. This hypothesized mechanism of action cannot however explain the increase of the total amount of Spo0B demonstrated in .
This insensitivity poses an interesting question, which is, based on the experimental data reported by Eswaramoorthy et al. , why should the cell expend energy to increase the abundance of Spo0B during sporulation? A conjecture could be that Spo0B abundance increase, which is in accord to Figure 10, does not significantly affect the average value of Spo0A~P but may affect its variability, which we intend to examine by additional running stochastic simulation of sporulation models in future work.
In this paper we distilled the current molecular modelling knowledge about the sporulation initiation process in B. subtilis. In line with previous studies we built a model that includes an artificial sporulation initiation mechanism based on IPTG induction, to allow for tuneable regulation of the sporulation signal, and we also included a mechanism to model induction occurring in B. subtilis wild type. We included in the model, the positive feedback exerted by Spo0A~P on KinA, Spo0F and Spo0A transcription.
Our model predictions are validated against experimental data reported in several studies ,,. We found that it is difficult to reconcile the predicted dynamics determined by a model with the IPTG stimulus with those observed in the wild type and reported by Eswaramoorthy et al. in . Indeed, while the model fits the experimental data for inducible systems, it cannot predict the amounts of phosphorelay proteins that are observed experimentally in wild-type B. subtilis.
This important discrepancy indicates that the reaction kinetics used in our model, which have been extracted from the literature, may not provide a faithful representation of the molecular mechanism for the wild type behaviour. To the best of our knowledge, there is a complete lack of experimental data describing the concentration of phosphorylated forms of the phosphorelay species. Under these circumstances it is impossible to determine what would be a more adequate modelling choice to elucidate the details of the phosphorelay dynamics. From our modelling results we can infer that some of the phosphorylated species, in particular Spo0F~P and Spo0B~P, would not be easily detected by experimental means as the phosphorelay dynamics reveals these to be parsimonious in their accumulation (see Figure 7). However for Spo0A~P our model predicts much higher concentration levels of the order of ?M. The experimental data provided by Fujita et al.  indirectly suggests that Spo0A~P concentration needs to increase substantially in the cell to activate the downstream genes, given the low affinity that the species has for some of the gene promoters. It would be useful to conduct experiments that determine the dynamics of phosphorylated Spo0A in the wild type.
Finally, we would like to draw attention to the experimental data in , which shows that the total concentration of Spo0B, the central phosphorelay species, increases linearly over time during sporulation initiation. This highlights a gap in the knowledge of the biology of the phosphorelay as there is no evidence in the literature relating to Spo0B regulation. Moreover, we showed with sensitivity analyses (see Figure 10), that the steady-state response of the phosphorelay, including that of the main effector Spo0A~P, is quite insensitive to an increase of Spo0B concentration. This has led us to formulate an intriguing question concerning the increase in cellular levels of Spo0B which should have little effect on the downstream gene activation. To find possible answers we plan to explore stochastic simulations to ascertain whether variations in Spo0B concentration have an effect on the variability of the main effector Spo0A~P. The stochastic analysis of models of the phosphorelay of B. subtilis is a task that has been partially tackled by Csik sz-Nagy et al.  and further investigations, may reveal behaviour that is additional to that established through continuous deterministic modelling.
The modelling of the molecular network regulating sporulation in B. subtilis is a task that has so far been tackled by various studies, all based on IPTG inducible systems. In this paper we summarized the current published modelling understanding of sporulation initiation process into a model built with the COPASI modelling and simulation software tool. We found that while the model fitted the experimental data for inducible systems, it did not predict the amounts of phosphorelay proteins observed experimentally in wild-type B. subtilis. This important discrepancy emphasized the necessity for additional experimental data. Specifically, we found no measurements for the abundance of the phosphorylated forms of phosphorelay proteins in wild-type B. subtilis, no knowledge available concerning the regulation of the phosphorelay protein Spo0B during sporulation, and variable measurements of reaction kinetics. In summary, we believe that the aforementioned issues may be some of the aspects currently impairing our ability in building mechanistic models at the molecular level of sporulation regulation in B. subtilis wild type.
The authors gratefully acknowledge the support of the Biotechnology and Biological Sciences Research Council (BBSRC). This research was funded by the BBSRC Institute Strategic Programme BB/J004529/1: The Gut Health and Food Safety ISP. We extend our thanks to Professor Mike W. Peck and his group for help and critical discussions during the course of this work.
- Opinion of the scientific panel on biological hazards (BIOHAZ) related to clostridium spp in foodstuffs. EFSA J. 2005, 199: 1-65.Google Scholar
- Opinion of the scientific panel on biological hazards (BIOHAZ) on bacillus cereus and other bacillus spp in foodstuffs. EFSA J. 2005, 175: 1-48.Google Scholar
- Carlin F: Origin of bacterial spores contaminating foods. Food Microbiol. 2011, 28 (2): 177-182. 10.1016/j.fm.2010.07.008.View ArticlePubMedGoogle Scholar
- Heyndrickx M: The importance of endospore-forming bacteria originating from soil for contamination of industrial food processing. Appl Environ Soil Sci. 2011, 2011: 11-10.1155/2011/561975.View ArticleGoogle Scholar
- Stragier P, Losick R: Molecular genetics of sporulation in Bacillus subtilis. Annu Rev Genet. 1996, 30: 297-241. 10.1146/annurev.genet.30.1.297.View ArticlePubMedGoogle Scholar
- De Jong H, Geiselmann J, Batt G, Hernandez C, Page M: Qualitative simulation of the initiation of sporulation in Bacillus subtilis. Bull Math Biol. 2004, 66 (2): 261-299. 10.1016/j.bulm.2003.08.009.View ArticlePubMedGoogle Scholar
- Fawcett P, Eichenberger P, Losick R, Youngman P: The transcriptional profile of early to middle sporulation in Bacillus subtilis. Proc Natl Acad Sci U S A. 2000, 97 (14): 8063-8068. 10.1073/pnas.140209597.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoffmann A, Levchenko A, Scott ML, Baltimore D: The IkappaB-NF-kappaB signaling module: temporal control and selective gene activation. Science. 2002, 298 (5596): 1241-1245. 10.1126/science.1071914.View ArticlePubMedGoogle Scholar
- Ihekwaba AE, Broomhead DS, Grimley RL, Benson N, Kell DB: Sensitivity analysis of parameters controlling oscillatory signalling in the NF-kappaB pathway: the roles of IKK and IkappaBalpha. Syst Biol. 2004, 1 (1): 93-103. 10.1049/sb:20045009.View ArticleGoogle Scholar
- Nelson DE, Ihekwaba AE, Elliott M, Johnson JR, Gibney CA, Foreman BE, Nelson G, See V, Horton CA, Spiller DG, Edwards SW, McDowell HP, Unitt JF, Sullivan E, Grimley R, Benson N, Broomhead D, Kell DB, White MR: Oscillations in NF-kappaB signaling control the dynamics of gene expression. Science. 2004, 306 (5696): 704-708. 10.1126/science.1099962.View ArticlePubMedGoogle Scholar
- Bischofs IB, Hug JA, Liu AW, Wolf DM, Arkin AP: Complexity in bacterial cell-cell communication: quorum signal integration and subpopulation signaling in the Bacillus subtilis phosphorelay. Proc Natl Acad Sci. 2009, 106 (16): 6459-6464. 10.1073/pnas.0810878106.PubMed CentralView ArticlePubMedGoogle Scholar
- Chastanet A, Vitkup D, Yuan GC, Norman TM, Liu JS, Losick RM: Broadly heterogeneous activation of the master regulator for sporulation in Bacillus subtilis. Proc Natl Acad Sci U S A. 2010, 107 (18): 8486-8491. 10.1073/pnas.1002499107.PubMed CentralView ArticlePubMedGoogle Scholar
- Jabbari S, Heap JT, King JR: Mathematical modelling of the sporulation-initiation network in Bacillus subtilis revealing the dual role of the putative quorum-sensing signal molecule PhrA. Bull Math Biol. 2011, 73 (1): 181-211. 10.1007/s11538-010-9530-7.View ArticlePubMedGoogle Scholar
- Kothamachu VB, Feliu E, Wiuf C, Cardelli L, Soyer OS: Phosphorelays provide tunable signal processing capabilities for the cell. PLoS Comput Biol. 2013, 9 (11): e1003322-10.1371/journal.pcbi.1003322.PubMed CentralView ArticlePubMedGoogle Scholar
- Kuchina A, Espinar L, Garcia-Ojalvo J, Suel GM: Reversible and noisy progression towards a commitment point enables adaptable and reliable cellular decision-making. PLoS Comput Biol. 2011, 7 (11): e1002273-10.1371/journal.pcbi.1002273.PubMed CentralView ArticlePubMedGoogle Scholar
- Narula J, Devi SN, Fujita M, Igoshin OA: Ultrasensitivity of the Bacillus subtilis sporulation decision. Proc Natl Acad Sci U S A. 2012, 109 (50): E3513-E3522. 10.1073/pnas.1213974109.PubMed CentralView ArticlePubMedGoogle Scholar
- Sen S, Garcia-Ojalvo J, Elowitz MB: Dynamical consequences of bandpass feedback loops in a bacterial phosphorelay. PLoS One. 2011, 6 (9): e25102-10.1371/journal.pone.0025102.PubMed CentralView ArticlePubMedGoogle Scholar
- Vishnoi M, Narula J, Devi SN, Dao HA, Igoshin OA, Fujita M: Triggering sporulation in Bacillus subtilis with artificial two-component systems reveals the importance of proper Spo0A activation dynamics. Mol Microbiol. 2013, 90 (1): 181-194.PubMedGoogle Scholar
- Levine JH, Fontes ME, Dworkin J, Elowitz MB: Pulsed feedback defers cellular differentiation. PLoS Biol. 2012, 10 (1): e1001252-10.1371/journal.pbio.1001252.PubMed CentralView ArticlePubMedGoogle Scholar
- Fujita M, Losick R: Evidence that entry into sporulation in Bacillus subtilis is governed by a gradual increase in the level and activity of the master regulator Spo0A. Genes Dev. 2005, 19 (18): 2236-2244. 10.1101/gad.1335705.PubMed CentralView ArticlePubMedGoogle Scholar
- Muchova K, Lewis RJ, Perecko D, Brannigan JA, Ladds JC, Leech A, Wilkinson AJ, Barak I: Dimer-induced signal propagation in Spo0A. Mol Microbiol. 2004, 53 (3): 829-842. 10.1111/j.1365-2958.2004.04171.x.View ArticlePubMedGoogle Scholar
- Varughese KI: Molecular recognition of bacterial phosphorelay proteins. Curr Opin Microbiol. 2002, 5 (2): 142-148. 10.1016/S1369-5274(02)00305-3.View ArticlePubMedGoogle Scholar
- Seredick S, Spiegelman GB: Lessons and questions from the structure of the Spo0A activation domain. Trends Microbiol. 2001, 9 (4): 148-151. 10.1016/S0966-842X(01)01981-3.View ArticlePubMedGoogle Scholar
- Fujita M, Gonzalez-Pastor JE, Losick R: High- and low-threshold genes in the Spo0A regulon of Bacillus subtilis. J Bacteriol. 2005, 187 (4): 1357-1368. 10.1128/JB.187.4.1357-1368.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Perego M, Wu JJ, Spiegelman GB, Hoch JA: Mutational dissociation of the positive and negative regulatory properties of the Spo0A sporulation transcription factor of Bacillus subtilis. Gene. 1991, 100: 207-212. 10.1016/0378-1119(91)90368-L.View ArticlePubMedGoogle Scholar
- Trach K, Burbulys D, Strauch M, Wu JJ, Dhillon N, Jonas R, Hanstein C, Kallio P, Perego M, Bird T, Spiegelman G, Fogher C, Hoch JA: Control of the initiation of sporulation in Bacillus subtilis by a phosphorelay. Res Microbiol. 1991, 142 (78): 815-823. 10.1016/0923-2508(91)90060-N.View ArticlePubMedGoogle Scholar
- York K, Kenney TJ, Satola S, Moran CP, Poth H, Youngman P: Spo0A controls the sigma A-dependent activation of Bacillus subtilis sporulation-specific transcription unit spoIIE. J Bacteriol. 1992, 174 (8): 2648-2658.PubMed CentralPubMedGoogle Scholar
- Satola SW, Baldus JM, Moran CP: Binding of Spo0A stimulates spoIIG promoter activity in Bacillus subtilis. J Bacteriol. 1992, 174 (5): 1448-1453.PubMed CentralPubMedGoogle Scholar
- Satola S, Kirchman PA, Moran CP: Spo0A binds to a promoter used by sigma A RNA polymerase during sporulation in Bacillus subtilis. Proc Natl Acad Sci U S A. 1991, 88 (10): 4533-4537. 10.1073/pnas.88.10.4533.PubMed CentralView ArticlePubMedGoogle Scholar
- Lewis PJ, Wu LJ, Errington J: Establishment of prespore-specific gene expression in Bacillus subtilis: localization of SpoIIE phosphatase and initiation of compartment-specific proteolysis. J Bacteriol. 1998, 180 (13): 3276-3284.PubMed CentralPubMedGoogle Scholar
- Strauch M, Webb V, Spiegelman G, Hoch JA: The SpoOA protein of Bacillus subtilis is a repressor of the abrB gene. Proc Natl Acad Sci U S A. 1990, 87 (5): 1801-1805. 10.1073/pnas.87.5.1801.PubMed CentralView ArticlePubMedGoogle Scholar
- Strauch MA, Wu JJ, Jonas RH, Hoch JA: A positive feedback loop controls transcription of the spoOF gene, a component of the sporulation phosphorelay in Bacillus subtilis. Mol Microbiol. 1993, 7 (6): 967-974. 10.1111/j.1365-2958.1993.tb01188.x.View ArticlePubMedGoogle Scholar
- Tojo S, Hirooka K, Fujita Y: Expression of kinA and kinB of Bacillus subtilis, necessary for sporulation initiation, is under positive stringent transcription control. J Bacteriol. 2013, 195 (8): 1656-1665. 10.1128/JB.02131-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Jiang M, Shao W, Perego M, Hoch JA: Multiple histidine kinases regulate entry into stationary phase and sporulation in Bacillus subtilis. Mol Microbiol. 2000, 38 (3): 535-542. 10.1046/j.1365-2958.2000.02148.x.View ArticlePubMedGoogle Scholar
- Rowland SL, Burkholder WF, Cunningham KA, Maciejewski MW, Grossman AD, King GF: Structure and mechanism of action of Sda, an inhibitor of the histidine kinases that regulate initiation of sporulation in Bacillus subtilis. Mol Cell. 2004, 13 (5): 689-701. 10.1016/S1097-2765(04)00084-X.View ArticlePubMedGoogle Scholar
- Molle V, Nakaura Y, Shivers RP, Yamaguchi H, Losick R, Fujita Y, Sonenshein AL: Additional targets of the Bacillus subtilis global regulator CodY identified by chromatin immunoprecipitation and genome-wide transcript analysis. J Bacteriol. 2003, 185 (6): 1911-1922. 10.1128/JB.185.6.1911-1922.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in Escherichia coli. Nature. 2000, 403 (6767): 339-342. 10.1038/35002131.View ArticlePubMedGoogle Scholar
- Eswaramoorthy P, Duan D, Dinh J, Dravis A, Devi SN, Fujita M: The threshold level of the sensor histidine kinase KinA governs entry into sporulation in bacillus subtilis. J Bacteriol. 2010, 192 (15): 3870-3882. 10.1128/JB.00466-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Eswaramoorthy P, Guo T, Fujita M: In vivo domain-based functional analysis of the major sporulation sensor kinase, KinA, in Bacillus subtilis. J Bacteriol. 2009, 191 (17): 5358-5368. 10.1128/JB.00503-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Winnen B, Anderson E, Cole JL, King GF, Rowland SL: Role of the PAS sensor domains in the Bacillus subtilis sporulation kinase KinA. J Bacteriol. 2013, 195 (10): 2349-2358. 10.1128/JB.00096-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee J, Tomchick DR, Brautigam CA, Machius M, Kort R, Hellingwerf KJ, Gardner KH: Changes at the KinA PAS-A dimerization interface influence histidine kinase function. Biochemistry. 2008, 47 (13): 4051-4064. 10.1021/bi7021156.View ArticlePubMedGoogle Scholar
- Wang L, Fabret C, Kanamaru K, Stephenson K, Dartois V, Perego M, Hoch JA: Dissection of the functional and structural domains of phosphorelay histidine kinase A of Bacillus subtilis. J Bacteriol. 2001, 183 (9): 2795-2802. 10.1128/JB.183.9.2795-2802.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Lewis M: The lac repressor. Comptes Rendus Biol. 2005, 328 (6): 521-548. 10.1016/j.crvi.2005.04.004.View ArticleGoogle Scholar
- Dunaway M, Olson JS, Rosenberg JM, Kallai OB, Dickerson RE, Matthews KS: Kinetic studies of inducer binding to lac repressor.operator complex. J Biol Chem. 1980, 255 (21): 10115-10119.PubMedGoogle Scholar
- Bell CE, Lewis M: A closer view of the conformation of the Lac repressor bound to operator. Nat Struct Biol. 2000, 7 (3): 209-214. 10.1038/73317.View ArticlePubMedGoogle Scholar
- Burbulys D, Trach KA, Hoch JA: Initiation of sporulation in B. subtilis is controlled by a multicomponent phosphorelay. Cell. 1991, 64 (3): 545-552. 10.1016/0092-8674(91)90238-T.View ArticlePubMedGoogle Scholar
- Mandic-Mulec I, Gaur N, Bai U, Smith I: Sin, a stage-specific repressor of cellular differentiation. J Bacteriol. 1992, 174 (11): 3561-3569.PubMed CentralPubMedGoogle Scholar
- Shafikhani SH, Mandic-Mulec I, Strauch MA, Smith I, Leighton T: Postexponential regulation of sin operon expression in Bacillus subtilis. J Bacteriol. 2002, 184 (2): 564-571. 10.1128/JB.184.2.564-571.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Cervin MA, Lewis RJ, Brannigan JA, Spiegelman GB: The Bacillus subtilis regulator SinR inhibits spoIIG promoter transcription in vitro without displacing RNA polymerase. Nucleic Acids Res. 1998, 26 (16): 3806-3812. 10.1093/nar/26.16.3806.PubMed CentralView ArticlePubMedGoogle Scholar
- Perego M, Hanstein C, Welsh KM, Djavakhishvili T, Glaser P, Hoch JA: Multiple protein-aspartate phosphatases provide a mechanism for the integration of diverse signals in the control of development in B. subtilis. Cell. 1994, 79 (6): 1047-1055. 10.1016/0092-8674(94)90035-3.View ArticlePubMedGoogle Scholar
- Perego M, Hoch JA: Negative regulation of Bacillus subtilis sporulation by the spo0E gene product. J Bacteriol. 1991, 173 (8): 2514-2520.PubMed CentralPubMedGoogle Scholar
- McQuade RS, Comella N, Grossman AD: Control of a family of phosphatase regulatory genes (phr) by the alternate sigma factor sigma-H of Bacillus subtilis. J Bacteriol. 2001, 183 (16): 4905-4909. 10.1128/JB.183.16.4905-4909.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Perego M, Glaser P, Hoch JA: Aspartyl-phosphate phosphatases deactivate the response regulator components of the sporulation signal transduction system in Bacillus subtilis. Mol Microbiol. 1996, 19 (6): 1151-1157. 10.1111/j.1365-2958.1996.tb02460.x.View ArticlePubMedGoogle Scholar
- Diaz AR, Stephenson S, Green JM, Levdikov VM, Wilkinson AJ, Perego M: Functional role for a conserved aspartate in the Spo0E signature motif involved in the dephosphorylation of the Bacillus subtilis sporulation regulator Spo0A. J Biol Chem. 2008, 283 (5): 2962-2972. 10.1074/jbc.M709032200.View ArticlePubMedGoogle Scholar
- Stephenson SJ, Perego M: Interaction surface of the Spo0A response regulator with the Spo0E phosphatase. Mol Microbiol. 2002, 44 (6): 1455-1467. 10.1046/j.1365-2958.2002.02974.x.View ArticlePubMedGoogle Scholar
- Eswaramoorthy P, Dinh J, Duan D, Igoshin OA, Fujita M: Single-cell measurement of the levels and distributions of the phosphorelay components in a population of sporulating Bacillus subtilis cells. Microbiology. 2010, 156 (8): 2294-2304. 10.1099/mic.0.038497-0.View ArticlePubMedGoogle Scholar
- Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, Singhal M, Xu L, Mendes P, Kummer U: COPASI a COmplex PAthway SImulator. Bioinformatics. 2006, 22 (24): 3067-3074. 10.1093/bioinformatics/btl485.View ArticlePubMedGoogle Scholar
- Mendes P, Kell D: Non-linear optimization of biochemical pathways: applications to metabolic engineering and parameter estimation. Bioinformatics. 1998, 14 (10): 869-883. 10.1093/bioinformatics/14.10.869.View ArticlePubMedGoogle Scholar
- Ingalls BP, Sauro HM: Sensitivity analysis of stoichiometric networks: an extension of metabolic control analysis to non-steady state trajectories. J Theor Biol. 2003, 222 (1): 23-36. 10.1016/S0022-5193(03)00011-0.View ArticlePubMedGoogle Scholar
- Kuo BC: Automatic Control Systems. 1987, Prentice Hall, Upper Saddle River, New Jersey, USAGoogle Scholar
- Cho KH, Shin SY, Kolch W, Wolkenhauer O: Experimental design in systems biology based on parameter sensitivity analysis with monte carlo simulation: a case study for the TNFalpha mediated NF-kappaB signal transduction pathway. Simulation. 2003, 79: 11-12.Google Scholar
- Gentleman R: R Programming for Bioinformatics. 2008, Chapman & Hall/CRC, LondonView ArticleGoogle Scholar
- Carabetta VJ, Tanner AW, Greco TM, Defrancesco M, Cristea IM, Dubnau D: A complex of YlbF, YmcA and YaaT regulates sporulation, competence and biofilm formation by accelerating the phosphorylation of Spo0A. Mol Microbiol. 2013, 88 (2): 283-300. 10.1111/mmi.12186.PubMed CentralView ArticlePubMedGoogle Scholar
- Grimshaw CE, Huang S, Hanstein CG, Strauch MA, Burbulys D, Wang L, Hoch JA, Whiteley JM: Synergistic kinetic interactions between components of the phosphorelay controlling sporulation in Bacillus subtilis. Biochemistry. 1998, 37 (5): 1365-1375. 10.1021/bi971917m.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.