Skip to main content

Serotonin synthesis, release and reuptake in terminals: a mathematical model



Serotonin is a neurotransmitter that has been linked to a wide variety of behaviors including feeding and body-weight regulation, social hierarchies, aggression and suicidality, obsessive compulsive disorder, alcoholism, anxiety, and affective disorders. Full understanding of serotonergic systems in the central nervous system involves genomics, neurochemistry, electrophysiology, and behavior. Though associations have been found between functions at these different levels, in most cases the causal mechanisms are unknown. The scientific issues are daunting but important for human health because of the use of selective serotonin reuptake inhibitors and other pharmacological agents to treat disorders in the serotonergic signaling system.


We construct a mathematical model of serotonin synthesis, release, and reuptake in a single serotonergic neuron terminal. The model includes the effects of autoreceptors, the transport of tryptophan into the terminal, and the metabolism of serotonin, as well as the dependence of release on the firing rate. The model is based on real physiology determined experimentally and is compared to experimental data.


We compare the variations in serotonin and dopamine synthesis due to meals and find that dopamine synthesis is insensitive to the availability of tyrosine but serotonin synthesis is sensitive to the availability of tryptophan. We conduct in silico experiments on the clearance of extracellular serotonin, normally and in the presence of fluoxetine, and compare to experimental data. We study the effects of various polymorphisms in the genes for the serotonin transporter and for tryptophan hydroxylase on synthesis, release, and reuptake. We find that, because of the homeostatic feedback mechanisms of the autoreceptors, the polymorphisms have smaller effects than one expects. We compute the expected steady concentrations of serotonin transporter knockout mice and compare to experimental data. Finally, we study how the properties of the the serotonin transporter and the autoreceptors give rise to the time courses of extracellular serotonin in various projection regions after a dose of fluoxetine.


Serotonergic systems must respond robustly to important biological signals, while at the same time maintaining homeostasis in the face of normal biological fluctuations in inputs, expression levels, and firing rates. This is accomplished through the cooperative effect of many different homeostatic mechanisms including special properties of the serotonin transporters and the serotonin autoreceptors. Many difficult questions remain in order to fully understand how serotonin biochemistry affects serotonin electrophysiology and vice versa, and how both are changed in the presence of selective serotonin reuptake inhibitors. Mathematical models are useful tools for investigating some of these questions.


Traditionally, serotonin (5-HT) has been associated to a wide variety of behaviors including feeding and body-weight regulation, social hierarchies, aggression and suicidality, obsessive compulsive disorder, alcoholism, anxiety, and affective disorders[1]. In addition, 5-HT has been linked to motor system function[2], sleep-wake cycles[3], circadian rhythms[4], respiratory stability[5], embryonic development[6], and reward processing[7]. Not surprisingly, the 5-HT neurons in the nuclei originally classified by Dalhstrom and Fuxe[8] project to a large variety of regions of the central nervous system including spinal cord, cerebellum, frontal cortex, hypothalamus, hippocampus, striatum, and a bewildering variety of 5-HT receptors have been identified [9]. A huge body of research on genomics, anatomy, neurochemistry, electrophysiology, and behavior has provided a wealth of information on serotonergic systems, but the causal mechanisms of serotonergic function, both normal and in the presence of various disorders and pharmacological agents, remain largely unknown.

Polymorphisms in the serotonin reuptake transporter (SERT) gene have been associated with depression and other mood disorders[1013] and may be associated with anxiety[14], autism[15], and suicidality[16, 17]. Polymorphisms in the tryptophan hydroxylase gene have been associated with unipolar[18] and bipolar disorder[19]. Furthermore, variations in gene expression very likely play a role in the regulation of serotonergic systems both normally and in response to selective serotonin reuptake in-hibitors(SSRIs). SERTs are downregulated in the presence of SSRIs [20, 21], 5-HT1A autoreceptor expression levels differ in different brain regions[22], and 5-HT1A mRNA levels are affected by gonadal hormones [23].

Because of efforts to understand the modes of action of SSRIs, the neurochemistry of serotonin has received much attention. Serotonin is synthesized in serotonergic terminals from tryptophan, which competes with tyrosine and the branched chain amino acids for transport across the blood-brain barrier[1, 24]. Autoreceptors play important roles in the regulation of 5-HT chemistry. For example, 5-HT1B autoreceptors on terminals decrease synthesis and release when extracellular 5-HT rises and 5-HT1A autoreceptors affect firing rates in the dorsal raphe nucleus[25]. In addition, these regulatory mechanisms are themselves regulated by dynamic changes in autoreceptor expression levels[26]. Serotonin acts both in one-to-one neural signaling and as a neuromodulator, via volume transmission, of the effects of other neurotransmitters[1, 27]. Each of these facts plays an important role in neuropsychiatry and neuropharmacology.

The electrophysiology of serotonergic signaling is related both to neurochemistry and to behavior. The classical experiments of Jacobs on cats[28] showed that the patterns of firing of nucleus centralis superior serotonergic neurons correspond to different sleep-wake states. 5-HT modulates motor firing patterns[2] and motor behavior[29, 30]. Autoreceptors affect the inhibition of CA3 hippocampal pyramidal neurons caused by stimulating the ascending serotonergic pathways[31, 32]. 5-HT increases the firing rates of histaminergic neurons in the hypothalamic tuberomammillary nucleus[33], inhibits the firing of somatosensory cortical neurons[34], and can inhibit or excite neurons in the ventromedial nucleus of the hypothalamus[35]. It has been proposed that 5-HT activates the hypothalamic-pituitary-adrenal axis by stimulating production of corticotropin-releasing hormone[36]. 5-HT influences dopaminergic signaling[37, 38] and may affect firing in the cerebral cortex by causing the release of glutamate[39]. Traditionally, dopamine was thought to be the primary neurotransmitter involved in reward processing, but recent work suggests a strong role for 5-HT[7]. Thus, the neurochemistry and electrophysiology affect each other, both affect behavior, and both are affected, of course, by neuronal morphology, which is itself changeable.

Even this brief discussion shows why understanding the casual mechanisms in serotonergic signaling is a challenging problem. Not only does one have to understand mechanism and function on four different levels, genomic, biochemical, electrophysiological, and behavioral, but changes on each level affect function on the other three levels, and this makes the interpretation of experimental and clinical results very difficult. In addition, the brain is not fixed, but dynamical changes on different time scales are happening at all four levels. Mathematical models can play an important role because they allow one to study explicitly the simultaneous effects of all the interactions in a large complex system. Ideas and hypotheses can then be tested by in silico experimentation, that is, by computer simulations of the mathematical model. Our main interest is to understand how the biochemistry of 5-HT (synthesis, release, reuptake) is regulated and how the biochemistry affects the electrophysiology and vice versa. As a first step, we present in this paper a model of 5-HT biochemistry in serotonergic terminals.

The model includes (see Figure 1): uptake of tryptophan across the blood-brain barrier and transport into terminals; synthesis of 5-HT by tryptophan hydroxylase (THP) and aromatic amino acid decarboxylase (AADC); transport of 5-HT into a vesicular compartment by the monamine transporter (MAT); release of 5-HT into the extracellular space depending on firing rate; reuptake via the SERTs; regulation by the autoreceptors. As much as possible, the model is based on real physiology that has been determined experimentally. It is worthwhile to say at the outset that there is no such thing as "the serotonergic terminal"; important parameters (like SERT and autoreceptor densities) vary in different projection regions and this variation is likely to be related to function. Our main purpose is to use the model as a platform for in silico experimentation that sheds light on the complex regulatory mechanisms of serotonergic signaling. Some results of some simulations with the model have previously appeared elsewhere [40].

Figure 1

Steady state concentrations and fluxes. The figure shows the reactions in the model. The rectangular boxes indicate substrates and blue ellipses contain the acronyms of enzymes, transporters, and autoreceptors; steady state values in the model are indicated. Full names of the substrates are given in Table 1. Names of enzymes and transporters are as follows: Trpin, neutral amino acid transporter; DRR, dihydrobiopterin reductase; TPH, tryptophan hydroxylase; AADC, aromatic amino acid decarboxylase; MAT, vesicular monoamine transporter; SERT, 5-HT reuptake transporter; auto, 5-HT autoreceptors; MAO monoamine oxidase; ALDH, aldehyde dehydrogenase. Removal means uptake by capillaries or glial cells or diffusion out of the system.

Mathematical methods have been used by a variety of authors to understand serotonergic signaling. The serotonergic model presented in this paper is conceptually similar to the dopaminergic model presented in [41]; both models were inspired by the original model of Justice et al.[42] for a dopaminergic terminal. Many studies use statistical methods to identify associations between variables on different levels of the serotonergic system. Cohen and colleagues used theoretical and experimental methods to show how 5-HT modulates the frequency and phase lag of bursting in lamprey spinal cord[2, 43]. Butera showed by modeling how 5-HT affects the bursting behavior of neuron R15 in Aplysia[44]. Waggoner and colleagues introduced a three state stochastic model for the serotonin dependence of egg laying in a nematode[45]. Bunin et al.[46] and Daws et al.[47] used mathematical models and data to compute apparent values of the Michaelis-Menten constants K m and V max for the SERTs in different projection regions. Venton et al.[48] used experiments and mathematical models to show that the extracellular space is well-mixed during tonic firing but not during burst firing. Kim et al.[4] used a mathematical model to explain why the rhythmic degradation of the mRNA of serotonin N-acetyltransferase is essential for its circadian rhythm. Tanaka et al.[49] used a mathematical model to show that 5-HT controls the time scale of reward prediction by differentially regulating activities in the striatum. Dayan and Huys[50] used a Markov model to study the effects of 5-HT on how the predictions of future outcomes lead to behavioral inhibition, suppression, and withdrawal and created a computational model to investigate 5-HT in affective control[51]. Stoltenberg and Nag[52] used a dynamical systems model to go directly from genes to behavior.


The mathematical model consists of nine differential equations for the variables listed in Table 1. The differential equations corresponding to the reactions diagrammed in Figure 1 follow in Table 2. Reaction velocities or transport velocities begin with a capital V followed by the name of the enzyme, the transporter, or the process as a subscript. For example, VTPH(trp, bh 4, e 5ht) is the velocity of the tryptophan hydroxylase reaction and it depends on the concentrations of its substrates, trp and bh 4, as well as e 5ht (via the autoreceptors). Below we discuss in detail the more difficult modeling issues and reactions with non-standard kinetics. Table 3 gives the parameter choices and references for reactions that have Michaelis-Menten kinetics in any of the following standard forms:

Table 1 Names used for Variables
Table 2 The Differential Equations
Table 3 Kinetic Parameters (μ M, μ M/hr,/hr).
V = V m a x [ S ] K m + [ S ]
V = V m a x [ S 1 ] [ S 2 ] ( K S 1 + [ S 1 ] ) ( K S 2 + [ S 2 ] )
V = V m a x f [ S 1 ] [ S 2 ] ( K S 1 + [ S 1 ] ) ( K S 2 + [ S 2 ] ) V m a x b [ P 1 ] [ P 2 ] ( K P 1 + [ P 1 ] ) ( K P 2 + [ P 2 ] )

for unidirectional, one substrate, unidirectional, two substrates, and bidirectional, two substrates, two products, respectively.

Table 1 gives the abbreviations used for the variables throughout. We use lower case italic abbreviations in the differential equations and other formulas so that they are easier to read. Full names for the enzymes appear in the legend to Figure 1.

