Latest revision as of 03:27, 3 October 2015

model

3 Signaling module

3.1 Introduction

In our laboratory work, we performed CPP tests to explore the impact of downregulating MOR protein on mouse behavior after morphine administration, which is the ultimate goal of our project. In this module, computational and systems biology approaches were applied to examine the root of behavior changes quantitatively at the molecular level. The most important brain reward circuit involves dopamine-containing neurons in the VTA of the midbrain. Morphine can cause indirect excitation of VTA dopamine neurons by reducing inhibitory synaptic transmission mediated by GABAergic neurons [1,2].

We modeled the signaling network to investigate the emergent properties of the reward pathway. By comparing the activation degree of the reward pathway before and after downregulating MOR protein levels, we could have a better mechanistic understanding of drug effects. Although we did not perform any experiment to support this modeling module, the methods and parameters we chose are grounded in literature reports.

Figure 8. Reward pathway of acute morphine administration. We focused on activation of MOR, inhibition of AC and release of GABA vesicles in this module. The reference pathway and figure are adapted from Kyoto Encyclopedia of Genes and Genomes database (KEGG).

3.2 Model methods

We used both deterministic and stochastic models to describe the activation of GPCR and release of GABA. In biological systems, signal transmission occurs primarily through two mechanisms: (i) mass-action laws governing protein synthesis, degradation and interactions; and (ii) standard Michaelis-Menten formulation for reactions catalyzed by enzymes [3].

Broadly, mathematical models of biochemical reactions can be divided into two categories: deterministic systems and stochastic systems [3]. In deterministic models, the change in time of the components’ concentrations is completely determined by specifying the initial and boundary conditions; by contrast, the changes in concentrations of components with respect to time cannot be fully predicted in stochastic models [3]. In the previous two modules, we modeled the delivery device and RNA interference using deterministic models.

3.2.1 Modeling the activation of MOR

MOR belongs to the class A (Rhodopsin) family of heterotrimeric Gi/o protein-coupled receptors [4]. The binding of opioids to MOR activates the G protein, upon which both G-protein α and βγ subunits interact with multiple cellular effector systems. As the first step of signal transmission, the degree of activation of MOR in response to opioid has a direct and far-reaching influence on the behavior of mice.

Deterministic models were applied to describe the biochemical reactions occurring in the diagram below. We used the Matlab Simbiology package to draw the diagram and to design the equation, the details of which are accessible on the uploaded files. This model was created on the basis of work by Bhalla and Iyengar on the activation of glutamate receptor [5].

Figure 9. Reaction schemes for the activation of MOR in the simulation. Reversible reactions are represented as bidirectional arrows; irreversible reactions, as unidirectional arrows. This figure is adapted from the literature [5].

3.2.2 Modeling adenylate cyclase inhibtion

The concentration of second messenger is a significant indicator of excitability of GABAergic neurons. Thus, we chose to simulate cAMP levels and adenylate cyclase (AC) activity to determine the effect of downregulating MOR protein levels on morphine reward signaling networks.

AC1/8 is a type of adenylate cyclases involved in the signaling of the acute morphine reward pathway [6]. When MOR is activated, the disassociated Gα subunit reacts with AC1/8 and subsequently inhibits its activity, leading to a decrease in cellular cAMP levels. The parameters of this model were primarily derived from the literature [5] with slight modifications to fit to the data presented in the literature [7].

Figure 10. Reaction schemes for inhibition of AC in simulation. Reversible reactions are represented as bidirectional arrows, and enzyme reactions are drawn as an arrow with two bends. AC: adenylate cyclase; PDE: phosphodiesterase.

3.2.3 Modeling GABA vesicle releases

A stochastic model was applied to describe the random behavior of neurotransmitter vesicles release [8]. GABA is an important inhibitory neurotransmitter, the level of which directly determines the firing rate of dopamine neurons and other physiological and behavioral statuses. The GABA synaptic vesicle cycle consists of three discrete processes: synthesis of GABA vesicles, docking of GABA vesicles at the inner membrane of presynapses and release of GABA vesicles reacting to a certain signal. The release of GABA vesicles is strictly regulated by cellular signaling networks. When Gi/o is activated and the cellular cAMP level drops, the release of GABA is inhibited. Many complicated mechanisms are involved in the inhibition of GABA release due to activation of Gi/o. Here, we simply studied the action potential-independent pathway of GABA release, through which the release of GABA is directly inhibited by activated Gβγ subunits [9].

