Team:Heidelberg/Templates/Modelling/Ind-Production

\begin{align} \mathrm{d}\mathrm{[Bac]}/\mathrm{d}t &= \frac{\mathrm{[Bac]} \cdot \left(\mathrm{[Bac]} - \mathrm{Bacmax}\right) \cdot \mathrm{beta}}{\mathrm{Bacmax}}\label{bacgrowth}\\ \mathrm{d}\mathrm{[Glu]}/\mathrm{d}t &= - \mathrm{[Bac]} \cdot \mathrm{[Glu]} \cdot \mathrm{ksyn}\label{glu}\\ \mathrm{d}\mathrm{[cGlu]}/\mathrm{d}t &= - \mathrm{kdim} \cdot {\mathrm{[cGlu]}}^2 + \mathrm{[Bac]} \cdot \mathrm{[Glu]} \cdot \mathrm{ksyn}\label{cglu}\\ \mathrm{d}\mathrm{[Ind]}/\mathrm{d}t &= {\mathrm{[cGlu]}}^2 \cdot \mathrm{kdim} - \mathrm{[Ind]} \cdot \mathrm{kdegi}\label{ind} \end{align}

Initially, inhibition of bacterial growth by indigoidine and degradation of cyclic glutamine were also described by the model, using equations \eqref{bacgrowthinhib} and \eqref{cgludeg} instead of \eqref{bacgrowth} and \eqref{cglu}, respectively, however the degradation rate turned out to be non-identifiable for all data sets, converging to the lower bound, whereas the growth inhibition term led to frequent convergence failures, the inhibition constant being very low in cases where convergence was achieved. These two terms were thus removed from the model. \begin{align} \mathrm{d}\mathrm{[Bac]}/\mathrm{d}t &= \frac{\mathrm{[Bac]} \cdot \left(\mathrm{[Bac]} - \mathrm{Bacmax}\right) \cdot \left(\mathrm{beta} - \mathrm{ki} \cdot \mathrm{[Ind]}\right)}{\mathrm{Bacmax}}\label{bacgrowthinhib}\\ \mathrm{d}\mathrm{[cGlu]}/\mathrm{d}t &= - \mathrm{kdim} \cdot {\mathrm{[cGlu]}}^2 - \mathrm{kdegg} \cdot \mathrm{[cGlu]} + \mathrm{[Bac]} \cdot \mathrm{[Glu]} \cdot \mathrm{ksyn}\label{cgludeg} \end{align}

• Bacmax: maximum capacity for bacterial growth
• beta: maximum attainable growth rate
• ksyn: synthesis rate of cyclized glutamine from glutamine
• kdim: dimerization rate of two cyclized glutamines to an Indigoidine dimer
• kdegi: degradation rate of indigoidine

Framework

as well as identifiability analysis yielding valuable insight into bacterial indigoidine production.

Results

