Next Article in Journal
Reductive Stress of Sulfur-Containing Amino Acids within Proteins and Implication of Tandem Protein–Lipid Damage
Next Article in Special Issue
A Lumenal Loop Associated with Catalytic Asymmetry in Plant Vacuolar H+-Translocating Pyrophosphatase
Previous Article in Journal
Structural Modelling of KCNQ1 and KCNH2 Double Mutant Proteins, Identified in Two Severe Long QT Syndrome Cases, Reveals New Insights into Cardiac Channelopathies
Previous Article in Special Issue
Ensemble of Template-Free and Template-Based Classifiers for Protein Secondary Structure Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DOPA Homeostasis by Dopamine: A Control-Theoretic View

1
Norwegian Center for Maritime and Diving Medicine, Haukeland University Hospital, 5021 Bergen, Norway
2
Department of Chemistry, Bioscience and Environmental Engineering, University of Stavanger, 4021 Stavanger, Norway
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2021, 22(23), 12862; https://doi.org/10.3390/ijms222312862
Submission received: 11 October 2021 / Revised: 23 November 2021 / Accepted: 23 November 2021 / Published: 28 November 2021
(This article belongs to the Collection Computational Studies of Biomolecules)

Abstract

:
Dopamine (DA) is an important signal mediator in the brain as well as in the periphery. The term “dopamine homeostasis” occasionally found in the literature refers to the fact that abnormal DA levels can be associated with a variety of neuropsychiatric disorders. An analysis of the negative feedback inhibition of tyrosine hydroxylase (TH) by DA indicates, with support from the experimental data, that the TH-DA negative feedback loop has developed to exhibit 3,4-dihydroxyphenylalanine (DOPA) homeostasis by using DA as a derepression regulator. DA levels generally decline when DOPA is removed, for example, by increased oxidative stress. Robust DOPA regulation by DA further implies that maximum vesicular DA levels are established, which appear necessary for a reliable translation of neural activity into a corresponding chemical transmitter signal. An uncontrolled continuous rise (windup) in DA occurs when Levodopa treatment exceeds a critical dose. Increased oxidative stress leads to the successive breakdown of DOPA homeostasis and to a corresponding reduction in DA levels. To keep DOPA regulation robust, the vesicular DA loading requires close to zero-order kinetics combined with a sufficiently high compensatory flux provided by TH. The protection of DOPA and DA due to a channeling complex is discussed.

1. Introduction

Dopamine (DA) is an important neurotransmitter. Its dysfunction is associated with a variety of neuropsychiatric disorders. Low activity levels of DA are related to movement disorders, irregular sleep patterns, or attention deficit hyperacticity disorder (ADHD). On the other hand, increased or dysregulated DA activity in different brain areas have been related to addiction [1] and schizophrenia [2,3] but also to improved cognitive functions [4]. Thus, the literature often refers to DA homeostasis in order to address the need for optimum DA levels to assure proper nervous functions.
The enzyme tyrosine hydroxylase (TH) converts tyrosine (Tyr) to the precursor 3,4-dihydroxy-phenylalanine (DOPA), which becomes DA by decarboxylation using dopa decarboxylase (DDC) (Figure 1). Cellular DA is moved into vesicles by vesicular monoamine transporter 2 (VMAT2); the vesicles are transported to the active site possibly by kinesins [5]. There, vesicles are emptied into the synaptic cleft upon the arrival of action potentials and the inflow of Ca 2 + ions [6]. There is evidence that the reaction chain Tyr → DOPA → vesicular DA inside the presynaptic cell is organized in terms of “metabolite channeling” [7], where the enzymes TH, DDC, and the transport VMAT2 form a complex attached to the vesicular membrane [8,9,10,11,12]. In the synaptic cleft DA is spread by diffusion and by reuptake into the presynaptic terminals with DA transporters (DAT) or by uptake by neighboring glial cells where it is metabolized [13,14,15].
There are several regulatory interactions which act on TH [23,31,32]. In addition, DOPA and dopaquinone have been reported to act as inhibitors [23,33], although the TH inhibition by quinones appears to be related to irreversible modifications of TH’s sulfhydryl groups [33]. In addition, TH has been found to possess DOPA oxidase activity [34] similar to the enzyme Tyrosinase (TYR) [28,29]. TYR, though expressed at low levels in dopaminergic neurons, uses DOPA as a substrate and seems to form neuromelanin [35].
DA synthesis and metabolism have been modeled in several studies addressing different regulatory aspects and model approaches. One of the earlier models of DA synthesis, release, and metabolism used the assumption that DOPA-derived DA enters the vesicular pool without mixing with cytosolic DA [36]. Three DA pools were considered, i.e., inactive bound DA, releasable bound DA, and free cytosolic DA. The TH regulation of DA synthesis was modeled [37], with an emphasis on regulation by phosphorylation/dephosphorylation, the cofactor BH4, and the cell’s redox status. Multiple activity states of TH were considered [37] together with influences by iron and α -synuclein. These enzymatic models are based on Michaelis–Menten (MM) kinetics. Others models described DA synthesis and metabolism by a power-law approach (biochemical systems theory) developed by Savageau [38]. The calculations by Qi et al. [30,39,40] contain an analysis of presynaptic DA homeostasis including details on DA metabolites, catechol auto-oxidation, and melanin formation. Best et al. [27] modeled DA synthesis and metabolism including TH regulation by autoreceptors. Cullen and Wong-Li [41] analyzed a reduced version of the model by Best et al. [27], which is considered to be computationally more efficient and addressing more clearly underlying key mechanisms. Véronneau-Veilleux et al. [42] developed an integrative model for Parkinson’s disease by integrating a pharmacokinetic Levodopa/DOPA treatment with DA dynamics and a neurocomputational model of basal ganglia. This illustrates nicely how different models can be combined to bridge interactions between different subsystems and to obtain clinically relevant predictions. Muddapu et al. [43,44] used a multiscale approach to model a network between the brain regions Substantia Nigra Pars Compacta (SNc), Subthalamic Nucleus (STN), Globus Pallidus externa (GPe ) with a spiking neuron model. To study the loss of dopaminergic cells [44] and Levodopa-induced toxicity [45], the Muddapu et al. model contains details about ion channels, calcium buffering mechanisms, energy metabolism pathways, DA turnover processes, and calcium-induced apoptosis. Recently, a detailed single-cell model of SNc was used to study the influence of energy deficiency on its subcellular processes to obtain insights into the neurodegeneration in Parkinson’s disease [46]. Practically all modeling approaches are based on solving the rate equations in the form of ordinary differential equations (ODEs). Among the most detailed models in terms of metabolites are probably the studies by Qi et al. [30].
In the present study, we use a minimal model approach to understand possible key (and ideal) features of DA regulation, which might not be immediately visible from a more complex model. Focus is taken on the presynaptic TH-DA negative feedback loop from a robust homeostatic viewpoint. As described in more detail below, we show that low K M values for the transport of DA into vesicles (as indicated by experimental data) allow for robust homeostasis of DOPA by using DA as the regulator. An intact homeostasis of DOPA by DA leads, in addition, to a maximum DA loading of vesicles. The feedback regulation of DOPA by DA is explored with respect to influences of oxidative stress, aging, and Levodopa treatment on DOPA and DA levels.

Mechanisms Leading to Robust Homeostasis

Before we go into more detailed properties of the TH-DA negative feedback and to DOPA regulation, we describe, by taking a simple negative feedback loop as an example, the requirements for robust homeostasis.
The term homeostasis was introduced by Walter Cannon in 1929 as “the sum of mechanisms which keep the steady states in organisms within certain but narrow limits” [47,48]. With the development of cybernetics [49,50,51] and systems biology [52,53,54], feedback mechanisms have become central in our understanding of how homeostasis can be achieved in biological systems [55,56].
During the past 20 years, evidence has accumulated that organisms can employ regulation methods, which originally were discovered/invented in control engineering [57,58,59,60,61,62,63,64,65,66,67,68]. Integral control, as one of these engineering methods, was invented in the beginning of the 20th century for the robust steering of ships [69,70] but is now an essential element in practically all industrial control applications [57].
In integral control, the difference between the controlled variable A and its set-point A s e t is integrated (Figure 2). The integrated error ϵ is proportional to the concentration of the manipulated (controller) variable, which is used to oppose perturbations of the controlled variable A. It can be proven that for step-wise perturbations the controlled variable A always converges toward its set-point A s e t [57]. This is generally referred to as the ”robustness” of the controller.
How can error integration be performed in a biochemical system? There are presently three approaches to achieve integral control in reaction kinetics: one is based on zero-order kinetics in the removal of the manipulated variable E [58,60], while another is based on a first-order autocatalytic formation of E combined with its first-order removal kinetics [65,71,72]. A third approach uses so-called antithetic control, where two controller variables are involved [65,66,68]. We focus here on how robust control can be achieved by zero-order kinetics.
Figure 3 gives a simple example of a negative feedback loop between A and E. Component P is a precursor of the controlled variable A. This feedback loop can be considered as a minimum model of TH inhibition by DA where P, A, and E are, respectively, Tyr, DOPA, and DA.
The rate equations are:
d P d t = k 1 k 2 · P k 3 + P · k 4 k 4 + E k 10 · P
d A d t = k 2 · P k 3 + P · k 4 k 4 + E k 9 · A perturbation k 5 · A k 6 + A
d E d t = k 5 · A k 6 + A k 7 · E k 8 + E
The rate equation for E includes an MM-type transport/removal term k 7 · E / ( k 8 + E ) . When this term becomes zero order with respect to E, i.e., k 8 E such that E / ( k 8 + E ) 1 , the steady-state value of A is under robust control with set-point A s e t :
d E d t = k 5 · A k 6 + A k 7 = ! 0 A s s = A s e t = k 6 k 7 k 5 k 7
In relation to integral control (Figure 2), we can rewrite the rate equation for E in the following form (see Appendix A):
d E d t = γ A ( A s e t A )
with γ A = ( k 5 k 7 ) ( k 6 + A ) . ( A s e t A ) represents the error between the actual value of A and A s e t . However, unlike to the presentations of integral control for linear or simpler chemical models [57,64]; here, γ A is not a constant but depends upon A.
To maintain robust control in A the flux j 2 ( = k 2 P / ( k 3 + P ) ) needs to be large enough to be able to compensate for the perturbing flux j 9 = k 9 · A ; otherwise, the concentrations of A are below A s e t . This is illustrated in Figure 4, where step-wise perturbations by k 9 are applied during four phases. In the third phase, j 2 is too low to compensate for the increased perturbation ( k 9 = 5.0). In this phase, E is practically zero, and j 2 cannot be increased further by E-depression. However, a larger P concentration in phase 4 by an increased k 1 value allows one to reestablish homeostasis in A. E is increased and can again act as a derepressing regulator E and compensate for the outflow perturbation j 9 = k 9 · A .

2. Materials and Methods

2.1. Method of Calculation

Rate equations were solved by using the Fortran subroutine LSODE [73]. Plots and curve fits were generated with gnuplot (version 5.4.2, www.gnuplot.info) and edited with Adobe Illustrator CS6 (version 16.0.0, adobe.com). To make notations simpler, concentrations of compounds are generally denoted by compound names without square brackets. Concentrations of the Tyr → DOPA → DA reaction chain are given in μ M, and time units are in minutes.

2.2. Model of Tyrosine Hydroxylase Regulation

