Targeted Metabolomic Profiling Reveals Association Between Altered Amino Acids and Poor Functional Recovery After Stroke

Amino acids have been shown to be among the most important metabolites to be altered following stroke; however, they are a double-edged sword with regard to regulating hemostasis. In this study, we conducted a targeted metabolomic study to examine the association between serum levels of amino acids and functional recovery after stroke. Three hundred and fifty-one patients with stroke admitted to an acute rehabilitation hospital were screened, and 106 patients were selected based on inclusion and exclusion criteria. Recruited patients were stratified using Montebello Rehabilitation Factor Score (MRFS) efficiency. We selected the top (n = 20, 19%) and bottom (n = 20, 19%) of MRFS efficiency for metabolomic analysis. A total of 21 serum amino acids levels were measured using ultra high performance liquid chromatography and mass spectrometry. The normalized data were analyzed by multivariate approaches, and the selected potential biomarkers were combined in different combinations for prediction of stroke functional recovery. The results demonstrated that there were significant differences in leucine-isoleucine, proline, threonine, glutamic acid, and arginine levels between good and poor recovery groups. In the training (0.952) and test (0.835) sets, metabolite biomarker panels composed of proline, glutamic acid, and arginine had the highest sensitivity and specificity in distinguishing good recovery from poor. In particular, arginine was present in the top 10 combinations of the average area under the receiver operating characteristic curve (AUC) test set. Our findings suggest that amino acids related to energy metabolism and excitotoxicity may play an important role in functional recovery after stroke. Therefore, the level of serum arginine has predictive value for the recovery rate after stroke.

Introduction

Increasing evidence suggests that amino acids, including homocysteine and branched-chain amino acids (BCAA), are one of the most important disturbed metabolites after stroke (1, 2). Importantly, recent studies suggest that amino acids could have beneficial and detrimental effects. For example, glutamate plays an important role in maintaining the normal signal transduction of nerve cells, which is beneficial to the synaptic plasticity of neurons and to the recovery of stroke (3). However, elevated levels of glutamate can trigger oxidative stress, inflammation, and endothelial damage (4). BCAAs were also found to be low in stroke patients compared with normal controls, and lower BCAA levels correlated with poor neurological outcome in stroke patients (5). On the other hand, higher concentrations of baseline BCAA were associated with increased risk of stroke in a high cardiovascular risk population (6). Therefore, elucidating the effects of dysregulated amino acid levels on stroke recovery will contribute to identifying prognostic biomarkers and formulating effective therapeutic interventions. However, there are few studies on the relationship between the changes in amino acid levels and the rate of functional recovery after stroke (7, 8).

Metabolomics is a new strategy that can detect changes in amino acids, vitamins, organic acids, and other small-molecule metabolites in biofluids (e.g., plasma or serum) of patients in real time (9, 10). As it is difficult to detect the metabolites directly in the brain, metabolites in serum are usually used as alternative indicators to reflect biological and pathological functions of the brain (11, 12). Metabolic alterations in the brain can result in changes in the metabolome of biofluids (12, 13), especially those metabolites with low molecular weight, which may be easily exchanged through the meningo between the cerebrospinal fluid (CSF) and the blood (5, 14). As a result, potential biomarkers associated with stroke recovery can be detected in serum by metabolomics, and signaling pathways involved in stroke recovery can be drawn.

Currently, most studies of metabolomic biomarkers in patients with stroke compare metabolomic profiles of stroke patients with the normal population (9, 13). Few studies explore the differences in metabolomic biomarkers among stroke patients with good or poor functional recovery (4). Therefore, in this study, we performed a metabolomic analysis of serum amino acid levels in stroke patients undergoing acute inpatient rehabilitation. Our initial hypothesis was that alterations in amino acids can affect the rate of functional recovery in stroke patients. The aim of this study was to find new biomarkers with high sensitivity and specificity, and to provide a basis for further investigation into the rehabilitation mechanism of stroke, as well as the development of targeted treatment methods and improvement of prognosis.

Materials and Methods

Study Participants

This is a retrospective study approved by the ethics committee of Spaulding Rehabilitation Hospital. Three hundred and fifty-one patients with stroke admitted to an acute rehabilitation hospital in Boston between January 2015 and December 2016 were screened with the following inclusion criteria: first ischemic stroke, confirmed by CT or MRI, no tPA (tissue plasminogen activator) treatment, age between 50 and 85 years, length of stay >6 days, admission total Functional Independence Measure (FIM) score between 36 and 71 (12, 15, 16), and has serum samples available upon admission. Exclusion criteria included: continuing with gastric tube feeding, active cancer, HIV (human acquired immunodeficiency) carrier, or severe liver and kidney dysfunction (12, 15, 16). All patients received standard inpatient rehabilitation including physical therapy, occupational therapy, and speech therapy. One hundred and six patients met the inclusion/exclusion criteria and were stratified at the top 20 of MRFS efficiency (top 19%) and the bottom 20 of MRFS efficiency (bottom 19%), and they were defined as good recovery group (GR) and poor recovery group (PR) (Figure 1). The MRFS efficiency formula is described in the following section.

FIGURE 1

Figure 1. Patients with stroke were screened according to the inclusion criteria and the exclusion criteria.

Functional Independence Measure Scale