Figure 11. Schematic representation of GABA release in which four steps are modeled using mass action law and the stochastic method.

3.2.4 Gillespie’s algorithm

When spatially restricted reactions, such as the release of neurotransmitter vesicles, are studied, the traditional deterministic model is no longer effective for ignoring the discrete nature of the problem [3]. Stochastic models convert reaction rates to probability, which allows users to explore the noise and randomness of signaling networks. A standard algorithm dealing with stochastic model is Gillespie’s algorithm. This algorithm starts with the initial condition for each molecule type in the reaction network. Then, Monte Carlo simulation is applied to generate some random variables and to calculate the smallest time interval in which the reaction will occur [3,10]. Finally, the number of molecules in the reaction network is updated, and the process is repeated.

3.3 Results

The simulation results revealed the kinetics of MOR activation in case and control studies. In the CPP test, the Western blot result demonstrated that the relative level of MOR protein after MOR-siRNA injection was 0.5. Thus, the concentration of MOR protein was set at half of the level in the case study.

The results indicated that almost all the MOR protein is activated in response to morphine. The quantity and action of Gα and βγ subunits highly correlates with the quantity of MOR protein. By downregulating the MOR protein to half of its initial level, we also inhibit approximately half of activated Gα and βγ subunits.

Figure 12. Concentration-time curves for the activation of MOR in response to morphine. A: Control study with the concentration of MOR set at 1 mM. B: Case study with the concentration of MOR set at 0.5 mM due to downregulation by MOR-siRNA. Ga_GTP and Gbg represents activated Gα and βγ subunit, respectively.

The primary effector of activated Gα subunit is AC. The activation degree of AC influences its product cAMP—an important second messenger that indicates the excitability of GABAergic neurons. We now theoretically predicted and compared the inhibition of Gα subunit on AC and the subsequent decrease in cellular cAMP levels in control (wild type) and case (MOR-siRNA injected) studies.

Figure 13. Effect of downregulating MOR protein on AC activity (A) and cellular cAMP levels (B) in response to morphine. The input level of MOR protein is based on the result shown in Figure 10.

Activation of wild type MOR protein inhibited over 25% of AC, and relative cellular cAMP levels dropped below 70%, which is consistent with findings in the literature[7]. The injection of MOR-siRNA reduces the activation quantity of MOR and significantly attenuates the inhibition of AC and decrease in cAMP levels. Maintaining the cellular cAMP level induced by the drug plays a crucial role in blocking reward pathways.

Finally, we explored the relationship between MOR activation and GABA release. The wild type study revealed significant inhibition of GABA vesicles due to activated G βγ subunits. MOR-siRNA counteracted this trend by downregulating MOR protein and activated G βγ subunit levels as depicted in the case study. Maintaining GABA release reduces the excitability and firing rate of dopamine neurons, which is consistent with the expected drug effect on blockage of the reward pathway and explain the behavioral changes observed in the CPP tests.

Figure 14. Stochastic modeling of GABA release. A: Control study with the establishment of mass balance between synthesized, docked and released GABA vesicles. B: Case study with MOR-siRNA injected to attenuate the inhibition of GABA release. C: Wild type study with a normal level of MOR protein activation resulting in inhibition of GABA release. D: Summary of numbers of released and inhibited GABA vesicles in different treatments. The results are presented as the mean±S.D.

3.4 Conclusion and remarks

In this module, we used deterministic and stochastic methods to model the cell signaling network and to predict the blockage of the reward pathway by injecting MOR- siRNA. The simulation results could somewhat explain the behavioral changes observed in the CPP tests (function tests) mechanistically.

3.5 Model equations, variables and parameters

The modeling details of the activation of MOR protein and inhibition of AC are truncated here because we use software to help us design the model and there are too many parameters and equations. We have uploaded relevant source code and files for those individuals interested in exploring the models. However, we want to emphasize that our parameters are all derived from the literature. The modeling of GABA release is inspired by the literature [8] and the parameters were estimated from literature [2]. These parameters, as well as initial conditions, can be accessed in our uploaded files and we selectively list part of them below.

These parameters, as well as initial conditions, can be accessed in our uploaded files and we selectively list part of them below.

3.5.1 Activation of MOR

Model Parameters:

