Cooperative transport by small teams of molecular motors

Molecular motors power directed transport of cargoes within cells. Even if a single motor is sufficient to transport a cargo, motors often cooperate in small teams. We discuss the cooperative cargo transport by several motors theoretically and explore some of its properties. In particular we emphasize how motor teams can drag cargoes through a viscous environment.Comment: 9 pages, 1 figure, uses ws-brl.cls, presented at Bio-Systems conference, Berlin, June 200

Tug-of-war as a cooperative mechanism for bidirectional cargo transport by molecular motors

Intracellular transport is based on molecular motors that pull cargos along cytoskeletal filaments. One motor species always moves in one direction, e.g. conventional kinesin moves to the microtubule plus end, while cytoplasmic dynein moves to the microtubule minus end. However, many cellular cargos are observed to move bidirectionally, involving both plus-end and minus-end directed motors. The presumably simplest mechanism for such bidirectional transport is provided by a tug-of-war between the two motor species. This mechanism is studied theoretically using the load-dependent transport properties of individual motors as measured in single-molecule experiments. In contrast to previous expectations, such a tug-of-war is found to be highly cooperative and to exhibit seven different motility regimes depending on the precise values of the single motor parameters. The sensitivity of the transport process to small parameter changes can be used by the cell to regulate its cargo traffic.Comment: 17 pages, latex, 11 figures, 4 tables, includes Supporting Informatio

Structural Analysis and Stochastic Modelling Suggest a Mechanism for Calmodulin Trapping by CaMKII

Activation of CaMKII by calmodulin and the subsequent maintenance of constitutive activity through autophosphorylation at threonine residue 286 (Thr286) are thought to play a major role in synaptic plasticity. One of the effects of autophosphorylation at Thr286 is to increase the apparent affinity of CaMKII for calmodulin, a phenomenon known as “calmodulin trapping”. It has previously been suggested that two binding sites for calmodulin exist on CaMKII, with high and low affinities, respectively. We built structural models of calmodulin bound to both of these sites. Molecular dynamics simulation showed that while binding of calmodulin to the supposed low-affinity binding site on CaMKII is compatible with closing (and hence, inactivation) of the kinase, and could even favour it, binding to the high-affinity site is not. Stochastic simulations of a biochemical model showed that the existence of two such binding sites, one of them accessible only in the active, open conformation, would be sufficient to explain calmodulin trapping by CaMKII. We can explain the effect of CaMKII autophosphorylation at Thr286 on calmodulin trapping: It stabilises the active state and therefore makes the high-affinity binding site accessible. Crucially, a model with only one binding site where calmodulin binding and CaMKII inactivation are strictly mutually exclusive cannot reproduce calmodulin trapping. One of the predictions of our study is that calmodulin binding in itself is not sufficient for CaMKII activation, although high-affinity binding of calmodulin is

Methods for simulating actin filament dynamics in MCell

Learning is mediated by activity-dependent changes in synaptic strength that rely on Calcium-dependent signalling in the postsynaptic neuron. In the hippocampus, synapses are located on dendritic spines, small mushroom-like structures that grow from the dendrite. Dendritic spines form micro-compartments that are - to some extent - chemically isolated from the rest of the dendrite due to the long and narrow spine neck limiting diffusion in and out. The shape and the physical structure of a dendritic spine are determined by the actin filaments that form the cellular cytoskeleton. Long-term potentiation (LTP) of the synapse leads to remodelling of the polymeric actin cytoskeleton, which increases the size of the dendritic spine. This process is regulated by interactions between actin and other postsynaptic proteins [1]. Here, addressing the lack of spatiotemporal models that combine structure with biochemical signalling, we present a set of methods that allows the modelling of the biochemistry and spatial dynamics of actin filament remodelling using the spatial stochastic simulator MCell [2, 3, 4]. These methods allow us to recapitulate the main events relevant for actin filament remodelling, including polymerisation, depolymerisation, branching and severing followed by filament displacement in space. They rely on the ability of MCell to model multi-state complex molecules [5] and on a system of virtual tags to label states of actin subunits. The ”complex molecule” feature of MCell 3.3 was utilised by defining an immobile three-dimensional matrix with 8000 subunits, each belonging to one of the defined states. This provides a frame in which some of the subunits are occupied by actin monomers representing the cytoskeleton, while the rest are vacant. In response to biochemical reactions with freely diffusing modifier proteins, the state of the subunits can be changed allowing the rearrangement of the filaments. In addition, with the help of a tagging system for subunits, the monomers forming a single filament can change their location inside the immobile matrix in corpore by occupying previously vacant slots. In the future, the model could be extended by incorporating more elements from signalling pathways and modelling filament bundles present in non-stimulated synapses. As such, the model can provide insights into how the actin cytoskeleton interacts with postsynaptic proteins that mediate LTP