FIM scale is widely used to evaluate the functional abilities of stroke patients undergoing rehabilitation (15, 17). The scale includes 18 items graded on a 7-point ordinal scale, with a maximum total score of 126 where a lower score means less functional independence. The total FIM scores were recorded at admission and discharge. The gain of the total FIM score, which is commonly used to evaluate functional recovery after stroke, was determined for every patient by calculating the difference of the total FIM score from admission to discharge (15, 17).

Because length of stay and total FIM scores at admission varied considerably, Montebello Rehabilitation Factor Score (MRFS) efficiency was used to evaluate functional recovery after stroke (17, 18). The MRFS evaluates relative gain, and this method depends on the validated FIM score. According to this method, the basis for calculating relative gain is a patient's specific potential for improvement (maximal possible FIM – actual admission FIM). The actual score ranges from 0 to 1, and the MRFS can overcome the misinterpretation of the ceiling effect. The MRFS was calculated using the following formula: MRFS = (discharge FIM – admission FIM)/(maximum FIM score – admission FIM). MRFS efficiency reflects recovery and functional outcomes of stroke more appropriately and precisely because it is measured relative to the potential for change and overcomes the fact that different patients have different admission FIM scores (15, 18). MRFS efficiency = MRFS/length of stay (15, 17, 18).

Clinical Characteristics

Demographic data and clinical characteristics including age, gender, education, body mass index (BMI), diagnosis, comobilities, and discharge destination were collected from electronic medical records. Stroke features including site (supratentorial and infratentorial) and side were collected from MRI or CT reports. Lesion size was not measured due to limited imaging available for the measurement.

Sample Preparation

Non-fasting venous blood was obtained from all stroke patients within 1 week after hospitalization. Studies showed that there were no significant difference of amino acids between fasting serum samples and non-fasting serum samples (19, 20). The blood sample testing process follows standard protocol of ultra high performance liquid chromatography—mass spectrometry (UHPLC-MS), and the standard protocol are as follow.

Blood samples were centrifuged at 13,000 g for 30 min at 4°C, after which time serum was removed and aliquotted before storage at −80°C (not more than 2 years) until ultra high performance liquid chromatography—mass spectrometry (UHPLC-MS) analysis. Prior to analysis, serum samples were thawed and centrifuged at 14,000 g for 10 min in a cold room (4°C). Then 200 μl supernatant was transferred into a new 1.5 ml microcentrifuge tube, and 800 μl of cool methanol (−80°C) (Fisher Scientific, cat. no. A452SK1) was added to make a final 80% (vol/vol) methanol solution. This mixture was incubated for 8 h at −80°C and then centrifuged at 14,000 g for 10 min (4°C). Subsequently, the supernatant was transferred into a new 1.5 ml microcentrifuge tube, dried in a SpeedVac (Savant AS160, Farmingdale, NY), and stored at −80°C until analysis (21, 22).

Each sample was resuspended in 20 μl of UHPLC-MS grade water (Fisher Scientific, cat. no. MWX00016), and 10 μl per sample was analyzed with UHPLC-MS using the selected reaction monitoring (SRM) method with positive/negative ion polarity switching on the hybrid triple quadrupole/linear ion trap mass spectrometer (AB/SCIEX). A total of 21 amino acids were monitored and detected for each sample. Peak areas from the total ion current for each metabolite SRM transition were integrated using MultiQuant v2.0 software (AB/SCIEX) (21, 22).

Statistical Analysis

Statistical analysis was carried out in SPSS 14.0. The T-test was used to analyze the difference in the continuous variables between the two groups, and the chi-square test was used to evaluate clinical categorical measures (23, 24). The Pearson correlation coefficient was used to calculate the correlations between age and identified metabolites in the two groups. The differences between the two groups were considered significant at p < 0.05.

Metabolomics Data Analysis

As glycine could not be detected in most samples, only 20 amino acids' data were analyzed. All peak areas were aligned and normalized using the median of all amino acids from each sample before further analysis. The normalized data were imported into SIMCA-P version 14.1 (Umetrics Inc., Umea, Sweden) for multivariate analysis, including principal component analysis (PCA) and partial least-squares discriminant analysis (PLS-DA) after mean-centering and unit variance (UV) scaling. The potential biomarkers were filtered and confirmed when their Variable Importance of Projection (VIP) scores are >1 (VIP > 1) (5, 12).

The quality of the PLS-DA model was determined based on a goodness of fit parameter (R2Y) and a goodness of prediction parameter (Q2Y). In addition, the PLS-DA model and the reliability models were further validated using a rigorous permutation test (n = 200). The parameters of the models, such as the R2 and Q2 intercepts, were investigated to ensure the quality of the models and to avoid over-fitting (5, 11, 25).

Through analysis of PLS-DA loadings, the metabolites contributing to sample discrimination were identified by Variable Importance of Projection (VIP) scores. The potential biomarkers were filtered and confirmed when their VIP scores are >1 (VIP > 1) (5, 11, 25).

Matlab R2014a (The MathWorks Inc., Natick, MA, USA) was used to perform variable selection of potential biomarkers. The selected potential biomarkers were combined in different combinations for prediction of stroke functional recovery. All prediction combinations were examined separately, using 10-fold cross-validation (26, 27).

In 10-fold cross-validation, nine-tenths of the serum samples from all 40 samples were randomly assigned to the training set. The metabolite profile of this training set was used to diagnose for this prediction task. The remaining one-tenth of the serum samples from all 40 samples formed the test set. This test set was used to validate the metabolite profile diagnostic for the feature of interest. This was repeated 10 times, so that each one-tenth split of the data set acts as the testing set once. Areas under the curve (AUCs) with 95% confidence intervals (CIs) were calculated for sensitivity and specificity values. Mean AUC in the training set, mean AUC in the test set, standard deviation (SD) of AUC in the test set, and 95% CI of mean AUC in the test set were analyzed (26, 27).

Predictive performance results for each prediction combination were compared using area under the receiver operating characteristic curve (AUC of ROC). We note that this predictive performance is for the stratified data set, so our primary interest is in the predictive factors. The values of mean AUC in the test set are regarded as the criteria for selecting the best combination of predictive biomarkers (26, 27). Since the mean test-set AUC scores were used for ranking models and not for formal hypothesis testing, multiple testing corrections were not needed.

Results

Clinical Characteristics of Stroke Patients

The clinical characteristics of stroke patients are shown in the Table 1. According to the effectiveness of MRFS, stroke patients were divided into two subgroups: good recovery (GR) group and poor recovery (PR) group. The average age of GR group was significantly lower than that of PR group (61.25 ± 7.84 vs. 71.55 ± 10.39 years) (p < 0.001). The time of hospitalization, destination of discharge, total score of FIM at admission and discharge, and effective rate of MRFS and MRFS in the GR group were significantly better than those in the PR group (p < 0.001). There was no significant difference in other clinical indexes (including education, BMI, medical history, laboratory items, stroke features, and side of hemiparesis) between the two groups (P > 0.05) (Detailed P-values can be seen in the Table 1).

TABLE 1

Table 1. The clinical characteristics of stroke patients.

Serum Metabolic Profile of Stroke Patients With Good Recovery and Poor Recovery

A PLS-DA model was performed to explore the correlation between the GR and PR groups. According to a PLS-DA score plot (Figure 2A), there was a significant separation between stroke patients with good recovery and poor recovery (R2Y = 0.495, Q2Y = 0.345) (Figure 2A), indicating that there is a difference in serum metabolite levels between the GR and PR groups. The PLS-DA model validation was performed using the number of permutations equaling 200 generated and the intercepts of Q2 (fewer than 0), which meant that the PLS-DA model was non-overfitting and reliable (Figure 2B).

FIGURE 2

Figure 2. (A) PLS-DA models—Two-component plots of GR and PR (green = GR, blue = PR) from UHPLC-MS metabolic profiles, demonstrate clear separation between the GR and PR groups. (B) Validation plots for the PLS-DA model obtained from 200 permutation tests (green = R2, blue = Q2): the intercepts of Q2 were <0.

Differences in Metabolites Between Stroke Patients With Good Recovery and Poor Recovery

The PLS-DA model was further analyzed to identify the serum metabolites associated with functional recovery of stroke. Through analysis of PLS-DA loadings, the metabolites contributing to sample discrimination were identified by VIP scores. The potential biomarkers were filtered and confirmed when their VIP scores are >1 (VIP > 1). Through analysis, five serum metabolites [threonine (VIP predicted values: 1.04), arginine (VIP predicted values: 1.47), glutamate (VIP predicted values: 1.89), proline (VIP predicted values: 1.97), and leucine-isoleucine (VIP predicted values: 2.46)] were screened out, which were closely related to the recovery level of stroke (VIP > 1) (Table 2). In contrast to those in the GR group, levels of glutamate (KEGG:C00025) and arginine (KEGG:C00062) were increased in the PR group, whereas levels of leucine-isoleucine (KEGG:C00123), proline (KEGG:C00148), and threonine (KEGG:C00188) were markedly decreased. The average normalized quantities of the differential metabolites in the GR and PR groups are shown in the heat map (Figure 3).

TABLE 2

Table 2. Results of the VIP predicted values of differential metabolites between stroke patients with GR and PR.

FIGURE 3

Figure 3. The heat map of the five different metabolites in the GR and PR groups. The colors changing from blue to red indicate the changes of characteristic metabolite content relative to the average normalized quantities of the metabolites in the GR group, blue plots indicate down-regulated metabolites, and red plots indicate up-regulated metabolites.

Identify Potential Predictive Biomarkers in Serum

Models were fit where five metabolites were considered in all possible combinations, including one single metabolite combination, two different metabolites in combination, and up to the combination of all five metabolites (Table 3). There were 31 combinations, and all combinations of the selected five metabolites for prediction of functional recovery after stroke were analyzed by cross-validation (Table 3).

TABLE 3

Table 3. Results of 10 independent replicates of 10-fold cross-validation of all combinations of the selected five metabolites for prediction of functional recovery after stroke.

Cross-validation is a well-known technique to choose tuning parameters of a model, while limiting the risk of overfitting (26, 27). As shown in Table 3, standard deviations (SD) of the AUC in the testing set are between 0.09 and 0.13, which indicate that the model is steady. Single metabolite did not offer good predictive value but displayed low sensitivity and specificity in both the training and test sets. The panel including proline, glutamate, and arginine had the highest values of mean AUC in the test set [0.835, 95% CI: (0.756, 0.913)] and mean AUC in the training set: [0.952, 95% CI: (0.951, 0.953)] of the 31 combinations (Table 3, Figure 4). Therefore, this panel was regarded as the best combination of predictive biomarkers out of the 31 combinations.

