Next Article in Journal
Evaluation of Reliability of the Recomputed Nutrient Intake Data in the National Heart, Lung, and Blood Institute Twin Study
Next Article in Special Issue
Influence of Ellagitannins Extracted by Pomegranate Fruit on Disulfide Isomerase PDIA3 Activity
Previous Article in Journal
Identification of Six Phytochemical Compounds from Asparagus officinalis L. Root Cultivars from New Zealand and China Using UAE-SPE-UPLC-MS/MS: Effects of Extracts on H2O2-Induced Oxidative Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Synergistic Effects of Sanghuang–Danshen Bioactives on Arterial Stiffness in a Randomized Clinical Trial of Healthy Smokers: An Integrative Approach to in silico Network Analysis

1
Department of Nutritional Science and Food Management, Ewha Womans University, Seoul 03760, Korea
2
Department of Neurology, Ewha Womans University School of Medicine, Seoul 07985, Korea
3
Department of Bio and Brain Engineering, KAIST, Daejeon 34141, Korea
4
Department of Food Science and Technology, Seoul National University of Science and Technology, Seoul 01811, Korea
5
Department of Neurology, Eunpyeong St. Mary’s Hospital, The Catholic University of Korea, Seoul 06591, Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Nutrients 2019, 11(1), 108; https://doi.org/10.3390/nu11010108
Submission received: 18 December 2018 / Revised: 31 December 2018 / Accepted: 2 January 2019 / Published: 7 January 2019
(This article belongs to the Special Issue Dietary Bioactive Compounds and Human Health and Disease)

Abstract

:
The vascular endothelium is a favorite early target of cardiovascular risk factors, including cigarette smoking. Here, we investigated the synergistic effects of Sanghuang–Danshen (SD) bioactives on vascular stiffness in a controlled clinical trial of healthy chronic smokers (n = 72). Relative to placebo, 4-week SD consumption at 900 mg/day improves pulse wave velocity (p = 0.0497), reduces systolic blood pressure (peripheral, p = 0.0008; brachial, p = 0.0046; and ankle, p = 0.0066), and increases endothelial nitric oxide synthase activation (p < 0.0001). We then mapped all differential markers obtained from the clinical data, Affymetrix microarray, and 1H NMR metabolomics, together with 12 SD bioactives, onto the network platform termed the context-oriented directed associations. The resulting vascular subnetwork demonstrates that ellagic acid, caffeic acid, protocatechuic acid, cryptotanshinone, tanshinone I, and tanshinone IIA are linked to NOS3, ARG2, and EDN1 for vascular dilation, implicated with arginine/proline metabolism. They are also linked to SUCLG1, CYP1A1, and succinate related to the mitochondrial metabolism and detoxification, implicated with various metabolic pathways. These results could explain the synergistic action mechanisms of SD bioactives in the regulation of vascular endothelial dilation and metabolism, confirming the potential of SD in improving vascular stiffness and blood pressure in healthy smokers.

1. Introduction

The vascular endothelial cells can actively and reactively participate in responding to numerous environmental stimuli and behavioral risk factors, by changing structure and function [1,2]. Typically, alterations of vascular endothelial structure and function will result in aggravating already-existent stress and promoting atherothrombotic diseases, if not correctly managed [3]. Among various risk factors for cardiovascular disease, cigarette smoking stands out as a significant behavioral health concern [4]. Cigarette smoking is mainly related to the impairment of vascular endothelium-dependent dilation, and thereby alterations in arterial stiffness and blood pressure [4,5,6,7]. Importantly, however, smoking-induced vascular endothelial change is an early and reversible phenomenon [8], providing a unique chance to evaluate the beneficial effects of lifestyle modification or drug therapies against vascular endothelium-dependent impairments.
Over the last two decades, remarkable progress has been achieved in the development of drugs that regulate vascular endothelial dysfunction and atherothrombosis, but negative results have often transpired, due to severe and unintended adverse events [9]. Alternatively, increasing attention has been paid to a range of functional foods that help to maintain or regain vascular endothelial structure and function. Multiple components in foods are expected to be working in concert with various targets in various physiological processes to modify the entire vascular system. In our previous animal and human umbilical vein endothelial cell studies, we demonstrated that a pair of herbs, including an edible mushroom “Sanghuang” (Phellinus baumii) and a highly valuable herb rhizome “Danshen” (Salvia miltiorrhiza Bunge) (i.e., Sanghuang–Danshen (SD)), contained 12 signature bioactives and ameliorated vascular endothelial dysfunction [10]. Phenolics in Sanghuang [11] and tanshinones in Danshen [12] were considered as bioactive ingredients having various biological activities. One of the mechanisms proposed in our previous study was the augmentation of endothelium-derived nitric oxide (NO) availability through the activation of endothelial NO synthase (eNOS) phosphorylation.
To advance the previous findings, in the present study, we explored the hypothesis that daily consumption of SD would be useful in improving arterial stiffness, due to the synergistic mechanisms of SD bioactives on endothelium-dependent vasodilation. To prove this hypothesis, we assessed vasomotor responses to SD consumption, by measuring blood pressure and pulse wave velocity (PWV) in a clinical trial of healthy chronic smokers. We also measured eNOS phosphorylation in erythrocytes, 1H nuclear magnetic resonance (NMR)-based metabolomics in urine, and microarray gene expressions in peripheral blood mononuclear cells (PBMCs). We then mapped all effects observed in the clinical trial and the analytical data of SD bioactives onto an anatomical context-specific network, termed context-oriented directed associations (CODA), which was recently developed in our previous study [13].

2. Subjects and Methods

2.1. Test Material

