ORIGINAL RESEARCH article

Front. Nucl. Med., 01 December 2022

Sec. PET and SPECT

Volume 2 - 2022 | https://doi.org/10.3389/fnume.2022.1030387

Transient equilibrium determination of dopamine D2/D3 receptor densities and affinities in brain

  • 1. Department of Neurology, Gødstrup Hospital, Herning, Denmark

  • 2. NIDO - Centre for Research and Education, Gødstrup Hospital, Herning, Denmark

  • 3. Department of Neurology, Aarhus University, Aarhus, Denmark

  • 4. Department of Nuclear Medicine and PET Centre, Aarhus University Hospital, Aarhus, Denmark

  • 5. Johns Hopkins Medical Institutions, Department of Radiology and Radiological Sciences, Division of Nuclear Medicine PET Center, Baltimore MD, United States

  • 6. School of Medicine, Washington University in St. Louis, St. Louis, MO, United States

  • 7. Institute of Regional Health Research, University of Southern Denmark, Odense, Denmark

  • 8. Medical Spinal Research Unit, Spine Centre of Southern Denmark, University Hospital of Southern Denmark, Odense, Denmark

  • 9. Department of Diagnostic Radiology and Nuclear Medicine, Saitama Medical Center, Saitama Medical University, Moroyama, Japan

  • 10. Translational Neuropsychiatry Unit, Department of Clinical Research, Aarhus University, Aarhus, Denmark

  • 11. Departments of Nuclear Medicine and Clinical Research, University of Southern Denmark and Odense University Hospital, Odense, Denmark

  • 12. Department of Neuroscience, University of Copenhagen, Copenhagen, Denmark

  • 13. Department of Neurology and Neurosurgery, McGill University, Montréal, QC, Canada

  • 14. Neuroscience Center, Tabriz University of Medical Sciences, Tabriz, Iran

Article metrics

View details

1

Citations

1,6k

Views

928

Downloads

Abstract

Long-term alteration of dopaminergic neurotransmission is known to modulate the D2/D3 receptor expression in the brain. The modulation can occur as a response to pathological processes or pharmacological intervention. The receptor density can be monitored by in vivo positron emission tomography (PET) of raclopride. To obtain accurate measurements of receptor-ligand interaction, it is essential to estimate binding parameters at true (if transient) equilibrium of bound and unbound ligand quantities. We designed this study as a comparison of two quantitative approaches to transient equilibrium, the TRansient EquilibriuM BoLus Estimation (TREMBLE) method and the Transient Equilibrium Model (TEM) method, to determine binding parameters at transient equilibrium with bolus injection of the radioligand. The data demonstrates that TREMBLE unlike TEM identified the time at which equilibrium existed. TREMBLE revealed that equilibrium prevailed at one or more times after bolus injection and identified differences of receptor density among regions such as putamen and caudate nucleus. We demonstrated that TREMBLE is a quantitative approach suitable for the study of pathophysiological conditions of certain types of neurotransmission the brain.

Highlights

  • 1.

    We invented a method of quantitation of radioligand binding at transient but true equilibrium of binding (TREMBLE).

  • 2.

    We compared the results of the new method of D2/D3 receptor density in striatum of humans with a standard method (TEM).

  • 3.

    The method of TREMBLE yielded superior equilibrium binding estimates of the radioligand raclopride compared to TEM.

Introduction

The D2/D3 dopaminergic receptor antagonist raclopride has been widely used to study the dopaminergic system by means of in vivo PET since 1985, when the tracer was introduced for the first time for human brain studies (1). The radioligand raclopride binds with high selectivity to / receptors, and the binding is inhibited by challenge from endogenous dopamine (DA) (2,3) and other dopamine receptor agonists and antagonists (4). Thus, raclopride satisfies pharmacological criteria of application to studies of dopaminergic neurotransmission and receptor occupancy by drugs.

Dopamine mediates its physiological action through five subtypes of G protein-coupled receptors, D1-D5 with subtypes. Particularly, the / receptors are highly expressed in brain areas critical for motor control, mesolimbic function, and memory processing, such as striatum, nucleus accumbens, olfactory tubercle, ventral tegmental area (VTA), and hippocampus [reviewed in Beaulieu and Gainetdinov (5)]. Therefore, dopaminergic signalling at these receptors is of special interest to neurological disorders with dysregulated dopaminergic transmission, including Parkinson’s disease (PD) and schizophrenia among others (6,7). Deficient dopaminergic transmission is one of the main underlying causes of motor impairment in PD, and post-mortem examinations have shown that the loss of dopamine is regionally heterogeneous, with greater loss in putamen than in caudate nucleus (8). Patients with PD display correspondingly greater / receptor binding in putamen than in caudate nucleus, suggesting compensatory upregulation in response to dopamine deficiency (9).

Contrary to the dopaminergic deficits observed in Parkinson’s disease, increased dopaminergic activity is held to be an underlying pathophysiological characteristic of schizophrenia, as shown both by single-photon emission computed tomography (SPECT) and PET studies, in which amphetamine administration evokes greater dopamine release in patients with schizophrenia than in healthy control subjects (3,1013). Interestingly, the magnitude of dopamine release in striatum evoked by amphetamine correlates significantly with changes of symptom severity on the Brief Psychiatric Rating Scale (BPRS) (3), and the rate of dopamine synthesis is higher in patients with schizophrenia than in healthy control subjects, as determined by means of PET with 3,4-dihydroxy-6-(18)F-fluoro-l-phenylalanine (FDOPA) (14).

To examine pathophysiology and effect of treatment by PET, it is necessary to estimate receptor density () and affinity () separately rather than binding potential. The necessity arises from the long-term modulation of dopaminergic neurotransmission that is known to lead to adaptive changes of the receptor system, as originally shown by early reports of autoradiography revealing increased receptor density in striatum by experimental administration of a selective antagonist (15,16).

The binding potential term originally introduced to PET quantification by Mintun et al. (17) as an equilibrium parameter reflects the ratio of bound to unbound radioligand that equals the product of the density of receptors without ligand and the affinity of the ligand (). Therefore, long-term modulation that alters the receptor density is reflected in the binding potential estimate. Of the several formulations of the binding potential, the term for the binding potential relative to non-displaceable binding in a region () was introduced by Innis et al. (18).

