User:TheLarry/Notebook/Larrys Notebook/2009/08/04

Project name

Main project page Previous entry      Next entry

Agent Base Modeling After talking with Koch yesterday it seems like taking an agent based modeling approach will be more interesting. We'll treat each head as a different object and watch as their interact. To do this i'll need a sub.vi that figures out the rate constants for each situation which is encountered. I thought about making a giant matrix and just looking up terms but this would be difficult and i want to see it work before i make it more complicated. So each head can be in a state of |position; state; nucleotide> which means it is defined with where it is right now, bound/unbound to MT, and if it has a nucleotide currently bound to it. So i'll need to calculate each rate constant for each state that comes across for both heads. There seems to be 12 rate constants. + means that it is moving forward, and - is backwards. {{#widget:Google Spreadsheet |key=t3AqLiWABJ4PwuXi5Zrhc7g|width=300|height=300}} So there it is 12 different rate constants. I think i have to calculate each one for each time the state changes. Most of these will be zero of course. I forgot the kinesin coming off the MT in general, but that should happen if the other foot disassociates while the other is in the air. So there are 16 different states possible (2 positions * 2 states * 4 nucleotide = 16 possible). This is why i thought about a matrix. But the elements in the matrix would have to be arrays of the 12 k's? I don't know and am not gonna think about it now. So let's get started ATP Binding Strain needs to be relieved. This is a bit controversial. Block says strain is relieved when phosphate is released in Kinesin Motor Mechanics: Binding, Stepping, Tracking, Gating and Limping Also said strain is relieved when the other foot unbinds like in Kinesin: World's Tiniest Biped by Charles Asbury. Asbury uses an earlier Block paper as a source for this but since Block says it later that phosphate release is needed i'll go with him. Head needs to be empty Head needs to be in front Head needs to be bound ATP released ATP needs to be bound not sure what else have to get back to this ATP Hydrolysis Foot must be behind ATP must be attached ADP-P phosphorylation Foot needs to be in ADP-P state Pi release Needs to be in ADP-P state There needs to be another state. like this needs to be behind or something but i can't find anything that states that explicitly Pi enter Needs to be in ADP state Unbind from MT Needs to be behind ATP needs to be in front This creates a forward strain that pulls the foot forward This might not be so much a pre-requisite as it increases the rate constant Needs to have ADP attached Needs to be bound Bind to MT needs to be unbound Diffuse forward needs to be unbound needs to be behind Diffuse backward needs to be unbound needs to be in front ADP release ADP needs to be bound to foot Needs to be in front needs to be bound other foot needs to be bound ADP enter foot needs to be empty Tomorrow I have completed what Koch suggested, and a base program is now completed. I am gonna put in some fake values to see if it is actually working tomorrow. I'll try to make it fit the flow. it should be fun. and by fun i mean crazy annoying to put in numbers in like 20 3 element arrays. but i am curious to see it work. Just read Theoretical model for motility and processivity of two-headed molecular motors. It is a model based on 3 object diffusion with limiting factors that they have to be close and some rate constants have to be small and shit like that. But this goes into long discussion about rate constants into what i call weird states in my program. Most of the time he just says that they should be small. Kinesin Motors: No Strain, No Gain: this article is a quick review of some work Yildiz did. It basically is a good review of how strain affects processivity.