The SD and a color-matched placebo were kindly provided by Pulmuone Foods Co., Ltd. (Seoul, Korea). The study product was prepared, as detailed in our previous publication [10]. Briefly, the dried fruit body of Sanghuang and the root of Danshen were extracted using water and 60% ethanol in sequential order and 75% ethanol, respectively. Then, each extract was combined, filtered, concentrated, and spray-dried to yield a powder. A total of 12 bioactives were identified from Sanghuang (caffeic acid, ellagic acid, fumaric acid, hispidin, and protocatechuic acid) and Danshen (cryptotanshinone, danshensu, rosmarinic acid, salvianolic acid A, salvianolic acid B, tanshinone I, and tanshinone IIA) in our previous study, using a high-performance liquid chromatography (HPLC) equipped with a diode-array detector coupled to an ion trap mass spectrometer (Agilent Technologies, Santa Clara, CA, USA) [10]. For the standardization purpose, protocatechuic acid and tanshinone IIA were chosen as marker components at the levels of 0.45–1.5 and 50 μg/g, respectively. The study product was packed in a gelatin capsule to give a dose of 0, 600, or 900 mg SD per day.

2.2. Study Population

Subjects were recruited through poster advertisements. Healthy smokers (≥10 cigarettes/day), aged between 20 and 70 years with smoking history for more than 1 year, were eligible for participation. Exclusion criteria were body mass index of ≤18.5 or ≥35 kg/m2; a history of body weight change of ≥10% in the previous 8 weeks; a history of alcohol drinking of ≥140 g/week; the use of any medication or dietary supplements; a history of platelet dysfunction, hypertension, stroke, diabetes, thyroid disease, and any other chronic diseases likely to interfere with study participation; a history of hypersensitivity to the test material; and pregnancy or breastfeeding. The study was approved by the Institutional Review Boards of Ewha Womans University Medical Center (Seoul, Korea) and conducted according to the Declaration of Helsinki. The study was also registered in the WHO International Clinical Trials Registry Platform under the identification number KCT0001642.

2.3. Study Design

The study followed a randomized, double-blinded, placebo-controlled design with placebo (P), low-dose SD (L), and high-dose SD (H) groups for 4 weeks. The sample size was estimated at 24 subjects/group to detect a difference in the thrombotic tendency among groups with a statistical power of 80%, based on our previous unpublished studies with a two-sided alpha of 0.05, allowing for an attrition rate of 20%. Before entering the run-in period, a list of flavonoid-rich foods and beverages were handed out to each subject to guide restricted diets, which might be a major confounder of the assessment of arterial wall stiffness. The subjects were instructed to maintain their usual diet and lifestyle, but to refrain from eating or drinking flavonoid-rich foods and beverages. Dietary intakes and physical activities were recorded 3 days a week, including 2 weekdays and 1 weekend day, using a smartphone application. The level of compliance was monitored at the end of the study by counting the remaining capsules. At the beginning and end of the intervention period, blood pressure and arterial wall stiffness were measured; and blood and urine samples were taken after an overnight fast of ≥8 h, for western blotting, and transcriptome and metabolome analyses.

2.4. Measurements of Blood Pressure and Arterial Wall Stiffness

Participants rested for 15 min in a seated position, and three consecutive blood pressure readings were recorded on the non-dominant upper arm, using a noninvasive oscillometric device (Dinamap ProCare 200, GE Medical Systems, Milwaukee, WI, USA). The brachial and ankle pulse wave velocity (baPWV) was measured using a volume-plethysmographic apparatus (VP-1000 Plus system, Omron Healthcare, Tokyo, Japan). Measurements of baPWV were repeated twice, and the mean value of the right and left sides was used as the baPWV value. The radial artery pressure waveform was recorded simultaneously, using an automated tonometry system (HEM-9000AI, Omron Healthcare, Tokyo, Japan), to estimate the central arterial pressure. All measurements were determined by one trained staff, who was blinded to the study design and objectives. The average value of two measures was considered for statistical analysis.

2.5. Western Blotting of eNOS Phosphorylation

In erythrocytes, western blotting analysis was performed to determine eNOS and eNOS phosphorylation at Ser1177. After lying erythrocytes in a lysis buffer, the supernatant was analyzed for protein concentration by the Bradford method (Bio-Rad, Hercules, CA, USA). Equal amounts of proteins were separated by 10% sodium dodecyl sulfate–polyacrylamide gel electrophoresis and transferred to a nitrocellulose membrane (Bio-Rad). Following incubation in 5% bovine serum albumin, membranes were probed with antibodies against eNOS, phosphorylated eNOS (Ser1177) (Cell Signaling Technology, Danvers, MA, USA), and β-actin (Santa Cruz Biotechnology, Santa Cruz, CA, USA). Membrane-bound primary antibodies were visualized by horseradish peroxidase-conjugated secondary antibodies (Santa Cruz Biotechnology), and immunoreactive band intensities were quantified using a ChemiDoc XRS imager equipped with Quantity One software (Bio-Rad).

2.6. Statistical Analyses

All statistical analyses were performed using SAS 9.4 (SAS Institute, Cary, NC, USA). Logarithmic transformation was performed on skewed variables. Differences in the baseline characteristics were examined by one-way analysis of variance (ANOVA) or the chi-square test. The effect of the SD on each parameter was expressed as the relative change between measurements obtained at baseline and those obtained at the 4-week follow-up within each group, to consider the possible differences at baseline, among groups. Statistical comparisons were based on the one-way analysis of covariance (ANCOVA) analysis of adjusted least squares (LS) means. p < 0.05 was considered significant.

2.7. Analysis of Transcriptome and Metabolome

The PBMCs of the P and H groups were separated from blood samples using Histopaque (Sigma, St Louis, MO, USA) density gradient centrifugation. Samples were sent to the DNA Link (Seoul, Korea) for transcriptome analysis, using Affymetrix GeneChip Human Gene 2.0 ST Arrays (Affymetrix, Santa Clara, CA, USA). RNA extraction, cDNA synthesis, hybridization, and scanning of the 2.0 GeneChip were performed according to standard Affymetrix protocols. The DAT files were converted into CEL files by Affymetrix Command Console software.
The NMR spectra for urine samples of the P and H groups were acquired using a Varian Unity Inova 600 MHz spectrometer (Varian, Palo Alto, CA, USA) at Pusan National University (Pusan, Korea). Acquisition parameters were as follows: a spectral width of 24,038.5 Hz; acquisition time of 3 s, and 128 transients. Additional conditions of a relaxation delay time of 1 s and saturation power of 4 were set, to suppress massive water peaks. NMR spectra were reduced to data using Chenomx NMR Suite program ver. 4.6 (Chenomx, Inc., Edmonton, Alberta, Canada). Metabolite concentrations were annotated and quantified manually in the NMR spectra, using the Chenomx NMR Suite Professional software package ver. 4.6 (Chenomx, Inc.). The metabolomics data were imported into SIMCA-P ver. 14.1 (Umetrics, Umeå, Sweden) for multivariate analysis. Principal component analysis (PCA), orthogonal partial least squares-discriminant analysis (OPLS-DA), heat map, and metabolic pathway analysis were performed using the MetaboAnalyst 4.0 software (http://metaboanalyst.ca).

2.8. Integration of All Relevant Data onto the CODA Network

To construct a comprehensive component–target–pathway network, we first assembled the component–target–phenotype by searching: (1) targets of signature bioactives in SD, and (2) related phenotypes and biological markers of endothelium-dependent vasodilation and arterial stiffness. We then linked the components, targets, and phenotypes, using the shortest paths between SD bioactives and markers in the vascular subnetwork [14]. The biological marker–phenotype network is constructed by manual curation. Next, we assembled the component–target–pathway by: (1) collecting all genes and metabolites in the vascular subnetwork, (2) performing pathway enrichment analysis in a group of genes and metabolites, by the Fisher’s exact test combined with a false discovery rate of <0.01, and (3) linking the component, targets, metabolites, and enriched KEGG pathways (http://www.genome.jp/kegg/pathway.html). Finally, the component–target–pathway networks were visualized using Cytoscape 3.6.1 (http://cytoscape.org/).

3. Results

3.1. Characteristics of Study Subjects

The Consolidated Standards of Reporting Trials (CONSORT) flow diagram for this study is shown in Figure 1. Of the 72 subjects who entered the study, 64 subjects completed the study without experiencing any serious adverse events. Eight subjects were lost to follow-up (n = 2) or withdrew from the study due to personal reasons (n = 6). Data were analyzed according to both intention-to-treat and per-protocol principles. Since the results of both analyses were comparable, only the results of the per-protocol analysis of 64 samples are presented. Demographic, anthropometric, and nutritional characteristics were similar among groups, but the recommended food score value was significantly lower in the L group than in the P group (p = 0.0496) (Table 1). The recommended food score is a simple dietary quality score modified for a Korean diet [15].

3.2. Changes in Arterial Wall Stiffness, Blood Pressure, and eNOS Activation

As shown in Figure 2a, there was a significant treatment effect on arterial wall stiffness assessed by baPWV (p = 0.0497). The L and H groups showed 21% and 46% reductions in baPWV relative to the P group, respectively. Consistent results were found with the H group for peripheral systolic blood pressure (SBP) (p = 0.0008), brachial SBP (p = 0.0046), and ankle SBP (p = 0.0066), although the difference in the central SBP was not significant among groups (Figure 2b). Consistent with these findings, when the phosphorylated-eNOS/eNOS ratio was normalized to 1 for the P group, the H group showed a significantly higher level of eNOS phosphorylation than the P group at 4-week post-treatment (p < 0.0001) (Figure 2c).

3.3. Changes in Transcriptomics and Metabolomics

Microarray analysis and 1H NMR metabolomics analysis were carried out to show the impact of the high-dose SD consumption on gene expressions in PBMCs and urinary metabolite profiles, respectively. One hundred and sixty-nine genes met a 1.3-fold cutoff and significance threshold of p < 0.05, and a total of 51 genes were identified. The heat map in Figure 3a represents an overall variation of gene expressions in the two groups (placebo versus SD), indicating clear differences by SD consumption. In the meantime, PCA and OPLS-DA plots of urinary metabolomics displayed a marked separation between placebo and SD groups (Figure 3b). The OPLS-DA model in the present study was well-fitted and displayed an acceptable predictive ability (R2Y = 0.793 and Q2Y = 0.447), suggesting that SD consumption led to changes in endogenous metabolites. The following variable importance in the projection analysis identified 21 metabolites that allowed for a distinct separation between placebo and SD groups (Table S1), as visualized in the heat map (Figure 3c). Of these, the levels of 8 metabolites (cis-aconitate, malonate, N-acetylglycine, O-acetylcholine, succinate, urea, valproate, and 2-oxoglutarate) were statistically different between the two groups (Table S2). The results were further analyzed to perform metabolic pathway analysis using MetaboAnalyst, identifying alanine/aspartate/glutamate metabolism, arginine/proline metabolism, D-glutamine/D-glutamate metabolism, the tricarboxylic acid (TCA) cycle, propanoate metabolism, and glycine/serine/threonine metabolism as the most differential pathways (Figure 3d).

3.4. Synergistic Mechanisms of SD Bioactives on Vascular Endothelial Dilation

To facilitate the visualization and interpretation of the complex relationships between SD bioactives, target proteins, and phenotypes, we mapped the 12 signature bioactives of SD and all data obtained from the clinical trial, metabolomics, and microarray analysis on the CODA network platform. The overall network consisted of 218 nodes and 497 edges (Table S3). The further grouping by vascular endothelial dilation demonstrated that six bioactives played essential roles in this subnetwork, by interacting with the 15 targets. The target protein names are detailed in Table S4. Notably, bioactives in Sanghuang (ellagic acid, caffeic acid, and protocatechuic acid) were recognized to be highly influential components, showing the highest number of target proteins relative to those in Danshen (cryptotanshinone, tanshinone I, and tanshinone IIA). Ellagic acid had 9 targets (ARG2, CALR, COX2, CYP1A1, JUN, NO53, PLA2G2A, SUCLG1, and succinate); caffeic acid (ALOX15, ALOX5, COX2, CYP2E1, EDN1, NO53, SUCLG1, and succinate) and protocatechuic acid (ALDH7A1, COX1, COX2, CYP1A1, CYP2E1, NO53, SUCLG1, and succinate) had 8 targets. In contrast, cryptotanshinone and tanshinone I shared 5 targets (ALDH7A1, CYP1A1, NO53, SUCLG1, and succinate) and tanshinone IIA had 4 targets (EDN1, NO53, SUCLG1, and succinate). Conversely, each target had a different number of associations with bioactives in SD. The NOS3, SUCLG1, and succinate received links with all six bioactives in SD. The ALDH7A1, CYP1A1, and EDN1 also had associations with bioactives from both Sanghuang and Danshen, although fewer in number. In contrast, the ALOX5, ALOX15, ARG2, CALR, COX1, COX2, CYP2E1, and JUN received links with bioactives exclusively from Sanghuang. To understand the underlying mechanisms, we also constructed in silico metabolic pathway network corresponding to vascular endothelial dilation. The results identified that the NOS3 was connected to the arginine and proline metabolism; SUCLG1 and succinate to the TCA cycle, alanine/aspartate/glutamate metabolism, and D-glutamate/D-glutamine metabolism; and the ALDH7A1 to the glycine/serine/threonine metabolism and propanoate metabolism (Figure 4a). The complete sequence of these six major metabolic pathways is visualized in Figure 4b, showing a general increase of measured metabolites in the SD consumption group compared with the placebo group.

4. Discussion

In this study, we hypothesized that SD exerts beneficial effects on improving vascular wall stiffness by the synergistic mechanisms of multiple bioactives on endothelium-dependent vasodilation in the human body. In support of this hypothesis, we demonstrated a decrease in baPWV and SBP and an increase in eNOS phosphorylation in healthy chronic smokers who consumed SD for 4 weeks relative to placebo consumption. Additionally, more fundamentally, we identified that six out of twelve bioactives in SD had direct interactions with the 15 targets in the vascular endothelial subnetwork, offering insight into the synergistic actions. To date, accumulating evidence, both in vitro and in vivo, has demonstrated that a complex mixture of bioactives in plants works in concert with various molecular targets in the body, thereby exerting synergistic effects for promoting health [16]. This concept is, however, hardly employed in clinical trials, due to the inherent limitation of traditional clinical assessments. To our knowledge, we are the first group to explain how multiple bioactives in SD exert synergistic effects on vascular relaxation in human clinical trials, by adopting an integrative approach to in silico database.
While the underlying molecular mechanisms remain unclear, there are data from human studies that convincingly demonstrate the associations between chronic smoking-induced hypertension and arterial stiffness with the increased oxidative stress [17,18] and impairment of vascular endothelium-dependent dilation [4]. Of the indices that have been introduced to quantify the elastic properties of arteries, the PWV is considered as a noninvasive and simple method in the clinical practice. The PWV can be calculated based on the pulse transit time and distance traveled by a pulse between two pressure catheters placed at a known distance from one another. The aortic PWV, which measures stiffness from the carotid to the femoral artery, is the most accurate. However, it is not ideal for routine use, because some subjects may feel uncomfortable when exposing the inguinal area [19]. Therefore, in this study, we determined the baPWV, which requires placing the blood pressure cuffs on the four extremities. Sugawara et al. [19] validated that baPWV provides qualitatively similar information to aortic PWV. Moreover, Hung et al. [20] reported that baPWV is useful to predict early atherosclerotic changes in the vascular tree. Recently, an increasing number of studies have identified the promising functional foods or bioactives for vascular endothelial health using this methodology [21,22,23].
Furthermore, we also measured changes in SBP, demonstrating a statistically lowered peripheral SBP level in the high-dose SD consumption group compared with the others. However, changes in central SBP were not statistically significant. This phenomenon may be explained by the so-called SBP amplification, which means that peripheral SBP level is higher than the central SBP level, due to the reduced diameter and also the increased stiffness in the distal vessels relative to the proximal vessels [24].
To support our hypothesis of the mechanism responsible for the observed arterial stiffness and blood pressure-lowering effects of SD, we further analyzed eNOS phosphorylation in erythrocytes. NO generated from erythrocytes contributes to the intravascular NO pool, playing a prominent role in the regulation of blood flow and pathogenesis of hypertension and stroke [25]. In the present study, there was a significant increase in eNOS activation by the high-dose SD consumption relative to the placebo or low-dose SD consumption. This result is consistent with our previous findings, showing a significant increase in eNOS phosphorylation and NO production in SD-treated human umbilical vein endothelial cells [10]. Moreover, we employed an integrative approach, by merging all data obtained from this study and the 12 signature bioactives in SD with the computational network platform termed CODA, to explore the synergistic actions of SD bioactives for vascular function. The resulting network demonstrated that ellagic acid, caffeic acid, protocatechuic acid, cryptotanshinone, tanshinone I, and tanshinone IIA were the responsible components, affecting the changes in vascular function in chronic smokers. These 6 bioactives were all linked to NOS3 encoding human eNOS, suggesting that there were convergent effects of bioactives on the endothelial NO production [26], which is implicated in the arginine/proline metabolic pathway. Ellagic acid in Sanghuang was also involved in the hydrolysis of arginine to ornithine and urea by regulating ARG2, which plays a critical role in NO metabolism [27]. This finding is compatible with other independent experiments that reported the modulatory effects of protocatechuic acid [28] and ellagic acid [29] on NO availability through eNOS phosphorylation. Meanwhile, the network analysis also identified that caffeic acid in Sanghuang, and tanshinone IIA in Danshen were linked to EDN1, which is another critical target related to vascular control, as a peptide with potent vasoactive properties [30]. Collectively, our data identify SD as exerting a positive role in improving arterial stiffness and blood pressure, which is likely to be attributed to increasing NO availability through stimulation of eNOS phosphorylation. Besides, recent studies have demonstrated that proline residues play a key role as gluten epitopes recognized by the immune system in the intestine [31]. Therefore, it was suggested that the degradation of proline residues might be important in reducing gluten intolerance [32]. Given that SD is implicated in arginine/proline metabolic pathway, we advance the hypothesis that SD bioactives may play a role in alleviating symptoms of gluten intolerance.
Our network analysis also captured the stimulation of mitochondrial metabolism and detoxification as additional potential biological events related to the effects of SD against chronic smoking. The above six bioactives were directly connected with SUCLG1, the alpha subunit of succinyl-coenzyme A ligase, implicated in the TCA cycle in mitochondria, alanine/aspartate/glutamate metabolism, and D-glutamate/D-glutamine metabolism, by catalyzing the conversion of succinyl CoA to succinate and free CoA [33]. Protocatechuic acid in Sanghuang, and cryptotanshinone and tanshinone I in Danshen were recognized to bind to ALDH7A1, implicated in the glycine/serine/threonine metabolism and propanoate metabolism. ALDH7A1 is an aldehyde dehydrogenase that is thought to play a crucial role in protecting against reactive aldehydes generated during oxidative stress and the metabolism of xenobiotics [34,35]. It is also involved in protecting cells from hyperosmotic stress, by the generation of osmolytes [34,36]. Elevated oxidative stress can be caused by the coupling of hyperosmotic stress to the production of reactive oxygen species within the cell [37]. Concurrently, ellagic acid and protocatechuic acid in Sanghuang, and cryptotanshinone and tanshinone I in Danshen were linked to CYP1A1, a detoxifying cytochrome P450 enzyme that is induced by polycyclic aromatic hydrocarbons found in cigarette smoking [38].
In the present study, we also recognized that caffeic acid, ellagic acid, and protocatechuic acid in Sanghuang, but no bioactives in Danshen were connected to COX and ALOX. Both COX and LOX are the key enzymes involved in synthesizing eicosanoids, including prostaglandins, leukotrienes, and thromboxanes that are mediators of some inflammatory and allergic conditions, and platelet aggregation [39]. Eicosanoids are generated primarily through an oxidative pathway from arachidonic acid [40]. However, in the current study, arachidonic acid metabolism was not recognized as significant in the metabolic pathway analysis. This result is probably, at least partly, due to the limitation of this study of using a sole analytical platform for measuring the metabolites. Although NMR has many advantages, no single analysis technique can achieve detection of all metabolites in a biological sample. The combined use of gas chromatography–mass spectrometry might be needed for the analysis of fatty acids that are thermostable and can be derivatized to create volatile compounds [41].

5. Conclusions

Here, we reported for the first time that SD consumption exhibits improvement of vascular stiffness and blood pressure in chronic smokers. We also demonstrated that subsequent mapping of all data obtained in a clinical trial, Affymetrix microarray, and 1H NMR metabolomics, together with signature bioactives in SD, allows understanding the complex synergistic actions of SD in the field of vascular dilation, mitochondrial metabolism, and detoxification. Mainly, in this study, the integration of metabolomics into systems biology was useful to explain the underlying mechanisms. It is a limitation of this research that only NMR was used for the metabolomics analysis. Utilizing a multiplatform analysis may accelerate the study of metabolomics, providing wider information to support the synergistic effects of bioactives.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6643/11/1/108/s1, Table S1: List of the 21 metabolites with the VIP > 1, Table S2: Significantly altered metabolites and their related metabolic pathways, Table S3: Overall network consisting of 218 nodes and 497 edges, and Table S4: Subnetwork comprising 6 bioactives and 15 targets.

Author Contributions

O.K., J.Y.K., T.-J.S., Y.-J.K., and D.L. conceived and designed the study; Y.L., T.-J.S., and Y.-J.K. conducted the clinical research; Y.L. and W.H. analyzed the data; O.K., Y.L., W.H., and D.L. interpreted the data; Y.L. and O.K. drafted the manuscript.

Funding

This research was funded by the Bio & Medical Technology Development Program of the National Research Foundation (NRF) funded by the Ministry of Science & ICT (NRF-2012M3A9C4048761).

Acknowledgments

We are grateful to the subjects in this study.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Galley, H.F.; Webster, N.R. Physiology of the endothelium. Br. J. Anaesth. 2004, 93, 105–113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Thuillez, C.; Richard, V. Targeting endothelial dysfunction in hypertensive subjects. J. Hum. Hypertens. 2005, 19 (Suppl. 1), S21–S25. [Google Scholar] [CrossRef] [PubMed]
  3. Viles-Gonzalez, J.F.; Fuster, V.; Badimon, J.J. Atherothrombosis: A widespread disease with unpredictable and life-threatening consequences. Eur. Heart J. 2004, 25, 1197–1207. [Google Scholar] [CrossRef] [PubMed]
  4. Ambrose, J.A.; Barua, R.S. The pathophysiology of cigarette smoking and cardiovascular disease. J. Am. Coll. Cardiol. 2004, 43, 1731–1737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Virdis, A.; Giannarelli, C.; Neves, M.F.; Taddei, S.; Ghiadoni, L. Cigarette smoking and hypertension. Curr. Pharm. Des. 2010, 16, 2518–2525. [Google Scholar] [CrossRef]
  6. Mazzone, A.; Cusa, C.; Mazzucchelli, I.; Vezzoli, M.; Ottini, E.; Ghio, S.; Tossini, G.; Pacifici, R.; Zuccaro, P. Cigarette smoking and hypertension influence nitric oxide release and plasma levels of adhesion molecules. Clin. Chem. Lab. Med. 2001, 39, 822–826. [Google Scholar] [CrossRef]
  7. Rhee, M.Y.; Na, S.H.; Kim, Y.K.; Lee, M.M.; Kim, H.Y. Acute effects of cigarette smoking on arterial stiffness and blood pressure in male smokers with hypertension. Am. J. Hypertens. 2007, 20, 637–641. [Google Scholar] [CrossRef]
  8. Kato, T.; Inoue, T.; Morooka, T.; Yoshimoto, N.; Node, K. Short-term passive smoking causes endothelial dysfunction via oxidative stress in nonsmokers. Can. J. Physiol. Pharmacol. 2006, 84, 523–529. [Google Scholar] [CrossRef]
  9. Johnson, S. Known knowns and known unknowns: Risks associated with combination antithrombotic therapy. Thromb. Res. 2008, 123 (Suppl. 1), S7–S11. [Google Scholar] [CrossRef]
  10. Lim, Y.; Lee, S.; Kim, J.Y.; Shin, J.-H.; Kwon, O. A phellinus baumii–based supplement containing salvia miltiorrhiza bunge improves atherothrombotic profiles through endothelial nitric oxide synthase and cyclooxygenase pathways in vitro and in vivo. J. Funct. Foods 2016, 24, 231–243. [Google Scholar] [CrossRef]
  11. Shon, M.Y.; Kim, T.H.; Sung, N.J. Antioxidants and free radical scavenging activity of Phellinus baumii (Phellinus of Hymenochaetaceae) extracts. Food Chem. 2003, 82, 593–597. [Google Scholar] [CrossRef]
  12. Yang, Q.; Zhang, X.L.; Li, X.Y.; Tang, W.K.; Zhang, J.X.; Fang, C.X.; Zheng, C.Y. Coupling continuous ultrasound-assisted extraction with ultrasonic probe, solid-phase extraction and high-performance liquid chromatography for the determination of sodium danshensu and four tanshinones in salvia miltiorrhiza bunge. Anal. Chim. Acta 2007, 589, 231–238. [Google Scholar] [CrossRef] [PubMed]
  13. Yu, H.; Jung, J.; Yoon, S.; Kwon, M.; Bae, S.; Yim, S.; Lee, J.; Kim, S.; Kang, Y.; Lee, D. Coda: Integrating multi-level context-oriented directed associations for analysis of drug effects. Sci. Rep. 2017, 7, 7519. [Google Scholar] [CrossRef] [PubMed]
  14. Yosef, N.; Ungar, L.; Zalckvar, E.; Kimchi, A.; Kupiec, M.; Ruppin, E.; Sharan, R. Toward accurate reconstruction of functional protein networks. Mol. Syst. Biol. 2009, 5, 248. [Google Scholar] [CrossRef] [PubMed]
  15. Kim, J.Y.; Yang, Y.J.; Yang, Y.K.; Oh, S.Y.; Hong, Y.C.; Lee, E.K.; Kwon, O. Diet quality scores and oxidative stress in korean adults. Eur. J. Clin. Nutr. 2011, 65, 1271–1278. [Google Scholar] [CrossRef] [PubMed]
  16. Schmidt, B.; Ribnicky, D.M.; Poulev, A.; Logendra, S.; Cefalu, W.T.; Raskin, I. A natural history of botanical therapeutics. Metabolism 2008, 57, S3–S9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Siasos, G.; Tousoulis, D.; Vlachopoulos, C.; Antoniades, C.; Stefanadi, E.; Ioakeimidis, N.; Andreou, I.; Zisimos, K.; Papavassiliou, A.G.; Stefanadis, C. Short-term treatment with l-arginine prevents the smoking-induced impairment of endothelial function and vascular elastic properties in young individuals. Int. J. Cardiol. 2008, 126, 394–399. [Google Scholar] [CrossRef]
  18. Wilkinson, I.B.; Hall, I.R.; MacCallum, H.; Mackenzie, I.S.; McEniery, C.M.; van der Arend, B.J.; Shu, Y.E.; MacKay, L.S.; Webb, D.J.; Cockcroft, J.R. Pulse-wave analysis: Clinical evaluation of a noninvasive, widely applicable method for assessing endothelial function. Arterioscler. Thromb. Vasc. Biol. 2002, 22, 147–152. [Google Scholar] [CrossRef]
  19. Sugawara, J.; Hayashi, K.; Yokoi, T.; Cortez-Cooper, M.Y.; DeVan, A.E.; Anton, M.A.; Tanaka, H. Brachial-ankle pulse wave velocity: An index of central arterial stiffness? J. Hum. Hypertens. 2005, 19, 401–406. [Google Scholar] [CrossRef]
  20. Hung, C.S.; Lin, J.W.; Hsu, C.N.; Chen, H.M.; Tsai, R.Y.; Chien, Y.F.; Hwang, J.J. Using brachial-ankle pulse wave velocity to associate arterial stiffness with cardiovascular risks. Nutr. Metab. Cardiovasc. Dis. 2009, 19, 241–246. [Google Scholar] [CrossRef]
  21. Böhm, V. Lycopene and heart health. Mol. Nutr. Food Res. 2012, 56, 296–303. [Google Scholar] [CrossRef] [PubMed]
  22. Grassi, D.; Desideri, G.; Necozione, S.; di Giosia, P.; Barnabei, R.; Allegaert, L.; Bernaert, H.; Ferri, C. Cocoa consumption dose-dependently improves flow-mediated dilation and arterial stiffness decreasing blood pressure in healthy individuals. J. Hypertens. 2015, 33, 294–303. [Google Scholar] [CrossRef] [PubMed]
  23. Siasos, G.; Tousoulis, D.; Kokkou, E.; Oikonomou, E.; Kollia, M.E.; Verveniotis, A.; Gouliopoulos, N.; Zisimos, K.; Plastiras, A.; Maniatis, K.; et al. Favorable effects of concord grape juice on endothelial function and arterial stiffness in healthy smokers. Am. J. Hypertens. 2014, 27, 38–45. [Google Scholar] [CrossRef] [PubMed]
  24. Herbert, A.; Cruickshank, J.K.; Laurent, S.; Boutouyrie, P.; Collaboration, R.V.f.A.M. Establishing reference values for central blood pressure and its amplification in a general healthy population and according to cardiovascular risk factors. Eur. Heart J. 2014, 35, 3122–3133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kolluru, G.K.; Siamwala, J.H.; Chatterjee, S. Enos phosphorylation in health and disease. Biochimie 2010, 92, 1186–1198. [Google Scholar] [CrossRef]
  26. Bredt, D.S. Endogenous nitric oxide synthesis: Biological functions and pathophysiology. Free Radic Res. 1999, 31, 577–596. [Google Scholar] [CrossRef]
  27. Wu, G.; Morris, S.M., Jr. Arginine metabolism: Nitric oxide and beyond. Biochem. J. 1998, 336, 1–17. [Google Scholar] [CrossRef]
  28. Kang, Z.; Zhu, H.; Jiang, W.; Zhang, S. Protocatechuic acid induces angiogenesis through pi3k-akt-enos-vegf signalling pathway. Basic. Clin. Pharmacol. Toxicol. 2013, 113, 221–227. [Google Scholar] [CrossRef]
  29. Berkban, T.; Boonprom, P.; Bunbupha, S.; Welbat, J.U.; Kukongviriyapan, U.; Kukongviriyapan, V.; Pakdeechote, P.; Prachaney, P. Ellagic acid prevents l-name-induced hypertension via restoration of enos and p47phox expression in rats. Nutrients 2015, 7, 5265–5280. [Google Scholar] [CrossRef]
  30. Li, H.; Louey, J.W.; Choy, K.W.; Liu, D.T.; Chan, W.M.; Chan, Y.M.; Fung, N.S.; Fan, B.J.; Baum, L.; Chan, J.C.; et al. Edn1 lys198asn is associated with diabetic retinopathy in type 2 diabetes. Mol. Vis. 2008, 14, 1698–1704. [Google Scholar]
  31. Malalgoda, M.; Simsek, S. Celiac disease and cereal proteins. Food Hydrocoll. 2017, 68, 108–113. [Google Scholar] [CrossRef]
  32. Albanese, L.; Ciriminna, R.; Meneguzzo, F.; Pagliaro, M. Gluten reduction in beer by hydrodynamic cavitation assisted brewing of barley malts. Lwt-Food Sci. Technol. 2017, 82, 342–353. [Google Scholar] [CrossRef]
  33. Landsverk, M.L.; Zhang, V.W.; Wong, L.C.; Andersson, H.C. A SUCLG1 mutation in a patient with mitochondrial DNA depletion and congenital anomalies. Mol. Genet. Metab. Rep. 2014, 1, 451–454. [Google Scholar] [CrossRef] [PubMed]
  34. Brocker, C.; Cantore, M.; Failli, P.; Vasiliou, V. Aldehyde dehydrogenase 7a1 (aldh7a1) attenuates reactive aldehyde and oxidative stress induced cytotoxicity. Chem. Biol. Interact. 2011, 191, 269–277. [Google Scholar] [CrossRef] [PubMed]
  35. O’Brien, P.J.; Siraki, A.G.; Shangari, N. Aldehyde sources, metabolism, molecular toxicity mechanisms, and possible effects on human health. Crit. Rev. Toxicol. 2005, 35, 609–662. [Google Scholar] [CrossRef]
  36. Brocker, C.; Lassen, N.; Estey, T.; Pappa, A.; Cantore, M.; Orlova, V.V.; Chavakis, T.; Kavanagh, K.L.; Oppermann, U.; Vasiliou, V. Aldehyde dehydrogenase 7a1 (aldh7a1) is a novel enzyme involved in cellular defense against hyperosmotic stress. J. Biol. Chem. 2010, 285, 18452–18463. [Google Scholar] [CrossRef]
  37. Schliess, F.; Görg, B.; Häussinger, D. Pathogenetic interplay between osmotic and oxidative stress: The hepatic encephalopathy paradigm. Biol. Chem. 2006, 387, 1363–1370. [Google Scholar] [CrossRef]
  38. Boström, C.E.; Gerde, P.; Hanberg, A.; Jernström, B.; Johansson, C.; Kyrklund, T.; Rannug, A.; Törnqvist, M.; Victorin, K.; Westerholm, R. Cancer risk assessment, indicators, and guidelines for polycyclic aromatic hydrocarbons in the ambient air. Environ. Health Perspect. 2002, 110 (Suppl. 3), 451–488. [Google Scholar]
  39. Funk, C.D. Prostaglandins and leukotrienes: Advances in eicosanoid biology. Science 2001, 294, 1871–1875. [Google Scholar] [CrossRef]
  40. Calder, P.C. Omega-3 fatty acids and inflammatory processes. Nutrients 2010, 2, 355–374. [Google Scholar] [CrossRef]
  41. Zhang, A.; Sun, H.; Wang, P.; Han, Y.; Wang, X. Modern analytical techniques in metabolomics analysis. Analyst 2012, 137, 293–300. [Google Scholar] [CrossRef] [PubMed]
Figure 1. CONSORT diagram of the study. The flow of subjects from enrollment to data analysis is shown, together with significant reasons for exclusion. All subjects who completed the study were analyzed. CONSORT, Consolidated Standards of Reporting Trials.
Figure 1. CONSORT diagram of the study. The flow of subjects from enrollment to data analysis is shown, together with significant reasons for exclusion. All subjects who completed the study were analyzed. CONSORT, Consolidated Standards of Reporting Trials.
Nutrients 11 00108 g001
Figure 2. Changes in vasomotor responses to SD consumption for 4 weeks (wk): (a) baPWV; (b) SBP (peripheral, ankle, brachial, and central); (c) Western blotting of eNOS phosphorylation at Ser1177 in erythrocytes. SD, Sanghuang–Danshen; baPWV, brachial and ankle pulse wave velocity; SBP, systolic blood pressure; eNOS, endothelial nitric oxide synthase. Values are expressed as mean ± SE (standard error). The different letters indicate significant differences at p < 0.05.
Figure 2. Changes in vasomotor responses to SD consumption for 4 weeks (wk): (a) baPWV; (b) SBP (peripheral, ankle, brachial, and central); (c) Western blotting of eNOS phosphorylation at Ser1177 in erythrocytes. SD, Sanghuang–Danshen; baPWV, brachial and ankle pulse wave velocity; SBP, systolic blood pressure; eNOS, endothelial nitric oxide synthase. Values are expressed as mean ± SE (standard error). The different letters indicate significant differences at p < 0.05.
Nutrients 11 00108 g002
Figure 3. Changes in the PBMC transcriptome and urinary metabolome after placebo or high-dose SD (900 mg) consumption for 4 weeks: (a) heat map visualization for the PBMC transcriptome; (b) PCA and OPLS-DA score plots derived from urinary NMR metabolites; (c) heat map visualization for the urinary 1H NMR metabolites; (d) metabolic pathway analysis overview showing altered urinary 1H NMR metabolic pathways by SD consumption. PBMC, peripheral blood mononuclear cell; SD, Sanghuang–Danshen; PCA, principal component analysis; OPLS-DA, orthogonal partial least squares-discriminant analysis; NMR, nuclear magnetic resonance.
Figure 3. Changes in the PBMC transcriptome and urinary metabolome after placebo or high-dose SD (900 mg) consumption for 4 weeks: (a) heat map visualization for the PBMC transcriptome; (b) PCA and OPLS-DA score plots derived from urinary NMR metabolites; (c) heat map visualization for the urinary 1H NMR metabolites; (d) metabolic pathway analysis overview showing altered urinary 1H NMR metabolic pathways by SD consumption. PBMC, peripheral blood mononuclear cell; SD, Sanghuang–Danshen; PCA, principal component analysis; OPLS-DA, orthogonal partial least squares-discriminant analysis; NMR, nuclear magnetic resonance.
Nutrients 11 00108 g003
Figure 4. Summary of synergistic mechanisms of SD bioactives in the vascular subnetwork: (a) component–target protein (round) or metabolite (triangle)–pathway network ( Nutrients 11 00108 i001); (b) complete sequence of six major metabolic pathways. The numbers in the box indicate the level of measured metabolites for the placebo (left) and SD (right) consumption group. SD, Sanghuang–Danshen; CA, caffeic acid; PCA, protocatechuic acid; EA, ellagic acid; CT, cryptotanshinone; TI, tanshinone I; and TIIA, tanshinone IIA. The target protein names are detailed in Table S4 (Supplementary Materials).
Figure 4. Summary of synergistic mechanisms of SD bioactives in the vascular subnetwork: (a) component–target protein (round) or metabolite (triangle)–pathway network ( Nutrients 11 00108 i001); (b) complete sequence of six major metabolic pathways. The numbers in the box indicate the level of measured metabolites for the placebo (left) and SD (right) consumption group. SD, Sanghuang–Danshen; CA, caffeic acid; PCA, protocatechuic acid; EA, ellagic acid; CT, cryptotanshinone; TI, tanshinone I; and TIIA, tanshinone IIA. The target protein names are detailed in Table S4 (Supplementary Materials).
Nutrients 11 00108 g004
Table 1. General characteristics of the participants 1.
Table 1. General characteristics of the participants 1.
CharacteristicPlacebo
(n = 21)
Low-Dose SD
(n = 23)
High-Dose SD
(n = 20)
p-Value
Gender (male/female)19/221/219/10.9204
Age (years)34.9 ± 2.533.0 ± 1.532.9 ± 2.80.7858
Cigarette smoking (per day)15.6 ± 1.615.3 ± 1.314.3 ±1.10.7909
Alcohol consumption (yes/no)14/713/1014/60.6262
Total physical activity (kcal/day)1995 ± 3392351 ± 6943223 ± 8680.4281
Waist (cm)87.8 ± 1.389.5 ± 2.284.2 ± 1.70.1234
Height (cm)173.9 ± 1.4174.1 ± 1.3172.6 ± 1.40.6906
Body weight (kg)73.4 ± 2.077.1 ± 2.769.6 ± 2.10.0855
Body mass index (kg/m2)24.2 ± 0.525.4 ± 0.923.3 ± 0.60.0934
Recommended food score22.4 ± 1.7 a16.7 ± 1.6 b20.8 ± 1.7 ab0.0496
Energy intake (kcal/day)1492 ± 831421 ± 941488 ± 910.8181
Carbohydrate (g/day)215 ± 10195 ± 14199 ± 130.5214
Protein (g/day)51.7 ± 6.053.2 ± 6.756.0 ± 6.10.7676
Fat (g/day)43.5 ± 7.545.7 ± 5.949.2 ± 6.20.6543
Sodium (mg/day)2770 ± 2372770 ± 2372675 ± 2780.3949
1 Data are shown as means ± SE (standard error). SD, Sanghuang–Danshen. The ANOVA and chi-square test were used to compare differences between groups. The different letters indicate significant differences at p < 0.05.

Share and Cite

MDPI and ACS Style

Lim, Y.; Song, T.-J.; Hwang, W.; Kim, J.Y.; Lee, D.; Kim, Y.-J.; Kwon, O. Synergistic Effects of Sanghuang–Danshen Bioactives on Arterial Stiffness in a Randomized Clinical Trial of Healthy Smokers: An Integrative Approach to in silico Network Analysis. Nutrients 2019, 11, 108. https://doi.org/10.3390/nu11010108

AMA Style

Lim Y, Song T-J, Hwang W, Kim JY, Lee D, Kim Y-J, Kwon O. Synergistic Effects of Sanghuang–Danshen Bioactives on Arterial Stiffness in a Randomized Clinical Trial of Healthy Smokers: An Integrative Approach to in silico Network Analysis. Nutrients. 2019; 11(1):108. https://doi.org/10.3390/nu11010108

Chicago/Turabian Style

Lim, Yeni, Tae-Jin Song, Woochang Hwang, Ji Yeon Kim, Doheon Lee, Yong-Jae Kim, and Oran Kwon. 2019. "Synergistic Effects of Sanghuang–Danshen Bioactives on Arterial Stiffness in a Randomized Clinical Trial of Healthy Smokers: An Integrative Approach to in silico Network Analysis" Nutrients 11, no. 1: 108. https://doi.org/10.3390/nu11010108

APA Style

Lim, Y., Song, T. -J., Hwang, W., Kim, J. Y., Lee, D., Kim, Y. -J., & Kwon, O. (2019). Synergistic Effects of Sanghuang–Danshen Bioactives on Arterial Stiffness in a Randomized Clinical Trial of Healthy Smokers: An Integrative Approach to in silico Network Analysis. Nutrients, 11(1), 108. https://doi.org/10.3390/nu11010108

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop