Next Article in Journal
Prevalence and Characteristics of Metabolic Hyperferritinemia in a Population-Based Central-European Cohort
Next Article in Special Issue
The Effects of Chronic Psychostimulant Administration on Bone Health: A Review
Previous Article in Journal
Structural and Biomechanical Properties of Supramolecular Nanofiber-Based Hydrogels in Biomedicine
Previous Article in Special Issue
Prefrontal Dopamine in Flexible Adaptation to Environmental Changes: A Game for Two Players
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modulation of D3R Splicing, Signaling, and Expression by D1R through PKA→PTB Phosphorylation

by
Orlando Casados-Delgado
1,
José Arturo Avalos-Fuentes
1,
Manuel Lara-Lozano
1,
Gisela Tovar-Medina
1,
Carla Daniela Florán-Hernández
1,
Karla Gisela Martínez-Nolasco
1,
Hernán Cortes
2,
Ricardo Felix
3,
José Segovia
1 and
Benjamín Florán
1,*
1
Departamento de Fisiología, Biofísica y Neurociencias, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Mexico City 07360, Mexico
2
Laboratorio de Medicina Genómica, Departamento de Genómica, Instituto Nacional de Rehabilitación Luis Guillermo Ibarra Ibarra, Mexico City 14389, Mexico
3
Departamento de Biología Celular, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Mexico City 07360, Mexico
*
Author to whom correspondence should be addressed.
Biomedicines 2024, 12(1), 206; https://doi.org/10.3390/biomedicines12010206
Submission received: 22 December 2023 / Revised: 11 January 2024 / Accepted: 11 January 2024 / Published: 17 January 2024
(This article belongs to the Special Issue Dopamine Signaling Pathway in Health and Disease)

Abstract

:
The D1R and D3R receptors functionally and synergistically interact in striatonigral neurons. Dopaminergic denervation turns this interaction antagonistic, which is correlated with a decrement in D3nf isoform and an increment in D3R membranal expression. The mechanisms of such changes in D3R are attributed to the dysregulation of the expression of their isoforms. The cause and mechanism of this phenomenon remain unknown. Dopaminergic denervation produces a decrement in D1R and PKA activity; we propose that the lack of phosphorylation of PTB (regulator of alternative splicing) by PKA produces the dysregulation of D3R splicing and changes D3R functionality. By using in silico analysis, we found that D3R mRNA has motifs for PTB binding and, by RIP, co-precipitates with PTB. Moreover, D1R activation via PKA promotes PTB phosphorylation. Acute and 5-day D1R blockade decreases the expression of D3nf mRNA. The 5-day treatment reduces D3R, D3nf, and PTB protein in the cytoplasm and increases D3R in the membrane and PTB in the nucleus. Finally, the blockade of D1R mimics the effect of dopaminergic denervation in D1R and D3R signaling. Thus, our data indicate that through PKA→PTB, D1R modulates D3R splicing, expression, and signaling, which are altered during D1R blockade or the lack of stimulation in dopaminergic denervation.

Graphical Abstract

1. Introduction

