1744-8069-8-25 1744-8069 Research <p>Serotonin-Induced Hypersensitivity via Inhibition of Catechol O-Methyltransferase Activity</p> TsaoDouglastsaodc11@wfu.edu WieskopfSJeffreyjeffrey.wieskopf@mail.mcgill.ca RashidNaimnaim@unc.edu SorgeERobertrobert.sorge@mail.mcgill.ca RedlerLRachelrachelredler@gmail.com SegallKSamanthasegall@DENTISTRY.UNC.EDU MogilSJeffreyjeffrey.mogil@mcgill.ca MaixnerWilliammaixnerw@DENTISTRY.UNC.EDU DokholyanVNikolaydokh@med.unc.edu DiatchenkoLudaluda@unc.edu

Department of Chemistry University of North Carolina, Chapel Hill, NC 27599, USA

Center for Neurosensory Disorders, University of North Carolina, Chapel Hill, NC 27599, USA

Department of Psychology, McGill University, 1205 Dr. Penfield Avenue, Montreal, QC, H3A 1B1, Canada

Department of Biostatistics Gillings School of Global Public Health University of North Carolina, Chapel Hill, NC 27599, USA

Department of Biochemistry and Biophysics, School of Medicine University of North Carolina, Chapel Hill, NC 27599, USA

Molecular Pain 1744-8069 2012 8 1 25 http://www.molecularpain.com/content/8/1/25 10.1186/1744-8069-8-2522500608
121201227320121342012 2012Tsao et al.; licensee BioMed Central Ltd.This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The subcutaneous and systemic injection of serotonin reduces cutaneous and visceral pain thresholds and increases responses to noxious stimuli. Different subtypes of 5-hydroxytryptamine (5-HT) receptors are suggested to be associated with different types of pain responses. Here we show that serotonin also inhibits catechol O-methyltransferase (COMT), an enzyme that contributes to modultion the perception of pain, via non-competitive binding to the site bound by catechol substrates with a binding affinity comparable to the binding affinity of catechol itself (K i  = 44 μM). Using computational modeling, biochemical tests and cellular assays we show that serotonin actively competes with the methyl donor S-adenosyl-L-methionine (SAM) within the catalytic site. Binding of serotonin to the catalytic site inhibits the access of SAM, thus preventing methylation of COMT substrates. The results of in vivo animal studies show that serotonin-induced pain hypersensitivity in mice is reduced by either SAM pretreatment or by the combined administration of selective antagonists for β2- and β3-adrenergic receptors, which have been previously shown to mediate COMT-dependent pain signaling. Our results suggest that inhibition of COMT via serotonin binding contributes to pain hypersensitivity, providing additional strategies for the treatment of clinical pain conditions.

Introduction

Catechol O-methyltransferase (COMT) is an enzyme that has been implicated in the perception of mechanical, thermal, and inflammatory pain in both humans and rodents 1 2 3 . Its association with a broad range of noxious stimuli suggests its critical role as an underlying factor for the fundamental processes of pain perception. The enzymatic activity of COMT is reversely correlated with the intensity of perceived pain; higher levels of COMT activity correspond to lower pain sensitivity 1 2 3 4 . Thus, compounds or cellular regulatory factors that lower COMT activity enhance pain sensation.

While commonly associated with antinociception 5 6 , serotonin (5-hydroxytryptamine, or 5-HT) is also known to produce a hyperalgesic response when injected subcutaneously or into deep tissue 7 8 . Different subtypes of 5-HT receptors are associated with hypersensitive responses to various noxious stimuli. For example, mechanical hypersensitivity (allodynia) induced by serotonin has been linked to 5-HT1A, 5-HT2B, and 5-HT2C receptors 7 9 , whereas the 5-HT2A receptor appears to be responsible for 5-HT-induced thermal hypersensitivity 10 .

A co-factor necessary for methylation of COMT substrates is S-adenosyl-L-methionine (SAM). SAM is responsible for donating the methyl group to a single hydroxyl group of the catechol substrate. Without SAM, COMT is unable to metabolize substrates such as norepinephrine and epinephrine, thereby leading to a potentiation of pain signaling through the downstream stimulation of β2 – and β3–adrenergic receptor pathways 2 . In comparing the chemical structures of SAM and serotonin (Figure 1A), we noted that the indole ring of serotonin is similar to the adenosine motif of SAM. Thus, we hypothesized that serotonin may play a role as a competitive inhibitor of SAM, which could provide a novel avenue for COMT-dependent pain modulation. Specifically, the inhibition of COMT activity through serotonin might reduce the ability of COMT to methylate endogenous catecholamine substrates, including several neurotransmitters, such as dopamine, epinephrine, and norepinephrine 11 . Here, we investigate the binding properties of serotonin to COMT using a combination of in silico, in vitro, and behavioral studies.

<p>Figure 1</p>

In silico and in vitro binding of serotonin and COMT.

In silico and in vitro binding of serotonin and COMT. ( A) Structures of serotonin and S-adenosyl-L-methionine (SAM). ( B) SPR data is fitted to association and dissociation curves to determine the kinetics of serotonin binding ( kon = 0.0012 s-1, koff = 13.7404 M-1 s-1). (C and D) Structural modeling of serotonin within COMT active site. ( C) Surface representation of serotonin binding to a pocket within COMT. Serotonin structure is shown in cyan; SAM structure is shown in green. The overlapped structures indicate that the amine sidechain of serotonin prevents SAM from actually binding to COMT. On the COMT surface, blue represents positively-charged, red represents negatively-charged, and white represents neutrally-charged regions. ( D) Interactions of serotonin with residues inside COMT active site. COMT structure and residues are shown in green; serotonin is shown in blue. Hydrogen-bonding interactions are highlighted with yellow-dashed lines.

Results

Serotonin Binds to COMT

Binding interactions between serotonin and COMT were determined through surface plasmon resonance (Figure 1B) whereby both the rate of association (k on) and dissociation (k off) can be determined. We find that serotonin dissociates from COMT with a k off rate of 0.0012 s-1, and binds to COMT with a k on rate of 13.7404 M-1 s-1. The equilibrium dissociation constant (K D) was determined by dividing the dissociation rate from the association rate, and calculated to be 87 μM.

Serotonin Computationally Docks to the COMT Active Site

Next, we investigated whether serotonin could be accommodated within the catalytic site of COMT. Using computational docking tools 12 13 , we dock serotonin to an apo-COMT protein to predict its native pose in the active site. We find that serotonin occupies a site within the pocket that sterically overlaps with the binding site of SAM (Figure 1)C. In the crystal structure of holo-COMT, His142 participates in an edge-to-face interaction with the adenosine of SAM 14 . However, during serotonin binding, His142 donates a hydrogen bond with the amine side chain of serotonin (Figure 1D). Another hydrogen bonding interaction takes place with the serotonin amine and the side chain carboxyl group of Asp90. These two residues are responsible for positioning serotonin near the binding pocket of SAM, thereby preventing SAM from accessing its site. The carbonyl backbone of His142 also accepts a hydrogen bond from the hydroxyl group of serotonin.

Two more residues form hydrogen-bonding interactions with serotonin. The carbonyl of Met40 forms a hydrogen bond with the imidazole group of serotonin, while Lys144 forms a hydrogen bond with the hydroxyl group of serotonin. The side chain of Met40 also packs against the indole ring of serotonin. Overall, the results of this docking simulation suggest that serotonin may inhibit COMT activity by preventing the binding of the critical cofactor, SAM.

Serotonin Inhibits COMT Activity

To test whether the binding of serotonin leads to inhibition of COMT activity, we performed a set of kinetics experiments to measure the reaction rates in either the presence or the absence of serotonin (Figure 2A,B). We measured the variance in reaction velocity with respect to changes in catechol concentration. A double-reciprocal plot comparing COMT velocities in the presence or absence of serotonin shows that addition of 100 μM of serotonin decreases the reaction velocity of COMT (Figure 2C). Thus, serotonin binding to COMT inhibits COMT activity.

<p>Figure 2</p>

Inhibition of COMT with serotonin.

Inhibition of COMT with serotonin. ( A) and ( B) show the reaction rates of COMT at different catechol concentrations in the absence and presence of 100 μM serotonin, respectively. Legend: 0.5 μM (black), 1 μM (red), 3.3 μM (green), 5 μM (blue), 10 μM (orange), 33 μM catechol (brown). ( C) Lineweaver-Burk plot comparing the reciprocal velocities of activity against the reciprocal concentration of catechol in the absence (black) and presence (red) of 100 μM serotonin. The higher slope observed for serotonin indicates some level of inhibition. ( D) Michaelis-Menten determination of equilibrium inhibition constant. Activity assays in the absence (black) and presence (red) of 100 μM serotonin were fit to a two-site Michaelis-Menten model. While the KM for catechol remains relatively the same, the Vmax shows some noticeable differences thereby suggesting a non-competitive mechanism of inhibition between serotonin and catechol.

We examined the kinetics of the observed inhibition in a set of experiments using a two-site Michaelis-Menten model (see Materials and Methods; Figure 2D). The Michaelis-Menten equilibrium constants for catechol binding were determined when serotonin was absent or present (100 μM) to provide the calculated values or the apparent K M and K M,app (see Materials and Methods). A comparison of the calculated K M (16 μM) and K M,app (18.5 μM) values shows a negligible difference between their affinities towards the catechol substrate. However, a comparison of the maximum velocities reveals a greater difference. The V max in the absence of serotonin was 16 nM/min while the apparent V max,app in the presence of serotonin was 7 nM/min, suggesting that serotonin inhibition is non-competitive with respect to catechol.