True receptor density and affinity must be determined at equilibrium with Eadie-Hofstee or Scatchard graphical analysis based on at least two levels of receptor occupancy. The gold standard of true albeit transient equilibrium in PET was achieved by continuous infusion of a radioligand that maintained a constant concentration in brain (19). However, continuous delivery of ligand experimentally is more difficult and any infusion conditions may not be more generalizable to all subjects in any given population than bolus injection approach. Therefore, quantitative kinetic models have been applied to bolus injection where the possible departure from the steady-state by continuous infusions are not an experimental limitation. The TRansient EquilibriuM BoLus Estimation (TREMBLE) method is an approach that yields the quantity of specifically bound radioligand on the assumption that the bound ligand reaches an equilibrium state at specific time points after bolus injection (20,21). The approach uses the radiolabelled metabolite-corrected plasma concentration of radioligand to specifically differentiate displaceable from non-displaceable receptor bound quantities of the radioligand, as derived in subsection I of the Appendix below.

To avoid the demand for arterial blood sampling, several approaches were presented in the past that apply an alternative reference region-derived input. One example is the Transient Equilibrium Model (TEM) that uses a reference region, e.g., the cerebellum, as an approximation of the tracer input from the activity in a region of non-displaceable binding (22), as derived in subsection II in the Appendix below. Both TREMBLE and TEM were designed to identify time points of transient equilibrium, as derived below in the Appendix. The aim of the present study was to compare the two approaches for quantification of receptor density using raclopride binding in healthy human brain.

Material and methods

Subjects

PET data of 21 healthy subjects were included in this study as a subset of a previous cohort (23). The criteria of selection of subjects for this study were based on the availability of data with blood sampling because the analysis with TREMBLE requires knowledge of tracer concentrations in arterial plasma. At the time of enrolment, all participants underwent physical examination, with no display of any abnormal neurological findings. The subjects had no neurological or psychiatric diseases in the past history, nor gave any evidence of substance abuse. All subjects signed informed consent forms prior to participation. The experiments were conducted at the Johns Hopkins Hospital in accordance with the Declaration of Helsinki as approved by the Institutional Review Board (IRB) of the Johns Hopkins Hospital.

Three subjects were excluded from the analysis because they displayed maximal receptor blockade at 75–97% in the challenge condition (as shown in subsection III in the Appendix with details of the basis of exclusion shown in Figure A4 and listed in Table A1). Complete or close-to-complete receptor blockade is not applicable to estimation of receptor density, because the Eadie-Hofstee linearisation requires evidence of receptor occupancy at a minimum of two different degrees of occupancy. Therefore, the remaining 18 subjects of the original 21 subjects had an average age of (). The eleven women and seven men included subjects of Caucasian (13), African-American (4), and Asian (1) ethnicity.

PET acquisition and analysis

All subjects underwent dual raclopride PET acquisitions at baseline followed by challenge with unlabeled raclopride. Molar activities were  GBq/mol () and  GBq/mol, respectively. Individual values are also listed in Table A2. Each acquisition lasted 90 min, during which we drew arterial blood samples to determine tracer input concentrations.

We analyzed the data by two kinetic models, described in subsections I and II in the Appendix, using a GUI custom-built in MATLAB (Mathworks), available at MATLAB central (24). Both models were designed to identify the instances of true but transient equilibrium of the bound quantity of tracer as the time(s) at which , e.g., at the peak of the binding curve (as indicated with arrows in Figures 1B,C). We automatically identified The time of steady-state by means of a function that pinpointed the numeric maximum of an array. A few subjects displayed a flat curve in the challenge condition of the TREMBLE analysis. In those cases, we determined the bound quantity as the mean value in the interval 20–60 min of the curve.

Figure 1

Statistics

To test if and estimates obtained by TREMBLE differed from the stimates obtained by TEM, we did one-way ANOVA, followed by Tukey’s correction of multiple comparisons. Paired t-tests were completed to test whether the estimates of the challenge condition differed from baseline. The applied statistical methods are specified in the respective figure legends. For all tests, we considered a -value of less than 0.05 to be indicative of significance. Statistical tests were performed by Graphpad Prism v. 7.00.

Compliance with ethical standards

The authors have no conflict of interest. Twenty-one healthy human participants were enrolled in the study, of whom 18 met the present criteria of moderate degrees of receptor occupation (i.e., 75% or less) by unlabeled raclopride. All subjects signed informed consent forms prior to participation. The experiments were conducted at Johns Hopkins Hospital in accordance with the Declaration of Helsinki, as approved by the Institutional Review Board (IRB) of the Johns Hopkins Hospital.

Results

The time-activity curves from a representative healthy subject demonstrate that accumulation of raclopride was markedly reduced upon challenge with unlabeled ligand in both putamen and caudate nucleus when compared to the baseline condition (Figure 1A). The magnitude and profile of the plasma input concentrations were similar in the two conditions, in support of the interpretation that the lower tracer binding in the challenge conditions in target regions of dopaminergic neurotransmission was due to competition from unlabeled ligand and not from differences of the arterial input concentrations. Cerebellum, known to be devoid of dopamine binding sites, had identical time-activity curves in the two conditions, indicating that no displacement of bound ligand occurred in the reference region.

The binding curves obtained by the TREMBLE and TEM analyses of this representative subject are shown in Figures 1B and 1C, respectively. For the TEM analysis, we used the activity of the reference region () as an approximation of the unbound and non-displaceable tracer quantity in brain. For the TREMBLE analysis, we computed the equivalent exchangeable quantity () from the measured arterial plasma input concentrations. The computed curve of initially was higher in the challenge condition than at baseline, suggesting a greater exchangeable quantity of raclopride as a result of displacement by unlabeled raclopride. This phenomenon was not observed with the TEM analysis.

The population average of the contents of respective compartments in putamen are shown in Figure 2. The population averaged curves consistently displayed the same time course pattern as curves of the single representative subject. In the challenge condition, at 0–20 min, the accumulation of ligand in the exchangeable compartment () exceeded that of the baseline condition, after which time the curves approached similar levels in both conditions (Figure 2C). We considered the initially higher quantity during the challenge condition a result of the lower receptor availability of receptors blocked by unlabeled ligand. This behavior of the exchangeable quantity was detected by the TREMBLE analysis but not by the TEM analysis.

Figure 2