Tryptophan and the tryptophan pool

Serum tryptophan concentrations have been measured in humans and other mammals both before and after meals with different protein composition. A range of 53-85 μ M was found in [53] and a range of 61-173 μ M was found in [54]. We take as our baseline the value of 96 μ M found by Fernstrom in fasted rats [24]. During the experiments with our model in Results A, the serum values of tryptophan were varied corresponding to meals.

Tryptophan is transported across the blood-brain barrier by the L-transporter and is then taken up by serotonergic neuron terminals [55]. We simplify these two steps into a single step with the kinetics of the L-transporter. Choosing the right K m for the L-transporter is complicated by two issues. First, the majority of tryptophan in the serum is not free but bound to albumin. Second, the other neutral and branched chain amino acids compete for the same transporter, so the effective K m depends on the concentrations of these other amino acids. Partridge [56] measured a K m = 190 μ M with respect to total serum tryptophan and Smith [57] measured K m = 15 μ M with respect to free serum tryptophan. We will use the effective K m = 330 μ M in the presence of other amino acids given in Kilberg [55]. We choose V max = 700 μ M/hr so that, in our model, the rate of transport into the brain (159 μ M/hr) closely matches that found by Kilberg (159 μ M/hr).

Intracellular tryptophan is used in a large number of biochemical pathways and, of course, in protein synthesis, which accounts for about half the use of tryptophan [58]. Protein breakdown and a variety of biochemical pathways are intracellular sources of tryptophan. Overall, about 2% of ingested tryptophan is used for the synthesis of serotonin [59, 60]. These numbers give some crude upper and lower bounds for the percentage of intracellular tryptophan that goes to the synthesis of serotonin, but accurate estimates are not known. In dopaminergic neurons about 90% of tyrosine goes to protein synthesis and other pathways and about 10% to dopamine synthesis [6163], so it seems reasonable to make a similar estimate for tryptophan. We let the variable trp-pool represent all the other intracellular sinks and sources of tryptophan and assume that intracellular tryptophan, trp, and trp-pool can be interconverted into each other:

t r p k 1 k 1 t r p p o o l .

We choose the rate constants k1 = 6 μ M/hr and k-1 = .6 μ M/hr so that trp-pool is approximately 10 times as large as trp:

Tryptophan hydroxylase

Tryptophan (trp) and tetrahydrobiopterin (bh 4) are converted by tryptophan hydroxylase (TPH) into 5-hydroxytryptamine (5htp) and dihyro-biopterin (bh2). The velocity of the reaction, VTPH, depends on trp, bh 4, and extracellular 5-HT (e 5ht) via the autoreceptors. We take the basic kinetics from [64] with K trp = 40 μ M, Kbh 4= 20 μ M. TPH exhibits substrate inhibition but it is quite weak, K i = 1000. The second term in the velocity equation below, which represents the effect of extracellular 5-HT on synthesis rate, is discussed in detail below under "autoreceptors." The constants are chosen so that at the normal steady state (e 5ht = .000768 μ M) this factor is equal to one, so the normal steady state is the same with and without the autoreceptors. This allows us to compare how the system changes with and without the autoreceptors when we perturb the system by changing enzyme properties, neuron firing rates, or transporter properties.

V T P H = V m a x ( t r p ) ( b h 4 ) ( K t r p + ( t r p ) + ( t r p ) 2 K i ) ( K b h 4 + ( b h 4 ) ) ( 1.5 ( e 5 h t ) 2 ( ( .000768 ) 2 + ( e 5 h t ) 2 ) )

Storage, release, and reuptake of serotonin

The 5-HTP produced by the TPH reaction is rapidly decarboxylated by the aromatic amino acid decarboxylase (AADC) to produce cytosolic serotonin. We take the parameters of AADC from the literature; see Table 3. The monoamine transporter, MAT, rapidly transports 5-HT into vesicles. We take the K m of the transporter to be 0.198 μ M as found in [65], which is consistent with the values in [66]. We choose the V max so that the concentration of cytosolic serotonin is very low. The experiments in [67] and the calculations in [68] in the case of dopamine suggest strongly that there is transport from the vesicles back into the cytosol, either dependent or independent of the MAT and it is likely that the same is true of serotonin [69]. We assume this transport is linear with rate constant, k out , chosen so that the vast majority (i.e., 98%) of the cellular serotonin is in the vesicular compartment. For simplicity we are assuming that the vesicular compartment is the same size as the non-vesicular cytosolic compartment. This assumption is unimportant since we take the cytosol to be well-mixed and we are not investigating vesicle creation, movement toward the synapic cleft, and recyling where geometry and volume considerations would be crucial. Of course, if we took the volume of the vesicular compartment to be much smaller than the volume of the cytosolic compartment, say 1 to 100, then the ratio of vesicular 5-HT concentration to cytosolic 5-HT concentration would approach the value of 104 suggested in [69].

In our model, vesicular 5-HT (v 5ht in the equations) is removed from the vesicles and put into the synaptic cleft, where it becomes e 5ht, by the term release(e 5ht) fire(t) vda(t) in the differential equations for v 5ht and e 5ht (see the differential equations above). fire is a function of time in some of our in silico experiments, for example in Results B and C where we investigate pulse experiments and in Results E where we consider the effects SSRIs. However, for determining our baseline steady state we take fire = 1 μ M/hr, which means that vesicular serotonin is released at a constant rate such that the entire pool turns over once per hour. The term release(e 5ht) represents the effect of e 5ht on release via the autoreceptors and is discussed below. The processes by which vesicles are created, move to the synapse, and release their serotonin are complicated and interesting [67, 7072], but are not included in this model.

Extracellular serotonin has three fates. It is pumped back into the cytosol by the SERTs; it is catabolized; it is removed from the system. The K m = .17 μ M for the SERTs is taken from [46]. As we will discuss later, the V max will vary considerably from one projection region to another because the density of SERTs varies by at least a factor of 5. For our baseline case, we take V max = 4700 μ M/hr which is in the middle of the range, 2052-6480 μ M/hr, found in [46]. The function fluox(t) that multiplies the term VSERT in the differential equations for the variables v 5ht and e 5ht is the fraction of SERTs that remain unblocked in the presence of an SSRI. In the absence of SSRIs, fluox(t) = 1. Catabolism and removal are discussed below.


It has been understood since the 1970s and 1980s that terminal autoreceptors (5-HT1B) sense the extracellular 5-HT concentration (e 5ht in the equations). When e 5ht goes up, they inhibit both the synthesis of 5-HT and the release of 5-HT from the vesicles into the synaptic cleft and when e 5HT goes down they facilitate synthesis and release [9, 25, 73]. Thus e 5ht provides a kind of end-point feedback for the entire serotonergic system from tryptophan in the serum to e 5HT in the extracellular space. It is also known [74] that autoreceptors modulate reuptake, but this effect is not included in the model. Extracellular 5-HT or autoreceptor agonists can decrease synthesis by 50% [7577] and by perhaps as much as 80-90% [78]. And, autoreceptor antagonists can increase synthesis by as much as 40-60% [77, 79]. These and many other experiments are often conducted with large amounts of agonists or antagonists, which leaves open the question of what range of extracellular 5-HT causes these effects. Experiments on rats [76, 80] showed that cocaine administration elevates extracellular 5-HT by factors of 2 to 5 and that such elevation has a large depressive effect on 5-HT synthesis, so it is reasonable to assume that synthesis is significantly affected by changes in e 5ht over less than an order of magnitude. The second term in the formula for VTPH above contains the effect of e 5ht on synthesis. When extracellular 5-HT has its steady state value of 0.768 nM the factor is equal to 1. As extracellular 5-HT declines towards 0, the factor increases to 1.5 and as extracellular 5-HT increases the factor declines to 0.5 (almost reaching that level when e 5ht = 3 nM). Thus facilitation of synthesis can be as much as 50% and inhibition of synthesis can be as much as 50% and most of the effect is between 0-3 nM of extracellular 5-HT.

Similarly, many experiments have shown that release of vesicular serotonin can be inhibited by increased e 5ht via the autoreceptors or facilitated if e 5ht goes down. For example, Gothert found that release can be inhibited 65% and facilitated by 50-60% [77]. It is not certain from the experiments over what range of e 5ht this effect takes place. We will assume a modest effect over a relatively small range. The factor release(e 5ht) descends linearly from 1.5 at e 5ht = 0 to 1.0 at e 5ht = .000768 μ M, the normal steady state. Then the factor descends linearly from 1.0 at e 5ht = .000768 μ M to 0.4 at e 5ht = .0023 μ M. For e 5ht > .0023, release(e 5ht) remains constant at 0.4. Thus, the maximal facilitation is 50% and the maximal inhibition is 60% and the effect takes place over the range 0 - 2.3 nM of extracellular 5-HT.

Metabolism and removal of serotonin

Serotonin is metabolized by monoamine oxidase (MAO) and aldehyde dehydrogenase (ALDH) to 5-hydroxyindoleacetic acid (5 - hiaa). In our simple model we are not investigating the details of catabolism, only in how c 5ht and e 5ht are removed from the system, so we combine these two steps into one and use the K m = 95 μ M determined in [81, 82]. The rate constant for the removal of 5hiaa was measured to be 0.82 ± .06 in [83]; we take it to be 1/hr. This results in a model steady state concentration of 5hiaa = 5.22 μ M. The ratio of 5-HIAA to 5-HT was measured to be around 1 in [84] and in the range 1-3 in [85]. Since tissue content of 5-HT in different brain regions is roughly 2-3 μ M [8688], the concentration 5hiaa = 5.22 μ M is reasonable.

In our model the extracellular space is a single compartment. One should think of it as the part of the entire extracellular space corresponding to this particular synapse. Of course, if we had many model synapses, the e 5ht from one will diffuse into the extracellular compartment of another (volume transmission). We are assuming for simplicity that the extracellular space is well-mixed, that is, we are ignoring diffusion gradients between different parts of the extracellular space. Venton et al.[48] have shown in the case of dopamine, using a combination of experiments and modeling, that the extracellular space is well-mixed during tonic firing but that substantial gradients exists between "hot spots" of release and reuptake and the rest of the extracellular space during and just after episodes of burst firing. The term k rem (e 5ht) in the differential equation for e 5ht represents removal of e 5ht though uptake by glial cells, uptake by the blood, and diffusion out of the tissue. This will vary from tissue to tissue; for the base case, we chose k rem = 400/hr. At steady state in the base case (see Figure 1), 21.4 μ M/hr of v 5ht is put into the extracellular space, 21.1 μ M/hr is put back into the cytosol by the SERTs, 0.1 μ M/hr is catabolized in the extracellular space, and 0.3 μ M/hr is removed. Thus the effect of the removal term is small at steady state but it plays a bigger role when large amounts of v 5ht are dumped into the extracellular space, for example in the pulse experiments described in Results B.

Figure 2

The effect of meals on brain DA and 5-HT. Panel A shows the blood concentration over a 24 hour period of either tyrosine or tryptophan due to three meals at 7 am, 12 pm and 6 pm. Panel B shows the tyrosine and tryptophan concentrations in dopaminergic and serotonergic synaptic terminals. Panel C shows the velocities of the TH and TPH reactions over the same 24 hour period. Panel D shows the extracellular DA and 5-HT concentrations. The vesicular stores of DA and 5-HT (not shown) vary like the extracellular concentrations in Panel D. All calculations for DA were done using the mathematical model described in [41] and the calculations for 5-HT were done using the mathematical model in this paper.

