Activity of Bradykinin B2 Receptor Is Regulated by Long-Chain Polyunsaturated Fatty Acids

The molecular and cellular mechanisms by which long-chain polyunsaturated fatty acids (LCPUFA) exert their beneficial effects on cardiovascular health remain obscure. While both LCPUFA and bradykinin (BK) signaling pathway play a role in the cardiovascular system, any direct link between the two is yet to be established. Using picosecond time-resolved fluorescence microscopy and a genetically engineered bradykinin B2 receptor (B2R) sensor (B2K-CC), we detected LCPUFA-induced conformational responses in the B2R similar to those caused by its cognate ligand, BK. The selective B2R antagonist (HOE-140) blocked the eicosapentaenoic acid (EPA, C20∶5, n-3) induced conformational response of the B2K-CC. Further analysis suggests that LCPUFA are capable of direct, B2R-dependent activation of extracellular ligand-regulated kinases (ERK). From a wide range of fatty acids studied, varying in chain length, saturation, and position of double bonds, EPA, docosahexaenoic (DHA, C22∶6, n-3), docosadienoic (DDA, C22∶2, n-6), and dihomo-gamma linoleic (DGLA, C20∶3, n-6) fatty acids caused the highest ERK phosphorylation. EPA or DHA dependent ERK phosphorylation was inhibited by the selective B2R antagonist. We show that LCPUFA stimulates downstream signaling by B2R such as B2R-dependent phosphorylation and expression regulation of endothelial nitric-oxide synthase (eNOS). Further analysis indicated that LCPUFA also alters levels of the eNOS transcription factor, kruppel-like factor 2 (KLF2). Moreover we show that EPA increases membrane fluidity on the same time scale as B2R conformational response, suggesting that partitioning of LCPUFA into bilayer is a primary step required for receptor activation. In summary our data show that LCPUFA activate B2R receptor at nanomolar concentrations suggesting a novel molecular mechanism by which fatty acids may affect the cardiovascular system.