Given the differences between V max and V max,app, we can calculate the inhibition constant for serotonin (apparent K i, or K i,app; see Materials and Methods). We find that K i,app is 44 μM, which is within the vicinity of our approximate K D determined via SPR (~87 μM). Overall, our results show that serotonin binds to COMT and inhibits its activity at a concentration that is similar to endogenous COMT substrates (~0.5 M intravesicular concentration; 25 μM reported for dopamine and ~10 μM for norepinephrine) 15 16 17 18 19 .

Serotonin-Induced Hypersensitivity is Diminished by Pretreatment with SAM or β2/3 Antagonists

If serotonin inhibition of COMT occurs through the competitive binding of the active site with SAM to produce pain hypersensitivity, we would predict that the administration of SAM should antagonize serotonin-induced hypersensitivity. We tested this hypothesis by performing a series of standard in vivo nociceptive assays in mice featuring the intraplantar (i.pl.) administration of serotonin 9 . We chose to use mechanical allodynia to assess the 5-HT-COMT interaction over the thermal pain, as thermal hyperalgesia is a far less important clinical symptom than mechanical allodynia, being observed in a much smaller percentage of chronic pain patients, being rated as much less bothersome, and having lower correlations with global pain ratings 20 21 22 .

Behavioral responsiveness to mechanical stimuli did not differ between groups prior to pharmacological manipulations, nor after SAM injection alone. I.pl. administration of serotonin (10 μM in a 20 μL injection volume) produced robust ipsilateral tactile hypersensitivity (i.e., mechanical allodynia) compared to saline (p < 0.001) (Figure 3A), consistent with previous observations 9 . No changes in mechanical withdrawal thresholds were observed in the contralateral hind paw (data not shown). Mice pretreated with SAM displayed modestly but significantly less serotonin-induced hypersensitivity compared to saline-pretreated mice (Figure 3A; p = 0.02).

<p>Figure 3</p>

Reduction of serotonin-induced hypersensitivity by SAM and β23 antagonism.

Reduction of serotonin-induced hypersensitivity by SAM and β23antagonism. In all graphs, symbols represent mean ± SEM withdrawal threshold (g). In both experiments, serotonin (5-HT) was injected into one hind paw and withdrawal thresholds were assessed over a 4–h period. ( A) Pretreatment with 80 mg/kg SAM (5-HT + SAM) lowers 5-HT-induced hypersensitivity compared to saline-pretreated mice (5-HT + Saline). SAM alone (Saline + SAM) produced no changes in withdrawal threshold at any time point. ( B) Administration of the selective β2 and β3 antagonists, ICI118,551 and SR59230A (5-HT + ICI/SR), reduces 5-HT-induced hypersensitivity compared to saline controls (5-HT + Saline). ICI/SR alone (Saline + ICI/SR) produced no changes in withdrawal threshold at any time point.

As an additional approach to examine whether a COMT is associated with serotonin-induced hypersensitivity, we determined whether the blockade of downstream receptors known to mediate COMT–dependent changes diminishes serotonin-induced hypersensitivity. Since we have previously shown that low COMT activity leads to increased pain sensitivity via a β2/3-adrenergic mechanism, blocking β2 and β3 adrenergic receptors should reduce the serotonin-induced hypersensitivity if that sensitivity is COMT-mediated 2 . Again, serotonin produced robust hypersensitivity compared to saline-treated controls (p < 0.001; Figure 3B). Consistent with our hypothesis, we find that relative to serotonin administration alone, co-administration of ICI118,551 and SR59230A—compounds that respectively block β2 and β3 adrenergic receptors—significantly reduces serotonin hypersensitivity (p = 0.01; Figure 3B).

Thus, the results of our behavioral experiments are consistent with our in vitro and in silico results, and suggest that at least a portion of serotonin-induced hypersensitivity is mediated by a COMT → β2/β3 adrenergic receptor mechanism. Furthermore, our results indicate that serotonin can actively compete with SAM for the catalytic site of COMT.

Discussion

Biochemical Interactions of Serotonin and COMT

