Cyclic AMP (cAMP) and its main effector Protein Kinase A (PKA) are critical for several aspects of neuronal function including synaptic plasticity. Specificity of synaptic plasticity requires that cAMP activates PKA in a highly localized manner despite the speed with which cAMP diffuses. Two mechanisms have been proposed to produce localized elevations in cAMP, known as microdomains: impeded diffusion, and high phosphodiesterase (PDE) activity. This paper investigates the mechanism of localized cAMP signaling using a computational model of the biochemical network in the HEK293 cell, which is a subset of pathways involved in PKA-dependent synaptic plasticity. This biochemical network includes cAMP production, PKA activation, and cAMP degradation by PDE activity. The model is implemented in NeuroRD: novel, computationally efficient, stochastic reaction-diffusion software, and is constrained by intracellular cAMP dynamics that were determined experimentally by real-time imaging using an Epac-based FRET sensor (H30). The model reproduces the high concentration cAMP microdomain in the submembrane region, distinct from the lower concentration of cAMP in the cytosol. Simulations further demonstrate that generation of the cAMP microdomain requires a pool of PDE4D anchored in the cytosol and also requires PKA-mediated phosphorylation of PDE4D which increases its activity. The microdomain does not require impeded diffusion of cAMP, confirming that barriers are not required for microdomains. The simulations reported here further demonstrate the utility of the new stochastic reaction-diffusion algorithm for exploring signaling pathways in spatially complex structures such as neurons.

Reinforcement learning theorizes that strengthening of synaptic connections in medium spiny neurons of the striatum occurs when glutamatergic input (from cortex) and dopaminergic input (from substantia nigra) are received simultaneously. Subsequent to learning, medium spiny neurons with strengthened synapses are more likely to fire in response to cortical input alone. This synaptic plasticity is produced by phosphorylation of AMPA receptors, caused by phosphorylation of various signalling molecules. A key signalling molecule is the phosphoprotein DARPP-32, highly expressed in striatal medium spiny neurons. DARPP-32 is regulated by several neurotransmitters through a complex network of intracellular signalling pathways involving cAMP (increased through dopamine stimulation) and calcium (increased through glutamate stimulation). Since DARPP-32 controls several kinases and phosphatases involved in striatal synaptic plasticity, understanding the interactions between cAMP and calcium, in particular the effect of transient stimuli on DARPP-32 phosphorylation, has major implications for understanding reinforcement learning. We developed a computer model of the biochemical reaction pathways involved in the phosphorylation of DARPP-32 on Thr34 and Thr75. Ordinary differential equations describing the biochemical reactions were implemented in a single compartment model using the software XPPAUT. Reaction rate constants were obtained from the biochemical literature. The first set of simulations using sustained elevations of dopamine and calcium produced phosphorylation levels of DARPP-32 similar to that measured experimentally, thereby validating the model. The second set of simulations, using the validated model, showed that transient dopamine elevations increased the phosphorylation of Thr34 as expected, but transient calcium elevations also increased the phosphorylation of Thr34, contrary to what is believed. When transient calcium and dopamine stimuli were paired, PKA activation and Thr34 phosphorylation increased compared with dopamine alone. This result, which is robust to variation in model parameters, supports reinforcement learning theories in which activity-dependent long-term synaptic plasticity requires paired glutamate and dopamine inputs. Reinforcement learning, based on the association of a stimulus-triggered movement with a reward, involves changes in connection strength between neurons. Memory storage occurs in the striatum, the input stage of the basal ganglia, when a stimulus or movement signal originating from the cortex and a reward signal originating from the midbrain reach the target striatal cells together. Repetitive pairing of these two signals strengthens the connection between cortical and striatal cells. The strengthening of the connections is caused by activation of biochemical signalling pathways inside the striatal cells. These intracellular signalling pathways are explored in a quantitative computational model describing the biochemical pathways important for reinforcement learning. Lindskog et al.'s study shows that when brief reward and stimuli signals are paired, a stronger response in the intracellular signalling occurs compared with the situation when each signal is given alone. This result illustrates mechanisms whereby paired stimuli, but not unpaired stimuli, can cause learning. Furthermore, the model predicts that the biochemical responses are different after brief stimulation as compared with prolonged stimulation. This result highlights the difficulties in predicting the nonlinear interactions within signalling cascades based on prolonged stimulations, which often are used in biochemical experiments.