Challenge with non-radiolabeled raclopride significantly reduced the estimates of by the two methods (Figure 3). The average reductions of the estimates averaged and () in putamen by TREMBLE and TEM, respectively. As the two models yielded similar reductions of , we asked whether a scalable difference existed between the results of the two methods. We compared the percentage declines in each subject by TREMBLE and TEM (Figure 3C) and found no consistent pattern, implying that the outcomes of TREMBLE and TEM yield random differences.

Figure 3

With the estimates of specific binding (, ) and binding potential () at two occupancy levels in the absence and presence of unlabelled raclopride, we applied Eadie-Hofstee plots to obtain the receptor density, , and half-saturation constant, , as presented in Figures 4A,B and subsection IV in the Appendix below. Comparison of receptor densities in Figure 4C revealed that TREMBLE detected a regional difference of with significantly higher density in putamen ( pmol/) than in caudate nucleus ( pmol/). In contrast, analysis with TEM yielded no difference of receptor density between the two parts of the striatum. Figure 4D shows the same values of in the two brain regions, regardless of analysis method.

Figure 4

We examined the dynamic time courses of the apparent estimates over time of both models, because the result of the Eadie-Hofstee plot reflects the magnitude of . Interestingly, TREMBLE revealed that estimates of reached a plateau at approximately 20 minutes after bolus injection when (Figures 5A and B). This finding is consistent with the hypothesis that peaks or troughs identified by TREMBLE mark the times of true transient equilibrium. The apparent constant magnitudes of reflect the times when the ratios of bound-to-free ligand quantities are sustained for longer intervals. In contrast, the apparent values of continuously rose for 40 minutes when quantified by TEM.

Figure 5

Because the estimations of and by TREMBLE is influenced by the magnitudes of two constants, and (Eqs. A1 and A5), we examined if the constants were affected by challenge with unlabeled raclopride. As shown in Figures 6A and B, the constants were not affected by the unlabelled raclopride.

Figure 6

Discussion

We determined / dopamine receptor densities in putamen and caudate nucleus, using two different kinetic models, TREMBLE and TEM, to identify the binding at transient equilibrium. Although, the averages of decreased significantly upon blocking with unlabeled raclopride according to both methods, the individual estimates of , and by the two models were not related.

Values of and obtained by the Eadie-Hofstee plot of TEM data match the results of previous PET studies that included TEM analysis (see Table 1 for a summary). In the present study, the TEM analysis yielded no significant difference of receptor density between the regions of putamen and caudate nucleus. In contrast, TREMBLE analysis of the same dataset showed a significantly greater receptor density of putamen than caudate nucleus.

Table 1

Brain regionKinetic% DiffAgeReference
model(pmol/)Mean (SD)Mean (SD)
Mean (SD)
Putamen18TEM29.3 (8.9)810.4 (2.9)34.9 (14.0)This study
Caudate nucleusTEM27.2 (9.4)10.6 (3.6)
PutamenTREMBLE26.7 (11.9)4812.2 (6.7)
Caudate nucleusTREMBLE18 (8.9)9.1 (6.7)
Putamen L20TEM27.7 (7.9)258.9 (2.5)27.5 (4.9)Farde et al. (25)
Caudate nucleus L22.1 (5.3)8.1 (2.0)
Putamen R27.6 (8.0)239.1 (2.3)
Caudate nucleus R22.5 (7.7)7.7 (2.4)
Putamen14TEM29.6 (6.8)1011.9 (3.6)59.2 (7.4)Rinne et al. (26)
Caudate nucleus26.9 (7.5)11.0 (3.2)
Striatum R8TEM28.7 (4.9)n/a8.2 (1.1)36.9 (6.4)Hietela et al. (27)
Striatum L28.4 (5.7)9 (2.2)
Striatum R (males)33TEM26.7 (7.6)n/a9.9 (2.6)40.2 (16.7)Pohjalainen et al. (28)
Striatum L (males)25.4 (7.2)9.4 (2.7)
Striatum R (females)21TEM27.2 (6.3)10.6 (2.0)
Striatum L (females)27.3 (7.9)11.0 (2.9)
Striatum (test)4TEM19.4 (2.3)n/a7.8 (0.9)28.0 (7.0)Hietela et al. (29)
Striatum (re-test, short term follow-up)21.8 (2.2)8.3 (0.8)
Striatum (test)4TEM28.7 (3.5)8.6 (1.2)
Striatum (re-test, long term follow-up)26.3 (2.3)8.2 (0.9)
Striatum3TREMBLE26 (4.0)n/a9.6 (4.0)n/aWong et al. (21)

Comparison of and of / receptors with other raclopride studies in healthy human subjects %Diff denotes the percentage difference in between putamen and caudate nucleus.

From the PET studies listed in Table 1, Farde et al. (25) and Rinne et al. (26) reported the / dopamine receptor density of putamen and caudate nucleus in healthy subjects obtained by application of TEM analysis and Scatchard plots. Farde et al. (25) reported 25% higher density in putamen than caudate nucleus, while the difference reported by Rinne et al. (26) was as low as 10%. In other studies of the receptor density with TEM analysis, the authors merged the subregions and reported densities in striatum that are comparable with the results of the present TEM analysis.

It was not clear from previous PET studies whether the density of / receptors is higher in putamen than in caudate nucleus. Therefore, we compared the present results with data from in vitro autoradiography in healthy humans, presented in (Table 2). To make the comparison with in vitro studies possible, we converted the density per unit wet weight to density per unit dry weight of protein. Under the assumption that brain tissue contains 10% solid material, conversion of wet to dry weight units uses division by 10 (36). The autoradiography studies listed in Table 2 support considerably higher density in putamen than in caudate nucleus. As the studies were performed in vitro and with other ligands of the / receptors, the results are not numerically translatable. However, the comparison of the relative regional difference shows that results from TREMBLE are consistent with data obtained by autoradiography.

Table 2

