Team:Leeds/Modeling

After attending YSB 1.0, we signed on with Manchester to help develop a modelling standard - this also interlinks with our work with Purdue. We hope that by being more tightly involved in these processes, we will better understand how to approach our project in terms of Characterisation and Modelling. Additional thanks to Newcastle for showing us BioNetGen.

Modelling The Cpx Pathway

The key to making MicroBeagle successful hinges on proper integration with the Cpx pathway. As such, it is essential we develop a god working model, not only to predict how much fluorescence we can expect in different environments, but also to prototype and test potential methods for controlling this; in turn reducing our false positive rate.
We are currently using [http://bionetgen.org/index.php/Main_Page BioNetGenLanguage] to write the model code, as this provides a rule-based coding language that is easier to use than SBML. Rule-based modelling allows mechanistic simulation of the pertinent reactants, reducing the need to know the exact kinetics of reactions, as sub-processes can be absorbed into rate constants. Additionally, adaptation of the model should be easier after a basic framework has been produced, as new processes simply require insertion of new rules into the code. Finally, BNGL offers easy integration with [http://www.ibiblio.org/virtualcell/ Virtual Cell], a graphical simulator of cell systems - ideal for outreach with the public and teaching.

Cpx Pathway Model

Due to the high complexity of the Cpx Pathway, it was necessary to graphically plot the processes taking place within the cell by hand, allowing us to then simplify and compartmentalize aspects which, while having an effect on the overall results, were of no interest to us. The above image shows the 3^rd iteration of this, with colour-coded interactions, products and proteins. The main enzymes of interest were specifically named, while a general product labelled as 'denaturants' was listed to encompass the more than 100 other Cpx-mediated protein syntheses. This would be inaccurate, if not for the use of Rules-Based software, such as RuleBender - for the model to work, one must list the starting concentrations of each molecular species, hence by setting the 'denaturant' promoter concentration to be at least 100 times greater than that of pGFP or pCpxP, the relative interaction mechanics are preserved.
So far, the model is still only at version 3, due to the difficulty in confirming the vast number of reaction constants needed - unlike Physics or Chemistry, biological factors will frequently differ by orders of magnitude, dependent upon the conditions they were measured in etc. This is the main downside to using a rule-based approach, as our good friends Manchester and York point out in their modelling support videos. However, the model currently describes the behaviour of CpxP, CpxA, CpxR and the phosphorylation processes these use, as well as the effect of misfolded protein sub-units (caused by membrane stress) and the expected relative production of GFP versus other species produced via CpxR promotion. This should provide a workable, skeletal basis for characterisation of our first Bio-System as GFP concentration scales directly with fluorescence. As we take measurements, we will then feed this back into the model to fine tune it, and thus develop it further for characterisation of our second system. Beyond this, it is hoped the model source code can then be submitted as a beta-stream to the modelling database. It is otherwise available as a raw text document here for future teams to make use of and improve upon.