Fluoxetine dosing

In part E of the results we use the model to investigate the results of giving a dose of the SSRI fluoxetine. The dose is represented in the model by changing the fraction of SERTs that are unblocked at any given time. The resulting function, fluox(t), multiplies VSERT in the differential equations for c 5ht and e 5ht. Since we give the dose at 1 hour, fluox(t) = 1 if t ≤ 1, and for t ≥ 1,

f l u o x ( t ) = 1 ( .95 ) ( t 1 ) 2 .04 + ( t 1 ) 2 e ( t 1 ) / 37 .

The half-life of fluoxetine is quite long; 1-4 days is reported in [89]. The number 37 appearing in the exponential corresponds to a half-life of fluoxetine of a little more than a day.

Steady state concentrations and velocities

Figure 1 shows the concentrations of all of the variables and the reaction and transport velocities at steady state.

The rate of tryptophan uptake from the serum is 159 μ M as found by [55]. The cellular tryptophan concentration, 20.6 μ M, is in the range found in most studies [86, 90]. The rate of the TPH and AADC reactions, which must be equal at steady state, is in the middle of the range of values for 5-HT synthesis reported in the literature, 2.4 μ M/hr [91] to approximately 13 μ M/hr [83, 92].

It is known that the cytosolic concentrations of 5-HTP and 5-HT are quite low. Fernstrom [86] measured 2 μM for 5-HTP and it is 2.26 μ M in the model. It is very diffcult to get reliable estimates for cytosolic 5-HT since tissue measurements include the vesicles where the concentration is known to be very high [69]. In our model cytosolic 5-HT = 0.5 μ M at steady state, a very low value since cytosolic 5-HT is rapidly pumped into the vesicular compartment by the monoamine transporter where the concentration is 21.45 μ M. The vesicular compartment turns over once an hour at steady state releasing 21.45 μ M/hr into the extracellular space. There, 21.13 μ M/hr is put back into the cytosol by the SERTs, 0.3 μ M/hr is removed from the system by uptake by glial cells or blood vessels or simply diffusion out of the tissue, and 0.008 μ M/hr is catabolized.

Most catabolism of 5-HT happens in the cytosol since the extracellular concentrations of 5-HT are so low. Of course, the cytosolic and extracellular catabolism rates plus the rate of removal must add up to the synthesis rate 5.6 μ M/hr. As explained above, the concentration of 5hiaa = 5.22 μ M is reasonable [84, 85].

In all cases, steady states or curves showing the variables as functions of time were computed using the stiff ODE solver in MATLAB.


A. The effect of meals on dopamine and serotonin

Since the early work of Fernstrom [93, 94] it has been generally thought that dopamine synthesis is not very sensitive to tyrosine availability but that serotonin synthesis is sensitive to the availability of tryptophan [1]. We have previously constructed a model of dopamine (DA) synthesis, release, and reuptake in dopaminergic terminals [41], so we can compare and contrast the effects of meals on synthesis in dopaminergic and serotonergic terminals. The model results are displayed in Figure 2, where Panel A shows the blood amino acid concentrations, Panel B the cellular DA and 5-HT concentrations, Panel C the rates of the TH and TPH reactions, and Panel D the concentrations of DA and 5-HT in the extracellular space, during a 24 time period with three meals. In each case the y-axis indicates percent of normal. The overall conclusion is clear: there is much more variation in the 5-HT than the DA concentrations and synthesis rates. To explain why the curves look the way they do, we shall discuss each panel in turn.

It is known that blood amino acid concentrations vary dramatically depending on meal content [5355] and also on the sequence of meals with different content [24]. During a 24 hour period the plasma amino acid concentration can vary as much as a factor of 6 but more typically varies by a factor of 2 to 4 [54]. The amino acid curves in the blood in Panel A of Figure 2 were produced by a simple model that assumes three hours of input corresponding to each meal and a relaxation time back to normal of about 6 hours after the beginning of each meal [53]. For the purpose of these model experiments we assume that the amino acid in the blood is either tyrosine or tryptophan.

Panel B shows the intracellular tyrosine and tryptophan concentrations in the dopaminergic and serotonergic terminals. These large swings in substrate availability correspond to what is seen experimentally; for example, Fernstrom found [95] that brain tyrosine can double after a meal. But why are the oscillations of tryptophan larger than the oscillations of tyrosine? In our model, the tyrosine input into the DA terminal is 241 μ M/hr and the tryptophan input into the 5-HT terminal is 159 μ M/hr corresponding to the experiments reported in [55]. However, the steady concentrations of tyrosine and tryptophan are 126 μ M and 21 μ M, respectively. Thus the tryptophan concentration is much smaller and has a much larger input and removal rate relative to its concentration than does tyrosine. This is why the percentage change due to meals is much larger in the case of tryptophan. This also explains why the tyrosine peaks increase from meal to meal, while the tryptophan peaks are all the same height because tryptophan returns almost to baseline before the next meal.

Panel C shows the velocities of the TH and TPH reactions during the 24 hour period. Despite large swings in tyrosine availability, the TH velocity remains almost constant over the 24 hour period because the reaction is running at near saturation due to the fact that the normal tyrosine concentration is well above the K m for tyrosine, which is 46 μ M in our DA model. In fact, as is clear from the graphs, the rate of the TH reaction actually goes down as tyrosine rises beyond about 100 μ M. This is because TH shows substrate inhibition [96]. In contrast, the oscillations in the TPH curve are large because the normal tryptophan concentration (21 μ M in the model) lies well below the K m of TPH for tryptophan (40 μ M in the model). Thus the rate of the TPH reaction is very sensitive to tryptophan availability. TPH also shows substrate inhibition but it is quite weak and only has an effect at very large (perhaps unphysiological) tryptophan concentrations. We have discussed the functional significance of substrate inhibition elsewhere [40, 97].

Panel D shows the extracellular concentrations of DA and 5-HT over the 24 hour period. The DA concentration varies very little while the 5-HT concentration varies by about 10%. The vesicular stores of DA and 5-HT (not shown) vary similarly in the terminals; since most DA and 5-HT is in these stores, this is what one would see if one measured brain DA or brain 5-HT. Note that we assume in these model calculations that both the dopaminergic and the serotonergic neuron are firing at their tonic rates. Thus these curves give the background concentrations due to tonic firing; burst firing for short periods of time will give significant temporary deviations. Since 5-HT is thought to be an appetite suppressant [1], it makes sense that extracellular 5-HT should rise during and after meals.

B. Release and Reuptake

A number of studies have examined the release of serotonin and reuptake from the extracellular space. Bunin et al.[46] use cyclic voltammetry and Daws et al.[47] use high speed chronoamperometry. Bunin et al. studied release and reuptake of 5-HT in rat brain slices after stimulation with electrical pulses at different frequencies. After stimulation for 1/5 of a second at 100 Hz, the extracellular 5-HT concentration in the dorsal raphe nucleus (DRN) rises to about 1.8 μ M and then declines rapidly back to baseline with half life of about 1 second (their Figure six). In the presence of 10 μ M fluoxetine, the extracellular 5-HT rises slightly higher and declines to baseline with a half-life of 2 seconds. We represent their stimulation in our model neuron by raising fire from 1 to 5000 for 1/5 of a second. The results can be seen in Figure 3. Extracellular 5-HT rises to 2 μ M and then decays back to baseline with a half-life of 1 second (blue curve). We modeled the presence of fluoxetine by blocking half the SERTs. In that case the extracellular 5-HT curve rises slightly higher and decays back to baseline with a half-life of about 2 seconds. Notice that the decay curves do not look exponential. In fact, they are linear until the concentration of extracellular 5-HT gets fairly low. This is because at 1-2 μ M the extracellular 5-HT concentration is well above the K m of the SERTs so the SERTs are saturated and pumping at a constant rate. The same effect can be seen in Figure six of [46]. The decay time back to baseline depends of course on the V max of the SERTs, which in turn depends on the SERT density that is quite different in different brain regions. Daws et al.[47] examined several different brain regions and found much longer half-lives than reported in [46], probably because those regions have much lower SERT densities. This difference may also be due to the fact that 5-HT is applied exogenously in the brains of anesthetized rats in [47] while in [46] tissue slices are stimulated by electrical pulses.

Figure 3

Release and reuptake. The time course of extracellular 5-HT is shown for a model experiment where the neuron was stimulated for 1/5 of a second (blue curve). In the presence of fluoxetine, the time course goes slightly higher and the decay time back to baseline doubles (green curve). We modeled the presence of fluoxetine by blocking half the SERTs. The curves are very similar to those in Figure six of [46].

C. SERT Knockouts

A large number of studies have examined the pharmacological and behavioral characteristics of mice that have the SERT gene knocked out. Such knockouts are of particular interest because they are an (extreme) model of what one could expect with high doses of SSRIs that block the SERTs. Table 4 shows steady state concentration and velocities for WT mice (left column) and steady state concentrations and velocities for SERT knockout mice (right column) in the model. Each column shows the steady state values if a certain fraction, f, of SERTs are functional. Thus, for WT mice, f = 1, and for SERT knockout mice f = 0. The intervening columns corresponds to the effects of progressively higher doses of SSRIs as one moves from left to right.

Table 4 Steady State Values from WT to SERT KO

It is known that 5-HT tissue levels are down both in knockouts [98] and in WT mice treated with fluoxetine [99]. Homberg et al.[98] found that tissue levels of 5-HT drop by 50-70% and Bengel et al.[100] found decreases of 60-80%. In our model, vesicular 5-HT (the main determinant of tissue 5-HT) drops from 21.5 μ M in WT to 6.41 μ M in SERT knockouts corresponding well to these experimental results. Homberg et al. also found that 5-HIAA tissue levels decrease 45-55%; in our model, where we have a very simple model of 5-HIAA metabolism, levels decrease by 88%.

Extracellular 5-HT rises as more and more SERTs are blocked, as expected. Gainetdinov and Caron [12] report that extracellular 5-HT rises 5-6 fold in SERT knockouts and Homberg et al. report a 9-fold increase [98]. In our model, extracellular 5-HT rises to 6.2 μ M from .768 nM, a 9-fold increase. We note that, in the model, extracellular 5-HT already rises substantially when only 50% of the SERTs are active. On the other hand, vesicular 5-HT remains almost normal when only 50% of the SERTs are active, decreasing from 21.5 μ M to 19.9 μ M. This corresponds well with the finding of Bengel et al.[100] that SERT knockout heterozygotes had almost normal tissue levels of 5-HT.

In our model of a 5-HT terminal, we made the SERT knockout terminal by simply setting the V max of VSERT equal to zero. The reality is much more complicated. SERT knockout mice have been that way all their lives and one would expect that other aspects of their serotonergic systems have also changed. SERT knockouts show developmental changes in neurons and brain, an impaired hypothalamic-pituitary-adrenal axis, and desensitization of 5-HT1A and other receptors [101, 102].

D. Homeostatic effects of the autoreceptors

It has been suggested that autoreceptors provide a kind of end product inhibition that tends to stabilize extracellular 5-HT [1, 103]. If extracellular 5-HT goes up, synthesis and release go down; if extracellular 5-HT goes down, synthesis and release go up. Indeed, Panel A of Figure 4 shows that extracellular 5-HT increases and decreases with the tonic firing rate, but the increase and decrease is much less in the presence of the autoreceptors. Thus the autoreceptors help to stabilize extracellular 5-HT in individuals against changes in inputs to the system like changes in firing rate or changes in mean blood tryptophan level (not shown).