The data sets could be grouped into two types, those with a significant indigoidine production and those without. Although both types of data sets resulted in seemingly good fits (Fig. 4), identifiability analysis using a profile likelihood exploiting approach showed that for all data sets without indigoidine production at least one parameter could not be reliably estimated (an exemplary plot is shown in Fig. 5B), whereas all parameters were identifiable for data sets showing indigoidine production (Fig. 5A). As our model only permits relative quantification of synthesis rates due to the lack of glutamine production and the arbitrarily chosen initial glutamine concentration, we normalized the indigoidine synthesis rates of the fully identifiable data sets to the range [0, 1] and visualized them in a heatmap (Fig. 3). The synthetic T domain T10 seems to be the most effective T domain for the entD and sfp PPTases, which corresponds very well to the preliminary analysis performed by the indigoidine team. Our mathematical model was able to describe most of our experimental data and was in those cases structurally and practically fully identifiable. We thus challenged the model to deduce additional facts on indigoidine synthesis. First, we had a look at the parameters that were estimated to describe our experimental data sets. For each parameter estimate, the lower boundary was set to 1e-5 and the upper boundary was set to 1e+5. Initial guesses for parameters during parameter estimations were generated via latin hypercube sampling to ensure broad coverage of all regions in the high-dimensional parameter space (http://www.jstor.org/stable/2346140). Nevertheless, some parameters exhibited broad variation among the best fits of the different data sets while other parameters were estimated in a narrow range (Fig. XXX). Interestingly, the bacterial growth rate beta and the maximum growth capacity Bacmax were very consistent. The parameters ksyn, kdim and kdegi contributed directly to indogoidine synthesis as it can be inferred from the above mentioned equations (2), (3) and (4). The synthesis rate ksyn characterizes the cyclization of glutamine. The dimerization rate kdim reflects how two cyclized glutamine molecules form one indgoidine molecule. As such, ksyn and kdim represent the efficiency of the T domain of indC. The degradation rate of indigoidine, kdegi, should be independent from the domain structure of the indigoidine synthetase as the stability of the peptide is not affected by the T domains or the PPTases that were tested. The rate did not vary widely and the diverse data sets could still be explained by the mathematical model. At first glance, kdegi seemed less important for our experimentally assessed conditions. To further investigate the role of kdegi for indigoidine production kinetics, we systematically varied this parameter and compared simulation results. The synthesis of indigoidine was clearly affected by variation of kdegi (Fig. XXX+1). The [Ind] trajectory collapsed more than 3-fold increase of the degradation rate. On the contrary, improved peptide stability could yield much more indigoidine. For very low degradation rates (i.e. kdegi < 1e-3), synthesis would not even saturate within the observed time frame.

Sensitivity analysis for optimized indigoidine production

To perceive the role of kdegi for indigoidine synthesis, we have varied the parameter over a broad range and observed qualitatively distinct model trajectories. For a more quantitative understanding of how indigoidine yield changes with altering parameters of the system, we conducted a sensitivity analysis (REVIEW NEEDED). Therefore, we first defined the variable IntInd, the integrated indigoidine reflecting production yield. The change of the indigoidine concentration with time is determined by the ODE (4) and its integral with time can be depicted as the area under the [Ind] curve, i.e. the amount of indigoidine produced by the bacteria with time. Now we were interested in the change of IntInd with the change of the different parameters of our system. Mathematically, this is expressed as the so called sensitivity $\frac{\partial \mathrm{X}}{\partial \mathrm{p}}, with X as our variable of interest (IntInd) and p as the parameter (beta, Bacmax, ksyn, kdim or kdegi, respectively). This partial derivative provides a measure of how much X changes with the change of p or how sensitive X is against small variations in p. Sensitivities can be scaled with the actual values of X and the respective p to fall in the range from -1 to +1. The higher the absolute value of the sensitivity, the more important is the parameter for the variable. A positive sensitivity means positive impact of the parameter on the variable whereas a sensitivity < 0 indicates reduction of X for subtle decrease in p. We utilised our mathematical model and computed the sensitivities of IntInd against the change of all five kinetic parameters of our system and scaled them accordingly. The analysis confirmed that the kdegi was the only parameter of the system that exhibited a negative impact on the indigoidine production (Tab. 1). Slight increase of all the other parameters also augmented additional yield. Sensitivity analysis for native_sfp:

Parameter	$$\left.\mathrm{[Bac]}\right|_{t=0}$$	kdegi	kdim	beta	ksyn	Bacmax
$$\frac{\partial\int\!\mathrm{Ind}\,\mathrm{d}t}{\partial p}$$	0.784	-0.7024	0.5117	0.3338	0.2919	0.1036

The initial concentration of bacteria represented the most crucial parameter for indigoidine production, which was intuitive since the amount of produced indigoidine depends on the amount of bacteria in culture that synthesise the peptide. The parameter kdegi exhibited a high impact on IntInd but did not vary so much among the experimental conditions. One would assume that the half-life of the peptide is not affected by different T domains or PPTases. Possibly, the arrangement of modules influences peptide stability as inter modular communication is a precondition for robust peptide assembly (REFERENCE). This links the peptide synthesis optimisation to our module shuffling project where modules of different hosts were combined for custom peptide synthesis. Indigoidine yield was also sensitive against small variations in beta and Bacmax. Culture conditions could be improved and metabolically more active strains could be chosen to further optimize performance in this respect. The parameters ksyn and kdim also determine peptide production. The question is whether those rates could be altered by engineering the NRPS to advanced performance. To answer this question, we performed a multiple regression analysis (Fig. XXX+2). We regressed the alignment score of the synthetic T domains of the indigoidine synthase with the rates from the respective best fit. The score was calculated (http://www.ebi.ac.uk/Tools/msa/clustalw2/) by pairwise multiple sequence alignment of the amino acid sequences of the synthetic T domains with the native T domain of the indigoidine synthase as a template. The parameter values were obtained from parameter estimations of our mathematical model that has been calibrated with experimental data as described above.

Conclusion and Outlook

constraints and a-priori assumptions indigoidine to learn lessons on optimizing production of synthetic peptides of course dependency of yield on rates does not necessarily scale linearly with peptide length but the approach is established