- Open Access
Serotonin synthesis, release and reuptake in terminals: a mathematical model
© Best et al; licensee BioMed Central Ltd. 2010
- Received: 25 June 2010
- Accepted: 19 August 2010
- Published: 19 August 2010
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.
- Firing Rate
- Dorsal Raphe Nucleus
- Serotonergic System
- Projection Region
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. In addition, 5-HT has been linked to motor system function, sleep-wake cycles, circadian rhythms, respiratory stability, embryonic development, and reward processing. Not surprisingly, the 5-HT neurons in the nuclei originally classified by Dalhstrom and Fuxe 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 . 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[10–13] and may be associated with anxiety, autism, and suicidality[16, 17]. Polymorphisms in the tryptophan hydroxylase gene have been associated with unipolar and bipolar disorder. 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, and 5-HT1A mRNA levels are affected by gonadal hormones .
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. In addition, these regulatory mechanisms are themselves regulated by dynamic changes in autoreceptor expression levels. 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 showed that the patterns of firing of nucleus centralis superior serotonergic neurons correspond to different sleep-wake states. 5-HT modulates motor firing patterns 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, inhibits the firing of somatosensory cortical neurons, and can inhibit or excite neurons in the ventromedial nucleus of the hypothalamus. It has been proposed that 5-HT activates the hypothalamic-pituitary-adrenal axis by stimulating production of corticotropin-releasing hormone. 5-HT influences dopaminergic signaling[37, 38] and may affect firing in the cerebral cortex by causing the release of glutamate. Traditionally, dopamine was thought to be the primary neurotransmitter involved in reward processing, but recent work suggests a strong role for 5-HT. 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.
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 ; both models were inspired by the original model of Justice et al. 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. Waggoner and colleagues introduced a three state stochastic model for the serotonin dependence of egg laying in a nematode. Bunin et al. and Daws et al. 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. 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. 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. 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 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. Stoltenberg and Nag used a dynamical systems model to go directly from genes to behavior.
Names used for Variables
the tryptophan pool
the tryptophan pool
The Differential Equations
Kinetic Parameters (μ M, μ M/hr,/hr).
aromatic amino acid decarboxylase
K bh 4
vesicular monoamine transporter
K i (substrate inhibition)
neutral amino acid transporter
catabolism and diffusion
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  and a range of 61-173 μ M was found in . We take as our baseline the value of 96 μ M found by Fernstrom in fasted rats . 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 . 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  measured a K m = 190 μ M with respect to total serum tryptophan and Smith  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 . 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).
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:
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 , which is consistent with the values in . We choose the V max so that the concentration of cytosolic serotonin is very low. The experiments in  and the calculations in  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 . 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 .
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, 70–72], 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 . 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 . 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  that autoreceptors modulate reuptake, but this effect is not included in the model. Extracellular 5-HT or autoreceptor agonists can decrease synthesis by 50% [75–77] and by perhaps as much as 80-90% . 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% . 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 ; 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  and in the range 1-3 in . Since tissue content of 5-HT in different brain regions is roughly 2-3 μ M [86–88], the concentration 5hiaa = 5.22 μ M is reasonable.
The half-life of fluoxetine is quite long; 1-4 days is reported in . 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 . 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  to approximately 13 μ M/hr [83, 92].
It is known that the cytosolic concentrations of 5-HTP and 5-HT are quite low. Fernstrom  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 . 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 . We have previously constructed a model of dopamine (DA) synthesis, release, and reuptake in dopaminergic terminals , 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 [53–55] and also on the sequence of meals with different content . 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 . 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 . 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  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 . 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 . 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 , it makes sense that extracellular 5-HT should rise during and after meals.
B. Release and Reuptake
C. SERT Knockouts
Steady State Values from WT to SERT KO
f* = 1(WT)
f= 0(SERT KO)
c 5 - HT
v 5 - HT
e 5 - HT **
5 - hiaa
It is known that 5-HT tissue levels are down both in knockouts  and in WT mice treated with fluoxetine . Homberg et al. found that tissue levels of 5-HT drop by 50-70% and Bengel et al. 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  report that extracellular 5-HT rises 5-6 fold in SERT knockouts and Homberg et al. report a 9-fold increase . 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. 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
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 . 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  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 . 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% . A genetic polymorphism of the promoter, rs11178997, reduces TPH2 transcriptional activity by 22% . 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 . 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 [106–109].
Casanovas et al. 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 , and therefore firing is reduced much more in the DRN than in the MRN.
Malagie et al. 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. 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 . 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, 70–72], 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 . 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  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 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  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 . Thus, much slower uptake was found in the dentate gyrus, the corpus callosum and the CA3 region of the hippocampus . 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 , 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  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 . In  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; 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.
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.
- Feldman R, Meyer J, Quenzer L: Principles of Neuropharmacology. 1997, Sunderland, MA: Sinauer Associates, IncGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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 CentralPubMedView ArticleGoogle Scholar
- Kranz GS, Kasper S, Lanzenberger R: Reward and the serotonergic system. Neuroscience. 2010, 166 (4): 1023-1035. 10.1016/j.neuroscience.2010.01.036.PubMedView ArticleGoogle Scholar
- 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
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- White KJ, Walline CC, Barker EL: Serotonin transporters: implications for antidepressant drug development. AAPS J. 2005, 7: E421-E433. 10.1208/aapsj070242.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Murphy DL, Lesch KP: Targeting the murine serotonin transporter: insights into human neurobiology. Nature Rev Neurosci. 2008, 9: 85-96. 10.1038/nrn2284.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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/sj.mp.4001336.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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
- 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 CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Jacobs BL, Fornal CA: 5-HT and motor control: a hypothesis. Trends Neurosci. 1993, 16: 346-352. 10.1016/0166-2236(93)90090-9.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Best JA, Nijhout HF, Reed MC: Models of dopaminergic and serotonergic signaling. Pharmacopsychiatry. 2010, 43 (Supp 1): 561-566.Google Scholar
- 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 CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Butera RJ: Implementation and analysis of neuromodulatory mechanisms in a mathematical model of neuron R15 in Aplysia. Master's thesis, Rice University. 1984Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Dayan P, Huys QJM: Serotonin, inhibition, and negative mood. PLoS Comp Biol. 2008, 4 (2): 1-11. 10.1371/journal.pcbi.0040004.View ArticleGoogle Scholar
- Dayan P, Huys QJM: Serotonin in affective control. Ann Rev Neurosci. 2009Google Scholar
- 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 CentralPubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Kilberg M, Haussinger D: Mammalian Amino Acid Transport, Mechanisms and Control. 1992, New York: Plenum PrView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Sirek A, Sirek OV: Serotonin: A Review. Canad Med Ass J. 1970, 102: 846-849.PubMed CentralPubMedGoogle Scholar
- 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
- 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 CentralPubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Schuldiner S, Shirvan A, Linial M: Vesicular neurotransmitter transporters: from bacteria to humans. Physiological Reviews. 1995, 75 (2): 369-392.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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 CentralPubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Starke K, Gothert M, Kilbinger H: Modulation of neurotransmitter release by presynaptic autoreceptors. Physiol Rev. 1989, 69: 864-989.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Galloway M: Regulation of dopamine and serotonin synthesis by acute administration of cocaine. Synapse. 1990, 6: 63-72. 10.1002/syn.890060108.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- Parsons LH, Weiss F, Koob GF: Serotonin 1B receptor stimulation enhances dopamine-mediated reinforcement. Psycopharmacology. 1996, 128: 150-160. 10.1007/s002130050120.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Hiemke C, Hartter S: Pharmacokinetics of selective serotonin reuptake inhibitors. Pharmaco Therap. 2000, 85: 11-28. 10.1016/S0163-7258(99)00048-0.View ArticleGoogle Scholar
- Leathwood PD: Tryptophan availability and serotonin synthesis. Proc Nutr Soc. 1987, 46: 143-156. 10.1079/PNS19870018.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Fernstrom J, Wurtman RD: Brain serotonin content: physiological dependence on plasma tryptophan levels. Science. 1971, 173: 149-152. 10.1126/science.173.3992.149.PubMedView ArticleGoogle Scholar
- Fernstrom J: Role of precursor availability in control of monoamine biosynthesis in brain. Physiol Rev. 1983, 63: 484-546.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Cooper J, Bloom F, Roth R: The Biochemical Basis of Neuropharmacology. 2003, New York, NY: Oxford U PressGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Winge I, McKinney J, Knappskog P, Haavik J: Characterization of wild-type and mutant forms of human tryptophan hydroxylase 2. J Neurochem. 2007Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Fernstrom J: Branched chain amino acids and brain function. J Nutr. 2005Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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 CentralPubMedView ArticleGoogle Scholar
- 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 CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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 CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Armarego W, Ohnishi A, Taguchi H: New pteridine substrates for dihydropteridine reductase and horseradish peroxidase. Biochem J. 1986, 234: 335-342.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Firgaira F, Cotton R, Danks D: Isolation and characterization of dihydropteridine reductase from human liver. Biochem J. 1981, 97: 31-43.View ArticleGoogle Scholar
- Schomburg D, Schomburg I: 6,7-Dihydropteridine Reductase. Springer Hand-book of Enzymes. 2005, 23: 248-272. full_text.View ArticleGoogle Scholar
- Craine J, Hall E, Kaufman S: The isolation and characterization of dihydropteridine reductase from sheep liver. J Biol Chem. 1972, 247: 6082-6091.PubMedGoogle Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.