Multi-Stage Modeling of the Kinetics of Activation of CaMKII

Ca 2+ /calmodulin-dependent protein kinase 2 (CaMKII) plays an important role in induction of long-term potentiation and formation of memory. It is abundant in dendritic spines, and is activated when Ca 2+ flows into the postsynaptic cytosol through open NMDA-type glutamate receptors. Its function is fine-tuned through interaction with other proteins as well as through subunit interactions and regulatory autophosphorylation. We have undertaken a multi-stage project to study the critical kinetics of activation of CaMKII in the spine by combining modeling and experimental studies. We are using computational modeling and simulations on various platforms, coupled with biochemical experiments in vitro, and eventually in vivo, to understand CaMKII regulation. The project includes the following steps: 1. Determining the parameters governing activation of a monomeric subunit. The CaMKII holoenzyme is a large dodecamer of similar, homologous subunits held together by interactions between the association domains located at the carboxyl end of each subunit. Individual, monomeric subunits can be expressed recombinantly by removing the association domain. Computer simulations of activation of monomeric CaMKII by Ca 2+ /calmodulin at both saturating and non-saturating concentrations in a test tube have helped to identify the binding parameters that are most crucial for modeling of regulation of CaMKII and thus have indicated the most useful biochemical assays to measure those parameters (Pepke et al., 2010). We are using these measurements to fine-tune our model of activation of individual catalytic subunits. 2. Building a model of the holoenzyme. Because a CaMKII holoenzyme contains 12 similar subunits, each of which can exist in several states, the holoenzyme can have a large number of state combinations. Thus, modeling the entire holoenzyme requires a computational framework that avoids the ensuing combinatorial complexity. The stochastic simulator MCell provides an elegant, rule-based way of modeling state changes in the CaMKII holoenzyme. 3. Modeling cooperativity that arises from the dodecameric structure of CaMKII. Autophosphorylation at threonine-286, which activates CaMKII subunits, is an inter-subunit event. Thus, it is greatly facilitated by the close proximity of subunits in the holoenzyme. In addition, the subunits within the holoenzyme are arranged as dimers which appears to result in cooperativity in the binding of Ca 2+ /CaM to individual subunits of the dimer (Chao et al., 2010). An accurate model of activation of subunits in the holoenzyme and their autophosphorylation will allow us to explore the effects of cooperativity on CaMKII activation on various time scales. 4. Modeling CaMKII within the context of a postsynaptic spine CaMKII interacts with a variety of other proteins, both in the postsynaptic density (PSD), close to major sources of Ca 2+ influx, and in other parts of the spine. In the fourth stage of this project we plan to implement kinetic models of activation of CaMKII in the context of an MCell model of Ca 2+ influx into a spine upon activation of NMDA-type glutamate receptors (Keller et al., 2008; Keller et al., 2011, submitted). We will explore the effects of different localization and numbers of CaMKII holoenzymes in the spine on CaMKII activation. References: Pepke, S., Kinzer-Ursem, T., Mihalas, S., and Kennedy, M.B. (2010). A dynamic model of interactions of Ca 2+ , calmodulin, and catalytic subunits of Ca 2+ /calmodulin-dependent protein kinase II. PLoS Comput Biol 6, e1000675. Chao, L.H., Pellicena, P., Deindl, S., Barclay, L.A., Schulman, H., and Kuriyan, J. (2010). Intersubunit capture of regulatory segments is a component of cooperative CaMKII activation. Nat Struct Mol Biol 17, 264-272. Keller, D.X., Franks, K.M., Bartol, T.M., Jr., and Sejnowski, T.J. (2008). Calmodulin activation by calcium transients in the postsynaptic density of dendritic spines. PLoS ONE 3, e2045. Keller, D.X., Bartol, T.M., Kinney, J.P, Kennedy, M.B., Bajaj, C., Harris, K.M., and Sejnowski, T.J. Regulation of synaptic calcium transients in reconstructed dendritic spines of hippocampal CA1 pyramidal neurons, submitted

