BACKGROUND
lactobacillus delbrueckii subsp.
bulgaricus has been widely applied in diary industry, especially as a starter for yogurt production.
this microorganism is facultative anaerobic and can produce lactic acid.
during the later period of fermentation and storage of yogurt, the production of lactic acid is dominated by l. bulgaricus  <cit> .
the excessive production of lactic acid during storage makes the yogurt taste too sour and this phenomenon is called post-acidification.
at industrial level, several attempts have been made to slow down post-acidification, such as using devoid of l. bulgaricus, selecting the weak post-acidification l. bulgaricus, pasteurization after fermentation and so on  <cit> .
however, these methods either control the post-acidification inefficiently or kill the probiotics in yogurt.
using genetic technology to modify the bacteria has become another option to decline post-acidification.
due to lacking molecular tools, mainly caused by the absence of a reliable transformation procedure, our understanding of the physiology and genetics of l. bulgaricus is still limited  <cit> .
in l. bulgaricus, the absorption of lactose is processed by lactose/galactose antiport transport system  <cit> .
firstly, lactose is transported into cells by lactose permease encoded by lacs gene.
then the lactose is decomposed into glucose and galactose by β-galactosidase encoded by lacz gene.
the free galactose is pumped out of cells or stored in cells in the form of macromolecule i.e.
carbohydrate gum.
glucose turns into pyruvate through glycolysis, and the glycolysis is transformed into lactic acid by the catalysis of lactate dehydrogenase  <cit> .
the lacs gene and the lacz gene constitute a lacsz operon.
in the downstream of lacz gene, there is a lacr gene encoding a repressor which makes the lacsz operon induced by lactose.
however, in l. bulgaricus, the lacr gene has lost regulatory function due to the insertion of some gene fragments, resulting in constitutive expression of lacsz operon  <cit> .
due to our limited knowledge, most attempts at the genetic level focused on changing the activity of one significant enzyme at a time in the metabolism of lactic acid.
druesne altered the histidine codon to alanine codon at the 552th locus of lactose permease, resulting in the reduction of the enzyme's activity  <cit> .
early in  <dig>  mollet obtained the mutant l. bulgaricus without β-galactosidase activity by spontaneous deletion  <cit> .
the new aforementioned organisms did produce less lactic acid but they had trouble in growing in milk independently due to the lack of energy.
adams selected two cold-sensitive mutants of the β-galactosidase from l. bulgaricus by using the expression system of e. coli  <cit> .
however, their later research turned back to directly mutate and screen l. bulgaricus since it was difficult to transform with l. bulgaricus.
so far, lactobacillus delbrueckii subsp.
bulgaricus atcc  <dig> has been reported successfully transformed by electroporation, although with very low reproducibility and efficiency  <cit> .
several ph-induced promoters from lactococcus lactis have been demonstrated, such as rcfb  <cit> , p <dig> and p <dig>  <cit> .
combining the two points mentioned above, we came up with a new method to build a ph-induced promoter with a repressor gene controlling the production of lactic acid.
thus, we could turn on the switch at appropriate time.
here we analyse the parameters in the kinetic model to investigate the regulation effect of ph-induced promoter on lactic acid.
methods
the reduced metabolic network
here we present a model of the production of lactic acid, which includes the important enzymes, lactose permease, β-galactosidase, lactate dehydrogenase and related regulation genes lacs, lacz, lacr, as shown in figure  <dig>  the process of glycolysis is reflected in one reaction, adopting the assumption that there is a constant flux from glucose to pyruvate.
pyruvate is partly catalysed to lactic acid, and the other two products, acetyl-coa and butanedione, are not taken into account.
the mathematical model
the mathematical model describes how the extracellular ph value influences the production of lactic acid.
the ordinary differential equations in the model are based on michaelis-menten equations  <cit>  and fundamental kinetic principles.
the variables are described in table  <dig>  the model assumptions are as follows.
each entity, described in table  <dig> is given as total amount in the population.
the model does not take energy metabolism into consideration.
all reactions are modelled by mass action principles, except for enzyme reactions which obey michaelis-menten equation, and transcription which obeys saturation kinetics.
the level of translation product is assumed to be proportional to mrna levels according to dockery and keener  <cit> .
there is no delay in synthesis of either component or delay because of protein transportations.
all substances are subject to degradation.
the metabolism of lactose
as mentioned before, the metabolism of lactose is described by michaelis-menten equation.
since lactose transport is reversible, a term was included to count for lactose efflux dependent on the internal lactose concentration  <cit> .
the vmax for lactose permease and β-galactosidase is associated with the concentration of enzymes, so we use kcat and the concentration of enzymes to estimate the vmax.
since the lac operon is usually induced by lactose, we assume that the repressor lacr protein could react with lactose.
thus we need to take the reduction into consideration when evaluating the concentration of intracellular lactose and lacr protein.
as for the utilization of sugar, l. bulgaricus belongs to the type of homofermentation.
theoretically, this microorganism could totally convert glucose to lactic acid, but the conversion of glucose is about 80%-90% in practice  <cit> .
therefore, the estimated proportion of pyruvate transformed into lactic acid is 80%.
the gene regulation
the transcription and decomposition of lactose are tightly regulated at the genetic level: production of the enzymes can be decreased at the transcriptional level by regulatory protein lacr binding at appropriate dna sites.
we set m to stand for the transcript concentration of lacr gene.
since lacs and lacz are located in the same operon, it is assumed that both genes are transcribed at same rate.
we used n to describe the concentration of the two genes' mrna.
when describing the equations, the mrnas are assumed to be transcribed at two distinct rates: basal one when there is no regulation and a higher/lower rate when being regulated  <cit> .
the transcription of lacr gene
for m, it is described as function .
the basal rate is given by v <dig> and increased at a constant rate v <dig> which is regulated by the ph value.
gr stands for the repressor gene concentration.
the degradation of the mrna occurs at the rate of dm.
it is also diluted due to cell growth rate at μ. we assumed the decrease of mrna is first-order proportion to the concentration of m.
dmdt=v0+v1⋅fph⋅gr-dm+μ⋅m
as the mechanisms of ph-induced promoters have not been fully elucidated, we included a generic ph-dependent switch f which turns on when ph value is below the threshold.
the switch takes the form of the smoothed step function  <cit> , where  <dig>  represents the threshold ph level around which the switch occurs  <cit> .
n indicates the steepness of the smooth switch function.
f=1-tanhn⋅ph- <dig>
since lactic acid is the main product during fermentation, we assumed ph value is all controlled by the concentration of lactic acid.
besides, we needed to take the buffer capacity of milk medium into account.
furthermore, we assumed that all the lactic acid could spread to the medium even when the extracellular concentration was pretty high.
we tested the ph value of milk medium when different amount of lactic acid was added.
three parallel experiments were taken and we fitted a function as the form of function .
to the ph data, c <dig>  c <dig> and c <dig> are constants.
this function was then fed into the model via the switch function f.
ph=c0-c1⋅la+c2⋅ <dig>
transcription of lacs and lacz gene
for n, it is described as function .
the basal rate is given by v <dig> and the transcription is inhibited by lacr.
for the part to simulate the regulation rate of lacr, tmax stands for the maximal rate of lacs and lacz gene transcription, s represents the sensitivity of lacs and lacz gene transcription to lactose permease and β-galactosidase  <cit> .
g is the concentration of transcription gene.
again, we assumed first-order degradation rate at dn and dilution rate at μ of the mrna.
dndt=v2-tmax⋅s⋅lr1+s⋅lr⋅g-dn+μ⋅n
cell growth
in milk medium, lactose is the main carbon source, so it is assumed that the cell growth rate is dependent on the extracellular lactose concentration.
the product, lactic acid, which declines the ph value of milk medium performs inhibition on cell growth.
the function to describe the production inhibition is the same as concepcion presented  <cit> .
in function , ks is the monod constant for growth in extracellular lactose, and μmax is the maximum specific growth rate.
kla is the maximum initial lactic acid concentration in which the microorganism growth is completely inhibited.
μ=μmax⋅lacoutlacout+ks⋅1-lakla
combining the metabolism of lactose with gene regulation
all the definitions of the parameters are shown in table  <dig>  combining with the enzyme reactions, gene regulations and cell growth discussed above, the process described in figure  <dig> could be represented by:
r1=kcat⋅lp⋅ls⋅lacoutlacout+km⋅lpr2=kcat⋅lp⋅ls⋅lacinlacin+km⋅lpr3=kcat⋅gal⋅lz⋅lacinlacin+km⋅galr4=vmax⋅glyc⋅gluglu+km⋅glycr5=vmax⋅ldh⋅pyrpyr+km⋅ldhr6=v0⋅grr7=v1⋅f⋅grr8=dm+μ⋅mr9=klr⋅mr10=kb⋅lr⋅lacinr11=dlr+μ⋅lrr12=v2⋅gr13=tmax⋅s⋅lr1+s⋅lr⋅gr14=dn+μ⋅nr15=kls⋅nr16=dls+μ⋅lsr17=klz⋅nr18=dlz+μ⋅lzr19=μ⋅x
the final metabolic model is given by:
dlacoutdt=-r1+r2dlacindt=r1-r2-r3-r10dgludt=r3-r4dpyrdt=2r4- <dig> r5dladt= <dig> r5dmdt=r6+r7-r8dlrdt=r9-r10-r11dndt=r12-r13-r14dlsdt=r15-r16dlzdt=r17-r18dxdt=r <dig>
we set stoichiometric constant of r <dig> with two since two molecules of pyruvate are formed from one molecule of glucose.
for the stoichiometric constant of r <dig>  it stands for the proportion of pyruvate which is converted to lactic acid.
data estimated
firstly, we did not take lacr gene into consideration.
this is just the case of wild type l. bulgaricus which lacsz operon is constitutive expression.
since most of the parameters in this reduced model have been reported, v <dig> is the only one needs to be modified.
we used sbtoolbox in matlab to set up the reduced model and simulated the metabolism process  <cit> .
we compared the concentration changes of extracellular lactose and lactic acid within  <dig> hours with the data from fatama, where the conditions used were  <dig> °c, 8% of skim milk concentration  and 4% of inoculum ratio  <cit> .
after that, we added lacr gene.
this time we did not have sufficient data for full parameterization of the model.
however, we did have enough information to estimate their relative sizes and give the qualitative nature of our investigations.
the parameter values are listed in table  <dig>
sensitivity and significance analysis of parameters
we still used sbtoolbox to perform local sensitivity analysis and sobol's method for global sensitivity analysis.
the local sensitivity analysis investigates how small changes in a single parameter value could affect the model output.
the method is based on the partial differentiation of the output with respect to the input parameters  <cit> .
herein, the partial differentiation is evaluated numerically by introducing a 1% increment from the specific parameter value.
we chose sobol sensitivity analysis method to calculate global sensitivity since the sensitivity index quantified the overall effects of a parameter, in combination with any other parameters, on the model output  <cit> .
the number of simulation to carry out is  <dig> times.
we first set the number as  <dig> according to schmidt  <cit> .
since the results are varied among different simulations, we increased the number to  <dig> and got stable outcomes.
during the above sensitivity analysis, we chose lactic acid concentration as the model output.
after figuring out the sensitivity parameters to the model output, we carried out different simulations with each sensitivity parameter at two levels, a relative low value and a relative high value.
the metabolism of lactic acid and changing of ph values was carefully considered.
those results indicate how ph values affect lactic acid production under different situations and help us to understand how to design lacr regulation in biology experiments in future.
RESULTS
reduced model without lacr gene
the results of modelled concentration changing of extracellular lactose and lactic acid are shown in figure  <dig>  along with the experimental data.
the model fits the experimental data well.
we have obtained the value of v <dig>  which indicates the basal rate of lacsz gene transcription.
the relative size of fitting value  is determined according to gustafsson's  <cit>  report.
sensitivity analysis of parameters
the results of local sensitivity analysis and sobol sensitivity analysis are shown in figure  <dig>  the parameters with sensitivity from high to low in local sensitivity analysis are followed as tmax, gr, klr, s, v <dig>  dm, v <dig>  umax, dlr and v <dig>  the other parameters perform little effect on the production of lactic acid.
the reason is that when one of those parameters increased, it would only affect the production rate of the output but have no effect on the total amount of lactic acid.
thus, they are not sensitive parameters in this method.
as for sobol sensitivity analysis, the parameters with sensitivity from high to low is tmax, km·lp, klz, klr, gr, km·lyc, v <dig>  s, dls, v <dig>  dn, v <dig> and dlr.
the parameter km·lp shows the second most sensitivity in sobol's method, just as druesne has proven, the reduction of lactose permease's activity would greatly decline the production of lactic acid  <cit> .
significant analysis of parameters
in sensitivity analysis, sensitivity does not mean importance, since sensitive parameters are always fixed inherently and unadjustable.
we need to combine the sensitivity results with operability, especially for the design of biology experiments.
although tmax, klz, v <dig>  dls, dlr, gr, dm, dn and umax perform high sensitivity in both methods, these parameters are usually inherently.
km·lp value is adjustable, but to decrease its value would generate fatal disadvantage on the growth of the bacteria.
so we do not include this parameter into consideration.
due to the same reason, we exclude km·glyc, which is the michaelis constant for glycolysis.
as for the parameter s, which means sensitivity of lacsz gene transcription to lacs and lacz proteins, it is difficult to measure and adjust in experiments.
therefore, the parameters which exhibit great influence on the production of lactic acid in both methods are klr, v <dig> and v <dig>  with klr representing rbs strength of lacr gene, and v <dig> together with v <dig> reflecting characters of the ph-induced promoter.
for the three parameters, the relative high levels of klr, v <dig> and v <dig> are  <dig> h- <dig>   <dig> h- <dig> and  <dig> h- <dig> respectively.
the relative low levels of klr, v <dig> and v <dig> are  <dig> h- <dig>   <dig>  h- <dig>   <dig>  h- <dig> respectively.
the results are shown in figure  <dig>  from picture a, we could figure out that with the increase of promoter strength, the production rate and amount of lactic acid would be greatly decreased.
however, for the relative high value of v <dig>  the finial ph value only declines to  <dig>  which is higher than casein isoelectric point.
thus we need to find an appropriate promoter strength between the relative high and low values when it is turned on by ph value.
from picture b, we could draw the conclusion that if the ph-induced promoter could be efficiently turned on when ph decreases to its threshold, leaking expression of the promoter would have little effect on the production of lactic acid.
in picture c, the increase of klr leads to strong decline of both the production rate and total amount of lactic acid.
by comparing the results of relative high level klr with the relative high level v <dig> in picture a, the two parameters almost perform the same influence on the production of lactic acid.
this is due to the reason that the two parameters both affect the concentration of lacr.
therefore, we could adjust the strength of ph-induced promoter and/or rbs strength of lacr gene to make our switch perform well in controlling lactic acid production.
other two important features for the promoter are the threshold ph level around which it would start and whether the promoter would keep turning on when the ph value decreased to  <dig>  or lower.
akyol  <cit>  and madsen  <cit>  only described the status of ph-induced promoters along with the ph declined to  <dig> .
further experiments are needed to test the characteristics of the promoters.
this model neglect the effect of lower ph on the promoter, assuming that it would still turn on with the decrease of ph value.
CONCLUSIONS
in conclusion, we propose a new method to control the production of lactic acid by l. bulgaricus which is to build a ph-induced switch  to turn on repressor gene.
the proposed method overcomes the disadvantage in bacterial growth by directly changing one of the enzyme related to lactic acid metabolism.
then, we build a reduced model of ph-induced promoter regulation on the production of lactic acid and adjust the model with data from wild type l. bulgaricus.
after that we carry out sensitivity analysis of the parameters and figure out three significant ones in the model.
to make our switch work well, we need to find an appropriate strength of promoter and/or an appropriate strength of rbs in lacr gene.
the values of those two parameters should between the high level and low level we have set in the analysis, so that when the ph value declines to the promoter's starting point, the lacr gene could express moderately to inhibit the production of lactic acid.
competing interests
the authors declare that they have no competing interests.
authors' contributions
kl, lq, xl, yy proposed the project and designed metabolic and gene regulation network.
kl and xz wrote the model and performed sensitivity analysis.
cd, ah and dx guided the model network design and significance analysis of the parameters.
kl wrote the manuscript.
all authors have read and approved the final manuscript.
declarations
this work and its publication is partially supported by the china natural science foundation , the national high-tech r&d program of china  , heilongjiang province science foundation , heilongjiang province key scientific and technological project , hit science foundation , national most special fund .
this article has been published as part of bmc bioinformatics volume  <dig> supplement  <dig>  2014: selected articles from the 9th international symposium on bioinformatics research and applications : bioinformatics.
the full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/15/s <dig>

 lactobacillus delbrueckii subsp. bulgaricus has been widely applied in diary industry, especially as a starter for yogurt production. this microorganism is facultative anaerobic and can produce lactic acid. during the later period of fermentation and storage of yogurt, the production of lactic acid is dominated by l. bulgaricus. the excessive production of lactic acid during storage makes the yogurt taste too sour and this phenomenon is called post-acidification. at industrial level, several attempts have been made to slow down post-acidification, such as using devoid of l. bulgaricus, selecting the weak post-acidification l. bulgaricus, pasteurization after fermentation and so on. however, these methods either control the post-acidification inefficiently or kill the probiotics in yogurt. using genetic technology to modify the bacteria has become another option to decline post-acidification. due to lacking molecular tools, mainly caused by the absence of a reliable transformation procedure, our understanding of the physiology and genetics of l. bulgaricus is still limited. in l. bulgaricus, the absorption of lactose is processed by lactose/galactose antiport transport system. firstly, lactose is transported into cells by lactose permease encoded by lacs gene. then the lactose is decomposed into glucose and galactose by β-galactosidase encoded by lacz gene. the free galactose is pumped out of cells or stored in cells in the form of macromolecule i.e. carbohydrate gum. glucose turns into pyruvate through glycolysis, and the glycolysis is transformed into lactic acid by the catalysis of lactate dehydrogenase. the lacs gene and the lacz gene constitute a lacsz operon. in the downstream of lacz gene, there is a lacr gene encoding a repressor which makes the lacsz operon induced by lactose. however, in l. bulgaricus, the lacr gene has lost regulatory function due to the insertion of some gene fragments, resulting in constitutive expression of lacsz operon. due to our limited knowledge, most attempts at the genetic level focused on changing the activity of one significant enzyme at a time in the metabolism of lactic acid. druesne altered the histidine codon to alanine codon at the 552th locus of lactose permease, resulting in the reduction of the enzyme's activity. early in mollet obtained the mutant l. bulgaricus without β-galactosidase activity by spontaneous deletion. the new aforementioned organisms did produce less lactic acid but they had trouble in growing in milk independently due to the lack of energy. adams selected two cold-sensitive mutants of the β-galactosidase from l. bulgaricus by using the expression system of e. coli. however, their later research turned back to directly mutate and screen l. bulgaricus since it was difficult to transform with l. bulgaricus. so far, lactobacillus delbrueckii subsp. bulgaricus atcc has been reported successfully transformed by electroporation, although with very low reproducibility and efficiency. several ph-induced promoters from lactococcus lactis have been demonstrated, such as rcfb, p and p. combining the two points mentioned above, we came up with a new method to build a ph-induced promoter with a repressor gene controlling the production of lactic acid. thus, we could turn on the switch at appropriate time. here we analyse the parameters in the kinetic model to investigate the regulation effect of ph-induced promoter on lactic acid. methods the reduced metabolic network here we present a model of the production of lactic acid, which includes the important enzymes, lactose permease, β-galactosidase, lactate dehydrogenase and related regulation genes lacs, lacz, lacr, as shown in figure the process of glycolysis is reflected in one reaction, adopting the assumption that there is a constant flux from glucose to pyruvate. pyruvate is partly catalysed to lactic acid, and the other two products, acetyl-coa and butanedione, are not taken into account. the mathematical model the mathematical model describes how the extracellular ph value influences the production of lactic acid. the ordinary differential equations in the model are based on michaelis-menten equations and fundamental kinetic principles. the variables are described in table the model assumptions are as follows. each entity, described in table is given as total amount in the population. the model does not take energy metabolism into consideration. all reactions are modelled by mass action principles, except for enzyme reactions which obey michaelis-menten equation, and transcription which obeys saturation kinetics. the level of translation product is assumed to be proportional to mrna levels according to dockery and keener. there is no delay in synthesis of either component or delay because of protein transportations. all substances are subject to degradation. the metabolism of lactose as mentioned before, the metabolism of lactose is described by michaelis-menten equation. since lactose transport is reversible, a term was included to count for lactose efflux dependent on the internal lactose concentration. the vmax for lactose permease and β-galactosidase is associated with the concentration of enzymes, so we use kcat and the concentration of enzymes to estimate the vmax. since the lac operon is usually induced by lactose, we assume that the repressor lacr protein could react with lactose. thus we need to take the reduction into consideration when evaluating the concentration of intracellular lactose and lacr protein. as for the utilization of sugar, l. bulgaricus belongs to the type of homofermentation. theoretically, this microorganism could totally convert glucose to lactic acid, but the conversion of glucose is about 80%-90% in practice. therefore, the estimated proportion of pyruvate transformed into lactic acid is 80%. the gene regulation the transcription and decomposition of lactose are tightly regulated at the genetic level: production of the enzymes can be decreased at the transcriptional level by regulatory protein lacr binding at appropriate dna sites. we set m to stand for the transcript concentration of lacr gene. since lacs and lacz are located in the same operon, it is assumed that both genes are transcribed at same rate. we used n to describe the concentration of the two genes' mrna. when describing the equations, the mrnas are assumed to be transcribed at two distinct rates: basal one when there is no regulation and a higher/lower rate when being regulated. the transcription of lacr gene for m, it is described as function. the basal rate is given by v and increased at a constant rate v which is regulated by the ph value. gr stands for the repressor gene concentration. the degradation of the mrna occurs at the rate of dm. it is also diluted due to cell growth rate at μ. we assumed the decrease of mrna is first-order proportion to the concentration of m. dmdt=v0+v1⋅fph⋅gr-dm+μ⋅m as the mechanisms of ph-induced promoters have not been fully elucidated, we included a generic ph-dependent switch f which turns on when ph value is below the threshold. the switch takes the form of the smoothed step function, where represents the threshold ph level around which the switch occurs. n indicates the steepness of the smooth switch function. f=1-tanhn⋅ph- since lactic acid is the main product during fermentation, we assumed ph value is all controlled by the concentration of lactic acid. besides, we needed to take the buffer capacity of milk medium into account. furthermore, we assumed that all the lactic acid could spread to the medium even when the extracellular concentration was pretty high. we tested the ph value of milk medium when different amount of lactic acid was added. three parallel experiments were taken and we fitted a function as the form of function. to the ph data, c c and c are constants. this function was then fed into the model via the switch function f. ph=c0-c1⋅la+c2⋅ transcription of lacs and lacz gene for n, it is described as function. the basal rate is given by v and the transcription is inhibited by lacr. for the part to simulate the regulation rate of lacr, tmax stands for the maximal rate of lacs and lacz gene transcription, s represents the sensitivity of lacs and lacz gene transcription to lactose permease and β-galactosidase. g is the concentration of transcription gene. again, we assumed first-order degradation rate at dn and dilution rate at μ of the mrna. dndt=v2-tmax⋅s⋅lr1+s⋅lr⋅g-dn+μ⋅n cell growth in milk medium, lactose is the main carbon source, so it is assumed that the cell growth rate is dependent on the extracellular lactose concentration. the product, lactic acid, which declines the ph value of milk medium performs inhibition on cell growth. the function to describe the production inhibition is the same as concepcion presented. in function, ks is the monod constant for growth in extracellular lactose, and μmax is the maximum specific growth rate. kla is the maximum initial lactic acid concentration in which the microorganism growth is completely inhibited. μ=μmax⋅lacoutlacout+ks⋅1-lakla combining the metabolism of lactose with gene regulation all the definitions of the parameters are shown in table combining with the enzyme reactions, gene regulations and cell growth discussed above, the process described in figure could be represented by: r1=kcat⋅lp⋅ls⋅lacoutlacout+km⋅lpr2=kcat⋅lp⋅ls⋅lacinlacin+km⋅lpr3=kcat⋅gal⋅lz⋅lacinlacin+km⋅galr4=vmax⋅glyc⋅gluglu+km⋅glycr5=vmax⋅ldh⋅pyrpyr+km⋅ldhr6=v0⋅grr7=v1⋅f⋅grr8=dm+μ⋅mr9=klr⋅mr10=kb⋅lr⋅lacinr11=dlr+μ⋅lrr12=v2⋅gr13=tmax⋅s⋅lr1+s⋅lr⋅gr14=dn+μ⋅nr15=kls⋅nr16=dls+μ⋅lsr17=klz⋅nr18=dlz+μ⋅lzr19=μ⋅x the final metabolic model is given by: dlacoutdt=-r1+r2dlacindt=r1-r2-r3-r10dgludt=r3-r4dpyrdt=2r4- r5dladt= r5dmdt=r6+r7-r8dlrdt=r9-r10-r11dndt=r12-r13-r14dlsdt=r15-r16dlzdt=r17-r18dxdt=r we set stoichiometric constant of r with two since two molecules of pyruvate are formed from one molecule of glucose. for the stoichiometric constant of r it stands for the proportion of pyruvate which is converted to lactic acid. data estimated firstly, we did not take lacr gene into consideration. this is just the case of wild type l. bulgaricus which lacsz operon is constitutive expression. since most of the parameters in this reduced model have been reported, v is the only one needs to be modified. we used sbtoolbox in matlab to set up the reduced model and simulated the metabolism process. we compared the concentration changes of extracellular lactose and lactic acid within hours with the data from fatama, where the conditions used were °c, 8% of skim milk concentration and 4% of inoculum ratio. after that, we added lacr gene. this time we did not have sufficient data for full parameterization of the model. however, we did have enough information to estimate their relative sizes and give the qualitative nature of our investigations. the parameter values are listed in table sensitivity and significance analysis of parameters we still used sbtoolbox to perform local sensitivity analysis and sobol's method for global sensitivity analysis. the local sensitivity analysis investigates how small changes in a single parameter value could affect the model output. the method is based on the partial differentiation of the output with respect to the input parameters. herein, the partial differentiation is evaluated numerically by introducing a 1% increment from the specific parameter value. we chose sobol sensitivity analysis method to calculate global sensitivity since the sensitivity index quantified the overall effects of a parameter, in combination with any other parameters, on the model output. the number of simulation to carry out is times. we first set the number as according to schmidt. since the results are varied among different simulations, we increased the number to and got stable outcomes. during the above sensitivity analysis, we chose lactic acid concentration as the model output. after figuring out the sensitivity parameters to the model output, we carried out different simulations with each sensitivity parameter at two levels, a relative low value and a relative high value. the metabolism of lactic acid and changing of ph values was carefully considered. those results indicate how ph values affect lactic acid production under different situations and help us to understand how to design lacr regulation in biology experiments in future.  reduced model without lacr gene the results of modelled concentration changing of extracellular lactose and lactic acid are shown in figure along with the experimental data. the model fits the experimental data well. we have obtained the value of v which indicates the basal rate of lacsz gene transcription. the relative size of fitting value is determined according to gustafsson's report. sensitivity analysis of parameters the results of local sensitivity analysis and sobol sensitivity analysis are shown in figure the parameters with sensitivity from high to low in local sensitivity analysis are followed as tmax, gr, klr, s, v dm, v umax, dlr and v the other parameters perform little effect on the production of lactic acid. the reason is that when one of those parameters increased, it would only affect the production rate of the output but have no effect on the total amount of lactic acid. thus, they are not sensitive parameters in this method. as for sobol sensitivity analysis, the parameters with sensitivity from high to low is tmax, km·lp, klz, klr, gr, km·lyc, v s, dls, v dn, v and dlr. the parameter km·lp shows the second most sensitivity in sobol's method, just as druesne has proven, the reduction of lactose permease's activity would greatly decline the production of lactic acid. significant analysis of parameters in sensitivity analysis, sensitivity does not mean importance, since sensitive parameters are always fixed inherently and unadjustable. we need to combine the sensitivity results with operability, especially for the design of biology experiments. although tmax, klz, v dls, dlr, gr, dm, dn and umax perform high sensitivity in both methods, these parameters are usually inherently. km·lp value is adjustable, but to decrease its value would generate fatal disadvantage on the growth of the bacteria. so we do not include this parameter into consideration. due to the same reason, we exclude km·glyc, which is the michaelis constant for glycolysis. as for the parameter s, which means sensitivity of lacsz gene transcription to lacs and lacz proteins, it is difficult to measure and adjust in experiments. therefore, the parameters which exhibit great influence on the production of lactic acid in both methods are klr, v and v with klr representing rbs strength of lacr gene, and v together with v reflecting characters of the ph-induced promoter. for the three parameters, the relative high levels of klr, v and v are h- h- and h- respectively. the relative low levels of klr, v and v are h- h- h- respectively. the results are shown in figure from picture a, we could figure out that with the increase of promoter strength, the production rate and amount of lactic acid would be greatly decreased. however, for the relative high value of v the finial ph value only declines to which is higher than casein isoelectric point. thus we need to find an appropriate promoter strength between the relative high and low values when it is turned on by ph value. from picture b, we could draw the conclusion that if the ph-induced promoter could be efficiently turned on when ph decreases to its threshold, leaking expression of the promoter would have little effect on the production of lactic acid. in picture c, the increase of klr leads to strong decline of both the production rate and total amount of lactic acid. by comparing the results of relative high level klr with the relative high level v in picture a, the two parameters almost perform the same influence on the production of lactic acid. this is due to the reason that the two parameters both affect the concentration of lacr. therefore, we could adjust the strength of ph-induced promoter and/or rbs strength of lacr gene to make our switch perform well in controlling lactic acid production. other two important features for the promoter are the threshold ph level around which it would start and whether the promoter would keep turning on when the ph value decreased to or lower. akyol and madsen only described the status of ph-induced promoters along with the ph declined to. further experiments are needed to test the characteristics of the promoters. this model neglect the effect of lower ph on the promoter, assuming that it would still turn on with the decrease of ph value.  in conclusion, we propose a new method to control the production of lactic acid by l. bulgaricus which is to build a ph-induced switch to turn on repressor gene. the proposed method overcomes the disadvantage in bacterial growth by directly changing one of the enzyme related to lactic acid metabolism. then, we build a reduced model of ph-induced promoter regulation on the production of lactic acid and adjust the model with data from wild type l. bulgaricus. after that we carry out sensitivity analysis of the parameters and figure out three significant ones in the model. to make our switch work well, we need to find an appropriate strength of promoter and/or an appropriate strength of rbs in lacr gene. the values of those two parameters should between the high level and low level we have set in the analysis, so that when the ph value declines to the promoter's starting point, the lacr gene could express moderately to inhibit the production of lactic acid. competing interests the authors declare that they have no competing interests. authors' contributions kl, lq, xl, yy proposed the project and designed metabolic and gene regulation network. kl and xz wrote the model and performed sensitivity analysis. cd, ah and dx guided the model network design and significance analysis of the parameters. kl wrote the manuscript. all authors have read and approved the final manuscript. declarations this work and its publication is partially supported by the china natural science foundation, the national high-tech r&d program of china, heilongjiang province science foundation, heilongjiang province key scientific and technological project, hit science foundation, national most special fund. this article has been published as part of bmc bioinformatics volume supplement 2014: selected articles from the 9th international symposium on bioinformatics research and applications: bioinformatics. the full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/15/s 