Figure 5 shows the considered model:
The DAergic midbrain neurons form highly extended and complex structures, in particular, their axonal arbors [79]. We can therefore expect considerable compartmentalization in these neurons. Our modeling approach does however not rely on compartmentalization or model this explicitly. Although metabolite intermediates are free in solution and act through mass action kinetics, we expect that local concentration differences of the metabolites can occur within the axonal structure. We tentatively assume that the enzymes TH, DDC and the transporter VMAT2 form a channeling complex as indicated in Figure 1. When using the term ’cytosolic DOPA’ or ’cytosolic DA’ later in the paper, we do not necessarily consider that DOPA and DA are homogeneously distributed in the cytoplasm, but due to the local distribution of the channeling complex at the vesicular membrane, we assume that the highest concentrations of free DA and DOPA are physically close to the complex (this situation is termed by Cornish-Bowden “dynamic channeling”; see Figure 13.30 in [7]). Despite the additional inhibitions acting on TH, the key element for a possible robust homeostasis is here the same as for the model in Figure 3, i.e., the occurrence of zero-order kinetics in the removal of the TH-inhibiting species DA. Figure 5 shows that this step is related to the loading of DA into the vesicles by transporter VMAT2. DOPA is thereby the potentially regulated species, and DA is the regulator. The transport of DA into the vesicles is ATP dependent, which is needed to build a necessary proton gradient by ATPase [10,80]. The ATP dependence is not included in the model.
Experimentally determined rate parameters and inhibition constants (see next sections) are used and include substrate (Tyr) inhibition of TH [24], the inhibition of TH by DA [22,23,31], and the inhibition of TH by extracellular DA via D2 autoreceptors [27]. Furthermore, DOPA is enzymatically removed both by TH [34] as well as by TYR. The experimentally determined K M value for the transport of DA into vesicles, which is critical for a robust regulation of DOPA by DA has indeed been found to be relatively low (0.29 μ M [81]; see also [82] for another low K M value with respect to vesicle loading). These K M values indicate that the TH-DA negative feedback may prove to be close to robust behavior.
ROS, such as OH · radicals or superoxide O 2 , play an additional factor for removing DOPA (leading to neuromelanin [30]), thereby perturbing DOPA homeostasis. The influence of external DOPA addition by mimicking Levodopa medication (with rate constant k 22 ) is included in the model together with the re-entry of extracellular DA into the cell. More detailed descriptions of the experimental parameter/rate constant values are given in the next sections.
The rate equations are:
d ( T y r ) d t = k 1 k 2 · T y r k 3 1 + D A k 4 · 1 + D O P A k 16 + T y r k 12 k 12 + T y r k 20 k 20 + D A e x k 10 · T y r
d ( D O P A ) d t = k 2 · T y r k 3 1 + D A k 4 · 1 + D O P A k 16 + T y r k 12 k 12 + T y r k 20 k 20 + D A e x k 9 · D O P A k 13 + D O P A k 5 · D O P A k 6 + D O P A k 14 · D O P A k 15 · D O P A k 16 1 + T y r k 3 + D O P A + k 22
d ( D A ) d t = k 11 · D A e x k 21 + D A e x + k 5 · D O P A k 6 + D O P A k 7 · D A k 8 + D A α · k 18 · D A k 19 + D A + k 23 · ( D A v e s )
d ( D A v e s ) d t = k 7 · D A k 8 + D A k 17 · D A v e s k 23 · ( D A v e s )
d ( D A e x ) d t = k 17 · D A v e s k 18 · D A e x k 19 + D A e x
The grayed term in Equation (8) indicates the removal of DA at increased MAO concentrations with α being an adjustable factor. The set-point for DOPA ( D O P A s e t ) is calculated, in an analogous manner as for the P-A-E model (Figure 3), i.e., from the steady-state condition of Equation (8) with α = 0. The influence of α is shown below when considering combination aging, oxidative stress, and DA auto-oxidation.
In calculating the set-point D O P A s e t , we consider that the transport velocity of DA into the vesicles, given by k 7 · D A / ( k 8 + D A ) , is an ideal zero-order process with respect to DA. In other words, at low enough k 8 values, we have that k 7 · D A / ( k 8 + D A ) k 7 .
Then, inserting the DA re-entry flux
j 11 = k 11 · D A e x k 21 + D A e x
and the vesicular DA leakage (for details see next section and Appendix B)
j 23 = k 23 · ( D A v e s )
into Equation (8), setting d ( D A ) / d t = 0 , and solving for DOPA, gives the expression for the DOPA set-point:
D O P A s e t = k 6 ( k 7 j 11 j 23 ) k 5 + j 11 + j 23 k 7
Due to the re-entry and leakage terms j 11 and j 23 , D O P A s e t depends on the DA level in the synaptic cleft and on the amount of DA stored in vesicles. The maximum set-point of DOPA is reached when j 11 = 0 and j 23 = 0 with D O P A s e t m a x = k 6 k 7 / ( k 5 k 7 ) . D O P A s e t goes to zero when fluxes j 11 and j 23 balance k 7 , i.e., when there is no net-uptake of DA into vesicles. Thus, under ideal controller conditions (low k 8 ), D O P A s e t reflects the net-uptake of DA into vesicles.

2.3. Parameter Values

To the extent these were found in the literature, experimentally determined rate parameter values were used in the calculations. While reported K M values can be directly applied, V m a x values depend on the actual enzyme concentrations, which are mostly unknown. Rate constants (see the scheme of the model, Figure 5) k 1 , k 10 , k 11 , k 14 , k 17 k 20 , and k 22 were used as adjustable parameters. In the following section, we give an overview of found the literature data.

2.3.1. Dopamine Transport into Vesicles and Leakage

Robust homeostasis of DOPA is based on zero-order kinetics for the transport of DA into vesicles (i.e., a low transporter K M value compared to the local DA concentration transporter VMAT2 is exposed to). Volz et al. [81] studied the DA uptake kinetics by VMAT2 using voltammetry. The V m a x and K M values were determined, respectively, to 1.9 fmol/(s·mg protein) and 0.289 μ M. Experimental determination of cytosolic DA levels were measured by capillary electrophoresis and voltammetry with values between 2 and 3 μ M [83,84,85] and 47–140 μ M [86]. Near [82] studied the binding of [2-3H]Dihydrotetrabenazine to bovine striatal synaptic vesicles. The kinetic determination of the equilibrium constant resulted in the value of 5.4 nM ( 5.4 × 10 3   μ M), which relates to the K M for vesicle loading under a rapid-equilibrium assumption.
There is evidence that catecholamine transmitters stored in vesicles leak into the cytoplasm [16,17,18,19,20]. An analysis of the kinetic data by Fried [16] and Schonn et al. [17] (Appendix B) show that the vesicular leakage of DA can be described as a first-order process with respect to D A v e s (Equation (12)). While both Fried and Schonn et al. data sets give k 23 values close to each other, we use k 23 = (0.0158 ± 0.0006) min 1 from the Schonn et al. data.

2.3.2. Dopamine Re-Entry by DAT

DA transport by DAT was reported to follow Michaelis–Menten kinetics with K M estimates varying between 0.2 to 2 μ M ( k 21 in the model) [77,78]. In our calculations, we treated the re-entry of DA by DAT as a zero-order process (low k 21 ). Although different K M values of DAT change D O P A s e t by changed j 11 fluxes (as the leakage flux j 23 of DA out from vesicles also does (see Equation (13)), we found that different D O P A s e t values due to different DAT K M ’s are still defended.

2.3.3. DOPA Decarboxylase (DCC) EC 4.1.1.28

BRENDA [87] lists K M values for human DDC of 0.028 [88] and 4.27 mM [89]. We take the multiplicative average of these rather different values and use k 6 = 0.028 · 4.27 mM = 346 μ M. The turnover number of the enzyme was determined to be 5.1 s 1 or 306 min 1 [88]. Assuming (rather arbitrarily) a DCC concentration of 10 8 M, we obtain a V m a x ( k 5 ) value of 3.06 μ M/min, which is used in the calculations.

2.3.4. Tyrosinase (TYR) EC 1.14.18.1

Three human K M entries for DOPA as a substrate are reported: 0.34 [28], and 0.48 and 0.49 mM [29]. We take the average and assign for k 13 a value of 437 μ M. The turnover number for this enzyme is reported to be 38.11 s 1 or 2286 min 1 [28]. Using an enzyme concentration of 10 9 M, the used V m a x ( k 9 ) value in the calculations is 2.3 μ M/min.

2.3.5. Tyrosine Hydroxylase (TH) EC 1.14.16.2. Inhibition by Tyr

BRENDA has 5 entries for the inhibition constant k 12 of TH by Tyrosine. These are data of four mutant enzymes and one wild-type enzyme with the values [24]: 0.037, 0.044, 0.046 (wild type), 0.048, and 0.073 mM. The average value is 0.0496 mM. Fossbakk [25] reports a similar inhibition constant as [24] for the wild-type enzyme (59 ± 18 μ M) and, in addition, for a variety of TH point mutations. In the calculations, we use a value of 50 μ M for k 12 .

2.3.6. Inhibition of TH by DA

Using bovine TH, Lazar et al. [22] reports an inhibition (dissociation) constant for DA, k 4 in our model, in dependence of the enzyme’s phosphorylation status. At pH 7.0, the inhibition constant is 6 ± 2 μ M for the nonphosphorylated form and 45 ± 7 μ M when the enzyme is phosphorylated. Sura et al. [23] report for human TH isoforms inhibition constants for DA in the range 78–208 μ M for S40 phosphorylated TH. For the unphosphorylated TH forms, the inhibition constants are significantly lower, i.e., lower than 3.5 nM.

2.3.7. Inhibition of TH by DOPA

Sura et al. [23] determined the binding of DOPA to human TH isoforms. For the unphosphorylated enzyme, the K d is reported to be between 2.7 and 4.5 μ M. For phosphorylated TH, the K d values are slightly higher and lie in the range 4.0–7.4 μ M. Unlike the binding of DA to phosphorylated and unphosphorylated TH, the binding of DOPA shows only minor changes in dependence of TH’s phosphorylation state on TH activity. Sura et al. [23] conclude that DOPA appears not to be important for regulation of TH in vivo, whereas TH inhibition by DA is.

2.3.8. TH Turnover Number/ V m a x for Tyr as Substrate

The turnover number of Tyr for the wild-type TH is reported to be 3.33 s 1 ≈ 200 min 1 [24]. Sura et al. [31] report a similar value for the Tyr turnover number of 3.03 s 1 . We use the 200 min 1 value for TH’s turnover number and show how V m a x ( k 2 ) in dependence of TH concentration influences the regulation of DOPA by DA.

2.3.9. TH K M (Tyr) Values, k 3

BRENDA gives an overview of 16 K M values for human TH (using L-Tyrosine as substrate), ranging from 0.0081 to 0.166 mM. The values cover 11 mutant enzymes and five recombinant isoforms of TH. The isoform K M values are 0.011 mM [90] and 0.055, 0.066, 0.074, and 0.166 mM [91] with the average of 0.0744 mM = 74.4 μ M. The K M values for the different mutant and wild-type enzymes are: 0.0081, 0.014, 0.015, 0.016, 0.019, 0.02, 0.025 [92], 0.017 [93], 0.039, 0.042, 0.043, 0.046 (’wild-type’), 0.066 [24], and 0.08 mM (’wild type’) [23]. Sura et al. [23] also provide an overview of K M values for different mutants and the wild-type enzyme. Fossbakk et al. [25] reports a K M of 5.1 ± 2.2 μ M, while Morgenroth et al. [94] had reported earlier a K M for Tyr of 55.3 ± 4.2 μ M similar to one of the above TH isoform values by Nasrin et al. [91]. For k 3 , we use here the average value by Nasrin et al. [91] of 74.4 μ M.

2.3.10. TH-Mediated Conversion of DOPA

Haavik [34] reports a K M ( k 16 ) for the DOPA oxidation by TH of 56 ± 12 μ M. The rate of disappearance was found to be very low. In our calculations, we used a V m a x ( k 15 ) for the removal of DOPA by TH to be the same as the V m a x ( k 9 ) for TYR, i.e., 2.3 μ M/min.

2.3.11. Monoamine Oxidase (MAO) EC 1.4.3.4

Li et al. [95] studied the liver enzyme at 11 C. They determined the K M ( k 19 in our model) to 111 μ M and the k c a t (turnover number) ( k 19 ) to 0.249 s 1 ≈ 15 min 1 . Making a rough estimate by using Van ’t Hoff’s rule that for each increase by 10 °C, the rate constants doubles, we estimate k c a t roughly to 90 min 1 ( = 2 × 2 × 1.5 × 15 min 1 ). Assuming a MAO concentration of 1 μ M, we obtain 90 μ M/min for V m a x ( k 18 ).

3. Results

3.1. Factors Influencing DOPA Homeostasis

Since experimental parameter values may vary depending on the studying conditions/methods, we investigated the system’s homeostatic behavior when certain parameter values are varied.

3.1.1. K M of Dopamine Loading into Vesicles

As indicated in the derivation of D O P A s e t above, an important parameter for DOPA homeostasis is the K M value for the transport of DA into vesicles ( k 8 , Equation (8)). A necessary, but not sufficient condition for robust homeostasis, is that k 8 D A , i.e., the kinetics of DA loading should be zero-order or close to zero order with respect to DA. To illustrate how the zero-order condition influences DOPA homeostasis, Figure 6 shows the controller’s behavior toward ROS perturbations (oxidative stress) by applying step-wise increased k 14 values during three successive phases (phase 1: k 14 = 0 min 1 ; phase 2: k 14 = 0.01 min 1 ; phase 3: k 14 = 0.1 min 1 ) and by varying k 8 around the experimentally determined value of 0.29 μ M [81].
The framed panel in Figure 6 shows the calculations using the experimentally determined value of k 8 = 0.29 μ M, while in panels (a) and (b), the k 8 ’s are lower by two and one orders of magnitude, respectively. In panel (d), the k 8 is one order of magnitude higher than the value determined by Volz et al. [81]. Figure 6 clearly shows that the controller’s accuracy to defend D O P A s e t is dependent on k 8 .

3.1.2. Influence of Compensatory Flux j 2

While a low K M value for the DA loading into vesicles ( k 8 ) is important for ensuring a good controller accuracy, a sufficiently high compensatory flux j 2
j 2 = k 2 · T y r k 3 1 + D A k 4 · 1 + D O P A k 16 + T y r k 12 k 12 + T y r k 20 k 20 + D A e x
is also required to allow a workable derepression by DA.
This is shown in Figure 7b where k 2 ( V m a x of TH) is increased by one order of magnitude. Figure 7a shows the response before the increase in k 2 . However, the improved controller accuracy comes at the expense of a slower response time and higher DA levels.
From Equation (14), we see that an increase in Tyr, for example by increasing the inflow of Tyr by k 1 , can have negative effects on the controller, due to the Tyr inhibition term k 12 / ( k 12 + T y r ) , which decreases j 2 with increasing Tyr concentrations. Taking Figure 6c as a reference and increasing the k 1 by one order of magnitude (from 10 to 100 μ M/min), Tyr increases by approximately one order of magnitude, which leads to an 11% reduction of j 2 in phase 3 and a significant offset of DOPA from D O P A s e t (Figure 8).

3.1.3. TH Inhibition by DA

Unphosphorylated TH is subject to a much stronger DA inhibition (lower inhibition constant) than phosphorylated TH. Since the reported values of the inhibition constants ( k 4 in the model) vary significantly (see Material and Methods), we were interested in the feedback loops’s behavior when k 4 is varied. Figure 9 shows, in comparison with the varied k 8 and k 14 values from Figure 6, the corresponding results when k 4 is decreased by one order of magnitude. Clearly, a tighter inhibition of TH by DA results in a much more rapid response by the system, but due to a decreased compensatory flux, DOPA levels are below the controller’s set-point at larger k 8 values. In addition, DA steady-state levels are lower in comparison with Figure 6. As in Figure 6, the outlined panel c shows the results with k 8 having the value of 0.29 μ M determined by Volz et al. [81].

3.1.4. Levodopa Treatment

Since increasing oxidative perturbations k 14 leads to a decrease in DOPA and thereby to an additional decrease in DA, we investigated the effect when external DOPA is added, e.g., when medication by Levodopa is applied (described by rate constant k 22 ). Adding external DOPA to the system has two effects:
(i)
When applied such that the combined fluxes j 5 + j 11 + j 23 are lower than j 7 ( j 7 is the rate of DA loading of vesicles, Figure 5), DOPA inflow helps to maintain DOPA homeostasis and slightly improves the performance of the controller/negative feedback. However, the improvement by DOPA addition is dependent on the controller accuracy (i.e., k 8 values). This is shown in Figure 10 where controller performances in absence of DOPA addition (panels a and b), and in its presence, (panels c and d) are compared for two different k 8 values. When k 8 is low and controller accuracy is high, DOPA addition improves controller performance and raises DA levels. In this case, all components of the model, including DA, are in a steady state.
(ii)
However, if DOPA inflow results in the flux condition j 5 + j 11 + j 23 > j 7 , then the controller breaks down, and DA levels start to grow continuously with a DOPA steady-state level above D O P A s e t . The controller tries to oppose the increased DOPA levels by downregulating the compensatory flux j 2 to zero with a continuous increase in DA. This behavior, also termed integral wind-up [57], is shown in Figure 11, when the DOPA inflow rate in phase 3 was increased to k 22 = 8 μ M/min. DOPA steady-state levels are now entirely uncontrolled. Since the compensatory flux j 2 is practically zero, the steady-state level of DOPA is now solely determined by the DOPA inflow rate ( k 22 ) and by the rates of DOPA removal.
As a result of the increasing DA levels during the wind-up, the rate of vesicular DA loading j 7 (Equation (9)) becomes zero-order ( k 8 D A ) with respect to DA, which has the effect that the vesicular concentration of DA ( D A v e s ) reaches its maximum steady-state value k 7 / ( k 17 + k 23 ) :
D A v e s = k 7 k 17 + k 23 D A k 8 + D A = D A v e s m a x D A k 8 + D A high DA D A v e s m a x
DOPA additions can restore DOPA levels when the activity of TH is low, such as in TH deficiency (THD) [96] or in DOPA-responsive dystonia (DRD) [97,98] patients. In the case that TH activity is not fully absent, DOPA additions can restore the homeostatic regulation of DOPA in presence of oxidative stress ( k 14 ). This is indicated in Figure 12 where the amount of TH in phases 2 and 3 (both panels) was decreased by two-orders of magnitude to k 2 = 1 × 10 2   μ M/min with k 14 = 0.1 μ M/min. By comparison, phases 2 do not have any DOPA addition, while DOPA addition occurs in phase 3 ( k 22 = 7 μ M/min). The breakdown of the controller is seen in phase 2 showing low DOPA and DA levels due to the insufficient TH activity. As a reference, in phase 1, sufficient TH is present ( k 2 = 1 × 10 4   μ M/min), and oxidative stress is absent ( k 14 = 0.0 μ M/min). Panels a and b differ by controller accuracy k 8 , showing that a lower k 8 (panel b) gives a better controller accuracy.
Since in the complete absence of TH activity the negative feedback is lacking, levodopa titrations need to be used to compensate for the missing levels of DOPA and DA.

3.1.5. Robust DOPA Homeostasis Implies Maximum Vesicular DA Loading

As DA decreases in order to counteract for an outflow perturbation in DOPA, the vesicular DA concentration D A v e s also decreases (Equation (15)). The accuracy of the DOPA controller, i.e., the value of k 8 (compared with Figure 6), then, according to Equation (15), determines how much the steady-state values of vesicle-loaded DA ( D A v e s ) change in relation to DA levels. Figure 13 shows D A v e s steady-state levels (in terms of percentage of its maximum value D A v e s m a x ) in relation to DA concentrations and k 8 . The black line 3 shows the D A v e s m a x D A relationship when k 8 = 0.29 μ M [81]. The other curves show the changes when k 8 is increased or decreased by one and two orders of magnitude.
Thus, a robust controller with low k 8 values keeps D A v e s levels close at its maximum k 7 / ( k 17 + k 23 ) (Equation (15)) even when cytosolic DA levels change. Figure 14 shows the decrease in vesicular DA concentration in dependence of the changed controller accuracy k 8 and applied oxidative stress. The same perturbations by k 14 are applied as in Figure 6. Curves 1, 2, and 3 refer to the conditions in Figure 6a,c,d, respectively.
One may argue that maximum filled DA vesicles also occur when external DOPA is applied in excess, as indicated by Figure 11 and Equation (15). However, when DOPA is in excess, DA shows the wind-up and becomes uncontrollably high. The resulting high DA and DOPA levels can lead to increased and unfavorable oxidative stress by reactions with MAO/ROS leading to cell toxicity and neuropsychiatric disorders [30,45,99,100].

3.1.6. Deteriorated DOPA Homeostasis by DA Removal/Auto-Oxidation

A removal of DA by MAO and by DA auto-oxidation disturbs DOPA homeostasis and may lead to the breakdown of the DOPA-DA control loop. Although the channeling complex may keep DOPA and DA protected against possible molecular attack by ROS, some “capturing” by ROS and/or MAO may occur.
Figure 15 shows how DA removal leads to controller breakdown. When α = 0, no DA removal occurs, and DOPA control by DA is operative. When MAO/auto-oxidation has access to DA ( α = 1), DOPA homeostasis breaks down leading to high DOPA and low DA levels (Figure 15, dashed lines).
The corresponding decrease in D A v e s levels is shown in Figure 16. In phase 3 when the perturbation k 14 is relative high ( k 14 = 0.1 μ M/min), there is a significant reduction in D A v e s when DA oxidation is present ( α = 1 ), compared with the situation DA oxidation does not occur ( α = 0 ).

4. Discussion

4.1. On DOPA Regulation by DA Derepression in Cells

Using PC12 cells, Goldstein et al. [101] applied different MAO inhibitors and analyzed their influence on a variety of metabolites, including DA and DOPA. The study clearly shows the ”inverse” relationship between DA and DOPA: by applying increasing amounts of each inhibitor, DA levels rise (because less DA is oxidized by MAO), while DOPA levels correspondingly decline. Interestingly, for higher clorgyline concentrations (larger than 10 nM), the trend is reversed, i.e., DA levels decline with increasing inhibitor concentration and DOPA levels correspondingly rise. These results support the view of DOPA regulation by DA repression/derepression. However, we were not able to find the literature data which would answer the question to what extent the TH-DA negative feedback leads to more or less robust DOPA homeostasis. To address this, more specific tests are required, such as varying expression of TH or applying perturbations, which specifically would only reduce DOPA levels.
Thus, we analyzed the TH-DA negative feedback loop (Figure 5) from an ideal working perspective: under what conditions would an ideal control be possible, and which of the species within the negative feedback loop would be controlled? We see one kinetic and one structural argument which would favor ideal control.
The kinetic argument for an ideal robust control of DOPA is a zero-order uptake of DA into the vesicles. It implies, as seen in Figure 6, that DOPA is a regulated species and DA acts, by derepression kinetics, as the regulator. The experimental data by Volz et al. [81] for the K M ( k 8 in the model) of the VMAT2 transporter indeed support this view, although the behavior of the controller when using the experimental K M value by Volz et al. only tangents an ideal DOPA control (Figure 6c). On the other hand, increased TH activity levels due to increased enzyme concentrations or a changed phosphorylation status can lead to an improved DOPA homeostasis (Figure 7).
The structural aspect is the indication that the enzymes TH, DCC, and the transporter VMAT2 appear to form a channeling complex [9,10,11]. The occurrence of metabolite channeling is well recognized [7] and is considered to be a central control mechanism of cellular metabolism [102]. There are two aspects of metabolic channeling which appear of importance in the context of DOPA regulation by DA. One is an effective DA filling mechanism of vesicles. The other is the minimization of environmental perturbations (for example, at increased oxidative stress) acting on DOPA and DA (Figure 15).

4.2. Why DOPA Homeostasis?

An argument suggested by Equation (15) is that under ideal controller conditions (low k 8 ), DOPA homeostasis allows an optimum DA filling of the vesicular stores, which seems to be of importance for a reliable translation of the electrical signal into a chemical one. In this respect, it is interesting that TH localization to VMAT2 is found to be regulated by its phosphorylation on Ser31 [12]. To what extent this regulation is also related to DA synthesis and DOPA homeostasis remains to be discovered. In addition, DOPA homeostasis could also be related to the regulation of flux control within a biosynthesis pathway originating from DOPA.
Figure 7 shows that TH is important for providing a sufficiently high compensatory flux j 2 for maintaining a functional DOPA controller. Since many cofactors are involved in the activity of TH, the dysregulation of these factors/genes has negative effects on DOPA control and DA levels. For example, the point mutation P301A and several others [25] show a decrease in TH’s activity ( V m a x ( T y r ) ) by two orders of magnitude, which certainly has an impact on DOPA homeostasis and the regulation of DA, as indicated above.

4.3. Oxidative Stress and Age

Reactive oxygen species are implicated in a variety of age-related diseases [103]. With regard to aging, DA and serotonin are those neurotransmitters that are mostly discussed. From early adulthood, DA levels decline steadily by around 10% per 10 years [104], which has been related to the increase in MAO activity with age [74]. Other factors that are implicated in the aging brain is the production of ROS by mitochondria [105]. ROS, which decrease DOPA levels, cause a reduction in DA, as the feedback loop tries to oppose the declining DOPA levels. When oxidative stress increases with age and removes DOPA, DA levels therefore also lead to a steady decrease. With an additional removal of DA by MAO, cytosolic and vesicular DA further decrease. As MAO produces hydrogen peroxide and radicals that may further decrease DOPA, a self-amplifying decrease in DA is expected to occur with age. On the other hand, there is a complex network of antioxidant molecules to keep an intracellular redox homeostasis [106], which needs to be considered in understanding the mechanisms of aging.
ATP is an important component which is required to maintain various homeostatic regulations, including protein (i.e., α -synuclein) homeostasis [107], the transport of DA into vesicles [10], and vesicle transport [5]. As ATP levels decline during aging [108], many of the homeostatic mechanisms decline with increased age.

4.4. Why Only DOPA Inflow Control?

Most of the homeostatic regulators in metabolism come in antagonistic pairs. Some examples are insulin and glucagon in blood sugar homeostasis or parathyroid hormone and calcitonin in blood calcium homeostasis (see the Supporting Material in [64] for descriptions of other controller pairs). Each individual controller in such a pair acts either as an inflow controller to allow minimum acceptable levels of the controlled metabolite or as an outflow controller to avoid that metabolite concentrations become too high [64]. As seen by the calculations above, the control of DOPA by DA opposes the removal of DOPA, for example by oxidative stress. The TH-DA negative feedback acts therefore as an inflow controller [64] to avoid that DOPA levels become too low. The reason for this is apparently that one needs DOPA for the making of DA, where the latter is needed in providing the chemical transmitter signal. However, this inflow control by DA does not oppose any DOPA levels above D O P A s e t . This is seen in Figure 10b, when Levodopa is applied in such a high dose leading to the controller shutdown and DA wind-up. The absence of an outflow controller which would restrict maximum DOPA levels seems to indicate that higher DOPA levels can be better tolerated than lower ones.

4.5. Role of Other TH Regulators

The presence of Tyr inhibition on TH has in our calculations a clear influence on DOPA homeostasis by reducing the compensatory flux j 2 when Tyr increases. This is seen in Figure 8 when the Tyr concentration is increased by a larger k 1 value. However, this does not mean that robust homeostasis in DOPA cannot be achieved in the presence of Tyr inhibition of TH. Interestingly, calculations by Best et al. [27] indicate that the substrate inhibition of TH by Tyr can work as an attenuator, i.e., the Tyr inhibition of TH ’flattens’ the reaction rate catalyzed by TH and smooths D A v e s concentration when Tyr levels fluctuate, for example, by the intake of meals.

5. Conclusions and Outlook

The results above suggest the view that DA regulates the homeostasis of DOPA by derepression. The critical parameter to achieve robust DOPA control is the K M ( k 8 ) for the VMAT2 transporter moving DA into vesicles. The results by Volz et al. [81] for this K M (0.289 μ M), and the value by Near ( 5.4 × 10 3 μ M, for the binding of [2-3H] dihydrotetrabenazine to bovine striatal synaptic vesicles) [82], indicate that these transporter K M s appear to be relative low and may work under zero-order conditions with respect to the uptaken compound. While k 8 is important for the accuracy of DOPA homeostasis (Figure 6), an increase in TH’s V m a x ( k 2 ) is able to improve controller accuracy even when k 8 is relatively high (Figure 7). However, a sufficiently high compensatory flux j 2 is needed to maintain DOPA homeostasis. Factors which decrease j 2 may lead to a poorer controller performance, as for example an increase in TH inhibition by Tyr. The dephosphorylation of TH results in a stronger TH inhibition by DA, which leads to a more rapid and stronger controller response (Figure 9).
An argument for robust DOPA homeostasis is the observation that a zero-order kinetic loading of DA into vesicles not only results in robust DOPA control but also in a maximum DA loading of vesicles. This appears to be necessary in order to reliably transform an incoming train of electrical signals into a chemical transmitter message. Although we presently cannot answer whether the in vivo system shows robust DOPA control, we feel that our findings points to a possible “functional intention” of DOPA control during evolution. This regulatory feature of DA synthesis could be of interest for treatment procedures using MAO inhibitors or DOPA supplementation.
We have so far not included the (oscillatory) dynamics of repeated Ca 2 + inflow, DA release into the synaptic cleft, and the influences on TH activity via auto-receptors and the re-entry of extracellular DA into the cell. Even under such oscillatory conditions, DOPA homeostasis may still be operative as calculations on related homeostatic controllers with derepression kinetics have shown [109,110]. We will study these oscillatory aspects in relation to DOPA/DA regulation in subsequent work.

Author Contributions

Conceptualization, R.K., Q.W. and P.R.; methodology, P.R.; software, P.R.; validation, R.K., Q.W. and P.R.; formal analysis, R.K., Q.W. and P.R.; investigation, R.K., Q.W. and P.R.; writing—original draft preparation, P.R.; writing—review and editing, R.K., Q.W. and P.R. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by the Western Norway Health Authorities (RK, project no. 912264) and by a postdoctoral grant to QW from the University of Stavanger.

Acknowledgments

We thank Aurora Martinez for comments.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

The Fortran programs with binaries or related Matlab versions (Windows 10, Mac OSX) are available from the corresponding author.

Abbreviations

The following abbreviations are used in this manuscript:
DAdopamine
DAvesvesicular dopamine
DAexextracellular dopamine
DATdopamine transporter
DDCDOPA dehydroxylase
DOPAL-3,4-dihydroxyphenylalanine
LevodopaL-3,4-dihydroxyphenylalanine
MAOmonoamine oxidase
ROSreactive oxygen species
ttime (min)
Tyrtyrosine
TYRtyrosinase
THtyrosine hydroxylase
VMAT2dopamine transporter located at vesicle membrane

Appendix A. Derivation of Equation (5)

We start with the dot notation of Equation (4) and rearrange:
E ˙ = k 5 · A k 6 + A k 7 = k 5 · A k 6 + A k 7 ( k 6 + A ) k 6 + A = k 5 · A k 6 k 7 k 7 · A k 6 + A
Multiplying the left side and last term of Equation (A1) by 1 k 6 k 7 gives:
E ˙ k 6 k 7 = ( k 5 k 7 ) · A k 6 k 7 1 k 6 + A
Noticing that
A s e t = k 6 k 7 k 5 k 7
and multiplying the left and right side of Equation (A2) by k 6 k 7 k 5 k 7 leads to:
E ˙ k 5 k 7 = ( A A s e t ) k 6 + A = ( A s e t A ) k 6 + A
Rearranging Equation (A4) then gives Equation (5), i.e.:
E ˙ = ( k 5 k 7 ) k 6 + A γ A · ( A s e t A ) error
If E ˙ 0 , then A A s e t .

Appendix B. Determination of Rate Constant k23

We analyzed the kinetic data by Schonn et al. [17] (Figure 2a, incubation with clomipramine and reserpine). The analyzed data (Figure A1) show an excellent fit to a first-order kinetic decrease in the vesicular transmitter.
Figure A1. The data by Schonn et al. [17] of Serotonin release in the presence of clomipramine and reserpine inhibitors (green dots). The data were fitted to the function f ( t ) = 100 · exp ( k 23 · t ) resulting in k 23 = (0.0158 ± 0.0006) min 1 . The resulting fit is shown by the red line.
Figure A1. The data by Schonn et al. [17] of Serotonin release in the presence of clomipramine and reserpine inhibitors (green dots). The data were fitted to the function f ( t ) = 100 · exp ( k 23 · t ) resulting in k 23 = (0.0158 ± 0.0006) min 1 . The resulting fit is shown by the red line.
Ijms 22 12862 g0a1
Another analysis of the earlier data by Fried (Figure 1 in [16]) of the vesicular release of Noradrenalin (Norepinephrine) in the presence of reserpine and Mg-ATP gives a very similar k 23 value (≈ 0.0136 min 1 ) as for the Schonn et al. data. In this case, we estimated the half-life ( t 1 / 2 ) of the assumed first-order kinetic release of Noradrenaline and calculated k 23 by the relationship
k 23 = ln 2 t 1 / 2
Figure A2. Replot of the Fried data (Figure 1 in [16]). The release of Noradrenalin was monitored in the presence of Mg-ATP and the inhibitor reserpine. k 23 is estimated from t 1 / 2 to be approximately 0.0136 min 1 .
Figure A2. Replot of the Fried data (Figure 1 in [16]). The release of Noradrenalin was monitored in the presence of Mg-ATP and the inhibitor reserpine. k 23 is estimated from t 1 / 2 to be approximately 0.0136 min 1 .
Ijms 22 12862 g0a2

References

  1. Wise, R.A.; Robble, M.A. Dopamine and addiction. Annu. Rev. Psychol. 2020, 71, 79–106. [Google Scholar] [CrossRef]
  2. Davis, K.L.; Kahn, R.S.; Ko, G.; Davidson, M. Dopamine in schizophrenia: A review and reconceptualization. Am. J. Psychiatry 1991, 148, 1474–1486. [Google Scholar]
  3. Grace, A.A. Dysregulation of the dopamine system in the pathophysiology of schizophrenia and depression. Nat. Rev. Neurosci. 2016, 17, 524–532. [Google Scholar] [CrossRef] [PubMed]
  4. Nieoullon, A. Dopamine and the regulation of cognition and attention. Prog. Neurobiol. 2002, 67, 53–83. [Google Scholar] [CrossRef]
  5. Friedman, M.H. Principles and Models of Biological Transport; Springer Science & Business Media: New York, NY, USA, 2008. [Google Scholar]
  6. Kandel, E.R.; Schwartz, J.H.; Jessell, T.M.; Siegelbaum, S.; Hudspeth, A.J.; Mack, S. Principles of Neural Science; McGraw-Hill: New York, NY, USA, 2000; Volume 4. [Google Scholar]
  7. Cornish-Bowden, A. Fundamentals of Enzyme Kinetics, 4th ed.; Wiley-VCH Verlag: Weinheim, Germany, 2012. [Google Scholar]
  8. Bowling, K.M.; Huang, Z.; Xu, D.; Ferdousy, F.; Funderburk, C.D.; Karnik, N.; Neckameyer, W.; O’Donnell, J.M. Direct binding of GTP cyclohydrolase and tyrosine hydroxylase: Regulatory interactions between key enzymes in dopamine biosynthesis. J. Biol. Chem. 2008, 283, 31449–31459. [Google Scholar] [CrossRef] [Green Version]
  9. Cartier, E.A.; Parra, L.A.; Baust, T.B.; Quiroz, M.; Salazar, G.; Faundez, V.; Egaña, L.; Torres, G.E. A biochemical and functional protein complex involving dopamine synthesis and transport into synaptic vesicles. J. Biol. Chem. 2010, 285, 1957–1966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Munoz, P.; Huenchuguala, S.; Paris, I.; Segura-Aguilar, J. Dopamine oxidation and autophagy. Park. Dis. 2012, 2012, 920953. [Google Scholar] [CrossRef]
  11. Chen, R.; Wei, J.; Fowler, S.C.; Wu, J.Y. Demonstration of functional coupling between dopamine synthesis and its packaging into synaptic vesicles. J. Biomed. Sci. 2003, 10, 774–781. [Google Scholar] [CrossRef]
  12. Jorge-Finnigan, A.; Kleppe, R.; Jung-Kc, K.; Ying, M.; Marie, M.; Rios-Mondragon, I.; Salvatore, M.F.; Saraste, J.; Martinez, A. Phosphorylation at serine 31 targets tyrosine hydroxylase to vesicles for transport along microtubules. J. Biol. Chem. 2017, 292, 14092–14107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mulvihill, K.G. Presynaptic regulation of dopamine release: Role of the DAT and VMAT2 transporters. Neurochem. Int. 2019, 122, 94–105. [Google Scholar] [CrossRef] [PubMed]
  14. Wayment, H.K.; Schenk, J.O.; Sorg, B.A. Characterization of extracellular dopamine clearance in the medial prefrontal cortex: Role of monoamine uptake and monoamine oxidase inhibition. J. Neurosci. 2001, 21, 35–44. [Google Scholar] [CrossRef] [Green Version]
  15. Sulzer, D.; Cragg, S.J.; Rice, M.E. Striatal dopamine neurotransmission: Regulation of release and uptake. Basal Ganglia 2016, 6, 123–148. [Google Scholar] [CrossRef] [Green Version]
  16. Fried, G. Noradrenaline release and uptake in isolated small dense cored vesicles from rat seminal ducts. Acta Physiol. Scand. 1981, 112, 41–46. [Google Scholar] [CrossRef] [PubMed]
  17. Schonn, J.S.; Desnos, C.; Henry, J.P.; Darchen, F. Transmitter uptake and release in PC12 cells overexpressing plasma membrane monoamine transporters. J. Neurochem. 2003, 84, 669–677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Eisenhofer, G.; Kopin, I.J.; Goldstein, D.S. Catecholamine metabolism: A contemporary view with implications for physiology and medicine. Pharmacol. Rev. 2004, 56, 331–349. [Google Scholar] [CrossRef] [PubMed]
  19. Wallace, L.J. A small dopamine permeability of storage vesicle membranes and end product inhibition of tyrosine hydroxylase are sufficient to explain changes occurring in dopamine synthesis and storage after inhibition of neuron firing. Synapse 2007, 61, 715–723. [Google Scholar] [CrossRef]
  20. Hondebrink, L.; Meulenbelt, J.; Timmerman, J.G.; Van Den Berg, M.; Westerink, R.H. Amphetamine reduces vesicular dopamine content in dexamethasone-differentiated PC12 cells only following l-DOPA exposure. J. Neurochem. 2009, 111, 624–633. [Google Scholar] [CrossRef]
  21. Meiser, J.; Weindl, D.; Hiller, K. Complexity of dopamine metabolism. Cell Commun. Signal. 2013, 11, 34. [Google Scholar] [CrossRef] [Green Version]
  22. Lazar, M.A.; Lockfeld, A.J.; Truscott, R.J.; Barchas, J.D. Tyrosine hydroxylase from bovine striatum: Catalytic properties of the phosphorylated and nonphosphorylated forms of the purified enzyme. J. Neurochem. 1982, 39, 409–422. [Google Scholar] [CrossRef] [PubMed]
  23. Sura, G.R.; Daubner, S.C.; Fitzpatrick, P.F. Effects of phosphorylation by protein kinase A on binding of catecholamines to the human tyrosine hydroxylase isoforms. J. Neurochem. 2004, 90, 970–978. [Google Scholar] [CrossRef] [PubMed]
  24. Royo, M.; Daubner, S.C.; Fitzpatrick, P.F. Effects of mutations in tyrosine hydroxylase associated with progressive dystonia on the activity and stability of the protein. Proteins 2005, 58, 14–21. [Google Scholar] [CrossRef] [Green Version]
  25. Fossbakk, A.; Kleppe, R.; Knappskog, P.M.; Martinez, A.; Haavik, J. Functional studies of tyrosine hydroxylase missense variants reveal distinct patterns of molecular defects in Dopa-responsive dystonia. Hum. Mutat. 2014, 35, 880–890. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Ford, C.P. The role of D2-autoreceptors in regulating dopamine neuron activity and transmission. Neuroscience 2014, 282, 13–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Best, J.A.; Nijhout, H.F.; Reed, M.C. Homeostatic mechanisms in dopamine synthesis and release: A mathematical model. Theor. Biol. Med. Model. 2009, 6, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Fogal, S.; Carotti, M.; Giaretta, L.; Lanciai, F.; Nogara, L.; Bubacco, L.; Bergantino, E. Human tyrosinase produced in insect cells: A landmark for the screening of new drugs addressing its activity. Mol. Biotechnol. 2015, 57, 45–57. [Google Scholar] [CrossRef] [PubMed]
  29. Han, H.Y.; Lee, J.R.; Xu, W.A.; Hahn, M.J.; Yang, J.M.; Park, Y.D. Effect of Cl on tyrosinase: Complex inhibition kinetics and biochemical implication. J. Biomol. Struct. Dyn. 2007, 25, 165–171. [Google Scholar] [CrossRef] [PubMed]
  30. Qi, Z.; Miller, G.W.; Voit, E.O. Computational systems analysis of dopamine metabolism. PLoS ONE 2008, 3, e2444. [Google Scholar] [CrossRef] [PubMed]
  31. Sura, G.R.; Lasagna, M.; Gawandi, V.; Reinhart, G.D.; Fitzpatrick, P.F. Effects of ligands on the mobility of an active-site loop in tyrosine hydroxylase as monitored by fluorescence anisotropy. Biochemistry 2006, 45, 9632–9638. [Google Scholar] [CrossRef] [Green Version]
  32. Quinsey, N.S.; Luong, A.Q.; Dickson, P.W. Mutational analysis of substrate inhibition in tyrosine hydroxylase. J. Neurochem. 1998, 71, 2132–2138. [Google Scholar] [CrossRef] [PubMed]
  33. Kuhn, D.M.; Arthur Jr, R.E.; Thomas, D.M.; Elferink, L.A. Tyrosine hydroxylase is inactivated by catechol-quinones and converted to a redox-cycling quinoprotein: Possible relevance to Parkinson’s disease. J. Neurochem. 1999, 73, 1309–1317. [Google Scholar] [CrossRef]
  34. Haavik, J. L-DOPA Is a Substrate for Tyrosine Hydroxylase. J. Neurochem. 1997, 69, 1720–1728. [Google Scholar] [CrossRef]
  35. Carballo-Carbajal, I.; Laguna, A.; Romero-Giménez, J.; Cuadros, T.; Bové, J.; Martinez-Vicente, M.; Parent, A.; Gonzalez-Sepulveda, M.; Peñuelas, N.; Torra, A.; et al. Brain tyrosinase overexpression implicates age-dependent neuromelanin production in Parkinson’s disease pathogenesis. Nat. Commun. 2019, 10, 973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Justice, J.B.; Nicolaysen, L.C.; Michael, A.C. Modeling the dopaminergic nerve terminal. J. Neurosci. Methods 1988, 22, 239–252. [Google Scholar] [CrossRef]
  37. Kaushik, P.; Gorin, F.; Vali, S. Dynamics of tyrosine hydroxylase mediated regulation of dopamine synthesis. J. Comput. Neurosci. 2007, 22, 147–160. [Google Scholar] [CrossRef]
  38. Savageau, M.A. Biochemical Systems Analysis: A Study of Function and Design in Molecular Biology; Addison-Wesley: Reading, MA, USA, 1976. [Google Scholar]
  39. Qi, Z.; Miller, G.W.; Voit, E.O. A mathematical model of presynaptic dopamine homeostasis: Implications for schizophrenia. Pharmacopsychiatry 2008, 41 (Suppl. S1), S89–S98. [Google Scholar] [CrossRef] [PubMed]
  40. Qi, Z.; Miller, G.W.; Voit, E.O. Computational modeling of synaptic neurotransmission as a tool for assessing dopamine hypotheses of schizophrenia. Pharmacopsychiatry 2010, 43 (Suppl. 1), S50–S60. [Google Scholar] [CrossRef] [PubMed]
  41. Cullen, M.; Wong-Lin, K. Integrated dopaminergic neuronal model with reduced intracellular processes and inhibitory autoreceptors. IET Syst. Biol. 2015, 9, 245–258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Véronneau-Veilleux, F.; Robaey, P.; Ursino, M.; Nekka, F. An integrative model of Parkinson’s disease treatment including levodopa pharmacokinetics, dopamine kinetics, basal ganglia neurotransmission and motor action throughout disease progression. J. Pharmacokinet. Pharmacodyn. 2021, 48, 133–148. [Google Scholar] [CrossRef]
  43. Muddapu, V.R.; Mandali, A.; Chakravarthy, V.S.; Ramaswamy, S. A Computational Model of Loss of Dopaminergic Cells in Parkinson’s Disease Due to Glutamate-Induced Excitotoxicity. Front. Neural Circuits 2019, 13, 11. [Google Scholar] [CrossRef] [Green Version]
  44. Muddapu, V.R.; Chakravarthy, V.S. A Multi-Scale Computational Model of Excitotoxic Loss of Dopaminergic Cells in Parkinson’s Disease. Front. Neuroinform. 2020, 14, 34. [Google Scholar] [CrossRef]
  45. Muddapu, V.R.; Vijayakumar, K.; Ramakrishnan, K.; Chakravarthy, V.S. A Computational Model of Levodopa-Induced Toxicity in Substantia Nigra Pars Compacta in Parkinson’s Disease. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  46. Muddapu, V.R.; Chakravarthy, V.S. Influence of energy deficiency on the subcellular processes of Substantia Nigra Pars Compacta cell for understanding Parkinsonian neurodegeneration. Sci. Rep. 2021, 11, 1754. [Google Scholar] [CrossRef] [PubMed]
  47. Cannon, B. Organization for Physiological Homeostasis. Physiol. Rev. 1929, 9, 399. [Google Scholar] [CrossRef]
  48. Langley, L.L. (Ed.) Homeostasis. Origins of the Concept; Dowden, Hutchinson & Ross, Inc.: Stroudsbourg, PA, USA, 1973. [Google Scholar]
  49. Ashby, W.R. An Introduction to Cybernetics; John Wiley: New York, NY, USA, 1956. [Google Scholar]
  50. Wiener, N. Cybernetics or Control and Communication in the Animal and the Machine; MIT Press: Cambridge, MA, USA, 1961. [Google Scholar]
  51. Cariani, P.A. The homeostat as embodiment of adaptive control. Int. J. Gen. Syst. 2009, 38, 139–154. [Google Scholar] [CrossRef]
  52. Kitano, H. Foundations of Systems Biology; The MIT Press: Cambridge, MA, USA; London, UK, 2001. [Google Scholar]
  53. Iglesias, P.A.; Ingalls, B.P. Control Theory and Systems Biology; MIT Press: Cambridge, MA, USA, 2010. [Google Scholar]
  54. Ingalls, B.P. Mathematical Modeling in Systems Biology: An Introduction; MIT Press: Cambridge, MA, USA, 2013. [Google Scholar]
  55. Hughes, G. (Ed.) Homeostasis and Feedback Mechanisms. In Symposia of the Society for Experimental Biology. Number XVIII; Academic Press: New York, NY, USA, 1964. [Google Scholar]
  56. Nijhout, H.F.; Best, J.A.; Reed, M.C. Systems biology of robustness and homeostatic mechanisms. Wiley Interdiscip. Rev. Syst. Biol. Med. 2019, 11, e1440. [Google Scholar] [CrossRef] [PubMed]
  57. Wilkie, J.; Johnson, M.A.; Katebi, R. Control Engineering: An Introductory Course; Palgrave Macmillan: Basingstoke, UK, 2002. [Google Scholar]
  58. Yi, T.M.; Huang, Y.; Simon, M.I.; Doyle, J. Robust perfect adaptation in bacterial chemotaxis through integral feedback control. Proc. Natl. Acad. Sci. USA 2000, 97, 4649–4653. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Yi, T.M. Constructing Mathematical Models of Biological Signal Transduction Pathways: An Analysis of Robustness. In Foundations of Systems Biology; Kitano, H., Ed.; The MIT Press: Cambridge, MA, USA; London, UK, 2001; Chapter 12; pp. 251–260. [Google Scholar]
  60. Ni, X.Y.; Drengstig, T.; Ruoff, P. The control of the controller: Molecular mechanisms for robust perfect adaptation and temperature compensation. Biophys. J. 2009, 97, 1244–1253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Jolma, I.W.; Ni, X.Y.; Rensing, L.; Ruoff, P. Harmonic oscillations in homeostatic controllers: Dynamics of the p53 regulatory system. Biophys. J. 2010, 98, 743–752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Ang, J.; Bagh, S.; Ingalls, B.P.; McMillen, D.R. Considerations for using integral feedback control to construct a perfectly adapting synthetic gene network. J. Theor. Biol. 2010, 266, 723–738. [Google Scholar] [CrossRef] [PubMed]
  63. Huang, Y.; Drengstig, T.; Ruoff, P. Integrating fluctuating nitrate uptake and assimilation to robust homeostasis. Plant Cell Environ. 2012, 35, 917–928. [Google Scholar] [CrossRef] [PubMed]
  64. Drengstig, T.; Jolma, I.W.; Ni, X.Y.; Thorsen, K.; Xu, X.M.; Ruoff, P. A basic set of homeostatic controller motifs. Biophys. J. 2012, 103, 2000–2010. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Briat, C.; Gupta, A.; Khammash, M. Antithetic Integral Feedback Ensures Robust Perfect Adaptation in Noisy Biomolecular Networks. Cell Syst. 2016, 2, 15–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Briat, C.; Zechner, C.; Khammash, M. Design of a Synthetic Integral Feedback Circuit: Dynamic Analysis and DNA Implementation. ACS Synth. Biol. 2016, 5, 1108–1116. [Google Scholar] [CrossRef] [PubMed]
  67. Krishnan, J.; Floros, I. Adaptive information processing of network modules to dynamic and spatial stimuli. BMC Syst. Biol. 2019, 13, 32. [Google Scholar] [CrossRef] [PubMed]
  68. Khammash, M.H. Perfect adaptation in biology. Cell Syst. 2021, 12, 509–521. [Google Scholar] [CrossRef] [PubMed]
  69. Lewis, F.L. Applied Optimal Control & Estimation: Digital Design & Implementation; Prentice Hall: Englewood Cliffs, NJ, USA, 1992. [Google Scholar]
  70. Bennett, S. A brief history of automatic control. IEEE Control Syst. Mag. 1996, 16, 17–25. [Google Scholar] [CrossRef]
  71. Drengstig, T.; Ni, X.Y.; Thorsen, K.; Jolma, I.W.; Ruoff, P. Robust Adaptation and Homeostasis by Autocatalysis. J. Phys. Chem. B 2012, 116, 5355–5363. [Google Scholar] [CrossRef] [PubMed]
  72. Shoval, O.; Goentoro, L.; Hart, Y.; Mayo, A.; Sontag, E.; Alon, U. Fold-change detection and scalar symmetry of sensory input fields. Proc. Natl. Acad. Sci. USA 2010, 107, 15995–16000. [Google Scholar] [CrossRef] [Green Version]
  73. Radhakrishnan, K.; Hindmarsh, A.C. Description and Use of LSODE, the Livermore Solver for Ordinary Differential Equations. NASA Reference Publication 1327, Lawrence Livermore National Laboratory Report UCRL-ID-113855; National Aeronautics and Space Administration, Lewis Research Center: Cleveland, OH, USA, 1993.
  74. Volchegorskii, I.A.; Shemyakov, S.E.; Turygin, V.V.; Malinovskaya, N.V. The age dynamics of monoamine oxidase activity and levels of lipid peroxidation products in the human brain. Neurosci. Behav. Physiol. 2004, 34, 303–305. [Google Scholar] [CrossRef] [PubMed]
  75. Serra, P.A.; Esposito, G.; Enrico, P.; Mura, M.A.; Migheli, R.; Delogu, M.R.; Miele, M.; Desole, M.S.; Grella, G.; Miele, E. Manganese increases L-DOPA auto-oxidation in the striatum of the freely moving rat: Potential implications to L-DOPA long-term therapy of Parkinson’s disease. Br. J. Pharmacol. 2000, 130, 937–945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Umek, N.; Geršak, B.; Vintar, N.; Šoštarič, M.; Mavri, J. Dopamine autoxidation is controlled by acidic pH. Front. Mol. Neurosci. 2018, 11, 467. [Google Scholar] [CrossRef] [Green Version]
  77. Jones, S.R.; Garris, P.A.; Kilts, C.D.; Wightman, R.M. Comparison of dopamine uptake in the basolateral amygdaloid nucleus, caudate-putamen, and nucleus accumbens of the rat. J. Neurochem. 1995, 64, 2581–2589. [Google Scholar] [CrossRef] [PubMed]
  78. Schmitz, Y.; Benoit-Marand, M.; Gonon, F.; Sulzer, D. Presynaptic regulation of dopaminergic neurotransmission. J. Neurochem. 2003, 87, 273–289. [Google Scholar] [CrossRef]
  79. Matsuda, W.; Furuta, T.; Nakamura, K.C.; Hioki, H.; Fujiyama, F.; Arai, R.; Kaneko, T. Single nigrostriatal dopaminergic neurons form widely spread and highly dense axonal arborizations in the neostriatum. J. Neurosci. 2009, 29, 444–453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Beyenbach, K.W.; Wieczorek, H. The V-type H+ ATPase: Molecular structure and function, physiological roles and regulation. J. Exp. Biol. 2006, 209, 577–589. [Google Scholar] [CrossRef] [Green Version]
  81. Volz, T.J.; Hanson, G.R.; Fleckenstein, A.E. Measurement of kinetically resolved vesicular dopamine uptake and efflux using rotating disk electrode voltammetry. J. Neurosci. Methods 2006, 155, 109–115. [Google Scholar] [CrossRef] [PubMed]
  82. Near, J.A. [3H]Dihydrotetrabenazine binding to bovine striatal synaptic vesicles. Mol. Pharmacol. 1986, 30, 252–257. [Google Scholar] [PubMed]
  83. Olefirowicz, T.M.; Ewing, A.G. Capillary electrophoresis in 2 and 5 microns diameter capillaries: Application to cytoplasmic analysis. Anal. Chem. 1990, 62, 1872–1876. [Google Scholar] [CrossRef] [PubMed]
  84. Olefirowicz, T.M.; Ewing, A.G. Dopamine concentration in the cytoplasmic compartment of single neurons determined by capillary electrophoresis. J. Neurosci. Methods 1990, 34, 11–15. [Google Scholar] [CrossRef]
  85. Olefirowicz, T.M.; Ewing, A.G. capillary electrophoresis for sampling signal nerve cell. Chimia 1991, 45, 106–108. [Google Scholar]
  86. Chien, J.B.; Wallingford, R.A.; Ewing, A.G. Estimation of Free Dopamine in the Cytoplasm of the Giant Dopamine Cell of Planorbis corneus by Voltammetry and Capillary Electrophoresis. J. Neurochem. 1990, 54, 633–638. [Google Scholar] [CrossRef] [PubMed]
  87. Chang, A.; Jeske, L.; Ulbrich, S.; Hofmann, J.; Koblitz, J.; Schomburg, I.; Neumann-Schaal, M.; Jahn, D.; Schomburg, D. BRENDA, the ELIXIR core data resource in 2021: New developments and updates. Nucleic Acids Res. 2021, 49, D498–D508. [Google Scholar] [CrossRef] [PubMed]
  88. Bertoldi, M. Mammalian Dopa decarboxylase: Structure, catalytic activity and inhibition. Arch. Biochem. Biophys. 2014, 546, 1–7. [Google Scholar] [CrossRef] [PubMed]
  89. Bertoldi, M.; Voltattorni, C.B. Multiple roles of the active site lysine of Dopa decarboxylase. Arch. Biochem. Biophys. 2009, 488, 130–139. [Google Scholar] [CrossRef]
  90. Martinez, A.; Abeygunawardana, C.; Haavik, J.; Flatmark, T.; Mildvan, A.S. Interaction of substrate and pterin cofactor with the metal of human tyrosine hydroxylase as determined by 1H-NMR. Adv. Exp. Med. Biol. 1993, 338, 77–80. [Google Scholar] [CrossRef] [PubMed]
  91. Nasrin, S.; Ichinose, H.; Hidaka, H.; Nagatsu, T. Recombinant human tyrosine hydroxylase types 1-4 show regulatory kinetic properties for the natural (6R)-tetrahydrobiopterin cofactor. J. Biochem. 1994, 116, 393–398. [Google Scholar] [CrossRef] [PubMed]
  92. Briggs, G.D.; Bulley, J.; Dickson, P.W. Catalytic domain surface residues mediating catecholamine inhibition in tyrosine hydroxylase. J. Biochem. 2014, 155, 183–193. [Google Scholar] [CrossRef] [PubMed]
  93. Ogawa, S.; Ichinose, H. Effect of metals and phenylalanine on the activity of human tryptophan hydroxylase-2: Comparison with that on tyrosine hydroxylase activity. Neurosci. Lett. 2006, 401, 261–265. [Google Scholar] [CrossRef] [PubMed]
  94. Morgenroth, V.H., III; Walters, J.R.; Roth, R.H. Dopaminergic neurons—Alteration in the kinetic properties of tyrosine hydroxylase after cessation of impulse flow. Biochem. Pharmacol. 1976, 25, 655–661. [Google Scholar] [CrossRef]
  95. Li, M.; Hubálek, F.; Newton-Vinson, P.; Edmondson, D.E. High-level expression of human liver monoamine oxidase A in Pichia pastoris: Comparison with the enzyme expressed in Saccharomyces cerevisiae. Protein Expr. Purif. 2002, 24, 152–162. [Google Scholar] [CrossRef] [PubMed]
  96. Ng, J.; Papandreou, A.; Heales, S.J.; Kurian, M.A. Monoamine neurotransmitter disorders—Clinical advances and future perspectives. Nat. Rev. Neurol. 2015, 11, 567–584. [Google Scholar] [CrossRef]
  97. Nygaard, T.G.; Marsden, C.; Duvoisin, R. Dopa-responsive dystonia. Adv. Neurol. 1988, 50, 377–384. [Google Scholar]
  98. Wijemanne, S.; Jankovic, J. Dopa-responsive dystonia—Clinical and genetic heterogeneity. Nat. Rev. Neurol. 2015, 11, 414–424. [Google Scholar] [CrossRef] [PubMed]
  99. Giannopoulos, S.; Samardzic, K.; Raymond, B.B.A.; Djordjevic, S.P.; Rodgers, K.J. L-DOPA causes mitochondrial dysfunction in vitro: A novel mechanism of L-DOPA toxicity uncovered. Int. J. Biochem. Cell Biol. 2019, 117, 105624. [Google Scholar] [CrossRef] [PubMed]
  100. Hörmann, P.; Delcambre, S.; Hanke, J.; Geffers, R.; Leist, M.; Hiller, K. Impairment of neuronal mitochondrial function by l-DOPA in the absence of oxygen-dependent auto-oxidation and oxidative cell damage. Cell Death Discov. 2021, 7, 151. [Google Scholar] [CrossRef]
  101. Goldstein, D.S.; Jinsmaa, Y.; Sullivan, P.; Holmes, C.; Kopin, I.J.; Sharabi, Y. Comparison of monoamine oxidase inhibitors in decreasing production of the autotoxic dopamine metabolite 3, 4-dihydroxyphenylacetaldehyde in PC12 cells. J. Pharmacol. Exp. Ther. 2016, 356, 483–492. [Google Scholar] [CrossRef] [Green Version]
  102. Obata, T. Toward an evaluation of metabolite channeling in vivo. Curr. Opin. Biotechnol. 2020, 64, 55–61. [Google Scholar] [CrossRef] [PubMed]
  103. Liguori, I.; Russo, G.; Curcio, F.; Bulli, G.; Aran, L.; Della-Morte, D.; Gargiulo, G.; Testa, G.; Cacciatore, F.; Bonaduce, D.; et al. Oxidative stress, aging, and diseases. Clin. Interv. Aging 2018, 13, 757–772. [Google Scholar] [CrossRef] [Green Version]
  104. Peters, R. Ageing and the brain. Postgrad. Med. J. 2006, 82, 84–88. [Google Scholar] [CrossRef]
  105. Melov, S. Modeling mitochondrial function in aging neurons. Trends Neurosci. 2004, 27, 601–606. [Google Scholar] [CrossRef]
  106. Finkel, T. Signal transduction by reactive oxygen species. J. Cell Biol. 2011, 194, 7–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Morimoto, R.I.; Selkoe, D.J.; Kelly, J.W. Protein Homeostasis; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 2011. [Google Scholar]
  108. Chaudhari, S.N.; Kipreos, E.T. The Energy Maintenance Theory of Aging: Maintaining Energy Metabolism to Allow Longevity. Bioessays 2018, 40, e1800005. [Google Scholar] [CrossRef] [PubMed]
  109. Drobac, G.; Waheed, Q.; Heidari, B.; Ruoff, P. An amplified derepression controller with multisite inhibition and positive feedback. PLoS ONE 2021, 16, e0241654. [Google Scholar] [CrossRef] [PubMed]
  110. Thorsen, K.; Agafonov, O.; Selsto, C.H.; Jolma, I.W.; Ni, X.Y.; Drengstig, T.; Ruoff, P. Robust concentration and frequency control in oscillatory homeostats. PLoS ONE 2014, 9, e107766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Overview of the synthesis and regulation of DA. There is evidence that DA is synthesized from tyrosine (Tyr) by a channeling mechanism with the enzymes TH, DDC, and the transporter VMAT2, forming a complex at the vesicular membrane [9,10,11]. The complex, whose details are presently not known, possibly minimizes the oxidations of DA and DOPA by MAO and reactive oxygen species (ROS), respectively [10]. It should however be noted that our analysis does not rely on the existence of such a protein complex. The synthesized DA is stored in synaptic vesicles for release. Significant leakage of transmitter molecules out from vesicles has been observed [16,17,18,19,20], and we tentatively assume that leakage also occurs in vivo. Upon arrival of an action, potential vesicles are emptied, and DA is released into the synaptic cleft. There, DA diffuses to DA-responsive target sites (postsynaptic neurons) or taken up by DA transporters (DAT) at neighboring dopaminergic terminals or metabolizing glial cells [10,14,21]. TH is inhibited both by DA [22,23] and its substrate Tyr [24,25]. There is also an inhibition of TH from extracellular DA via D2 autoreceptors that inhibits stimulatory cAMP/PKA phosphorylation of TH [26,27]. Tyrosinase (TYR) [28,29], converts DOPA to dopaquinone with the final formation of Neuromelanin [30]. Monoamine oxidase (MAO) forms hydrogen peroxide (a reactive oxygen species (ROS)) during the metabolization of DA. ROS are able to oxidize DOPA using a similar pathway as TYR, leading to neuromelanin (see [30] and references therein).
Figure 1. Overview of the synthesis and regulation of DA. There is evidence that DA is synthesized from tyrosine (Tyr) by a channeling mechanism with the enzymes TH, DDC, and the transporter VMAT2, forming a complex at the vesicular membrane [9,10,11]. The complex, whose details are presently not known, possibly minimizes the oxidations of DA and DOPA by MAO and reactive oxygen species (ROS), respectively [10]. It should however be noted that our analysis does not rely on the existence of such a protein complex. The synthesized DA is stored in synaptic vesicles for release. Significant leakage of transmitter molecules out from vesicles has been observed [16,17,18,19,20], and we tentatively assume that leakage also occurs in vivo. Upon arrival of an action, potential vesicles are emptied, and DA is released into the synaptic cleft. There, DA diffuses to DA-responsive target sites (postsynaptic neurons) or taken up by DA transporters (DAT) at neighboring dopaminergic terminals or metabolizing glial cells [10,14,21]. TH is inhibited both by DA [22,23] and its substrate Tyr [24,25]. There is also an inhibition of TH from extracellular DA via D2 autoreceptors that inhibits stimulatory cAMP/PKA phosphorylation of TH [26,27]. Tyrosinase (TYR) [28,29], converts DOPA to dopaquinone with the final formation of Neuromelanin [30]. Monoamine oxidase (MAO) forms hydrogen peroxide (a reactive oxygen species (ROS)) during the metabolization of DA. ROS are able to oxidize DOPA using a similar pathway as TYR, leading to neuromelanin (see [30] and references therein).
Ijms 22 12862 g001
Figure 2. The principle of integral control in negative feedback regulation. Variable A is the controlled variable with set-point A s e t . The error ϵ = A s e t A is integrated in time, which leads to the manipulated variable E, which can oppose perturbations affecting A. For step-wise perturbations, integral control leads to the precise regulation of A to A s e t [57].
Figure 2. The principle of integral control in negative feedback regulation. Variable A is the controlled variable with set-point A s e t . The error ϵ = A s e t A is integrated in time, which leads to the manipulated variable E, which can oppose perturbations affecting A. For step-wise perturbations, integral control leads to the precise regulation of A to A s e t [57].
Ijms 22 12862 g002
Figure 3. Negative feedback loop with robust homeostasis of controlled variable A. Same color code as in Figure 2 indicating the parameters which contribute to A s e t (red), the controlled variable A (blue) and the manipulated variable E (green). Outlined in brown is the parameter which leads to the integration of the error ϵ in Figure 2. The error integration is achieved when the manipulated variable E is removed by zero-order kinetics with respect to E, i.e., k 8 E (see Equations (3)–(5) in main text). The pair of rate constants k i , k j represent the respective V m a x and K M values of applied Michaelis–Menten kinetics. Compound P is the source of the compensatory flux j 2 . j 2 needs to have a sufficiently high value in order to compensate for the step-wise perturbations k 9 which remove A.
Figure 3. Negative feedback loop with robust homeostasis of controlled variable A. Same color code as in Figure 2 indicating the parameters which contribute to A s e t (red), the controlled variable A (blue) and the manipulated variable E (green). Outlined in brown is the parameter which leads to the integration of the error ϵ in Figure 2. The error integration is achieved when the manipulated variable E is removed by zero-order kinetics with respect to E, i.e., k 8 E (see Equations (3)–(5) in main text). The pair of rate constants k i , k j represent the respective V m a x and K M values of applied Michaelis–Menten kinetics. Compound P is the source of the compensatory flux j 2 . j 2 needs to have a sufficiently high value in order to compensate for the step-wise perturbations k 9 which remove A.
Ijms 22 12862 g003
Figure 4. Demonstration of robust A homeostasis in response to step-wise changes of k 9 . To keep A at its set-point at increased k 9 values, the concentration of E has to decrease by derepression under zero-order kinetics with respect to E. Upper panel: concentration of P as a function of time. Lower panel: concentrations of regulated species A (in blue) and controller species E (in green) as functions of time. For illustration purposes, the rate constants were chosen such, that A s e t = 1.0. Phase 1: Perturbation k 9 = 0.0. Rate constant k 1 represents the inflow of P, which is changed from 5.0 to 10.0 at the transition from phase 3 to phase 4. Phase 2: k 9 = 2.0; phases 3 and 4: k 9 = 5.0. Robust homeostasis in A breaks down in phase 3, because there is not sufficient P required to generate a sufficiently large compensatory flux j 2 . When k 1 is increased in phase 4, homeostasis of A is reestablished by the increase in P and j 2 . Other rate constants are: k 2 = 10.0, k 3 = 10.0, k 4 = 1.0, k 5 = 2.0, k 6 = 1.0, k 7 = 1.0, k 8 = 1 × 10 6 , and k 10 = 1 × 10 2 . Initial concentrations: P 0 = 400.0, A 0 = 1.0, and E 0 = 8.76. All values are in arbitrary units (au).
Figure 4. Demonstration of robust A homeostasis in response to step-wise changes of k 9 . To keep A at its set-point at increased k 9 values, the concentration of E has to decrease by derepression under zero-order kinetics with respect to E. Upper panel: concentration of P as a function of time. Lower panel: concentrations of regulated species A (in blue) and controller species E (in green) as functions of time. For illustration purposes, the rate constants were chosen such, that A s e t = 1.0. Phase 1: Perturbation k 9 = 0.0. Rate constant k 1 represents the inflow of P, which is changed from 5.0 to 10.0 at the transition from phase 3 to phase 4. Phase 2: k 9 = 2.0; phases 3 and 4: k 9 = 5.0. Robust homeostasis in A breaks down in phase 3, because there is not sufficient P required to generate a sufficiently large compensatory flux j 2 . When k 1 is increased in phase 4, homeostasis of A is reestablished by the increase in P and j 2 . Other rate constants are: k 2 = 10.0, k 3 = 10.0, k 4 = 1.0, k 5 = 2.0, k 6 = 1.0, k 7 = 1.0, k 8 = 1 × 10 6 , and k 10 = 1 × 10 2 . Initial concentrations: P 0 = 400.0, A 0 = 1.0, and E 0 = 8.76. All values are in arbitrary units (au).
Ijms 22 12862 g004
Figure 5. Reaction scheme of the TH-DA negative feedback model with added rate constants and fluxes. Not all DOPA and DA are necessarily bound to the channeling complex, but since the complex is not evenly distributed within the cytoplasm, both cytosolic DOPA and DA are considered to be located proximate to the complex. With increasing age and increased MAO levels [74] MAO may oxidize DA and perturb the feedback control of DOPA (path outlined in gray). In addition, both DOPA and DA are subject to auto-oxidation [75,76], which are included in the model by k 14 (together with ROS) and α ( k 18 ,   k 19 ) (together with MAO). As in Figure 3, pairs of rate constants ( k i , k j ) represent the respective ( V m a x , K M ) values of applied Michaelis–Menten kinetics. With increased MAO concentrations, oxidative stress and ROS also increase due to MAO-generated hydrogen peroxide. ROS can oxidize DOPA [30] (outlined in red). Cytosolic DA concentration is increased due to the re-entry of DA (flux j 11 ) and by vesicular leakage [16,17,18,19,20] (flux j 23 ). Inflow of Ca 2 + by neuronal stimulation lead to release of vesicular DA into the synaptic cleft. Clearance of DA in the synaptic cleft occurs by diffusion, glial cell metabolization [14,21], or uptake of DA by DAT [77,78].
Figure 5. Reaction scheme of the TH-DA negative feedback model with added rate constants and fluxes. Not all DOPA and DA are necessarily bound to the channeling complex, but since the complex is not evenly distributed within the cytoplasm, both cytosolic DOPA and DA are considered to be located proximate to the complex. With increasing age and increased MAO levels [74] MAO may oxidize DA and perturb the feedback control of DOPA (path outlined in gray). In addition, both DOPA and DA are subject to auto-oxidation [75,76], which are included in the model by k 14 (together with ROS) and α ( k 18 ,   k 19 ) (together with MAO). As in Figure 3, pairs of rate constants ( k i , k j ) represent the respective ( V m a x , K M ) values of applied Michaelis–Menten kinetics. With increased MAO concentrations, oxidative stress and ROS also increase due to MAO-generated hydrogen peroxide. ROS can oxidize DOPA [30] (outlined in red). Cytosolic DA concentration is increased due to the re-entry of DA (flux j 11 ) and by vesicular leakage [16,17,18,19,20] (flux j 23 ). Inflow of Ca 2 + by neuronal stimulation lead to release of vesicular DA into the synaptic cleft. Clearance of DA in the synaptic cleft occurs by diffusion, glial cell metabolization [14,21], or uptake of DA by DAT [77,78].
Ijms 22 12862 g005
Figure 6. The K M ( = k 8 ) of the transport of DA into vesicles has a major influence on the controller accuracy for DOPA homeostasis toward step-wise increased perturbations by oxidative stress ( k 14 ). (a) k 8 = 2.9 × 10 3   μ M; (b) k 8 = 2.9 × 10 2   μ M; (c) k 8 = 2.9 × 10 1   μ M; (d) 2.9 μ M. Other rate constants (ad): k 1 = 10 μ M/min, k 2 = 1 × 10 4 μ M/min, k 3 = 74.4 μ M, k 4 = 6 μ M, k 5 = 10.0 μ M/min, k 6 = 346 μ M, k 7 = 6.0 μ M/min, k 9 = 2.3 μ M/min, k 10 = 0.025 min 1 , k 11 = 1 μ M/min, k 12 = 50 μ M, k 13 = 437 μ M, k 14 phase 1: 0.0 min 1 , k 14 phase 2: 0.01 min 1 , k 14 phase 3: 0.1 min 1 , k 15 = 2.3 μ M/min, k 16 = 56 μ M, k 17 = 0.01 min 1 , k 18 = 90 μ M/min, k 19 = 111 μ M, k 20 = 0.1 μ M, k 21 = 1 × 10 4   μ M, k 22 = 0.0 min 1 , k 23 = 0.0158 min 1 . Initial concentrations, (ad): T y r = 323.2 μ M, D O P A = 52.9 μ M, D A = 512.5 μ M, D A v e s = 232.6 μ M, D A e x = 1.7 μ M. Framed panel c shows the calculations with the experimentally determined value [81] of k 8 .
Figure 6. The K M ( = k 8 ) of the transport of DA into vesicles has a major influence on the controller accuracy for DOPA homeostasis toward step-wise increased perturbations by oxidative stress ( k 14 ). (a) k 8 = 2.9 × 10 3   μ M; (b) k 8 = 2.9 × 10 2   μ M; (c) k 8 = 2.9 × 10 1   μ M; (d) 2.9 μ M. Other rate constants (ad): k 1 = 10 μ M/min, k 2 = 1 × 10 4 μ M/min, k 3 = 74.4 μ M, k 4 = 6 μ M, k 5 = 10.0 μ M/min, k 6 = 346 μ M, k 7 = 6.0 μ M/min, k 9 = 2.3 μ M/min, k 10 = 0.025 min 1 , k 11 = 1 μ M/min, k 12 = 50 μ M, k 13 = 437 μ M, k 14 phase 1: 0.0 min 1 , k 14 phase 2: 0.01 min 1 , k 14 phase 3: 0.1 min 1 , k 15 = 2.3 μ M/min, k 16 = 56 μ M, k 17 = 0.01 min 1 , k 18 = 90 μ M/min, k 19 = 111 μ M, k 20 = 0.1 μ M, k 21 = 1 × 10 4   μ M, k 22 = 0.0 min 1 , k 23 = 0.0158 min 1 . Initial concentrations, (ad): T y r = 323.2 μ M, D O P A = 52.9 μ M, D A = 512.5 μ M, D A v e s = 232.6 μ M, D A e x = 1.7 μ M. Framed panel c shows the calculations with the experimentally determined value [81] of k 8 .
Ijms 22 12862 g006
Figure 7. Increase in V m a x ( k 2 ) improves the accuracy of DOPA homeostasis but leads to a slower response of the controller. Panel (a) is a recalculation of Figure 6d. Panel (b) shows the controller’s response when k 2 has increased from k 2 = 1 × 10 4 μ to k 2 = 1 × 10 5 μ M/min. Other rate constants as in Figure 6d. Initial concentrations, panel (b): T y r = 323.1 μ M, D O P A = 53.0 μ M, D A = 5.3 × 10 3 μ M, D A v e s = 232.4 μ M, D A e x = 1.7 μ M. Note that the controller’s accuracy can be further improved by having k 8 values lower than 2.9 μ M.
Figure 7. Increase in V m a x ( k 2 ) improves the accuracy of DOPA homeostasis but leads to a slower response of the controller. Panel (a) is a recalculation of Figure 6d. Panel (b) shows the controller’s response when k 2 has increased from k 2 = 1 × 10 4 μ to k 2 = 1 × 10 5 μ M/min. Other rate constants as in Figure 6d. Initial concentrations, panel (b): T y r = 323.1 μ M, D O P A = 53.0 μ M, D A = 5.3 × 10 3 μ M, D A v e s = 232.4 μ M, D A e x = 1.7 μ M. Note that the controller’s accuracy can be further improved by having k 8 values lower than 2.9 μ M.
Ijms 22 12862 g007
Figure 8. Increase in k 1 leads to a poorer homeostatic regulation of DOPA due to TH inhibition by increased Tyr levels. The graph shows DOPA and DA concentrations when k 1 = 100 μ M/min. Other rate constant values as in Figure 6c. Initial concentrations: Tyr = 3.7 × 10 3   μ M, DOPA = 51.2 μ M, DA = 18.1 μ M, D A v e s = 228.9 μ M, D A e x = 1.6 μ M.
Figure 8. Increase in k 1 leads to a poorer homeostatic regulation of DOPA due to TH inhibition by increased Tyr levels. The graph shows DOPA and DA concentrations when k 1 = 100 μ M/min. Other rate constant values as in Figure 6c. Initial concentrations: Tyr = 3.7 × 10 3   μ M, DOPA = 51.2 μ M, DA = 18.1 μ M, D A v e s = 228.9 μ M, D A e x = 1.6 μ M.
Ijms 22 12862 g008
Figure 9. Stronger inhibition of TH by DA (lower k 4 ), for example by dephosphorylation [22], leads to more rapid response kinetics (in comparison with Figure 6). However, the stronger inhibition also leads to lower DA steady-state levels and to a larger offset of DOPA from its set-point. Rate constants as in Figure 6, except for k 4 , which was reduced by one order of magnitude to 0.6 μ M. Initial concentrations: (a) Tyr = 323.2 μ M, DOPA = 52.9 μ M, DA = 51.3 μ M, D A v e s = 232.5 μ M, D A e x = 1.7 μ M; (b) Tyr = 323.3 μ M, DOPA = 52.8 μ M, DA = 51.4 μ M, D A v e s = 232.4 μ M, D A e x = 1.7 μ M; (c) Tyr = 324.0 μ M, DOPA = 52.3 μ M, DA = 52.6 μ M, D A v e s = 231.3 μ M, D A e x = 1.6 μ M; (d) Tyr = 329.3 μ M, DOPA = 48.2 μ M, DA = 63.3 μ M, D A v e s = 222.4 μ M, D A e x = 1.5 μ M. Framed panel c shows the calculations with k 8 = 0.29 μ M [81].
Figure 9. Stronger inhibition of TH by DA (lower k 4 ), for example by dephosphorylation [22], leads to more rapid response kinetics (in comparison with Figure 6). However, the stronger inhibition also leads to lower DA steady-state levels and to a larger offset of DOPA from its set-point. Rate constants as in Figure 6, except for k 4 , which was reduced by one order of magnitude to 0.6 μ M. Initial concentrations: (a) Tyr = 323.2 μ M, DOPA = 52.9 μ M, DA = 51.3 μ M, D A v e s = 232.5 μ M, D A e x = 1.7 μ M; (b) Tyr = 323.3 μ M, DOPA = 52.8 μ M, DA = 51.4 μ M, D A v e s = 232.4 μ M, D A e x = 1.7 μ M; (c) Tyr = 324.0 μ M, DOPA = 52.3 μ M, DA = 52.6 μ M, D A v e s = 231.3 μ M, D A e x = 1.6 μ M; (d) Tyr = 329.3 μ M, DOPA = 48.2 μ M, DA = 63.3 μ M, D A v e s = 222.4 μ M, D A e x = 1.5 μ M. Framed panel c shows the calculations with k 8 = 0.29 μ M [81].
Ijms 22 12862 g009
Figure 10. DOPA addition can improve DOPA homeostasis when j 5 + j 11 + j 23 < j 7 . Panels a and b show reference calculations with two different k 8 values (representing different controller accuracies) when DOPA is not added ( k 22 = 0 μ M/min). Panels c and d show the effect of DOPA additions. For all panels, (ad): phase 1, no DOPA addition; phase 2, k 22 = 2 μ M/min addition; phase 3, k 22 = 4 μ M/min addition. Other rate constant values are as in Figure 9. Initial concentrations in panels a and c are as in Figure 9c. Initial concentrations in panels b and d are as in Figure 9b.
Figure 10. DOPA addition can improve DOPA homeostasis when j 5 + j 11 + j 23 < j 7 . Panels a and b show reference calculations with two different k 8 values (representing different controller accuracies) when DOPA is not added ( k 22 = 0 μ M/min). Panels c and d show the effect of DOPA additions. For all panels, (ad): phase 1, no DOPA addition; phase 2, k 22 = 2 μ M/min addition; phase 3, k 22 = 4 μ M/min addition. Other rate constant values are as in Figure 9. Initial concentrations in panels a and c are as in Figure 9c. Initial concentrations in panels b and d are as in Figure 9b.
Ijms 22 12862 g010
Figure 11. DA integral wind-up occurs when, here in phase 3, k 22 is set to 8 μ M/min, which results in the condition j 5 + j 11 + j 23 > j 7 . Other rate constant values and initial concentrations as for Figure 10c.
Figure 11. DA integral wind-up occurs when, here in phase 3, k 22 is set to 8 μ M/min, which results in the condition j 5 + j 11 + j 23 > j 7 . Other rate constant values and initial concentrations as for Figure 10c.
Ijms 22 12862 g011
Figure 12. Re-establishing robust DOPA homeostasis at reduced TH activity levels by DOPA addition. In both panels: phase 1 has a normal TH activity ( k 2 = 1 × 10 4   μ M/min) and no oxidative stress ( k 14 = 0.0 μ M/min); in phase 2, oxidative stress is present ( k 14 = 0.1 μ M/min), and TH activity level is reduced by two orders of magnitude ( k 2 = 1 × 10 2   μ M/min). As a result, DOPA homeostasis breaks down in phase 2 with low DOPA and DA levels. In phase 3, DOPA is added ( k 22 = 7 μ M/min) in the presence of low TH and oxidative stress as in phase 2. In phase 3, DOPA homeostasis is restored. Panel (a): k 8 = 0.29 μ M; panel (b): k 8 = 0.029 μ M. Other rate constant values and initial concentrations, as in Figure 10a.
Figure 12. Re-establishing robust DOPA homeostasis at reduced TH activity levels by DOPA addition. In both panels: phase 1 has a normal TH activity ( k 2 = 1 × 10 4   μ M/min) and no oxidative stress ( k 14 = 0.0 μ M/min); in phase 2, oxidative stress is present ( k 14 = 0.1 μ M/min), and TH activity level is reduced by two orders of magnitude ( k 2 = 1 × 10 2   μ M/min). As a result, DOPA homeostasis breaks down in phase 2 with low DOPA and DA levels. In phase 3, DOPA is added ( k 22 = 7 μ M/min) in the presence of low TH and oxidative stress as in phase 2. In phase 3, DOPA homeostasis is restored. Panel (a): k 8 = 0.29 μ M; panel (b): k 8 = 0.029 μ M. Other rate constant values and initial concentrations, as in Figure 10a.
Ijms 22 12862 g012
Figure 13. Calculated vesicular steady-state DA concentrations D A v e s as a function of cytosolic DA levels and DOPA controller accuracy k 8 (Equation (15)). At high controller accuracy (low k 8 , line no. 5), D A v e s levels are less affected when DA concentrations decrease. At low controller accuracy (high k 8 , line no. 1), D A v e s levels decrease rapidly with increased DOPA perturbations and with decreasing DA concentrations. k 8 values used: 1, 29 μ M; 2, 2.9 μ M; 3, 0.29 μ M; 4, 0.029 μ M; and 5, 0.0029 μ M.
Figure 13. Calculated vesicular steady-state DA concentrations D A v e s as a function of cytosolic DA levels and DOPA controller accuracy k 8 (Equation (15)). At high controller accuracy (low k 8 , line no. 5), D A v e s levels are less affected when DA concentrations decrease. At low controller accuracy (high k 8 , line no. 1), D A v e s levels decrease rapidly with increased DOPA perturbations and with decreasing DA concentrations. k 8 values used: 1, 29 μ M; 2, 2.9 μ M; 3, 0.29 μ M; 4, 0.029 μ M; and 5, 0.0029 μ M.
Ijms 22 12862 g013
Figure 14. Vesicular DA concentration as a function of perturbation k 14 and controller accuracy k 8 . The step-wise changes in k 14 for phases 2 and 3, as well as the other rate constants and initial concentrations, are the same as in Figure 6. Curves 1, 2, and 3 refer to the conditions of Figure 6a,c,d, respectively.
Figure 14. Vesicular DA concentration as a function of perturbation k 14 and controller accuracy k 8 . The step-wise changes in k 14 for phases 2 and 3, as well as the other rate constants and initial concentrations, are the same as in Figure 6. Curves 1, 2, and 3 refer to the conditions of Figure 6a,c,d, respectively.
Ijms 22 12862 g014
Figure 15. Deteriorated DOPA homeostasis when DA is subject to removal by MAO/auto-oxidation. Solid lines show the performance when α = 0 (Figure 6c). Dashed lines show the concentrations when α = 1 and DA is removed. Rate constant values as for Figure 6c with k 18 = 90 μ M/min and k 19 = 111 μ M. Initial concentrations (dashed curves, α = 1): Tyr = 210.7 μ M, DOPA = 518.3 μ M, DA = 1.705 μ M, D A v e s = 149.7 μ M, D A e x = 0.62 μ M. Red and blue lines show DOPA and DA concentrations, respectively.
Figure 15. Deteriorated DOPA homeostasis when DA is subject to removal by MAO/auto-oxidation. Solid lines show the performance when α = 0 (Figure 6c). Dashed lines show the concentrations when α = 1 and DA is removed. Rate constant values as for Figure 6c with k 18 = 90 μ M/min and k 19 = 111 μ M. Initial concentrations (dashed curves, α = 1): Tyr = 210.7 μ M, DOPA = 518.3 μ M, DA = 1.705 μ M, D A v e s = 149.7 μ M, D A e x = 0.62 μ M. Red and blue lines show DOPA and DA concentrations, respectively.
Ijms 22 12862 g015
Figure 16. Effect of cytosolic DA oxidation/auto-oxidation on vesicular DA levels. α = 0 indicates the results when using the conditions as in Figure 6c, i.e., no DA removal occurs. α = 1 allows oxidation of DA to occur. Same rate parameters as in Figure 6c are applied. Initial concentrations for α = 1 are the same as in Figure 15.
Figure 16. Effect of cytosolic DA oxidation/auto-oxidation on vesicular DA levels. α = 0 indicates the results when using the conditions as in Figure 6c, i.e., no DA removal occurs. α = 1 allows oxidation of DA to occur. Same rate parameters as in Figure 6c are applied. Initial concentrations for α = 1 are the same as in Figure 15.
Ijms 22 12862 g016
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kleppe, R.; Waheed, Q.; Ruoff, P. DOPA Homeostasis by Dopamine: A Control-Theoretic View. Int. J. Mol. Sci. 2021, 22, 12862. https://doi.org/10.3390/ijms222312862

AMA Style

Kleppe R, Waheed Q, Ruoff P. DOPA Homeostasis by Dopamine: A Control-Theoretic View. International Journal of Molecular Sciences. 2021; 22(23):12862. https://doi.org/10.3390/ijms222312862

Chicago/Turabian Style

Kleppe, Rune, Qaiser Waheed, and Peter Ruoff. 2021. "DOPA Homeostasis by Dopamine: A Control-Theoretic View" International Journal of Molecular Sciences 22, no. 23: 12862. https://doi.org/10.3390/ijms222312862

APA Style

Kleppe, R., Waheed, Q., & Ruoff, P. (2021). DOPA Homeostasis by Dopamine: A Control-Theoretic View. International Journal of Molecular Sciences, 22(23), 12862. https://doi.org/10.3390/ijms222312862

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop