Reversal of autism-like behaviors and metabolism in adult mice with single-dose antipurinergic therapy.

Autism spectrum disorders (ASDs) now affect 1-2% of the children born in the United States. Hundreds of genetic, metabolic and environmental factors are known to increase the risk of ASD. Similar factors are known to influence the risk of schizophrenia and bipolar disorder; however, a unifying mechanistic explanation has remained elusive. Here we used the maternal immune activation (MIA) mouse model of neurodevelopmental and neuropsychiatric disorders to study the effects of a single dose of the antipurinergic drug suramin on the behavior and metabolism of adult animals. We found that disturbances in social behavior, novelty preference and metabolism are not permanent but are treatable with antipurinergic therapy (APT) in this model of ASD and schizophrenia. A single dose of suramin (20 mg kg(-1) intraperitoneally (i.p.)) given to 6-month-old adults restored normal social behavior, novelty preference and metabolism. Comprehensive metabolomic analysis identified purine metabolism as the key regulatory pathway. Correction of purine metabolism normalized 17 of 18 metabolic pathways that were disturbed in the MIA model. Two days after treatment, the suramin concentration in the plasma and brainstem was 7.64 μM pmol μl(-1) (±0.50) and 5.15 pmol mg(-1) (±0.49), respectively. These data show good uptake of suramin into the central nervous system at the level of the brainstem. Most of the improvements associated with APT were lost after 5 weeks of drug washout, consistent with the 1-week plasma half-life of suramin in mice. Our results show that purine metabolism is a master regulator of behavior and metabolism in the MIA model, and that single-dose APT with suramin acutely reverses these abnormalities, even in adults.

Introduction

Genetic,[1, 2, 3] environmental[4,5] and metabolic[6] factors can contribute to the risk of autism to different extents in each affected child. Despite this etiologic heterogeneity, and the well-known clinical variations that make each child unique, clinical studies suggest that a common denominator may underlie the shared behavioral and cognitive features that define autism spectrum disorders (ASDs) as a group. For example, in a prospective study conducted by the Zimmerman group at the Kennedy Krieger Institute in 2007, 83% of children with autism spectrum disorders were found to improve transiently in association with fever.[7] Improvements were lost with the resolution of the fever. This study showed that, despite the many different causes of ASD, the symptoms were not permanent and could be improved in a substantial fraction of children.

Transient improvements with fever have also been found in patients with certain forms of post-infection brain syndromes, movement disorders, dementia and schizophrenia in the early 1900s, although those early studies were made complicated by the use of live malarial parasites to produce the fevers.[8,9] In contrast to these beneficial effects, when the exposure to serious infection happens before the onset of disease—during early development, and particularly during pregnancy—the metabolic changes associated with significant fever or infection are known to increase the risk of neurodevelopmental disorders in the offspring. These disorders include schizophrenia,[10] ASDs,[11] attention deficit/hyperactivity disorder,[12] bipolar disorder,[13] epilepsy[14] and cerebral palsy.[15] The nature and developmental timing of the exposure are important. Metabolism and mitochondrial function change adaptively during and after infection, and are well-known regulators of neurotransmission[16] and synaptic plasticity.[17] Collectively, these studies suggest that, despite many different causes, the symptoms of several neurodevelopmental disorders such as ASD, schizophrenia and bipolar disorder may have a metabolic basis and be acutely responsive to treatment using the right metabolic intervention.

The maternal immune activation (MIA) mouse model of neurodevelopmental disorders produces symptoms that are biologically similar to those of ASD[18] and schizophrenia.[19] Pregnant females that are exposed to a simulated viral infection by injection of the double-stranded RNA poly(Inosine:Cytosine) produce offspring with features of ASD[20] and schizophrenia.[21] Exposure to poly(IC) activates an evolutionarily conserved metabolic response to a threat called the cell danger response (CDR).[22] Pathological persistence of the CDR, beyond the physical presence of the threat, has been observed in a variety of chronic disorders including ASDs.[22] Purinergic signaling has been hypothesized to be a key regulator of the CDR;[22] however, this has not yet been proven. In support of this hypothesis, we recently showed that antipurinergic therapy (APT) in the MIA mouse model corrected all of the behavioral, molecular and neuropathological abnormalities when weekly treatment with the antipurinergic drug suramin was begun at 1.5 months of age, near the age of reproductive maturity for mice.[23] Significant reductions in mitochondrial oxygen consumption and body temperature were also found. However, comprehensive metabolomic analysis was not reported in that study.[23]

In the present study, we tested the hypothesis that the behavioral manifestations of the MIA model are a consequence of pathological persistence of the evolutionarily conserved CDR,[22] and that the CDR is maintained by dysregulated purine metabolism and secondary abnormalities in purinergic signaling. We found that a single dose of the antipurinergic drug suramin given to adult animals about 6 months of age (21–27 weeks) produced the concerted correction of over 90% of the metabolic pathway disturbances, and all of the behavioral abnormalities that we tested in the MIA model. Six-month-old mice are the human biological age equivalents of about 30 years[24] (see Materials and methods). After washout of the drug, these improvements were lost and the former abnormalities returned. These data show that purine metabolism and purinergic signaling represent a novel neurochemical switch that regulates both behavior and metabolism in the MIA model of neurodevelopmental disorders such as ASD and schizophrenia.

Materials and methods

Animals and husbandry

All studies were conducted at the University of California, San Diego (UCSD) in facilities accredited by the Association for Assessment and Accreditation of Laboratory Animal Care International (AAALAC) under the UCSD Institutional Animal Care and Use Committee-approved animal subjects protocols, and followed the National Institutes of Health Guidelines for the use of animals in research. Six- to eight-week-old C57BL/6J (strain no. 000664) mice were obtained from Jackson Laboratories (Bar Harbor, ME, USA), given food and water ad libitum, identified by ear tags, and used to produce the timed matings. Animals were housed in a temperature- (22–24 °C) and humidity (40–55%)-controlled vivarium with a 12-h light–dark cycle (lights on at 0700 hours). Nulliparous dams were mated at 9–10 weeks of age. The sires were also 9–10 weeks of age. The human biological age equivalent for the C57BL/6J strain of laboratory mouse (Mus musculus) can be estimated from the following equation: 12 years for the first month, 6 years for the second month, 3 years for months 3–6 and 2.5 years for each month thereafter.[24] Therefore, a 6-month-old mouse would be the biological equivalent of 30 years old (=12+6+3 × 4) on a human timeline.

Poly(IC) preparation and gestational exposure

To initiate the MIA model, pregnant dams were given two intraperitoneal injections of Poly(I:C) (Potassium salt; Sigma-Aldrich, St. Louis, MO, USA, Cat no. P9582; >99% pure; <1% mononucleotide content). These were quantified by UV spectrophotometry. One unit (U) of poly(IC) was defined as 1 absorbance unit at 260 nm. Typically, 1U=12 μg of RNA. 0.25 U/g [3 mg kg−1] of poly(IC) was given on E12.5 and 0.125 U g−1 (1.5 mg kg−1) on E17.5 as previously described.[23] Contemporaneous control pregnancies were produced by timed matings and randomized assignment of pregnant dams to saline injection (5 μl g−1 intraperitoneally (i.p.)) on E12.5 and E17.5.

Postnatal handling and antipurinergic therapy (APT)

Offspring of timed matings were weaned at 3–4 weeks of age into cages of two to four animals. No mice were housed in isolation. Only males were evaluated in these studies. Littermates were identified by ear tags and distributed into different cages in order to minimize litter and dam effects. To avoid chance differences in groups selected for single-dose treatment, the saline and poly(IC) exposure groups were each balanced according to their social approach scores at 2.25 months. At 5.25 or 6.5 months of age, half the animals received a single injection of either saline (5 μl g−1 i.p.) or suramin (hexasodium salt, 20 mg kg−1 i.p.; Tocris Bioscience, Bristol, UK, Cat no. 1472). Beginning 2 days later, behaviors were evaluated as described below. After completing the behavioral measurements, half of the subjects were killed after a 5-week-washout period for measurement of suramin tissue levels. For acute suramin levels, the other half was injected at 7.75 months of age and killed 2 days later for tissue level determinations.

Behavioral testing

Behavioral testing began at 2.25 months (9 weeks) of age. Mice were tested in social approach, rotarod, t-maze test of spontaneous alternation and light–dark box test. If abnormalities were found, treatment with suramin or saline was given at 5.25 months (21 weeks) or 6.5–6.75 months (26–27 weeks) and the testing was repeated. Only male animals were tested.

Social approach

Social behavior was tested as social preference as previously described[23] with minor modifications (see Full Methods in Supplementary Information; N=19–25, 2.25-month-old males per group before adult treatment with suramin. N=8–13, 6.5-month-old males per group).

T-Maze

Novelty preference was tested as spontaneous alternation behavior in the T-maze by a modification of the methods of Frye and Walf[25] (see Full Methods Supplementary Information). N=19–25, 4-month-old males per group before adult treatment with suramin. N=8–13, 5.25-month-old males per group.

Rotarod

Sensorimotor coordination was tested as latency to fall on the rotarod as previously described[23] (see Full Methods Supplementary Information; N=19–25, 2.5-month-old males per group before adult treatment with suramin. N=8–13, 6.75-month-old males per group).

Light–dark box

Certain anxiety-related and light-avoidance behaviors were tested in the light–dark box paradigm as previously described[26] (see Full Methods Supplementary Information; N=19–25, 3.5-month-old males per group).

Absence of abnormal behaviors produced by suramin

This was assessed in the non-MIA control animals (indicated as the ‘Saline' group in the pretreatment figures) that were injected with suramin as adults (indicated as the ‘Sal-Sur' groups in the single-dose treatment figures) using each of the above behavioral paradigms.

Suramin quantitation

Tissue samples (brainstem, cerebrum and cerebellum) were ground into powder under liquid nitrogen in a pre-cooled mortar. Powdered tissue (15–50 mg) was weighed and mixed with the internal standard trypan blue to a final concentration of 5 μM (pmol mg−1) and incubated at room temperature for 10 min to permit metabolite interaction with binding proteins. Nine volumes of methanol:acetonitrile:H2O (43:43:16) pre-chilled to −20 °C was added to produce a final solvent ratio of 40:40:20, and the samples were deproteinated and macromolecules removed by precipitation on crushed ice for 30 min. The mixture was centrifuged at 16 000 g for 10 min at 4 °C and the supernatant was transferred to a new tube and kept at −80 °C for further LC-MS/MS (liquid chromatography-tandem mass spectrometry) analysis. For plasma, 90 μl was used, to which 10 μl of 50 μM stock of trypan blue was added to achieve an internal standard concentration of 5 μM. This was incubated at room temperature for 10 min to permit metabolite interaction with binding proteins, then extracted with 4 volumes (400 μl) of pre-chilled methanol:acetonitrile (50:50) to produce a final concentration of 40:40:20 (methanol:acetonitrile:H2O) and precipitated on ice for 10 min. Other steps were the same as for solid tissue extraction.

Suramin was analyzed on an AB SCIEX QTRAP 5500 triple quadrupole mass spectrometer equipped with a Turbo V electrospray ionization source, Shimadzu LC-20A UHPLC system, and a PAL CTC autosampler (AB SCIEX, Framingham, MA, USA). Ten microliters of extract were injected onto a Kinetix pentafluorophenyl column (150 × 2.1 mm, 2.6 μm; Phenomenex, Torrance, CA, USA) held at 30 °C for chromatographic separation. The mobile phase A was water with 20 mM ammonium acetate (NH4OAC; pH 7) and mobile phase B was methanol with 20 mM NH4OAC (pH 7). Elution was performed using the following gradient: 0 min—0% B, 15 min—100% B, 18 min—100% B, 18.1 min—0% B, 23 min—end. The flow rate was 300 μl min−1. All the samples were kept at 4 °C during analysis. Suramin and trypan blue were detected using scheduled multiple reaction monitoring (MRM) with a dwell time of 30 ms in negative mode and retention time window of 7.5–8.5 min for suramin and 8.4–9.4 min for trypan blue. MRM transitions for the doubly charged form of suramin were 647.0 m z−1 (Q1) precursor and 382.0 m z−1 (Q3) product. MRM transitions for trypan blue were 435.2 (Q1) and 185.0 (Q3). Absolute concentrations of suramin were determined for each tissue using a tissue-specific standard curve to account for matrix effects, and the peak area ratio of suramin to the internal standard trypan blue. The declustering potential, collision energy, entrance potential and collision exit potential were −104, −9.5,−32 and −16.9, and −144.58, −7, −57.8 and −20.94 for suramin and trypan blue, respectively. The electrospray ionization source parameters were set as follows: source temperature 500 °C; curtain gas 30; ion source gas 1, 35; ion source gas 2 35; spray voltage −4500 V. Analyst 1.6.1 was used for data acquisition and analysis. N=4–6 per tissue. Results are reported as means±s.e.m. in absolute μM (pmol μl−1) concentration for plasma, and pmol mg−1 wet weight for tissues.

Metabolomics

Broad-spectrum analysis of 478 targeted metabolites from 44 biochemical pathways in the plasma was performed by a modification of the methods described by Bajad and Shulaev.[27] Only male animals that had been behaviorally evaluated were tested. Samples were analyzed on an AB SCIEX QTRAP 5500 triple quadrupole mass spectrometer equipped with a Turbo V electrospray ionization source, Shimadzu LC-20A UHPLC system and a PAL CTC autosampler (AB SCIEX). Whole blood was collected 2 daysafter a single dose of suramin (20 mg kg−1 i.p.) or saline (5 μl g−1 i.p.) from animals that were lightly anesthetized with isoflurane (Med-Vet International, Mettawa, IL, USA, Cat no. RXISO-250) in a drop jar into BD Microtainer tubes containing lithium heparin (Becton Dickinson, San Diego, CA, USA, Ref no. 365971) by submandibular vein lancet.[28] Plasma was separated by centrifugation at 600 g × 5 min at 20 °C within 1 h of collection. Fresh lithium–heparin plasma was transferred to labeled tubes for storage at −80 °C for analysis. Typically, 45 μl of plasma was thawed on ice and transferred to a 1.7-ml Eppendorf tube. Two and one-half (2.5) microliters of a cocktail containing 35 commercial stable isotope internal standards (Supplementary Table S3) and 2.5 μl of 310 stable isotope internal standards that were custom-synthesized in Escherichia coli and Saccharomyces cerevisiae by metabolic labeling with 13C-glucose and 13C-bicarbonate were added, mixed and incubated for 10 min at 20 °C to permit small molecules and vitamins in the internal standards to associate with plasma-binding proteins. Macromolecules (protein, DNA, RNA and so on) were precipitated by extraction with 4 volumes (200 μl) of cold (−20 °C), acetonitrile:methanol (50:50) (LCMS grade, Cat no. LC015-2.5 and GC230-4, Burdick & Jackson, Honeywell, Muskegon, MI, USA), vortexed vigorously and incubated on crushed ice for 10 min, and then removed with centrifugation at 16 000 g × 10 min at 4 °C. The supernatants containing the extracted metabolites and internal standards in the resulting 40:40:20 solvent mix of acetonitrile:methanol:water were transferred to labeled cryotubes and stored at −80 °C for LC-MS/MS (liquid chromatography-tandem mass spectrometry) analysis.

LC-MS/MS analysis was performed by MRM under the Analyst v1.6.1 software control in both negative and positive modes with rapid polarity switching (50 ms). Nitrogen was used for curtain gas (set to 30), collision gas (set to high) and ion source gases 1 and 2 (set to 35). The source temperature was 500 °C. Spray voltage was set to −4500 V in negative mode and to 5500 V in positive mode. The values for Q1 and Q3 mass-to-charge ratios (m z−1), declustering potential, entrance potential, collision energy and collision cell exit potential were determined and optimized for each MRM for each metabolite. Ten microliters of extract were injected with PAL CTC autosampler into a 250 mm × 2.1 mm, 5-μm Luna NH2 aminopropyl HPLC column (Phenomenex) held at 25 °C for chromatographic separation. The mobile phase was solvent A: 95% water with 23.18 mM NH4OH (Sigma, Fluka Cat no. 17837-100ML), 20 mM formic acid (Sigma, Fluka Cat no. 09676-100ML) and 5% acetonitrile (pH 9.44); solvent B: 100% acetonitrile. Separation was achieved using the following gradient: 0 min—95% B, 4 min—B, 19min—2% B, 22 min—2% B, 23 min—95% B, 28 min—end. The flow rate was 300 μl min−1. All the samples were kept at 4 °C during analysis. The chromatographic peaks were identified using MultiQuant v2.1.1 (AB SCIEX), confirmed by manual inspection and the peak areas were integrated. The median of the peak area of stable isotope internal standards was calculated and used for the normalization of metabolite concentration across the samples and batches. N=6, 6.5-month-old males per group. Metabolite data were log-transformed before multivariate and univariate analyses.