LigandBrain region (fmol/mg protein) Mean (SD)% DiffAge Mean (SD)Reference
PET[11C]raclopridePutamen18267 (119)2534.9 (14.0)This study
Caudate nucleus180 (89)
Autoradiography[3H]raclopridePutamen20249 (112)1477.2 (11.0)Piggott et al. (30)
Caudate nucleus218 (98)
[3H]spiroperidolPutamen4450 (6)750.2 (14.5)Joyce et al. (31)
Caudate nucleus420 (11)
Putamen10268 (28)37n/aCamps et al. (32)
Caudate nucleus195 (30)
[125I]epidepridePutamen4152 (25)20n/aJoyce et al. (33)
Caudate nucleus5127 (18)
Putamena24174 (19)3068.0 (14.0)Murray et al. (34)
Caudate nucleusa134 (21)
Putamen4100b (2.7)22n/aHall et al. (35)
Caudate nucleus82b (0.3)

Comparison of with autoradiography studies in healthy human brain.

The studies listed in the table consistently showed that putamen has markedly higher receptor density than caudate nucleus. %Diff denotes the percentage difference in between putamen and caudate nucleus.

aIndicates rostral part of striatum.

bIndicates the mean % of the density in putamen.

We also compared the values of from the TEM and TREMBLE analyses of this study with the results reported in the literature. The Eadie-Hofstee regression is driven by the magnitudes of and the bound quantity of tracer ( or , depending on method). Table 3 lists the publications that included values of obtained by linear regression methods, including the reference region version of the Logan plot (42) and the simplified reference tissue method (SRTM) (43). The publications listed in Table 3 consistently report greater values of in putamen than in caudate nucleus, with a difference of 20–30% that matches the results of TREMBLE. The lack of a difference between estimates of for the two striatal regions with TEM analysis can be attributed to the use of cerebellum activity in TEM as an approximation of non-specifically bound ligand in regions of binding. Instead, the correct estimate of unbound ligand in a region of specific binding depends on the number of receptors available. When receptors are blocked, the quantity of unbound ligand increases as the receptor availability declines. Thus, the use of cerebellum activity as an estimate of unbound and non-displaceable ligand in the computation of specific binding is erroneous. Ito and colleagues (44) simulated a study to characterise the error associated with the use of cerebellum activity in TEM instead of the real plasma input function. The authors demonstrated that TEM caused larger errors and led to overestimation of the values of when the assumed magnitudes of and both were low, and the authors also demonstrated mathematically that the concentration in cerebellum equals the free unbound and non-specifically bound ligand at the time of steady-state, but is not equal to zero and the bound ligand therefore is not at equilibrium at the time of steady-state identified by the TEM analysis.

Table 3

Brain regionKinetic model (pmol/) Mean (SD)% DiffAge Mean (SD)Reference
Putamen18TEM2.9 (0.3)434.9 (14.0)This study
Caudate nucleusTEM2.8 (0.3)
PutamenTREMBLE2.0 (0.2)25
Caudate nucleusTREMBLE1.6 (0.3)
Putamen (right)8Logan2.1 (0.3)2425.0 (5.8)Yoder et al. (37)
Caudate nucleus (right)1.7 (0.2)
Putamen (left)2.3 (0.3)41
Caudate nucleus (left)1.6 (0.3)
Putamen110Logan3.0 (0.4)30n/aZhou et al. (38)
Caudate nucleus2.3 (0.3)
Putamen10Logan2.8 (0.4)3332.2 (11.1)Black et al. (39)
Caudate nucleus2.1 (0.5)
Putamen12Logan2.6 (0.2)1363.0 (7.0)Politis et al. (40)
Caudate nucleus2.3 (0.2)
Putamen10SRTM3.0 (0.3)36n/aShotbolt et al. (41)
Caudate nucleus2.2 (0.3)

Comparison of of raclopride binding in healthy human subjects.

Logan indicates graphical anlysis using Logan plot with cerebellum as reference input. SRTM indicates simplified reference tissue model. %Diff denotes the percentage difference in between putamen and caudate nucleus.

Notably, the apparent values of computed with TREMBLE in the present study (Figures 5A and B) approached a constant value at the time when . The prolonged constant ratio between bound and free ligand quantities with TREMBLE analysis strongly supports the contention that TREMBLE identifies the ultimately transient equilibrium that persisted for some time after bolus injection. Determination of at equilibrium is crucial to obtain valid estimates of the binding potential that is defined only at instances of true equilibrium (17). Thus, identification of the time of equilibrium is critical to accurate and reproducible quantification of the values of and .

In conclusion, we demonstrated that the TREMBLE method yields valid estimates of receptor-ligand interaction at true (albeit transient) equilibrium after bolus injection of tracer. This conclusion would also apply to to experiments based on bolus-plus-infusion of the tracer, The reason is that the bolus-plus-infusion experiments focus on the establishment of constant concentrations of all tracer molecules in the tissue, rather than on the requirement of constant levels of bound tracer associated with true equilibrium, a goal that is very difficult to reach. A disadvantage of TREMBLE is the requirement for arterial plasma concentration sampling that makes clinical studies labor intensive, but the approach is important to the detection of moderate changes as found especially in neuropsychiatric disorders, for example.

Key points

QUESTION: A true (albeit often transient) equilibrium of the binding of radioligands to neuroreceptors must be present for binding parameters to be valid. To idnetify instances of true equilibrium, we designed a method of quantitation of radioligand binding that allowed us to identify instances of transient but true equilibrium of binding.

PERTINENT FINDINGS: We compared the results of the new method of identification of instances of transient equilibrium (TREMBLE) of / receptor binding of labeled raclopride in striatum of humans with the standard transient equilibrium method (TEM) that fails to identify such instances of true but transient equilibrium (TEM).

IMPLICATIONS FOR PATIENT CARE: The method of TREMBLE allows measurements of parameters of receptor binding to reach a level of accuracy that assures more precise and therefore relevant estimates of binding potentials on which to base the therapeutic engagement of dopamine receptors during treatment of pathological conditions involving dopaminergic neurotranismission of human brain.

Statements

Data availability statement

The datasets presented in this article are not readily available because the data is a few years old and has only the TAC curves and not the original images. Requests to access the datasets should be directed to Dean F. Wong and Albert Gjedde .

Ethics statement

The studies involving human participants were reviewed and approved by the Johns Hopkins Medicine Investigational Review Board. The patients/participants provided their written informed consent to participate in this study.

Author contributions

All authors contributed to the data analysis and review of the paper. DFW acquired the data and AG led the idea of the algorithm with DFW. All authors contributed to the article and approved the submitted version.

Funding

This work was supported in part by Public Health Service Grant RO1 NIH MH42821 (DFW and AG). J-AP completed the work in the course of an MD-PhD fellowship from Aarhus University.

Acknowledgments

We are grateful for the administrative and technical support we received from Johns Hopkins University. Special thanks to Ayon Nandi and Andrew Crabb for excellent technical support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1.

    FardeLEhrinEErikssonLGreitzTHallHHedströmC, et al. Substituted benzamides as ligands for visualization of dopamine receptor binding in the human brain by positron emission tomography. Proc Natl Acad Sci. (1985) 82:38637. 10.1073/pnas.82.11.3863

  • 2.

    EndresCJKolachanaBSSaundersRCSuTWeinbergerDBreierA, et al. Kinetic modeling of [] raclopride: combined pet-microdialysis studies. J Cereb Blood Flow Metab. (1997) 17:93242. 10.1097/00004647-199709000-00002

  • 3.

    BreierASuTPSaundersRCarsonRKolachanaBDe BartolomeisA, et al. Schizophrenia is associated with elevated amphetamine-induced synaptic dopamine concentrations: evidence from a novel positron emission tomography method. Proc Natl Acad Sci. (1997) 94:256974. 10.1073/pnas.94.6.2569

  • 4.

    TauscherJHussainTAgidOVerhoeffNPLWilsonAAHouleS, et al. Equivalent occupancy of dopamine D1, D2 receptors with clozapine: differentiation from other atypical antipsychotics. Am J Psychiatry. (2004) 161:16205. 10.1176/appi.ajp.161.9.1620

  • 5.

    BeaulieuJMGainetdinovRR. The physiology, signaling, pharmacology of dopamine receptors. Pharmacol Rev. (2011) 63(1):182217.

  • 6.

    ElsingaPHHatanoKIshiwataK. Pet tracers for imaging of the dopaminergic system. Curr Med Chem. (2006) 13:213953. 10.2174/092986706777935258

  • 7.

    FazioPSchainMMrzljakLAminiNNagSAl-TawilN, et al. Patterns of age related changes for phosphodiesterase type-10A in comparison with dopamine receptors, sub-cortical volumes in the human basal ganglia: A pet study with -MNI-659, -raclopride with correction for partial volume effect. Neuroimage. (2017) 152:3309. 10.1016/j.neuroimage.2017.02.047

  • 8.

    KishSJShannakKHornykiewiczO. Uneven pattern of dopamine loss in the striatum of patients with idiopathic Parkinson’s disease. N Engl J Med. (1988) 318:87680. 10.1056/NEJM198804073181402

  • 9.

    AntoniniASchwarzJOertelWHPogarellOLeendersKL. Long-term changes of striatal dopamine receptors in patients with Parkinson’s disease: a study with positron emission tomography and [] raclopride. Mov Disord. (1997) 12:338. 10.1002/mds.870120107

  • 10.

    McCutcheonRAAbi-DarghamAHowesOD. Schizophrenia, dopamine and the striatum: from biology to symptoms. Trends Neurosci. (2019) 42(3):20520.

  • 11.

    Abi-DarghamAvan de GiessenESlifsteinMKegelesLSLaruelleM. Baseline and amphetamine-stimulated dopamine activity are related in drug-naive schizophrenic subjects. Biol Psychiatry. (2009) 65:10913. 10.1016/j.biopsych.2008.12.007

  • 12.

    LaruelleMAbi-DarghamAVan DyckCHGilRD’SouzaCDErdosJ, et al. Single photon emission computerized tomography imaging of amphetamine-induced dopamine release in drug-free schizophrenic subjects. Proc Natl Acad Sci. (1996) 93:923540. 10.1073/pnas.93.17.9235

  • 13.

    LaruelleMAbi-DarghamAGilRKegelesLInnisR. Increased dopamine transmission in schizophrenia: relationship to illness phases. Biol Psychiatry. (1999) 46:5672. 10.1016/S0006-3223(99)00067-0

  • 14.

    ReithJBenkelfatCSherwinAYasuharaYKuwabaraHAndermannF, et al. Elevated dopa decarboxylase activity in living brain of patients with psychosis. Proc Natl Acad Sci. (1994) 91:116514. 10.1073/pnas.91.24.11651

  • 15.

    McGoniglePBoysonSJReuterSMolinoffPB. Effects of chronic treatment with selective and nonselective antagonists on the subtypes of dopamine receptors. Synapse. (1989) 3:7482. 10.1002/syn.890030111

  • 16.

    SubramaniamSLuckiIMcGonigleP. Effects of chronic treatment with selective agonists on the subtypes of dopamine receptors. Brain Res. (1992) 571:31322. 10.1016/0006-8993(92)90670-5

  • 17.

    MintunMARaichleMEKilbournMRWootenGFWelchMJ. A quantitative model for the in vivo assessment of drug binding sites with positron emission tomography. Ann Neurol. (1984) 15:21727. 10.1002/ana.410150302

  • 18.

    InnisRBCunninghamVJDelforgeJFujitaMGjeddeAGunnRN, et al. Consensus nomenclature for in vivo imaging of reversibly binding radioligands. J Cereb Blood Flow Metab. (2007) 27:15339. 10.1038/sj.jcbfm.9600493

  • 19.

    PatlakCSPettigrewKD. A method to obtain infusion schedules for prescribed blood concentration time courses. J Appl Physiol. (1976) 40:45863. 10.1152/jappl.1976.40.3.458

  • 20.

    SøllingTBrustPCunninghamVWongDGjeddeA. True equilibrium bolus estimation (tremble) confirms rapid transient equilibrium. Neuroimage. (1997) 5:A29.

  • 21.

    WongDSøllingTYokoiFGjeddeA. Quantification of extracellular dopamine release in schizophrenia and cocaine use by means of tremble. Quantitative functional brain imaging with positron emission tomography. Academic Press, Published by Elsevier Inc. (1998). p. 463–8.

  • 22.

    FardeLErikssonLBlomquistGHalldinC. Kinetic analysis of central [] raclopride binding to -dopamine receptors studied by pet—a comparison to the equilibrium analysis. J Cereb Blood Flow Metab. (1989) 9:696708. 10.1038/jcbfm.1989.98

  • 23.

    KuwabaraHMcCaulMEWandGSEarleyCJAllenRPWeertsEM, et al. Dissociative changes in the , of dopamine / receptors with aging observed in functional subdivisions of the striatum: a revisit with an improved data analysis method. J Nucl Med. (2012) 53:80512. 10.2967/jnumed.111.098186

  • 24.

    [Dataset] PhanJA. Tremble, TEM analysis (2019).

  • 25.

    FardeLWieselFAStone-ElanderSHalldinCNordströmALHallH, et al. dopamine receptors in neuroleptic-naive schizophrenic patients: a positron emission tomography study with [] raclopride. Arch Gen Psychiatry. (1990) 47:2139. 10.1001/archpsyc.1990.01810150013003

  • 26.

    RinneJOLaihinenARuottinenHRuotsalainenUNågrenKLehikoinenP, et al. Increased density of dopamine receptors in the putamen, but not in the caudate nucleus in early Parkinson’s disease: a pet study with [] raclopride. J Neurol Sci. (1995) 132:15661. 10.1016/0022-510X(95)00137-Q

  • 27.

    HietalaJWestCSyvälahtiENågrenKLehikoinenPSonninenP, et al. Striatal dopamine receptor binding characteristics in vivo in patients with alcohol dependence. Psychopharmacology. (1994) 116:28590. 10.1007/BF02245330

  • 28.

    PohjalainenTRinneJONågrenKSyvälahtiEHietalaJ. Sex differences in the striatal dopamine receptor binding characteristics in vivo. Am J Psychiatry. (1998) 155:76873. 10.1176/ajp.155.6.76

  • 29.

    HietalaJNågrenKLehikoinenPRuotsalainenUSyvälahtiJ. Measurement of striatal dopamine receptor density and affinity with []-raclopride in vivo: a test-retest analysis. J Cereb Blood Flow Metab. (1999) 19:2107. 10.1097/00004647-199902000-00012

  • 30.

    PiggottMMarshallEThomasNLloydSCourtJJarosE, et al. Dopaminergic activities in the human striatum: rostrocaudal gradients of uptake sites and of and but not of receptor binding or dopamine. Neuroscience. (1999) 90:43345. 10.1016/S0306-4522(98)00465-5

  • 31.

    JoyceJNSappDWMarshallJF. Human striatal dopamine receptors are organized in compartments. Proc Natl Acad Sci. (1986) 83:80026. 10.1073/pnas.83.20.8002

  • 32.

    CampsMCortesRGueyeBProbstAPalaciosJ. Dopamine receptors in human brain: autoradiographic distribution of sites. Neuroscience. (1989) 28:27590. 10.1016/0306-4522(89)90179-6

  • 33.

    JoyceJJanowskyANeveK. Characterization and distribution of [125I] epidepride binding to dopamine receptors in basal ganglia and cortex of human brain. J Pharmacol Exp Ther. (1991) 257:125363.

  • 34.

    MurrayAMRyooHLGurevichEJoyceJN. Localization of dopamine receptors to mesolimbic and receptors to mesostriatal regions of human forebrain. Proc Natl Acad Sci. (1994) 91:112715. 10.1073/pnas.91.23.11271

  • 35.

    HallHFardeLHalldinCHurdYLPauliSSedvallG. Autoradiographic localization of extrastriatal -dopamine receptors in the human brain using [125I] epidepride. Synapse. (1996) 23:11523. 10.1002/(SICI)1098-2396(199606)23:2%3C115::AID-SYN7%3E3.0.CO;2-C

  • 36.

    EricssonCPeredoINistérM. Optimized protein extraction from cryopreserved brain tissue samples. Acta Oncol. (2007) 46:1020. 10.1080/02841860600847061

  • 37.

    YoderKKKarekenDAMorrisED. What were they thinking?: Cognitive states may influence [] raclopride binding potential in the striatum. Neurosci Lett. (2008) 430:3842. 10.1016/j.neulet.2007.10.017

  • 38.

    ZhouYYeWBrašićJRCrabbAHHiltonJWongDF. A consistent and efficient graphical analysis method to improve the quantification of reversible tracer binding in radioligand receptor dynamic pet studies. Neuroimage. (2009) 44:66170. 10.1016/j.neuroimage.2008.09.021

  • 39.

    BlackKJPiccirilloMLKollerJMHseihTWangLMintunMA. Levodopa effects on [] raclopride binding in the resting human brain. F1000Research. (2015) 4. 10.12688/f1000research.5672.1

  • 40.

    PolitisMWilsonHWuKBrooksDJPicciniP. Chronic exposure to dopamine agonists affects the integrity of striatal d2 receptors in Parkinson’s patients. NeuroImage: Clin. (2017) 16:45560. 10.1016/j.nicl.2017.08.013

  • 41.

    ShotboltPTziortziACSearleGEColasantiAVan Der AartJAbanadesS, et al. Within-subject comparison of []-()-phno and [] raclopride sensitivity to acute amphetamine challenge in healthy humans. J Cereb Blood Flow Metab. (2012) 32:12736. 10.1038/jcbfm.2011.115

  • 42.

    LoganJFowlerJSVolkowNDWolfAPDeweySLSchlyerDJ, et al. Graphical analysis of reversible radioligand binding from time—activity measurements applied to [N--methyl]-()-cocaine pet studies in human subjects. J Cereb Blood Flow Metab. (1990) 10:7407. 10.1038/jcbfm.1990.127

  • 43.

    LammertsmaAAHumeSP. Simplified reference tissue model for pet receptor studies. Neuroimage. (1996) 4:1538. 10.1006/nimg.1996.0066

  • 44.

    ItoHHietalaJBlomqvistGHalldinCFardeL. Comparison of the transient equilibrium and continuous infusion method for quantitative pet analysis of [] raclopride binding. J Cereb Blood Flow Metab. (1998) 18:94150. 10.1097/00004647-199809000-00003

  • 45.

    GjeddeAWongDFRosa-NetoPCummingP. Mapping neuroreceptors at work: on the definition and interpretation of binding potentials after 20 years of progress. Int Rev Neurobiol. (2005) 63:120. 10.1016/S0074-7742(05)63001-2

  • 46.

    GjeddeAWongDF. Four decades of mapping and quantifying neuroreceptors at work in vivo by positron emission tomography. Front Neurosci. (2022) 16. 10.3389/fnins.2022.943512

  • 47.

    GjeddeA. Calculation of cerebral glucose phosphorylation from brain uptake of glucose analogs in vivo: a re-examination. Brain Res Rev. (1982) 4:23774. 10.1016/0165-0173(82)90018-2

  • 48.

    PatlakCSBlasbergRGFenstermacherJD. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab. (1983) 3:17. 10.1038/jcbfm.1983.1

  • 49.

    KhodaiiJAraj-KhodaeiMVafaeeMSWongDFGjeddeA. Relative strengths of three linearizations of receptor availability: Saturation, inhibition, and occupancy plots. J Nucl Med. (2022) 63:294301. 10.2967/jnumed.117.204453

  • 50.

    PhanJALandauAMJakobsenSGjeddeA. Radioligand binding analysis of adrenoceptors with [] yohimbine in brain in vivo: Extended inhibition plot correction for plasma protein binding. Sci Rep. (2017) 7:15979. 10.1038/s41598-017-16020-1

  • 51.

    GjeddeAWongD. Receptor occupancy in absence of reference region. Neuroimage. (2000) 11:S48.

  • 52.

    PhanJALandauAMWongDFJakobsenSNahimiADoudetDJ, et al. Quantification of [] yohimbine binding to adrenoceptors in rat brain in vivo. J Cereb Blood Flow Metab. (2015) 35:50111. 10.1038/jcbfm.2014.225

  • 53.

    LassenNBartensteinPLammertsmaAPrevettMTurtonDLuthraS, et al. Benzodiazepine receptor quantification in vivo in humans using [] flumazenil and pet: application of the steady-state principle. J Cereb Blood Flow Metab. (1995) 15:15265. 10.1038/jcbfm.1995.17

  • 54.

    CunninghamVJRabinerEASlifsteinMLaruelleMGunnRN. Measuring drug occupancy in the absence of a reference region: the Lassen plot re-visited. J Cereb Blood Flow Metab. (2010) 30:4650. 10.1038/jcbfm.2009.190

  • 55.

    EadieG. On the evaluation of the constants Vm and KM in enzyme reactions. Science. (1952) 116:688. 10.1126/science.116.3025.688.a

