A Dynamic Model of Interactions of Ca^(2+), Calmodulin, and Catalytic Subunits of Ca^(2+)/Calmodulin-Dependent Protein Kinase II

During the acquisition of memories, influx of Ca^(2+) into the postsynaptic spine through the pores of activated N-methyl-D-aspartate-type glutamate receptors triggers processes that change the strength of excitatory synapses. The pattern of Ca^(2+) influx during the first few seconds of activity is interpreted within the Ca^(2+)-dependent signaling network such that synaptic strength is eventually either potentiated or depressed. Many of the critical signaling enzymes that control synaptic plasticity, including Ca^(2+)/calmodulin-dependent protein kinase II (CaMKII), are regulated by calmodulin, a small protein that can bind up to 4 Ca^(2+) ions. As a first step toward clarifying how the Ca^(2+)-signaling network decides between potentiation or depression, we have created a kinetic model of the interactions of Ca^(2+), calmodulin, and CaMKII that represents our best understanding of the dynamics of these interactions under conditions that resemble those in a postsynaptic spine. We constrained parameters of the model from data in the literature, or from our own measurements, and then predicted time courses of activation and autophosphorylation of CaMKII under a variety of conditions. Simulations showed that species of calmodulin with fewer than four bound Ca^(2+) play a significant role in activation of CaMKII in the physiological regime, supporting the notion that processing ofCa^(2+) signals in a spine involves competition among target enzymes for binding to unsaturated species of CaM in an environment in which the concentration of Ca^(2+) is fluctuating rapidly. Indeed, we showed that dependence of activation on the frequency of Ca^(2+) transients arises from the kinetics of interaction of fluctuating Ca^(2+) with calmodulin/CaMKII complexes. We used parameter sensitivity analysis to identify which parameters will be most beneficial to measure more carefully to improve the accuracy of predictions. This model provides a quantitative base from which to build more complex dynamic models of postsynaptic signal transduction during learning

Phase-Locked Signals Elucidate Circuit Architecture of an Oscillatory Pathway

This paper introduces the concept of phase-locking analysis of oscillatory cellular signaling systems to elucidate biochemical circuit architecture. Phase-locking is a physical phenomenon that refers to a response mode in which system output is synchronized to a periodic stimulus; in some instances, the number of responses can be fewer than the number of inputs, indicative of skipped beats. While the observation of phase-locking alone is largely independent of detailed mechanism, we find that the properties of phase-locking are useful for discriminating circuit architectures because they reflect not only the activation but also the recovery characteristics of biochemical circuits. Here, this principle is demonstrated for analysis of a G-protein coupled receptor system, the M3 muscarinic receptor-calcium signaling pathway, using microfluidic-mediated periodic chemical stimulation of the M3 receptor with carbachol and real-time imaging of resulting calcium transients. Using this approach we uncovered the potential importance of basal IP3 production, a finding that has important implications on calcium response fidelity to periodic stimulation. Based upon our analysis, we also negated the notion that the Gq-PLC interaction is switch-like, which has a strong influence upon how extracellular signals are filtered and interpreted downstream. Phase-locking analysis is a new and useful tool for model revision and mechanism elucidation; the method complements conventional genetic and chemical tools for analysis of cellular signaling circuitry and should be broadly applicable to other oscillatory pathways

Simulation of Postsynaptic Glutamate Receptors Reveals Critical Features of Glutamatergic Transmission

Activation of several subtypes of glutamate receptors contributes to changes in postsynaptic calcium concentration at hippocampal synapses, resulting in various types of changes in synaptic strength. Thus, while activation of NMDA receptors has been shown to be critical for long-term potentiation (LTP) and long term depression (LTD) of synaptic transmission, activation of metabotropic glutamate receptors (mGluRs) has been linked to either LTP or LTD. While it is generally admitted that dynamic changes in postsynaptic calcium concentration represent the critical elements to determine the direction and amplitude of the changes in synaptic strength, it has been difficult to quantitatively estimate the relative contribution of the different types of glutamate receptors to these changes under different experimental conditions. Here we present a detailed model of a postsynaptic glutamatergic synapse that incorporates ionotropic and mGluR type I receptors, and we use this model to determine the role of the different receptors to the dynamics of postsynaptic calcium with different patterns of presynaptic activation. Our modeling framework includes glutamate vesicular release and diffusion in the cleft and a glutamate transporter that modulates extracellular glutamate concentration. Our results indicate that the contribution of mGluRs to changes in postsynaptic calcium concentration is minimal under basal stimulation conditions and becomes apparent only at high frequency of stimulation. Furthermore, the location of mGluRs in the postsynaptic membrane is also a critical factor, as activation of distant receptors contributes significantly less to calcium dynamics than more centrally located ones. These results confirm the important role of glutamate transporters and of the localization of mGluRs in postsynaptic sites in their signaling properties, and further strengthen the notion that mGluR activation significantly contributes to postsynaptic calcium dynamics only following high-frequency stimulation. They also provide a new tool to analyze the interactions between metabotropic and ionotropic glutamate receptors

Identification of Key Processes that Control Tumor Necrosis Factor Availability in a Tuberculosis Granuloma

Tuberculosis (TB) granulomas are organized collections of immune cells comprised of macrophages, lymphocytes and other cells that form in the lung as a result of immune response to Mycobacterium tuberculosis (Mtb) infection. Formation and maintenance of granulomas are essential for control of Mtb infection and are regulated in part by a pro-inflammatory cytokine, tumor necrosis factor-α (TNF). To characterize mechanisms that control TNF availability within a TB granuloma, we developed a multi-scale two compartment partial differential equation model that describes a granuloma as a collection of immune cells forming concentric layers and includes TNF/TNF receptor binding and trafficking processes. We used the results of sensitivity analysis as a tool to identify experiments to measure critical model parameters in an artificial experimental model of a TB granuloma induced in the lungs of mice following injection of mycobacterial antigen-coated beads. Using our model, we then demonstrated that the organization of immune cells within a TB granuloma as well as TNF/TNF receptor binding and intracellular trafficking are two important factors that control TNF availability and may spatially coordinate TNF-induced immunological functions within a granuloma. Further, we showed that the neutralization power of TNF-neutralizing drugs depends on their TNF binding characteristics, including TNF binding kinetics, ability to bind to membrane-bound TNF and TNF binding stoichiometry. To further elucidate the role of TNF in the process of granuloma development, our modeling and experimental findings on TNF-associated molecular scale aspects of the granuloma can be incorporated into larger scale models describing the immune response to TB infection. Ultimately, these modeling and experimental results can help identify new strategies for TB disease control/therapy

Akinetic view of GPCR allostery and biased agonism

G-protein-coupled receptors (GPCRs) are one of the most tractable classes of drug targets. These dynamic proteins can adopt multiple active states that are linked to distinct functional outcomes. Such states can be differentially stabilized by ligands interacting with the endogenous agonist-binding orthosteric site and/or by ligands acting via spatially distinct allosteric sites, leading to the phenomena of 'biased agonism' or 'biased modulation'. These paradigms are having a major impact on modern drug discovery, but it is becoming increasingly apparent that 'kinetic context', at the level of both ligand-receptor and receptor-signal pathway kinetics, can have a profound impact on the observation and quantification of these phenomena. The concept of kinetic context thus represents an important new consideration that should be routinely incorporated into contemporary chemical biology and drug discovery studies of GPCR bias and allostery