FIGURE 4

Figure 4. The ROC curve of the identified panel—The values of mean AUC of the panel including proline, glutamate, and arginine are 0.952 (95% CI: 0.756, 0.913) in the training set (A), and 0.835 (95% CI: 0.951, 0.953) in the test set (B).

As shown in Table 3, arginine was in all of the top 10 combinations ranked by values of mean AUC in the test set. Four out of the top five combinations ranked by values of mean AUC in the test set contained both glutamate and arginine. When leucine-isoleucine, threonine, and proline were combined with arginine, these new panels all had high predictive value for stroke functional recovery (Table 3). The models suggest that arginine is important in stroke recovery.

The Correlations Between Age and Metabolites

There were significant differences in age between the GR and PR groups. So, we investigated the correlation between age and identified metabolites in the two groups. As shown in Table 4 and Figure 5, there was little significant correlation between identified metabolites and age in the two groups.

TABLE 4

Table 4. The Pearson's correlation coefficients between age and identified metabolite.

FIGURE 5

Figure 5. Pearson's correlation coefficients graph between age and identified metabolites in the two groups (A) Arginine vs. Age, (B) Glutamate vs. Age, (C) Proline vs. Age, (D) Threonine vs. Age, (E) Leucine-isoleucine vs. Age.

Based on the above results, we recombined the age and the panel of proline, glutamate, and arginine. We analyzed the 4-variable panel, respectively, in the training and test sets. The values of mean AUC of the 4-variable panel was 0.962 (95% CI: 0.961, 0.963) in the training set, and 0.871 (95% CI: 0.819, 0.924) in the test set, which had higher sensitivity and specificity distinguished between patients with GR and PR compared with the identified panel from 31 combinations (Figure 6).

FIGURE 6

Figure 6. The ROC curve of the 4-variable panel including age, proline, glutamate, and arginine—The values of mean AUC of the panel including age, proline, glutamate, and arginine in the training set (A) is 0.962 (95% CI:0.961, 0.963) and 0.871 (95% CI:0.819, 0.924) in the test set (B).

Discussion

Amino acids are among the most important disturbed metabolites after stroke, but the association of amino acid levels with level of stroke recovery is not clear. Our main findings were that there was a significant difference in leucine-isoleucine, proline, threonine, glutamic acid, and arginine levels between the GR and PR groups. The panel of combined proline, glutamate, and arginine provided high sensitivity and specificity in prediction of functional recovery.

Our results showed that arginine levels in the PR group were significantly higher than those in the GR group. This is in line with evidence suggesting that high arginine concentration can induce neurotoxic substances (12, 14, 28). Studies have shown that arginine in plasma can cross the blood-brain barrier (BBB) into the brain (16, 28), and that there is a high correlation between serum arginine and arginine in cerebrospinal fluid (CSF) in stroke patients (16). As shown in the arginine metabolic pathway (Figure 7A), arginine is the precursor of nitric oxide (NO), and NO is synthesized from L-arginine by nitric oxide synthase (NOS), which includes neuronal NOS (nNOS), endothelial NOS (eNOS), and inducible NOS (iNOS). In the human brain, the synthesis of NO is mainly related to nNOS that lies in neurons (13). At physiological concentrations (EC₅₀ 1–4 nM) (29, 30), NO can regulate blood flow, relax blood vessels, and inhibit platelet aggregation, which are beneficial for the recovery of stroke dysfunction (28, 31). Several studies have shown that arginine supplementation contributes to stroke recovery, which is related to the PKC-mediated NO signaling pathway (32). However, recent studies indicate that high concentrations of arginine can increase oxidative stress in the general population (28). Arginase (also present in the brain) is an enzyme that catalyzes the conversion of arginine to ornithine, competing with eNOS (in the choroid plexus and vascular endothelium) for arginine (28). A high concentration of arginine could stimulate the expression and activation of arginase, and the maximal catalytic activity of arginase is higher than that of eNOS (15, 28). The increase of arginase activity will decrease the catalytic activity of arginine for NO production and lead to the uncoupling of eNOS, while uncoupled eNOS can induce the transfer of electrons from NADPH to oxygen molecules to form superoxide anions (${\text{O}}_2^{-}$). In vitro studies also found that a high concentration of arginine decreased the antioxidant capacity of brain tissue, which in turn increased oxidative stress, inhibited the activity of glutathione peroxidase in the brain tissue, induced the production of neurotoxic substances, and finally decreased the recovery of nerve function (15, 28).

FIGURE 7

Figure 7. (A) A higher glutamate concentration in the brain can induce excitotoxicity, and the increase of Ca2+ concentrations in the cell, which may cause mitochondrial dysfunction and increase oxidative stress and neuron death. A high arginine concentration can increase arginase activity, which leads to an uncoupling of eNOS and induces more neurotoxic substance (imaginary line). Excessive glutamate has excitotoxicity, and arginine can also increase oxidative stress and induce more neurotoxins by itself. Arginine and glutamate can be converted into each other and linked in their metabolism. When higher concentrations of arginine and glutamate come together, more neurotoxins such as peroxynitrite are produced. Adenosine monophosphate activated protein kinase (AMPK) also can be activated by peroxynitrite, so nNOS, peroxynitrite, and AMPK become one vicious cycle. All these will ultimately lead to neuronal death and functional impairment. (B) Both leucine-isoleucine and threonine can be converted to acetyl-coA and come into the TCA). TCA can provide ATP and AKG for the brain, which have an important effect on brain function recovery after stroke.