Appendix: Theories of transient equilibrium approach

Clarification on chosen nomenclature for binding potential equations

The nomenclature used here is based on equations reported by Gjedde et al. (45), Innis et al. (18), and Gjedde & Wong (46). However, the first-order differential equations reported by Gjedde et al. (45) and Gjedde & Wong (46) are not included in the formalism presented by Innis et al. (18) that focused on steady-state solutions. The first-order differential equations are shown below as the equivalent basis for the transient solutions on which TREMBLE depends. The steady-state solutions reported by Innis et al. (18) are equivalent to the steady-state solutions reported by Gjedde et al. (45) and will be presented as such at the appropriate place in the Appendix. The main distinction concerns the difference between the units of mass per unit volume used by Gjedde et al. (45) and Gjedde & Wong (46), and the units of concentration used by Innis et al. (18) that refer to spaces of tissue rather than volumes of fluid. It is the contention of Phan et al. authors that the units of mass per unit volume are more appropriate. In the retrospect we will suggest that the Innis et al (18) be modified in future extensions of nomenclature of Innis et al. (18).

I. TRansient EquilibriuM BoLus Estimation

The TRansient EquilibriuM Bolus Estimation (TREMBLE) method estimates the quantity of bound ligand to receptors when an equilibrium of bound and unbound tracer quantities is assumed to occur transiently. The method was introduced by Sølling et al. (20) and Wong et al. (21). As illustrated in Figure A1, the kinetic model describes two simultaneous dynamics, the tracer transfer from the plasma compartment () to the exchangeable tissue compartment () and to a separate compartment of specifically bound tracer (). The first order differential equation describes the distribution in the vascular compartment,where is the mass of ligand detectable by the tomograph and denotes the volume of the vascular bed, and is the time-variable concentration in plasma. The exchange from plasma to the exchangeable compartment can be described as follows,and to the binding compartment as follows,where denotes the plasma clearance, is the rate constant of exchange from tissue to plasma, and and denote the respective quantities in the exchangeable and specific bound compartments.

Figure A1

Multilinear analysis can be applied to fit multiple unknown variables, but with three or more compartments the number of unknown parameters and degrees of freedom become excessive to the extent that the analysis fails to yield accurate rate constants. To avoid the uncertainty of many variables, the fundamental aim of the analysis is to identify the quantity in an intermediate compartment that controls the exchange of ligands with subsequent compartments. Only in the initial phase of an tomography session is the quantity of ligand bound to receptors is negligible, and only then is the flux of ligand is controlled by the concentration gradient between plasma and brain tissue. At all time, the measured activity in a region of interest obeys the equation,such that rearrangement of Eq. A4 yields the quantity of ligand in the exchangeable compartment,where is the partition volume that equals . To solve Eq. A5, it is necessary to estimate , and the derivatives of and . The estimation of and using the Gjedde-Patlak plot and estimation of in a reference region with the two-compartment model is shown in separate sections below. When the quantity of ligand in the exchangeable compartment is known, the bound quantity is obtained by subtraction,where at transient equilibrium, the first derivative of the bound quantity, , equals zero. Innis et al. (18) correctly cite Mintun et al. (17) for the definition of binding potential as “the equilibrium ratio of specifically bound ligand (B) to its free concentration (F),” which we here interpret as the quantities of bound and unbound ligand at equilibrium, as expressed in Eq. A7 below. The time of the transient equilibrium corresponds to a peak of the curve of bound ligand, and at this time the steady-state binding potential is defined as the ratio of bound to free ligand, indicated by upper case symbols, and ,

Estimation of from a reference region

Assuming that the volume of distribution, equal to the ratio of the rate constants and , is identical in a reference region and in the regions of specific binding, the ratio of to in cerebellum can be obtained from the radioactivity in plasma and cerebellum.

When the specific binding of the ligand is negligible in cerebellum, as assumed, the kinetic model of three compartments is reduced to two compartments, as illustrated in Figure A2. The operational equation is then,where refers to the magnitude of in cerebellum. Then, the total quantity of ligand in the reference region, is obtained by integration of the differential Eq. A8,

Figure A2

Because is unknown, is expressed in terms of , and substitution of with the expression yields,