Using binding and activity assays, we show that serotonin interacts with COMT in a manner that inhibits enzymatic activity. We estimate that the COMT’s KD for serotonin is approximately 87 μM, as determined with surface plasmon resonance experiments. Based on previous literature, our calculated K D of 87 μM for serotonin is a reasonable 15 16 17 18 19 . SAM has the lowest K D reported among COMT ligands, with a K D of 20.2 μM 17 . In contrast, many common COMT catechol substrates exhibited dissociation constants in the range of 40 – 400 μM 17 . The higher dissociation constants observed for COMT substrates is most likely a consequence of its low substrate specificity. The binding site for catechol is very solvent accessible and thus can accommodate a large array of ligands 14 . Although COMT does not methylate most phenols or gallols 23 , virtually all catechols are susceptible to methylation including bulkier steroid compounds. Thus, given a K D for serotonin that is of the same order of magnitude observed for other known substrates, the interaction between serotonin and COMT is very plausible within cells that contain COMT.

We also find that the activity of COMT is reduced in the presence of serotonin, with a K i of approximately 44 μM. Furthermore, our data suggest that serotonin is a non-competitive inhibitor with respect to catechol substrates. These results are in line with our structural model because serotonin occupies a region within the active site that appears to hinder SAM binding and still allows catechol binding from a steric viewpoint. Moreover, serotonin does not possess a methyl group to donate that would substitute for SAM binding.

COMT Contributes to Serotonin-Induced Hypersensitivity

Our in silico models and in vitro results suggest that serotonin is an inhibitor of COMT. However, whether this inhibition has biological significance cannot be ascertained from these data alone. Our structural model suggests that serotonin inhibits COMT activity by actively competing with SAM at the active site. This mechanism is further supported by our kinetics studies. We performed behavioral experiments to determine if our in silico model and in vitro data are predictive of in vivo effects on pain behaviors. Mice pretreated with SAM displayed diminishes serotonin-induced mechanical hypersensitivity (Figure 3A). We hypothesize that the additional SAM produces this attenuation by decreasing the probability of serotonin occupying the active site of COMT. The relatively modest anti-allodynic effect is most likely due to the limited half-life of SAM 24 . The donating methyl group resides within a sulfonium center that is unstable due to its cationic charge. In addition, each SAM molecule can only donate one methyl group to a substrate. Once methylation occurs, SAM is converts into S-adenosyl-L-homocysteine (SAH). SAH can further serve as an inhibitor of COMT activity due to the high structural similarity to SAM but like serotonin, SAH is without a methyl donor 25 .

SAM is a cofactor that is responsible for many methylation reactions within the cell. The behavioral effects of SAM found presently could be non-specific, and not mediated only by COMT. Thus, we examined whether serotonin enhances pain sensitivity via the downstream stimulation of a β2/3-adrenergic receptors as occurs in response to the pharmacological inhibition of COMT 2 . Inhibiting COMT, and thereby preventing methylation of epinephrine and norepinephrine, promotes pain signaling through the stimulation of β2 and β3 receptors 2 . If serotonin inhibits COMT in vivo, we would thus expect that the observed increases in pain sensitivity results, at least in part, from the activation of β2/3 receptors. Consistent with this view, blockade of β2/3 receptors inhibited serotonin-induced pain hypersensitivity in mice (Figure 3B).

Implications for Serotonin-Induced Hypersensitivity and SSRIs

Our current results are of potentially substantial clinical significance. Increasing the bioavailability of serotonin with selective serotonin reuptake inhibitors (SSRIs) produces analgesia in some animal models 6 26 and can be used to treat clinical pain conditions 27 28 . However, SSRIs show relatively low analgesic efficacy in patients suffering from persistent pain conditions such as fibromyalgia, rheumatoid arthritis, and migraines 27 29 30 31 and have not gained widespread use for the treatment of persistent pain conditions 32 33 34 35 . In contrast, agents that have greater effects on the inhibition of norepinephrine relative serotonin reuptake (SNRIs), like duloxetine and milnacipran, show greater effect sizes and are more widely used for the treatment of a number of persistent pain states 36 37 38 .

The clinical analgesic effects of SSRIs are very modest and do not match the expectations seen in animal models. Furthermore, several animal studies note a dual role for serotonin in both analgesia and hypersensitivity, which appears to be dependent on the model and site of administration. Thus, the unraveling of the neural mechanisms that underlie the dual action of serotonin on pain perception is of importance. Mechanisms proposed to date that explain the dual action of serotonin primarily involve the activation of different subpopulations of 5-HT receptors that are distributed at different anatomical locations 7 8 9 10 . Our in silico, in vitro and in vivo results suggest that serotonin’s pain-promoting effect can also at least in part be attributed to serotonin-dependent COMT inhibition. This finding opens a new avenue for increasing the analgesic efficacy of SSRIs by co-administrating SAMe and/or non-selective beta-blockers like propranolol.

Conclusions