In this study, in addition to arginine, serum glutamate levels in patients with PR were also significantly higher than those in patients with GR. Glutamate is one of the most abundant free amino acids in the mammalian central nervous system (CNS) and is at the intersection of multiple metabolic pathways (33). In its physiological concentration, glutamate is crucial for various physiological processes, particularly synaptic transmission. But a higher glutamate concentration in the brain may trigger secondary brain injury following acute ischemic stroke, including neuron death, axonal injury, and mitochondrial dysfunction (34, 35). First of all, neuronal death, impaired energy supply, and increased oxidative stress caused by mitochondrial deregulation are not conducive to functional recovery after stroke. Secondly, glutamate receptors expressed on brain endothelial cells play an important role in regulating the function of the BBB (36). The excessive activation of glutamate receptors may lead to the abnormal expression and distribution of tight junction proteins in endothelial cells, resulting in the destruction of the BBB (36). In this case, the harmful substances may easily pass through the BBB and aggravate the death and injury of neurons. Thirdly, abnormal synaptic transmission induced by a high concentration of glutamate may decrease the synaptic plasticity and inhibit the recovery of neural function after stroke (Figure 6A) (36).

High concentrations of arginine and glutamate may lead to stronger neuroexcitotoxicity through synergistic action (37, 38). nNOS is Ca²⁺-dependent and will be activated to generate NO when calmodulin forms Ca²⁺/calmodulin complex with calcium (37). During the ischemia reperfusion (IR) phase of ischemic stroke, N-methyl-D-aspartic acid (NMDA) receptor-mediated excitotoxicity caused by increased glutamate levels may lead to calcium dysregulation and stimulate nNOS to produce more NO in the neurons (35). At the same time, this excitotoxicity may also trigger mitochondrial dysfunction and reduce the ability of mitochondria to resist oxidative stress, resulting in more reactive oxygen species (38, 39). Under oxidative stress, the production of S-nitrosoglutathione (GSNO) will be decreased, and the superoxide anion will react with NO to form a peroxynitrite anion (ONOO⁻), which is of high neurotoxicity (15, 40). ONOO⁻ can inhibit the activity of cytochrome c oxidase in the respiratory chain and destroy the electron transport-associated proteins in mitochondria, which leads to energetic failure and neuron death (15, 40). With the death of neurons and mitochondria, the oxidative stress levels in the brain may increase and lead to the production of more neurotoxic substances, such as nitrotyrosine (Figure 7A).

It was found that NMDA-receptor-mediated calcium dysregulation can activate nNOS by re-modifying its chemical structure, such as phosphorylation of Ser1412 (40). Adenosine monophosphate activated protein kinase (AMPK) is one energy sensor with a high expression of neurons (41). AMPK is activated when cellular energy is decreasing. The activation of AMPK may keep nNOS in a hyperactivated state via sustained phosphorylation of Ser1412 (40, 41). However, AMPK can also be activated by peroxynitrite, resulting in a vicious circle of nNOS, peroxynitrite, and AMPK (Figure 6A) (40). All of these will ultimately lead to neuronal death and functional impairment. Therefore, glutamate and arginine, both of which have high sensitivity and specificity to functional results, appeared four times in the top five combinations sorted according to the average value of AUC. This result is consistent with the synergistic effect of excessive glutamate and arginine on neurotoxicity. As shown in Figure 6A, arginine is also the precursor to glutamate. Therefore, an increase in arginine can cause an increase in glutamate, and arginine itself may also trigger this synergistic effect. This may partly explain why arginine could be found in the top ten combinations in the test set.

The levels of serum leucine-isoleucine and threonine in the PR group were significantly lower than those in the GR group. Studies have shown that leucine-isoleucine plays an active role in inhibiting an excessive glutamate concentration and excitotoxicity induced by stroke (5), and it is also an important component of energy metabolism (42). Threonine is also a fuel substrate that can be converted to pyruvic acid and become part of the tricarboxylic acid cycle (TCA) (43). TCA can provide adenosine triphosphate (ATP) and æ-ketoglutarate to the brain, while reducing leucine-isoleucine and threonine leads to a reduction in ATP and æ-ketoglutarate. The brain is an organ that consumes a lot of energy without storing any, so an abnormal energy metabolism may cause a brain disorder. æ–ketoglutarate is an important intermediate metabolite in the TCA cycle that can also improve ischemia-induced brain disorders (Figure 7B) (44).

The relationship between proline and stroke is unclear. It has been reported that the level of serum proline in patients with acute ischemic stroke is significantly lower than that in normal controls (7). He et al. further proved that the level of serum proline in the PR group was significantly lower than that in the GR group, which may be related to the fact that proline can enhance the stability of proteins and cell membranes (45).

Our results revealed that the serum levels of five amino acids (leucine-isoleucine, proline, threonine, glutamic acid, and arginine) differed significantly between GR and PR groups. The panel grouped by proline, glutamate, and arginine had the highest sensitivity and specificity in predicting recovery and functional outcomes of stroke. This finding suggests that the metabolomic process combined high level of neuroexcitoxicity with low level of cell membrane stability may worsen stroke recovery.

Study Limitations

The stratification method used in this study enhanced the power of detection by targeting the extremes; however, this study is still limited by relatively small sample size. Second, subjects from this study were from one single, urban rehabilitation hospital; therefore, the findings cannot be generalized to the general population. Third, in this study, the average age of the GR group was significantly younger than the average age of the PR group (p < 0.05). The age difference may contribute to the difference in metabolimic profiling. In addition, the age difference could also cause bias of the results. Further study with age-matched groups, multi-center and large sample size is warranted to confirm the findings.

Conclusion

1. High levels of glutamate and arginine are associated with PR after stroke, which is related to oxidative stress and excitotoxicity.

2. Leucine-isoleucine, threonine, and proline involved in energy metabolism are positively related to functional recovery after stroke.

3. Reduced leucine-isoleucine, threonine, and proline, as well as all or only one or two of them, were combined with elevated arginine into new panels, and these new panels illustrate a high predictive value for stroke functional recovery. The panel grouped by proline, glutamate, and arginine has the highest sensitivity and specificity on predicting recovery and functional outcomes of stroke of all 31 combinations. Age can increase sensitivity and specificity of this identified panel for stroke functional recovery.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ethics Statement

The study was approved by the ethics committee of Spaulding Rehabilitation Hospital. Written informed consent was not required as per local legislation and national guidelines.

Author Contributions

XW, TL, and QW contributed to the inception and design of the paper and writing of the manuscript. HS, SC, GL, AC, and PF contributed to patient recruitment and data collection. XW, QW, XL, and LW contributed to the statistical analysis and research design. All authors approved the final version of the manuscript.

Funding

This study was supported by the Harvard Catalyst/The Harvard Clinical and Translational Science Center (National Center for Research Resources and the National Center for Advancing Translational Sciences, National Institutes of Health award no. UL1 TR001102); Harvard University and its affiliated academic health care centers; and Fundamental Research Funds for the Central Universities (grant no. 12ykpy39). The content is solely the responsibility of the authors and does not necessarily represent the official views of Harvard Catalyst, Harvard University, and its affiliated academic health care centers, or the National Institutes of Health.

Conflict of Interest

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

Abbreviations

AMPK, adenosine monophosphate activated protein kinase; AUC, area under the receiver operating characteristic curve; CSF, cerebrospinal fluid; FIM, Functional Independence Measure; GR, good recovery; IR, ischemia reperfusion; MRFS, Montebello Rehabilitation Factor Score; NO, nitric oxide; NOS, nitric oxide synthase; PCA, principal component analysis; PLS-DA, partial least-squares discriminant analysis; PR, poor recovery; ROC, receiver operating characteristic; TCA, tricarboxylic acid cycle; UHPLC-MS, ultra high performance liquid chromatography—mass spectrometry; UV, unit variance; VIP, variable importance of projection.

References

1. Fan YL, Zhan R, Dong YF, Huang L, Ji XX, Lu P, et al. Significant interaction of hypertension and homocysteine on neurological severity in first-ever ischemic stroke patients. J Am Soc Hypertens. (2018) 12:534–41. doi: 10.1016/j.jash.2018.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Tobias DK, Lawler PR, Harada PH, Demler OV, Ridker PM, Manson JE, et al. Circulating branched-chain amino acids and incident cardiovascular disease in a prospective cohort of US women. Circ Genom Precis Med. (2018) 11:22–41. doi: 10.1161/CIRCGEN.118.002157

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Duan GF, Ye Y, Xu S, Tao W, Zhao S, Jin T, et al. Signal peptide represses GluK1 surface and synaptic trafficking through binding to amino-terminal domain. Nat Commun. (2018) 9:1–11. doi: 10.1038/s41467-018-07403-7

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Meng XE, Li N, Guo DZ, Pan SY, Li H, Yang C. High plasma glutamate levels are associated with poor functional outcome in acute ischemic stroke. Cell Mol Neurobiol. (2015) 35:159–65. doi: 10.1007/s10571-014-0107-0

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kimberly WT, Wang Y, Pham L, Furie KL, Gerszten RE. Metabolite profiling identifies a branched chain amino acid signature in acute cardioembolic stroke. Stroke. (2013) 44:1389–95. doi: 10.1161/STROKEAHA.111.000397

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ruiz-Canela M, Toledo E, Clish CB, Hruby A, Liang L, Salas-Salvadó J, et al. Plasma branched-chain amino acids and incident cardiovascular disease in the PREDIMED trial. Clin Chem. (2016) 62:582–92. doi: 10.1373/clinchem.2015.251710

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Goldstein LB. Effects of amphetamines and small related molecules on recovery after stroke in animals and man. Neuropharmacology. (2000) 39:852–9. doi: 10.1016/S0028-3908(99)00249-X

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Chen SH, Chang FM, Tsai YC, Huang KF, Lin CL, Lin MT. Infusion of human umbilical cord blood cells protect against cerebral ischemia and damage during heatstroke in the rat. Exp Neurol. (2006) 199:67–76. doi: 10.1016/j.expneurol.2005.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Wang D, Kong J, Wu J, Wang X, Lai M. GC–MS-based metabolomics identifies an amino acid signature of acute ischemic stroke. Neurosci Lett. (2017) 2:7–13. doi: 10.1016/j.neulet.2017.01.039