Figure 4

Homeostatic effects of the autoreceptors. Panel A shows hows extracellular 5-HT (e 5ht) changes as the firing rate of the neuron varies above and below normal both with and without the autoreceptors. Panel B shows how extracellular 5-HT changes with the expression level of the SERTs both with and without the autoreceptors. s/s and l/l indicate the activities of the corresponding genotypes. Panel C shows how extracellular 5-HT changes with the activity level of TPH both with and without the autoreceptors. The activities of the R441 H and P449R polymorphims are indicated. In all cases, the autoreceptors reduce the effect of changes in firing rate and polymorphisms on extracellular 5-HT.

However, the autoreceptors provide another kind of homeostasis, too. The genes for many of the enzymes and transporters in the serotonergic system have common polymorphisms in different human populations. Many of these polymorphisms are known to be functional in that they change the activity of the corresponding enzymes or the efficacy of the transporters. The autoreceptors tend to keep these serotonergic systems functioning normally despite the polymorphisms.

Polymorphisms in the SERT gene have been associated with depression and other mood disorders [10, 11, 13]. The SERT gene has a polymorphic regulatory region (the 5-HTT gene-linked polymorphic region, 5-HTTLPR), which consists of a variable tandem repeat: the short allele has 14 repeats, whereas the long allele has 16 repeats. The short allele reduces transcriptional activity of the gene and results in decreased expression of the serotonin transporter. The transcriptional activity of the short (s) allele is about 1/3 of that of the long (l) allele [14]. Although the level of transcription of a gene does not necessarily correspond to the activity of its product, we will assume that SERT activity in a s/s homozygote is 1/3 that in the l/l genotype. A study of 505 subjects [14] revealed that in a population sample of 505 individuals, 19% were of the s/s genotype, 49% were l/s, and 32% were l/l. Thus heterozygotes are the most common genotype, and if we assume their SERT activity is 1.0, then the activity of the s/s genotype would be 0.5 and of the l/l genotype 1.5. Panel B of Figure 4 shows that varying SERT activity over this range has a large effect (.5 μ M to 1.6 μ M) on the extracellular 5-HT concentration if the autoreceptors are turned off and a much smaller effect (.6 μ M to 1.1 μ M) if the autoreceptors are turned on.

There are several functional polymorphisms in the TPH2 gene and some are associated with the risk of bipolar disorder [19]. The SNP C2755A changes the amino acid from serine to tyrosine at peptide position 41; the tyrosine coding allele reduces the activity of TPH2 by about 35% [19]. A genetic polymorphism of the promoter, rs11178997, reduces TPH2 transcriptional activity by 22% [104]. The R441 H mutation of TPH2 codes for an enzyme that has only 19% of the wild type activity and the P449R mutation has an activity of 65% of wild type [105]. Thus genetic variation in human populations can cause variation of TPH2 activity between 0.19 and 1.0 of normal. Panel C of Figure 4 shows that varying TPH2 activity over this range has significant effect on the extracellular 5-HT concentration but the effect is less in the presence of the autoreceptors.

As can be seen, the autoreceptors significantly reduce the variation in extracellular serotonin caused by polymorphisms in TPH and SERT.

E. Interaction of autoreceptors and SERTs

Many investigators have studied the effects of doses of SSRIs on extracellular 5-HT in different brain regions. A particular focus of these studies has been to understand the role of the autoreceptors. We have conducted experiments with our model that correspond to some of the experiments in [106109].

Casanovas et al.[108] measured the extracellular 5-HT in the frontal cortex and the hippocampus in the rat after applying doses of 5-HT1A autoreceptor agonists. They found a rapid decline in the frontal cortex to about 30% of basal values and a decline in the hippocampus to about 70% of basal values. These effects are attributed to the stimulation of the 5-HT1A autoreceptors on cells in the raphe since it is known that such stimulation substantially decreases the firing rate of the serotonergic neurons in the raphe that project to the frontal cortex and the hippocampus. The dorsal raphe (DRN) projects to the frontal cortex and the median raphe (MRN) projects to the hippocampus. Casanovas et al. attribute the greater decline in the frontal cortex to the fact the density of 5-HT1A autoreceptors is higher in the DRN than the MRN [110], and therefore firing is reduced much more in the DRN than in the MRN.

In Figure 4, we see that, in the model, a reduction of fire to 58% of normal causes a reduction of extracellular 5-HT at steady state to 70% of normal. Therefore, to simulate the effect of a dose of a 5-HT1A agonist on the extracellular 5-HT in the hippocampus we lowered fire in a time-dependent manner to 58% of normal and then let it recover. Similarly, a reduction of fire to 20% of normal causes a reduction of extracellular 5-HT at steady state to 30% of normal. Therefore, to simulate the effect of a dose of a 5-HT1A agonist on the extracellular 5-HT in the frontal cortex we lowered fire in a time-dependent manner to 20% of normal and then let it recover. The effects on extracellular 5-HT in the frontal cortex and hippocampus can be seen in Panel A of Figure 5. These curves are very similar to those in Figure one (a, b) of Casanovas et al.[108]. The extracellular concentrations of 5-HT decrease in the frontal cortex and hippocampus because the firing rates in the DRN and the MRN are reduced due to the binding of the agonist to the autoreceptors on cell bodies. The concentrations of 5-HT in the frontal cortex and hippocampus begin to recover after the initial decline because the terminal autoreceptors in the frontal cortex and hippocampus increase synthesis and release of 5-HT.

Figure 5

Effects of 5-HT1A agonists and fluoxetine. Panel A shows the change in extracellular 5-HT in the hippocampus and the frontal cortex computed by the model after a 5-HT1A agonist is given. The curves are similar to those in [108]. Panel B shows model computations of the extracellular concentrations of 5-HT in the hippocampus and the frontal cortex after a dose of an SSRI (fluoxetine or paroxetine); the solid curves are wild type and the dashed curves are 5-HT1B knockouts. These curves should be compared to Figure one (10 mg/kg dose) in [106] and Figure one (a, b, c, d) in [109]. For discussion, see the text.

Malagie et al.[106] administered fluoxetine to anaesthetized rats and measured extracellular 5-HT in the frontal cortex and hippocampus. This is a very interesting experiment because fluoxetine blocks SERTs in the DRN and MRN and thus extracellular 5-HT will rise, stimulating the 5-HT1A autoreceptors and decreasing firing as in the experiments Casanovas et al.[107, 108]. This effect will tend to lower extracellular 5-HT in projection regions. However, fluoxetine will also block SERTs in the projection regions, which tends to raise extracellular 5-HT there. Thus, in the projection regions the level of extracellular 5-HT reflects a balance between these two effects. To see what the balance is in our model, we represent a dose of fluoxetine as described under "fluoxetine dosing" in Methods, and assume that fire drops as a function of time in the DRN and MRN as indicated above in the discussion of the experiments of Casanovas et al.. The results are shown by the blue (frontal cortex) and green (hippocampus) solid curves in Figure 5B. Hippocampal extracellular 5-HT rises by approximately 147% and frontal cortex extracellular 5-HT rises approximately 63%. These curves are very similar to the analogous curves in Malagie et al., Figure one (the 10 mg/kg dose), where extracellular 5-HT rises about 110% in the hippocampus and 60% in the frontal cortex. In a second study, Malagie et al.[109] performed similar experiments with paroxetine on mice whose 5-HT1B autoreceptors on the terminals had been knocked out. They saw a large increase in hippocampal extracellular 5-HT and a smaller rise in frontal cortex extracellular 5-HT compared to wild type. In the model we make a 5-HT1B knockout by simply turning off the terminal autoreceptors. The results can be seen in the blue (frontal cortex) and green (hippocampal) dashed curves in Figure 5B. The model hippocampal extracellular 5-HT concentration rises 369% in the knockout, compared to 147% in wild type mice. The model frontal cortex extracellular 5-HT rises 82% as compared to 63% for the wild type. The model curves should be compared to those in Malagie et al. (Figure one,a,b,c,d ). They gave two doses (1 mg/kg and 5 mg/kg). Our model results are simlar to their results but with larger changes than those induced by their 1 mg/kg dose and smaller changes than than those induced by their 5 mg/kg dose. This indicates that our "dose" is between their doses. As they saw, we also found that the knockout induced smaller changes in the frontal cortex than in the hippocampus. Notice that in all cases, there is an increase in extracellular 5-HT caused by the blockage of SERTs. But this effect is then moderated by the increased rate of removal and catabolism, which slowly bring the extracellular concentrations back to equilibria (that are higher than prior to the dose of fluoxetine).


We have presented a relatively simple and straightforward model of synthesis, release, and reuptake of 5-HT in a serotonergic terminal. The kinetics for individual reactions and the values of constants were chosen as much as possible from the experimental literature. The purpose is to create a model that can be used, here and in future investigations, as a platform for exploring various hypotheses about serotonergic homeostasis and serotonergic signaling. Some results and predictions of the model have already appeared in [40]. We note that we have not altered parameters and kinetics to fit any particular set of experimental data. The parameter values remain the same in all the model experiments in the Results sections, except as indicated for changes corresponding to the particular experimental situations that we were examining.

Any model includes many oversimplifications. We have not included the details of the use of tryptophan in other metabolic pathways. The processes by which vesicles are created, move to the synapse, and release their serotonin are complicated and interesting [67, 7072], but are not included in this model. In our model the SERTs put released serotonin back into the terminal, but we do not include leakage of cytosolic serotonin through the SERTs into the extracellular space. We include in the model the effects of the terminal autoreceptors on serotonin synthesis (via TPH) and on serotonin release, but we do not include effects of autoreceptors on reuptake [74]. In this first model we have not included the soma explicitly, so the effects of the somatic autoreceptors are modeled by directly affecting firing rate and thus release in the terminal.

In Section A of Results we use the model to give reasons for the well known observations that dopamine synthesis is relatively insensitive to tyrosine availability, but serotonin synthesis is quite sensitive to tryptophan availability[1, 93, 94]. First, at the normal intracellular concentration of DA, the TH reaction is already running close to saturation, however the normal intracellular concentration of tryptophan is well below the K m of TPH, so changes in availability cause big changes in synthesis rate. Second, the flux into and out of the intracellular pool is much larger(relative to the pool size) in the case of tryptophan than in the case of tyrosine. We showed (Figure 2) the consequences of these differences for the time-dependent behavior of extracellular DA and 5-HT due to meals. In our model calculations, for simplicity, we assumed that the transport of the amino acids tyrosine and tryptophan across the blood brain barrier are independent of each other. In fact, both tyrosine and tryptophan compete for the L-transporter [55] with many other amino acids including the branched chain amino acids (BCAA). The protein composition of meals affects how much tyrosine and tryptophan is imported into the brain and thus how much brain DA and 5-HT change [54, 95, 111, 112]. Even more interesting, Fernstrom[24] has shown that the order of meals affects how much tryptophan gets into the brain. The reason is that carbohydrate meals stimulate insulin production and this tends to transport amino acids into skeletal muscle, but tryptophan is partially protected from this transport because it is bound to serum albumin. A mathematical model for these competitive transport processes is in preparation.

In Section B of results we examine the time courses of extracellular 5-HT after an electric shock both with and without the presence of fluoxetine. Our time courses are very similar to those found in [46] for the DRN and substantia nigra reticulata. The shapes of the curves depend heavily on the density of SERTs, which is known to vary by a factor of 5 in different projection regions [113]. Thus, much slower uptake was found in the dentate gyrus, the corpus callosum and the CA3 region of the hippocampus [47]. This is a good reminder that there is no such thing as a single model of "the" serotonergic terminal. Parameters, both SERT density and also expression levels of 5-HT1A receptors [22], can vary by large amounts in different projection regions, presumably for important functional reasons.

In Section C we examined the steady state concentrations and velocities in the model corresponding to different densities of SERTs, or, equivalently, different doses of SSRIs that block the SERTs. The case where the SERTs are completely blocked corresponds to SERT knockout mice. The model concentrations of extracellular 5-HT and vesicular 5-HT are similar to those found in experiments of SERT knockout mice. Extracellular 5-HT is up 9-fold and vesicular 5-HT is down 70%. Interestingly, as more and more SERTs are blocked corresponding to higher and higher doses of fluoxetine, vesicular 5-HT decreases fairly slowly (Table 4). It is known [99] that tissue levels of 5-HT do decrease in the presence of SSRIs, but this decrease has not been remarked on very much in the literature, perhaps because it is a relatively moderate effect as predicted by our model.

In Section D we examined the steady state effects of the autoreceptors and showed that they produce two kinds of homeostasis. First, they moderate the effects of changes in the cell's environment on the concentration of extracellular 5-HT. We illustrate this by changing the firing rate (Figure 4, Panel A), but similar moderating effects are seen with changes in tryptophan availability or MAO activity. Thus, the autoreceptors allow serotonergic signaling to continue more or less as before in the face of a changed environment. Second, the autoreceptors partially compensate for the effects of various polymorphisms in the genes for TPH and SERTs (Figure 4, Panels B and C). Even though a polymorphism reduces the activity of TPH by 50%, the vesicular and extracellular 5-HT decrease by only 13% (Panel C). Thus, the autoreceptors give a kind of protection against the effects of polymorphisms. We have provided some of the first calculations that show quantitatively the importance of this aspect of autoreceptor function.

Finally, in Section E we conducted model experiments that correspond to experiments in which the time course of extracellular 5-HT was measured in different brain regions of animals after a dose of an SSRI. The purpose of many of these experiments was to investigate functional effects of the somatodendritic (5-HT1A) autoreceptors or the terminal (5-HT1B) autoreceptors, so the experiments were carried both in and without the presence of autoreceptor antagonists. In general, our model calculations correspond reasonably well to the experimental results and give some insight into the reasons why the experimental results look the way they do. We note that autoreceptor densities vary considerably in different brain regions [9, 25, 22], and this variation is likely to have important electrophysiological and behavioral consequences.

Other results of this model have been published previously. It is known that serotonergic neurons in the DRN fire tonically at frequencies of 0.4-2.5 spikes/second and that they also fire short bursts at higher frequencies that convey sensory or motor information [1]. In [40] we showed that the model rersponse to burst firing is very dependent on the density of SERTs on the terminal, which is proportional to the V max of VSERT. The size of this V max determines how long it takes to clear the extracellular space of excess 5-HT after a spike. If this clearance time is approximately the time between spikes during tonic ring, then even a short burst will raise extracellular 5-HT considerably. However, if this clearance time is very short, for example, 1/3 of the time between spikes in tonic firing, then a burst of three spikes at triple the tonic frequency raises extracellular 5-HT little. It is known that the density of SERTs varies by about a factor of five across different projection regions [46, 47, 113]. Interestingly, the frequency of tonic firing of serotonergic cells in the dorsal raphe nucleus also varies by about a factor of four or five [1, 114, 115]. This led us to predict that the SERT density in projection regions is tuned to the tonic firing rate of the DRN cells that project to that region, where "tuned" means that the clearance time is approximately the interspike interval for tonic firing. If it is possible to determine experimentally how the tonic firing rates of DRN cells relate to the region they project to, this prediction can be confirmed or refuted.

The fact that the mathematical model presented here is only for a serotonergic terminal limits our ability to address important issues involving mechanisms at the soma of serotonergic cells and their influence on extracellular 5-HT at terminals in projection regions. The serotonergic cells in the DRN and MRN release 5-HT from both the soma and dendrites and only 70% of the release is related to firing[25, 116, 117]. SSRIs block SERTs on these cell bodies as well as on terminals in projection regions, raising extracellular 5-HT in the DRN and MRN and decreasing firing rate via the 5-HT1A autoreceptors[107, 118]. Thus, acute use of SSRIs can have two conflicting consequences in terminal regions: increased extracellular 5-HT because of SERT blockage on the terminal, and decreased extracellular 5-HT because of SERT blockage on the cell bodies. A common hypothesis is that chronic use of SSRIs does not have a therapeutic effect for several weeks because it takes that long for the 5-HT1A autoreceptors on cells bodies to desensitize[25, 31, 119]. While this hypothesis may have merit, additional factors such as SERTs are also likely to be involved. Studies have shown dramatic downregulation of SERT mRNA during chronic use of SSRIs[20, 21]. It is also possible that the dynamic time scale of such SERT downregulation may contribute to the delay period. An increase in extracellular DA rapidly recruits more DATs to the terminal membrane, but downregulation of DAT activity and density follows as the increase becomes chronic[120]; similar dynamic regulation of SERTs is possible, both at the soma and at terminals. In order to examine the interplay between 5HT1A receptors and dynamic SERT regulation in the presence of SSRIs, we plan to extend our model to include the cell body of the serotonergic cell.

It is worthwhile to keep in mind how difficult the study of the serotonergic system really is. Though much new information is available that gives associations between genotypes and behaviors, the causal mechanisms are mostly unknown. These casual mechanisms necessarily involve cell biochemistry and morphology and the connections between the biochemistry and morphology and the electrophysiology of neurons and networks of neurons. Even more daunting is the fact that the four levels, gene expression, biochemical, electrophysiological, and behavioral, influence each other, both chronically and dynamically. Experiments are often difficult to interpret because changes at more than one level may be involved. In this situation, mathematical models based on real physiology can contribute to understanding, for they provide a platform for testing hypotheses and investigating how changes, chronic or dynamic, at one level cause changes at the other levels.