Calcium Input Frequency, Duration and Amplitude Differentially Modulate the Relative Activation of Calcineurin and CaMKII

NMDA receptor dependent long-term potentiation (LTP) and long-term depression (LTD) are two prominent forms of synaptic plasticity, both of which are triggered by post-synaptic calcium elevation. To understand how calcium selectively stimulates two opposing processes, we developed a detailed computational model and performed simulations with different calcium input frequencies, amplitudes, and durations. We show that with a total amount of calcium ions kept constant, high frequencies of calcium pulses stimulate calmodulin more efficiently. Calcium input activates both calcineurin and Ca2+/calmodulin-dependent protein kinase II (CaMKII) at all frequencies, but increased frequencies shift the relative activation from calcineurin to CaMKII. Irrespective of amplitude and duration of the inputs, the total amount of calcium ions injected adjusts the sensitivity of the system to calcium input frequencies. At a given frequency, the quantity of CaMKII activated is proportional to the total amount of calcium. Thus, an input of a small amount of calcium at high frequencies can induce the same activation of CaMKII as a larger amount, at lower frequencies. Finally, the extent of activation of CaMKII signals with high calcium frequency is further controlled by other factors, including the availability of calmodulin, and by the potency of phosphatase inhibitors

Motility states of molecular motors engaged in a stochastic tug-of-war

Intracellular transport is mediated by molecular motors that pull cargos along cytoskeletal filaments. Many cargos move bidirectionally and are transported by two teams of motors which move into opposite directions along the filament. We have recently introduced a stochastic tug-of-war model for this situation. This model describes the motion of the cargo as a Markov process on a two-dimensional state space defined by the numbers of active plus and active minus motors. In spite of its simplicity, this tug-of-war model leads to a complex dependence of the cargo motility on the motor parameters. We present new numerical results for the dependence on the number of involved motors. In addition, we derive a simple and intuitive sharp maxima approximation, from which one obtains the cargo motility state from only four simple inequalities. This approach provides a fast and reliable method to determine the cargo motility for a given experimental system.Comment: 23 pages, latex, 12 figures, 1 tabl

Biophysical properties of presynaptic short-term plasticity in hippocampal neurons: insights from electrophysiology, imaging and mechanistic models

Hippocampal neurons show different types of short-term plasticity (STP). Some exhibit facilitation of their synaptic responses and others depression. In this review we discuss presynaptic biophysical properties behind heterogeneity in STP in hippocampal neurons such as alterations in vesicle priming and docking, fusion, neurotransmitter filling and vesicle replenishment. We look into what types of information electrophysiology, imaging and mechanistic models have given about the time scales and relative impact of the different properties on STP with an emphasis on the use of mechanistic models as complementary tools to experimental procedures. Taken together this tells us that it is possible for a multitude of different mechanisms to underlie the same STP pattern, even though some are more important in specific cases, and that mechanistic models can be used to integrate the biophysical properties to see which mechanisms are more important in specific cases of STP