At later times, the approaching equilibrium state is associated with substantial washout of ligand from plasma, and the ligand accumulates in brain tissue and other organs. At this time, the magnitudes of and the ratio become negligible. Thus, Eq. A10 reduces to,where division of Eq. A11 with yields a simple linearized solution as previously described by Gjedde et al. (47),We resolved the cerebellar rate constants and graphically from the y-intercept and slope, respectively, and is simply the the ratio of the cerebellar and .

Estimation of the vascular volume

Regional and values can be estimated by means of the two-compartment model applied to the estimation of . At the onset of the PET acquisition, the tracer predominantly exists in the vascular compartment, and the loss of tracer from brain to plasma can be ignored,where is the observed radioactivity in a region of interest. When divided by , Eq. A13 yields the equation of Gjedde-Patlak plot (47,48),where is solved from the y-intercept, and the value of is obtained from the initial acquisition frames. In this study, we solved regional and at 0-2 minutes post-injection.

II. Transient equilibrium model

The ligand and neuroreceptor interaction described by Farde and colleagues (22) as the Transient Equilibrium Model (TEM) method (22) also involves a three-compartment model and four first order rate constants (Figure A3). In contrast to TREMBLE, however, where the center compartment represents the quantity of exchangeable ligand, the center compartment of TEM represents unbound or “free” quantity of the tracer denoted as .

Figure A3

The exchange of ligand in the compartments can be described by two differential equations,andwhere denotes the free quantity in tissue and is the concentration in plasma. For convenience, the specifically bound quantity in tissue alleged by TEM is denoted with an asterisk (). Assuming that specific binding of raclopride is negligible in the reference region, the exchange of ligand in the reference region is reduced to,andwhere by rearrangement of Eq. A18, the specifically bound quantity yields,In TEM, the alleged specifically bound quantity is obtained by subtraction of the quantity in the reference region () from the total quantity (). In this calculation, it is assumed that the free quantities in tissue and in plasma have equal magnitudes in and . As in TREMBLE, quantity of specifically bound ligand then is obtained in TEM at the time of transient equilibrium when the first derivative equals zero. The equilibrium binding potential in TEM is determined from the ratio of the bound quantity in a region of binding and the total quantity of tracer in the reference region,where it is important to note that the binding potential obtained by TREMBLE in principle therefore is different from the binding potential obtained by TEM because assumptions of non-displaceable quantities have different origins and are determined differently.

III. Inhibition plots of receptor saturation

We used Inhibition Plot procedures (49) (see Figure A4 below) to determine the degree of competition between labeled and unlabeled tracer molecules, to rule out degrees of competition that would be inconsistent with active exchange of labeled ligand molecules, because of complete or close-to-complete saturationffo receptors (Figure A4 and Table A1).

Figure A4

Table A1

Ex 1Ex 2Ex3
ROI
Putamen1.250.580.671.150.420.731.2250.3670.858
Putamen left1.230.580.651.210.450.761.2110.3390.872
Putamen right1.220.580.641.100.390.711.2390.3880.851
Caudate nucleus0.950.530.421.000.320.681.1090.2720.837
Caudate nucleus left0.840.510.330.930.290.641.1170.2950.822
Caudate nucleus right1.080.560.521.040.360.681.0650.2430.822
Cerebellum0.340.350.010.340.350.010.3210.3170.004

Regional values of the three excluded subjects.

The values in the three excluded subjects are listed. and denote the distribution volume at baseline and inhibition (challenge) conditions, whereas indicates the difference in between the two conditions. The estimates are obtained from graphical analysis as previously described by Phan et al. (50).

IV. Eadie-Hofstee plots of receptor density and affinity

We jointly determined receptor density and affinity by graphical means of the Eadie-Hofstee plot (55), based on binding parameters obtained at two different levels of receptor occupancy. The steady-state estimates of were plotted as a function of the bound quantity () (Figure 4),where the ordinate intercept yielded the receptor density , and the slope yielded the half-saturation concentration, equal to the reciprocal of the affinity . Innis et al. (18) focus on the number of free receptor sites () that equals the the difference between the total number of sites and the sites holding bound ligand(s). It is possible to rephrase the equation with the term by substitution, such that the Eadie-Hofstee equation at equilibrium now readsorthat allows the calculation of the Michaelis constant with unit of concentration from values of , , and .

V. Radiochemistry Details in Low Molar Activity Scans

In low molar activity scans, nonradioactive raclopride in sterile saline solution was added to the raclopride sterile saline solution. The solution was prepared just prior to injection into the subject. The final molar activity adjusted to the injection adjusted to the injection time was used for the calculation of and . The molar activity at baseline (high molar activity) was always higher than the challenge (low molar activity). The individual values are listed in Table A2.

Table A2

SubjectHigh molar activity (HMA)Low molar activity (LMA)HMA/LMA
116.090.86135
411.010.61678
44.570.6569
924.870.731259
1480.970.682181
271.770.75364
706.780.431651
Ex 135.850.4286
140.690.74189
272.270.72377
Ex 232.200.7145
533.360.61873
1424.430.951498
75.580.8490
53.980.5696
90.960.72127
Ex 383.740.63134
76.170.48158
228.120.61373
282.890.80355
254.320.60425
Mean359.080.67
SD432.140.14

Individual molar activity.

The values of individual high and low molar activities (HMA and LMA) are listed in units of GBq/mol. The values in the three excluded subjects are indicated as Ex 1,2 and 3.

Summary

Keywords

receptor density, dopamine receptor, positron emission tomography, raclopride, bolus injection

Citation

Phan J-A, Wong DF, Chang NHS, Kumakura Y, Bauer WR and Gjedde A (2022) Transient equilibrium determination of dopamine D2/D3 receptor densities and affinities in brain. Front. Nucl. Med. 2:1030387. doi: 10.3389/fnume.2022.1030387

Received

28 August 2022

Accepted

27 October 2022

Published

01 December 2022

Volume

2 - 2022

Edited by

Nicolas A. Karakatsanis, Cornell University, United States

Reviewed by

Anton Lindberg, University of Toronto, Canada, Sridhar Goud Nerella, National Institutes of Health (NIH), United States

Updates

Copyright

*Correspondence: Dean F. Wong Albert Gjedde

These authors have contributed equally to this work and share first authorship.

Specialty Section: This article was submitted to PET and SPECT, a section of the journal Frontiers in Nuclear Medicine

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics