BACKGROUND
during anaerobic growth e. coli uses different electron acceptors in the electron transport chain, including nitrite .
given the high toxicity of nitrite, expression regulation of genes involved in the nitrite-associated electron transport chain is closely linked to the expression regulation of nitrite metabolism genes.
nitrite reductases nrfa  and nirb  are two main components of the nitrite-reducing system.
these enzymes have different metabolite activities and localizations in the cell.
nrfa reductase is a respiratory enzyme localized in the cell periplasm; it takes part in the proton gradient generation and uses no <dig> as an electron acceptor.
nirb reductase is restricted to the cytoplasm and, most probably, has excessive nitrite detoxification as its main role in the cell.
both enzymes catalyze the reduction of no <dig> to ammonium, but nrfa is most active at low nitrite levels, whereas nirb activity is observed only at high nitrite levels  <cit> .
differential expression of nrf and nir operons, encoding nrfa and nirb reductases, respectively, enables such a combination of nrfa and nirb reductase activities.
the nrf and nir operons expression is controlled by transcription factors, narl and narp, the activity of which depends on narq and narx kinases  <cit> .
unlike nir operon, the expression of which is always positively regulated by nitrite  <cit> , the nrf operon has regulatory region that allows its expression activation at low nitrite levels and expression inhibition at high nitrite levels  <cit> .
as a result, at nitrite concentrations in the media > 2 mm, genetic system that encodes nrfa reductase and allows nitrite metabolism in the periplasm of the cell switches to a more effective in terms of nitrite recycling genetic system, that encodes nirc transporter and nirb nitrite reductase and allows nitrite transport into the cell and its further recycle by the cytoplasmic reductase  <cit> .
due to a high toxicity of nitrite, cell transport system plays a crucial role in the nitrite level regulation in the cell.
three transport proteins involved in the nitrite transmembrane passage are known: nark, naru and nirc  <cit> , but the nitrite transport protein nirc is known to have the highest activity among three  <cit> .
localization of nirС gene, encoding the nirc protein, in the same operon as nirb gene demonstrates a tight connection between intracellular nitrite utilization and nitrite import/export systems at high nitrite levels; however, complex interplay between periplasmic and cytoplasmic systems, which determines the nitrite accumulation and utilization in the cell as well as its export from the cell, represents a significant knowledge gap.
the existing amount of qualitative and quantitative data on nrfa, nirb and nirС genes expression regulation and nitrite accumulation dynamics during e. coli stationary phase growth in chemostat  <cit> , allowed us to develop a mathematical model of nitrite intracellular utilization and investigate how nirb and nrfa reductases and nirc transporter contribute to this process.
using this model, mechanisms involved in regulation of these proteins’ activities were investigated for their roles in controlling the intracellular nitrite levels.
mathematical model of nitrite utilization in e. coli cells
the exact conditions of e. coli cells cultivation in a nitrite-supplied flow chemostat are reproduced in the model.
the model is adapted to the experimental data on nitrite accumulation dynamics in a chemostat and on the inducer levels dependent nrf and nir operons expression  <cit> .
a steady-state growth rate and cell culture density in the chemostat was achieved by fixating the glucose uptake  <cit> .
the model of nitrite utilization in e. coli cells cultured in a chemostat consists of several elementary subsystems.
subsystems are detailed in accordance with the process scheme depicted in fig.
<dig> and are described below.
when describing these processes, generalized hill functions were applied  <cit> .fig.
<dig> scheme of the processes occurring in the chemostat during nitrite utilization in e. coli cells.
– nitrite inflow and outflow at steady-state flow rate;  – reduction of extracellular nitrite to ammonium mediated by periplasmic nrfa reductase; , – import of the extracellular nitrite into the cell and its export from the cell into the chemostat mediated by the nirc transporter;  – intracellular nitrite utilization mediated by cytoplasmic nirb reductase;  – nrfa and nrfb proteins synthesis from the nrf operon mrna and nirb, nirc and nird synthesis from the nir operon mrna;  – transport of the nrf reductase a and b subunits from the cytoplasm to the periplasm and back;  – formation of the nirc transporter active pentameric form;  – formation of the nirb reductase active heteromeric form;  – formation of the nrf active heterotetrameric form;  – degradation of proteins and their complexes
elementary subsystems of the model
subsystem  describes the nitrite flow through the chemostat
nitrite inflow rate into the chemostat is a constant parameter throughout the experiment and nitrite outflow rate from the chemostat is proportional to its current concentration, so the total nitrite flow rate through the chemostat can be described with the equation dudt=vin−kflow⋅u, where u – extracellular nitrite concentration in the chemostat, vin  – nitrite inflow rate into the chemostat, kflow  – the rate constant for the nitrite outflow.
later on, for convenience, vin parameter will be presented as vin = kflow ⋅ s, where s is a nitrite concentration established in the chemostat at a given vin and kflow values in the absence of the cell culture.
later on we will name this parameter «added nitrite».
subsystem  describes the reduction of extracellular nitrite to ammonium mediated by the periplasmic formate-dependent nrfa reductase
the reduction of extracellular nitrite to ammonium by periplasmic formate-dependent nrfa reductase is described in accordance with the stoichiometric scheme: no2− + 6е− + 7h+ → nh3 + 2h2o with the simplest michaelis–menten model  dudt=−ckdis,nrfa2b2uu×nrfa2b2/km,nrf−nrfa2b2u,dnrfa2b2dt=kdis,nrfa2b2uu×nrfa2b2/km,nrf−nrfa2b2u−kcat,nrfnrfa2b2u,dnrfa2b2dt=−kdis,nrfa2b2uu×nrfa2b2/km,nrf−nrfa2b2u+kcat,nrfnrfa2b2u, where nrfa2b <dig> – nitrite-dependent nrfa reductase concentration in the cell periplasm ), c – cell volume fraction in the chemostat volume, kcat,nrf – catalytic constant of nrfa reductase turnover, km,nrf – the michaelis constant.
subsystem  describes the extracellular nitrite import into the cell
it is known that extracellular nitrite transport from the environment into the cell is mediated by the nirc, naru и nark transporter proteins  <cit> , but due to significantly higher nirc nirtrite transport activity, compared with naru and nark  <cit> , only nirc transporter function was considered in the model.
the structure of the e. coli nirc transporter catalytically active form is not known.
we assume that it should be similar to the salmonella typhimurium nirc structure because both proteins belong to one family  <cit> ; that is why in the model nirc is assumed to be a pentameric protein.
the mechanism of nitrite import by the nirc transporter is also not known, therefore we described this process rate with the simplest system, resulted from the simplest michaelis–menten model dudt=−ckdis,nirc5uu×nirc5/km,nircin−nirc5u,dnirc5dt=kdis,nirc5uu×nirc5/km,nircin−nirc5u−kcat,nircinnirc5u,dnirc5udt=−kdis,nirc5uu×nirc5/km,nircin−nirc5u+kcat,nircinnirc5u, where nirc <dig> – nitrite-dependent concentration of the nirc transporter active form in the cell periplasm, kdis,nirc5u – the rate constant for the u × nirc <dig> complex dissociation into subunits, kcat,nircin – the rate constant for the nitrite import by the nirc transporter, km,nircin – the michaelis constant.
subsystem  describes intracellular nitrite export from the cell to the chemostat
in the model intracellular nitrite export from the cell to the chemostat, as well as its import, is mediated only by the nirc transporter.
the naru and nark transporters contribution to this process was not considered  <cit> .
because the mechanism of nitrite import by the nirc transporter is not yet studied, the rate of this process is described with the simplest equation dwdt=−kdis,nirc5ww×nirc5/km,nircout−nirc5w,dnirc5dt=kdis,nirc5ww×nirc5/km,nircout−nirc5w−kcat,nircoutnirc5w,dnirc5wdt=−kdis,nirc5ww×nirc5/km,nircout−nirc5w+kcat,nircoutnirc5w, where w –intracellular nitrite concentration, kdis,nirc5w – the rate constant for the w × nirc <dig> complex dissociation into subunits, kcat,nircout – the rate constant for the nitrite export by the nirc transporter, km,nircout – the michaelis constant.
subsystem  describes the intracellular nitrite utilization mediated by the nadn-dependent nirb reductase
in the model this process is described in accordance with the following stoichiometric scheme: no2−+3nadh+5h+→nh4++3nad++2h2o.
the utilization rate is calculated with the michaelis–menten law  <cit>  dwdt=−kdis,nirb2dww×nirb2d/km,nirb−nirb2dw,dnirb2ddt=kdis,nirb2dww×nirb2d/km,nirb−nirb2dw−kcat,nirbnirb2dw,dnirb2dwdt=−kdis,nirc5ww×nirc5/km,nircout−nirc5w+kcat,nirbnirb2dw, where nirb2d – nitrite-dependent concentration of the nirb reductase active form, kdis,nirb2dw – the rate constant for the w × nirb2d molecular complex dissociation into subunits, kcat,nirb – the recycle catalytic constant, km,nirb – the michaelis constant.
subsystem  describes nrfa, nrfb, nirb, nirc и nird proteins synthesis
in e. coli cells nrfa and nrfb proteins are synthesized from mrna templates of the nrf operon and nirb, nirc and nird proteins – from mrna templates of the nir operon.
the nrf operon mrna concentration in the cell that exists in the chemostat at a current nitrite concentration u is assumed to be proportional to the relative nrfa-β-gal chimeric protein activity observed at the same nitrite concentration ; and the nir operon mrna concentration in the same cell is assumed to be proportional to the relative nirb-β-gal chimeric protein activity , measured by wang and co-authors  <cit> .fig.
<dig> effect of nitrite concentration on nrfabcd and nirbdc operons expression during anaerobic cell growth.
dots represent experimental values of activities of the nrfa-β-gal and nirb-β-gal chimeric proteins measured in  <cit>  and approximating their theoretical curves m
nrf   and m
nir .
the x-axis – steady-state nitrite concentration  in the chemostat ; the y-axis  – the nrfa-β-gal chimeric protein relative activity,  – the nrfb-β-gal chimeric protein relative activity.
the steady-state nitrite concentration in the chemostat used as a scale bar on the x-axis was calculated earlier  <cit> , and based on the data from the previous research  <cit>
activities of nrfa-β-gal and nirb-β-gal chimeric proteins are dependent on the nitrite concentration and were approximated with generalized hill functions, mnrfu=1+δnrf,1×uknrf,1hnrf,11+uknrf,1hnrf,1×1+δnrf,2×uknrf,2hnrf,21+uknrf,2hnrf, <dig>  mniru=1+δnir,1×uknir,1hnir,1+δnir,2×uknir,2hnir,21+uknir,1hnir,1+uknir,2hnir, <dig>
optimal parameters of the mnrf  and mnir  functions were estimated with the gradient descent algorithm and are represented in the appendix: table  <dig>
as a result, synthesis of nrfa and nrfb proteins can be described with the following system of differential equations: dnrfacdt=ksnrf,ac⋅mnrfu,dnrfbcdt=ksnrf,bc⋅mnrfu, where nrfac and nrfbc – concentrations of nrfa and nrfb proteins in the cytoplasm, ksnrf – the rate constant for the synthesis of nrfa and nrfb proteins, and synthesis of nirb, nirc and nird proteins can be described with the following system of equations: dnirbdt=ksnir,b⋅mniru,dnircdt=ksnir,c⋅mniru,dnirddt=ksnir,d⋅mniru, where nirb, nirc, nird – concentrations of nirb, nirc and nird proteins in the cytoplasm, ksnir – the rate constant for the protein synthesis.
it is considered that the nirc protein translation level is more than five times higher than such of the nirb protein in the millimolar range of nitrite concentration  <cit>
subsystem  describes transport of the nrfa and nrfb subunits of the nrf nitrite reductase from the cytoplasm to the periplasm
the nrf enzyme is assumed to be transported from the cytoplasm to the periplasm as monomer subunits nrfa and nrfb, which are synthezised in the cytoplasm.
we also assume that the specific rate constant for the nrf reductase a and b monomers transport from the cytoplasm, where their synthesis takes place, to the periplasmic space, where the active enzyme is being formed, is a variable that depends on the membrane potential: ktnrf ⋅ ku).
generalized hill functions is used for the ku description ktnrf,cp,uus=ktnrf,cp1+duus, whose value is made dependent on the u function, named the nominal potential.
in its turn, the u function mimics the membrane potential that changes in relation to the nitrite concentration in the media.
due to the lack of knowledge on the mechanism of e. coli natural membrane potential formation, the phenomenological function u, which qualitatively reproduces the experimental data received by motteram et al.
<cit> , was used in the model.
according to these data, the membrane potential highest value was observed in the  <dig> -1 mm nitrite concentration range and was stable starting from the  <dig>  mm nitrite concentrations, and decreased outside this range.
to reproduce this phenomenon, we used a modified formula  from the previous research  <cit> : us=skpmf,1hpmf,1+δ×skpmf,2hpmf,2+ω2×skpmf,3hpmf,31+skpmf,1hpmf,1+skpmf,2hpmf,2+ω1×skpmf,3hpmf, <dig>
further, we estimated parameter values for the given function  so that its behavior was qualitatively consistent with the above listed features.
qualitative behavior of the u function is presented in the fig.
<dig> fig.
<dig> the nominal membrane potential u in relation to the added nitrite level.
the membrane potential value is presented in arbitrary units
the system of differential equations describing the transport process is as follows: dnrfacdt=−ktnrf,cp,uus⋅nrfac−ktnrf,pc⋅nrfa,dnrfbcdt=−ktnrf,cp,uus⋅nrfbc−ktnrf,pc⋅nrfb,dnrfadt=δperiplktnrf,cp,uus⋅nrfac−ktnrf,pc⋅nrfa,dnrfbdt=δperiplktnrf,cp,uus⋅nrfbc−ktnrf,pc⋅nrfb, where nrfac and nrfbc – concentrations of nrfА and nrfВ proteins in the cytoplasm, nrfa and nrfb – concentrations of nrfА and nrfВ proteins in the periplasm, ktnrf,c
p – the rate constant for the nrfА and nrfВ proteins transport from the cytoplasm to the periplasm, ktnrf,p
c – the rate constant for the nrfА and nrfВ proteins transport from the periplasm to the cytoplasm, δperipl – ratio of cytoplasmic volume to periplasmic volume.
subsystem  describes the nirc transporter active form formation
as stated above, the structure of the e. coli nirc transporter catalytically active form is not known.
because both e. coli nirc and s. typhimurium nirc proteins belong to the same family, pentameric nirc protein of s. typhimurium was used for modeling the e. coli nirc  <cit> .
accordingly, the complex formation is described based on a fifth-order reaction nirc+nirc+nirc+nirc+nirc⇄kdis,nirc5kform,nirc5nirc <dig> which implies the following system of differential equations dnircdt=5kdis,c5nirc5−kdis,c5−4nirc <dig> dnirc5dt=−kdis,c5nirc5−kdis,c5−4nirc <dig>  where nirc and nirc <dig> – concentrations of nirc monomeric protein and nirc <dig> pentameric protein, kform,c <dig> kdic,c <dig> – the rate constants for the forward and reverse reactions, kdic,c5=kdic,c5kform,c <dig> – the dissociation constant of the nirc <dig> complex.
subsystem  describes the periplasmic nrf reductase active form formation
it is assumed in the model that under physiological conditions nrfa operates as part of the heterotetrameric nrfa2–nrfb <dig> complex  <cit> .
genes encoding nrfa and nrfb subunits of the enzyme are clustered in one operon and are translated in the e. coli cell in equimolar amounts  <cit> .
the nrf active form  concentration in the periplasm was calculated with the biochemical model  <cit> : nrfa+nrfb⇄kdis,abkform,abnrfab,nrfab+nrfab⇄kdis,a2b2kform,a2b2nrfa2b <dig>
this implies a system of differential equations that describes the nrfab and nrfa2b <dig> complexes formation and degradation in the periplasm, as well as nrfa and nrfb proteins degradation in the periplasm: dnrfadt=kdis,abnrfab−kdis,ab−1nrfa⋅nrfb,dnrfbdt=kdis,abnrfab−kdis,ab−1nrfa⋅nrfb,dnrfabdt=−kdis,abnrfab−kdis,ab−1nrfa⋅nrfb+2kdis,a2b2nrfa2b2−kdis,a2b2−1nrfab <dig> dnrfa2b2dt=−kdis,a2b2nrfa2b2−kdis,a2b2−1nrfab <dig>  where nrfa, nrfb – concentrations of free nrfa and nrfb monomeric proteins in the periplasm, nrfab –nrfab complex concentration in the periplasm, nrfa2b <dig> – concentration of the active heterotetrameric form  of the nrf enzyme, kform,ab, kform,a2b <dig> and kdis,ab, kdis,a2b2– the rate constants for the forward and reverse reactions, kdis,ab=kdis,abkform,ab, kdis,a2b2=kdis,a2b2kform,a2b2– the dissociation constants for the nrfab and nrfa2b <dig> complexes.
subsystem  describes the active form formation of the cytoplasmic nirb reductase
according to  <cit> , it is assumed in the model that under physiological conditions nirb nitrite reductase operates as part of nirb2nird the heterotrimeric complex.
the concentration of the nirb reductase active form was evaluated with the biochemical model: nirb+nirb⇄kdis,b2kform,b2nirb <dig> nirb2+nird⇄kdis,b2dkform,b2dnirb2d
on this basis, a system of differential equations describing the nirb <dig> and nirb2nird complexes formation and degradation as well as nirb and nird proteins degradation in the cytoplasm is as follows: dnirbdt=2kdis,b2nirb2−kdis,b2−1nirb <dig> dnirb2dt=−kdis,b2nirb2−kdis,b2−1nirb2+kdis,b2dnirb2d−kdis,b2d−1nirb2×nird,dnirddt=kdis,b2dnirb2d−kdis,b2d−1nirb2×nird,dnirb2ddt=−kdis,b2dnirb2d−kdis,b2d−1nirb2×nird, where nirb, nird –intracellular concentrations of the nirb and nird monomeic proteins, nirb <dig> – intracellular concentration of the nirb <dig> dimer, nirb2d – intracellular concentration of the cytoplasmic nitrite reductase active form, which is a nirb2d trimer, kdis,b <dig> kform,b <dig> kdis,b <dig> kform,b2– the rate constants for the forward and reverse reactions, kdis,b2=kdis,b2kform,b <dig>  kdis,b2d=kdis,b2dkform,b2d– the dissociation constants for the nirb <dig> and nirb2nird complexes.
subsystem  describes the degradation of proteins and their complexes as well as their dilution during cell growth
degradation processes of proteins and their complexes are described by the linear monomolecular reaction.
the dilution rate of the substance intracellular concentration due to the cell growth, when concentration is x, in general equals x, where v – current cell volume and v’ – cell volume growth rate.
because the growth rate in the chemostat is synchronized with the outflow rate, the average cell growth rate constant is equal to the outflow rate constant  <cit> ; as a result, we have the following system dxdt=kd,x+kflowx,x=w,nrfacnrfc,nrfa,nrfb,nirc,nirb,nird,nrfab,nrfa2b <dig> nirb <dig> nirb2d,nirc <dig>
disequilibrium model assemblage
to assemble a disequilibrium model from the equations describing – subsystems, the rate law was used, according to which the total change rate of a particular substance concentration equals to the change rates sum of a substance concentrations in each local process.
in this case, the local processes are the – subsystems, and the corresponding differential equations predetermine the local rates.
as a result, we have a system of differential equations , which is listed below.
<dig> dudt=kflows−u+c⋅kdis,nrfa2b2unrfa2b2u−kdis,nrfa2b2u−1nrfa2b2×u++kdis,nirc5nirc5−kdis,nirc5u−1nirc5×u+kcat,nircoutnirc5w,dwdt=kcat,nircinnirc5u+kdis,nirc5wnirc5w−kdis,nirc5w−1nirc5×w++kdis,nirb2dwnirb2dw−kdis,nirb2dw−1nirb2d×w−kflow⋅w,dnircdt=ksnirmniru+5kdis,nirc5nirc5−kdis,nirc5−4nirc5−kd,nirc+kflownirc,dnirc5dt=kcat,nircinnirc5u+kcat,nircoutnirc5w−kdis,nirc5nirc5−kdis,nirc5−4nirc5−−kdis,nirc5unirc5u−kdis,nirc5u−1nirc5×u−kdis,nirc5wnirc5w−kdis,nirc5w−1nirc5×w−kd,nirc5+kflownirc <dig> dnirc5udt=−kdis,nirc5wnirc5u−kdis,nirc5u−1nirc5×u−kcat,nircin+kflownirc5u,dnirc5wdt=−kdis,nirc5wnirc5w−kdis,nirc5w−1nirc5×w−kcat,nircout+kflownirc5w,dnirbdt=ksnirmniru+2kdis,nirb2nirb2−kdis,nirb2−1nirb2−kd,nirb+kflownirb,dnirddt=ksnirmniru+kdis,nirb2dnirb2d−kdis,nirb2d−1nirb2×nird−kd,nird+kflownird,dnirb2dt=kdis,nirb2nirb2−kdis,nirb2−1nirb2++kdis,nirb2dnirb2d−kdis,nirb2d−1nirb2×nird−kd,nirb2+kflownirb <dig> dnirb2ddt=−kdis,nirb2dnirb2d−kdis,nirb2d−1nirb2×nird+kcat,nirbnirb2dw++kdis,nirb2dwnirb2dw−kdis,nirb2dw−1nirb2d×w−kd,nirb2d+kflownirb2d,dnirb2dwdt=−kdis,nirb2dwnirb2dw−kdis,nirb2dw−1nirb2d×w−kcat,nirb+kflownirb2dw,dnrfacdt=ksnrfmnrfu−ktnrf,cp,uusnrfac−ktnrf,pcnrfa−kd,nrfac+kflownrfac,dnrfbcdt=ksnrfmnrfu−ktnrf,cp,uusnrfbc−ktin,pcnrfb−kd,nrfbc+kflownrfbc,dnrfadt=δperipl⋅ktnrf,cp,uusnrfac−ktnrf,pcnrfa+kform,nrfabkdis,nrfabnrfab−nrfa⋅nrfb−kd,nrfa+kflownrfa,dnrfbdt=δperipl⋅ktnrf,cp,uusnrfbc−ktnrf,pcnrfb+kform,nrfabkdis,nrfabnrfab−nrfa⋅nrfb−kd,nrfb+kflownrfb,dnrfabdt=2kdis,nrfa2b2nrfa2b2−kdis,nrfa2b2−1nrfab2−kdis,nrfabnrfab−kdis,nrfab−1nrfa⋅nrfb−kd,nrfab+kflownrfab,dnrfa2b2dt=kdis,nrfa2b2unrfa2b2u−kdis,nrfa2b2u−1nrfa2b2×u+kcat,nrfa2b2unrfa2b2u.−−kdis,nrfa2b2nrfa2b2−kdis,nrfa2b2−1nrfab2−kd,nrfa2b2+kflownrfa2b <dig> dnrfa2b2udt=−kdis,nrfa2b2unrfa2b2u−kdis,nrfa2b2u−1nrfa2b2×u−kcat,nrfa2b2u+kflownrfa2b2u.
the first equation in the  model describes the total change rate of the nitrite concentration in the chemostat, which is a sum of nitrite inflow and outflow rates into/from the chemostat, the rate of the nrfa, nirc and nitrite enzyme-substrate complexes formation and the rate of the nitrite export from e. coli cells.
since c is the total cell volume fraction in the chemostat volume, in the first equation rates of processes that occur within an individual cell are multiplied by the c value.
other equations describe processes that occur within an individual cell.
the second equation describes the change rate of the intracellular nitrite concentration, and other equations–formation processes of the nrfa and nirb reductases and nirc transporter active forms, transport and catalytic reactions.
estimation of the model parameters
the complete list of  model parameters and their values used in the calculations are represented in the appendix: table  <dig>  the procedure for obtaining parameter values for the mnrf, mnir and u functions was described above  and ).
in this section we will consider the  model parameters, the values of which were obtained on the basis of the existing experimental data on the dynamics of the e. coli cell culture growth under glucose-limiting conditions in the flow chemostat, the kinetics of the nitrite utilization by nrfa and nrfb nitrite reductases and parameters of the nirc protein mediated nitrite transport .
parameters of the protein passage from the cytoplasm to the periplasm under the membrane potential influence were estimated on the basis of the assumption that nrf enzyme is transported from the cytoplasm to the periplasm as nrfa and nrfb monomer subunits, which are synthesized in the cytoplasm ).
the ratio of the cytoplasmic volume to the periplasmic volume δperipl was assessed based on the experimental data, according to which the periplasmic space proportion makes 8–40 % of the whole e. coli cell volume .
on this basis, it is determined that the δperipl value may vary in the  <dig> – <dig>  range.
in the model δperipl is equal to six .
when evaluating the degradation rates of nrfА, nrfВ, nirb, nird, nirc proteins and their complexes, we considered the mean half-life of a pool of e. coli cytoplasmic proteins   <cit>  and stability of the periplasmic fraction of protein complexes, which is significantly higher  <cit> .
the rate constant values for the nrf reductase dimeric and trimeric forms dissociation were estimated based on the experimental data  <cit> .
the rate constant values for the nrfab, nrfa2b <dig>  nirc <dig>  nirb <dig> and nirb2d complexes formation were evaluated on the basis of indirect data on the protein-protein interaction kinetics  <cit> .
the e. coli cell culture volume fraction in the chemostat volume was assessed based on the previous data  <cit> , according to which when glucose concentration in the medium is fixated  and glucose inflow rate is constant, the cell growth rate is given and constant.
the cell cycle duration under such conditions is 70 min.
accordingly, the c parameter in the  model is assumed constant.
considering the described above parameters for the culture development in the chemostat and the rate of glucose uptake by the cells under similar limiting conditions  <cit> , the range of the c parameter possible variation was evaluated:  <dig> ≤c≤ <dig> .
the c =  <dig>  parameter value was taken in the model analysis.
values of other parameters that have not been evaluated on the basis of experimental data were estimated in the process of model numerical adaptation to the experimental data on the nitrite accumulation in the chemostat  <cit> .
the model analysis approach
the model analysis approach was selected based on the purpose of the study, which was to explain the kinetic data on the nitrite accumulation in the chemostat according to molecular-genetic mechanisms of nitrite utilization in e. coli cell cultured in a flow chemostat.
since the target experiment performed in  <cit>  consisted of a series of measurements of nitrite equilibrium concentrations, which were established in the chemostat at every given constant nitrite inflow rate, calculations of the positive stationary points were also performed in the  model.
positive stationary point for each given s value of the added nitrate was found numerically by solving the system of algebraic equations, which was obtained by equating the right-hand side of the  model to zero.
the mathematica program was used for all calculations.
positive solutions found in all calculations were the only ones in a positive change range of the system arguments.
RESULTS
assessment of the contribution of mechanisms regulating the activity of enzymes involved in no <dig> metabolism and transport to the control of no <dig> intracellular levels
as discussed above, the existing amount of experimental data on nrfa, nirb and nirС gene expression regulation and nitrite accumulation dynamics during e. coli stationary phase growth in the chemostat, presented by wang et al  <cit> , allowed us to develop a mathematical model of nitrite intracellular utilization and investigate the role of various components of the nitrite utilization system in controlling nitrite intracellular level.
after adapting the model to the experimental data, the  model calculation results, presented in fig.
<dig>  demonstrated a good agreement between theoretical curve and quantitative data on the nitrite accumulation dynamics in a chemostat .
the maximum value of the intracellular nitrite concentration in the studied range of the added nitrite did not exceed 14 μm , which is in qualitative agreement with the experimental data on the maintenance of sufficiently lower, less than  <dig>  mm, nitrite concentration in the cell under e. coli cell cultures anaerobic growth at 20 mm nitrate, when the nitrite concentration in the growth medium reaches 5 mm  <cit> .fig.
<dig> the steady-state nitrite concentration in a chemostat and inside a cell.
curve –  model calculation, dots – experimental data from wang et al.
<cit> , x-axis – added nitrite concentration , y-axis: а – steady-state nitrite concentration in the chemostat , b – nitrite concentration in a cell .
parameter values are presented in the appendix: table 1
using the  model calculations, various components of the nitrite utilization system in e. coli cells cultured in a flow chemostat at different concentrations of the added nitrite were evaluated for their contribution to the nitrite accumulation dynamics in the chemostat and inside the cell .fig.
<dig> contribution of the no <dig> utilization and transport system various components to the nitrite accumulation dynamics.
dots – experimental data on the nitrite utilization by the cell culture population, calculated with the kflow/c law based on the data from  <cit> ; curves 1- <dig> –  model calculation; curves 6- <dig> –  model calculation; curve  <dig> – the rate of nitrite utilization by the perilasmic nrfa reductase; curve  <dig>  – the rate of nitrite import by the nirc transporter; curve  <dig> – the rate of nitrite export by the nirc transporter, measured in negative units; curve  <dig> – the rate of nitrite utilization by the cytoplasmic nirb reductase; curve  <dig> – the total rate of nitrite utilization by nrfa and nirb reductases; curve  <dig> – the rate of nitrite utilization by the nrfА reductase, implemented in the  model, excluding the effect of the membrane potential  on the nira and nirb proteins diffusion rate kt
nrf,cp =  <dig>  s−1; curve  <dig> – difference between curves  <dig> and 6; curve  <dig> – the total rate of nitrite assimilation by the cell, calculated with .
other parameters values are presented in the appendix: table 1
we see the periplasmic nrfa reductase  and cytoplasmic nirb reductase  contributions to the nitrite utilization dynamics in the chemostat , the experimental parameters of which were calculated based on the data from  <cit>  and are shown with dots.
the nirc transporter activity dynamics in relation to the nitrite import into the cell  and its export from the cell  was demonstrated.
analysis of results has shown that the range of the added nitrite concentrations can be divided into two nominal range areas: the area of low values 0 ≤ s ≤ ~ <dig>  mm and the area of high values s > ~ <dig>  mm, in which nitrite is utilized only by the periplasmic or the cytoplasmic nitrite reductase, respectively.
this is due to the functional peculiarities of the nrf and nir operons, encoding these enzymes, the activity of which depends on the nitrite concentration  <cit> .
let us consider the area of high values of the added to the chemostat nitrite .
in this range there is a high nitrite processing efficiency of the nirb reductase , which is sufficient for describing the experimentally observed nitrite utilization rate in e. coli cells , and low nitrate processing efficiency of the nrfa reductase , which is consistent with the expression dynamics of encoding these enzymes operons  <cit> .
moreover, according to the model, the nitrite import rate  strictly corresponds to the rate of its processing by the nirb reductase , indicating the close association of the intracellular nitrite utilization system with the transport system of nitrite import into the cell at high substrate levels, and is a result of the arrangement of genes encoding their structure in one operon.
however, as seen in fig.
<dig>  with the parameter values taken , nitrite export from the cell into the medium by the nirc transporter is practically absent .
such property of the model, in our opinion, is a result of some parameter values taken during calculations  and is associated with the transport system, the evaluation of which was carried out based on the s. typhimurium nirc catalytic property data  <cit> .
according to the  <cit> , s. typhimurium nirc is active in the millimolar range of the nitrite concentration.
these data transfer to the e. coli cell system  puts practically no restrictions on the nitrite import process as steady-state nitrite concentration in the chemostat in the range of added nitrite >  <dig>  mm is higher than  <dig>  mm.
however, when using those parameters of nitrite utilization in the cell by the cytoplasmic nirb reductase, taken from  <cit> , in the model calculations, when s >  <dig>  mm, the intracellular nitrite concentration was in the range of  <dig> to 15 μm .
it is not surprising that at such nitrite concentrations in the cell the export activity of the enzyme is practically absent.
it should be noted that the parameters of the nitrite export from the cell  indirectly indicate the concentration threshold of nitrite toxicity for cells.
therefore, if the value kmnirc,out = 1 mm is valid for e. coli cells, the intracellular nitrite becomes toxic for e. coli cells only at concentrations of about 1 mm.
based on this, under experimental conditions, implemented in  <cit> , the cell functions far from the toxicity threshold and has no need for nitrite export from the cell.
however, experimental data indicate the significant nirc transporter export activity  <cit> , so the need for kmnirc,out =1мМ value to explain the data on the nitrite accumulation in the chemostat  <cit>  in the range of high added nitrite concentration remains an open question.
we will consider this issue below.
let us now consider the range of low added nitrite concentrations .
as stated above, in this range, the inflowing nitrite is utilized only by the periplasmic nrfa reductase .
high nrf reductase activity is ensured by the combined effect, which is achieved through the activation of nrf operon expression , the de novo synthesis of nrfa and nrfb proteins, and through the increased rate of nrfa and nrfb proteins passage from the cytoplasm to the periplasm as a result of the membrane potential u formation .
let us consider the potential contribution of each mechanism in the formation of the nitrite utilization rate in e. coli cells cultured in the chemostat.
contribution evaluation of the genetic component of the nrfa periplasmic reductase activity regulation mechanism to the no <dig> utilization in e. coli cells
to assess the contribution of the nrf reductase encoding nrf operon expression regulating mechanisms to the nitrite utilization in the chemostat at micromolar concentrations of the substrate we explored the  model functioning dynamics, ignoring mechanisms of membrane potential formation.
to do this, we put value du =  <dig>  excluding subsystem  <dig> from the model .
in this version of the model, to achieve best agreement between calculations and experimental data, we had to increase the rate constant value of the constitutive nrfa and nrfb proteins transport from the cytoplasm to the periplasm, compared to the baseline, by  <dig>  times .
the results are shown in fig.
<dig> .
it is well seen that the nrfa reductase activity obtained by calculating the  model, in which the nrfa activity is determined only by the nitrite-dependent genetic mechanisms of the nrf operon expression regulation  at 1 mm point, is significantly lower than that obtained by calculating the total  model .
the difference between these activities  is a bell curve, which corresponds to the added activity, which provides an adequate description of the experimental data on nitrite utilization in the chemostat by the  model  as shown in fig.
<dig> with dots.
this added activity for describing the same data is missing in the  model .
thus, an adequate description of the nitrite accumulation kinetics in the chemostat in the micromolar range of the added nitrite concentration  cannot be achieved on the basis of the genetic data of nitrite-dependent nrf operon expression activity  <cit> .
from here, we conclude that membrane potential plays an important role in the nitrite utilization dynamics in the micromolar range of the substrate concentration in the chemostat.
above, we have given calculations that showed that it is sufficient to consider the membrane potential impact on the rate of the enzyme redistribution between the cytoplasm and the periplasm for explaining the experimental data.
however, since there is evidence of other possibilities of membrane potential influence on the nrfa reductase activity, it is necessary to consider them as alternative mechanisms for explaining the data on nitrite consumption in the chemostat.
these mechanisms are discussed in the next section.
mechanisms of influence of the membrane potential on the periplasmic nrfa reductase activity
we should note that, according to the existing data, membrane potential plays a significant role in the formation of catalytically active nrf enzyme molecule, affecting its passage into the periplasm, its stability and catalytic properties.
it is known that protein passage into the periplasm is initiated by the electric membrane potential advent as a result of the activity of respiratory enzymes involved in the proton gradient formation  <cit> .
nrf reductase refers specifically to such enzymes and nitrite is an electron acceptor  <cit> .
proteins passage into the periplasm affects their stability  <cit> .
enzyme assembling and correct orientation in the periplasmic space is also dependent on the presence of membrane potential  <cit> .
moreover, it has been shown that michaelis constant for the nrf enzyme is dependent on the membrane potential value  <cit> , and the membrane potential value is dependent on the nitrite concentration  <cit> .
to answer the question of which of the above mechanisms is essential to explain the experimental data, we evaluated the each mechanism contribution to the nitrite utilization at low added nitrite levels in the chemostat.
first, we tested if the enzyme catalytic properties change under the influence of the membrane potential can affect the nitrite utilization kinetics in the chemostat.
it is known that michaelis constant for the nrfa enzyme varies depending on the membrane potential   <cit> .
analysis of the  model have shown that the km,nrf value change by an order of magnitude or more, that goes beyond the experimentally observed differences  <cit> , leads to minor changes  of the nrfa mediated nitrite utilization rate in the examined range and cannot be a source of additional enzyme activity at low nitrite concentrations, sufficient to explain the experiments  <cit> .
introduction to the  model a positive non-linear relationship between the nrfa reductase catalytic activity and nitrite concentration in the medium allows to achieve the desired effect, however, this hypothesis is not supported by the data on a simple reaction kinetics described by the michaelis-menten equation  <cit> .
regarding the membrane potential impact on the rate of the enzyme passage into the periplasm and catalytically active enzyme molecule assembly, the consideration of this process is consistent with experimental data and, according to the  model analysis presented above, is sufficient to describe the experimental data on nitrite accumulation dynamics in the chemostat  <cit> .
it should be noted that passage of enzyme subunits into the periplasm under the membrane potential influence is also accompanied by an increase in their local concentrations and stability.
how important are these processes to explain the experimental data?
the local concentration coefficient of variation is equal to the ratio of cytoplasmic volume to periplasmic volume, which, based on the experimental data , is equal to six in the model .
in essence, the mechanism of membrane potential is to increase the rate of the nrfa and nrfb proteins transport from the cytoplasm to the periplasm and provide the accumulation of proteins synthesized in the cytoplasm, in a smaller volume of periplasm.
the δperipl =  <dig> ratio used in a model is sufficient to ensure that the concentration of nrf enzyme in the periplasm reached a level ensuring adequate nitrite consumption rate in the chemostat, meeting the experimental data.
note that the enzyme concentration increase is equivalent to the enzyme catalytic activity value increase , as during the catalytic process the absolute reaction rate is determined by multiplying the catalytic constants by the concentration of the enzyme active molecules.
theoretically, it is also possible that the increase of protein stability in the periplasm by almost an order of magnitude  <cit> , which is also considered in the model, may be a crucial factor.
however, according to the calculations, this factor is of secondary importance, and it is not required to explain the experimental data.
thus, the model analysis demonstrates that membrane potential significantly contributes to the regulation of the periplasmic nrfa reductase activity and it must be considered in order to explain the nitrite utilization dynamics at substrate concentrations <1 mm.
the most probable mechanism is the rate change of nrfa and nrfb proteins passage from the cytoplasm to the periplasm, depending on the membrane potential value.
in light of these theoretical results, attention is drawn to the fact that both mechanisms: the nitrite-dependent membrane potential change and the nitrite-dependent effectiveness of nrf operon expression, have qualitatively identical unimodal activity curves .
as a result, the necessity of addressing both mechanisms to explain the experimentally observed nitrite accumulation curve in the low range of the added nitrite concentration emerges.
this question is discussed in the next section.
on the functional redundancy of the no <dig> utilization mechanisms in e. coli cell
in this section, we analyze the possibility of  model functional redundancy in the utilization mechanisms and related issues.
the essence of the matter derives from the fact that subsystem –, described above and including a certain number of genetic, transport, complex-forming, and enzymatic reactions, make up the structural basis of the model.
we used two approaches to describe the subsystems.
regulation of the nrf and nir operons expression efficiency, as well as the membrane potential, we described with phenomenological functions of the class of generalized hill functions  <cit> .
other processes  were described with non-equilibrium biochemical reactions.
meaning the right-hand side of  system includes members that describe both direct and reverse non-equilibrium processes.
at the same time, to describe the dynamic characteristics of the nitrite accumulation in the chemostat, the  system equilibrium state has to be calculated.
it is well known from the general theory of dynamical systems that the equilibrium state specific the value is not determined by the values of the parameters, but by values of their specific combinations, the number of which is usually less than the number of the parameters themselves.
the simplest example is the bimolecular reaction.
to describe the reaction far from equilibrium, it is necessary to know the rate constant values of the forward and reverse reactions.
however, at equilibrium, it is sufficient to know the value of one equilibrium constant, which is the ratio of non-equilibrium constants.
as a result, we get a traditional problem of determining the degree of functional and parametric redundancy of mechanisms reflected in , which are necessary to describe the nitrite consumption kinetics in the chemostat.
in this context, we consider the redundancy issue of mechanisms, by which the nitrite accumulation kinetics in the chemostat at concentrations of the added nitrite s < 2 mm is explained in the model.
as stated above, the apparent similarity of the curves describing the nrf operon expression and membrane potential influence on the rate of the nrfa and nrfb proteins redistribution in relation to the added nitrite concentration  suggest that in order to explain the experimental data from  <cit>  it is sufficient to use only one of these mechanisms and the other can be omitted.
is it so?
as for the need to consider the membrane potential, we showed above that without considering the influence of membrane-dependent potential on the rate of the nrfa and nrfb proteins redistribution between the cytoplasm and the periplasm it is impossible to adequately describe the experimental data on nitrite utilization in the chemostat  <cit> .
this output, as shown above, is considerably based on the data on nitrite-dependent genetic regulation of nrf operon expression .
however, let’s assume that nrf operon is constitutively expressed.
it appears that for this version of  model it is not difficult to pick up a set of parameter values  ≡  <dig>  mnir ≡  <dig>  ks,nrfa = ks,nrfb =  <dig>  mm/s, knrf1 =  <dig> , hnrf,1 =  <dig>  δnrf,2 =  <dig> , knr,f2 =  <dig> , hnrf,2 =  <dig>  ω2 =  <dig> , knrf,3 =  <dig> , hnrf,3 =  <dig>  values of other parameters are unchanged and presented in the appendix: table 1), for which the model will adequately describe the nitrite accumulation curve.
that is, the data on the genetic regulation of the nrf operon expression are redundant for explaining the nitrite accumulation curve in relation to the added nitrite concentration.
it is sufficient to assume the nrf operon constitutive expression and take into account the membrane potential impact on the nira and nirb proteins redistribution between the cytoplasm and the periplasm.
however, this conclusion is theoretical, since the native e. coli cells reveal a strong nitrite-dependent regulation of the nrf operon expression  <cit> .
and, nevertheless, the revealed redundancy raises important theoretical question of evolutionary expediency of such a complex mechanism of nrf operon expression regulation in e. coli cells: because the cell could easily utilize nitrite in simpler and more economical way!
this issue requires further study and more detailed analysis.
the second issue considered in this section, is the redundancy of the nitrite export mechanism.
it is indeed known that e. coli cells have a certain threshold value of the intracellular nitrite concentration, above which it becomes toxic to the bacteria.
it is believed that in this case the cell has a protective system based on the nirc nitrite export activity, which is basic during nitrite breathing  <cit> .
however, the parameters of such activity were not measured.
in the model, we used values of the parameters that have been measured for the s. typhimurium nirc transporter.
but then, it implies that  the nitrite concentration toxic for e. coli cells is comparable to millimolar concentrations, and  in the  model nitrite export from cells is practically absent .
however, there is reason to believe that in e. coli cells nitrite toxic concentration is in the micromolar range  <cit> , i.e., michaelis constant for the nirc mediated nitrite export can be much smaller than the one given in appendix: table  <dig>  on this basis, we have considered a model version, in which the value of the michaelis constant for the nirc mediated nitrite export from the cell is  <dig> times smaller: km,nirc,out =  <dig> μМ.
the corresponding calculation is presented in fig.
<dig> fig.
<dig> nitrite utilization and transport rate in relation to the concentration of the added nitrite.
dots – the rate of nitrite consumption by the cell, calculated based on the experimental data from  <cit> ; curve  <dig> – the rate of nitrite utilization by the periplasmic nrfa reductase; curve  <dig> – the rate of nitrite import into the cell by the nirc transporter; curve  <dig> – the rate of nitrite export from the cell by the nirc transporter, measured in negative units; curve  <dig> – the rate of nitrite utilization by the cytoplasmic nirb reductase; curve  <dig> – the total rate of nitrite utilization by nrfa and nirb reductases.
calculations performed with the  model version, in which k
М,nirc,out =  <dig>  mМ, k
s,nirc,0 =  <dig>  mm/s; other parameter values are from the appendix: table 1
calculations have shown that such parametric version of the model does not impair the data reproduction quality.
at the same time, under a given set of parameter values there is a high rate of nitrite export from the cell into the chemostat .
consequently, the ratio of nirc transporter import activity functional parameter to its export activity functional parameter is a free parameter in the  model, and to determine its value, additional information that is not available in the scientific literature is needed.
in conclusion, the nitrite concentration curve at different values of the added nitrite concentration, shown in fig.
4b, in a sense, is nominal.
it is easy to choose sets of parameters for which the model calculations of the intracellular nitrite concentration will be different from the above calculation in orders of magnitude, and, at the same time, the accuracy of the nitrite accumulation approximation curve will not be altered , so this characteristic is also free in the model and for its more accurate characterization it is necessary to have additional data, which is not yet available.
CONCLUSIONS
in the present work we introduce the model of nitrite utilization in e. coli cells cultured in the flow chemostat.
the parameter values were chosen in such a way that the model adequately describes the experimental curve of nitrite accumulation in the chemostat  <cit> .
the model analysis revealed the existence of two ranges of the added no <dig> concentration: low  and high .
according to the model, in the high concentration range the nitrite is substantially utilized by the nitrite transport/utilization system .
in this range, the nitrite-dependent genetic regulation of the nir operon expression is sufficient for explaining experimental data on the nitrite accumulation in the chemostat  <cit> .
on the contrary, in the low concentration range of the added nitrite the known genetic mechanisms of the nrf operon expression are not enough for explaining experimental data on the nitrite accumulation in the chemostat  <cit> .
analysis of different hypothesis has shown that the most probable additional mechanism of the periplasmic nrf reductase activity regulation in the micromolar concentration range of the added nitrite, that allows co-ordination of physiological and genetic data  <cit> , is a local enzyme concentration change due to its passage from the cytoplasm to the periplasm under the influence of nitrite concentration-dependent membrane potential.
in the framework of our model, the potential-dependent mechanism of the nrfa and nrfb proteins redistribution between cytoplasm and periplasm is the essential element required to explain observed nitrite accumulation dynamics in the chemostat in the low concentration range of the substrate.
analysis of different nitrite utilization and transport mechanisms contribution to the nitrite accumulation in the chemostat and inside the cell revealed some redundancy of molecular-genetic mechanisms of nrfa, nrfb reductases and nirc transporter activities regulation.
thus, to describe the nitrite accumulation curve in chemostat in the entire concentration range of the added nitrite it is sufficient to assume a constitutive nrf operon expression, whereas the nirc transporter and nirb nitrite reductase can be ignored.
in the high concentration range of the added nitrite the ratio of nirc transporter import activity functional parameter to its export activity functional parameter as well as intracellular nitrite concentration are free parameters of the model and may vary over a wide range without affecting the nitrite accumulation curve description accuracy of the model.
appendix
k
flow
δnrf,1
k
nrf,1
h
nrf,1
δnrf,2
k
nrf,2
h
nrf,2
ks
nrf
ks
nirc
ksnirb==ksnird
ksnirb==ksnirc/5
k
d,nrf
kt
nrf,cp
kt
nrf,pc
k
dis,nrfab
k
dis,nrfab
kdis,nrfa2b2
kdis,nrfa2b2
δ
peripl
d
u
δ
k
pmf,1
h
pmf,1
k
pmf,2
h
pmf,2
ω
<dig>  ω
2
с
k
pmf,3
с
h
pmf,3
с
k
d,nrfa
sec –1
k
d,nrfa
k
d,nrfab
kd,nrfa2b2
kcat,nrfa2b2u
kdis,nrfa2b2u
δnir,1
k
nir,1
h
nir,1
δnir,2
k
nir,2
h
nir,2
kdis,nirb2
kdis,nirb2d
kcat,nirb2dw
kdis,nirb2dw
kdis,nirc5
k
cat,nircin
kdis,nirc5u
k
cat,nircout
kdis,nirc5w
k
d,nirc
k
d,nirc5
k
d,nirb
k
d,nird
k
d,nirb2
k
d,nirb2d
c
a – mm , sec , if dimension is not specified, the quantity is dimensionless
b – parameter value was estimated based on the experimental data
c – members, before which the ω
<dig> and ω
<dig> coefficients are standing in the membrane potential formula, are zeroed in  model, the formula used in the model calculation is with nrf operon constitutive expression and enabled nirc and nirb enzymes functions
competing interests
the authors declare that they have no competing interests.
authors’ contributions
tmk conceived the study.
val and nar developed the model and analyzed the data.
tmk and nar evaluated the model parameters.
val and nar performed the model calculations.
tmk wrote the manuscript with editing by val and nar.
all authors read and approved the final manuscript.

 during anaerobic growth e. coli uses different electron acceptors in the electron transport chain, including nitrite. given the high toxicity of nitrite, expression regulation of genes involved in the nitrite-associated electron transport chain is closely linked to the expression regulation of nitrite metabolism genes. nitrite reductases nrfa and nirb are two main components of the nitrite-reducing system. these enzymes have different metabolite activities and localizations in the cell. nrfa reductase is a respiratory enzyme localized in the cell periplasm; it takes part in the proton gradient generation and uses no as an electron acceptor. nirb reductase is restricted to the cytoplasm and, most probably, has excessive nitrite detoxification as its main role in the cell. both enzymes catalyze the reduction of no to ammonium, but nrfa is most active at low nitrite levels, whereas nirb activity is observed only at high nitrite levels. differential expression of nrf and nir operons, encoding nrfa and nirb reductases, respectively, enables such a combination of nrfa and nirb reductase activities. the nrf and nir operons expression is controlled by transcription factors, narl and narp, the activity of which depends on narq and narx kinases. unlike nir operon, the expression of which is always positively regulated by nitrite, the nrf operon has regulatory region that allows its expression activation at low nitrite levels and expression inhibition at high nitrite levels. as a result, at nitrite concentrations in the media > 2 mm, genetic system that encodes nrfa reductase and allows nitrite metabolism in the periplasm of the cell switches to a more effective in terms of nitrite recycling genetic system, that encodes nirc transporter and nirb nitrite reductase and allows nitrite transport into the cell and its further recycle by the cytoplasmic reductase. due to a high toxicity of nitrite, cell transport system plays a crucial role in the nitrite level regulation in the cell. three transport proteins involved in the nitrite transmembrane passage are known: nark, naru and nirc, but the nitrite transport protein nirc is known to have the highest activity among three. localization of nirС gene, encoding the nirc protein, in the same operon as nirb gene demonstrates a tight connection between intracellular nitrite utilization and nitrite import/export systems at high nitrite levels; however, complex interplay between periplasmic and cytoplasmic systems, which determines the nitrite accumulation and utilization in the cell as well as its export from the cell, represents a significant knowledge gap. the existing amount of qualitative and quantitative data on nrfa, nirb and nirС genes expression regulation and nitrite accumulation dynamics during e. coli stationary phase growth in chemostat, allowed us to develop a mathematical model of nitrite intracellular utilization and investigate how nirb and nrfa reductases and nirc transporter contribute to this process. using this model, mechanisms involved in regulation of these proteins’ activities were investigated for their roles in controlling the intracellular nitrite levels. mathematical model of nitrite utilization in e. coli cells the exact conditions of e. coli cells cultivation in a nitrite-supplied flow chemostat are reproduced in the model. the model is adapted to the experimental data on nitrite accumulation dynamics in a chemostat and on the inducer levels dependent nrf and nir operons expression. a steady-state growth rate and cell culture density in the chemostat was achieved by fixating the glucose uptake. the model of nitrite utilization in e. coli cells cultured in a chemostat consists of several elementary subsystems. subsystems are detailed in accordance with the process scheme depicted in fig.  and are described below. when describing these processes, generalized hill functions were applied.fig. scheme of the processes occurring in the chemostat during nitrite utilization in e. coli cells. – nitrite inflow and outflow at steady-state flow rate; – reduction of extracellular nitrite to ammonium mediated by periplasmic nrfa reductase;, – import of the extracellular nitrite into the cell and its export from the cell into the chemostat mediated by the nirc transporter; – intracellular nitrite utilization mediated by cytoplasmic nirb reductase; – nrfa and nrfb proteins synthesis from the nrf operon mrna and nirb, nirc and nird synthesis from the nir operon mrna; – transport of the nrf reductase a and b subunits from the cytoplasm to the periplasm and back; – formation of the nirc transporter active pentameric form; – formation of the nirb reductase active heteromeric form; – formation of the nrf active heterotetrameric form; – degradation of proteins and their complexes elementary subsystems of the model subsystem describes the nitrite flow through the chemostat nitrite inflow rate into the chemostat is a constant parameter throughout the experiment and nitrite outflow rate from the chemostat is proportional to its current concentration, so the total nitrite flow rate through the chemostat can be described with the equation dudt=vin−kflow⋅u, where u – extracellular nitrite concentration in the chemostat, vin – nitrite inflow rate into the chemostat, kflow – the rate constant for the nitrite outflow. later on, for convenience, vin parameter will be presented as vin = kflow ⋅ s, where s is a nitrite concentration established in the chemostat at a given vin and kflow values in the absence of the cell culture. later on we will name this parameter «added nitrite». subsystem describes the reduction of extracellular nitrite to ammonium mediated by the periplasmic formate-dependent nrfa reductase the reduction of extracellular nitrite to ammonium by periplasmic formate-dependent nrfa reductase is described in accordance with the stoichiometric scheme: no2− + 6е− + 7h+ → nh3 + 2h2o with the simplest michaelis–menten model dudt=−ckdis,nrfa2b2uu×nrfa2b2/km,nrf−nrfa2b2u,dnrfa2b2dt=kdis,nrfa2b2uu×nrfa2b2/km,nrf−nrfa2b2u−kcat,nrfnrfa2b2u,dnrfa2b2dt=−kdis,nrfa2b2uu×nrfa2b2/km,nrf−nrfa2b2u+kcat,nrfnrfa2b2u, where nrfa2b – nitrite-dependent nrfa reductase concentration in the cell periplasm ), c – cell volume fraction in the chemostat volume, kcat,nrf – catalytic constant of nrfa reductase turnover, km,nrf – the michaelis constant. subsystem describes the extracellular nitrite import into the cell it is known that extracellular nitrite transport from the environment into the cell is mediated by the nirc, naru и nark transporter proteins, but due to significantly higher nirc nirtrite transport activity, compared with naru and nark, only nirc transporter function was considered in the model. the structure of the e. coli nirc transporter catalytically active form is not known. we assume that it should be similar to the salmonella typhimurium nirc structure because both proteins belong to one family; that is why in the model nirc is assumed to be a pentameric protein. the mechanism of nitrite import by the nirc transporter is also not known, therefore we described this process rate with the simplest system, resulted from the simplest michaelis–menten model dudt=−ckdis,nirc5uu×nirc5/km,nircin−nirc5u,dnirc5dt=kdis,nirc5uu×nirc5/km,nircin−nirc5u−kcat,nircinnirc5u,dnirc5udt=−kdis,nirc5uu×nirc5/km,nircin−nirc5u+kcat,nircinnirc5u, where nirc – nitrite-dependent concentration of the nirc transporter active form in the cell periplasm, kdis,nirc5u – the rate constant for the u × nirc complex dissociation into subunits, kcat,nircin – the rate constant for the nitrite import by the nirc transporter, km,nircin – the michaelis constant. subsystem describes intracellular nitrite export from the cell to the chemostat in the model intracellular nitrite export from the cell to the chemostat, as well as its import, is mediated only by the nirc transporter. the naru and nark transporters contribution to this process was not considered. because the mechanism of nitrite import by the nirc transporter is not yet studied, the rate of this process is described with the simplest equation dwdt=−kdis,nirc5ww×nirc5/km,nircout−nirc5w,dnirc5dt=kdis,nirc5ww×nirc5/km,nircout−nirc5w−kcat,nircoutnirc5w,dnirc5wdt=−kdis,nirc5ww×nirc5/km,nircout−nirc5w+kcat,nircoutnirc5w, where w –intracellular nitrite concentration, kdis,nirc5w – the rate constant for the w × nirc complex dissociation into subunits, kcat,nircout – the rate constant for the nitrite export by the nirc transporter, km,nircout – the michaelis constant. subsystem describes the intracellular nitrite utilization mediated by the nadn-dependent nirb reductase in the model this process is described in accordance with the following stoichiometric scheme: no2−+3nadh+5h+→nh4++3nad++2h2o. the utilization rate is calculated with the michaelis–menten law dwdt=−kdis,nirb2dww×nirb2d/km,nirb−nirb2dw,dnirb2ddt=kdis,nirb2dww×nirb2d/km,nirb−nirb2dw−kcat,nirbnirb2dw,dnirb2dwdt=−kdis,nirc5ww×nirc5/km,nircout−nirc5w+kcat,nirbnirb2dw, where nirb2d – nitrite-dependent concentration of the nirb reductase active form, kdis,nirb2dw – the rate constant for the w × nirb2d molecular complex dissociation into subunits, kcat,nirb – the recycle catalytic constant, km,nirb – the michaelis constant. subsystem describes nrfa, nrfb, nirb, nirc и nird proteins synthesis in e. coli cells nrfa and nrfb proteins are synthesized from mrna templates of the nrf operon and nirb, nirc and nird proteins – from mrna templates of the nir operon. the nrf operon mrna concentration in the cell that exists in the chemostat at a current nitrite concentration u is assumed to be proportional to the relative nrfa-β-gal chimeric protein activity observed at the same nitrite concentration; and the nir operon mrna concentration in the same cell is assumed to be proportional to the relative nirb-β-gal chimeric protein activity, measured by wang and co-authors.fig. effect of nitrite concentration on nrfabcd and nirbdc operons expression during anaerobic cell growth. dots represent experimental values of activities of the nrfa-β-gal and nirb-β-gal chimeric proteins measured in and approximating their theoretical curves m nrf and m nir. the x-axis – steady-state nitrite concentration in the chemostat; the y-axis – the nrfa-β-gal chimeric protein relative activity, – the nrfb-β-gal chimeric protein relative activity. the steady-state nitrite concentration in the chemostat used as a scale bar on the x-axis was calculated earlier, and based on the data from the previous research activities of nrfa-β-gal and nirb-β-gal chimeric proteins are dependent on the nitrite concentration and were approximated with generalized hill functions, mnrfu=1+δnrf,1×uknrf,1hnrf,11+uknrf,1hnrf,1×1+δnrf,2×uknrf,2hnrf,21+uknrf,2hnrf, mniru=1+δnir,1×uknir,1hnir,1+δnir,2×uknir,2hnir,21+uknir,1hnir,1+uknir,2hnir, optimal parameters of the mnrf and mnir functions were estimated with the gradient descent algorithm and are represented in the appendix: table as a result, synthesis of nrfa and nrfb proteins can be described with the following system of differential equations: dnrfacdt=ksnrf,ac⋅mnrfu,dnrfbcdt=ksnrf,bc⋅mnrfu, where nrfac and nrfbc – concentrations of nrfa and nrfb proteins in the cytoplasm, ksnrf – the rate constant for the synthesis of nrfa and nrfb proteins, and synthesis of nirb, nirc and nird proteins can be described with the following system of equations: dnirbdt=ksnir,b⋅mniru,dnircdt=ksnir,c⋅mniru,dnirddt=ksnir,d⋅mniru, where nirb, nirc, nird – concentrations of nirb, nirc and nird proteins in the cytoplasm, ksnir – the rate constant for the protein synthesis. it is considered that the nirc protein translation level is more than five times higher than such of the nirb protein in the millimolar range of nitrite concentration subsystem describes transport of the nrfa and nrfb subunits of the nrf nitrite reductase from the cytoplasm to the periplasm the nrf enzyme is assumed to be transported from the cytoplasm to the periplasm as monomer subunits nrfa and nrfb, which are synthezised in the cytoplasm. we also assume that the specific rate constant for the nrf reductase a and b monomers transport from the cytoplasm, where their synthesis takes place, to the periplasmic space, where the active enzyme is being formed, is a variable that depends on the membrane potential: ktnrf ⋅ ku). generalized hill functions is used for the ku description ktnrf,cp,uus=ktnrf,cp1+duus, whose value is made dependent on the u function, named the nominal potential. in its turn, the u function mimics the membrane potential that changes in relation to the nitrite concentration in the media. due to the lack of knowledge on the mechanism of e. coli natural membrane potential formation, the phenomenological function u, which qualitatively reproduces the experimental data received by motteram et al., was used in the model. according to these data, the membrane potential highest value was observed in the -1 mm nitrite concentration range and was stable starting from the  mm nitrite concentrations, and decreased outside this range. to reproduce this phenomenon, we used a modified formula from the previous research: us=skpmf,1hpmf,1+δ×skpmf,2hpmf,2+ω2×skpmf,3hpmf,31+skpmf,1hpmf,1+skpmf,2hpmf,2+ω1×skpmf,3hpmf, further, we estimated parameter values for the given function so that its behavior was qualitatively consistent with the above listed features. qualitative behavior of the u function is presented in the fig.  fig. the nominal membrane potential u in relation to the added nitrite level. the membrane potential value is presented in arbitrary units the system of differential equations describing the transport process is as follows: dnrfacdt=−ktnrf,cp,uus⋅nrfac−ktnrf,pc⋅nrfa,dnrfbcdt=−ktnrf,cp,uus⋅nrfbc−ktnrf,pc⋅nrfb,dnrfadt=δperiplktnrf,cp,uus⋅nrfac−ktnrf,pc⋅nrfa,dnrfbdt=δperiplktnrf,cp,uus⋅nrfbc−ktnrf,pc⋅nrfb, where nrfac and nrfbc – concentrations of nrfА and nrfВ proteins in the cytoplasm, nrfa and nrfb – concentrations of nrfА and nrfВ proteins in the periplasm, ktnrf,c  p – the rate constant for the nrfА and nrfВ proteins transport from the cytoplasm to the periplasm, ktnrf,p  c – the rate constant for the nrfА and nrfВ proteins transport from the periplasm to the cytoplasm, δperipl – ratio of cytoplasmic volume to periplasmic volume. subsystem describes the nirc transporter active form formation as stated above, the structure of the e. coli nirc transporter catalytically active form is not known. because both e. coli nirc and s. typhimurium nirc proteins belong to the same family, pentameric nirc protein of s. typhimurium was used for modeling the e. coli nirc. accordingly, the complex formation is described based on a fifth-order reaction nirc+nirc+nirc+nirc+nirc⇄kdis,nirc5kform,nirc5nirc which implies the following system of differential equations dnircdt=5kdis,c5nirc5−kdis,c5−4nirc dnirc5dt=−kdis,c5nirc5−kdis,c5−4nirc where nirc and nirc – concentrations of nirc monomeric protein and nirc pentameric protein, kform,c kdic,c – the rate constants for the forward and reverse reactions, kdic,c5=kdic,c5kform,c – the dissociation constant of the nirc complex. subsystem describes the periplasmic nrf reductase active form formation it is assumed in the model that under physiological conditions nrfa operates as part of the heterotetrameric nrfa2–nrfb complex. genes encoding nrfa and nrfb subunits of the enzyme are clustered in one operon and are translated in the e. coli cell in equimolar amounts. the nrf active form concentration in the periplasm was calculated with the biochemical model: nrfa+nrfb⇄kdis,abkform,abnrfab,nrfab+nrfab⇄kdis,a2b2kform,a2b2nrfa2b this implies a system of differential equations that describes the nrfab and nrfa2b complexes formation and degradation in the periplasm, as well as nrfa and nrfb proteins degradation in the periplasm: dnrfadt=kdis,abnrfab−kdis,ab−1nrfa⋅nrfb,dnrfbdt=kdis,abnrfab−kdis,ab−1nrfa⋅nrfb,dnrfabdt=−kdis,abnrfab−kdis,ab−1nrfa⋅nrfb+2kdis,a2b2nrfa2b2−kdis,a2b2−1nrfab dnrfa2b2dt=−kdis,a2b2nrfa2b2−kdis,a2b2−1nrfab where nrfa, nrfb – concentrations of free nrfa and nrfb monomeric proteins in the periplasm, nrfab –nrfab complex concentration in the periplasm, nrfa2b – concentration of the active heterotetrameric form of the nrf enzyme, kform,ab, kform,a2b and kdis,ab, kdis,a2b2– the rate constants for the forward and reverse reactions, kdis,ab=kdis,abkform,ab, kdis,a2b2=kdis,a2b2kform,a2b2– the dissociation constants for the nrfab and nrfa2b complexes. subsystem describes the active form formation of the cytoplasmic nirb reductase according to, it is assumed in the model that under physiological conditions nirb nitrite reductase operates as part of nirb2nird the heterotrimeric complex. the concentration of the nirb reductase active form was evaluated with the biochemical model: nirb+nirb⇄kdis,b2kform,b2nirb nirb2+nird⇄kdis,b2dkform,b2dnirb2d on this basis, a system of differential equations describing the nirb and nirb2nird complexes formation and degradation as well as nirb and nird proteins degradation in the cytoplasm is as follows: dnirbdt=2kdis,b2nirb2−kdis,b2−1nirb dnirb2dt=−kdis,b2nirb2−kdis,b2−1nirb2+kdis,b2dnirb2d−kdis,b2d−1nirb2×nird,dnirddt=kdis,b2dnirb2d−kdis,b2d−1nirb2×nird,dnirb2ddt=−kdis,b2dnirb2d−kdis,b2d−1nirb2×nird, where nirb, nird –intracellular concentrations of the nirb and nird monomeic proteins, nirb – intracellular concentration of the nirb dimer, nirb2d – intracellular concentration of the cytoplasmic nitrite reductase active form, which is a nirb2d trimer, kdis,b kform,b kdis,b kform,b2– the rate constants for the forward and reverse reactions, kdis,b2=kdis,b2kform,b kdis,b2d=kdis,b2dkform,b2d– the dissociation constants for the nirb and nirb2nird complexes. subsystem describes the degradation of proteins and their complexes as well as their dilution during cell growth degradation processes of proteins and their complexes are described by the linear monomolecular reaction. the dilution rate of the substance intracellular concentration due to the cell growth, when concentration is x, in general equals x, where v – current cell volume and v’ – cell volume growth rate. because the growth rate in the chemostat is synchronized with the outflow rate, the average cell growth rate constant is equal to the outflow rate constant; as a result, we have the following system dxdt=kd,x+kflowx,x=w,nrfacnrfc,nrfa,nrfb,nirc,nirb,nird,nrfab,nrfa2b nirb nirb2d,nirc disequilibrium model assemblage to assemble a disequilibrium model from the equations describing – subsystems, the rate law was used, according to which the total change rate of a particular substance concentration equals to the change rates sum of a substance concentrations in each local process. in this case, the local processes are the – subsystems, and the corresponding differential equations predetermine the local rates. as a result, we have a system of differential equations, which is listed below. dudt=kflows−u+c⋅kdis,nrfa2b2unrfa2b2u−kdis,nrfa2b2u−1nrfa2b2×u++kdis,nirc5nirc5−kdis,nirc5u−1nirc5×u+kcat,nircoutnirc5w,dwdt=kcat,nircinnirc5u+kdis,nirc5wnirc5w−kdis,nirc5w−1nirc5×w++kdis,nirb2dwnirb2dw−kdis,nirb2dw−1nirb2d×w−kflow⋅w,dnircdt=ksnirmniru+5kdis,nirc5nirc5−kdis,nirc5−4nirc5−kd,nirc+kflownirc,dnirc5dt=kcat,nircinnirc5u+kcat,nircoutnirc5w−kdis,nirc5nirc5−kdis,nirc5−4nirc5−−kdis,nirc5unirc5u−kdis,nirc5u−1nirc5×u−kdis,nirc5wnirc5w−kdis,nirc5w−1nirc5×w−kd,nirc5+kflownirc dnirc5udt=−kdis,nirc5wnirc5u−kdis,nirc5u−1nirc5×u−kcat,nircin+kflownirc5u,dnirc5wdt=−kdis,nirc5wnirc5w−kdis,nirc5w−1nirc5×w−kcat,nircout+kflownirc5w,dnirbdt=ksnirmniru+2kdis,nirb2nirb2−kdis,nirb2−1nirb2−kd,nirb+kflownirb,dnirddt=ksnirmniru+kdis,nirb2dnirb2d−kdis,nirb2d−1nirb2×nird−kd,nird+kflownird,dnirb2dt=kdis,nirb2nirb2−kdis,nirb2−1nirb2++kdis,nirb2dnirb2d−kdis,nirb2d−1nirb2×nird−kd,nirb2+kflownirb dnirb2ddt=−kdis,nirb2dnirb2d−kdis,nirb2d−1nirb2×nird+kcat,nirbnirb2dw++kdis,nirb2dwnirb2dw−kdis,nirb2dw−1nirb2d×w−kd,nirb2d+kflownirb2d,dnirb2dwdt=−kdis,nirb2dwnirb2dw−kdis,nirb2dw−1nirb2d×w−kcat,nirb+kflownirb2dw,dnrfacdt=ksnrfmnrfu−ktnrf,cp,uusnrfac−ktnrf,pcnrfa−kd,nrfac+kflownrfac,dnrfbcdt=ksnrfmnrfu−ktnrf,cp,uusnrfbc−ktin,pcnrfb−kd,nrfbc+kflownrfbc,dnrfadt=δperipl⋅ktnrf,cp,uusnrfac−ktnrf,pcnrfa+kform,nrfabkdis,nrfabnrfab−nrfa⋅nrfb−kd,nrfa+kflownrfa,dnrfbdt=δperipl⋅ktnrf,cp,uusnrfbc−ktnrf,pcnrfb+kform,nrfabkdis,nrfabnrfab−nrfa⋅nrfb−kd,nrfb+kflownrfb,dnrfabdt=2kdis,nrfa2b2nrfa2b2−kdis,nrfa2b2−1nrfab2−kdis,nrfabnrfab−kdis,nrfab−1nrfa⋅nrfb−kd,nrfab+kflownrfab,dnrfa2b2dt=kdis,nrfa2b2unrfa2b2u−kdis,nrfa2b2u−1nrfa2b2×u+kcat,nrfa2b2unrfa2b2u.−−kdis,nrfa2b2nrfa2b2−kdis,nrfa2b2−1nrfab2−kd,nrfa2b2+kflownrfa2b dnrfa2b2udt=−kdis,nrfa2b2unrfa2b2u−kdis,nrfa2b2u−1nrfa2b2×u−kcat,nrfa2b2u+kflownrfa2b2u. the first equation in the model describes the total change rate of the nitrite concentration in the chemostat, which is a sum of nitrite inflow and outflow rates into/from the chemostat, the rate of the nrfa, nirc and nitrite enzyme-substrate complexes formation and the rate of the nitrite export from e. coli cells. since c is the total cell volume fraction in the chemostat volume, in the first equation rates of processes that occur within an individual cell are multiplied by the c value. other equations describe processes that occur within an individual cell. the second equation describes the change rate of the intracellular nitrite concentration, and other equations–formation processes of the nrfa and nirb reductases and nirc transporter active forms, transport and catalytic reactions. estimation of the model parameters the complete list of model parameters and their values used in the calculations are represented in the appendix: table the procedure for obtaining parameter values for the mnrf, mnir and u functions was described above and ). in this section we will consider the model parameters, the values of which were obtained on the basis of the existing experimental data on the dynamics of the e. coli cell culture growth under glucose-limiting conditions in the flow chemostat, the kinetics of the nitrite utilization by nrfa and nrfb nitrite reductases and parameters of the nirc protein mediated nitrite transport. parameters of the protein passage from the cytoplasm to the periplasm under the membrane potential influence were estimated on the basis of the assumption that nrf enzyme is transported from the cytoplasm to the periplasm as nrfa and nrfb monomer subunits, which are synthesized in the cytoplasm ). the ratio of the cytoplasmic volume to the periplasmic volume δperipl was assessed based on the experimental data, according to which the periplasmic space proportion makes 8–40 % of the whole e. coli cell volume. on this basis, it is determined that the δperipl value may vary in the – range. in the model δperipl is equal to six. when evaluating the degradation rates of nrfА, nrfВ, nirb, nird, nirc proteins and their complexes, we considered the mean half-life of a pool of e. coli cytoplasmic proteins and stability of the periplasmic fraction of protein complexes, which is significantly higher. the rate constant values for the nrf reductase dimeric and trimeric forms dissociation were estimated based on the experimental data. the rate constant values for the nrfab, nrfa2b nirc nirb and nirb2d complexes formation were evaluated on the basis of indirect data on the protein-protein interaction kinetics. the e. coli cell culture volume fraction in the chemostat volume was assessed based on the previous data, according to which when glucose concentration in the medium is fixated and glucose inflow rate is constant, the cell growth rate is given and constant. the cell cycle duration under such conditions is 70 min. accordingly, the c parameter in the model is assumed constant. considering the described above parameters for the culture development in the chemostat and the rate of glucose uptake by the cells under similar limiting conditions, the range of the c parameter possible variation was evaluated: ≤c≤. the c =  parameter value was taken in the model analysis. values of other parameters that have not been evaluated on the basis of experimental data were estimated in the process of model numerical adaptation to the experimental data on the nitrite accumulation in the chemostat. the model analysis approach the model analysis approach was selected based on the purpose of the study, which was to explain the kinetic data on the nitrite accumulation in the chemostat according to molecular-genetic mechanisms of nitrite utilization in e. coli cell cultured in a flow chemostat. since the target experiment performed in consisted of a series of measurements of nitrite equilibrium concentrations, which were established in the chemostat at every given constant nitrite inflow rate, calculations of the positive stationary points were also performed in the model. positive stationary point for each given s value of the added nitrate was found numerically by solving the system of algebraic equations, which was obtained by equating the right-hand side of the model to zero. the mathematica program was used for all calculations. positive solutions found in all calculations were the only ones in a positive change range of the system arguments.  assessment of the contribution of mechanisms regulating the activity of enzymes involved in no metabolism and transport to the control of no intracellular levels as discussed above, the existing amount of experimental data on nrfa, nirb and nirС gene expression regulation and nitrite accumulation dynamics during e. coli stationary phase growth in the chemostat, presented by wang et al, allowed us to develop a mathematical model of nitrite intracellular utilization and investigate the role of various components of the nitrite utilization system in controlling nitrite intracellular level. after adapting the model to the experimental data, the model calculation results, presented in fig.  demonstrated a good agreement between theoretical curve and quantitative data on the nitrite accumulation dynamics in a chemostat. the maximum value of the intracellular nitrite concentration in the studied range of the added nitrite did not exceed 14 μm, which is in qualitative agreement with the experimental data on the maintenance of sufficiently lower, less than  mm, nitrite concentration in the cell under e. coli cell cultures anaerobic growth at 20 mm nitrate, when the nitrite concentration in the growth medium reaches 5 mm.fig. the steady-state nitrite concentration in a chemostat and inside a cell. curve – model calculation, dots – experimental data from wang et al., x-axis – added nitrite concentration, y-axis: а – steady-state nitrite concentration in the chemostat, b – nitrite concentration in a cell. parameter values are presented in the appendix: table 1 using the model calculations, various components of the nitrite utilization system in e. coli cells cultured in a flow chemostat at different concentrations of the added nitrite were evaluated for their contribution to the nitrite accumulation dynamics in the chemostat and inside the cell.fig. contribution of the no utilization and transport system various components to the nitrite accumulation dynamics. dots – experimental data on the nitrite utilization by the cell culture population, calculated with the kflow/c law based on the data from; curves 1- – model calculation; curves 6- – model calculation; curve – the rate of nitrite utilization by the perilasmic nrfa reductase; curve – the rate of nitrite import by the nirc transporter; curve – the rate of nitrite export by the nirc transporter, measured in negative units; curve – the rate of nitrite utilization by the cytoplasmic nirb reductase; curve – the total rate of nitrite utilization by nrfa and nirb reductases; curve – the rate of nitrite utilization by the nrfА reductase, implemented in the model, excluding the effect of the membrane potential on the nira and nirb proteins diffusion rate kt nrf,cp =   s−1; curve – difference between curves and 6; curve – the total rate of nitrite assimilation by the cell, calculated with. other parameters values are presented in the appendix: table 1 we see the periplasmic nrfa reductase and cytoplasmic nirb reductase contributions to the nitrite utilization dynamics in the chemostat, the experimental parameters of which were calculated based on the data from and are shown with dots. the nirc transporter activity dynamics in relation to the nitrite import into the cell and its export from the cell was demonstrated. analysis of results has shown that the range of the added nitrite concentrations can be divided into two nominal range areas: the area of low values 0 ≤ s ≤ ~  mm and the area of high values s > ~  mm, in which nitrite is utilized only by the periplasmic or the cytoplasmic nitrite reductase, respectively. this is due to the functional peculiarities of the nrf and nir operons, encoding these enzymes, the activity of which depends on the nitrite concentration. let us consider the area of high values of the added to the chemostat nitrite. in this range there is a high nitrite processing efficiency of the nirb reductase, which is sufficient for describing the experimentally observed nitrite utilization rate in e. coli cells, and low nitrate processing efficiency of the nrfa reductase, which is consistent with the expression dynamics of encoding these enzymes operons. moreover, according to the model, the nitrite import rate strictly corresponds to the rate of its processing by the nirb reductase, indicating the close association of the intracellular nitrite utilization system with the transport system of nitrite import into the cell at high substrate levels, and is a result of the arrangement of genes encoding their structure in one operon. however, as seen in fig.  with the parameter values taken, nitrite export from the cell into the medium by the nirc transporter is practically absent. such property of the model, in our opinion, is a result of some parameter values taken during calculations and is associated with the transport system, the evaluation of which was carried out based on the s. typhimurium nirc catalytic property data. according to the, s. typhimurium nirc is active in the millimolar range of the nitrite concentration. these data transfer to the e. coli cell system puts practically no restrictions on the nitrite import process as steady-state nitrite concentration in the chemostat in the range of added nitrite >   mm is higher than  mm. however, when using those parameters of nitrite utilization in the cell by the cytoplasmic nirb reductase, taken from, in the model calculations, when s >   mm, the intracellular nitrite concentration was in the range of to 15 μm. it is not surprising that at such nitrite concentrations in the cell the export activity of the enzyme is practically absent. it should be noted that the parameters of the nitrite export from the cell indirectly indicate the concentration threshold of nitrite toxicity for cells. therefore, if the value kmnirc,out = 1 mm is valid for e. coli cells, the intracellular nitrite becomes toxic for e. coli cells only at concentrations of about 1 mm. based on this, under experimental conditions, implemented in, the cell functions far from the toxicity threshold and has no need for nitrite export from the cell. however, experimental data indicate the significant nirc transporter export activity, so the need for kmnirc,out =1мМ value to explain the data on the nitrite accumulation in the chemostat in the range of high added nitrite concentration remains an open question. we will consider this issue below. let us now consider the range of low added nitrite concentrations. as stated above, in this range, the inflowing nitrite is utilized only by the periplasmic nrfa reductase. high nrf reductase activity is ensured by the combined effect, which is achieved through the activation of nrf operon expression, the de novo synthesis of nrfa and nrfb proteins, and through the increased rate of nrfa and nrfb proteins passage from the cytoplasm to the periplasm as a result of the membrane potential u formation. let us consider the potential contribution of each mechanism in the formation of the nitrite utilization rate in e. coli cells cultured in the chemostat. contribution evaluation of the genetic component of the nrfa periplasmic reductase activity regulation mechanism to the no utilization in e. coli cells to assess the contribution of the nrf reductase encoding nrf operon expression regulating mechanisms to the nitrite utilization in the chemostat at micromolar concentrations of the substrate we explored the model functioning dynamics, ignoring mechanisms of membrane potential formation. to do this, we put value du =  excluding subsystem from the model. in this version of the model, to achieve best agreement between calculations and experimental data, we had to increase the rate constant value of the constitutive nrfa and nrfb proteins transport from the cytoplasm to the periplasm, compared to the baseline, by times. the results are shown in fig. . it is well seen that the nrfa reductase activity obtained by calculating the model, in which the nrfa activity is determined only by the nitrite-dependent genetic mechanisms of the nrf operon expression regulation at 1 mm point, is significantly lower than that obtained by calculating the total model. the difference between these activities is a bell curve, which corresponds to the added activity, which provides an adequate description of the experimental data on nitrite utilization in the chemostat by the model as shown in fig.  with dots. this added activity for describing the same data is missing in the model. thus, an adequate description of the nitrite accumulation kinetics in the chemostat in the micromolar range of the added nitrite concentration cannot be achieved on the basis of the genetic data of nitrite-dependent nrf operon expression activity. from here, we conclude that membrane potential plays an important role in the nitrite utilization dynamics in the micromolar range of the substrate concentration in the chemostat. above, we have given calculations that showed that it is sufficient to consider the membrane potential impact on the rate of the enzyme redistribution between the cytoplasm and the periplasm for explaining the experimental data. however, since there is evidence of other possibilities of membrane potential influence on the nrfa reductase activity, it is necessary to consider them as alternative mechanisms for explaining the data on nitrite consumption in the chemostat. these mechanisms are discussed in the next section. mechanisms of influence of the membrane potential on the periplasmic nrfa reductase activity we should note that, according to the existing data, membrane potential plays a significant role in the formation of catalytically active nrf enzyme molecule, affecting its passage into the periplasm, its stability and catalytic properties. it is known that protein passage into the periplasm is initiated by the electric membrane potential advent as a result of the activity of respiratory enzymes involved in the proton gradient formation. nrf reductase refers specifically to such enzymes and nitrite is an electron acceptor. proteins passage into the periplasm affects their stability. enzyme assembling and correct orientation in the periplasmic space is also dependent on the presence of membrane potential. moreover, it has been shown that michaelis constant for the nrf enzyme is dependent on the membrane potential value, and the membrane potential value is dependent on the nitrite concentration. to answer the question of which of the above mechanisms is essential to explain the experimental data, we evaluated the each mechanism contribution to the nitrite utilization at low added nitrite levels in the chemostat. first, we tested if the enzyme catalytic properties change under the influence of the membrane potential can affect the nitrite utilization kinetics in the chemostat. it is known that michaelis constant for the nrfa enzyme varies depending on the membrane potential. analysis of the model have shown that the km,nrf value change by an order of magnitude or more, that goes beyond the experimentally observed differences, leads to minor changes of the nrfa mediated nitrite utilization rate in the examined range and cannot be a source of additional enzyme activity at low nitrite concentrations, sufficient to explain the experiments. introduction to the model a positive non-linear relationship between the nrfa reductase catalytic activity and nitrite concentration in the medium allows to achieve the desired effect, however, this hypothesis is not supported by the data on a simple reaction kinetics described by the michaelis-menten equation. regarding the membrane potential impact on the rate of the enzyme passage into the periplasm and catalytically active enzyme molecule assembly, the consideration of this process is consistent with experimental data and, according to the model analysis presented above, is sufficient to describe the experimental data on nitrite accumulation dynamics in the chemostat. it should be noted that passage of enzyme subunits into the periplasm under the membrane potential influence is also accompanied by an increase in their local concentrations and stability. how important are these processes to explain the experimental data? the local concentration coefficient of variation is equal to the ratio of cytoplasmic volume to periplasmic volume, which, based on the experimental data, is equal to six in the model. in essence, the mechanism of membrane potential is to increase the rate of the nrfa and nrfb proteins transport from the cytoplasm to the periplasm and provide the accumulation of proteins synthesized in the cytoplasm, in a smaller volume of periplasm. the δperipl = ratio used in a model is sufficient to ensure that the concentration of nrf enzyme in the periplasm reached a level ensuring adequate nitrite consumption rate in the chemostat, meeting the experimental data. note that the enzyme concentration increase is equivalent to the enzyme catalytic activity value increase, as during the catalytic process the absolute reaction rate is determined by multiplying the catalytic constants by the concentration of the enzyme active molecules. theoretically, it is also possible that the increase of protein stability in the periplasm by almost an order of magnitude, which is also considered in the model, may be a crucial factor. however, according to the calculations, this factor is of secondary importance, and it is not required to explain the experimental data. thus, the model analysis demonstrates that membrane potential significantly contributes to the regulation of the periplasmic nrfa reductase activity and it must be considered in order to explain the nitrite utilization dynamics at substrate concentrations <1 mm. the most probable mechanism is the rate change of nrfa and nrfb proteins passage from the cytoplasm to the periplasm, depending on the membrane potential value. in light of these theoretical results, attention is drawn to the fact that both mechanisms: the nitrite-dependent membrane potential change and the nitrite-dependent effectiveness of nrf operon expression, have qualitatively identical unimodal activity curves. as a result, the necessity of addressing both mechanisms to explain the experimentally observed nitrite accumulation curve in the low range of the added nitrite concentration emerges. this question is discussed in the next section. on the functional redundancy of the no utilization mechanisms in e. coli cell in this section, we analyze the possibility of model functional redundancy in the utilization mechanisms and related issues. the essence of the matter derives from the fact that subsystem –, described above and including a certain number of genetic, transport, complex-forming, and enzymatic reactions, make up the structural basis of the model. we used two approaches to describe the subsystems. regulation of the nrf and nir operons expression efficiency, as well as the membrane potential, we described with phenomenological functions of the class of generalized hill functions. other processes were described with non-equilibrium biochemical reactions. meaning the right-hand side of system includes members that describe both direct and reverse non-equilibrium processes. at the same time, to describe the dynamic characteristics of the nitrite accumulation in the chemostat, the system equilibrium state has to be calculated. it is well known from the general theory of dynamical systems that the equilibrium state specific the value is not determined by the values of the parameters, but by values of their specific combinations, the number of which is usually less than the number of the parameters themselves. the simplest example is the bimolecular reaction. to describe the reaction far from equilibrium, it is necessary to know the rate constant values of the forward and reverse reactions. however, at equilibrium, it is sufficient to know the value of one equilibrium constant, which is the ratio of non-equilibrium constants. as a result, we get a traditional problem of determining the degree of functional and parametric redundancy of mechanisms reflected in, which are necessary to describe the nitrite consumption kinetics in the chemostat. in this context, we consider the redundancy issue of mechanisms, by which the nitrite accumulation kinetics in the chemostat at concentrations of the added nitrite s < 2 mm is explained in the model. as stated above, the apparent similarity of the curves describing the nrf operon expression and membrane potential influence on the rate of the nrfa and nrfb proteins redistribution in relation to the added nitrite concentration suggest that in order to explain the experimental data from it is sufficient to use only one of these mechanisms and the other can be omitted. is it so? as for the need to consider the membrane potential, we showed above that without considering the influence of membrane-dependent potential on the rate of the nrfa and nrfb proteins redistribution between the cytoplasm and the periplasm it is impossible to adequately describe the experimental data on nitrite utilization in the chemostat. this output, as shown above, is considerably based on the data on nitrite-dependent genetic regulation of nrf operon expression. however, let’s assume that nrf operon is constitutively expressed. it appears that for this version of model it is not difficult to pick up a set of parameter values  ≡  mnir ≡  ks,nrfa = ks,nrfb =   mm/s, knrf1 = , hnrf,1 =  δnrf,2 = , knr,f2 = , hnrf,2 =  ω2 = , knrf,3 = , hnrf,3 =  values of other parameters are unchanged and presented in the appendix: table 1), for which the model will adequately describe the nitrite accumulation curve. that is, the data on the genetic regulation of the nrf operon expression are redundant for explaining the nitrite accumulation curve in relation to the added nitrite concentration. it is sufficient to assume the nrf operon constitutive expression and take into account the membrane potential impact on the nira and nirb proteins redistribution between the cytoplasm and the periplasm. however, this conclusion is theoretical, since the native e. coli cells reveal a strong nitrite-dependent regulation of the nrf operon expression. and, nevertheless, the revealed redundancy raises important theoretical question of evolutionary expediency of such a complex mechanism of nrf operon expression regulation in e. coli cells: because the cell could easily utilize nitrite in simpler and more economical way! this issue requires further study and more detailed analysis. the second issue considered in this section, is the redundancy of the nitrite export mechanism. it is indeed known that e. coli cells have a certain threshold value of the intracellular nitrite concentration, above which it becomes toxic to the bacteria. it is believed that in this case the cell has a protective system based on the nirc nitrite export activity, which is basic during nitrite breathing. however, the parameters of such activity were not measured. in the model, we used values of the parameters that have been measured for the s. typhimurium nirc transporter. but then, it implies that the nitrite concentration toxic for e. coli cells is comparable to millimolar concentrations, and in the model nitrite export from cells is practically absent. however, there is reason to believe that in e. coli cells nitrite toxic concentration is in the micromolar range, i.e., michaelis constant for the nirc mediated nitrite export can be much smaller than the one given in appendix: table on this basis, we have considered a model version, in which the value of the michaelis constant for the nirc mediated nitrite export from the cell is times smaller: km,nirc,out = μМ. the corresponding calculation is presented in fig.  fig. nitrite utilization and transport rate in relation to the concentration of the added nitrite. dots – the rate of nitrite consumption by the cell, calculated based on the experimental data from; curve – the rate of nitrite utilization by the periplasmic nrfa reductase; curve – the rate of nitrite import into the cell by the nirc transporter; curve – the rate of nitrite export from the cell by the nirc transporter, measured in negative units; curve – the rate of nitrite utilization by the cytoplasmic nirb reductase; curve – the total rate of nitrite utilization by nrfa and nirb reductases. calculations performed with the model version, in which k М,nirc,out =  mМ, k s,nirc,0 =   mm/s; other parameter values are from the appendix: table 1 calculations have shown that such parametric version of the model does not impair the data reproduction quality. at the same time, under a given set of parameter values there is a high rate of nitrite export from the cell into the chemostat. consequently, the ratio of nirc transporter import activity functional parameter to its export activity functional parameter is a free parameter in the model, and to determine its value, additional information that is not available in the scientific literature is needed. in conclusion, the nitrite concentration curve at different values of the added nitrite concentration, shown in fig. 4b, in a sense, is nominal. it is easy to choose sets of parameters for which the model calculations of the intracellular nitrite concentration will be different from the above calculation in orders of magnitude, and, at the same time, the accuracy of the nitrite accumulation approximation curve will not be altered, so this characteristic is also free in the model and for its more accurate characterization it is necessary to have additional data, which is not yet available.  in the present work we introduce the model of nitrite utilization in e. coli cells cultured in the flow chemostat. the parameter values were chosen in such a way that the model adequately describes the experimental curve of nitrite accumulation in the chemostat. the model analysis revealed the existence of two ranges of the added no concentration: low and high. according to the model, in the high concentration range the nitrite is substantially utilized by the nitrite transport/utilization system. in this range, the nitrite-dependent genetic regulation of the nir operon expression is sufficient for explaining experimental data on the nitrite accumulation in the chemostat. on the contrary, in the low concentration range of the added nitrite the known genetic mechanisms of the nrf operon expression are not enough for explaining experimental data on the nitrite accumulation in the chemostat. analysis of different hypothesis has shown that the most probable additional mechanism of the periplasmic nrf reductase activity regulation in the micromolar concentration range of the added nitrite, that allows co-ordination of physiological and genetic data, is a local enzyme concentration change due to its passage from the cytoplasm to the periplasm under the influence of nitrite concentration-dependent membrane potential. in the framework of our model, the potential-dependent mechanism of the nrfa and nrfb proteins redistribution between cytoplasm and periplasm is the essential element required to explain observed nitrite accumulation dynamics in the chemostat in the low concentration range of the substrate. analysis of different nitrite utilization and transport mechanisms contribution to the nitrite accumulation in the chemostat and inside the cell revealed some redundancy of molecular-genetic mechanisms of nrfa, nrfb reductases and nirc transporter activities regulation. thus, to describe the nitrite accumulation curve in chemostat in the entire concentration range of the added nitrite it is sufficient to assume a constitutive nrf operon expression, whereas the nirc transporter and nirb nitrite reductase can be ignored. in the high concentration range of the added nitrite the ratio of nirc transporter import activity functional parameter to its export activity functional parameter as well as intracellular nitrite concentration are free parameters of the model and may vary over a wide range without affecting the nitrite accumulation curve description accuracy of the model. appendix k flow δnrf,1 k nrf,1 h nrf,1 δnrf,2 k nrf,2 h nrf,2 ks nrf ks nirc ksnirb==ksnird ksnirb==ksnirc/5 k d,nrf kt nrf,cp kt nrf,pc k dis,nrfab k dis,nrfab kdis,nrfa2b2 kdis,nrfa2b2 δ peripl d u δ k pmf,1 h pmf,1 k pmf,2 h pmf,2 ω ω 2 с k pmf,3 с h pmf,3 с k d,nrfa sec –1 k d,nrfa k d,nrfab kd,nrfa2b2 kcat,nrfa2b2u kdis,nrfa2b2u δnir,1 k nir,1 h nir,1 δnir,2 k nir,2 h nir,2 kdis,nirb2 kdis,nirb2d kcat,nirb2dw kdis,nirb2dw kdis,nirc5 k cat,nircin kdis,nirc5u k cat,nircout kdis,nirc5w k d,nirc k d,nirc5 k d,nirb k d,nird k d,nirb2 k d,nirb2d c a – mm, sec, if dimension is not specified, the quantity is dimensionless b – parameter value was estimated based on the experimental data c – members, before which the ω and ω coefficients are standing in the membrane potential formula, are zeroed in model, the formula used in the model calculation is with nrf operon constitutive expression and enabled nirc and nirb enzymes functions competing interests the authors declare that they have no competing interests. authors’ contributions tmk conceived the study. val and nar developed the model and analyzed the data. tmk and nar evaluated the model parameters. val and nar performed the model calculations. tmk wrote the manuscript with editing by val and nar. all authors read and approved the final manuscript. 