Striatal fast-spiking (FS) interneurons are interconnected by gap junctions into sparsely connected networks. As demonstrated for cortical FS interneurons, these gap junctions in the striatum may cause synchronized spiking, which would increase the influence that FS neurons have on spiking by the striatal medium spiny (MS) neurons. Dysfunction of the basal ganglia is characterized by changes in synchrony or periodicity, thus gap junctions between FS interneurons may modulate synchrony and thereby influence behavior such as reward learning and motor control. To explore the roles of gap junctions on activity and spike synchronization in a striatal FS population, we built a network model of FS interneurons. Each FS connects to 30–40% of its neighbors, as found experimentally, and each FS interneuron in the network is activated by simulated corticostriatal synaptic inputs. Our simulations show that the proportion of synchronous spikes in FS networks with gap junctions increases with increased conductance of the electrical synapse; however, the synchronization effects are moderate for experimentally estimated conductances. Instead, the main tendency is that the presence of gap junctions reduces the total number of spikes generated in response to synaptic inputs in the network. The reduction in spike firing is due to shunting through the gap junctions; which is minimized or absent when the neurons receive coincident inputs. Together these findings suggest that a population of electrically coupled FS interneurons may function collectively as input detectors that are especially sensitive to synchronized synaptic inputs received from the cortex.

Calcium plays a role in long-term plasticity by triggering postsynaptic signaling pathways for both the strengthening (LTP) and weakening (LTD) of synapses. Since these are opposing processes, seve...

The ability of neurons to differentially respond to specific temporal and spatial input patterns underlies information storage in neural circuits. One means of achieving spatial specificity is to restrict signaling molecules to particular subcellular compartments using anchoring molecules such as A-Kinase Anchoring Proteins (AKAPs). Disruption of protein kinase A (PKA) anchoring to AKAPs impairs a PKA-dependent form of long term potentiation (LTP) in the hippocampus. To investigate the role of localized PKA signaling in LTP, we developed a stochastic reaction-diffusion model of the signaling pathways leading to PKA activation in CA1 pyramidal neurons. Simulations investigated whether the role of anchoring is to locate kinases near molecules that activate them, or near their target molecules. The results show that anchoring PKA with adenylyl cyclase (which produces cAMP that activates PKA) produces significantly greater PKA activity, and phosphorylation of both inhibitor-1 and AMPA receptor GluR1 subunit on S845, than when PKA is anchored apart from adenylyl cyclase. The spatial microdomain of cAMP was smaller than that of PKA suggesting that anchoring PKA near its source of cAMP is critical because inactivation by phosphodiesterase limits diffusion of cAMP. The prediction that the role of anchoring is to colocalize PKA near adenylyl cyclase was confirmed by experimentally rescuing the deficit in LTP produced by disruption of PKA anchoring using phosphodiesterase inhibitors. Additional experiments confirm the model prediction that disruption of anchoring impairs S845 phosphorylation produced by forskolin-induced synaptic potentiation. Collectively, these results show that locating PKA near adenylyl cyclase is a critical function of anchoring. The hippocampus is a part of the cerebral cortex involved in formation of certain types of long term memories. Activity-dependent change in the strength of neuronal connections in the hippocampus, known as synaptic plasticity, is one mechanism used to store memories. The ability to form crisp and distinguishable memories of different events implies that learning produces plasticity of specific and distinct subsets of synapses within each neuron. Synaptic activity leads to production of intracellular signaling molecules, which ultimately cause changes in the properties of the synapses. The requirement for synaptic specificity seems incompatible with the diffusibility of intracellular signaling molecules. Anchoring proteins restrict signaling molecules to particular subcellular compartments thereby combating the indiscriminate spread of intracellular signaling molecules. To investigate whether the critical function of anchoring proteins is to localize proteins near their activators or their targets, we developed a stochastic reaction-diffusion model of signaling pathways leading to synaptic plasticity in hippocampal neurons. Simulations demonstrate that colocalizing proteins with their activator molecules is more important due to inactivation mechanisms that limit the spatial extent of the activator molecules.