Data analysis

Animals were randomized into active (suramin) and mock (saline) treatment groups at ~6 months of age. Group means and s.e.m. are reported. Behavioral data involving more than two groups were analyzed by two-way analysis of variance (ANOVA) and one-way ANOVAs (GraphPad Prism 5.0d, GraphPad Software Inc., La Jolla, CA, USA). Pair-wise post hoc testing was performed by the method of Tukey. Repeated measures ANOVA with prenatal treatment and drug as between subject factors and stimulus (mouse/cup) on time spent with mouse or cup was used as an additional test of social preference. Student's t-test was used for comparisons involving the two groups. Significance was set at P<0.05. Bonferroni post hoc correction was used to control for multiple hypothesis testing when t-tests were used to test social preference in two or more experimental groups. Metabolomic data were analyzed using multivariate partial least squares discriminant analysis, Ward hierarchical clustering and univariate one-way ANOVA with pairwise comparisons and post hoc correction by Fisher's least significant difference test in MetaboAnalyst.[29]

Results

Restoration of normal social behavior

Social behavior in mice can be quantified as the time spent interacting with a novel (‘stranger') mouse compared with the total time spent interacting with either a mouse or a novel inanimate object.[30] MIA animals showed social deficits from an early age (Figure 1a, Supplementary Figure S1a). Single-dose APT with suramin completely reversed the social abnormalities in 6.5-month-old adults (Figure 1b, Supplementary Figure S1b). Five weeks (5 half-lives) after suramin washout, a small residual benefit to social behavior was still detectable (Figure 1c, Supplementary Figure S1c). The residual social benefit of APT even after 5 weeks following suramin was correlated with retained metabolomic benefits (see below). This phenomenon was not investigated further but may be related to the development of metabolic memory and/or somatic epigenetic DNA changes that lasted longer than the physical presence of the drug.[31]

Figure 1

Single-dose correction of behavioral abnormalities. (a) Social abnormalities in male MIA animals were found at the earliest ages of testing at 2.25 months of age. (Student's t-test ****P<0.0002; N=19 Saline and 25 Poly(IC)). (b) A single dose of suramin given to adult MIA mice restored normal social behavior (PIC-Sur). two-way ANOVA was first used to test for the presence of interaction between drug treatment and experimental groups. This revealed an interaction consistent with the observation that suramin benefited social behavior in the MIA animals but had no effect on normal controls (F(1,39)=13.48); P=0.0007). We then performed one-way ANOVA to test for suramin effects. A single treatment with suramin (20 mg kg−1 i.p.) given 2–4 days before testing restored normal social behavior (one-way ANOVA F(3,40)=8.95; P<0.0001; Tukey post hoc PIC-Sal versus PIC-Sur ****P<0.0001; N=8–13 per group). (c) After 5 weeks of suramin washout, the social behavior remained improved compared with saline-treated animals but was decreased from the first week after treatment. (F(3,40)=10.5; Tukey post hoc PIC-Sal versus PIC-Sur *P<0.05; N=8–13 per group). Values are expressed as means±s.e.m. (d) We estimated the strength of novelty preference[32] as spontaneous alternation in the T-maze. MIA mice showed deficits in spontaneous alternation from the age of earliest testing at 4 months of age (Student's t-test; ****P<0.0001; N=19 Saline and 25 PIC). (e) Two-way ANOVA was first used to test for the presence of interaction between drug treatment and experimental groups. This revealed an interaction consistent with the observation that suramin restored spontaneous alternation in the MIA animals but had no effect on normal controls (F(1,40)=7.609; P=0.0087). We then performed one-way ANOVA to test for suramin effects. A single dose of suramin (20 mg kg−1 i.p.) injected 2–4 days before testing corrected the deficits in young adult animals that were 5.25 months of age. (F(3,40)=9.46; ; Tukey post hoc Sal-Sal versus PIC-Sal **P<0.01; PIC-Sal versus PIC-Sur ***P<0.001); N=8–13 per group). (f) This benefit was lost after a drug washout period of 5 weeks, leaving a significant difference between control (Sal) and MIA (PIC) groups (F(3,39)=18.05; P<0.0001), but no remaining effect of suramin by post hoc testing. (Tukey post hoc PIC-Sal versus PIC-Sur P=ns; N=8–13 per group). Values are expressed as means±s.e.m. (g) Motor coordination abnormalities were quantified on the rotarod as latency to fall. Performance was abnormal from the earliest age of testing at 2.5 months of age (Student's t-test ****P<0.0001; N=19 Saline and 25 Poly(IC)). (h) Suramin did not improve performance after two doses (20 mg kg−1 i.p.) given at 6.5 and 6.75 months of age and tested 2–4 days after the second dose. (two-way ANOVA interaction F(1,39)=0.1227; P=0.728 (ns); Poly(IC) effect F(1,39)=25.06; ****P<0.0001; treatment effect F(1,39)=0.01; P=0.908 (ns)). Values are expressed as means±s.e.m.

Restoration of spontaneous alternation in the T-maze

Novelty preference is an innate feature of normal rodent[32] and human[33] behavior and a predictor of socialization and communication growth in children with ASD.[34] The loss or suppression of novelty preference in children with ASD is associated with the phenomenon known as insistence on sameness.[35] We estimated preference for novelty as spontaneous alternation behavior in the T-maze.[25] The T-maze can also be used to estimate spatial working memory, especially when food-motivated.[36] We did not use the food-motivated variation in our study. We found that MIA animals showed deficient novelty preference as reflected by chance (near 50%) spontaneous alternation behavior (Figure 1d). These deficits were normalized after a single dose of suramin (Figure 1e). Five weeks after suramin washout, no residual benefit remained (Figure 1f).

Failure to restore rotarod performance in adult animals

Previous studies have shown age-dependent, postnatal loss of cerebellar Purkinje cells in the MIA model. This can reach up to 60% of Purkinje cells lost by 4 months (16 weeks) of age.[23,37] Motor coordination measured by rotarod performance is deficient in the MIA model[23] (and Figure 1g) and is critically dependent on the integrity of Purkinje cell circuits in the cerebellum.[38] We hypothesized that since Purkinje cells are known to be lost in MIA animals by 4 months (16 weeks) of age, that APT given later in life would have no effect. Our results confirmed this. We found that a single injection of suramin given to 6-month-old adults failed to restore normal motor coordination (Figures 1g and h). Although cerebellar Purkinje cell density was not quantified in this study, our results are consistent with the notion that once Purkinje cells are lost, their function cannot be restored by APT in adult animals.

Other behaviors

Certain features of ASD and schizophrenia were not captured by our studies of the MIA model and therefore could not be interrogated for pharmacologic response to APT. For example, we did not find any abnormalities in the MIA mouse model using our protocol when we looked for certain types of anxiety-related behavior in the light–dark box paradigm (Supplementary Figure S2). Likewise, our earlier studies showed no stereotypic repetitive movements in the C57BL/6J mouse strain either by clinical observation or by testing in the hole board exploration beam break mouse behavioral pattern monitor.[23] Finally, no abnormal behaviors were produced by suramin treatment itself. This was shown by the absence of behavioral differences between control mice treated with saline (‘Sal-Sal') and those treated with suramin (‘Sal-Sur') in Figures 1b, e and h.

Brainstem uptake of suramin

Suramin is known not to pass the blood–brain barrier;[39] however, no studies have looked at suramin concentrations in areas of the brain similar to the area postrema in the brainstem that lack a blood–brain barrier.[40] After completing the behavioral studies described above, we used mass spectrometry to measure drug levels in plasma, cerebrum, cerebellum and brainstem following a 5-week period of drug washout. The plasma half-life of suramin after a single dose in mice is 1 week.[41] No suramin was detected in any tissue after 5 weeks of drug washout (data not shown). We next gave an acute injection of suramin (20 mg kg−1 i.p.) to the remaining subjects. After 2 days, plasma suramin was 7.64 μM±0.50, and brainstem suramin was 5.15 pmol mg−1±0.49 (Figure 2). No drug was detectable in the cerebrum or cerebellum (<0.10 pmol mg−1 wet weight) in either control (Sal-Sur) or MIA (PIC-Sur) animals, consistent with an intact blood–brain barrier that excluded suramin from these tissues. In contrast to the cerebrum and cerebellum, the brainstem showed significant suramin uptake (Figure 2). These results are consistent with the notion that nuclei in brainstem, or their projection targets in distant sites of the brain, may mediate the dramatic behavioral effects of acute and chronic APT in this model.[23]

Figure 2

Plasma and brainstem suramin quantitation. (a) Suramin is a polysulphonated napthylurea with a molecular weight of 1297 g mol−1. (b) Suramin was present in the plasma and brainstem but was not detectable in the cerebrum or cerebellum. Two days after a single 20 mg kg−1 i.p. dose of suramin, drug levels were measured in plasma (7.64 μM±0.50), brainstem (5.15 pmol mg−1±0.49), cerebrum (<0.1 pmol mg−1) and cerebellum (<0.1 pmol mg−1) in both controls (Sal) and maternal immune activation (PIC) animals. (N=4–6 per tissue). Values are expressed as means±s.e.m.

Restoration of normal purine metabolism rescues other metabolic disturbances

We analyzed the acute metabolomic effects in plasma 2 days after single-dose treatment with suramin or saline in the same animals studied behaviorally. We measured 478 metabolites from 44 pathways using mass spectrometry, analyzed the data by partial least squares discriminant analysis and visualized the results by projection in two dimensions (Figures 3a and b). This revealed sharp differences between control and MIA animals that were substantially normalized by a single treatment with suramin (Figure 3a). Figure 3b shows a similar analysis that illustrates the gradual return to disease-associated metabolism after 5 weeks of drug washout. Using hierarchical cluster analysis we found that the metabolic profiles of controls (Sal-Sal; light blue) and MIA animals that were treated with one dose of suramin (PIC-Sur; green) were more similar (major branch on the left of Figure 3c) than the metabolic profiles of saline-treated MIA animals (PIC-Sal; red) and the MIA animals tested 5 weeks after suramin washout (PIC-Sur W/O; dark blue; major branch on the right of Figure 3c). The reason that the metabolic profile had not returned completely to pretreatment conditions (to the position of the red triangles in Figure 3b) even after 5 weeks following a dose of suramin was not investigated but could be due to the development of metabolic memory and/or somatic epigenetic DNA changes that lasted longer than the physical presence of the drug.[31]

Figure 3

Metabolomic analysis. (a) APT rescues widespread metabolic abnormalities. Plasma samples were collected 2 days after a single dose of suramin (20 mg kg−1 i.p.) or saline (5 μl g−1 i.p.). This analysis shows that a single dose of suramin (PIC-Sur; green) drives the metabolism of MIA animals (PIC-Sal; red) strongly in the direction of controls (Sal-Sal; blue). Metabolomic profiles consisted of 478 metabolites from 44 biochemical pathways measured with LC-MS/MS. N=6, 6.5-month-old males per group. (b) Metabolic memory preserves metabolic rescue by APT. This analysis shows that 5 weeks after a single dose of suramin (PIC-Sur W/O; green) the metabolism of treated animals has drifted back toward that of untreated, MIA animals (PIC-Sal; red; N=6 males per group). (c) Hierarchical clustering of suramin-treated and suramin-washout metabotypes. This analysis illustrates the metabolic similarity of control (Sal-Sal; light blue) and MIA animals treated with one dose of suramin (PIC-Sur; green) compared with saline-treated MIA animals (PIC-Sal; red) and ASD-like animals tested 5 weeks after suramin washout (PIC-Sur W/O; dark blue). The numbers listed along the x axis are animal ID numbers. (d) Rank Order of metabolites disturbed in the MIA model. Multivariate analysis across the four treatment groups (PIC-Sal=MIA; PIC-Sur=acute suramin treatment; PIC-Sur w/o=5 weeks post-suramin washout; Sal-Sal=Controls). Biochemical pathway assignments are listed on the left. Relative magnitudes of each metabolite disturbance are listed on the right as high (red), intermediate (yellow or light green) and low (dark green). Variable importance in projection (VIP) scores are a multivariate statistic that reflects the impact of each metabolite on the partial least squares discriminant analysis model. VIP scores above 1.5 are significant.

Figure 3d shows the top 48 significant metabolites found in the untreated MIA animals, ranked according to their impact by variable importance in projection (VIP) score. The color-coded columns on the right of the figure indicate the direction of the change. In 43 of the 48 (90%) discriminating metabolites, suramin treatment (PIC-Sur) resulted in a metabolic shift in concentration that was either intermediate (coded yellow or light green) or in the direction of and beyond that found in control animals (Sal-Sal). The biochemical pathways represented by each metabolite are indicated on the left of Figure 3d.

Metabolic pathway analysis

The most influenced biochemical pathway in the MIA mouse was purine metabolism (Table 1). Eleven (23%) of the 48 discriminant metabolites were purines. Nine (82%) of the 11 purine metabolites were increased in the untreated MIA mice, consistent with hyperpurinergia. Only ATP and allantoin, the end product of purine metabolism in mice, were decreased in the plasma. A limitation of plasma metabolomics is that it cannot measure the effective concentration of nucleotides in the pericellular halo that defines the unstirred water layer near the cell surface where receptors and ligands meet.[23] The concentration of ATP in the unstirred water layer is regulated according to conditions of cell health and danger[22] in the range of 1–10 μM, which is near the EC50 of most purinergic receptors.[42] This is up to 1000-fold more concentrated than the 10–20 nM levels of ATP in compartments removed from the cell surface such as the plasma.[43] In the plasma we found that suramin restored 9 (82%) of the 11 purine metabolites to more normal levels, including ATP and allantoin (Figure 3d, right PIC-Sur column, coded yellow or light green) and increased inosine and deoxyinosine to above normal.

Table 1

Biochemical pathways with metabolites altered in the MIA mouse model of neurodevelopmental disorders

No.	Pathway	Measured metabolites in the pathway (N)	Expected pathway proportion (P=N/478)	Expected hits in a sample of 48 (P*48)	Observed hits in the top 48 metabolites	Fold- enrichment (Obs/Exp)	Impact (Σvip)	Fraction of VIP explained (% of 116.16)	Pathway normalized by single-dose suramin treatment
1	Purine metabolism	48	0.1004	4.8201	11	2.3	28.19	24.3%	Yes (9/11)
2	Microbiome metabolism	32	0.0669	3.2134	6	1.9	17.53	15.1%	Yes (6/6)
3	Phospholipid metabolism	88	0.1841	8.8368	4	0.5	9.76	8.4%	Yes (4/4)
4	Bile salt metabolism	4	0.0084	0.4017	3	7.5	9.23	7.9%	No (0/3)
5	Sphingolipid metabolism	72	0.1506	7.2301	4	0.6	8.28	7.1%	Yes (4/4)
6	Cholesterol, cortisol, steroid metabolism	19	0.0397	1.9079	4	2.1	8.08	7.0%	Yes (4/4)
7	Glycolysis and gluconeogenesis	17	0.0356	1.7071	3	1.8	6.25	5.4%	Yes (3/3)
8	Oxalate, glyxoylate metabolism	3	0.0063	0.3013	2	6.6	5.02	4.3%	Yes (2/2)
9	Tryptophan metabolism	11	0.0230	1.1046	1	0.9	4.11	3.5%	Yes (1/1)
10	Krebs cycle	18	0.0377	1.8075	2	1.1	3.58	3.1%	Yes (2/2)
11	Vitamin B3 (niacin/NAD) metabolism	7	0.0146	0.7029	1	1.4	3.19	2.7%	Yes (1/1)
12	GABA, glutamate, arginine, ornithine, proline metabolism	6	0.0126	0.6025	1	1.7	2.33	2.0%	Yes (1/1)
13	Pyrimidine metabolism	35	0.0732	3.5146	1	0.3	2.24	1.9%	Yes (1/1)
14	Vitamin B2 (riboflavin) metabolism	4	0.0084	0.4017	1	2.5	1.97	1.7%	Yes (1/1)
15	Thyroxine metabolism	1	0.0021	0.1004	1	10.0	1.66	1.4%	Yes (1/1)
16	Amino-sugar and galactose metabolism	10	0.0209	1.0042	1	1.0	1.61	1.4%	Yes (1/1)
17	SAM, SAH, methionine, cysteine, glutathione metabolism	22	0.0460	2.2092	1	0.5	1.57	1.3%	Yes (1/1)
18	Biopterin, neopterin, molybdopterin metabolism	1	0.0021	0.1004	1	10.0	1.56	1.3%	Yes (1/1)
		398 (0.8326 × 478)	0.8326	40 (0.8326 × 48)	48		116.16	100%	94% (17/18)

Abbreviation: VIP, variable importance in projection.

Pathways were ranked by their impact measured by summed VIP (ΣVIP) scores. A total of 48 metabolites were found to discriminate treatment, control, washout and MIA groups by multivariate partial least squares discriminant analysis (PLSDA). Significant metabolites had VIP scores of ⩾1.5. Eighteen (41%) of the 44 pathways interrogated had at least one metabolite with VIP scores ⩾1.5. The total impact of these 48 metabolites corresponded to a summed VIP score of 116.16. The fractional impact of each pathway is quantified as the percent of the summed VIP score and displayed in the final column on the right in the table. Single dose APT with suramin not only corrected purine metabolism but also normalized 17 (94%) of 18 metabolic pathway abnormalities that defined the MIA model of neurodevelopmental disorders.

Additional pathway analysis revealed a pattern of disturbances that was remarkably similar to metabolic disturbances that have been found in children with ASDs (Table 1). Eighteen of the 44 pathways were disturbed in the MIA model. The 44 pathways interrogated by this analysis are reported in Supplementary Table S1. After purine metabolism, the next most influenced pathway was the microbiome. Microbiome metabolites are molecules that are produced by biochemical pathways that are absent in mammalian cells but are present in bacteria that reside in the gut microbiome.[44] Together, purine and microbiome metabolism accounted for nearly 40% (ΣVIP=39.4%) of the impact measured by VIP scores. The two top discriminant metabolites were products of the microbiome (Figure 3d). A total of seven pathways each contributed 5% or more to the VIP pathway impact scores (Table 1). These top seven pathways were purines, microbiome metabolism, phospholipids, bile salt metabolism, sphingolipids, cholesterol, cortisol, and steroid metabolism and glycolysis. Seventy-five percent (75%) of the metabolite VIP score impact was accounted for by metabolites in these seven pathways (Table 1). Univariate statistical analysis was conducted by one-way ANOVA and pairwise group comparisons with post hoc correction (Supplementary Table S2). Forty-six (46) metabolites satisfied a false discovery rate threshold of less than 10% in this analysis. These were rank ordered by P-values. This univariate analysis identified 16 (35% of 46) metabolites (shaded yellow, Supplementary Table S3) that were also found by multivariate analysis across the four groups, and 30 (65%) additional metabolites (unshaded in Supplementary Table S2) that were discriminating only in pairwise group comparisons.

Restoration of normal purine metabolism by APT led to the concerted normalization of 17 (94%) of the 18 biochemical pathway disturbances that characterized the MIA model (Table 1; far right column). Only the bile salt pathway was not restored by suramin (Table 1,Figure 3d). The three bile salt metabolites were highest in the plasma of control animals (Figure 3d; Sal-Sal coded red in the columns on the right), lower in MIA animals (Figure 3d; PIC-Sal coded yellow) and made even lower by suramin (Figure 3d; PIC-Sur, coded dark green). Overall, we found that restoration of normal purine metabolism with APT led to the concerted improvement in both the behavioral and metabolic abnormalities in this model.

Discussion

Children with inborn errors in purine and pyrimidine metabolism have long been known to have neurodevelopmental and behavioral abnormalities.[45] However, the neurochemical basis for the brain and behavioral manifestations of these classic disorders such as Lesch-Nyhan syndrome[46] and adenosine deaminase deficiency[45] is not yet understood. Genetic disorders of purine metabolism were some of the first single-gene disorders found to be associated with ASD behaviors. In 1969, William Nyhan described a 3-year-old boy with autistic behavior and hearing impairment caused by a defect in the first committed step in purine biosynthesis, phosphoribosyl pyrophosphate synthase.[47] This mutation created an enzyme that was resistant to feedback inhibition by ATP, resulting in phosphoribosyl pyrophosphate synthase superactivity, and excess purine biosynthesis.[47, 48, 49] Other reports of purine[50] and pyrimidine[51] disorders linked to autism soon followed.

The discovery of ATP signaling by Geoffrey Burnstock just a few years later in 1972 showed for the first time that extracellular nucleotides could act as neurotransmitters.[52] However, the fields of inborn errors in purine and pyrimidine metabolism and purinergic signaling developed independently for the next 40 years. Our work in autism lies at the confluence of these two rivers of investigation. A large number of other single gene disorders,[53] genetic variants, chromosomal and copy number variations[2,54] and environmental factors[6,55] are known to increase the risk of ASD. Most of these have no apparent ties to purine and pyrimidine metabolism. Our results are consistent with the unifying hypothesis that each factor, or mixture of factors, that causes ASD, does so by triggering a conserved cellular response to stress that we have called the CDR.[22] In the MIA mouse model, the CDR is maintained by abnormalities in purinergic signaling that can be treated with antipurinergic drugs such as suramin.

It is well known that ATP inside the cell acts as an energy carrier. It is less well known that ATP and other nucleotides located outside the cell can bind to cell surface receptors and act as signaling molecules and neuromodulators that are important in inflammation,[56] neurotransmission[57] and many other biological processes. As such, ATP, other nucleotides and certain metabolites have been collectively considered as mitokines—signaling molecules traceable to mitochondrial function.[23] Fifteen different isoforms of phosphorylated purinergic nucleotide receptors are now known.[58] Four additional purinergic receptors are responsive to the unphosphorylated nucleoside adenosine, and called the P1, or adenosine receptors.[59] The 15 phosphorylated nucleotide receptors are divided into seven ionotropic P2X1-7 receptors, and eight metabotropic, P2Y G-protein-coupled receptors numbered 1–14. Six former P2Y receptors (nos 3, 5 and 7–10) have been reclassified. Purinergic receptors are expressed in tissue-specific patterns on every cell type in the body. Together, the 19 P1, P2X and P2Y receptors control a broad range of biological functions that have relevance to a number of neurodevelopmental disorders, including autism, schizophrenia, attention deficit/hyperactivity and bipolar disorder. These include synaptogenesis and brain development,[58] neuronal plasticity,[60,61] sleep,[62] regulation of the PI3K/AKT/GSK3β (glycogen synthase kinase-3β) pathway,[63] control of immune responses and chronic inflammation,[64] gut motility,[65] gut permeability,[66] taste chemosensory transduction,[67] sensitivity to food allergies,[68] hearing,[69] innate immune signaling,[70] neuroinflammation, antiviral signaling, microglial activation, neutrophil chemotaxis, autophagy, chronic pain syndromes,[58] cerebellar Purkinje cell signaling and motor coordination,[71] and the regulation of autonomic parasympathetic control of heart rate variability.[72] The importance of metabolism in regulating behavior has recently been highlighted by the discovery that acute inhibition of the metabolic control protein GSK3β restores normal spontaneous alternation in the Y-maze in the MIA model.[73]

Our results show that purine metabolism is a master regulatory pathway in the MIA model (Table 1, Figure 3d, Supplementary Table S1). Correction of purine metabolism with APT restored normal social behavior (Figure 1b) and novelty preference (Figure 1e). Comprehensive metabolomic analysis revealed disturbances in several other metabolic pathways relevant to children with ASDs. These included disturbances in microbiome, phospholipid, cholesterol/sterol, sphingolipid, glycolytic and bile salt metabolism (Table 1). The top, non-microbiome-associated metabolite was quinolinic acid (Figure 3d), which was decreased in the MIA model. Quinolinic acid is a product of the indoleamine 2,3-dioxygenase pathway of tryptophan metabolism.[74] Interestingly, abnormalities in purine,[47,75] tryptophan,[76,77] microbiome,[78,79] phospholipid,[80] cholesterol/sterol[81] and sphingolipid[82,83] metabolism have each been reported in children with ASDs. Abnormalities in purine metabolism,[84] tryptophan,[85] cholesterol/sterol,[86] sphingolipid[87] and phospholipid[88] metabolism have also been described in schizophrenia. Although the detailed metabolic features of ASD and schizophrenia are different, these disorders share biochemical pathway disturbances that reveal the persistent activation of the evolutionarily conserved CDR[22] in both ASD and schizophrenia. These data show that the metabolic disturbances in the MIA model and human ASD and schizophrenia are similar and provide strong support for the biochemical validity of this animal model.

Our results also reveal a novel mechanism by which drugs that cannot penetrate the blood–brain barrier can still have central nervous system effects. We show that suramin is taken up into the brainstem (Figure 2b), an area of the brain known to contain critical collections of chemosensory neurons such as the area postrema that is not protected by the blood–brain barrier.[40] Purinergic neurons in the area postrema are known to integrate blood-borne signals that activate the hypothalamic–pituitary–adrenal axis during stress,[89] coordinate parasympathetic and sympathetic autonomic balance[90] and regulate whole-body metabolism and sickness behavior during inflammation[91] in response to the evolutionarily conserved CDR associated with ASDs and other disorders.[22] These data provide a plausible new mechanism that could apply to many drugs that have central nervous system effects but are known not to pass the blood–brain barrier.

Although our results of single-dose correction of abnormal behaviors in an animal model of autism and schizophrenia are encouraging, there are several caveats that must be considered before extending the results to humans. First, while the MIA mouse model captures several features of ASD and schizophrenia, no animal model can fully capture the complexities of human behavior. Second, suramin is a poor drug choice for chronic use because of potentially toxic side effects that can occur with prolonged treatment.[92] Third, human forms of ASD and schizophrenia may occur by mechanisms not captured by the MIA model. Mechanisms that do not involve the CDR[22] may not be amenable to APT. Human clinical trials will be necessary to answer these questions.

In summary, this study identifies purinergic signaling as a novel neurochemical switch that regulates both behavior and metabolism in the MIA mouse model. Our results provide new tools for refining current concepts of pathogenesis in autism, schizophrenia and several other neurodevelopmental disorders, and create a fresh path forward for the development of newer and safer drugs. The data provide preclinical support for the hypothesis that some of this new class of antipurinergic medicines need not be given chronically. Rather, new drugs might be given only once, or intermittently, during sensitive windows to unblock metabolism, restore more normal neural network function, improve resilience and plasticity, and permit improved development in response to behavioral and interdisciplinary therapies, and to natural play.