???: Although activation of MOR has not been modeled yet, we use activation of Glutamate receptor, which has been modeled in the literature as an approximation.

Model Equations:

3.5.2 Activation of MOR

Model Parameters:

???: No literature has directly reported binding and disassociation constant of Gi to AC. Therefore, we use the binding and disassociation constant of Gs to AC as an approximation derived from literature[5].

Model Equations:

3.5.3 GABA release

Model Parameters

References:
1.Fields, H.L. and Margolis, E.B. (2015) Understanding opioid reward. Trends in neurosciences, 38, 217-225.
2.Sotomayor, R., Forray, M.I. and Gysling, K. (2005) Acute morphine administration increases extracellular DA levels in the rat lateral septum by decreasing the GABAergic inhibitory tone in the ventral tegmental area. Journal of neuroscience research, 81, 132-139.
3.Eungdamrong, N.J. and Iyengar, R. (2004) Computational approaches for modeling regulatory cellular networks. Trends in cell biology, 14, 661-669.
4.Waldhoer, M., Bartlett, S.E. and Whistler, J.L. (2004) Opioid receptors. Annual Review of Biochemistry, 73, 953-990.
5.Bhalla, U.S. and Iyengar, R. (1999) Emergent properties of networks of biological signaling pathways. Science, 283, 381-387.
6.Nestler, E.J. and Aghajanian, G.K. (1997) Molecular and cellular basis of addiction. Science, 278, 58-63.
7.Charalampous, K.D. and Askew, W.E. (1977) Cerebellar cAMP levels following acute and chronic morphine administration. Can J Physiol Pharmacol, 55, 117-120.
8.Ribrault, C., Sekimoto, K. and Triller, A. (2011) From the stochasticity of molecular processes to the variability of synaptic transmission. Nature reviews. Neuroscience, 12, 375-387.
9.Stephens, G.J. (2009) G-protein-coupled-receptor-mediated presynaptic inhibition in the cerebellum. Trends Pharmacol Sci, 30, 421-430.
10.Gillespie, D.T. (1977) Exact stochastic simulation of coupled chemical reactions. The Journal of Physical Chemistry, 81, 2340-2361.

@@ Line 200: / Line 200: @@
 <body class="index" data-view="index">
 <div class="content">
-<nav>
+<nav style="height:80px">
-<ul>
+<ul style="margin-top:10px">
 <li><a class="button" href="#detail" style="font-weight:bold;font-family:Microsoft YaHei">content</a></li>
 </ul>
@@ Line 241: / Line 241: @@
 <TABLE borderColor=#00ff99 height="100%" width="100%" border=0 style="table-layout:fixed">
      <TR height="100%">
-	<TD width="27%" bgColor=#EEE8AA style="vertical-align:top">
+	<TD width="27%" bgColor=#E6E8FA style="vertical-align:top">
- <ul id="nav" class="dyn" style="height:100px;line-height:100px;margin-top:20%;margin-left:10%;list-style-type:none">
+</br>
-<li><a href="https://2015.igem.org/Team:NJU-China" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Home</a></li></br>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/Team:NJU-China" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Home</a></li>
-<li><a href="https://2015.igem.org/NJU-China-background.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Background</a></li></br>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-background.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Background</a></li>
-<li><a href="https://2015.igem.org/NJU-China-project.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Project</a></li></br>
+<div style="line-height:250%;margin-left:10%" id="main1" onClick="document.all.child0.style.display=(document.all.child0.style.display =='none')?'':'none'" ><a href="#" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Project</a></div>
-<li><a href="#" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Model</a>
+<div id="child0" style="display:none">
-        <ul>
+<ul>
-            <li><a href="https://2015.igem.org/Team:NJU-China/Modeling.html#deliverymodel" style="font-weight:bold;font-family:幼圆;font-size:20px;color:black">Delivery model</a></li>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/Team:NJU-China/Design" style="font-weight:bold;font-family:幼圆;color:black">design</a></li>
-            <li><a href="https://2015.igem.org/Team:NJU-China/Modeling#" style="font-weight:bold;font-family:幼圆;font-size:20px;color:black">RNAi model</a></li>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-project/result.html" style="font-weight:bold;font-family:幼圆;color:black">results</a></li>
-        </ul>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-project/conclusion.html" style="font-weight:bold;font-family:幼圆;color:black">conclusion</a></li>
-</li></br>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-project/future work.html" style="font-weight:bold;font-family:幼圆;color:black">future work</a></li>
-<li style="line-height:250%"><a href="https://2015.igem.org/Team:NJU-China/Practices" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Human Practice</a></li></br>
+</ul>
-<li style="line-height:250%"><a href="https://2015.igem.org/NJU-China-parts.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Parts</a></li></br>
+</div>
-<li style="line-height:250%"><a href="https://2015.igem.org/NJU-China-team.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Team</a></li></br>
+<div style="line-height:250%;margin-left:10%" id="main1" onClick="document.all.child1.style.display=(document.all.child1.style.display =='none')?'':'none'" ><a href="#" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Model</a></div>
-<li style="line-height:250%"><a href="https://2015.igem.org/Team:NJU-China/Attributions" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Attribution</a></li></br>
+<div id="child1" style="display:none">
-<li style="line-height:250%"><a href="https://2015.igem.org/Team:NJU-China/Collaborations" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Colaborations</a></li></br>
+<ul>
-<li style="line-height:250%"><a href="https://2015.igem.org/NJU-China-safty.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Safety</a></li></br>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-model.html" style="font-weight:bold;font-family:幼圆;color:black">Delivery model</a></li>
-<li style="line-height:250%"><a href="#" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Notebook</a>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/Team:NJU-China/RNAi" style="font-weight:bold;font-family:幼圆;color:black">RNAi model</a></li>
-        <ul>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/Team:NJU-China/signaling" style="font-weight:bold;font-family:幼圆;color:black">Signaling</a></li>
-            <li><a href="https://2015.igem.org/NJU-China-notebook.html#methods" style="font-weight:bold;font-family:幼圆;font-size:20px;color:black">Methods</a></li>
+</ul>
-            <li><a href="https://2015.igem.org/NJU-China-notebook.html#protocal" style="font-weight:bold;font-family:幼圆;font-size:20px;color:black">Protocal</a></li>
+</div>
-            <li><a href="https://2015.igem.org/NJU-China-notebook.html#notebook" style="font-weight:bold;font-family:幼圆;font-size:20px;color:black">Notebook</a></li>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/Team:NJU-China/Practices" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Human Practice</a></li>
-        </ul>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-parts.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Parts</a></li>
-</li></br>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-team.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Team</a></li>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-attribution.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Attribution</a></li>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/Team:NJU-China/Collaborations" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Collaborations</a></li>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-safty.html" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Safety</a></li>
+<div style="line-height:250%;margin-left:10%" id="main2" onClick="document.all.child2.style.display=(document.all.child2.style.display =='none')?'':'none'" > <a href="#" style="font-weight:bold;font-family:幼圆;font-size:25px;color:black">Notebook </div>
+<div id="child2" style="display:none">
+<ul>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-notebook.html#methods" style="font-weight:bold;font-family:幼圆;font-size:20px;color:black">Methods</a></li>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-notebook.html#protocal" style="font-weight:bold;font-family:幼圆;font-size:20px;color:black">Protocal</a></li>
+<li style="line-height:250%;margin-left:10%"><a href="https://2015.igem.org/NJU-China-notebook.html#notebook" style="font-weight:bold;font-family:幼圆;font-size:20px;color:black">Notebook</a></li>
+</ul>
+</div>
-</ul>
 	</TD>
      <TD width="73%" bgColor=#FFFFFF style="vertical-align:top;padding-left:80px;padding-right:80px;padding-top:50px;padding-bottom:50px;word-wrap:break-word;">
-	你在火中涅槃，灵魂不死；你在花间长眠，记忆不消。绚烂的阳光是迎接你西去的光环，不知你独自一人去那方可会感到孤单？
+							<h1> 3 Signaling module </h1>
-[一]魂去
-      起初听到你的死，我是愕然的。毕竟听说你的病情好了不少，却不料那只是最后的回光返照。
-      我没有流泪，甚至没有悲伤，脑子里只是混混沌沌的一团。死？那是我们所有人最后的归宿吧。或许于你更是能够减少痛苦，未尝不是一件好事。
+	<h2> 3.1 Introduction </h2>
-      一切被我想得太简单，甚至太过理所当然了。理智，真是一个可怕的局。
-[二]席间
-      小号、大鼓敲个不停，那种声音不得不说是一种噪音。良久声止，开始一个个念悼词。
-      我这时才发现，原来我并没有那么冷漠，缘是城市的喧嚣麻木了我的心。想起那年你的温言软语、慈爱笑容，泪竟是就这样忍不住地从眼眶中溢出了。我撇过头，不愿让人看见，嘴角仍是挂着笑容。
+In our laboratory work, we performed CPP tests to explore the impact of downregulating
-      呵呵，我竟觉着在这葬礼上哭泣都是可耻的。望着周围人谈笑风生，我的心有点寒。到场的那么多人，究竟又有几个是真心的呢？可怜，还是可悲啊？
-      原来人可以铁石心肠这种程度，可毕竟他们中很多人比我和你的关系更近啊！听着那个花钱雇来哭的人呜咽，看着那人热泪纵横，我终是忍不住苦笑，竟是一个毫不相识的外人哭得最伤心。这是讽刺，还是笑话？
+MOR protein on mouse behavior after morphine administration, which is the ultimate goal
-	  [三]送葬
-       你说，为什么他们还能那样淡然，可我偏偏忍不住泪流满面？为什么他们那个惊讶的表情让我竟是感到羞愧？好像自己做了什么惊天动地的事一样？
+of our project. In this module, computational and systems biology approaches were
-       执着一根香，撒下一碗酒，三拜三叩，泪融于酒，我如若痴人。
-       照礼我是不必披麻戴孝的，只需戴着那顶白色的帽子即可。沿着幽僻的小路，我跟着吹着小号开路的队伍走在末尾。
+applied to examine the root of behavior changes quantitatively at the molecular level.
-       风景甚好，埋在那棵大树底下，你会不会成为那树的枝，树的叶？
-[四]感怀
+The most important brain reward circuit involves dopamine-containing neurons in the VTA
-      死生有命，无可更改。谋事在人，成事在天。任你生前浮华万千，到头来不过一抔尘土，一捧青灰。愿逝者安息，生者恬淡。
-      有时候我们看不到自己拥有的美好，直到失去才恍然叹息。逝者已矣，页已掀过，惜福、感恩，生活其实温暖如花。
+of the midbrain. Morphine can cause indirect excitation of VTA dopamine neurons by
-      花开半夏，倾尽繁华，愿你回眸半刹。你可听到我在同你夜话？
-	</TD>
+reducing inhibitory synaptic transmission mediated by GABAergic neurons [1,2].
+<br><br>
+<B>
+We modeled the signaling network to investigate the emergent properties of the reward
+pathway. By comparing the activation degree of the reward pathway before and after
+downregulating MOR protein levels, we could have a better mechanistic understanding of
+drug effects. Although we did not perform any experiment to support this modeling
+module, the methods and parameters we chose are grounded in literature reports.
+</B>
+<br><br>
+<!--插入第八张图--> <img src="https://static.igem.org/mediawiki/2015/e/ec/NJU-China-
+Model_Figure8.jpg" style="width:600px" > <br><br>
+Figure 8. Reward pathway of acute morphine administration. We focused on activation of
+MOR, inhibition of AC and release of GABA vesicles in this module. The reference
+pathway and figure are adapted from Kyoto Encyclopedia of Genes and Genomes database
+(KEGG).<br><br>
+	<h2> 3.2 Model methods </h2>
+<B> We used both deterministic and stochastic models to describe the activation of GPCR
+and release of GABA. </B>
+In biological systems, signal transmission occurs primarily through two mechanisms: (i)
+mass-action laws governing protein synthesis, degradation and interactions; and (ii)
+standard Michaelis-Menten formulation for reactions catalyzed by enzymes [3].
+<br><br>
+Broadly, mathematical models of biochemical reactions can be divided into two
+categories: deterministic systems and stochastic systems [3]. In deterministic models,
+the change in time of the components’ concentrations is completely determined by
+specifying the initial and boundary conditions; by contrast, the changes in
+concentrations of components with respect to time cannot be fully predicted in
+stochastic models [3]. In the previous two modules, we modeled the delivery device and
+RNA interference using deterministic models.
+<br><br>
+	<h3> 3.2.1 Modeling the activation of MOR </h3>
+MOR belongs to the class A (Rhodopsin) family of heterotrimeric Gi/o protein-coupled
+receptors [4]. The binding of opioids to MOR activates the G protein, upon which both
+G-protein α and βγ subunits interact with multiple cellular effector systems. As the
+first step of signal transmission, the degree of activation of MOR in response to
+opioid has a direct and far-reaching influence on the behavior of mice.
+<br><br>
+Deterministic models were applied to describe the biochemical reactions occurring in
+the diagram below. We used the Matlab Simbiology package to draw the diagram and to
+design the equation, the details of which are accessible on the uploaded files. This
+model was created on the basis of work by Bhalla and Iyengar on the activation of
+glutamate receptor [5].
+<br><br>
+<!-- 插入第九张图 -->   <img src="https://static.igem.org/mediawiki/2015/e/e3/NJU-China-
+Model_Figure9.jpg" style="width:600px" > <br><br>
+Figure 9. Reaction schemes for the activation of MOR in the simulation. Reversible
+reactions are represented as bidirectional arrows; irreversible reactions, as
+unidirectional arrows. This figure is adapted from the literature [5]. <br><br>
+	<h3> 3.2.2 Modeling adenylate cyclase inhibtion </h3>
+The concentration of second messenger is a significant indicator of excitability of
+GABAergic neurons. Thus, we chose to simulate cAMP levels and adenylate cyclase (AC)
+activity to determine the effect of downregulating MOR protein levels on morphine
+reward signaling networks.
+<br><br>
+AC1/8 is a type of adenylate cyclases involved in the signaling of the acute morphine
+reward pathway [6]. When MOR is activated, the disassociated Gα subunit reacts with
+AC1/8 and subsequently inhibits its activity, leading to a decrease in cellular cAMP
+levels. The parameters of this model were primarily derived from the literature [5]
+with slight modifications to fit to the data presented in the literature [7].
+<br><br>
+<!-- 插入第十张图 --> <img src="https://static.igem.org/mediawiki/2015/7/78/NJU-China-
+Model_Figure10.jpg" style="width:400px" > <br><br>
+Figure 10. Reaction schemes for inhibition of AC in simulation. Reversible reactions
+are represented as bidirectional arrows, and enzyme reactions are drawn as an arrow
+with two bends. AC: adenylate cyclase; PDE: phosphodiesterase.
+<br><br>
+	<h3> 3.2.3 Modeling GABA vesicle releases </h3>
+A stochastic model was applied to describe the random behavior of neurotransmitter
+vesicles release [8]. GABA is an important inhibitory neurotransmitter, the level of
+which directly determines the firing rate of dopamine neurons and other physiological
+and behavioral statuses. The GABA synaptic vesicle cycle consists of three discrete
+processes: synthesis of GABA vesicles, docking of GABA vesicles at the inner membrane
+of presynapses and release of GABA vesicles reacting to a certain signal. The release
+of GABA vesicles is strictly regulated by cellular signaling networks. When Gi/o is
+activated and the cellular cAMP level drops, the release of GABA is inhibited. Many
+complicated mechanisms are involved in the inhibition of GABA release due to activation
+of Gi/o. Here, we simply studied the action potential-independent pathway of GABA
+release, through which the release of GABA is directly inhibited by activated Gβγ
+subunits [9].
+<br><br>
+<!-- 插入第十一张图 --> <img src="https://static.igem.org/mediawiki/2015/1/16/NJU-China-
+Model_Figure11.jpg" style="width:600px" > <br><br>
+Figure 11. Schematic representation of GABA release in which four steps are modeled
+using mass action law and the stochastic method. <br><br>
+	<h3> 3.2.4 Gillespie’s algorithm </h3>
+When spatially restricted reactions, such as the release of neurotransmitter vesicles,
+are studied, the traditional deterministic model is no longer effective for ignoring
+the discrete nature of the problem [3]. Stochastic models convert reaction rates to
+probability, which allows users to explore the noise and randomness of signaling
+networks. A standard algorithm dealing with stochastic model is Gillespie’s algorithm.
+This algorithm starts with the initial condition for each molecule type in the reaction
+network. Then, Monte Carlo simulation is applied to generate some random variables and
+to calculate the smallest time interval in which the reaction will occur [3,10].
+Finally, the number of molecules in the reaction network is updated, and the process is
+repeated.
+<br><br>
+	<h2> 3.3 Results </h3>
+The simulation results revealed the kinetics of MOR activation in case and control
+studies. In the CPP test, the Western blot result demonstrated that the relative level
+of MOR protein after MOR-siRNA injection was 0.5.
+<B> Thus, the concentration of MOR protein was set at half of the level in the case
+study. </B>
+<br><br>
+<B> The results indicated that almost all the MOR protein is activated in response to
+morphine. The quantity and action of Gα and βγ subunits highly correlates with the
+quantity of MOR protein. By downregulating the MOR protein to half of its initial
+level, we also inhibit approximately half of activated Gα and βγ subunits. </B>
+<br><br>
+<!-- 插入第十二张图 --> <img src="https://static.igem.org/mediawiki/2015/b/b3/NJU-China-
+Model_Figure12.jpg"  > <br><br>
+Figure 12. Concentration-time curves for the activation of MOR in response to morphine.
+A: Control study with the concentration of MOR set at 1 mM. B: Case study with the
+concentration of MOR set at 0.5 mM due to downregulation by MOR-siRNA. Ga_GTP and Gbg
+represents activated Gα and βγ subunit, respectively.
+<br><br>
+The primary effector of activated Gα subunit is AC. The activation degree of AC
+influences its product cAMP—an important second messenger that indicates the
+excitability of GABAergic neurons. We now theoretically predicted and compared the
+inhibition of Gα subunit on AC and the subsequent decrease in cellular cAMP levels in
+control (wild type) and case (MOR-siRNA injected) studies.
+<br><br>
+<!-- 插入第十三张图 --> <img src="https://static.igem.org/mediawiki/2015/d/d4/NJU-China-
+Model_Figure13.jpg" > <br><br>
+Figure 13. Effect of downregulating MOR protein on AC activity (A) and cellular cAMP
+levels (B) in response to morphine. The input level of MOR protein is based on the
+result shown in Figure 10.
+<br><br>
+<B> Activation of wild type MOR protein inhibited over 25% of AC, and relative cellular
+cAMP levels dropped below 70%, which is consistent with findings in the literature[7].
+</B>
+The injection of MOR-siRNA reduces the activation quantity of MOR and significantly
+attenuates the inhibition of AC and decrease in cAMP levels. Maintaining the cellular
+cAMP level induced by the drug plays a crucial role in blocking reward pathways.
+<br><br>
+Finally, we explored the relationship between MOR activation and GABA release. The wild
+type study revealed significant inhibition of GABA vesicles due to activated G βγ
+subunits. MOR-siRNA counteracted this trend by downregulating MOR protein and activated
+G βγ subunit levels as depicted in the case study. Maintaining GABA release reduces
+the excitability and firing rate of dopamine neurons, which is consistent with the
+expected drug effect on blockage of the reward pathway and explain the behavioral
+changes observed in the CPP tests.
+<br><br>
+<!-- 插入第十四张图 --> <img src="https://static.igem.org/mediawiki/2015/5/51/NJU-China-
+Model_Figure14.jpg" style="width:600px" > <br><br>
+Figure 14. Stochastic modeling of GABA release. A: Control study with the establishment
+of mass balance between synthesized, docked and released GABA vesicles. B: Case study
+with MOR-siRNA injected to attenuate the inhibition of GABA release. C: Wild type study
+with a normal level of MOR protein activation resulting in inhibition of GABA release.
+D: Summary of numbers of released and inhibited GABA vesicles in different treatments.
+The results are presented as the mean±S.D.
+<br><br>
+	<h2> 3.4 Conclusion and remarks </h2>
+<B>
+In this module, we used deterministic and stochastic methods to model the cell
+signaling network and to predict the blockage of the reward pathway by injecting MOR-
+siRNA. The simulation results could somewhat explain the behavioral changes observed in
+the CPP tests (function tests) mechanistically.
+<br><br>
+</B>
+	<h2> 3.5 Model equations, variables and parameters </h2>
+The modeling details of the activation of MOR protein and inhibition of AC are
+truncated here because we use software to help us design the model and there are too
+many parameters and equations. We have uploaded relevant source code and files for
+those individuals interested in exploring the models. However, we want to emphasize
+that our parameters are all derived from the literature. The modeling of GABA release
+is inspired by the literature [8] and the parameters were estimated from literature
+[2].
+These parameters, as well as initial conditions, can be accessed in our uploaded files
+and we selectively list part of them below.
+<br><br>
+These parameters, as well as initial conditions, can be accessed in our uploaded files
+and we selectively list part of them below.
+<br><br>
+	<h3> 3.5.1 Activation of MOR </h3>
+Model Parameters:
+<br><br>
+<img src="https://static.igem.org/mediawiki/2015/d/de/NJU-China-model-sig-1.jpg" >
+<br><br>
+???: Although activation of MOR has not been modeled yet, we use activation of
+Glutamate receptor, which has been modeled in the literature as an approximation.
+<br><br>
+Model Equations:
+<table>
+<tr>
+<td>
+<!-- 插入第一张公式 --> <img src="https://static.igem.org/mediawiki/2015/b/ba/NJU-China-
+Equation_Sig_1.jpg" style="width:500px"> </td> <br><br>
+</tr>
+<tr>
+<td>
+<!-- 插入第二张公式 --> <img src="https://static.igem.org/mediawiki/2015/0/0b/NJU-China-
+Equation_Sig_2.jpg" style="width:500px" > </td> <br><br>
+</tr>
+</table>
+	<h3> 3.5.2 Activation of MOR </h3>
+Model Parameters:
+<br><br>
+<img src = "https://static.igem.org/mediawiki/2015/f/f2/NJU-China-model-sig-2.jpg" >
+<br><br>
+???: No literature has directly reported binding and disassociation constant of Gi to
+AC. Therefore, we use the binding and disassociation constant of Gs to AC as an
+approximation derived from literature[5].
+<br><br>
+Model Equations:
+<table>
+ <img src="https://static.igem.org/mediawiki/2015/a/a7/NJU-China-
+Equation_Sig_3.jpg", style="width:500px"> <br><br>
+</table>
+	<h3> 3.5.3 GABA release </h3>
+Model Parameters
+<br><br>
+<img src = "https://static.igem.org/mediawiki/2015/1/19/NJU-China-model-sig-3.jpg" >
+<br><br>
+References: <br>
+.Fields, H.L. and Margolis, E.B. (2015) Understanding opioid reward. Trends in
+neurosciences, 38, 217-225. <br>
+.Sotomayor, R., Forray, M.I. and Gysling, K. (2005) Acute morphine administration
+increases extracellular DA levels in the rat lateral septum by decreasing the GABAergic
+inhibitory tone in the ventral tegmental area. Journal of neuroscience research, 81,
+-139. <br>
+.Eungdamrong, N.J. and Iyengar, R. (2004) Computational approaches for modeling
+regulatory cellular networks. Trends in cell biology, 14, 661-669. <br>
+.Waldhoer, M., Bartlett, S.E. and Whistler, J.L. (2004) Opioid receptors. Annual
+Review of Biochemistry, 73, 953-990. <br>
+.Bhalla, U.S. and Iyengar, R. (1999) Emergent properties of networks of biological
+signaling pathways. Science, 283, 381-387. <br>
+.Nestler, E.J. and Aghajanian, G.K. (1997) Molecular and cellular basis of addiction.
+Science, 278, 58-63. <br>
+.Charalampous, K.D. and Askew, W.E. (1977) Cerebellar cAMP levels following acute and
+chronic morphine administration. Can J Physiol Pharmacol, 55, 117-120. <br>
+.Ribrault, C., Sekimoto, K. and Triller, A. (2011) From the stochasticity of molecular
+processes to the variability of synaptic transmission. Nature reviews. Neuroscience,
+, 375-387. <br>
+.Stephens, G.J. (2009) G-protein-coupled-receptor-mediated presynaptic inhibition in
+the cerebellum. Trends Pharmacol Sci, 30, 421-430. <br>
+.Gillespie, D.T. (1977) Exact stochastic simulation of coupled chemical reactions.
+The Journal of Physical Chemistry, 81, 2340-2361. <br>
      </TR>
 </TABLE>

Difference between revisions of "Team:NJU-China/signaling"

Latest revision as of 03:27, 3 October 2015

3 Signaling module

3.1 Introduction

3.2 Model methods

3.2.1 Modeling the activation of MOR

3.2.2 Modeling adenylate cyclase inhibtion

3.2.3 Modeling GABA vesicle releases

3.2.4 Gillespie’s algorithm

3.3 Results

3.4 Conclusion and remarks

3.5 Model equations, variables and parameters

3.5.1 Activation of MOR

3.5.2 Activation of MOR

3.5.3 GABA release