Take it to the limit (one more time)

The postdoc who's been working on the perl simulation of USS evolution now has everything working, including the equilibrium test that I said I would post about but didn't. I'll have to describe this test later (or maybe she will) because here I want to comment on a different issue.

One problem with simulating biological processes is that they need to be very unrealistic, otherwise they can't be made to work at all. One of the unrealistic features of our USS simulation is the need to use unrealistically high mutation rates, as much as a million times higher than real mutation rates. It's possible to simulate lower mutation rates, but each simulation would take years to run, rather than minutes or hours.

It's tempting to think that the most realistic simulations will be the most informative, but that's far from universally true. In most models the effects of changing variables, and the interactions between variables, are at least as important as the absolute values of the variables. In our USS model, we want to understand how changing each of the following affects the equilibrium numbers of perfect and imperfect USS sequences:
  1. the mutation rate;
  2. the strength of the uptake bias;
  3. the ratio of mutation rate to uptake bias.
The most realistic simulations would use very low mutation rates, and they'd take a lllooonnngg time to run. But we may be able to identify the patterns of change by doing simulations with much higher rates.

For convenience we could start with a mutation rate and uptake bias that let the simulation reach equilibrium in 20 minutes, and test higher and lower rates, keeping the bias constant. Does a higher rate mean more USS or fewer? How big is the effect of a 2-fold change in rate? A 4-fold change? A 10-fold change? What happens to the imperfect USS? Next we can keep the rate constant and vary the bias. Does a 2-fold higher bias give more USS? How many more? What about a 2-fold lower bias? What happens if we change both mutation rate and uptake bias, keeping the ratio constant?

Such experiments can suggest the fundamental properties of the system even though they use very unrealistic values. Once this is done, there are several more directions to go.

First, we don't need to just assume that the properties we've found are likely to hold at more realistic values of rate and bias. Once we understand what happens when mutation rates are unrealistically high, we can think clearly about whether this should change when mutation rates are much lower. We can test our thinking by making a prediction of what we should see in runs with lower mutations rates, and then doing a few runs that test our predictions. If we have also tracked how long each simulation took to reach equilibrium, we'll be able to predict how long these more realistic runs should take, letting us decide whether they're practical or not.

Second, we can devise new tests based on our observations. For example, we should start with simple simulations where cells have no bias towards imperfect USS. But once we understand the properties of these simulations, we can introduce this extra bias and see what happens. Better, we can first predict what we think should happen, then test our predictions. These predictions and tests can be done at high mutation rates so they run fast. Once the results are in, we can consider whether we expect the same effects at low mutation rates, and test key predictions.

Once we're fairly sure of the results, we can work on polishing the paper that will report them. And because paper-polishing takes a long time, we can have a few key simulations running at very low mutations rates while we do this. The final paper will either report that the very low rates gave the expected result, or consider why they didn't.

2 comments:

  1. Stacy said...
    Hey Rosie,
    I have tried to respond to you but the web says I am unauthorized to communicate with you. I agree, alot of "dog people" overlook the people factor in uncomfortable situations. The dog thing just led into a story I have told several times over. People from Mexico not only speak a romance language, they live from the heart. Yes, there are good and bad people EVERYWHERE. I just needed a lead in to write this observation on humanity and culture; now maybe I wont have to keep telling this story over and over. I left out the part about the men in suits in the Denver airport yelling at me to "control my kids" on our jouney back to the states. I really cried then; maybe some other Blog. You are a Scientist, eh? Do you read any Margaret Atwood? How about Oryx and Crake?

    10/08/2006 7:07 PM
    Also, what is the timing differntial in the accelerated process of mutation to what would be normal? Are we talking years, months or days? My totally right brained self is curious.

    ReplyDelete
  2. Hi Stacy!

    Yes I'm a scientist, and I've read just about all of Margaret Atwood's novels, including Oryx and Crake.

    The mutation rate thing is as follows. Bacteria make about one mistake (= one wrong base = one mutation) every time they copy their DNA, which is about once a day. But the DNA is made up of about 5,000,000 base, so each base has only a 1/5,000,000 chance of mutating each time. That's a very slow rate.

    In our model we pretend that each base has a 1/100 or 1/1000 chance of mutating each time - thats 5000 to 50,000 times faster than real bacteria. And our computer also copies the pretend DNA much faster than a real cell would copy its DNA, in a few seconds rather than a day. So rather than taking millions of years to evolve like real bacteria, our pretend bacteria can evolve in a few hours.

    Rosie

    ReplyDelete

Markup Key:
- <b>bold</b> = bold
- <i>italic</i> = italic
- <a href="http://www.fieldofscience.com/">FoS</a> = FoS