Team:Utah State/Modeling

There is currently limited information on how to properly select which secretion pathway should be used for a given protein (Choi, 2004). The construction of models that evaluate protein secretion may provide a useful framework for studying and attempting to predict secretion efficiency. Models of the Sec and Tat translocation pathways were made using MATLAB’s Simbiology toolbox. They were evaluated using the embedded ordinary differential equation solver in the software. Assumptions were made for both models due to lack of time scale information for individual steps for these translocation mechanisms. However, both provide a flexible framework that can become more detailed as additional information is found in literature or discovered in the laboratory. Within both models, a protein is first generated and then carried out of the cytoplasm to the periplasm. Only protein species were tracked, as species involved with making the protein are currently not as useful to monitor.

In both models, parameters for the manufacture of protein were the same and found using averaged values for E. coli. The average length of an mRNA for E. coli is 1100 nucleotides and an E. coli cell transcribes at an average rate of 70 nucleotides per second. A simple rate calculation determines that a strand of mRNA is made approximately every 15.7 seconds. This value was then used to yield a first order reaction rate of (1/15.7) 1/second. The translation rate was found in a similar manner. Given an average protein size in E. coli as 360 amino acids and the average translation rate is 40 amino acids per second, the first order translation rate constant is (1/9) 1/s (Institute for Biomolecular Design, 2008). To initiate the model, the concentration of the gene was set at 1, which was locked at that value over the course of the simulation. The median half life for mRNA was found to be 3.7 minutes, after which it will degrade (Milo, n. d.). For the Tat-dependent mechanism, the protein is fully folded prior to translocation (Mergulhão, 2005). For this folding time, a place holder value of 10 minutes was used as an approximate order of magnitude, although the value can easily be changed depending on parameters of the protein in question. Degradation of protein was modeled as the average time it takes the cell to reproduce, which is given as half an hour (Institute for Biomolecular Design, 2008).

In addition to the discussed values regarding the manufacture of protein, parameters for the Sec and Tat pathway are required. A literature review yielded no specific rate constants or time parameters for the individual steps involved with either process. However, one article stated that the Sec and Tat pathways take a few seconds and a few minutes, respectively, to translocate protein to the periplasmic space (Mergulhão, 2005). Accordingly, the Sec and Tat pathways were modeled to require 3 seconds and 3 minutes for protein translocation, respectively. First order rate constants were determined by taking the inverse of these values in seconds. Images of the model diagrams are seen below. The different pathways after the formation of protein can be activated and deactivated as more detailed information comes available about the translocation and secretion processes. The Sec secretion pathway also includes the SRP secretion pathway as a possible option. The constructed models for the Sec and Tat pathways are given as Figures 1 and 2, respectively.

Simulations for the Sec pathways were run for 5400 seconds (1.5 hours), which allowed the number of periplasmic to reach approximately steady state. The simulation for the Sec pathway secretion model is shown below as Figure 3.

Protein generated by the E. coli cell reaches steady state quickly, as it is both produced and translocated in to the periplasm rapidly. The periplasmic protein reaches approximately 90 proteins upon achieving steady state. The long settling time and high number of periplasmic protein are a result of the long protein degradation time.

Simulations for the Tat pathway were run for 3600 seconds (1 hour). The shorter length of simulation relative to the Sec model is used despite having significantly longer protein translocation time because the time-consuming folding step limits the reaction. The results of the ordinary differential equation simulation are shown below in Figure 4. The amount of periplasmic protein is lower than cytoplasmic and folded protein values upon reaching steady state. This is due to the rate-limiting folding step, which also leads to an expected high concentration of folded protein as well.

The accuracy of models, such as those constructed, will increase as more knowledge about the time step parameters for each of the mechanisms is determined. These improvements will allow for a better comparison of the two secretion methods, which can lead to better signal peptide selection (and corresponding pathway selection) for the protein in question. Future work could also include incorporating promoter strength into the model to affect transcription rates to provide a more accurate depiction of what is happening inside the cell. Better correlation between these parameters and the model will result in a more useful aid in when studying protein secretion.

• Institute for Biomolecular Design (2008, August 1) CyberCell Database. Retrieved October 2009.
• Choi JH, Lee SY (2004) Secretory and extracellular production of recombinant proteins using Escherichia coli. Appl Microbiol Biotechnol 64:625-635
• Mergulhão FJM, Summers DK, Montier GA (2005) Recombinant protein secretion in Escherichia coli. Biotechnol Adv 23:177-202
• Milo R (n.d.) BioNumbers: The Database of Useful Biological Numbers. Retrieved October 2009.