Team:Newcastle/Modelling/L-form Switch

L-form switch

Although it has been documented that L-forms can be produced by a specific mutation in the murE gene, we strongly suspected that preventing the expression of the the entire murE operon would have the same effect. However we were not certain, and therefore modelled the biosynthesis of peptidoglycan, the primary component of the cell wall, using BioNetGen. Modelling was based off the KEGG pathway for peptidoglycan synthesis in B. subtilis. The model put the operon under the control of a xylose inducable promoter, and demonstrated that in at low xylose concentrations (and thus no expression of the murE operon no peptidoglycan was produced. In the presence of xylose peptidoglycan synthesis did occur. As we were not worried about the quantity of peptidoglycan that was produced, just whether it was, arbritary unit were used for peptidoglycan production. This was also because quantifying peptidoglycan is non-trivial. These results gave the team confidence that the switch BioBrick we designed would work as desired. Experimental results later validated this.

The biosynthesis of peptidoglycan is quite complicated, with one initial pathway splitting into 3, each fork resulting in peotidoglycan being formed. This can be seen by the network below- don't worry if it looks confusing! Synthesis can be broadly defined in a few steps steps. In the first step fructose 6-phosphate is converted to glucosamine-6-phosphate. In step 2 a acetyl group is donated, creating N-acetyl-glucosamine-6-phosphate. This is subsequently isomerized to N-acetyl-glucosamine-1-phosphate in step 3. In step 4 UTP reacts with N-acetyl-glucosamine-1-phosphate forming UDP-N-acetylglucosamine. This is where the KEGG pathway, and our model, starts. UDP-N-acetylglucosamine is converted to UDP-N-acetylmuramic acid, and this undergoes multiple rounds of amino acid addition. It is these steps, indicated on the pathway below, that are governed by a few of the murE operon genes- murD, murE, murF,murG and mraY. The switching off of these genes prevents the addition to UDP-MurNAc-L-Ala of more amino-acids and thus prevent peptidoglycan biosynthesis.

With a functional murE UDP-MurNAc pentapeptide is formed. This precursor then passes from the cytosol to the cytoplasmic membrane, and is converted to the lipid PP-MurNac penta. After a few more reactions a disaccharide is formed and added to a long glycan chain. One transglycoslyation later and peptidoglycan has been increased.

This model is comprised of a vast number of reactions and substrates/enzymes, many of which have not been characterized. We were able to find parameters for the catalytic rates of murD, murE and murF, although not for B. subtilis. As these rates are similar, we extended these values to other enzymes in the pathway. Whilst this may not be correct we aren't attempting to intricately quantify the amount of peptidoglycan, and so should not be a problem. Substrates involved in but not formed by the peptidoglycan biosynthesis pathway were kept high so as not to be rate limiting, and were chosen as the concentrations at which the enzymes best functioned in vivo as described by sscacadckac. Enzyme concentrations were set to a few millimolar and the simulation was done as an ode, as we weren't worried about how precise the model was.