the proposed method of introducing a ph-induced promoter to regulate a repressor gene could restrain the synthesis of lactic acid if an appropriate strength of promoter and/or an appropriate strength of ribosome binding sequence in lacr gene has been designed. the excessive production of lactic acid by l. bulgaricus during yogurt storage is a phenomenon we are always tried to prevent. a model of ph-induced promoter regulation on the production of lactic acid by l. bulgaricus was built. the modelled lactic acid metabolism without ph-induced promoter regulation fitted well with wild type l. bulgaricus. 20- may 9th international symposium on bioinformatics reseaerch and applications charlotte, nc, usa klr and v both exhibited great influence on the producing of lactic acid. genetic methods of changing the activity of one enzyme related to lactic acid metabolism make the bacteria short of energy to growth, although they are efficient ways in controlling lactic acid production. v had little effect on lactic acid production if the ph-induced promoter could be well induced when ph decreased to its threshold. both the local sensitivity analysis and sobol sensitivity analysis indicated parameters tmax, gr, klr, s, v v and dlr were sensitive.
 the excessive production of lactic acid by l. bulgaricus during yogurt storage is a phenomenon we are always tried to prevent. the methods used in industry either control the post-acidification inefficiently or kill the probiotics in yogurt. genetic methods of changing the activity of one enzyme related to lactic acid metabolism make the bacteria short of energy to growth, although they are efficient ways in controlling lactic acid production.  a model of ph-induced promoter regulation on the production of lactic acid by l. bulgaricus was built. the modelled lactic acid metabolism without ph-induced promoter regulation fitted well with wild type l. bulgaricus. both the local sensitivity analysis and sobol sensitivity analysis indicated parameters tmax, gr, klr, s, v v and dlr were sensitive. in order to guide the future biology experiments, three adjustable parameters, klr, v and v were chosen for further simulations. v had little effect on lactic acid production if the ph-induced promoter could be well induced when ph decreased to its threshold. klr and v both exhibited great influence on the producing of lactic acid.  the proposed method of introducing a ph-induced promoter to regulate a repressor gene could restrain the synthesis of lactic acid if an appropriate strength of promoter and/or an appropriate strength of ribosome binding sequence in lacr gene has been designed. 20- may 9th international symposium on bioinformatics reseaerch and applications charlotte, nc, usa 