CrossRef Full Text | Google Scholar

10. Newgard CB. Metabolomics and metabolic diseases: where do we stand? Cell Metab. (2017) 25:43–56. doi: 10.1016/j.cmet.2016.09.018

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Liu P, Li R, Antonov AA, Wang L, Li W, Hua Y, et al. Discovery of metabolite biomarkers for acute ischemic stroke progression. J Proteome Res. (2017) 16:773–9. doi: 10.1021/acs.jproteome.6b00779

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Mapstone M, Lin F, Nalls MA, Cheema AK, Singleton AB, Fiandaca MS, et al. What success can teach us about failure: the plasma metabolome of older adults with superior memory and lessons for Alzheimer's disease. Neurobiol Aging. (2017) 51:148–55. doi: 10.1016/j.neurobiolaging.2016.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Hasan N, McColgan P, Bentley P, Edwards RJ, Sharma P. Towards the identification of blood biomarkers for acute stroke in humans: a comprehensive systematic review. Br J Clin Pharmacol. (2012)74:230–40. doi: 10.1111/j.1365-2125.2012.04212.x

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Armengou A, Hurtado O, Leira R, Obón M, Pascual C, Moro MA, et al. L-arginine levels in blood as a marker of nitric oxide–mediated brain damage in acute stroke: a clinical and experimental study. J Cereb Blood Flow Metab. (2003) 23:978–84. doi: 10.1097/01.WCB.0000080651.64357.C6

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Bottemiller KL, Bieber PL, Basford JR, Harris M. FIM scores, FIM efficiency, and discharge disposition following inpatient stroke rehabilitation. Rehabil Nurs. (2006) 31:22–5. doi: 10.1002/j.2048-7940.2006.tb00006.x

CrossRef Full Text | Google Scholar

16. Jové M, Mauri-Capdevila G, Suárez I, Cambray S, Sanahuja J, Quílez A, et al. Metabolomics predicts stroke recurrence after transient ischemic attack. Neurology. (2015) 84:36–45. doi: 10.1212/WNL.0000000000001093

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wen H, Alshikho MJ, Wang Y, Luo X, Zafonte R, Herbert MR, et al. Correlation of fractional anisotropy with motor recovery in patients with stroke after postacute rehabilitation. Arch Phys Med Rehabil. (2016) 97:1487–95. doi: 10.1016/j.apmr.2016.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Koh GC, Chen CH, Petrella R, Thind A. Rehabilitation impact indices and their independent predictors: a systematic review. BMJ Open. (2013) 3:1–12. doi: 10.1136/bmjopen-2013-003483

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Carayol M, Licaj I, Achaintre D, Sacerdote C, Vineis P, Key TJ, et al. Reliability of serum metabolites over a two-year period: A targeted metabolomic approach in fasting and non-fasting samples from EPIC. PLoS ONE. (2015) 10:1–10. doi: 10.1371/journal.pone.0135437

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Iwasaki M, Ishihara J, Takachi R, Todoriki H, Yamamoto H, Miyano H, et al. Validity of a self-administered food-frequency questionnaire for assessing amino acid intake in Japan: comparison with intake from 4-day weighed dietary records and plasma levels. J Epidemiol. (2016) 26:36–44. doi: 10.2188/jea.JE20150044

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Perl A, Hanczko R, Lai ZW, Oaks Z, Kelly R, Borsuk R, et al. Comprehensive metabolome analyses reveal N-acetylcysteine -responsive accumulation of kynurenine in systemic lupus erythematosus: implications for activation of the mechanistic target of rapamycin. Metabolomics. (2015) 11:1157–74. doi: 10.1007/s11306-015-0772-0

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Yuan M, Breitkopf SB, Yang X, Asara JM. A positive/negative ion–switching, targeted mass spectrometry–based metabolomics platform for bodily fluids, cells, and fresh and fixed tissue. Nat Protoc. (2012) 7:872–81. doi: 10.1038/nprot.2012.024

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Wang H, Liu W, Fang X, Wang H, Ma W, Dong H, et al. Effect of 405 nm low intensity irradiation on the absorption spectrum of in-vitro hyperlipidemia blood. Technol Health Care. (2018) 26:135–43. doi: 10.3233/THC-174302

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Karthik R, Hafila MIF, Saravanan C, Vivek N, Priyadarsini P, Ashwath B. Assessing prevalence of temporomandibular disorders among university students: a questionnaire study. J Int Soc Prev Community Dent. (2017) 7(Suppl 1):S24–9. doi: 10.4103/jispcd.JISPCD_146_17

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Sethi S, Brietzke E. Omics-based biomarkers: application of metabolomics in neuropsychiatric disorders. Int J Neuropsychopharmacol. (2015) 19:1–13. doi: 10.1093/ijnp/pyv096

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Dykstra MA, Switzer N, Eisner R, Tso V, Foshaug R, Ismond K, et al. Urine metabolomics as a predictor of patient tolerance and response to adjuvant chemotherapy in colorectal cancer. Mol Clin Oncol. (2017) 7:767–70. doi: 10.3892/mco.2017.1407

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Bouhaddani SE, Uh HW, Jongbloed G, Hayward C, Klarić L, Kiełbasa SM, et al. Integrating omics datasets with the OmicsPLS package. BMC Bioinformatics. (2018) 19:371–9. doi: 10.1186/s12859-018-2371-3

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Carvalho AM, Oliveira AA, Loureiro AP, Gattás GJ, Fisberg RM, Marchioni DM. Arginine intake is associated with oxidative stress in a general population. Nutrition. (2017) 33:211–5. doi: 10.1016/j.nut.2016.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Bellamy TC, Griffiths C, Garthwaite J. Differential sensitivity of guanylyl cyclase and mitochondrial respiration to nitric oxide measured using clamped concentrations. J Biol Chem. (2002) 277:31801–7. doi: 10.1074/jbc.M205936200

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Gibb BJ, Wykes V, Garthwaite J. Properties of NO-activated guanylyl cyclases expressed in cells. Br J Pharmacol. (2003) 139:1032–40. doi: 10.1038/sj.bjp.0705318

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Ma L, Hu L, Feng X, Wang S. Nitrate and nitrite in health and disease. Aging Dis. (2018) 9:938–45. doi: 10.14336/AD.2017.1207

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gulati P, Singh N. Neuroprotective mechanism of ischemic postconditioning in mice: a possible relationship between protein kinase C and nitric oxide pathways. J Surg Res. (2014) 189:174–83. doi: 10.1016/j.jss.2014.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Macrez R, Stys PK, Vivien D, Lipton SA, Docagne F. Mechanisms of glutamate toxicity in multiple sclerosis: biomarker and therapeutic opportunities. Lancet Neurol. (2016) 15:1089–102. doi: 10.1016/S1474-4422(16)30165-X

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Chamorro A, Dirnagl U, Urra X, Planas AM. Neuroprotection in acute stroke: targeting excitotoxicity, oxidative and nitrosative stress, and inflammation. Lancet Neurol. (2016) 15:869–81. doi: 10.1016/S1474-4422(16)00114-9

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Joshi DC, Tewari BP, Singh M, Joshi PG, Joshi NB. AMPA receptor activation causes preferential mitochondrial Ca2+ load and oxidative stress in motor neurons. Brain Res. (2015) 1616:1–9. doi: 10.1016/j.brainres.2015.04.042

