1 / 48

Michael Yudkin Kellogg College Oxford

Initiation of differential gene expression in sporulating Bacillus subtilis – a sceptical biochemist looks at mathematical modelling. Michael Yudkin Kellogg College Oxford. Part I. The biological system. Sporulation in Bacillus subtilis. Sporulation cycle. Vegetative cycle.  F.  E.

ace
Download Presentation

Michael Yudkin Kellogg College Oxford

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Initiation of differential gene expression in sporulating Bacillus subtilis –a sceptical biochemist looks at mathematical modelling Michael Yudkin Kellogg College Oxford

  2. Part I. The biological system

  3. Sporulation in Bacillus subtilis Sporulation cycle Vegetative cycle

  4. F E Establishing differential gene expression

  5. Three proteins act to regulate sF: • SpoIIAB (AB), a homodimer of 33 kDa; • SpoIIAA (AA), a monomer of 13 kDa; • SpoIIE (IIE), a multidomain protein of 91 kDa.

  6. AB can engage in three interactions: • It can make a complex with sF in the presence of ATP; • It can make a complex with AA in the presence of ADP; • It can use ATP to phosphorylate AA on Ser-58, yielding AA-P. AA is the substrate of the AB kinase. IIE hydrolyses AA-P back to AA.

  7. Default situation –pre-divisional cell and mother cell •sF is in a sF·AB·ATP complex, with one molecule of sF bound to an AB dimer; • AA is phosphorylated to AA-P, which cannot interact with AB; • IIE is absent or almost absent.

  8. Release of sF in the prespore • IIE is made and hydrolyses AA-P to AA. • The AA attacks the sF·AB·ATP complex. • sF is liberated (“induced release”). • AA is phosphorylated to AA-P; AB·ATP is converted to AB·ADP.

  9. F F F AA AA AA AA AB AB AB AB AB AB AB AB Cartoon of induced release

  10. Induced release and its consequences AA + sF·AB·ATP → sF + AA·AB·ATP AA·AB·ATP → AA-P + AB·ADP AB·ADP+ATP → AB·ATP + ADP sF + AB·ATP → sF·AB·ATP

  11. How is the activation of sF maintained?

  12. ATP AA Pi AB IIE ADP AA-P H2O Cycling of SpoIIAA Summary: ATP + H2O → AA + Pi A wasteful cycle? Can the prespore slow down the cycle to diminish the waste of ATP?

  13. Time course of phosphorylation of SpoIIAA

  14. The time course shows that, after one round of phosphorylation, the enzyme has to return from an inactive to an active form. This return clearly includes a slow step. Our experiments have identified the slow step as the loss of ADP from AB·ADP. But we know from previous work that the loss of ADP from AB·ADP is usually extremely rapid. Why is this case different?

  15. Effect of adding ADP on the phosphorylation of AA

  16. We can account for the results by suggesting that an interaction with AA changes the conformation of AB to AB*. • AB* binds exceptionally tightly to ADP. • AB*·ADP loses its ADP into the solution extremely slowly, and so is long-lived. • Unlike AB·ADP, AB*·ADP can interact with AA to make AA·AB*·ADP. • The formation of AA·AB*·ADP dramatically slows the phosphorylation of AA.

  17. Conclusions so far • While phosphorylating AA, AB accumulates in the form AA·AB*·ADP. • The formation of this complex greatly slows the phosphorylation reaction. • So long as AB is phosphorylating AA, it cannot interact with sF. • Prolonging the phosphorylation reaction therefore serves to maintain the activity of sF.

  18. Phosphorylation reaction scheme ATP AB AB·ATP AA ADP AB·ADP AA·AB*·ATP 1 sec-1 AA·AB*·ADP AB*·ADP AA-P·AB*·ADP AA AA-P

  19. Structure of SpoIIAB

  20. Conformational change in AB AB·ATP AA·AB*·ADP We suggest that the conformational change is a closing of the ATP-lid when AA binds

  21. ATP ADP 1 sec-1

  22. Source of dephosphorylated AA • Induced release of sF from the sF·AB·ATP complex depends on dephosphorylation of AA-P by IIE. • Maintenance of free sF over the necessary period also depends on IIE activity. • Where is the IIE activity located?

  23. IIE is a membrane-bound enzyme, found at the asymmetric septum. IIE Is IIE activity confined to the prespore? If so, how?

  24. Part II. The mathematical model

  25. Mathematical models(a sceptical biochemist’s view) • How do they work? • What are they for? • Why are biochemists sceptical of them? • How do we know if our model is any good?

  26. dA = k1·B + k−2·C− k−1·A− k2·A dt How do models work? • We translate on and off rates of individual interactions into a set of equations, e.g. A k1 k2 k−1 k−2 k−3 C B k3

  27. The kinetic constants (k1, k2 etc.) and the concentrations of the components ([A], [B] etc.) are the parameters. • These parameters are used to produce a set of linked differential equations. • These linked equations constitute the model.

  28. What are models for? I: Qualitative aspects • A mathematical model should include all the components of a system and all of their interactions. • If the model fails to simulate the known behaviour of the system, it is likely that one (or more than one) of the interactions has been omitted.

  29. What are models for? II: Quantitative aspects • Genetics and biochemistry identify the components in a system and show how these components interact. • But to get a quantitative view of the outcome of such interactions, a verbal description is not enough.

  30. The need for a quantitative view: IIE activity in the prespore • As we saw earlier, IIE activity is located on the asymmetric septum. • The prespore is ~ 5-fold smaller in volume than the mother cell. • So the effective concentration of IIE is 5-fold higher in the prespore. • Does this difference account for the specificity of gene expression?

  31. The need for a quantitative view: the sF/sAparadox • In the prespore sF has to displace sA from the core RNA polymerase. • But the affinity of core polymerase for sF is 25-fold less than for sA, and the concentration of sF is only 2-fold higher. So how can sF displace sA?

  32. Why are biochemists sceptical of mathematical models? The number of interactions involved in a regulatory scheme is often large. If we want to make the mathematics come out right, we may be tempted: • to change the parameters without regard to whether they make physical sense; • to increase the number of interactions without regard to the principle of parsimony.

  33. How do we know if our model is any good? • The values for all the parameters should be justifiable (because they have been measured explicitly, or because they are similar to those measured in analogous systems); • The model should make predictions that we can test experimentally.

  34. The sF regulatory interactions AB + ATP → AB·ATP sF + AB·ATP → sF·AB·ATP AA + sF·AB·ATP → sF + AA·AB*·ATP AA·AB*·ATP → AA-P + AB*·ADP AA-P + H2O → AA + Pi AA + AB*·ADP → AA·AB*·ADP AB·ADP+ATP → AB·ATP + ADP sF + core RNAP → sF·holoenzyme

  35. We have previously measured the kinetic constants for these interactions. • We have also measured the concentrations of the intermediates. • We can make a plausible estimate of the kinetic constants for the conformational changes in AB. • We thus have a set of parameters, with which we can construct a model – Model 1. • We now use Model 1 to make verifiable predictions, both qualitative and quantitative.

  36. Predicted formation of sF-holoenzyme(Model 1)

  37. Evidently some factor is missing from our model. What could it be? • AB is a dimer, and we know that it can undergo conformational changes; so maybe AB is an allosteric protein. • If AB is allosteric, it is possible that AA binds to it cooperatively. • Cooperativity can now be included in the model, to give Model 2.

  38. Sporulation Model 2

  39. Predicted formation of sF-holoenzyme (Models 1 and 2)

  40. Binding of AA to AB – experiment and simulation

  41. Other simulations from Model 2 The model can successfully simulate results obtained in vitro from experiments on: • Binding of AA to AB·ADP • Binding of sF to AB·ATP • Disruption of sF·AB·ATP complexes by AA • Rebinding of sF to AB·ATP as AA is phosphorylated • Response of this rebinding to IIE • Time course of phosphorylation of AA.

  42. Quantitative aspects of Model 2: IIE activity in the prespore • As we saw earlier, IIE activity is located on the asymmetric septum. • The prespore is ~ 5-fold smaller in volume than the mother cell. • So the effective concentration of IIE is 5-fold higher in the prespore. • Does this difference account for the specificity of gene expression?

  43. Predicted release of sF

  44. Predicted release of sF IIE+AA-P IIE

  45. Quantitative aspects of Model 2: the sF/sAparadox • In the prespore sF has to displace sA from the core RNA polymerase. • But the affinity of core polymerase for sF is 25-fold less than for sA, and the concentration of sF is only 2-fold higher. So how can sF displace sA?

  46. Predicted formation of sA-holoenzyme and sF-holoenzyme

  47. Conclusions of a convinced sceptic When used judiciously and critically, mathematical modelling is an extremely valuable technique. Mathematical modelling can supply insights that could not be reached by any other method now available.

  48. The model was developed by Joanna Clarkson and Dagmar Iber, using data obtained in my lab by the following: Kyung-Tai Min Mahmoud Najafi Thierry Magnin Matt Lord Daniela Barillà Brian Lee Isabelle Lucet Jwu-Ching Shu Joanna Clarkson Acknowledgments

More Related