While several 5-HT receptors subtypes are known to contribute to pain perception, we have demonstrated that 5-HT action on COMT activity is another mechanism underlying 5-HT induced hypersensitivity, through a non-competitive binding process between SAM and serotonin at COMT’s catalytic site. Elucidation of the mechanisms of serotonin-induced hypersensitivity will contribute to the further understanding of its pharmacodynamics, and may lead to the future development of new and effective serotonin-based treatments for human pain conditions.

Materials and methods

Docking Simulations of Serotonin to COMT Active Site

We employ the MedusaDock package to generate possible ligand conformations within the protein active site, and utilize the MedusaScore package to evaluate each conformation generated by MedusaDock 12 13 . MedusaDock enables flexible docking of both the ligand and the side-chain amino acids of the protein. We perform all ligand docking simulations using the crystal structure of human COMT (PDB: 3BWM) 14 .

Prior to all docking simulations, we first stripped all crystallized ligands bound to COMT. Crystallographic waters that were in the active site were retained for our docking simulations. We performed 200 docking simulations for serotonin, where each simulation began with a different seed number and each conformation generated from a simulation was subsequently minimized. All conformations were then ranked according to their free energy value. The lowest energy structure is determined to be the native pose.

Binding Measurements using Surface Plasmon Resonance

Purified COMT protein (derived from porcine liver) and serotonin ligand were obtained from Sigma-Aldrich. Protein concentration was verified by performing a BCA assay (Thermo Scientific). Purified COMT was biotinylated at its surface lysine residues using EZ-Link NHS-LC-LC Biotin (Thermo Scientific) at a 1:1 mole ratio of biotinylating reagent to protein and subsequently purified using a 1 mL Sephadex G-25 medium spin column. Biotinylated COMT was buffered with 10 mM Tris pH 7.4, 1 mM MgCl2, 1 μM DTT and stored at 4°C until use.

Surface plasmon resonance (SPR) measurements were conducted at 25°C using a Biacore 2000 instrument with a previously established protocol 39 . Biotinylated COMT was first loaded onto a streptavidin-coated flow cell (sensor chip SA or Biotin CAPture Kit, GE Healthcare) at a flow rate of 5 μL/min, followed by a buffer flow, introduction of serotonin (buffered in 10 mM Tris pH 7.4, 1 mM MgCl2, 1 μM DTT), and then a final wash with buffer (10 mM Tris pH 7.4, 1 mM MgCl2, 1 μM DTT). Here the flow rate was established at 5 μL/min.

The dissociation rate of the serotonin-COMT complex was determined by fitting the dissociation data to a single-exponential decay.

R = R o × e k d ( t t o )

Here, the response unit from the SPR machine is denoted as R. Once the dissociation rate is determined, the association rate of serotonin to COMT can be found by fitting to the following equation:

R = R eq [ 1 - e - ( k a Serotonin + k a ) ( t - t a ) ]

Radiolabeled SAM Activity Assay

Activity of COMT was measured through a series of experiments where the enzyme was incubated with the substrate catechol (dissolved in H2O, Fisher Scientific) and S-[methyl-3 H]-adenosyl-L-methionine (PerkinElmer, 10.0 Ci/mmol) 40 . Reactions were carried out in a buffer containing 10 mM Tris (pH 7.4), 1 mM MgCl2, and 1 μM DTT with a total reaction volume of 100 μL.

To determine the linear range of the activity assay, SAM concentration (1 μM) was held constant while COMT concentration was varied (10 – 80 ng/μL). Each reaction was placed in a clean PCR tube and incubated at 37°C for 30 minutes. After incubation, reactions were terminated by addition of 100 μL of 1 M HCl to each tube. Radiolabeled catechol products were extracted from the reaction by adding 10 mL of scintillation fluid (MonoFlow I, National Diagnostics) and quantified using a scintillation counter. Activities were normalized for each reaction by performing a duplicate reaction with the COMT inhibitor, OR-486 (0.5 mg of inhibitor dissolved in 30 μL DMSO and 20 μL of 0.9% saline; OR-486 buffer alone does not inhibit reaction), and subtracting its radioactivity. We found that the linear range was between 10 and 40 ng/μL, and thus for all subsequent experiments we used a COMT concentration of 20 ng/μL.

Kinetic experiments were performed at a variety of catechol concentrations (0, 0.5, 1, 3.3, 5, 10, and 33 μM) while maintaining constant SAM (1 μM) and COMT (20 ng/μL) concentrations. Reactions were set-up for a variety of time points: 0, 5, 10, 30, 60, and 120 minutes. Each time-point reaction is incubated at 37°C and then quenched with 1 M HCl at its specified time point. Products were quantified as discussed above. Serotonin-inhibition experiments were conducted for the same catechol concentrations, but with the addition of 100 μM of serotonin (Sigma) to each reaction.