nrfa reductase is localized in the cell periplasm and uses no as an electron acceptor to create a proton gradient; nirb reductase is restricted to the cytoplasm and metabolizes excessive nitrite inside the cell, the uptake of which is mediated by the transporter protein nirc. it is not sufficient to consider molecular-genetic mechanisms of nrfa reductase activity regulation to describe the nitrite accumulation dynamics in the chemostat in the micromolar range of added nitrite concentrations. using a mathematical model describing the nitrite utilization in e. coli cells cultured in a flow chemostat, the role of enzymes involved in nitrite metabolism and transport in controlling nitrite intracellular levels was investigated. at nitrite concentrations in the media more than 2 mm, the model adapted to the experimental data on nitrite utilization dynamics in e. coli cells cultured in the flow chemostat demonstrates the largest contribution of genetic mechanisms involved in nrf and nir operons activity regulation to the control of nitrite intracellular levels. due to a high toxicity of nitrite and its metabolites, it is of high interest to study mechanisms underlying the low no level maintenance in the cell. analysis of different hypotheses has shown that the mechanism of local enzyme concentration change due to membrane potential-induced diffusion from the cytoplasm to the periplasm at low nitrite levels is sufficient to explain the nitrite accumulation dynamics in the chemostat. escherichia colianaerobic respirationmathematical modelgene expression regulationnitrite reductasenitrite transportmembrane potentialthe 7th international young scientists school “systems biology and bioinformatics” novosibirsk, russia 22- june http://conf.nsc.ru/plantgen2015issue-copyright-statement© the author 2016 the model predicts a significant contribution of the membrane potential to the periplasmic nrfa nitrite reductase activity regulation and nitrite utilization dynamics at substrate concentrations ≤1 mm. it was demonstrated that the model adapted to the experimental data on expression of nrfa and nirb genes encoding nrfa and nirb nitrite reductases, can describe nitrite accumulation kinetics in the chemostat in the millimolar range of added substrate concentrations without any additional assumptions.
 due to a high toxicity of nitrite and its metabolites, it is of high interest to study mechanisms underlying the low no level maintenance in the cell. during anaerobic growth of escherichia coli the main nitrite-reducing enzymes are nrfa and nirb nitrite reductases. nrfa reductase is localized in the cell periplasm and uses no as an electron acceptor to create a proton gradient; nirb reductase is restricted to the cytoplasm and metabolizes excessive nitrite inside the cell, the uptake of which is mediated by the transporter protein nirc. while it is known that these three systems, periplasmic, cytoplasmic and transport, determine nitrite uptake and assimilation in the cell as well as its excretion, little is known about their co-ordination.  using a mathematical model describing the nitrite utilization in e. coli cells cultured in a flow chemostat, the role of enzymes involved in nitrite metabolism and transport in controlling nitrite intracellular levels was investigated. it was demonstrated that the model adapted to the experimental data on expression of nrfa and nirb genes encoding nrfa and nirb nitrite reductases, can describe nitrite accumulation kinetics in the chemostat in the millimolar range of added substrate concentrations without any additional assumptions. according to the model, in this range, low intracellular nitrite level, weakly dependent on its concentration in the growth media, is maintained. it is not sufficient to consider molecular-genetic mechanisms of nrfa reductase activity regulation to describe the nitrite accumulation dynamics in the chemostat in the micromolar range of added nitrite concentrations. analysis of different hypotheses has shown that the mechanism of local enzyme concentration change due to membrane potential-induced diffusion from the cytoplasm to the periplasm at low nitrite levels is sufficient to explain the nitrite accumulation dynamics in the chemostat.  at nitrite concentrations in the media more than 2 mm, the model adapted to the experimental data on nitrite utilization dynamics in e. coli cells cultured in the flow chemostat demonstrates the largest contribution of genetic mechanisms involved in nrf and nir operons activity regulation to the control of nitrite intracellular levels. the model predicts a significant contribution of the membrane potential to the periplasmic nrfa nitrite reductase activity regulation and nitrite utilization dynamics at substrate concentrations ≤1 mm. keywords escherichia colianaerobic respirationmathematical modelgene expression regulationnitrite reductasenitrite transportmembrane potentialthe 7th international young scientists school “systems biology and bioinformatics” novosibirsk, russia 22- june http://conf.nsc.ru/plantgen2015issue-copyright-statement© the author 2016 