Serotonergic systems must respond robustly to important biological signals, while at the same time maintaining homeostasis in the face of normal biological fluctuations in inputs, expression levels, and firing rates. Our mathematical model gives insight into how this homeostasis is accomplished through the cooperative effect of many different homeostatic mechanisms including the special properties of tryptophan hydroxylase, the serotonin reuptake transporters, and the serotonin autoreceptors. The model also shows how the autoreceptors moderate the effects of polymorphisms in the genes for the SERTs and TPH. The model calculations correspond quite well to a variety of experimental data. Thus, the model can be useful for testing hypotheses about the relationships between gene expression, biochemistry, and serotonergic signaling.


  1. 1.

    Feldman R, Meyer J, Quenzer L: Principles of Neuropharmacology. 1997, Sunderland, MA: Sinauer Associates, Inc

    Google Scholar 

  2. 2.

    Harris-Warrick RM, Cohen AH: Serotonin modulates the central pattern generator for locomotion in the isolated lamprey spinal cord. J Exper Biol. 1985, 116: 27-46.

    CAS  Google Scholar 

  3. 3.

    Strecker RE, Thackar MM, Porkka-Heiskanen T, Dauphin LJ, Bjorkum AA, McCarley RW: Behavioral state-related changes of extracellular serotonin concentration in the pedunculopontine tegmental nucleus: a microdialysis study in freely moving animals. Sleep Res Online. 1999, 2: 21-27.

    CAS  PubMed  Google Scholar 

  4. 4.

    Kim TD, Kim JS, Kim JH, Myung J, Chae HD, Woo KC, Jang SK, Koh DH, Kim KT: Rhythmic serotonin N-acetyltransferase mRNA degradation is essential for the maintenance of its circadian oscillation. Mole Cell Biol. 2005, 25: 3232-3246. 10.1128/MCB.25.8.3232-3246.2005.

    CAS  Article  Google Scholar 

  5. 5.

    Carley DW, Radulovacki M: Mirtazapine, a mixed-profile serotonin agonist/antagonist, suppresses sleep apnea in the rat. Am J Respir Crit Care Med. 1999, 160: 1824-1829.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Cote F, Fligny C, Bayard E, Launay JM, Gershon MD, Mallet J, Vodjdani G: Maternal serotonin is crucial for murine embryonic development. PNAS. 2007, 104: 329-334. 10.1073/pnas.0606722104.

    PubMed Central  CAS  PubMed  Article  Google Scholar 

  7. 7.

    Kranz GS, Kasper S, Lanzenberger R: Reward and the serotonergic system. Neuroscience. 2010, 166 (4): 1023-1035. 10.1016/j.neuroscience.2010.01.036.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Dahlstrom A, Fuxe K: Evidence for the existence of monoamine-containing neurons in the central nervous system, I: Demonstration of monoamines in the cell bodies of brainstem neurons. Acta Physiol Scand. 1964, 62 (suppl 232): 1-55.

    Google Scholar 

  9. 9.

    Roth BL, Lopez E, Patel S, Kroeze WK: The multiplicity of serotonin receptors: uselessly diverse molecules or an embarrassment of riches?. The Neuroscientist. 2000, 6: 252-262. 10.1177/107385840000600408.

    CAS  Article  Google Scholar 

  10. 10.

    Heils A, Teufel A, Petri S, Stöber G, Riederer P, Bengel D, Lesch KP: Allelic variation of human serotonin transporter gene expression. J Neurochem. 1996, 66: 2621-2624. 10.1046/j.1471-4159.1996.66062621.x.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Caspi A: Influence of life stress on depression: moderation by a polymorphism in the 5-HTT gene. Science. 2003, 301: 386-389. 10.1126/science.1083968.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Gainetdinov RR, Caron MG: Monoamine transporters: from genes to behavior. Ann Rev Pharma Toxicol. 2003, 43: 261-284. 10.1146/annurev.pharmtox.43.050802.112309.

    CAS  Article  Google Scholar 

  13. 13.

    White KJ, Walline CC, Barker EL: Serotonin transporters: implications for antidepressant drug development. AAPS J. 2005, 7: E421-E433. 10.1208/aapsj070242.

    PubMed Central  CAS  PubMed  Article  Google Scholar 

  14. 14.

    Lesch KP, Bengal D, Heils A, Sabol SZ, Greenberg BD: Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region. Science. 1996, 274: 1527-1531. 10.1126/science.274.5292.1527.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Murphy DL, Lesch KP: Targeting the murine serotonin transporter: insights into human neurobiology. Nature Rev Neurosci. 2008, 9: 85-96. 10.1038/nrn2284.

    CAS  Article  Google Scholar 

  16. 16.

    Arango V, Huang YY, Underwood MD, Mann JJ: Genetics of the serotonergic system in suicidal behavior. J Physchiatr Res. 2003, 37: 375-386. 10.1016/S0022-3956(03)00048-7.

    Article  Google Scholar 

  17. 17.

    Anguelova M, Benkelfat C, Turecki G: A systematic review of association studies investigating genes coding for serotonin receptors and the serotonin transporter: II Suicidal behavior. Mol Psychiatry. 2003, 8 (7): 646-653. 10.1038/

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Zhang X, Gainetdinov RD, Beaulieu JM, Sotnikova TD, Burch LH, Williams RB, Schwartz DA, Krishnan KRR, Caron MG: Loss-of-function mutation in tryptophan hydroxylase-2 identified in unipolar major depression. Neuron. 2005, 45: 11-16. 10.1016/j.neuron.2004.12.014.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Lin YMJ, Chao SC, Chen TM, Lai TJ, Chen JS, Sun HS: Association of functional polymorphisms of the human tryptophan hydroxylase 2 gene with risk for bipolar disorder in Han Chinese. Arch Gen Psychiatry. 2007, 64 (9): 1015-1024. 10.1001/archpsyc.64.9.1015.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Benmansour S, Owens WA, Cecchi M, Morilak D, Frazer A: Serotonin clearance in vivo is altered to a greater extent by antidepressant-induced downregulation of the serotonin transporter than by acute blockade of the transporter. J Neurosci. 2002, 22 (15): 6766-6772.

    CAS  PubMed  Google Scholar 

  21. 21.

    Mizra NR, Nielson EO, Troelsen KB: Serotonin transporter density and anxiolytic-like effects of antidepressants in mice. Prog Neuropsycho Biol Psych. 2007, 31: 858-866. 10.1016/j.pnpbp.2007.01.020.

    Article  CAS  Google Scholar 

  22. 22.

    Stein P, Savli M, Wadsakj W, Mitterhauser M, Spindelegger C, Mien LK, Moser U, Dudczak R, Kletter K, Kasper S, Lanzenberger R: The serotonin-1A receptor distribution in healthy men and women measured by PET and [carbonyl-11C]WAY-100635. Eur J Nucl Med. 2008, 35 (12): 2159-2168. 10.1007/s00259-008-0850-x.

    CAS  Article  Google Scholar 

  23. 23.

    Birzniece V, Johansson IM, Wang MD, Secki JR, Backstrom T, Olsson T: Serotonin 5-HT1A receptor mRNA expression in dorsal hippocampus and raphe nuclei after gonadal hormone manipulation in female rats. Neuroendocrinology. 2001, 74 (2): 135-142. 10.1159/000054679.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Fernstrom MH, Fernstrom JD: Brain tryptophan concentrations and serotonin synthesis remain responsive to food consumption after the ingestion of sequential meals. Am J Clin Nutr. 1995, 61: 312-319.

    CAS  PubMed  Google Scholar 

  25. 25.

    Adell A, Celada P, Abellán MT, Artigas F: Origin and functional role of the extracellular serotonin in the midbrain raphe nuclei. Brain Res Rev. 2002, 39: 154-180. 10.1016/S0165-0173(02)00182-0.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Jansson A, Descarries L, Cornea-Hebert V, Riad M, Verge D, Bancila M, Agnati LF, Fuxe K: Transmitter-Receptor Mismatches in Central Dopamine, Serotonin, and Neuropeptide Systems: Further Evidence for Volume Transmission. The Neuronal Environment: Brain Homeostasis in Health and Disease. Edited by: Walz W. 2002, Humana Press, 83-108.

    Google Scholar 

  27. 27.

    Robinson DR, Hermans A, Seipel AT, Wightman RM: Monitoring rapid chemical communication in the brain. Chem Rev. 2008, 108: 2554-2584. 10.1021/cr068081q.

    PubMed Central  CAS  PubMed  Article  Google Scholar 

  28. 28.

    Rasmussen K, Heym J, Jacobs BL: Activity of serotonin-containing neurons in nucleus centralis superior of freely moving cats. Exper Neurol. 1984, 83: 302-317.

    CAS  Article  Google Scholar 

  29. 29.

    Fletcher PJ, Tampakeras M, Yeomans JS: Median raphe injections of 8-OH-DPAT lower frequency thresholds for lateral hypothalamic self-stimulation. Pharma Biochem Beh. 1995, 52: 65-71. 10.1016/0091-3057(94)00441-K.

    CAS  Article  Google Scholar 

  30. 30.

    Jacobs BL, Fornal CA: 5-HT and motor control: a hypothesis. Trends Neurosci. 1993, 16: 346-352. 10.1016/0166-2236(93)90090-9.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Chaput Y, Blier P, d Montigny C: In vivo electrophysiological evidence for the regulatory role of autoreceptors on serotonergic terminals. J Neurosci. 1986, 6 (10): 2796-2801.

    CAS  PubMed  Google Scholar 

  32. 32.

    Blier P, Chaput Y, Montigny CD: Long-term 5-HT reuptake blockade, but not monoamine oxidase inhibition, decreases the function of terminal 5-HT autoreceptors: an electrophysiological study in the rat brain. Naun-Schmied Arch Pharmacol. 1988, 337: 246-254.

    CAS  Article  Google Scholar 

  33. 33.

    Eriksson KS, Stevens DR, Haas HL: Serotonin excites tuberomammillary neurons by activation of Na+/Ca2+ exchange. Neuropharmacology. 2001, 40: 345-351. 10.1016/S0028-3908(00)00175-1.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Yarbrough GG, Singh DK, Pettibone DJ: A comparative electrophysiological and biochemical assessment of serotonin (5-HT) and a novel 5-HT agonist (MK-212) on central serotonergic receptors. Neuropharmacology. 1984, 23 (11): 1271-1277. 10.1016/0028-3908(84)90044-3.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Kow LM, Tsai YF, Wang L, Pfaff DW: Electrophysiological analyses of serotonergic actions on neurons in hypothalamic ventromedial nucleus in vitro: receptor subtypes involved and implications for regulation of feeding and lordosis behaviors. Chin J Physiol. 1992, 35 (2): 105-121.

    CAS  PubMed  Google Scholar 

  36. 36.

    Heisler LK, pronchuk N, Nonogaki K, Zhou L, Raber J, Tung L, Yeo GSH, O'Rahilly S, Colmers WF, Elmquist JK, Tecott LH: Serotonin activates the hypothalamic-pituitary-adrenal axis via serotonin 2C receptor stimulation. J Neurosci. 2007, 27 (26): 6956-6964. 10.1523/JNEUROSCI.2584-06.2007.

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Mascio MD, Giovanni GD, Matteo VD, Esposito E: Decreased chaos of mid-brain dopaminergic neurons after serotonin denervation. Neuroscience. 1999, 92: 237-243. 10.1016/S0306-4522(98)00757-X.

    PubMed  Article  Google Scholar 

  38. 38.

    georgetti M, Tecott LH: Contributions of 5-HT2C receptors to multiple actions of central serotonin systems. Eur J Pharmacol. 2004, 488: 1-9. 10.1016/j.ejphar.2004.01.036.

    Article  CAS  Google Scholar 

  39. 39.

    Marek GJ, Aghajanian GK: The electrophysiology of prefrontal serotonin systems: therapeutic implications for mood and psychosis. Biolog Psychiatry. 1998, 44 (11): 1118-1127. 10.1016/S0006-3223(98)00036-5.

    CAS  Article  Google Scholar 

  40. 40.

    Best JA, Nijhout HF, Reed MC: Models of dopaminergic and serotonergic signaling. Pharmacopsychiatry. 2010, 43 (Supp 1): 561-566.

    Google Scholar 

  41. 41.

    Best JA, Nijhout HF, Reed MC: Homeostatic mechanisms in dopamine synthesis and release: a mathematical model. Theor Biol Med Model. 2009, 6: 21-10.1186/1742-4682-6-21.

    PubMed Central  PubMed  Article  CAS  Google Scholar 

  42. 42.

    Justice J, Nicolaysen L, Michael A: Modeling the dopaminergic nerve terminal. J Neurosci Meth. 1988, 22: 239-252. 10.1016/0165-0270(88)90045-3.

    Article  Google Scholar 

  43. 43.

    Cohen AH, Holmes PJ, Rand RH: The nature of the coupling between segmental oscillators of the lamprey spinal generator for locomotion:a mathematical model. J Math Biol. 1982, 13: 345-369. 10.1007/BF00276069.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Butera RJ: Implementation and analysis of neuromodulatory mechanisms in a mathematical model of neuron R15 in Aplysia. Master's thesis, Rice University. 1984

    Google Scholar 

  45. 45.

    Waggoner LE, Zhou GT, Schafer RW, Schafer WR: Control of alternative behavioral states by serotonin in Caenorhabditis elegans. Neuron. 1998, 21: 203-214. 10.1016/S0896-6273(00)80527-9.

    CAS  PubMed  Article  Google Scholar 

  46. 46.

    Bunin M, Prioleau C, Mailman R, Wightman R: Release and uptake rates of 5-hydroxytryptamine in the dorsal Raphe and substantia nigra of the rat brain. J Neurochem. 1998, 70: 1077-1087. 10.1046/j.1471-4159.1998.70031077.x.

    CAS  PubMed  Article  Google Scholar 

  47. 47.

    Daws L, Montenez S, Owens W, Gould G, Frazer A, Toney G, Gerhardt G: Transport mechanisms governing serotonin clearance in vivo revealed by high speed chronoamperometry. J Neurosci Meth. 2005, 143: 49-62. 10.1016/j.jneumeth.2004.09.011.

    CAS  Article  Google Scholar 

  48. 48.

    Venton B, Zhang H, Garris P, Sulzer D, Phillips P, Wightman R: Real-time decoding of dopamine neurotransmission in the caudate-putamen during tonic and phasic firing. J Neurochem. 2003, 87: 1284-1295. 10.1046/j.1471-4159.2003.02109.x.

    CAS  PubMed  Article  Google Scholar 

  49. 49.

    Tanaka SC, Schweighofer N, Asahi S, Shishida K, Okamoto Y, Yamawaki S, Doya K: Serotonin differentially regulates short- and long-term prediction of rewards in the ventral and dorsal striatum. PLoS ONE. 2007, 2 (12): 1-7. 10.1371/journal.pone.0001333.

    Article  CAS  Google Scholar 

  50. 50.

    Dayan P, Huys QJM: Serotonin, inhibition, and negative mood. PLoS Comp Biol. 2008, 4 (2): 1-11. 10.1371/journal.pcbi.0040004.

    Article  CAS  Google Scholar 

  51. 51.

    Dayan P, Huys QJM: Serotonin in affective control. Ann Rev Neurosci. 2009

    Google Scholar 

  52. 52.

    Stoltenberg SF, Nag P: Description and validation of a dynamical systems model of presynaptic serotonin function:genetic variation, brain activation and impulsivity. Behavior Genetics. 2010, 40: 262-279. 10.1007/s10519-010-9335-3.

    PubMed Central  PubMed  Article  Google Scholar 

  53. 53.

    Nasset ES, Heald FP, Calloway DH, Margen S, Scheeman P: Amino acids in human blood plasma after single meals of meat, oil, sucrose, and whiskey. J Nutr. 1979, 109: 621-630.

    CAS  PubMed  Google Scholar 

  54. 54.

    Glaeser BS, Maher TJ, Wurtman RJ: Changes in brain levels of acidic, basic, and neutral amino acids after consumption of single meals containing various proportions of protein. J Neurochem. 1983, 41: 1016-1021. 10.1111/j.1471-4159.1983.tb09045.x.

    CAS  PubMed  Article  Google Scholar 

  55. 55.

    Kilberg M, Haussinger D: Mammalian Amino Acid Transport, Mechanisms and Control. 1992, New York: Plenum Pr

    Book  Google Scholar 

  56. 56.

    Partridge WM, Oldendorf WH: Kinetic analysis of blood-brain barrier transport of amino acids. Biochenica Biophysica Acta. 1975, 401: 128-136. 10.1016/0005-2736(75)90347-8.

    Article  Google Scholar 

  57. 57.

    Smith QR, Momma S, Aoyagi M, Rapoport SI: Kinetics of neutral amino acid transport across the blood-brain barrier. J Neurchem. 1987, 49: 1651-1658. 10.1111/j.1471-4159.1987.tb01039.x.

    CAS  Article  Google Scholar 

  58. 58.

    Kalyanasundaramand S, Ramanamurthy PSV: Utilization of tyrosine and tryptophan for protein synthesis by undernourished developing rat brain. Neurochem Res. 1983, 8 (11): 1471-1480. 10.1007/BF00965002.

    Article  Google Scholar 

  59. 59.

    Sirek A, Sirek OV: Serotonin: A Review. Canad Med Ass J. 1970, 102: 846-849.

    PubMed Central  CAS  PubMed  Google Scholar 

  60. 60.

    Brown RR: Tryptophan Metabolism: A Review. L-Tryptophan: Current Prospects in Medicine and Drug Safety. Edited by: Kochen W, Steinhart H. 1994, Berlin: Walter de Gruyter, 17-30.

    Google Scholar 

  61. 61.

    Mains R, Patterson P: Primary cultures of dissociated sympathetic neurons. II. Initial studies on catecholamine metabolism. J Cell Biol. 1973, 59: 346-360. 10.1083/jcb.59.2.346.

    PubMed Central  CAS  PubMed  Article  Google Scholar 

  62. 62.

    Coenen H, Kling P, Stocklin G: Cerebral metabolism of L-[2-18F] fluorotyrosine, a new PET tracer of protein synthesis. J Nucl Med. 1989, 30: 1367-1372.

    CAS  PubMed  Google Scholar 

  63. 63.

    Cumming P, Asa A, Kuwabara H, Gjedde A: [3H]DOPA formed from [3H]tyrosine in living rat brain is not committed to dopamine synthesis. J Cereb Blood Flow Metab. 1998, 18: 491-499. 10.1097/00004647-199805000-00004.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    McKinney J, Knappskog PM, Haavik J: Different properties of the central and peripheral forms of human tryptophan hydroxylase. J Neurochem. 2005, 92: 311-320. 10.1111/j.1471-4159.2004.02850.x.

    CAS  PubMed  Article  Google Scholar 

  65. 65.

    Rau KS, Volz TJ, Riordan JA, BP Adair AJBI, Bitter R, Gibb JW, Hanson GR, Fleckenstein AE: Methamphetamine Administration Reduces Hippocampal Vesicular Monoamine Transporter-2 Uptake. J Pharmacol Exp Ther. 2006, 318: 676-682. 10.1124/jpet.105.099200.

    CAS  PubMed  Article  Google Scholar 

  66. 66.

    Slotkin TA, Seidler FJ, Whitmore WL, Lau C, Salvaggio M, Kirksey DK: Rat brain synaptic vesicles: uptake specificities of [3H]norepinephrine and [3H]serotonin in preparations from whole brain and brain regions. J Neurochem. 1978, 31: 961-968. 10.1111/j.1471-4159.1978.tb00134.x.

    CAS  PubMed  Article  Google Scholar 

  67. 67.

    Floor E, Leventhal P, Wang Y, Meng L, Chen W: Dynamic storage of dopamine in rat brain synaptic vesicles in vitro. J Neurochem. 1995, 64: 689-699. 10.1046/j.1471-4159.1995.64020689.x.

    CAS  PubMed  Article  Google Scholar 

  68. 68.

    Wallace L: A small dopamine permeability of storage vesicles membranes and end product inhibition of tyrosine hydroxylase are sufficient to explain changes occurring in dopamine synthesis and storage after inhibition of neuronal firing. Synapse. 2007, 61: 715-723. 10.1002/syn.20408.

    CAS  PubMed  Article  Google Scholar 

  69. 69.

    Schuldiner S, Shirvan A, Linial M: Vesicular neurotransmitter transporters: from bacteria to humans. Physiological Reviews. 1995, 75 (2): 369-392.

    CAS  PubMed  Google Scholar 

  70. 70.

    Reimer R, Fon E, Edwards R: Vesicular neurotransmitter transport and the presynaptic regulation of quantal size. Curr Opin Neurobiol. 1998, 8: 405-412. 10.1016/S0959-4388(98)80068-8.

    CAS  PubMed  Article  Google Scholar 

  71. 71.

    Kuhn DM, Wolf WA, Youdin MB: 5-Hydroxytryptamine release in vivo from a cytoplasmic pool: studies on the 5-HT behavioural syndrome in reserpinized rats. Br J Pharmacol. 1985, 84: 121-129.

    PubMed Central  CAS  PubMed  Google Scholar 

  72. 72.

    de Iraldi AP: Compartmentalization of monoaminergic synaptic vesicles in the storage and release of neurotransmitter. Mol Neurobiol. 1992, 6 (4): 323-337. 10.1007/BF02757939.

    Article  Google Scholar 

  73. 73.

    Starke K, Gothert M, Kilbinger H: Modulation of neurotransmitter release by presynaptic autoreceptors. Physiol Rev. 1989, 69: 864-989.

    CAS  PubMed  Google Scholar 

  74. 74.

    Daws LC, Gerhardt GA, Frazer A: 5-HT1B antagonists modulate clearance of extracellular serotonin in rat hippocampus. Neurosci Lett. 1999, 266 (3): 165-168. 10.1016/S0304-3940(99)00277-3.

    CAS  PubMed  Article  Google Scholar 

  75. 75.

    Estaban S, Llado J, Sastre-cool A, Garcia-Sevilla JA: Activation and desensitization by cyclic antidepressant drugs of 2-autoreceptors, 2-heteroreceptors and 5-HT1A-autoreceptors regulating monoamine synthesis in the rat brain in vivo. Naunyn-Schmiedeberg's Arch Pharmacol. 1999, 360: 135-143. 10.1007/s002109900045.

    Article  Google Scholar 

  76. 76.

    Galloway M: Regulation of dopamine and serotonin synthesis by acute administration of cocaine. Synapse. 1990, 6: 63-72. 10.1002/syn.890060108.

    CAS  PubMed  Article  Google Scholar 

  77. 77.

    Gothert M: Presynaptic serotonin receptors in the central nervous system. Ann N Y Acad Sci. 1990, 604: 102-112. 10.1111/j.1749-6632.1990.tb31986.x.

    CAS  PubMed  Article  Google Scholar 

  78. 78.

    Hjorth S, Suchowski CS, Galloway MP: Evidence for 5-HT autoreceptor-mediated, nerve impulse-independent, control of 5-HT synthesis in the rat brain. Synapse. 1995, 19: 170-176. 10.1002/syn.890190304.

    CAS  PubMed  Article  Google Scholar 

  79. 79.

    Jacoby JH, Shabshelowitz H, Fernstrom JD, Wurtman RJ: The mechanisms by which methiothepin, a putative serotonin receptor antagonist, increases brain 5-hydroxyindole levels. J Pharmacol Exp Ther. 1975, 195 (2): 257-264.

    CAS  PubMed  Google Scholar 

  80. 80.

    Parsons LH, Weiss F, Koob GF: Serotonin 1B receptor stimulation enhances dopamine-mediated reinforcement. Psycopharmacology. 1996, 128: 150-160. 10.1007/s002130050120.

    CAS  Article  Google Scholar 

  81. 81.

    Gottowik J, Cesura AM, Malherbe P, Lang G, Prada MD: Characterisation of wild-type and mutant forms of human monoamine oxidase A and B expressed in a mammalian cell line. Fed Eur Biochem Soc. 1993, 317 (1-2): 152-156.

    CAS  Article  Google Scholar 

  82. 82.

    Fowler CJ, Ross SB: Selective inhibitors of monoamine oxidase A and B: biochemical, pharmacological, and clinical properties. Med Res Rev. 1984, 4 (3): 323-358. 10.1002/med.2610040303.

    CAS  PubMed  Article  Google Scholar 

  83. 83.

    Echizen H, Freed CR: Measurement of serotonin turnover rate in rat dorsal raphe nucleus by in vivo electrochemistry. J Neurchem. 1984, 42: 1483-1486. 10.1111/j.1471-4159.1984.tb02815.x.

    CAS  Article  Google Scholar 

  84. 84.

    Schweiger U, Broocks A, Tuschel RJ, Pirke KM: Serotonin turnover in rat brain during semistarvation with high-protein and high-carbohydrate diets. J Neural Tranm. 1989, 77: 131-139. 10.1007/BF01248926.

    CAS  Article  Google Scholar 

  85. 85.

    Asikainen M, Deboer T, Porkka-Heiskanen T, Stenberg D, Tobler I: Sleep deprivation increases brain serotonin turnover in the Djungarian hamster. Neurosci Lett. 1995, 198: 21-24. 10.1016/0304-3940(95)11953-T.

    CAS  PubMed  Article  Google Scholar 

  86. 86.

    Fernstrom JD, Wurtman RJ: Brain serotonin content: physiological regulation by plasma neutral amino acids. Science. 1972, 178 (4059): 414-416. 10.1126/science.178.4059.414.

    CAS  PubMed  Article  Google Scholar 

  87. 87.

    Haubrich DR, Denzer JS: Simultaneous extraction of brain serotonin, and fluorometric measurement of brain serotonin, catecholamines, 5-hydroxyindoleacetic acid and homovanillic acid. Analytical Biochem. 1973, 55: 306-312. 10.1016/0003-2697(73)90320-5.

    CAS  Article  Google Scholar 

  88. 88.

    Nowak P, Bortel A, Dabrowska J, Oswiecimska J, Drosik M, Kwiecinski A, Opara J, Kostrzewa RM, Brus R: Amphetamine and mCPP effects on dopamine and serotonin striatal in vivo microdialysates in an animal model of hyperactivity. Neurotox Res. 2007, 11 (2): 131-144. 10.1007/BF03033391.

    CAS  PubMed  Article  Google Scholar 

  89. 89.

    Hiemke C, Hartter S: Pharmacokinetics of selective serotonin reuptake inhibitors. Pharmaco Therap. 2000, 85: 11-28. 10.1016/S0163-7258(99)00048-0.

    CAS  Article  Google Scholar 

  90. 90.

    Leathwood PD: Tryptophan availability and serotonin synthesis. Proc Nutr Soc. 1987, 46: 143-156. 10.1079/PNS19870018.

    CAS  PubMed  Article  Google Scholar 

  91. 91.

    Diksic M, Hagahiro S, Sourkes TL, Yamamoto YL: A new method to mea-sure brain serotonin synthesis in vivo. I. Theory and basic data for a biological model. J Cereb Blood Flow Metab. 1990, 10: 1-12.

    CAS  PubMed  Article  Google Scholar 

  92. 92.

    Sawada M, Nagatsu T: Stimulation of the serotonin autoreceptor prevents the calcium-calmodulin-dependent increase of serotonin biosynthesis in rat raphe slices. J Neurochem. 1986, 46: 963-967. 10.1111/j.1471-4159.1986.tb13063.x.

    CAS  PubMed  Article  Google Scholar 

  93. 93.

    Fernstrom J, Wurtman RD: Brain serotonin content: physiological dependence on plasma tryptophan levels. Science. 1971, 173: 149-152. 10.1126/science.173.3992.149.

    CAS  PubMed  Article  Google Scholar 

  94. 94.

    Fernstrom J: Role of precursor availability in control of monoamine biosynthesis in brain. Physiol Rev. 1983, 63: 484-546.

    CAS  PubMed  Google Scholar 

  95. 95.

    Fernstrom J, Fernstrom M: Dietary effects on tyrosine availability and cat-echolamine synthesis in the central nervous system: possible relevance to the control of protein intake. Proc Nutr Soc. 1994, 53: 419-429. 10.1079/PNS19940047.

    CAS  PubMed  Article  Google Scholar 

  96. 96.

    Nakashima A, Mori K, Suzuki T, Kurita H, Otani M, Nagatsu T, Ota A: Dopamine inhibition of human tyrosine hydroxylase type I is controlled by the specific portion of the N-terminus of the enzyme. J Neurchem. 1999, 72: 2145-2153. 10.1046/j.1471-4159.1999.0722145.x.

    CAS  Article  Google Scholar 

  97. 97.

    Reed M, Lieb A, Nijhout H: The biological significance of substrate inhibition: a mechanism with diverse functions. BioEssays. 2010, 32: 422-429. 10.1002/bies.200900167.

    CAS  PubMed  Article  Google Scholar 

  98. 98.

    Homberg JR, Olivier JDA, Smits BMG, Mul JD, Mudde J, Verheul M, Nieuwe-huizen OFM, Cools AR, Ronken E, Cremers T, Schoffelmeer ANM, Ellenbroek BA, Cuppen E: Characterization of the serotonin transporter knock out rat: a selective change in the functioning of the serotonergic system. Neurosci. 2007, 146: 1662-1676. 10.1016/j.neuroscience.2007.03.030.

    CAS  Article  Google Scholar 

  99. 99.

    Zusso M, Debetto P, Guidolin D, Barbierato M, Manev H, Giusti P: Fluoxetine-induced proliferation and differentiation of neural progenitor cells isolated from rat postnatal cerebellum. Biochemical Pharmacology. 2008, 76: 391-403. 10.1016/j.bcp.2008.05.014.

    CAS  PubMed  Article  Google Scholar 

  100. 100.

    Bengel D, Murphy DL, Andrews AM, Wichems CH, Feltner D, Heils A, Mossner R, Westphal H, Lesch KP: Altered brain serotonin homeostasis and locomotor insensitivity to 3,4-Methylenedioxymethamphetamine ("Ecstasy") in serotonin transporter-deficient mice. Amer Soc Pharma Exper Therap. 1998, 53: 649-655.

    CAS  Google Scholar 

  101. 101.

    Li Q: Cellular and molecular alterations in mice with deficient and reduced serotonin transporters. Molecular Neurobiology. 2006, 34: 51-65. 10.1385/MN:34:1:51.

    PubMed  Article  Google Scholar 

  102. 102.

    Gainetdinov RR, Sotnikova TD, Caron MG: Monoamine transporter pharmacology in mutant mice. Trends in Pharmacological Sciences. 2002, 23: 367-373. 10.1016/S0165-6147(02)02044-8.

    CAS  PubMed  Article  Google Scholar 

  103. 103.

    Cooper J, Bloom F, Roth R: The Biochemical Basis of Neuropharmacology. 2003, New York, NY: Oxford U Press

    Google Scholar 

  104. 104.

    Scheuch L, Lautenschlager M, Grohmann M, Stahlberg S, Kirchheiner J, Zill P, Heinz A, Walther D, Priller J: Characterization of a functional promoter polymorphism of the human tryptophan hydroxylase 2 gene in serotonergic raphe neurons. Biol Psychiatry. 2007, 62 (11): 1288-1294. 10.1016/j.biopsych.2007.01.015.

    CAS  PubMed  Article  Google Scholar 

  105. 105.

    Winge I, McKinney J, Knappskog P, Haavik J: Characterization of wild-type and mutant forms of human tryptophan hydroxylase 2. J Neurochem. 2007

    Google Scholar 

  106. 106.

    Malagie I, Trillat AC, Jacquot C, Gardier AM: Effects of acute fluoxetine on extracellular serotonin levels in the Raphe: an in vivo microdialysis study. Eur J Pharmacol. 1995, 286: 213-217. 10.1016/0014-2999(95)00573-4.

    CAS  PubMed  Article  Google Scholar 

  107. 107.

    Casanovas JM, Artigas F: Differential effects of ipsapirone on 5-hydroxytryptamine release in the dorsal and median Raphe neuronal pathways. J Neurochem. 1996, 67: 1945-1952. 10.1046/j.1471-4159.1996.67051945.x.

    CAS  PubMed  Article  Google Scholar 

  108. 108.

    Casanovas JM, Lesourd M, Artigas F: The effect of the selective 5-HT1A agonists alnespirone (S-20499) and 8-OH-DPAT on extracellular 5-hydroxytryptamine in different regions of rat brain. Brit J Pharmacol. 1997, 122: 733-741. 10.1038/sj.bjp.0701420.

    CAS  Article  Google Scholar 

  109. 109.

    Malagie I, Trillat AC, Bourin M, Jacquot C, Hen R, Gardier AM: 5-HT1B autoreceptors limit the effects of selective serotonin re-uptake inhibitors in mouse hippocampus and frontal cortex. J Neurochem. 2001, 76: 865-871. 10.1046/j.1471-4159.2001.00083.x.

    CAS  PubMed  Article  Google Scholar 

  110. 110.

    Weissmann-nanopoulos D, Mach E, Magre J, Demassey Y, Pujol JF: Evidence for the localization of 5HT1A binding sites on serotonin containing neurons in the raphe dorsalis and raphe centralis nuclei of the rat brain. Neurochem Int. 1985, 7: 1061-1072. 10.1016/0197-0186(85)90156-1.

    CAS  PubMed  Article  Google Scholar 

  111. 111.

    Fernstrom J: Branched chain amino acids and brain function. J Nutr. 2005

    Google Scholar 

  112. 112.

    Richard DM, Dawes MA, Mathias CW, Acheson A, Hill-Kapturczak N, Dougherty DM: L-tryptophan: basic metabolic functions, behavioral research, and therapeutic indications. Int J Tryp Res. 2009, 2: 45-60.

    CAS  Google Scholar 

  113. 113.

    Lin KJ, Yen TC, Wey SP, Hwang JJ, Ye XX, Tzen KY, Fu YK, Chen JC: Characterization of the binding sites 123I-ADAM and the relationship to the serotonin transporter in rat and mouse brains using quantitative autoradiography. J Nuc Med. 2004, 45: 673-681.

    CAS  Google Scholar 

  114. 114.

    Waterhouse BD, Devilbiss D, Seiple S, Markowitz R: Sensorimotor-related discharge of simultaneously recorded, single neurons in the dorsal raphe nucleus of the awake unrestrained rat. Brain Res. 2004, 1000: 183-191. 10.1016/j.brainres.2003.11.030.

    CAS  PubMed  Article  Google Scholar 

  115. 115.

    Aghajanian GK, Wang RY, Baraban J: Serotonergic and non-serotonergic neurons of the dorsal raphe: reciprocal changes in firing induced by peripheral nerve stimulation. Brain Res. 1978, 153: 169-175. 10.1016/0006-8993(78)91140-X.

    CAS  PubMed  Article  Google Scholar 

  116. 116.

    Becquet D, Faudon M, Hery F: The role of serotonin release and autoreceptors in the dorsal raphe nucleus in the control of serotonin release in the cat caudate nucleus. Neuroscience. 1990, 3: 639-647. 10.1016/0306-4522(90)90248-3.

    Article  Google Scholar 

  117. 117.

    O'Connor JJ, Kruk ZL: Pharmacological characteristics of 5-hydroxytryptamine autoreceptors in rat brain slices incorporating the dorsal raphe or the suprachiasmatic nucleus. Br J Pharmacol. 1992, 106: 524-532.

    PubMed Central  PubMed  Article  Google Scholar 

  118. 118.

    Gartside SE, Umbers V, Hajos M, Sharp T: Interaction between a selective 5-HT1A receptor antagonist and an SSRI in vivo: effects on 5-HT cell firing and extracellular 5-HT. Br J Pharmacol. 1995, 115: 1064-1070.

    PubMed Central  CAS  PubMed  Article  Google Scholar 

  119. 119.

    Stahl S: Mechanism of action of serotonin selective reuptake inhibitors: serotonin receptors and pathways mediate therapeutic effects and side effects. Affective Disorders. 1998, 51: 215-235. 10.1016/S0165-0327(98)00221-3.

    CAS  Article  Google Scholar 

  120. 120.

    Chen R, Furman CA, Gnegy ME: Dopamine transporter tracking: rapid response on demand. Future Neurology. 2010, 5: 123-134. 10.2217/fnl.09.76.

    PubMed Central  CAS  PubMed  Article  Google Scholar 

  121. 121.

    Sumi C, Ichinose H, Nagatsu T: Characterization of recombinant human aromatic L-amino acid decarboxylase expressed in COS cells. J Neurochem. 1990, 55: 1075-1078. 10.1111/j.1471-4159.1990.tb04601.x.

    CAS  PubMed  Article  Google Scholar 

  122. 122.

    Armarego W, Ohnishi A, Taguchi H: New pteridine substrates for dihydropteridine reductase and horseradish peroxidase. Biochem J. 1986, 234: 335-342.

    PubMed Central  CAS  PubMed  Article  Google Scholar 

  123. 123.

    Bailey S, Ayling J: 6,6-Dimethylpterins: stable quinoid dihydropterin substrate for dihydropteridine reductase and tetrahydropterin cofactor for phenylalanine hydroxylase. Biochem. 1983, 22: 1790-1798. 10.1021/bi00277a008.

    CAS  Article  Google Scholar 

  124. 124.

    Firgaira F, Cotton R, Jennings I, Danks D: Use of naphthoquinone adsorbant for the isolation of human dihydropteridine reductase. Meth Enzymol. 1987, 142: 116-126. full_text.

    CAS  PubMed  Article  Google Scholar 

  125. 125.

    Firgaira F, Cotton R, Danks D: Isolation and characterization of dihydropteridine reductase from human liver. Biochem J. 1981, 97: 31-43.

    Article  Google Scholar 

  126. 126.

    Schomburg D, Schomburg I: 6,7-Dihydropteridine Reductase. Springer Hand-book of Enzymes. 2005, 23: 248-272. full_text.

    Article  Google Scholar 

  127. 127.

    Craine J, Hall E, Kaufman S: The isolation and characterization of dihydropteridine reductase from sheep liver. J Biol Chem. 1972, 247: 6082-6091.

    CAS  PubMed  Google Scholar 

Download references


This work was supported by NSF grant DMS-061670 (MR, HFN), NSF agreement 0112050 through the Mathematical Biosciences Institute (JB, MR), and NSF CAREER Award DMS-0956057 (JB). JB is an Alfred P. Sloan Research Fellow.

Author information



Corresponding author

Correspondence to Janet Best.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All three authors (JB, MR, HFN) contributed equally to the formulation of the model, the estimation of parameters, experimentation with the model, the biological interpretations and conclusions, and the writing and editing of the manuscript. All authors have read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Best, J., Nijhout, H.F. & Reed, M. Serotonin synthesis, release and reuptake in terminals: a mathematical model. Theor Biol Med Model 7, 34 (2010).

Download citation


  • Fluoxetine
  • Firing Rate
  • Dorsal Raphe Nucleus
  • Serotonergic System
  • Projection Region