CrossRef Full Text | Google Scholar

36. Hamilton NB, Kolodziejczyk K, Kougioumtzidou E, Attwell D. Proton-gated Ca2+-permeable TRP channels damage myelin in conditions mimicking ischaemia. Nature. (2016) 529:523–7. doi: 10.1038/nature16519

CrossRef Full Text | Google Scholar

37. Garry PS, Ezra M, Rowland MJ, Westbrook J, Pattinson KT. The role of the nitric oxide pathway in brain injury and its treatment—from bench to bedside. Exp Neurol. (2015) 263:235–43. doi: 10.1016/j.expneurol.2014.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Kozlov AV, Bahrami S, Redl H, Szabo C. Alterations in nitric oxide homeostasis during traumatic brain injury. Biochim Biophys Acta. (2017) 1863(10 Pt B):2627–32. doi: 10.1016/j.bbadis.2016.12.020

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Liu F, Lu J, Manaenko A, Tang J, Hu Q. Mitochondria in ischemic stroke: new insight and implications. Aging Dis. (2018) 9:924–37. doi: 10.14336/AD.2017.1126

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Khan M, Dhammu TS, Matsuda F, Singh AK, Singh I. Blocking a vicious cycle nNOS/peroxynitrite/AMPK by S-nitrosoglutathione: implication for stroke therapy. BMC Neurosci. (2015) 16:42–53. doi: 10.1186/s12868-015-0179-x

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Apoorv TS, Karthik C, Babu PP. AMP-activated protein kinase (AMPK) is decreased in the mouse brain during experimental cerebral malaria. Neurosci Lett. (2018) 662:290–4. doi: 10.1016/j.neulet.2017.10.054

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Zhang S, Zeng X, Ren M, Mao X, Qiao S. Novel metabolic and physiological functions of branched chain amino acids: a review. J Anim Sci Biotechnol. (2017) 8:10–21. doi: 10.1186/s40104-016-0139-z

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Ye JZ, Lin XM, Cheng ZX, Su YB, Li WX, Ali FM, et al. Identification and efficacy of glycine, serine and threonine metabolism in potentiating kanamycin-mediated killing of Edwardsiella piscicida. J Proteomics. (2018) 183:34–44. doi: 10.1016/j.jprot.2018.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Kovalenko TN, Ushakova GA, Osadchenko I, Skibo GG, Pierzynowski SG. The neuroprotective effect of 2-oxoglutarate in the experimental ischemia of hippocampus. J Physiol Pharmacol. (2011) 62:239–46.

PubMed Abstract | Google Scholar

45. He G, Wu C, Hunag J, Zhou R. Effect of exogenous proline on metabolic response of Tetragenococcus halophilus under salt stress. J Microbiol Biotechnol. (2017) 27:1681–91. doi: 10.4014/jmb.1702.02060

PubMed Abstract | CrossRef Full Text | Google Scholar