The co-activation of D3 receptor (D3R) and dopamine D1 receptor (D1R) results in two possible signaling pathways. In typical signaling, the D3R Gi-mediated effects oppose D1R Gs/olf effects [1,2,3,4]; in atypical signaling, D3R synergizes and potentiates D1R effects via a dimeric interaction mediated by Gs/olf ([1,5,6,7,8,9]. In the nigral projections from the dorsal striatum, dopaminergic denervation promotes the apparition of D3R typical signaling that coexists with the normal atypical version [10]. The D3R typical signal (in which D3R opposes the D1R stimulation of cAMP accumulation and GABA release) masks the normal atypical signaling, with D3R potentiating both parameters. This phenomenon was related to an increase in the expression of D3R in the membrane [10]. The dopamine D3R isoform D3nf regulates the canonical D3R expression in the membrane [11,12,13,14]. During dopaminergic denervation, decreased D3nf isoform expression and protein–protein interaction with D3R allows D3R to be inserted into the membrane and signal [10,15]. These findings indicate the altered regulation of D3R membranal location and its mRNA alternative splicing during dopaminergic denervation, which can explain the appearance of functional D3R that opposes D1R. However, to date, the mechanism that generates these changes remains unknown.
Polypyrimidine tract-binding protein, PTB, is a splicing regulator [16]. PTB regulates dopamine D2 receptors (D2R) for alternative splicing in heterologous expression systems, decreasing D2Rshort isoform expression [17]. On the other hand, the nuclear location and function of PTB in splicing are regulated by phosphorylation from PKA [18]. When PKA phosphorylates PTB, it is exported from the nucleus, decreasing its function in splicing [18]. With these antecedents in the case of D1R and D3R interaction in striatonigral neurons, it can be proposed that the lack of D1R activation (activator of PKA) during denervation decreased the phosphorylation of PTB that remains in the nucleus, repressing splicing. This may produce a reduced expression of the D3nf isoform that lets D3R canonical move to the membrane and signal, producing the typical response. Our proposal assumes that PTB also regulates D3R splicing. Here, we study this proposal.
By using in silico analysis, we found that D3R mRNA has motifs for PTB binding and, by RIP, co-precipitates with PTB. Moreover, D1R activation via PKA promotes PTB phosphorylation. Acute and 5-day blockades of D1R decrease the expression of D3nf mRNA, which correlates with a decrease in D3R, D3nf, and PTB protein in the cytoplasm and an increment in D3R in the membrane and PTB in the nucleus. Finally, the blockade of D1R mimics the effect of dopaminergic denervation, producing the apparition of the D3R-mediated typical response. Thus, our data indicate that D1R modulates through PKA→PTB, D3R splicing, expression, and signaling, which is altered during D1R blockade or the lack of stimulation through dopaminergic denervation.

2. Materials and Methods

2.1. In Silico Analysis

We evaluated the presence of splicing regulatory motifs for the deleted sequence of 98 nucleotides that were previously reported for D3nf production in D3R mRNA using in silico analysis [11,19]. We used the Human Splicing Finder program, version 3.1. We searched for exon and intron regions, GC motifs for 5’ and GA for 3’ splicing sites, branch points (CURAY sequence), and ISS and ESS motifs on polypyrimidine tracts (TY). Delimited sequences were the sixth exon and the two flanking introns with 300 up and down nucleotides that were analyzed. We compared the literature reports for PTB motifs in different splicing sites with respect to the obtained results from our analysis to test whether PTB is a potential regulator of D3R splicing at this exon.

2.2. Animals

Male Wistar rats (200–250 g) housed together (five per cage) with food and water ad libitum and kept under a natural light cycle were used throughout. All procedures followed the National Institutes of Health Guide for Care and Use of Laboratory Animals. They were approved by the Institutional Animal Care Committee of the CINVESTAV, making all efforts to minimize animal suffering.

2.3. SCH23390 Treatment

Rats were injected in two ways; according to the experiment, some rats received a single dose of 0.5 mg/kg i.p. of D1R antagonist SCH 23390 and were sacrificed at 6 or 12 h after injection. Others received 0.5 mg/kg of SCH23390 or saline twice daily (every 12 h) for 5 days and, on the 6th day, were sacrificed for experiments. Rats for mRNA immunoprecipitation and PTB phosphorylation assay did not receive SCH 23390 treatment.

2.4. Whole Striatum Homogenates and Brain Slice Preparation

After rapid decapitation, the brains were removed and immersed in ice-cold artificial cerebrospinal fluid (aCSF) using the following composition (mM): 118.25 NaCl, 1.75 KCl, 1 MgSO4, 1.25 KH2PO4, 25 NaHCO3, 2 CaCl2, and 10 D-glucose. A rapid dissection of the whole striatum was performed, then the tissue was immersed in aCSF and processed for WB experiments. The cAMP accumulation and PTB phosphorylation experiments were performed in brain slices. The brain was glued to a metal cube mounted on a Petri dish filled with ice-cold aCSF, and brain slices (300 µm thick) containing the striatum were obtained with a vibroslicer (Campden Inc., Cambridge, UK). The slices were transferred to cold slides, and under a stereoscopic microscope, the striatum was microdissected and used in different experiments. The atlas of Paxinos and Watson (2006) [20] was utilized to identify individual nuclei. The slices or whole striatum tissue obtained from the two striata of each rat were pooled and used separately in each experiment so that n = 1 represents the data obtained from the two striatum of one rat.

2.5. Cellular Fractionation Procedures

For the WBs of D3R, D3nf, and PTB in subcellular fractions enriched with the nucleus, cytoplasm, and cell membranes, we used the method modified from Heydorn [21,22]. The whole striatal nuclei of three to five male Wistar rats (8 to 12 weeks old, weighing approximately 220 g) were used. Once dissected, the tissue was weighed and homogenized in seven volumes (w/v) of 0.32 M sucrose/5 mM HEPES pH 7.4 with protease inhibitors. The homogenate was centrifuged at 800× g for 10 min at 4 °C. The resulting precipitate (P1) is the crude nuclear fraction, resuspended in 0.32 M sucrose/5 mM HEPES pH 7.4 with protease inhibitors added; this was stored at −20/−70 °C until its use. The supernatant (S1) was centrifuged at 9200× g for 20 min at 4 °C. The precipitate 2 (P2) and supernatant 2 (S2) were obtained. S2 was centrifuged at 105,000× g for 90 min at 4 °C, and the supernatant 3 (S3) corresponding to the cytoplasmic fraction was obtained and stored. The precipitate corresponding to the microsomes was discarded. P2 was resuspended for lysis in I0 volumes of 5 mM HEPES pH 7.4 on ice for 30 min. Then, it was centrifuged at 25,000× g for 20 min at 4 °C, producing precipitate 3 (P3), and the supernatant containing the post-synaptic vesicles was discarded. P3 was resuspended in a small volume of 0.32 M sucrose/5 mM HEPES pH 7.4 and carefully placed on a stepwise gradient of 0.32, 0.9, and 1.2 M sucrose, respectively, which was centrifuged at 83,000× g for 120 min at 4 °C. Two bands were obtained from the previous step, and the interface between the 0.32 and 0.9 sucrose gradient, corresponding to myelin, was discarded. The second interface, between the 0.9 and 1.2 M sucrose gradient, corresponds to the membrane fraction, which was carefully extracted, added with protease inhibitors, and stored for later use.
Cellular fractionation was carried out according to [23,24] for RIP procedures to obtain the enriched nuclear fraction. In brief, striatal tissue was homogenized with 8 mL of cold lysis buffer (HEPES 20 mM, KCl 125 mM, MgCl2 4 mM, NP40 0.05%, one tablet of cOMPLETE (Cat. # 11873580001, Roche, Basel, Switzerland, Sigma-Aldrich, St. Louis, MO, USA) previously solved in 50 mL buffer, 40 μL of RNaseOUT (Cat # 10777-019 Invitrogene. Thermofisher Scientific, Waltham, MA, USA) and 80 μL of Ditiotreitol (DTT), pH 8). The homogenate was stored for 24 h at −80 °C to promote cell lysis. Then, it was centrifuged at 100× g for 5 min at 4 °C to remove cell debris. Then, another centrifugation was performed at 900× g for 10 min at 4 °C to separate the cytoplasmic and nuclear fractions. The pellet corresponding to the nuclear fraction was washed with 1 mL of sucrose buffer (Tris-HCl 20 mM, NaCl 60 mM, KCl 15 mM, Sucrose 0.34 M, two tablets of Complete previously solved in 20 mL buffer, 5 μL of RNaseOUT, 10 μL of DTT, pH 7.65). The pellet was resuspended with 2 mL of sucrose buffer with 10 μL of RNaseOUT, 20 μL of DTT, and 800 μL of high-salt buffer (Tris-HCl 20 mM, EDTA 0.2mM, Glycerol 25%, NaCl 900 mM, MgCl2 1.5 mM, half of a tablet of cOMPLETE previously solved in 5 mL buffer, 4 μL RNaseOUT and 8 μL DTT, pH 7.65) to promote nuclear membrane rupture; this was incubated in ice for 30 min, reversing every 5 min. For DNA digestion, 4.6 mL of sucrose buffer was added with 25 μL of RNaseOUT and 46 μL of DTT, supplemented with 7.5 μL of CaCl2 1M. Next, 75 μL of DNAsase (Cat # 18068015 DNAsase I, Invitrogen. Thermofisher Scientific) was added and incubated for 15 min at 37 °C, inverting the sample every 5 min. After this time, 60 μL of 0.5M EDTA was added. For soluble nuclear fraction collection, the sample was centrifuged at 16,000× g for 20 min at 4 °C, and the supernatant was transferred to a new cold tube.

2.6. PTB Phosphorylation Assay

In vitro phosphorylation assay was performed as described by Loya López (2020). Striatum (300 μM thick slices) were incubated for 30 min with the D1R agonist SKF 38393 (1 μM), the D1R antagonist SCH 23390 (100 nM), or PKA blocker H89 (10 μM), or a combination as indicated. After incubation, the slices were homogenized and centrifugated to obtain nuclei and cytoplasm, purified according to [25], and then separated by electrophoresis in 10% acrylamide gels. Gels were fixed in 50% methanol and 10% acetic acid overnight. Then, they were washed with ultrapure water three times and stained with Pro-Q Diamond phosphoprotein gel staining (Invitrogen) for 90 min. Gels were de-stained with a solution containing 20% acetonitrile, 50 mM sodium acetate, pH 4.0, 30 min, three times, and, finally, the phosphoprotein bands were detected using a GelDoc EZ from BioRad. The resulting bands were analyzed by using densitometry using ImageJ software, Version 1.46 (NIH, Bethesda, MD, USA).

2.7. mRNA Immunoprecipitation (RIP)

First, we prepared magnetic beads with PTB antibodies. A total of 100 μL of lysis buffer with 130 μL of magnetic beads and 5 μg of antibody against PTB were mixed. Then, they were incubated for 1 h at 4 °C to promote binding and were centrifuged at 5000× g for 15 s. Six washes were performed with 1 mL of lysis buffer to resuspend the beads in 900 μL of cold lysis buffer with 4.5 μL RNaseOUT and 9 μL DTT.
The supernatant of the nuclear fraction was transferred to the magnetic beads already attached to the antibody and was left to incubate in rotation at 4 °C overnight. The next day, four washes were performed with 500 μL cold lysis buffer, 2.5 μL RNaseOUT, and 5 μL DTT. The beads were resuspended with 100 μL lysis buffer, 0.5 μL of RNaseOUT, and 1 μL DTT. In order to elute the RNA bound to the beads, the protein was digested by adding 1.5 μL of Proteinase K solution (Cat 4333793 Invitrogene Thermofisher Scientific) 20 mg/mL) and incubated for 30 min at 55 °C in constant agitation. Finally, the supernatant was collected.
Total RNA was extracted by the TRIzol method. A total of 1 mL of TRIzol (Invitrogen, 15596026) was added to the supernatant collected from RNA immunoprecipitation. The sample was incubated for 5 min at room temperature, and then 200 μL of chloroform was added. The sample was shaken by inversion, incubated for 3 min at room temperature, and centrifuged for 15 min at 12,000 rpm at 4 °C. The aqueous phase containing the RNA was transferred to a new tube. Next, 500 μL of 100% isopropanol was added and incubated for 24 h at −20 °C. The next day, the sample was centrifuged for 10 min at 12,000 rpm at 4 °C, and the supernatant was removed. Subsequently, the RNA precipitate was washed two times with 1 mL of 75% ethanol. The supernatant was discarded, and the RNA pellet was dried for 10 min. The pellet was resuspended in 20 μL of RNase-free water (Gibco, Waltham, MA, USA) and incubated in a bath at 60 °C for 10 min.
Finally, cDNA was synthesized. The reaction was prepared with 2 μg of RNA, 1 μL of Random Primers (Random Primers Cat. 48190011 Invitrogene. Thermofisher Scientific), and 2 μL of 10 mM dNTPs. The reaction was incubated at 65 °C for 5 min and then placed on ice. Subsequently, 4 μL of 5X FS buffer, 2 μL of DTT 100 mM, and 1 μL of RNaseOUT were added to leave the sample incubating for 2 min at 37 °C. Then, 1 μL of the enzyme M-MLV RT (Thermofisher, 28025013) was added and incubated for 10 min at 25 °C. It was then incubated at 37 °C for 50 min. Finally, the sample was incubated at 70 °C for 15 min to inactivate the enzyme and was stored on ice.

2.8. Western Blot

For the Western Blot (WB) assay, 50 µg of protein (to detect PTB, D3nf, and D3R from coprecipitate samples) from subcellular fractionation were resolved by SDS-PAGE, transferred onto nitrocellulose membranes, and blotted for 2 h at room temperature in Tris-buffered saline containing 0.1% Tween 20 and 10% nonfat powdered milk. Membranes were incubated overnight at 4 °C with either antibody against PTB (1:1000), D3R (1:1000), or D3nf (1:1000). Polyclonal antibodies for D3R and D3nf were obtained from Santa Cruz Biotechnology Inc. numbers: D3DR (H-50): sc-9114, and D3nf (C-16): sc-9184. PTB antibody was purchased from Invitrogen, Thermo Fisher Scientific Inc., PTBP1 monoclonal antibody 32-4800. Antibodies were detected by chemiluminescence (ECL Plus Amersham) with HRP-conjugated secondary antibodies (1:20,000 dilution for antibodies against primary antibodies.

2.9. qRT-PCR Determinations

Following the treatment periods, the animals were euthanized, and both striatal regions were carefully dissected from brain slices, as described in Section 2.4. The dissected striatal tissues were collected in tubes containing 200 µL of TRIzol on ice. Subsequently, the tissues were macerated to achieve homogenization, and TRIzol was added until a final volume of 1 mL was reached. Total RNA extraction was then performed, as previously mentioned in Section 2.7. After total RNA purification, the sample’s purity and total RNA integrity were assessed before cDNA synthesis.
For tissue cDNA synthesis, 5 µg of total RNA was used in a final volume of 8 µL. This RNA was treated with 1 µL of DNase I and 1 µL of DNase I Buffer (10X). The reaction was incubated for 15 min at room temperature, and it was stopped by a heat shock and the addition of 1 µL of EDTA. In order to the previous reaction, 1 µL of Oligo dT (500 µg/µL) was added, followed by incubation at 65 °C for 5 min and 1 min on ice. The RT protocol was carried out using the parameters reported in Section 2.7.
For the qRT-PCR reactions from RIP cDNA and tissue cDNA, pre-designed probes from Thermofisher were used to evaluate the expression of the D3R gene (Rn00567568_m1) and D3nf gene (AIT96RG), and for the tissue samples, the Hprt1 gene (Rn01527840_m) was employed as an endogenous control [10]. The reaction mixtures contained 5 µL of 2x Master Mix (Applied Biosystems, Waltham, MA, USA), 2 µL of the sample, 0.5 µL of the respective probe, and 2.5 µL of water. The thermocycling conditions were as follows: an initial denaturation cycle at 94 °C for 5 min, followed by 55 cycles with two steps, 94 °C for 15 s, and 60 °C for 60 s. These qPCR reactions were carried out using the Bio-Rad CFX96 real-time equipment. Data were analyzed by using 2−ΔCt and 2−ΔΔCt calculations [26]. The fold change was obtained by comparing the SCH 23390-treated rats vs. the saline group and the ratio D3nf/D3R + D3nf to evaluate the relative expression of the isoforms in the saline vs. SCH 23390-treated condition.

2.10. cAMP Accumulation

We carried out [3H]cAMP accumulation assays according to the method described by Rangel-Barajas et al. [27]. Striatal slices obtained from one rat were pooled and incubated with [3H]-Adenine (130 nM) for 1 h at 37 °C. After this, the radioactive samples was washed thrice and then suspended in aSCF with 3-isobutyl-1-methylxanthine (IBMX) 1 mM. Incubation with dopamine receptor drugs was continued for 15 min and stopped by adding 100 μL of ice-cold trichloroacetic acid (15%) containing non-labeled ATP (2.5 mM) and cAMP (4.5 mM). N-ethylene maleimide (NEM) [10] was incubated with dopaminergic drugs to block Gi protein signaling. After 20 min on ice, the supernatant was loaded onto Dowex 50W-X4 (300 μL per column). A fraction containing [3H]-ATP was eluted with 3 mL of distilled water. A second eluent obtained with 5 mL of distilled water was directly loaded onto neutral alumina columns. Alumina columns were finally eluted with 4 mL of 50 mM Tris–HCl buffer (pH 7.4) to obtain [3H]-cAMP. ATP and cAMP eluents were transferred to vials, and radioactivity was determined by scintillation counting. The results were expressed as the ratio of [3H]- cAMP × 100/[3H]-cAMP + [3H]-ATP and then as a percent of the control condition.

2.11. Drugs

3-isobutyl-1-methylxanthine, IBMX; R(+)-7-Chloro-8-hydroxy-3-methyl-1-phenyl-2,3,4,5-tetrahydro-1H-3-benzazepine hydrochloride, SCH 23390; (+)-(4aR,10bR)-3,4,4a,10b-Tetrahydro-4-propyl-2H,5H-[1]benzopyrano [4,3-b]-1,4-oxazin-9-ol hydrochloride, PD 128,907; (±)-1-Phenyl-2,3,4,5-tetrahydro-(1H)-3-benzazepine-7,8-diol hydrochloride, SKF38393; were purchase from Merk KGaA, Mexico. 4′-Acetyl-N-[4-[4-(2-methoxyphenyl)-1-piperazinyl] butyl]-[1,1′-biphenyl]-4-carboxamide, GR103691 was obtained from Tocris Bioscience, CTR Scientific, Mexico.
Radiochemical: Adenine, [2,8-3H], >97%, 1 mCi (37MBq), [3H] Adenine; Aminobutyric Acid (GABA) γ-[2,3-3H(N)], Specific Activity: 25–40 Ci (925 GBq-1.48 TBq)/mmoL, 1 mCi (37 MBq) was purchased from Perkin Elmer (Waltham, MA, USA).

2.12. Data Analysis

Data are expressed as mean ± SEM. All data were analyzed using Graph Pad Software (San Diego, CA, USA), version 9.4.1. For the differences in three or more experimental conditions running in parallel, a one-way ANOVA combined with Tukey’s test was used to compare the differences in the same experiment. An unpaired t-test was employed to compare the two experimental conditions from different experiments or only two experimental conditions from the same experiment.

3. Results

3.1. In Silico Motifs for Splicing in the mRNA of D3R and Co-Immunoprecipitation with PTB

Silico analysis of the sixth exon of D3R mRNA showed that the 98 deleted nucleotides in D3nf are characterized as an intronic region. The GC motif 5’ position with a 78.48 score and GA motif 3’ position with a 73.04 score in the 98-nucleotide region are donor and acceptor regions; moreover, we found two branch points scores of 68.88 and 86.16, and finally, a critical ISS in the 26 position on a TY tract that has a score of 80.18. Furthermore, multiple motifs with high silencer activity were found (see Supplementary Material).
In order to assess whether the PTB protein binds to D3R mRNA, we conducted an RIP (RNA immunoprecipitation) procedure. After the RIP procedure, the presence of the D3R transcript was evaluated through qRT-PCR (see Section 2). The amplification of the D3R transcript indicates that D3R mRNA co-precipitates with PTB (Figure 1A). In order to confirm that the results of the qRT-PCR were not a nonspecific product, we evaluated the molecular weight of qRT-PCR products by agarose gel electrophoresis. We used cDNA from striatum tissue as a control. As observed, the molecular size of the amplicons from the RIP or striatum tissue presents the same molecular weight, corroborating the findings (Figure 1B).

3.2. D1R Activation Phosphorylates PTB

In order to test the phosphorylation of PTB promoted by D1R activation via PKA, we incubated striatal slices for 15 min with the D1R agonist SKF 38393 (1μM) [28] or SKF 38393 in the presence of the D1R selective receptor antagonist SCH 23390 (100 nM) [29] or the PKA blocker H89 (10μM) [30], and we determined the PTB phosphorylation level. Figure 2B,D show that D1R stimulation increased the basal phosphorylation levels of PTB by almost two times when compared to the basal state (Basal 1 vs. SKF38393 2.98 ± 0.32, F = 1.74 (2,6), p < 0.001, n = 3 determinations, ANOVA, followed by Tukey’s test). The addition of the antagonist SCH 23390 prevents this effect (Basal 1 vs. SKF38393+ SCH23390 1.30 ± 0.07, F = 1.74 (2,6), p = 0.52 ns, n = 3 determinations, ANOVA, followed by Tukey’s test). Similarly, the blockade of PKA with H89 prevents the SKF 38393 phosphorylation of PTB (Basal 1 vs. SKF38393+ H89 1.28 ± 0.10, F = 1.01 (2,6), p = 0.63 ns, n = 3 determinations, ANOVA, followed by Tukey’s test). Figure 2A,C show the representative blots from the three experiments.

3.3. D1R Blockade Decreases D3nf mRNA, Protein Expression, and Cytoplasmic Location

We tested if the lack of stimulation of D1R changes D3nf expression in a similar way to what is observed in dopaminergic denervation. First, in normal rats, we blockaded D1R via a single injection of SCH 23390 0.5 mgr./kg i.p.; after 6 and 12 h later, we determined the striatal expression of D3R and D3nf mRNA via qRT-PCR. Figure 3A shows a decrement of nearly 65% in the 2-ΔΔCT-fold change in D3nf mRNA 6 h post-injection, with a tendency to recover to saline control values 12 h post-injection (p = 0.014 at 6 h; p = 0.785 at 12 h; n = 6, ANOVA, followed by Tukey’s test); no significant changes in 2-ΔΔCT-fold change for D3R mRNA were found (Figure 3A). When evaluating the relative expression of D3nf mRNA compared to D3+D3nf in treated and saline control samples, a significant decrement was observed at 6 h. However, the mRNA levels tended to recover 12 h post-treatment (Figure 3B: ratio control 0.68 ± 0.08 vs. SCH 23390-treated, 6 h, 0.20 ± 0.07, p = 0.0007 vs. 12 h, 0.49 ± 0.06, ns, n = 6; ANOVA, followed by Tukey’s test).
The effect of 5 days of D1R blockade from a twice-a-day injection of D1R selective antagonist SCH 23390 on the mRNA and protein expression of D3nf in the whole striatum and the membranal and cytoplasmic subcellular fractions is shown in Figure 3C,D and Figure 4. In Figure 3C, it can be observed that D1R blockade decreases the expression of the mRNA of D3nf by nearly 75% with respect to the saline condition since the 2−ΔΔCT-fold change decreases to 0.24 ± 0.10 (p = 0.005, n = 4, t-test). In Figure 3D, it can be observed that the percent ratio of D3nf/D3nf+D3R also decreases significantly in the treated group (ratio control 0.67 ± 0.18 vs. SCH 23390-treated 0.12 ± 0.05, p = 0.021, t = 3.3, n = 4, and unpaired t-test). As expected from the mRNA decrement observed, an important decrement in D3nf protein expression was observed in the whole striatal tissue (Figure 4A,B) (D3nf relative expression with respect to actine: control 1.10 ± 0.11 vs. SCH23390-treated 0.34 ± 0.10, p = 0.003, t = 4.81, n = 4, and unpaired t-test). The decrement in protein expression occurs mainly in the cytoplasm, as can be observed in Figure 4E,F (D3nf relative expression with respect to the actine control 1.57 ± 0.36 vs. SCH23390-treated 0.41 ± 0.07, p = 0.041, t = 2.95, n = 3, and unpaired t-test); but this is not the case in the protein located in the cell membrane (Figure 4E,F; D3nf relative expression respect ATPase control 0.912 ± 0.18 vs. SCH23390-treated 0.79 ± 0.08, p = 0.92, t = 0.54, n = 3. Unpaired t-test). Figure 4A,C,E represent typical Western blots.

3.4. D1R Blockade Did Not Modify D3R mRNA but Increased Its Membrane Location

We also studied if D1R blockade affects the expression of D3R canonical isoform. In Figure 3A, it can be observed that no statistically significant change in D3R mRNA expression occurs in the striatum of the SCH 23390-treated rats since a 2-ΔΔCT-fold change concerning the saline condition did not differ from the expected value of 1 (1.31 ± 0.37 vs. 1, p = 0.46, n = 4; one sample t-test). Additionally, SCH 23390 did not affect D3R protein levels when evaluating the whole striatal protein content for both the monomers (45 kDa) and dimers (75 kDa) of the receptor [10,31,32] (Figure 5B: D3R 45 kDa relative expression with respect to actine control 1.03 ± 0.08 vs. SCH23390-treated 1.023 ± 0.23, p = 0.97, t = 0.04 and D3R 75 kDa relative expression control 1.036 ± 0.09 vs. SCH 23390-treated 1.027 ± 0.22, p = 0.97, ns, n = 3. Unpaired t-test). Interestingly, D3R protein expression increases in the membrane and decreases in the cytoplasm (Figure 5D: D3R 45 kDa relative expression with respect to ATPase, control 0.27 ± 0.04 vs. SCH23390-treated 1.027 ± 0.16, p = 0.01, t = 4.55 and D3R 75 kDa relative expression control 0.49 ± 0.12 vs. SCH 23390-treated 1.09 ± 0.02, p = 0.008; Figure 5F: D3R 45 kDa relative expression respect actine control 1.06 ± 0.11 vs. SCH23390-treated 0.396 ± 0.033, p = 0.005, t = 5.59 and D3R 75 kDa relative expression control 1.43 ± 0.11 vs. SCH 23390-treated 0.357 ± 0.09, p = 0.001 n = 3. Unpaired t-test). Figure 5A,C,E represent typical Western blots.

3.5. D1R Blockade Increased PTB Expression in the Nucleus

The effects on PTB protein expression via D1R blockade are shown in Figure 6. First, SCH 23390 treatment did not modify the protein expression in whole striatum tissue (Figure 6B: PTB relative expression with respect to actine control 1.64 ± 0.08 vs. SCH 23390-treated 1.64 ± 0.23, p = 0.98, t = 0.03, n = 3. Unpaired t-test). Nevertheless, it decreases protein in the cytoplasm and increases it in the nuclear fraction (Figure 6D: PTB relative expression with respect to actine control in the cytoplasm 1.64 ± 0.25 vs. SCH 23390-treated 0.59 ± 0.16, p = 0.025, t = 3.50, n = 3; Figure 6F: PTB relative expression with respect to actine control in nucleus 0.62 ± 0.06 vs. SCH 23390-treated 1.13 ± 0.06, p = 0.004, t = 5.70, n = 3. Unpaired t-test). Figure 6A,C,E represent representative Western blots.

3.6. Blockade of D1R Produces D3R Typical Signaling and Masks the Atypical One

Finally, we studied if D1R blockade produces the same changes in D3R signaling as dopaminergic denervation [8,10]. We studied the stimulation of cAMP accumulation via D1R and D1R+D3R co-activation in the striatal slices of the saline and SCH 23390-treated rats. D1R was stimulated with SKF38393 (1 μM) and D3R, with a PD of 128,907 (100 nM) [33]. In Figure 7A, the atypical response of D3R activation on the D1R stimulation of cAMP accumulation can be observed (potentiation), whereas the D3R agonist itself did not have an effect (Figure 7A: cAMP accumulation SKF38393 124 ± 3% vs. SKF38393+PD128,907 160 ± 5%, p = 0.0014, n = 6). In the SCH 23390-treated group, D3R activation antagonized the D1R stimulation of cAMP accumulation, which signifies the atypical response (Figure 7B: cAMP accumulation SKF38393 144 ± 1.85% vs. SKF38393+PD128,907 103 ± 3.83%, p < 0.0001, n = 3).
In order to test if the D3R atypical response masks the typical one, we blockaded Gi protein signaling by performing experiments in the presence of NEM (100 μM) [34]. In Figure 7C, it can be observed that in the saline-treated rats, the atypical response is preserved, even in the presence of NEM. On the contrary, in the SCH 23390-treated rats, the same atypical response as in the saline condition is observed (see Figure 7D).

4. Discussion

Several patho-physiological processes are related to changes in the splicing of D3R, which modify its signal and location [35]. Here, for the first time, we found evidence of the role of D1R, PKA, and PTB in regulating D3R splicing in striatal neurons. Our data indicate that D1R regulates D3R alternative splicing through PTB phosphorylation via PKA. They also suggest that PTB is a splicing regulator of D3R mRNA that promotes the expression of the D3nf isoform, which, as was previously suggested, regulates the membranal expression of the canonical receptor [14]. Finally, this mechanism impacts the functional effect of D1R and D3R interaction in striatal neurons.

4.1. PTB Modulates D3R mRNA Alternative Splicing and Produces D3nf Expression

Initially, the mRNA 98-nucleotide sequence deleted in D3nf mRNA was characterized in the human gene [11,36,37], but some differences in the mRNA of Rattus Norvegicus can be observed in the intron, based on the ensemble genome browser (see supplementary data), which indicates only an 87% analogy between both. In order to test the presence of PTB motifs in the rat sequence, we performed an in silico analysis of the reported rat sequence using Human Splicing Finder 3.1 software [11] (see Supplementary Material). In our analysis and agreement with the Schmauss group [11,36,37], the GC and GA sequences found are potential 3’ and 5’ splicing sites; GC is the more frequently atypic donor motif, and the GA site is considered a non-canonic acceptor site with a frequency of 1:19 [38]. In addition, our analysis shows that the deleted sequence has two potential ramification points with a high score. Thus, all in silico analysis data support the idea that the sequence is an atypical intronic region in the mRNA of the rat, as was reported for humans. Sironi’s motifs analysis of human deleted mRNA sequence indicates 15 silencer sequences with scores of more than 60, from which nine matched the rat mRNA [39]. The literature-reported union motifs for PTB: CTC [16], TCC [40], CCT [16], CTTA [41], and CCCT [16] are present in the deleted sequence as well as the CCTGC sequence of the polypyrimidine tract [42,43]. The presence of the 3’ splicing site near the PT tract, the presence of several Sironi’s motifs, the literature-reported repressor motifs for PTB, and our RIP data (Figure 1) strongly suggest that PTB can bind the deleted sequence in the mRNA and inhibit the intronic definition [44]. More research is needed to elucidate the specific sites involved in the intronic definition via PTB.

4.2. D1R Modulation of PTB Phosphorylation, Function, and Cellular Location

Our data indicate that D1R activation by SKF 38393 promotes the robust phosphorylation of PTB (almost 150%, Figure 2), which is D1R-specific and PKA-dependent. Knowing that PKA is the main regulator of PTB phosphorylation, and according to its nuclear location [18,45], this result is expected. Moreover, to our knowledge, this is the first report in a naïve system in which D1R phosphorylates PTB and modifies its splicing repressor activity during its blockade, as well as changing its cellular location and affecting the mRNA and protein expression of the D3nf isoform.
The effects of D1R on PTB functionality can be observed better in the long-time receptor blockade with SCH23390. That is because it produces chronic effects in the expression, location, and function of proteins processed by the regulated mRNA splicing, which is the case in dopaminergic denervation [10,15] and the 5-day blockade of D1R (this report). In the striatum, D1R activates adenylyl cyclase and incrementally effects cAMP production, which activates PKA [46]. When PKA is activated, the catalytic subunit translocates to the nucleus and phosphorylates PTB at Ser 16 [45]; then, PTB relocates from the nucleus to the cytoplasm, decreasing the repression of splicing and increasing the expression of isoforms: [47] in this case, D3nf. On the contrary, based on this proposal regarding the regulation of PTB by D1R, when D1R is chronically blockaded, it maintains decreased PKA activity, whereby no displacement of PTB occurs, remaining in the nucleus (Figure 6), repressing D3R splicing, and, as a consequence, decreasing D3nf isoform mRNA and protein expression (Figure 3).
Our model of a 5-day D1R blockade has two main characteristics that are required for producing the changes observed: D1R blockade and decreased PKA activity. It has been shown that injected SCH 23390 (at a dose of 0.1 mg/kg i.p.) has a peak concentration in the striatum for upwards of 3 h, whereas in plasma, it decreases over 1 h. This long-lasting concentration in the brain correlates with a decrement in Vmax for adenylyl cyclase and the locomotor activity (up to 80%) 12 h after injection concerning the control [48]. Additionally, the changes in locomotor activity stimulated by D1R in the striatum correlate with changes in PKA activity [49]. These antecedents indicate that in our model, the two doses a day of 0.5 mgr./kg of SCH 23390 for 5 days would produce a persistent blockade of D1R activity since this dose of antagonist occupies 95% of receptors [50]. This long-time receptor occupancy decreases cAMP formation via adenylyl cyclase and decreases PKA activity. Although D1R blockade could increase D1R membranal receptors [51], they remain occupied by SCH 23390. The need for the 5-day treatment is because, for the blockade of D1R using a single dose, the decrease in D3nf expression at 6 h tends to recover at 12 h, as can be seen in Figure 3A,B; thus, no significant accumulated changes in the expression of proteins was evident. On the other hand, studying the regulation of the splicing of D3R mRNA is important because the changes in D3R and D3nf isoform expression have important repercussions in chronic diseases in which D1R is not active or overactivated [10,14,35].
Interestingly, D1R blockade did not modify D3R mRNA or protein (Figure 3 and Figure 5), except for its cellular location. One can expect that an increment in canonical D3R could compensate for a decrease in the expression of D3nf, but that is not the case. It should be considered that the relative expression of D3nf mRNA concerning D3R mRNA in a normal condition is 0.6%, and this decreases up to 0.1% during blockade (Figure 3). The amount of D3R mRNA is too high compared to D3nf mRNA; thus, no significant change can be observed in our experimental conditions if such compensation occurs.
A fact that has claimed attention is the amount of D3nf protein expressed, which is almost comparable to that of D3R protein; this is essential since it has a physiological role: regulating D3R expression into the membrane [14,35]. How can this relatively low mRNA produce a significant amount of protein? This requires more experimental research; however, PTB has also been noted to be involved in maintaining mature mRNA stability and ribosomal transduction into protein [47,52]. Thus, PTB may be the factor that ensures enough levels of D3nf protein. On the other hand, in neurons, PTB is a factor capable of regulating the level of genic expression by inhibiting the exportation of incompletely spliced mRNA from the cytoplasm, a factor that triggers the nuclear degradation of the messenger (in this case, of D3R mRNA [53]), thus equalizing the translation into the protein of the D3R and D3nf isoforms [10,15]. More research is needed to clarify these points.
Finally, the changes in the location of the D3R protein are related to the decrement in the D3nf isoform since D3nf dimerizes the D3R protein, retaining it in the cytoplasm ([10,14]); the decrement in the expression of D3nf promotes D3R membranal location (Figure 5) [10]. In summary, our data indicate that through PTB phosphorylation, D1R regulates the expression of mRNA and the protein location of the PTB, D3R, and D3nf isoforms.

4.3. D1R Regulation of D3R Splicing Modifies Its Function

The central hypothesis for the role of D3nf in D1R and D3R interaction comes from Ritchand’s research, which proposes that the D3nf isoform regulates D3R membranal expression [14]. In an antagonistic interaction, the D3R opposes D1R effects, and by sequestering the canonical D3R into the cytoplasm, D3nf modulates dopamine’s effect on the receptors that favor the action in the D1R. According to a condition in which the D3nf isoform decreases, more D3R moves to the membrane and efficiently antagonizes D1R activity. Thus, Ritchand’s proposal explains the increased membrane expression of D3R observed in dopaminergic denervation and in the 5-day blockade of D1R [10,15,54], except when there is no favor for D1R activity because the normal response between D1R and D3R are synergistic. The most important finding in this work is that D1R regulates the expression of D3R in the membrane by regulating splicing. On the other hand, in the nucleus accumbens, D3R receptor activation also regulates D3R and D3nf mRNAs expression [55]; thus, in dopaminergic denervation, it is expected that the lack of stimulation of D3R will reduce both D3R and D3nf messengers. However, in our previous report in the denervated striatum, this is not the case; only D3nf messenger is decreased, and the fact that D1R blockade mimics the effects of denervation supports the notion of the role of D1R in regulating the membranal expression of D3R and decrements in the D3nf isoform [10,54]. It is interesting to note that the function of the nucleus accumbens regarding D3R is antagonistic with respect to D1R; thus, increments in D3R expression will produce less D1R activity in terms of competition for the dopamine available [14] in the case of the dorsal striatum when D3R is retained in the cytoplasm and favors the activity of the D1R-D3R dimmers; this regulation marks substantial differences in the function of D3R in both regions of the striatum.
In the dorsal striatum, D1R and D3R interact in proportion to medium-sized spine striatal neurons and their projections [6,7,8,9]. Therein, D3R potentiates D1R effects (atypical signaling) in a dimeric interaction; thus, during denervation or 5 days of SCH 23390 treatment, the new membranal functional D3R appears, which opposes (typical signaling) D1R-D3R synergistic effects [10]. This new membrane receptor is functionally coupled to the Gi protein, masking atypical signaling since Gi protein blockade unmasks a synergistic relationship that is not lost (Figure 7). In agreement with this, no significant changes in D3R mRNA occur (Figure 2) [10], nor do they happen in the number of D1R-D3R mRNA-expressing neurons during denervation [56,57]. These observations imply that in a normal condition, not all D3R protein is located in the membrane; that is, an amount of protein is retained in the cytoplasm by the D3nf isoform and the action regulated by D1R (Figure 5). Why does a neuron produce D3R protein and keep it in the cytoplasm but not in the membrane? This requires research; however, since D1R-D3R participates in regulating motor behavior [58], the regulation of membranal D3R expression is critical since its activation will antagonize the motor effects of the dimmer [10]. In this sense, it can be proposed that regulating D3R splicing is a normal mechanism to control the adequate function of the D1R-D3R dimmer. In summary, the data support the notion of D1R regulating D3R splicing and function, which is altered during dopaminergic denervation by a lack of activity of dopamine on D1R, leading to a modified functional response (Figure 8).

5. Conclusions

Our data indicate that through PKA→PTB, D1R modulates D3R splicing, expression, and signaling, which is altered during D1R blockade or a lack of stimulation, as is the case in dopaminergic denervation. The dysregulation of D3R expression implies the co-existence of synergistic and antagonistic effects on D1R function, altering the behavioral response to their activation.
Since several patho-physiological processes are related to the changes in the splicing of D3R that modify its signal and location [35], knowledge of those factors that determine the regulation of the splicing is essential to understand the changes in the pathology of Parkinson’s disease and to propose of new therapeutical approach.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines12010206/s1.

Author Contributions

Conceptualization, O.C.-D., J.S. and B.F.; methodology, O.C.-D., J.A.A.-F., M.L.-L., G.T.-M., C.D.F.-H., K.G.M.-N., H.C., R.F. and B.F.; formal analysis, O.C.-D., J.A.A.-F., M.L.-L., G.T.-M., C.D.F.-H., K.G.M.-N. and B.F.; investigation, O.C-D., M.L.-L., G.T.-M., C.D.F.-H., K.G.M.-N. and B.F.; resources, R.F., J.S. and B.F.; data curation, O.C-D., J.A.A.-F., M.L.-L., G.T.-M., C.D.F.-H., K.G.M.-N. and B.F.; writing—original draft preparation, O.C.-D., C.D.F.-H. and B.F.; writing—review and editing, O.C-D. and B.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Center for Research and Advanced Studies of the National Polytechnic Institute to BF.

Institutional Review Board Statement

The animal study was approved by the Institutional Animal Care Committee of the CINVESTAV. Approval Code: 0146-15 Approval Date: 1 September 2021, until 31 December 2024.

Informed Consent Statement

This study was not conduced in humans.

Data Availability Statement

Data are available upon request due to privacy/ethical considerations.

Acknowledgments

The authors thank Francisco Paz for technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ridray, S.; Griffon, N.; Mignon, V.; Souil, E.; Carboni, S.; Diaz, J.; Schwartz, J.C.; Sokoloff, P. Coexpression of dopamine D1 and D3 receptors in islands of Calleja and shell of nucleus accumbens of the rat: Opposite and synergistic functional interactions. Eur. J. Neurosci. 1998, 10, 1676–1686. [Google Scholar] [CrossRef]
  2. Schwartz, J.C.; Ridray, S.; Bordet, R.; Diaz, J.; Sokoloff, P. D1/D3 receptor relationships in brain coexpression, coactivation, and coregulation. Adv. Pharmacol. 1998, 42, 408–411. [Google Scholar] [PubMed]
  3. Schwartz, J.C.; Diaz, J.; Bordet, R.; Griffon, N.; Perachon, S.; Pilon, C.; Ridray, S.; Sokoloff, P. Functional implications of multiple dopamine receptor subtypes: The D1/D3 receptor coexistence. Brain Res. Brain Res. Rev. 1998, 26, 236–242. [Google Scholar] [CrossRef]
  4. Liu, X.Y.; Mao, L.M.; Zhang, G.C.; Papasian, C.J.; Fibuch, E.E.; Lan, H.X.; Zhou, H.F.; Xu, M.; Wang, J.Q. Activity-dependent modulation of limbic dopamine D3 receptors by CaMKII. Neuron 2009, 61, 425–438. [Google Scholar] [CrossRef] [PubMed]
  5. Barik, S.; De Beaurepaire, R. Hypothermic effects of dopamine D3 receptor agonists in the island of Calleja Magna. Potentiation by D1 activation. Pharmacol. Biochem. Behav. 1998, 60, 313–319. [Google Scholar] [CrossRef] [PubMed]
  6. Fiorentini, C.; Busi, C.; Gorruso, E.; Gotti, C.; Spano, P.; Missale, C. Reciprocal regulation of dopamine D1 and D3 receptor function and trafficking by heterodimerization. Mol. Pharmacol. 2008, 74, 59–69. [Google Scholar] [CrossRef]
  7. Marcellino, D.; Ferre, S.; Casado, V.; Cortes, A.; Le Foll, B.; Mazzola, C.; Drago, F.; Saur, O.; Stark, H.; Soriano, A.; et al. Identification of dopamine D1-D3 receptor heteromers. Indications for a role of synergistic D1-D3 receptor interactions in the striatum. J. Biol. Chem. 2008, 283, 26016–26025. [Google Scholar] [CrossRef]
  8. Avalos-Fuentes, A.; Loya-López, S.; Flores-Pérez, A.; Recillas-Morales, S.; Cortés, H.; Paz-Bermúdez, F.; Aceves, J.; Erlij, D.; Florán, B. Presynaptic CaMKIIα modulates dopamine D3 receptor activation in striatonigral terminals of the rat brain in a Ca2+ dependent manner. Neuropharmacology 2013, 71, 273–281. [Google Scholar] [CrossRef]
  9. Cruz-Trujillo, R.; Avalos-Fuentes, A.; Rangel-Barajas, C.; Paz-Bermudez, F.; Sierra, A.; Escartin-Perez, E.; Aceves, J.; Erlij, D.; Floran, B. D3 dopamine receptors interact with dopamine D1 but not D4 receptors in the GABAergic terminals of the SNr of the rat. Neuropharmacology 2013, 67, 370–378. [Google Scholar] [CrossRef]
  10. Campos Campos, B.; Ávalos-Fuentes, A.; Pina Leyva, C.; Sánchez-Zavaleta, R.; Loya-López, S.; Rangel-Barajas, C.; Leyva-Gómez, G.; Cortés, H.; Erlij, D.; Florán, B. Coexistence of D3R typical and atypical signaling in striatonigral neurons during dopaminergic denervation. Correlation with D3nf expression changes. Synapse 2020, 74, e22152. [Google Scholar] [CrossRef]
  11. Liu, K.; Bergson, C.; Levenson, R.; Schmauss, C. On the origin of mRNA encoding the truncated dopamine D3-type receptor D3nf and detection of D3nf-like immunoreactivity in human brain. J. Biol. Chem. 1994, 269, 29220–29226. [Google Scholar] [CrossRef] [PubMed]
  12. Elmhurst, J.L.; Xie, Z.; O’Dowd, B.F.; George, S.R. The splice variant D3nf reduces ligand binding to the D3 dopamine receptor: Evidence for heterooligomerization. Brain Res. Mol. Brain Res. 2000, 80, 63–74. [Google Scholar] [CrossRef] [PubMed]
  13. Karpa, K.D.; Lin, R.; Kabbani, N.; Levenson, R. The dopamine D3 receptor interacts with itself and the truncated D3 splice variant d3nf: D3-D3nf interaction causes mislocalization of D3 receptors. Mol. Pharmacol. 2000, 58, 677–683. [Google Scholar] [CrossRef]
  14. Richtand, N.M. Behavioral sensitization, alternative splicing, and d3 dopamine receptor-mediated inhibitory function. Neuropsychopharmacology 2006, 31, 2368–2375. [Google Scholar] [CrossRef] [PubMed]
  15. Prieto, G.A.; Perez-Burgos, A.; Palomero-Rivero, M.; Galarraga, E.; Drucker-Colin, R.; Bargas, J. Upregulation of D2-class signaling in dopamine-denervated striatum is in part mediated by D3 receptors acting on Ca V 2.1 channels via PIP2 depletion. J. Neurophysiol. 2011, 105, 2260–2274. [Google Scholar] [CrossRef]
  16. Oberstrass, F.C.; Auweter, S.D.; Erat, M.; Hargous, Y.; Henning, A.; Wenter, P.; Reymond, L.; Amir-Ahmady, B.; Pitsch, S.; Black, D.L. Structure of PTB bound to RNA: Specific binding and implications for splicing regulation. Science 2005, 309, 2054–2057. [Google Scholar] [CrossRef] [PubMed]
  17. Sasabe, T.; Futai, E.; Ishiura, S. Polypyrimidine tract-binding protein 1 regulates the alternative splicing of dopamine receptor D2. J. Neurochem. 2011, 116, 76–81. [Google Scholar] [CrossRef]
  18. Ma, S.; Liu, G.; Sun, Y.; Xie, J. Relocalization of the polypyrimidine tract-binding protein during PKA-induced neurite growth. Biochim. Biophys. Acta (BBA)-Mol. Cell Res. 2007, 1773, 912–923. [Google Scholar] [CrossRef]
  19. Pritchard, L.M.; Logue, A.D.; Taylor, B.C.; Ahlbrand, R.; Welge, J.A.; Tang, Y.; Sharp, F.R.; Richtand, N.M. Relative expression of D3 dopamine receptor and alternative splice variant D3nf mRNA in high and low responders to novelty. Brain Res. Bull. 2006, 70, 296–303. [Google Scholar] [CrossRef]
  20. Paxinos, G.; Watson, C. The Rat Brain in Stereotaxic Coordinates: Hard Cover Edition; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  21. Heydorn, W.E.; Creed, G.J.; Patel, J.; Jacobowitz, D.M. Distribution of proteins in different subcellular fractions of rat brain studied by two-dimensional gel electrophoresis. Neurochem. Int. 1986, 9, 357–370. [Google Scholar] [CrossRef]
  22. Guillemin, I.; Becker, M.; Ociepka, K.; Friauf, E.; Nothwang, H.G. A subcellular prefractionation protocol for minute amounts of mammalian cell cultures and tissue. Proteomics 2005, 5, 35–45. [Google Scholar] [CrossRef] [PubMed]
  23. Gagliardi, M.; Matarazzo, M.R. Rip: Rna Immunoprecipitation. In Polycomb Group Proteins: Methods and Protocols; Springer: Berlin/Heidelberg, Germany, 2016; pp. 73–86. [Google Scholar]
  24. Heim, M.; Blot, L.; Besse, F. An RNA-immunoprecipitation protocol to identify RNAs associated with RNA-binding proteins in cytoplasmic and nuclear Drosophila head fractions. STAR Protoc. 2022, 3, 101415. [Google Scholar] [CrossRef]
  25. Loya-López, S.; Sandoval, A.; González-Ramírez, R.; Calderón-Rivera, A.; Ávalos-Fuentes, A.; Rodríguez-Sánchez, M.; Caballero, R.; Tovar-Soto, D.; Felix, R.; Florán, B. Cdk5 phosphorylates CaV1. 3 channels and regulates GABAA-mediated miniature inhibitory post-synaptic currents in striato-nigral terminals. Biochem. Biophys. Res. Commun. 2020, 524, 255–261. [Google Scholar] [CrossRef] [PubMed]
  26. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  27. Rangel-Barajas, C.; Silva, I.; Lopéz-Santiago, L.M.; Aceves, J.; Erlij, D.; Florán, B. L-DOPA-induced dyskinesia in hemiparkinsonian rats is associated with up-regulation of adenylyl cyclase type V/VI and increased GABA release in the substantia nigra reticulata. Neurobiol. Dis. 2011, 41, 51–61. [Google Scholar] [CrossRef]
  28. Dubois, A.; Savasta, M.; Curet, O.; Scatton, B. Autoradiographic distribution of the D1 agonist [3H] SKF 38393, in the rat brain and spinal cord. Comparison with the distribution of D2 dopamine receptors. Neuroscience 1986, 19, 125–137. [Google Scholar] [CrossRef]
  29. Bourne, J.A. SCH 23390: The first selective dopamine D1-like receptor antagonist. CNS Drug Rev. 2001, 7, 399–414. [Google Scholar] [CrossRef] [PubMed]
  30. Chijiwa, T.; Mishima, A.; Hagiwara, M.; Sano, M.; Hayashi, K.; Inoue, T.; Naito, K.; Toshioka, T.; Hidaka, H. Inhibition of forskolin-induced neurite outgrowth and protein phosphorylation by a newly synthesized selective inhibitor of cyclic AMP-dependent protein kinase, N-[2-(p-bromocinnamylamino) ethyl]-5-isoquinolinesulfonamide (H-89), of PC12D pheochromocytoma cells. J. Biol. Chem. 1990, 265, 5267–5272. [Google Scholar]
  31. Nimchinsky, E.A.; Hof, P.R.; Janssen, W.G.; Morrison, J.H.; Schmauss, C. Expression of dopamine D3 receptor dimers and tetramers in brain and in transfected cells. J. Biol. Chem. 1997, 272, 29229–29237. [Google Scholar] [CrossRef]
  32. Ikeda, E.; Matsunaga, N.; Kakimoto, K.; Hamamura, K.; Hayashi, A.; Koyanagi, S.; Ohdo, S. Molecular mechanism regulating 24-hour rhythm of dopamine D3 receptor expression in mouse ventral striatum. Mol. Pharmacol. 2013, 83, 959–967. [Google Scholar] [CrossRef]
  33. Pugsley, T.A.; Davis, M.D.; Akunne, H.C.; MacKenzie, R.G.; Shih, Y.H.; Damsma, G.; Wikstrom, H.; Whetzel, S.Z.; Georgic, L.M.; Cooke, L.W.; et al. Neurochemical and functional characterization of the preferentially selective dopamine D3 agonist PD 128907. J. Pharmacol. Exp. Ther. 1995, 275, 1355–1366. [Google Scholar] [PubMed]
  34. Shapiro, M.S.; Wollmuth, L.P.; Hille, B. Modulation of Ca2+ channels by PTX-sensitive G-proteins is blocked by N-ethylmaleimide in rat sympathetic neurons. J. Neurosci. 1994, 14, 7109–7116. [Google Scholar] [CrossRef] [PubMed]
  35. Prieto, G.A. Abnormalities of Dopamine D3 Receptor Signaling in the Diseased Brain. J. Cent. Nerv. Syst. Dis. 2017, 9, 1179573517726335. [Google Scholar] [CrossRef] [PubMed]
  36. Schmauss, C.; Haroutunian, V.; Davis, K.L.; Davidson, M. Selective loss of dopamine D3-type receptor mRNA expression in parietal and motor cortices of patients with chronic schizophrenia. Proc. Natl. Acad. Sci. USA 1993, 90, 8942–8946. [Google Scholar] [CrossRef]
  37. Schmauss, C. Enhanced cleavage of an atypical intron of dopamine D3-receptor pre-mRNA in chronic schizophrenia. J. Neurosci. 1996, 16, 7902–7909. [Google Scholar] [CrossRef] [PubMed]
  38. Burset, M.; Seledtsov, I.A.; Solovyev, V.V. Analysis of canonical and non-canonical splice sites in mammalian genomes. Nucleic Acids Res. 2000, 28, 4364–4375. [Google Scholar] [CrossRef]
  39. Sironi, M.; Menozzi, G.; Riva, L.; Cagliani, R.; Comi, G.P.; Bresolin, N.; Giorda, R.; Pozzoli, U. Silencer elements as possible inhibitors of pseudoexon splicing. Nucleic Acids Res. 2004, 32, 1783–1791. [Google Scholar] [CrossRef]
  40. van Bergeijk, P.; Seneviratne, U.; Aparicio-Prat, E.; Stanton, R.; Hasson, S.A. SRSF1 and PTBP1 are trans-acting factors that suppress the formation of a CD33 splicing isoform linked to Alzheimer’s disease risk. Mol. Cell. Biol. 2019, 39, e00568–18. [Google Scholar] [CrossRef]
  41. Xue, Y.; Zhou, Y.; Wu, T.; Zhu, T.; Ji, X.; Kwon, Y.-S.; Zhang, C.; Yeo, G.; Black, D.L.; Sun, H. Genome-wide analysis of PTB-RNA interactions reveals a strategy used by the general splicing repressor to modulate exon inclusion or skipping. Mol. Cell 2009, 36, 996–1006. [Google Scholar] [CrossRef]
  42. Singh, R.; Valcarcel, J.; Green, M.R. Distinct binding specificities and functions of higher eukaryotic polypyrimidine tract-binding proteins. Science 1995, 268, 1173–1176. [Google Scholar] [CrossRef]
  43. Corioni, M.; Antih, N.; Tanackovic, G.; Zavolan, M.; Krämer, A. Analysis of in situ pre-mRNA targets of human splicing factor SF1 reveals a function in alternative splicing. Nucleic Acids Res. 2011, 39, 1868–1879. [Google Scholar] [CrossRef]
  44. Spellman, R.; Smith, C.W. Novel modes of splicing repression by PTB. Trends Biochem. Sci. 2006, 31, 73–76. [Google Scholar] [CrossRef]
  45. Xie, J.; Lee, J.-A.; Kress, T.L.; Mowry, K.L.; Black, D.L. Protein kinase A phosphorylation modulates transport of the polypyrimidine tract-binding protein. Proc. Natl. Acad. Sci. USA 2003, 100, 8776–8781. [Google Scholar] [CrossRef] [PubMed]
  46. Undieh, A.S. Pharmacology of signaling induced by dopamine D1-like receptor activation. Pharmacol. Ther. 2010, 128, 37–60. [Google Scholar] [CrossRef] [PubMed]
  47. Sawicka, K.; Bushell, M.; Spriggs, K.A.; Willis, A.E. Polypyrimidine-tract-binding protein: A multifunctional RNA-binding protein. Biochem. Soc. Trans. 2008, 36, 641–647. [Google Scholar] [CrossRef]
  48. Schulz, D.W.; Staples, L.; Mailman, R.B. SCH23390 causes persistent antidopaminergic effects in vivo: Evidence for longterm occupation of receptors. Life Sci. 1985, 36, 1941–1948. [Google Scholar] [CrossRef]
  49. Ma, L.; Day-Cooney, J.; Benavides, O.J.; Muniak, M.A.; Qin, M.; Ding, J.B.; Mao, T.; Zhong, H. Locomotion activates PKA through dopamine and adenosine in striatal neurons. Nature 2022, 611, 762–768. [Google Scholar] [CrossRef] [PubMed]
  50. Haraguchi, K.; Ito, K.; Kotaki, H.; Sawada, Y.; Iga, T. Prediction of drug-induced catalepsy based on dopamine D1, D2, and muscarinic acetylcholine receptor occupancies. Drug Metab. Dispos. 1997, 25, 675–684. [Google Scholar]
  51. Porceddu, M.L.; Ongini, E.; Biggio, G. [3H] SCH 23390 binding sites increase after chronic blockade of D-1 dopamine receptors. Eur. J. Pharmacol. 1985, 118, 367–370. [Google Scholar] [CrossRef]
  52. Mitchell, S.A.; Spriggs, K.A.; Bushell, M.; Evans, J.R.; Stoneley, M.; Le Quesne, J.P.; Spriggs, R.V.; Willis, A.E. Identification of a motif that mediates polypyrimidine tract-binding protein-dependent internal ribosome entry. Genes. Dev. 2005, 19, 1556–1571. [Google Scholar] [CrossRef]
  53. Yap, K.; Lim, Z.Q.; Khandelia, P.; Friedman, B.; Makeyev, E.V. Coordinated regulation of neuronal mRNA steady-state levels through developmentally controlled intron retention. Genes. Dev. 2012, 26, 1209–1223. [Google Scholar] [CrossRef] [PubMed]
  54. Avalos-Fuentes, A.; Albarran-Bravo, S.; Loya-Lopez, S.; Cortes, H.; Recillas-Morales, S.; Magana, J.J.; Paz-Bermudez, F.; Rangel-Barajas, C.; Aceves, J.; Erlij, D.; et al. Dopaminergic denervation switches dopamine D3 receptor signaling and disrupts its Ca(2+) dependent modulation by CaMKII and calmodulin in striatonigral projections of the rat. Neurobiol. Dis. 2015, 74, 336–346. [Google Scholar] [CrossRef] [PubMed]
  55. Richtand, N.M.; Liu, Y.; Ahlbrand, R.; Sullivan, J.R.; Newman, A.H.; McNamara, R.K. Dopaminergic regulation of dopamine D3 and D3nf receptor mRNA expression. Synapse 2010, 64, 634–643. [Google Scholar] [CrossRef] [PubMed]
  56. Lanza, K.; Centner, A.; Coyle, M.; Del Priore, I.; Manfredsson, F.P.; Bishop, C. Genetic suppression of the dopamine D3 receptor in striatal D1 cells reduces the development of L-DOPA-induced dyskinesia. Exp. Neurol. 2021, 336, 113534. [Google Scholar] [CrossRef]
  57. Solís, O.; Garcia-Montes, J.R.; González-Granillo, A.; Xu, M.; Moratalla, R. Dopamine D3 receptor modulates l-DOPA-induced dyskinesia by targeting D1 receptor-mediated striatal signaling. Cereb. Cortex 2017, 27, 435–446. [Google Scholar] [CrossRef]
  58. Moreno, E.; Casajuana-Martin, N.; Coyle, M.; Campos, B.C.; Galaj, E.; del Torrent, C.L.; Seyedian, A.; Rea, W.; Cai, N.-S.; Bonifazi, A. Pharmacological targeting of G protein-coupled receptor heteromers. Pharmacol. Res. 2022, 185, 106476. [Google Scholar] [CrossRef]
Figure 1. D3R mRNA co-precipitates with PTB. (A) A qRT-PCR amplification curve for the D3R gene from the extract of PTB immunoprecipitated (RIP) and the cDNA from striatal tissue (Str). (B) Agarose gel electrophoresis from the amplified cDNA recovered from the qRT-PCR.
Figure 1. D3R mRNA co-precipitates with PTB. (A) A qRT-PCR amplification curve for the D3R gene from the extract of PTB immunoprecipitated (RIP) and the cDNA from striatal tissue (Str). (B) Agarose gel electrophoresis from the amplified cDNA recovered from the qRT-PCR.
Biomedicines 12 00206 g001
Figure 2. Activation of D1R by SKF 38393 phosphorylates PTB by PKA. (A,C) Representative blots of PTB phosphorylation in homogenates from striatal slices treated with the D1R agonist SKF 38393 and the D1R antagonist SCH 23390 or the PKA blocker H89. (B,D) The densitometry analysis from three blots. ** p < 0.01; *** p < 0.001; ns: not significant differences. (B) F = 30.96, p = 0.0007; (D) F = 19.86, p = 0.0023, One-way ANOVA, followed by Tukey’s test, n = 3.
Figure 2. Activation of D1R by SKF 38393 phosphorylates PTB by PKA. (A,C) Representative blots of PTB phosphorylation in homogenates from striatal slices treated with the D1R agonist SKF 38393 and the D1R antagonist SCH 23390 or the PKA blocker H89. (B,D) The densitometry analysis from three blots. ** p < 0.01; *** p < 0.001; ns: not significant differences. (B) F = 30.96, p = 0.0007; (D) F = 19.86, p = 0.0023, One-way ANOVA, followed by Tukey’s test, n = 3.
Biomedicines 12 00206 g002
Figure 3. D1R receptor blockade with SCH 23390 decreases D3nf mRNA expression and relative expression concerning D3R mRNA. (A) The effect of a single dose of the D1R antagonist SCH 23390 with respect to saline in striatal D3nf and D3R mRNAs at 6 and 12 h after injection. (B) The D3nf/D3R+D3nf ratio of the relative expression of the isoforms. (C,D) The same determinations after a 5-day two-dose antagonist treatment. * p < 0.05; ** p < 0.01; *** p < 0.001; ns: not significant; ## p < 0.01 with respect to the indicated group. (A) F = 4.29, p = 0.0056; (B) F = 10.25, p = 0.0021; (C) F = 18.72, p < 0.0001. One-way ANOVA, followed by Tukey’s test, n = 3–6 determinations, (D) t = 4.91, df = 6, unpaired t-test n = 4.
Figure 3. D1R receptor blockade with SCH 23390 decreases D3nf mRNA expression and relative expression concerning D3R mRNA. (A) The effect of a single dose of the D1R antagonist SCH 23390 with respect to saline in striatal D3nf and D3R mRNAs at 6 and 12 h after injection. (B) The D3nf/D3R+D3nf ratio of the relative expression of the isoforms. (C,D) The same determinations after a 5-day two-dose antagonist treatment. * p < 0.05; ** p < 0.01; *** p < 0.001; ns: not significant; ## p < 0.01 with respect to the indicated group. (A) F = 4.29, p = 0.0056; (B) F = 10.25, p = 0.0021; (C) F = 18.72, p < 0.0001. One-way ANOVA, followed by Tukey’s test, n = 3–6 determinations, (D) t = 4.91, df = 6, unpaired t-test n = 4.
Biomedicines 12 00206 g003
Figure 4. D1R blockade with SCH 23390 decreases striatal D3nf protein expression in the cytoplasm. (A,C,E) WB from D3nf protein from total striatal homogenates and protein fraction from striatal membranes or cytoplasm, respectively. (B,D,F) The densitometry analysis. * p < 0.05; ** p < 0.01; ns: not significant: (B), t = 4.81, df = 6, n = 4; (D), t = 0.65, df = 4, n = 3; (F), t = 2.96, df = 4, n = 3. Unpaired t-test.
Figure 4. D1R blockade with SCH 23390 decreases striatal D3nf protein expression in the cytoplasm. (A,C,E) WB from D3nf protein from total striatal homogenates and protein fraction from striatal membranes or cytoplasm, respectively. (B,D,F) The densitometry analysis. * p < 0.05; ** p < 0.01; ns: not significant: (B), t = 4.81, df = 6, n = 4; (D), t = 0.65, df = 4, n = 3; (F), t = 2.96, df = 4, n = 3. Unpaired t-test.
Biomedicines 12 00206 g004
Figure 5. D1R blockade with SCH 23390 did not change striatal D3R protein expression but increased it in the membrane. (A,C,E) WB from D3R protein from total striatal homogenates and protein fraction from striatal membranes or cytoplasm, respectively; the 45 kD bands represent the D3R monomeric forms and the 75 kD bands represent the dimeric forms. (B,D,F) The densitometry analysis. * p < 0.05; ** p < 0.01; **** p < 0.0001; ns: not significant. (B), F = 0.001, p > 0.999; (D), F = 15.22, p = 0.0011; (C), F = 31.73, p < 0.0001. One-way ANOVA, followed by Tukey n = 3.
Figure 5. D1R blockade with SCH 23390 did not change striatal D3R protein expression but increased it in the membrane. (A,C,E) WB from D3R protein from total striatal homogenates and protein fraction from striatal membranes or cytoplasm, respectively; the 45 kD bands represent the D3R monomeric forms and the 75 kD bands represent the dimeric forms. (B,D,F) The densitometry analysis. * p < 0.05; ** p < 0.01; **** p < 0.0001; ns: not significant. (B), F = 0.001, p > 0.999; (D), F = 15.22, p = 0.0011; (C), F = 31.73, p < 0.0001. One-way ANOVA, followed by Tukey n = 3.
Biomedicines 12 00206 g005
Figure 6. D1R blockade with SCH 23390 did not change striatal PTB protein expression but decreased it in the cytoplasm and increased it in the nucleus. (A,C,E) WB for D3R protein from total striatal homogenates and protein fraction from striatal cytoplasm or nucleus, respectively. (B,D,F) The densitometry analysis. * p < 0.05; ** p < 0.01; ns: not significant differences, (B), t = 0.03, df = 4, n = 3; (D), t = 3.5, df = 4, n = 3; (F), t = 5.7, df = 4, n = 3. Unpaired t-test.
Figure 6. D1R blockade with SCH 23390 did not change striatal PTB protein expression but decreased it in the cytoplasm and increased it in the nucleus. (A,C,E) WB for D3R protein from total striatal homogenates and protein fraction from striatal cytoplasm or nucleus, respectively. (B,D,F) The densitometry analysis. * p < 0.05; ** p < 0.01; ns: not significant differences, (B), t = 0.03, df = 4, n = 3; (D), t = 3.5, df = 4, n = 3; (F), t = 5.7, df = 4, n = 3. Unpaired t-test.
Biomedicines 12 00206 g006
Figure 7. D1R blockade with SCH 23390 produces the typical D3R response that masks the atypical response. (A,B) cAMP accumulation experiments in control striatal slices and 5-day SCH23390 treatment. (A) D1R activation by SKF38393 increases cAMP, and the D3R activation by PD 128,907 potentiates it; (B) The response of D1R activation is prevented by D3R activation. (C,D) are shown in the same experiments but in the presence of Gi protein coupling blocked by NEM, which restores the D3R atypical response in (D). ** p < 0.01, *** p < 0.001; ns: not significant. (A) F = 37.79, **** p < 0.0001; (B) F = 116.7, p < 0.0001; (C) F = 53.55, p < 0.0001; (D) F = 76.88, p < 0.0001. One-way ANOVA, followed by Tukey n = 3.
Figure 7. D1R blockade with SCH 23390 produces the typical D3R response that masks the atypical response. (A,B) cAMP accumulation experiments in control striatal slices and 5-day SCH23390 treatment. (A) D1R activation by SKF38393 increases cAMP, and the D3R activation by PD 128,907 potentiates it; (B) The response of D1R activation is prevented by D3R activation. (C,D) are shown in the same experiments but in the presence of Gi protein coupling blocked by NEM, which restores the D3R atypical response in (D). ** p < 0.01, *** p < 0.001; ns: not significant. (A) F = 37.79, **** p < 0.0001; (B) F = 116.7, p < 0.0001; (C) F = 53.55, p < 0.0001; (D) F = 76.88, p < 0.0001. One-way ANOVA, followed by Tukey n = 3.
Biomedicines 12 00206 g007
Figure 8. Schematic drawing illustrating the D1R or D1R/D3R signaling to regulate the PTB modulation of the splicing of D3R in a normal condition and D1R blockade or dopaminergic denervation.
Figure 8. Schematic drawing illustrating the D1R or D1R/D3R signaling to regulate the PTB modulation of the splicing of D3R in a normal condition and D1R blockade or dopaminergic denervation.
Biomedicines 12 00206 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Casados-Delgado, O.; Avalos-Fuentes, J.A.; Lara-Lozano, M.; Tovar-Medina, G.; Florán-Hernández, C.D.; Martínez-Nolasco, K.G.; Cortes, H.; Felix, R.; Segovia, J.; Florán, B. Modulation of D3R Splicing, Signaling, and Expression by D1R through PKA→PTB Phosphorylation. Biomedicines 2024, 12, 206. https://doi.org/10.3390/biomedicines12010206

AMA Style

Casados-Delgado O, Avalos-Fuentes JA, Lara-Lozano M, Tovar-Medina G, Florán-Hernández CD, Martínez-Nolasco KG, Cortes H, Felix R, Segovia J, Florán B. Modulation of D3R Splicing, Signaling, and Expression by D1R through PKA→PTB Phosphorylation. Biomedicines. 2024; 12(1):206. https://doi.org/10.3390/biomedicines12010206

Chicago/Turabian Style

Casados-Delgado, Orlando, José Arturo Avalos-Fuentes, Manuel Lara-Lozano, Gisela Tovar-Medina, Carla Daniela Florán-Hernández, Karla Gisela Martínez-Nolasco, Hernán Cortes, Ricardo Felix, José Segovia, and Benjamín Florán. 2024. "Modulation of D3R Splicing, Signaling, and Expression by D1R through PKA→PTB Phosphorylation" Biomedicines 12, no. 1: 206. https://doi.org/10.3390/biomedicines12010206

APA Style

Casados-Delgado, O., Avalos-Fuentes, J. A., Lara-Lozano, M., Tovar-Medina, G., Florán-Hernández, C. D., Martínez-Nolasco, K. G., Cortes, H., Felix, R., Segovia, J., & Florán, B. (2024). Modulation of D3R Splicing, Signaling, and Expression by D1R through PKA→PTB Phosphorylation. Biomedicines, 12(1), 206. https://doi.org/10.3390/biomedicines12010206

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