Difference between revisions of "Team:Nanjing-China/Model"

This year our team created a mathematical model to optimize the arrangement of the nif gene cluster. This model helped we refined our design and provided some new perspectives of our nitrogen-fixation system attranscriptional level.

We developed this model with two goals in mind:
1. We want to achieve the putative best stoichiometric proportion of each nif gene, which is nifB: nifH: nifD: nifK: nifE: nifN: nifX: nifV=1: 3: 4: 4: 1: 1: 1: 1.
2. We want our system as simple as possible, that means minimizing numbers of promoters and each nif gene.

We made the following assumptions:
1. There are two kinds of promoters, both of which can successfully launch the expression of every nitrogen fixation gene involved in our system.
2. One promoter is stronger(called H) while the other is relatively weak(called L). Under promoter H, each gene’s transcription level is double that of under promoter L.
3. The order of genes has little influence on their transcription level.

We conducted Real-time Quantitative PCR to detect the transcription level of nif gene cluster and the experimental data we received became an important reference for our modeling.

Method:

To start with, we put all genes into two groups. One group is under the strong promoter while the other is under the weak one. We constructed two arrays, weak[i] and expected[i].

Then we did some necessary preprocessings. Firstly, we found the smallest data in weak[i] and called it “min”. We normalized all the other data accordingly by doing:

We constructed strong[i]:

Secondly, to guarantee the existence of a solution, we adjusted expected[i]* by examining whether it is no less than the corresponding weak[i]*, if not, we did:

After that, we began the organization. In order to minimize the total numbers of genes, we arranged the strong promoter group first, and considered the weak group later. Because each gene can be considered separately, here we only describe the organization of the i^th gene as an example.
For the i^th gene, we tried adding one copy of it under the strong promoter. If

we actually added it. Until we have added (m_i+1) i^th genes, and got

Then we stopped adding it and recorded that we have added m_i i^th genes under the strong promoter.

For the weak promoter group, we applied a similar method. For the i^th gene, we tried adding one copy of it under the weak promoter. If

we actually added it. Until we have added (n_i+1) i^th genes, and got

Then we stopped adding it and recorded that we have added n_i i^th genes under the weak promoter.

In that way, we were able to determine numbers of the i^th gene under the two promoters with which the deviation was the smallest.

Fig 1. A flow diagram describing the idea of our modeling process

According to this flow diagram, we programmed with Python and got the following results:

Fig 2. The best arrangement of nif genes according to our calculation

With this arrangement, the proportion of nifB: nifH: nifD: nifK: nifE: nifN: nifX: nifV= 15.44: 46.93: 71.88: 62.10: 16.44: 16.04: 16.0: 15.94, which is close enough to the ideal proportion among all the solutions.

Refinement of our model:

We modified the putative best expression level of nifB:nifH:nifD:nifK:nifE:nifN:nifX:nifV to 5:3:4:4:1:1:1:1. We believed in this way, we could better simulate the expression of nitrogenase in our engineered E.coli strains. We made this change because of three reasons.

Firstly, nifB is indispensable for nitrogenase assembly no matter in diazotrophs or engineered E. coli strains. Apart from the minimal nitrogen fixation gene cluster, the genomic DNA of wide type Paenibacillus polymyxa includes analogues of nifM, nifU, nifS and other genes which exist in other nitrogen-fixing microorganisms and are essential for the correct folding of nitrogenase iron protein. However, the E. coli genome doesn’t have such analogues. Nevertheless, it has been reported that the excessive expression of nifB can compensate for the absence of nifU and nifS. That is, if nifB is overexpressed in E. coli, these auxiliaries are not necessary. Therefore, the expression level of nifB should be the highest 5.

Secondly, compared with nitrogen-fixing microorganisms, E. coli also lacks some genes that provide electron transfer function, such as nifF and nifJ. So the intracellular reductive power of E. coli is insufficient to accomplish nitrogen fixation. Thus it is necessary to overexpress nifH(nitrogenase reductase) and the value is set to 3 instead of 5 because our semiconductor, the CdS part, can provide additional electrons.

Thirdly, we set the expression level of nifD and nifK to be 4 because molybdenum iron protein is an ɑ2β2 allotetramer and is the core of nitrogenase.

Based on the new ideal stoichiometric proportion, we adjusted the code and received a more accurate result.

Fig 3 The best arrangement of nif genes version 2.0.

The achieved stoichiometric proportion of nifB: nifH: nifD: nifK: nifE: nifN: nifX: nifV=77.23: 46.93: 71.88: 62.10: 16.44: 16.04: 16.0: 15.94, which is close enough to the ideal 5:3:4:4:1:1:1:1.

This model provided a potential strategy for the improvement of biological activity of nitrogenase expressed in our engineered E. coli strain and offered a great help to our further experiments.

References

1. Wang, X., et al., Using synthetic biology to distinguish and overcome regulatory and functional barriers related to nitrogen fixation. PLoS One,2013. 8(7):p.e68677.
2. Yang, J., et al., Modular electron-transport chains from eukaryotic organelles function to support nitrogenase activity. Proc Natl Acad Sci U S A, 2017. 114(12):p.E2460-E2465.
3. Yang, J., et al., Polyprotein strategy for stoichiometric assembly of nitrogen fixation components for synthetic biology. Proc Natl Acad Sci U S A, 2018. 115(36):p.E8509-E8517.
4. Yang, J.G., et al., Reconstruction and minimal gene requirements for the alternative iron-only nitrogenase in Escherichia coli. Proceedings of the National Academy of Sciences of the United States of America, 2014. 111(35):p.E3718-E3725.