Kinetic parameters were determined by first fitting the data sets to a single-site Michaelis-Menten model (for reference), followed by a two-site Michaelis-Menten model as given by 17 :

v = V max SAM Cat ( SAM Cat + K m , SAM Cat + K m,Cat + K D,SAM K m,Cat )

where catechol is abbreviated as “Cat.” To determine the mechanism of serotonin-inhibition, we calculated the apparent K i (K i,app ) of inhibition for both a competitive and non-competitive inhibitor, and compared which K i,app yielded the lowest value. The apparent K i for a competitive inhibitor is given by

K i,app = K m,Cat Serotonin K m,app K m,Cat

where K m,app is the apparent K m,Cat for the serotonin-inhibited curve. The apparent K i for a noncompetitive inhibitor is given by

K i,app = V max,spp Serotonin V max

where Vmax,app is the apparent Vmax for the serotonin-inhibited curve.

Mouse Behavioral Experiments

Subjects were naïve, young adult (6–14 weeks old) outbred CD-1® (ICR:CrI) mice of both sexes (Charles River, Boucherville, QC). Animals were housed in groups of 2 to 6 with same-sex littermates, with food (Harlan Teklad 8604) and tap water available ad lib. Approximately equal numbers of male and female mice were used; no sex differences were noted so data from both sexes were pooled. All procedures followed the guidelines and regulations of the Canadian Council on Animal Care, and were approved by the McGill Downtown Animal Care and Use Committee.

Mechanical sensitivity was measured using the automated von Frey test. Mice were placed in individual transparent Plexiglas cubicles set atop a perforated metal floor (with 5-mm diameter holes 7 mm apart), and separated from each other by opaque dividers. Prior to testing, the mice were habituated to their surroundings for 1 h. A von Frey fiber with automatically increasing force (Ugo Basile Dynamic Plantar Aesthesiometer) was applied to the mid-plantar hind paw. Three separate withdrawal threshold determinations (on each hind paw) were taken at baseline (pre-injection) and then averaged. One withdrawal threshold determination (on each hind paw) was taken at every post-injection time point.

For the SAM pre-treatment study, after measuring baseline thresholds, mice received i.p. injections of S- (5’-Adenosyl)-L-methionine chloride (80 mg/kg) (Sigma), followed, 15 min later, by 20 μl i.pl. injections of serotonin hydrochloride (10 uM) (Sigma). Mice were tested at 15, 30, 45, 60, 90, 120, 180 and 240 min following the serotonin injection.

For the β-adrenergic antagonist study, mice received 20 μL i.pl. injections of serotonin hydrochloride (5 uM), followed, 15 min later, by i.p. injection of a cocktail of ICI118,551 (5 mg/kg) and SR59230A (50 mg/kg) (Tocris). Mice were tested at 15, 30, 45, 60, 90, 120, 180 and 240 min following the antagonist injection.

Pain hypersensitivity was quantified as the area over the time-threshold curve (post-injection), using the trapezoidal method.

GEE Analysis

We utilize Generalized Estimating Equations (GEE) to determine whether SAM pretreatment or co-administration of SR and ICI treatment versus serotonin alone resulted in significantly different paw withdrawal thresholds 41 across the repeated measurements. The response was the repeated measurements of paw withdrawal threshold at fixed intervals post-administration of treatment. This model included an intercept and a binary covariate reflecting treatment status (e.g., ICI118,551 and SR59230A treatment versus serotonin alone) as the only covariate, where parameter estimates for this covariate can be interpreted as the overall mean change over measurements in PWT due to treatment. We assume an AR 1 correlation structure between repeated measurements, gaussian distribution of the data, and the identity link. We tested whether the estimate for our binary covariate for treatment was significantly different from zero (indicating a significant change in paw withdrawal threshold due to treatment effect) using Proc Genmod in SAS version 9.2, and report p-values from this test. We also report the estimated mean change in paw withdrawal threshold due to treatment effect ± the standard error of the mean change.

Competing interests

The authors declare that they have no competing interests.

Author’s contributions

DT performed the in silico simulations, SPR experiments, COMT activity assays, and wrote the paper. J.S.W. performed the mouse model studies and wrote the paper. N.R. performed the GEE analysis and wrote the paper. R.E.S. performed mouse model studies. R.L.R. performed SPR experiments. S.K.S. developed the COMT activity assay. J.S.M. designed the mouse model experiments and wrote the paper. W.M. conceived of the hypothesis, designed the activity assay experiments, and wrote the paper. N.V.D. designed the in silico experiments, developed the hypothesis, and wrote the paper. L.D. designed the activity assay experiments, developed the hypothesis, and wrote the paper.

Acknowledgments

This work was supported by the US National Institutes of Health grant (R01GM080742 to N.V.D.); American Recovery and Reinvestment Act supplements (GM080742-03 S1, GM066940-06 S1 to N.V.D.x; National Institute of Dental and Craniofacial Research and National Institute of Neurological Disorders and Stroke grants (RO1-DE16558, UO1-DE017018, PO1 NS045685 to L.D. and W.M.).

<p>Genetic basis for individual variations in pain perception and the development of a chronic pain condition</p>DiatchenkoLSladeGDNackleyAGBhalangKSigurdssonABelferIGoldmanDXuKShabalinaSAShaginDHum Mol Genet20051413514315537663<p>Catechol-O-methyltransferase inhibition increases pain sensitivity through activation of both beta2- and beta3-adrenergic receptors</p>NackleyAGTanKSFechoKFloodPDiatchenkoLMaixnerWPain200712819920810.1016/j.pain.2006.09.022190586117084978<p>Comt1 genotype and expression predicts anxiety and nociceptive sensitivity in inbred strains of mice</p>SegallSKNackleyAGDiatchenkoLLariviereWRLuXMarronJSGrabowski-BoaseLWalkerJRSladeGGauthierJGenes Brain Behav2010993394610.1111/j.1601-183X.2010.00633.x297580520659173<p>Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure</p>NackleyAGShabalinaSATchivilevaIESatterfieldKKorchynskyiOMakarovSSMaixnerWDiatchenkoLScience20063141930193310.1126/science.113126217185601<p>The role of serotonin in analgesia elicited by morphine in the periaqueductal gray matter (PAG)</p>SchulRFrenkHBrain Res199155335335710.1016/0006-8993(91)90849-Q1681985<p>5-Hydroxytryptamine and antinociception</p>RobertsMHNeuropharmacology1984231529153610.1016/0028-3908(84)90097-26098854<p>Serotonin is a directly-acting hyperalgesic agent in the rat</p>TaiwoYOLevineJDNeuroscience19924848549010.1016/0306-4522(92)90508-Y1534874<p>Receptor mediation of 5-HT-induced inflammation and nociception in rats</p>SufkaKJSchomburgFMGiordanoJPharmacol Biochem Behav199241535610.1016/0091-3057(92)90058-N1531705<p>Serotonin receptor 5-HT2B mediates serotonin-induced mechanical hyperalgesia</p>LinSYChangWJLinCSHuangCYWangHFSunWHJ Neurosci2011311410141810.1523/JNEUROSCI.4682-10.201121273425<p>5-HT2A receptor subtype is involved in the thermal hyperalgesic mechanism of serotonin in the periphery</p>TokunagaASaikaMSenbaEPain19987634935510.1016/S0304-3959(98)00066-99718253<p>Catechol-O-methyltransferase (COMT): biochemistry, molecular biology, pharmacology, and clinical efficacy of the new selective COMT inhibitors</p>MannistoPTKaakkolaSPharmacol Rev19995159362810581325<p>MedusaScore: an accurate force field-based scoring function for virtual drug screening</p>YinSBiedermannovaLVondrasekJDokholyanNVJ Chem Inf Model2008481656166210.1021/ci8001167266500018672869<p>Rapid flexible docking using a stochastic rotamer library of ligands</p>DingFYinSDokholyanNVJ Chem Inf Model2010501623163210.1021/ci100218t294761820712341<p>Crystal structures of human 108 V and 108 M catechol O-methyltransferase</p>RutherfordKLe TrongIStenkampREParsonWWJ. Mol. Biol.200838012013010.1016/j.jmb.2008.04.04018486144<p>Estimation of intrasynaptic norepinephrine concentrations at vascular neuroeffector junctions in vivo</p>KopinIJZukowska-GrojecZBayorhMAGoldsteinDSNaunyn Schmiedebergs Arch Pharmacol198432529830510.1007/BF005043726728040<p>Role of norepinephrine in the pathophysiology and treatment of mood disorders</p>ResslerKJNemeroffCBBiol Psychiatry1999461219123310.1016/S0006-3223(99)00127-410560027<p>Kinetics of human soluble and membrane-bound catechol O-methyltransferase: a revised mechanism and description of the thermolabile variant of the enzyme</p>LottaTVidgrenJTilgmannCUlmanenIMelenKJulkunenITaskinenJBiochemistry1995344202421010.1021/bi00013a0087703232<p>Efflux of dopamine from the synaptic cleft in the nucleus accumbens of the rat brain</p>GarrisPACiolkowskiELPastorePWightmanRMJ Neurosci199414608460937931564<p>Estimation of intrasynaptic norepinephrine concentrations in humans</p>GoldsteinDSZimlichmanRStullRKeiserHRKopinIJHypertension1986847147510.1161/01.HYP.8.6.4713710555<p>What should we be measuring in behavioral studies of chronic pain in animals?</p>MogilJSCragerSEPain2004112121510.1016/j.pain.2004.09.02815494180<p>Animal models of pain: progress and challenges</p>MogilJSNat Rev Neurosci20091028329419259101<p>Pain sensitivity risk factors for chronic TMD: descriptive data and empirically identified domains from the OPPERA case control study</p>GreenspanJDSladeGDBairEDubnerRFillingimRBOhrbachRKnottCMulkeyFRothwellRMaixnerWJ Pain201112T617410.1016/j.jpain.2011.08.00622074753<p>Enzymatic O-methylation of epinephrine and other catechols</p>AxelrodJTomchickRJ Biol Chem195823370270513575440<p>Chromatographic analysis of the chiral and covalent instability of S-adenosyl-L-methionine</p>HoffmanJLBiochemistry1986254444444910.1021/bi00363a0413530324<p>S-Aristeromycinyl-L-homocysteine, a Potent Inhibitor of S-Adenosylmethionine-Dependent Transmethylations</p>BorchardtRTWuYJ. Med. Chem.19761919719810.1021/jm00223a0431246047<p>Effects of increasing serotonergic receptor activity in brain on analgesic activity in rats</p>LinMTChandraAChiMLKauCLExp Neurol19806854855410.1016/0014-4886(80)90108-96966578<p>The analgesic effects of selective serotonin reuptake inhibitors</p>SmithAJJ Psychopharmacol19981240741310.1177/02698811980120041310065917<p>Serotonin-containing neurons: their possible role in pain and analgesia</p>MessingRBLytleLDPain1977412122060<p>Antidepressants as analgesics: a review of randomized controlled trials</p>LynchMEJ Psychiatry Neurosci2001263036140804011212591<p>Mechanism-based treatment in chronic neuropathic pain: the role of antidepressants</p>ColuzziFMattiaCCurr Pharm Des2005112945296010.2174/138161205486499316178755<p>Antidepressants in the management of cancer pain</p>MaysTACurr Pain Headache Rep2001522723610.1007/s11916-001-0036-811309210<p>A double-blind placebo controlled trial of fluoxetine in fibromyalgia</p>WolfeFCatheyMAHawleyDJScand J Rheumatol19942325525910.3109/030097494091037257973479<p>A randomized, double-blind crossover trial of fluoxetine and amitriptyline in the treatment of fibromyalgia</p>GoldenbergDMayskiyMMosseyCRuthazerRSchmidCArthritis Rheum1996391852185910.1002/art.17803911118912507<p>A randomized, placebo-controlled, double-blind, flexible-dose study of fluoxetine in the treatment of women with fibromyalgia</p>ArnoldLMHessEVHudsonJIWelgeJABernoSEKeckPEAm J Med200211219119710.1016/S0002-9343(01)01089-011893345<p>A randomized controlled trial of citalopram in the treatment of fibromyalgia</p>NorregaardJVolkmannHDanneskiold-SamsoeBPain19956144544910.1016/0304-3959(94)00218-47478688<p>A double-blind, multicenter trial comparing duloxetine with placebo in the treatment of fibromyalgia patients with or without major depressive disorder</p>ArnoldLMLuYCroffordLJWohlreichMDetkeMJIyengarSGoldsteinDJArthritis Rheum2004502974298410.1002/art.2048515457467<p>Efficacy and safety of duloxetine for treatment of fibromyalgia in patients with or without major depressive disorder: Results from a 6-month, randomized, double-blind, placebo-controlled, fixed-dose trial</p>RussellIJMeasePJSmithTRKajdaszDKWohlreichMMDetkeMJWalkerDJChappellASArnoldLMPain200813643244410.1016/j.pain.2008.02.02418395345<p>A review of SSRIs and SNRIs in neuropathic pain</p>LeeYCChenPPExpert Opin Pharmacother2010112813282510.1517/14656566.2010.50719220642317<p>The rate and equilibrium constants for a multistep reaction sequence for the aggregation of superoxide dismutase in amyotrophic lateral sclerosis</p>KhareSDCaplowMDokholyanNVProc Natl Acad Sci U S A2004101150941509910.1073/pnas.040665010152406815475574<p>Functional analysis of genetic variation in catechol-O-methyltransferase (COMT): effects on mRNA, protein, and enzyme activity in postmortem human brain</p>ChenJLipskaBKHalimNMaQDMatsumotoMMelhemSKolachanaBSHydeTMHermanMMApudJAm J Hum Genet20047580782110.1086/425589118211015457404LiangKZegerSLLongitudinal data analysis using generalized linear models Biometrika1986731322