Introduction
We have recently reported that LCPUFAs stimulate conformational changes and activate parathyroid hormone receptor (PTH1R) in a dose dependent manner [1]. PTH1R is a G protein coupled receptor (GPCR) involved in bone homeostasis [2]. There has been increasing consensus that many GPCRs are ''allosteric machines'' controlled by various membrane components such as lipids, ions, cholesterol or mechanical perturbation of the cell membrane [3][4][5][6][7][8][9][10]. For example, we have shown that bradykinin B 2 R (a GPCR expressed in endothelial cells [11]) responds to mechanical perturbation and changes in fluidity of the cell membrane [12]. B 2 R activates multiple signaling cascades but generally is reported to signal through G q -protein -phospholipase C pathway stimulating phosphoinositide hydrolysis and transient increase in intracellular free Ca 2+ . B 2 R also signals through the G s -cAMP-protein kinase A and mitogen-activated protein kinases (MAPKs) leading to various biological effects [11]. B 2 R signaling pathway is involved in activation of eNOS, the enzyme that catalyze the production of nitric oxide (NO) [13], which is essential in maintaining normal vessel homeostasis [14]. The coronary dilation induced by bradykinin is mediated by the activation of the B 2 R. A number of earlier studies have demonstrated that omega-3 polyunsaturated fatty acids induce in-vivo vasodilation [15][16][17][18][19][20]. These observations along with established effects of fatty acids on cardiovascular health provide a strong rationale for the study of the signal transduction pathway of LCPUFAs in endothelial cells. Since endothelial cells are exposed to fatty acids at high concentrations, in this study we sought to determine if LCPUFA can also regulate function of the B 2 R. The venous blood concentration of fatty acids is known to vary from 250 mM to 3 mM depending on the nutritional state [21]; most of these fatty acids are bound to serum albumin while a small percent is unbound in the plasma. The amounts incorporated into plasma membranes (PM) of cells can be enhanced by up to 10 times through dietary supplements enriched in LCPUFA [22][23][24]. Further examination of the fatty acid composition in blood reported that roughly 5.3% of total serum fatty acids are LCPUFA (2.4% docosahexaenoic (DHA, C22:6, n-3), 1.8% eicosapentaenoic (EPA, C20:5, n-3), 1% docosapentaenoic (DPA, C22:5, n-3 ) [25]; similar fatty acid concentration/composition has been reported in mice [26,27].
It is speculated that due to its intrinsic ability to affect so many seemingly unrelated conditions, LCPUFA must function at a fundamental level common to most cells. There is strong evidence suggesting LCPUFA act by altering cell plasma membrane properties [28,29] that affect GPCR function, however a direct interaction between LCPUFA and B 2 R has to be considered.
To gain insights into the mechanisms by which LCPUFAs influence endothelial cells, we have chosen to study B 2 R-dependent activation of the MAPK/ERK pathway, which is known to have an important role in the cardiovascular system [30]. In this study we tested if selected LCPUFAs can: (1) cause conformational changes in B 2 R, (2) activate ERKs via B 2 R, (3) modulate expression and activate eNOS; and (4) induce changes in expression of KLF2. Additionally we aimed to determine if a selective B 2 R antagonist can block EPA-mediated activation of B 2 R conformational changes and signaling. Furthermore, we investigated the dynamics of LCPUFA incorporation into cell membrane and changes in plasma membrane fluidity induced by LCPUFA.

Detection of B2K-CC Conformational Change with LCPUFA
To test if LCPUFAs induce conformational changes in B 2 R, we used the previously described B 2 R FRET sensor (B2K-CC) containing an intramolecular FRET pair that enables detection of conformational activity upon ligand stimulation [12]; a number of other GPCR FRET sensors has been constructed and used to study conformational GPCR activity using similar approach [9,12,[31][32][33][34][35][36]. B2K-CC was expressed in human embryonic kidney 293 (HEK293) cells. In the presence of EPA, a decrease in the FRET ratio was observed with time constants (t) of 4.960.1 min for 200 mM and 5.260.3 min for 10 mM treatments, respectively. In contrast Oleic acid had no effect on the B2K-CC FRET ratio (Fig. 1A). HEK293 cells expressing the control FRET sensor (PM-CC) showed no response to EPA (Fig. 1A). To test specificity of EPA-induced B 2 R conformational change, preincubation of the B2K-CC transfected HEK293 cells with the selective B 2 R antagonist, HOE-140 [37], blocked the conformational response of the B 2 R to EPA. Additionally, treatment with 100 mM of DDA (C22:2, n-6) and DGLA (C20:3, n-6) also caused a decrease in the FRET ratio observed with time constants of t = 2.4960.3 min and 5.4160.4 min, respectively (Fig. 1B). Although GPCR FRET constructs of this type have been shown to have a reduced efficacy of cognate ligands in stimulating conformational changes (due to the insertion of fluorescing proteins into the native receptor structure [9,38]), we were able to observe significant FRET response to LCPUFA.
Fatty Acid Induced ERK Phosphorylation in HEK293 Cells Transfected with B 2 R To examine fatty acid activation of the ERK1/2 cascade via B 2 R, we treated B 2 R transfected HEK293 cells with distinct fatty acids for 5 minutes as shown on Fig. 2A. Untransfected cells were also treated with fatty acids since HEK293 express low levels of endogenous B 2 R [39,40]; no significant increase in ERK1/2 phosphorylation was found compared to non-treated transfected cells, indicating that ERK1/2 activation was B 2 R dependent.
Note, that simultaneous treatment with BK (100 nM) and fatty acid (10 mM EPA) did not produce any additional increase in ERK phosphorylation beyond that caused by BK alone.
To exclude the possibility that LCPUFA metabolites are activating the ERK pathway, we treated B 2 R transfected HEK293 cells with EPA or DHA for 15 seconds (Fig. 2B). The short treatment significantly increased ERK1/2 phosphorylation but to a lower magnitude than after the 5 minute treatment; lower activation levels were also observed for the BK treated cells suggesting that the reduction in ERK activation after 15 second treatment is due to the intrinsic time scale needed for activation of ERK and not due to the finite production rate of LCPUFA metabolites. Furthermore, the majority of fatty acids have been reported to remain in the plasma membrane after 10 minutes [41] while the LCPUFA metabolizing enzyme, cyclooxygenase is localized in the endoplasmic reticulum [42] making it unlikely that the LCPUFAs metabolites are activating B 2 R in our experiments (in cells the prostaglandin production rate from EPA is of 4 to 5 hours [43], whereas in vitro the cyclooxygenase reaction kinetics with EPA at nanomolar concentrations is expected to be 7 minutes [44]).

B 2 R Antagonist Blocks EPA Induced Activation of the ERK Cascade
To determine if a B 2 R selective antagonist can inhibit fatty acidinduced ERK phosphorylation, HEK293 cells transfected with B 2 R were treated with EPA in the presence of B 2 R inhibitor, HOE-140 (Fig. 3A). Phospho-ERK immunoblot analysis showed that the antagonist is able to significantly inhibit ERK phosphorylation of cells treated with bradykinin (positive control) and EPA. To determine if these results were reproducible in a different cell line, we used bovine aortic endothelial cells (BAEC) expressing B 2 R at endogenous levels (no transfection was used). The cells were treated with BK (positive control), EPA, or DHA in the presence of the inhibitor. Phospho-ERK immunoblot analysis demonstrated that the inhibitor significantly blocked the ERK phosphorylation of the treated cells (Fig. 3B). To determine inhibitor specificity, we treated BAECs with noradrenaline (NA) to activate the endogenous adrenergic receptor and determine if the HOE-140 inhibitor would unspecifically inhibit ERK phosphorylation. The HOE-140 did not have any effect in inhibiting adrenergic receptor stimulation by NA (Fig. 3C) suggesting that HOE-140 selectively antagonizes B 2 receptor.

Measurement of Antagonist Potency of HOE-140 on B 2 R Expressed in HEK293 Cells
To further examine the mechanism of action of antagonist HOE-140 on inhibiting B 2 R activation by EPA, we have performed a Schild analysis of the antagonist-induced inhibition of ERK activation. HOE-140 produced a parallel rightward shift in the dependence of the B 2 R-mediated ERK activation on EPA concentration (Fig. 4A). The antagonist did not significantly affect the maximum possible effect (E max ) of EPA and did not detectably affect ERK phosphorylation in the absence of the agonist (Fig. 3). The Schild slope was 1.1360.08 (Fig. 4B). These observations suggest that HOE-140 acts as a competitive antagonist of EPA stimulated ERK phosphorylation at the B 2 R over the range of antagonist concentrations tested.

B 2 R Antagonist Blocks EPA-induced Activation and Expression of eNOS
To investigate if EPA-dependent B 2 R activation has a direct biological effect, we investigated eNOS activation via phosphorylation at serine 1177 which is known to activate the eNOS enzyme leading to production of nitric oxide (NO) [45]. Treatment of BAECs with BK or EPA for 5 minutes induced an increase in eNOS phosphorylation (Fig. 5A). To determine if eNOS activation is a result of EPA-induced B 2 R activation, we repeated the experiment in the presence of the antagonist HOE-140. The antagonist blocked the B 2 R activation by BK and EPA.

NOS3 Expression is Modulated by EPA-dependent B 2 R Activation
To further characterize the association between the fatty acids, B 2 R and eNOS activation, alterations in NOS3 gene expression were studied. BAECs were treated with BK, EPA, or both for 1 hour prior to gene expression analysis via RT-PCR. NOS3 expression levels were increased as a result of BK or EPA stimulation while simultaneous treatment with BK and EPA did not increase the expression further (Fig. 5B). Specificity to B 2 R stimulation was determined by the ability of the B 2 R antagonist, HOE-140, to significantly inhibit the increase in expression. To support these findings, the experiment was repeated in human coronary artery endothelial cells (HCAEC) and similar results were obtained (Fig. 5C).

Expression of NOS3 Transcription Factor, KLF2, is Modulated by EPA-dependent B 2 R Activation
To gain insights into the possible upstream regulation of EPAinduced B 2 R stimulation that causes alterations in NOS3 expression, we investigated expression of Kruppel-like factor 2 (KLF2; which is a known transcription factor that induces eNOS expression [46]). Treatment of BAECs with BK or EPA alone, or in combination for 1 hour was performed prior to analysis of gene expression via RT-PCR. KLF2 expression levels were increased as a result of BK or EPA stimulation while simultaneous treatment with both did not increase the expression further (Fig. 5D). B 2 R dependent regulation of KLF2 expression was confirmed by the ability of the selective B 2 R antagonist to significantly inhibit the expression increase.

EPA Increases Plasma Membrane Fluidity of HEK293 Cells
To gain insights into possible non-specific mechanisms by which LCPUFA is activating the B 2 R and to determine how fast EPA partitions into the cell membrane, we investigated effects of fatty acids on membrane viscosity. Using the molecular rotor farnesyl-(2-carboxy-2-cyanovinyl)-julolidine (FCVJ) and picosecond timeresolved fluorescent microscopy, we recorded time-dependent changes in membrane fluidity in response to EPA treatment. Since FCVJ has a viscosity-dependent quantum yield [47], we stained HEK293 cells with FCVJ and recorded its signal intensity. Treatment with 100 mM of EPA caused a time-dependent decrease in fluorescence intensity compared to the control solution, indicating an increase in membrane fluidity with a time scale of 5.460.3 min (Fig. 6). A treatment with control solution of 0.1% Ethanol and 3.5610 25% BSA in DPBS did not induce significant changes in fluidity.

Discussion
Analysis of the FRET data ( Fig. 1) indicates that EPA-treated B 2 R FRET sensor undergoes the conformational change (negative change in FRET ratio) that is similar to the one induced by its cognate ligand, BK [12] but occurs on a longer time scale. We have previously shown that the time scale of B2K-CC FRET sensor response to BK stimulation is a few seconds (depending on the concentration) [12]. In contrast Fig. 1 shows that LCPUFA induce conformational changes in B 2 R on a time scale of 2 to 5 minutes. Note, that the time scale of the B2K-CC response to LCPUFA is similar to the time scale of LCPUFA incorporation into the membrane bilayer (see below). The kinetics of FRET response (Fig. 2B) is expected to be different between DDA and DGLA due to differences in the structure of these fatty acids. Compared to the untreated receptor, our data show that EPA, DDA, and DGLA induce a decrease in FRET signal that ultimately saturates and stabilizes at a different FRET level, indicating that B 2 R changes conformational state (Fig. 1).
The maximal effect (E max ) of fatty acids on ERK1/2 phosphorylation is lower as compared to stimulation by BK (Fig. 1A, Fig. 2A,B) suggesting that the interactions with fatty acids stabilize the receptor in an activated state that is different from the one stimulated by BK. Although our primary focus in this study was EPA, other LCPUFAs, both n-3 and n-6, also caused a significant, B 2 R-dependent ERK phosphorylation with varying potency (Fig. 2A). To decipher a possible trend in fatty acid induced B 2 R-dependent ERK activation, we tested fatty acids with different alkyl chain lengths, saturation level and positions of double bonds. Our data presented in Fig. 2A suggest that the length of the fatty acid and the number of carbon chains play a major role in B 2 R activation. More specifically, our data suggest that significant B 2 R activation requires fatty acids to have an alkyl chain length of no less than 20 carbons with 3 or more double bonds, except for alpha-linolenic acid (18:3 n-3) which showed significant potency in ERK activation. Note, EPA appears more potent in the FRET assay than in the ERK assay due to differences in treatment concentrations. Overall, our results are similar to those we earlier observed for the PTH1 receptor [1] and consistent with other studies in which specific GPCRs have been shown to become activated by fatty acids in a length-dependent manner [48][49][50][51].
The slope of the Schild plot [52,53] (Fig. 4B) suggests that EPA is interacting with the orthosteric binding site of B 2 R. Note, however, that the concentration of fatty acids in the membrane are expected to be significantly higher than concentrations in the suspension buffer due to natural tendency of LCPUFA to partition into bilayers [28]. The lipid bilayer effectively increases the concentration of LCPUFA in GPCR environment thereby  increasing the probability of their interaction. In contrast to our earlier report on PTH1 receptor [1] we did not observe synergism from simultaneous stimulation by BK and LCPUFA (Fig. 2 A, 5 B, 5 C) supporting the hypothesis that EPA acts through the orthosteric site of B 2 R. However due to the high concentrations of LCPUFA in the membrane, contributions from allosteric and nonspecific interactions cannot be excluded.
BK is a potent stimulator of endothelial nitric-oxide synthase (eNOS), the enzyme that catalyze the production of nitric oxide (NO) via B 2 R signaling pathway [13]. NO contributes to vessel homeostasis by inhibiting vascular smooth muscle contraction and growth, platelet aggregation, and leukocyte adhesion to the endothelium [14]. In our study, EPA had similar biological effects on eNOS as BK. EPA was able to increase phosphorylation of serine 1177 which activates the eNOS enzyme to produce nitric oxide (NO) [45] in a B 2 R-dependent manner. According to a recent study, EPA stimulates phosphorylation of AMP-activated protein kinase (AMPK)-Thr172 leading to an increase in eNOS phosphorylation and NO release [54], while our present results similarly show that EPA increases B 2 R-dependent expression and activation of eNOS. NO plays an important role in the protection against the onset and progression of cardiovascular disease [14] and this finding provides new insights into the putative mechanism by which LCPUFAs act at preventing and alleviating CVDs; these findings are in agreement with previous studies that reported EPAinduced NO production in endothelial cells [19,20].
KLF2 is an endothelial transcription factor that is uniquely induced by fluid shear stress [55]. Blood flow turbulence at vessel bifurcations and curvatures have been implicated as a factor in atherosclerosis [56]. The ability of LCPUFA to increase expression levels of KLF2 (Fig. 5 D), which was the first endothelial transcription factor shown to be uniquely induced by flow, might provide a molecular basis for beneficial effects of LCPUFA; notably it has been recently reported that LCPUFA supplemented diet plays a similar role in preventing bone loss as mechanical loading of bone (which is known to be fluid shear stress dependent) [57].
Changes in membrane-dependent functions of receptors and membrane associated proteins have been hypothesized to depend on cell membrane fluidity [9,10,12,47,58]. It has been shown that stimulation with BK decreases cell membrane fluidity [59] which is in contrast to the fluidizing effect of EPA on the cell membrane; therefore LCPUFA-induced changes in membrane fluidity are unlikely to be responsible for the observed response of B 2 R to LCPUFA. The time scale of LCPUFA incorporation into membrane (Fig. 6) is within the same range as that one of EPAinduced conformational changes in both B 2 R and PTH1R [1]. This finding suggest that fatty acids first partition into the cell membrane and then may either specifically bind to the transmembrane/intracellular domain of the GPCR or affect membrane properties in a non-specific way, stabilizing GPCR in a tertiary conformation required for efficient signaling. As an example of the effects of other membrane components, several GPCR structure determination efforts have shown that the addition of cholesterol analogs is often critical for maintaining protein stability [60]. The PTH1R and B 2 R show a similar change (decrease) in FRET in response to LCPUFA that may be suggestive of a general underlying mechanism. However, the existence of putative general mechanism by which GPCRs are activated by LCPUFA does not diminish the relevance to hemodynamic function or function in other biological tissues because cells in specialized tissues express specific GPCRs at relatively high levels as e.g. B 2 R in endothelial cells [61]; as noted above activation of B 2 R in endothelial cells is known to invoke a number of physiological hemodynamic responses.
While bradykinin treatment favors cardioprotective properties, continuous treatment induces side effects [11]. A similar profile has not yet been reported for EPA or other LCPUFAs. In contrast to bradykinin, one possible explanation for the differences in activity is the ability of EPA to alter other GPCRs that serve other functions. For example, Grp120, a GPCR that binds to EPA, mediates potent anti-inflammatory and insulin sensitizing effects [62]. EPA may also bind to other GPCRs important in the prevention of CVDs.
Currently, therapeutic medications are available for the prevention and/or treatment of CVDs but are expensive and have significant side effects. Studies have shown that supplementation with n-3 LCPUFA after myocardial infarction reduced the risk of recurrence, stroke, and death by 10% [63]. Many studies suggest n-3 LCPUFAs, derived from seafood and plant sources may reduce CVD risk [64,65]. n-3 fatty acids have been reported to contain antihypertensive effects, reducing features of inflammatory atherosclerotic plaques, and decrease in death rate over a two year post myocardial infarction [66]. The importance of BK signaling in functioning of cardiovascular system and the effects of LCPUFA on BK pathway described in this study call for more detailed studies of the role of LCPUFA on BK pathway. Further studies in both in in vitro and in vivo models may contribute to development of alternative approaches for CVD treatments.
Overall, we characterized the effect of LCPUFA on the B 2 R and our results indicate that EPA is a novel B 2 R agonist. More specifically, we demonstrated the ability of LCPUFA (EPA, DHA, DDA, or DGLA) to act as agonist to the B 2 R at nanomolar concentrations and to induce receptor conformational changes. EPA was able to induce a B 2 R response with an EC50 value of 214 nM 6107 nM. The Schild analysis indicated that EPA was competitively binding to the orthosteric site of the B 2 R. Furthermore, we showed that other LCPUFAs were also able to activate B 2 R at the concentrations of fatty acids typically found in blood. Experiments to determine activation of downstream signaling pathways revealed an increase in ERK phosphorylation, regulation of eNOS activity and expression. Collectively, these results provide a putative molecular basis for the involvement of the B 2 R in mediation of LCPUFAs effects in endothelial cells. We believe our study is a starting point for more detailed investigations of interactions between fatty acids and other types of GPCRs.

Materials and Methods
Cell Culture, Transfection, and Chemicals

Fatty Acid and Inhibitor Treatment
To examine fatty acid dependent ERK phosphorylation, HEK293 cells transfected with B 2 R were treated with various fatty acids at 10 mM concentration containing 0.1% Ethanol and 3.5610 25 % BSA in DPBS at 37uC for 5 minutes. The negative control treatment consisted of only 0.1% Ethanol and 3.5610 25 % BSA in DPBS at 37uC for 5 minutes. To characterize the effect of B 2 R inhibitor on fatty acid stimulation, HEK293 cells transfected with B 2 R were incubated with either 1 mM, 300 nm, 100 nm, or 10 nm of the B 2 R inhibitor, HOE-140, for 30 minutes prior to fatty acid treatment. A control FRET sensor that is localized to the plasma membrane of cells (PM-CC) was constructed as previously described [67]. The construct PM-CC encodes a protein that is comprised of Citrine and Cerulean linked together with a short and flexible GGGGPV (ProVal to encode Age1 restriction site) linker peptide to ensure Förster resonance energy transfer between the two fluorescent proteins, and is fused to another 10 residue leader peptide MGCINSKRKD to direct its translocation to plasma membrane [68].
B2K-CC was constructed and characterized as described in our previous study [69]. A primer extension method was used to construct B2K-CC by inserting EYFP into the third cytoplasmic loop at a position between F232 and K233 and fusing ECFP to the truncated C terminus at the position C329. The sequences of all constructs were confirmed by sequencing service.

FRET Measurements
FRET measurements in single living cells, 24 hours after being transfected, were performed using multichannel time resolved single photon counting as previously described [69]. Briefly, fluorescence emission kinetics and spectra were measured by using a multichannel, time-correlated single photon counting spectrograph (PML-16/SPC630; Becker & Hickl, Berlin, Germany) coupled to an inverted microscope (Axiovert 200 M; Zeiss, Thornwood, NY) via fiber optic link. A femtosecond Ti:Sapphire oscillator (Spectra-Physics, Irvine, CA, USA) was used as the excitation source. The repetition frequency of the light pulses from the oscillator was reduced to 8 MHz and the wavelength was doubled to 435 nm. The excitation light was defocused to a spot size of 20-50 mm to enable spatially homogenous excitation of a single cell. Single cell fluorescence spectra were obtained by integrating time-resolved fluorescence data. Presented data were recorded by detecting fluorescence emission polarized at the magic angle (54.7u) to the polarization of the excitation light at 435 nm to reduce contributions from potential reorientation of the sensor in the plasma membrane under shear stress. FRET ratio was defined as ratio of Citrine emission intensity at , 525 nm to Cerulean emission intensity at 475 nm. Duration FRET measurements were kept to a minimum to prevent photobleaching.

FCVJ Staining
A staining solution was prepared from 50 ml stock solution of 20 mM of the FCVJ in dimethyl sufoxide (DMSO), which was dispersed in 2 ml fetal cow serum under vigorous stirring as described [47]. 10 ml of DMEM was added. The slides with the cells were covered with the staining solution. Incubation took place over 10 min in the dark at 37uC. After the incubation period, the cells were rinsed with HBSS, and the fluorescence checked under an epifluorescent microscope (Diaphot TMD, Nikon, Garden City, NY, USA) using the G2B filter set.

Statistical Analyses
Statistical significance was evaluated for differences between groups from at least three independent experiments using Stu dent's t-test. A p value of ,0.05 was considered to be statistically significant. Calculations of the mean, standard error of the mean, and Student's t-test were performed using Microsoft Excel and OriginPro (OriginLab Co., Northampton, MA).