Introduction

Hypertension is a risk factor for the development of atherosclerosis and cardiovascular diseases as well as for renal failure. Obesity and hyperlipidemia increase the risk of cardiovascular disease as well, as well as of diabetes mellitus type 2, gout and cognitive decline, including Alzheimer’s disease. In addition to efforts to change life-style, antihypertensive drugs and lipid-lowering agents are prescribed. Due to the high prevalence of these risk factors, this medication ranks among the most frequently prescribed drugs in medical practice.

Beta-blockers are a heterogeneous class of antihypertensive drugs. Their main mode of action is the inhibition of noradrenaline and adrenaline β-adrenergic receptors, which diminishes the effect of the sympathetic nervous system on its target organ, especially on the smooth muscle tissue of the blood vessels and the heart. Most beta-blockers act on β1-, β2- and/or β3-adrenergic receptors, while some beta-blockers of the newer generation show an additional vasodilatory effect by, for example, blocking α1-adrenergic receptors [13]. Another important class of antihypertensive agents are angiotensin-I-converting enzyme (ACE) inhibitors. The ACE catalyzes the conversion of angiotensin-I to angiotensin-II which causes the muscles surrounding blood vessels to contract, thereby restricting the blood flow and increasing blood pressure. In addition, ACE cleaves and reduces the action of the vasodilator bradykinin and, thereby, increases the blood pressure. ACE inhibitors diminish the action of this enzyme and, thus, blood vessels dilate and blood pressure is reduced [1, 4]. A drug therapy with diuretics has long been one of the most important treatments of hypertension and heart failure. However, today, many diuretics have found their main role in areas outside of hypertension. Diuretics inhibit the reabsorption of sodium in the renal tubules and their ability to alter long-term sodium balance induces important hemodynamic changes that result in a reduction in peripheral resistance and a decrease in blood pressure [57].

Two important classes of lipid-lowering drugs are statins and fibrates. Statins lower the concentration of cholesterol in the blood by reducing the cholesterol biosynthesis. The rate-limiting enzyme of this biosynthesis is the 3-hydroxy-3-methyl-glutaryl-CoA reductase (HMG-CoA reductase; HMGCR), which is competitively inhibited by the statins [8, 9]. In contrast, fibrates decrease the blood levels of fatty acids and triglycerides. They activate the peroxisome proliferator-activated receptor (PPAR), mainly PPAR-α, which mediates the stimulation of the fatty acid β-oxidation [1013].

In this study we used a pharmacometabolomic approach. Pharmacometabolomics aims to identify metabolic traits that offer insights into the intended as well as the unintentional effects of drugs on the human organism. The power of pharmacometabolomics was already shown in several studies [1416]. One example is the study from Trupp et al. [14]. Examining the metabolic response to simvastatin treatment in 148 study participants they identified metabolites from multiple pathways not directly connected to the drug target which may explain the variation in therapeutic response.

For the identification of drug effects clinical drug testing is the standard approach. These trials are conducted under very controlled conditions and they are mainly focused on the objectives of the trials. The metabolomics approach for the analysis of medication effects in a big study population can provide a broader picture of the pathways affected by the respective drug by giving a snapshot of the physiological state of the organism. Using this kind of analysis, hypotheses about on-target as well as off-target effects can be generated, that reflect the drug action in the population and under everyday life conditions. A population based approach also allows the analysis of several drugs in the same study.

In this study, we used a mass spectrometry based non-targeted metabolomics approach to investigate the effect of antihypertensive and lipid-lowering agents. To this end, we analyzed the metabolic profiles of 1,762 participants from the KORA (Cooperative Health Research in the Region of Augsburg) study population for associations with the intake of these drugs (classification into cases/controls see Table 1).

Table 1 Characteristics of the population related to antihypertensives and lipid-lowering drugs

Materials and methods

Study population

The research platform KORA conducts population-based surveys and subsequent follow-up studies in the fields of epidemiology, health care research, health economics and genetics. A multitude of different parameters is provided, including medical history, life style factors, and socio-demographic variables. The dataset analyzed here was taken from the F4 study, which was conducted in 2006–2008 as a follow-up of the fourth KORA survey (S4; 1999–2001). For this follow-up study F4, 3,080 of the 4,261 participants of the S4 were reexamined 7 years after baseline examination. From this group, the metabolic profiles of 1,768 participants, aged between 32 and 77 years were measured. In total, 1,762 of these study participants provided information about medication with beta-blockers, ACE inhibitors, diuretics, statins and/or fibrates. To determine the use of medication within the last 7 days before the examination an instrument for database-assisted online collection of medication data (IDOM) was applied [17]. The drug classes were considered according to the recommendations of the German Hypertension Association [18]. The characteristics of the population related to these drug classes are shown in Table 1a–e. The usage of the agents within each drug class is given in the Supplemental Fig. 1. Informed consent was obtained from each study participant and all study protocols were approved by the local ethics committees.

Blood samples

For this study, we measured the metabolic profile in the blood serum of the study participants. To avoid variation due to circadian rhythm, blood was drawn in the morning between 8 and 10:30 am after overnight fasting (at least 8 h). Medication was taken in the morning as usual. Material was drawn into serum gel tubes, gently inverted twice and followed by 30 min resting at room temperature (18–25 °C) to obtain complete coagulation. The material was then centrifuged for 10 min (2750 g at 15 °C). Serum was aliquoted and stored at 4 °C, after which it was deep frozen to −80 °C on the same day until analysis of the metabolites.

Metabolite profiling and metabolite spectrum

The metabolites were measured by the US-company Metabolon Inc., a commercial supplier of metabolic analyses. Their platform integrates the chemical analysis, including identification and relative quantification, data reduction, and quality assurance components of the process. Two separate ultrahigh performance liquid chromatography/tandem mass spectrometry (UHPLC/MS/MS; positive and negative mode) injections and one gas chromatography/mass spectrometry (GC/MS) injection were done on this platform.

The UHPLC/MS/MS platform utilized a Waters Acquity UPLC and a ThermoFisher LTQ mass spectrometer, while for GC/MS analysis a Thermo-Finnigan Trace DSQ MS was used.

A standard library containing retention time, molecular mass to charge ratio (m/z), preferred in-source fragments and adducts as well as their associated MS/MS spectra for all molecules in the library, subsequently enables to identify a multitude of metabolites in resulting MS/MS data. A more detailed description can be found in the Supplementary Material.

In a recent publication, Metabolon showed that their analytical platform is able to perform relative quantitative analysis of analytical data in a high-throughput mode and that it identifies a broad spectrum of molecules with a high degree of confidence [19]. Further details of the platform technology are described in Suhre et al. [20]. The measured panel includes 295 metabolites from many relevant classes such as amino acids, carbohydrates, cofactors and vitamins, acylcarnitines, glyceropospholipids, lipids, nucleotides, small peptides and xenobiotics.

This platform has been successfully used in several studies, e.g. for the analysis of the adult human plasma metabolome [21] and for the identification of sarcosine as a biomarker for prostate cancer [22].

Statistical analysis

The statistical analysis system R (http://www.r-project.org/) and SPSS for Windows (Version 19.0. Chicago: SPSS Inc.) were used for the statistical analysis.

The analytic approach applied by Metabolon is semi-quantitative. This means that the added standards are mainly used for determining the retention time and not for the calculation of the metabolite concentrations. Thus, a relative intensity is measured and the measurement is sensitive to instrument parameters as well as fluctuations caused by maintenances like column change. These fluctuations are day dependent and the run day has to be considered for each individual. Thus, a run day normalization of the metabolic data was done: For each individual and each metabolite the data was first divided by the day median of the respective metabolite and then multiplied with the overall median of this metabolite.

To avoid false-positive associations, data points of metabolic traits that lay more than three standard deviations off the mean were excluded from further analysis. Metabolite concentrations were log-transformed since a test of normality showed that in most cases the log-transformed concentrations were closer to a normal distribution than the untransformed values.

For the identification of metabolites associated with medication use, a linear regression test was applied to model the association of each metabolite concentration and each of the five drug classes beta-blockers, ACE inhibitors, diuretics, statins and fibrates, respectively. As cofactors BMI, age, gender, HDL cholesterol, LDL cholesterol, total cholesterol, triglycerides, hypertension (>160/95 mmHg or known, medicated hypertension) and diabetes were chosen. To control for the effect of testing multiple hypotheses we used the Bonferroni correction. Thus, basing on a nominal level of 0.05, we used 3.39 × 10−5 [0.05/(295 metabolites * five drugs)] as an estimate of the significance level.

Results

For each drug class—beta-blockers, ACE inhibitors, diuretics, statins and fibrates—we analyzed the metabolic dataset to identify metabolites whose blood serum levels significantly differ between subjects taking and not taking the respective drug (p value smaller than 3.39 × 10−5, the estimated significance level after correction for multiple testing—see “Statistical analysis”). The results of the linear regression are given in Table 2.

Table 2 Results of the linear regression test with the strongest association to beta-blockers, ACE inhibitors, diuretics, statins or fibrates

For patients who took beta-blockers we observed increased concentrations of pyroglutamine, homocitrulline, salicylate, and acylcarnitins in the blood serum. In contrast, serotonin, fatty acids and 3-hydroxybutyrate [FA(4:0-OH)] were decreased. Their p values ranged between 1.9 × 10−9 and 2.2 × 10−5.

For the group of ACE inhibitors we identified four metabolites that significantly associated with the intake of these drugs. The p values showed a broad range from 2.4 × 10−80 up to 7.5 × 10−13. While levels of HWESASXX and des-arg(9)-bradykinin were higher in case of medication with ACE inhibitors, we found lower levels of phenylalanylphenylalanine and aspartylphenylalanine. Diuretics showed associations with increased serum levels of pseudouridine, C-glycosyltryptophan, glutaroylcarnitine [C5-DC] and urate. Additional metabolites with a p value smaller than 3.39 × 10−5, namely homocitrulline, HWESASXX (both increased) and phenylalanylphenylalanine (decreased) were already found to associate with beta-blockers and ACE inhibitors, respectively.

To identify associations between lipid-lowering drugs and metabolites we analyzed the effects of statins as well as fibrates. For the statins the resulting metabolites with the lowest p values were 1-arachidonoylglycerophosphocholine [LPC(20:4)], 1-arachidonoylglycerophosphoethanolamine [LPE(20:4)], isobutyrylcarnitine [C3-M], 1-docosahexaenoylglycerophosphocholine [LPC(22:6)], alpha-tocopherol, uridine (all increased), 7-alphahydroxy-3-oxo-4-cholestenoate, 1-palmitoylglycerophosphoinositol [LPI(16:0)], lathosterol and glycochenodeoxycholate (all decreased). For fibrates most of the significant metabolites showed a positive association: 2-hydroxyisobutyrate [FA(3:0-OH-M)], 3-dehydrocarnitine, riboflavin, pantothenate, indolelactate, carnitine, pipecolate and uridine. Only for one of the resulting metabolites—pyroglutamine—a significant negative association was detected. Pyroglutamine was already observed to associate with the intake of beta-blockers. However, in contrast to the intake of fibrates, the association between the beta-blockers and the concentration of pyroglutamine was positive.

Discussion

In this study we analyzed the effect of antihypertensive drugs and lipid-lowering drugs on the human metabolism. To this end, 295 metabolites were measured in the serum of 1,762 participants of the population-based KORA F4 study. We found hypothesis-generating associations with metabolites for four different drugs, however, not for diuretics. The results of the linear regression are given in Table 2. In the following we will discuss the main results for each drug class.

Beta-blockers associate with decreased serotonin and free fatty acid levels

Beta-blockers diminish the effect of the sympathetic nervous system on its target organ mainly by inhibiting the action of noradrenaline and adrenaline on β-adrenergic receptors. Among the metabolites with the lowest p values we found several fatty acids that were decreased with beta-blocker intake and some acylcarnitines that were increased. The increase in the concentration of acylcarnitines is also supported by the only nominally significant results. The decrease of free fatty acids agrees with the function of beta-blockers inhibiting the action of noradrenaline and adrenaline on β-adrenergic receptors. Since lipolysis is dependent on the action of these hormones, less triglycerides are broken down to free fatty acids [23, 24]. Vanhees et al. [25] also reported this effect for the beta-blocker bisoprolol which reduced the availability of plasma free fatty acids in healthy men.

Serotonin, a monoamine neurotransmitter, was significantly decreased in the blood of participants taking beta-blockers. Low plasma levels of serotonin were found in depressive patients [26, 27] and depression has been discussed as side-effect of beta-blockers [28]. Hence, the serotonin level variation might explain an adverse effect.

The metabolite with the strongest association to the intake of beta-blockers was pyroglutamine—a cyclic derivative of glutamine—that showed an increased concentration. In this context, it is noteworthy that an opposite behaviour, namely a decrease in the pyroglutamine concentration was associated with fibrate intake. However, further research on this metabolite is needed since too little is known so far about its physiological role.

In summary, our results for beta-blockers show associations that may indicate side effects of this medication.

Direct action of ACE inhibitors on its target ACE

The Angiotensin I-converting enzyme (ACE) cleaves not only angiotensin, but also several other polypeptides. One substrate of ACE is Cholecystokinin-8 (CCK-8) from which the dipeptide Asp-Phe-NH2 is cleaved [29]. In line with an inhibition of this process by ACE inhibitors, we found lower levels of aspartylphenylalanine (Asp-Phe) in users of ACE inhibitiors. Thus, one might expect that the substrate CCK-8 was increased in our patients on ACE inhibitors. Unfortunately, due to the used metabolite profiling technique CCK-8 was not in our metabolite panel and, therefore, we can only hypothesize. CCK-8 was shown to exert an antidiabetogenic action via increased insulin secretion [30]. Aguilar and Solomon suggested that ACE inhibitors may prevent the onset of diabetes [31] and Andraws and Brown could confirm this effect in a meta-analysis [32]. Thus, our observation of decreased levels of aspartylphenylalanine and the assumed elevated levels of CCK-8, suggest that this antidiabetogenic effect of ACE inhibitors is mediated by the inhibition of the cleavage of CCK-8. This hypothesis might be a starting point for further investigation.

The other dipeptide, phenylalanylphenylalanine, was also negatively associated with the intake of ACE inhibitors. It might be a product of the cleaving action of ACE as well, although the correspondent ACE substrate has not yet been identified. Candidate peptide substrates may now be searched for based on this information.

Another substrate of ACE is des-Arg(9)-bradykinin—the active metabolite of bradykinin—which causes blood vessels to dilate [33, 34]. Recent reports suggest an important role of lower bradykinin levels in the development of hypertension [35, 36]. Thus, our observation of an increased blood level of des-Arg(9)-bradykinin under ACE inhibitors indicates the direct action of this drug and strengthens the suggested role of bradykinin in the treatment with ACE inhibitors [37].

Analysis of diuretics shows a heterogeneous picture

The results from the analysis of metabolic profiles of patients on diuretics are heterogeneous and difficult to interpret. An increased level of pseudouridine might indicate RNA degradation and cell turnover as well as the remodelling process in heart failure [38]. Our findings of an increase of C-glycosyltryptophan in patients on diuretics might be explained by the role of tryptophan. In liver cirrhosis the concentration of tryptophan is increased and tryptophan plays an important role in chronic kidney disease [3941]. This is consistent with our observation of increased levels of glutaroylcarnitine [C5-DC], since glutaryl-CoA hydrogenase (GCDH) is essential in the degradation pathway of tryptophan [42]. A direct association between the metabolites pseudouridine and C-glycosyltryptophan and end stage renal disease (ESRD) was recently observed by Niewczas et al. [43]. They reported higher blood levels of these compounds in type 2 diabetes patients who progressed to ESRD within several years compared to diabetes patients without ESRD. Moreover, Yonemura et al. [44] observed a negative correlation between the C-glycosylated compound C-mannosyltryptophane and estimated glomerular filtration rate (eGFR; indicator of kidney function). An irreversible formation of homocitrulline (carbamoylation) is reported for end-stage renal disease and may explain the increased concentration of this metabolite [45].

Apart from their role as antihypertensives diuretics are used in cases of chronic kidney disease, heart failure and liver cirrhosis and, thus, we might see the metabolic traits of these diseases. Only the observed increased concentration of urate is a known side-effect of loop diuretics and thiazide [46].

Since more than 70 % of the participants on diuretics used the agent hydrochlorothiazide (HCT) we repeated the analysis only with HCT-users. The two resulting metabolites with significant associations with HCT are pseudouridine and C-glycosyltryptophan. These metabolites were also the top hits in the analysis of the diuretics which suggests that the results for diuretics mainly reflect the results for HCT.

Statins and their effects on the synthesis and degradation of cholesterol

Statins lower the concentration of cholesterol by inhibiting the HMG-CoA-reductase (HMGCR) which is the rate controlling enzyme of the biosynthesis of cholesterol. The levels of lathosterol and 7-alpha-hydroxy-3-oxo-4-cholestenoate were decreased in the blood from study participants treated with statins. Lathosterol is a precursor of cholesterol in the cholesterol biosynthesis, while 7-alpha-hydroxy-3-oxo-4-cholestenoate is a metabolite of cholesterol [47, 48]. So these two metabolites show direct drug effects. The increased concentration of 1-arachidonoylglycerophosphocholine [LPC(20:4)] and 1-arachidonoylglycerophosphoethanolamine [LPE(20:4)] and the decreased level of 1-palmitoylglycerophosphoinositol [LPI(16:0)] agree with previous reports of an increase in arachidonic acid [FA(20:4)] synthesis and a decrease of palmitic acid [FA(16:0)] levels in persons on statins. This effect is assumed to be due to a statin-induced enhancement of delta five desaturase (FADS1) activity [4951]. The benefit of the polyunsaturated fatty acid arachidonic acid is controversial. On the one hand arachidonic acid was reported as “beneficial in preventing and/or improving age-related declines in brain and cardiovascular system function” [52] as well as protective against oxidative stress in neurons [53]. On the other hand its role in inflammation is not fully clear, since anti- [54] as well as pro-inflammatory [55] effects have been described. A beneficial effect is achieved by the decrease of the 1-palmitoylglycerophosphoinositol [LPI(16:0)] carrying a saturated fatty acid which is considered as unhealthy.

In summary, for statins we observed mainly associations that are directly related to the action of this drug class.

Fibrates: degradation marker and increased carnitine levels

As a second group of lipid-lowering drugs, we analyzed the intake of fibrates for associations with metabolic traits. Though the group size for fibrate intake was quite small, associations with small p values were observed due to the large effect size of the association. However, the possibility of a type I error cannot fully be excluded. The metabolite with the strongest association to fibrate intake was the upregulated 2-hydroxyisobutyrate [FA(3:0-OH-M)]. 2-hydroxyisobutyrate probably is a metabolite of fenofibrate. Liu et al. [56] propose a degradation of fenofibrate where reduced fenofibric acid is metabolized to reduced 4-chloro-4′-hydroxybenzophenone. We presume that 2-hydroxyisobutyrate is released from reduced fenofibric acid during this reaction. If this is correct, 2-hydroxyisobutyrate might be used as a marker for the degradation of fenofibrate. In addition, 2-hydroxyisobutyrate was reported to be associated with prediagnostic gestational diabetes mellitus conditions [57] and obesity, and it might also be derived from the gut microbiome [58].

We also observed associations that might be related to side-effects of fibrates. 3-Dehydrocarnitine and free carnitine showed an increased concentration in the blood of patients on fibrates. This agrees with previous findings that fenofibrate increases the level of carnitine [59, 60]. The fact that 3-dehydrocarnitine is an intermediate in the carnitine biosynthesis may explain the increase of the blood concentration of this metabolite. Pyroglutamine was already the result with the lowest p value in the beta-blockers analysis, but it also showed a strong association with the intake of fibrates. While with beta-blockers the association was positive, the association with fibrates was negative. The question if beta-blockers and fibrates influence the same metabolic phenomenon in opposite directions needs further investigation.

Our results show that indications for drug related metabolic changes can be identified in population studies although the participants are more influenced by their environment (e.g. xenobiotics, smoking, food, physical activity) and much less controlled than during a clinical drug testing, e.g. in phase 1. BMI, age, gender, HDL cholesterol, LDL cholesterol, total cholesterol, triglycerides, hypertension and diabetes are known to influence the here measured metabolic profile. To avoid seeing their disturbing influence in the results we used them as cofactors for our statistical analysis. In addition, we checked if any of the metabolites we found to be associated with beta-blockers, ACE inhibitors, statins or fibrates (not diuretics) shows an association with hypertension or the total blood cholesterol level in the non-medicated (no antihypertensives and no lipid-lowering drugs) participants. Except for one metabolite associated with the intake of statins (alpha-tocopherol), we found no such association, meaning that the associations observed for those drugs are mainly drug- and not disease-related. Also, a calculation considering for each drug the remaining four drugs as cofactors yielded the same significant associations.

Limitations

Our sampling population was composed of Caucasians of European ancestry, so these results cannot be generalized to other ethnicities—effects of genetic variation on metabolism needs to be considered [20, 61, 62]. Because of the cross-sectional design of this study, we can only generate hypotheses based on associations—causal nature of the observed associations would require longitudinal data.

Major covariates, such as age, gender, BMI, HDL cholesterol, LDL cholesterol, total cholesterol, triglycerides, hypertension and diabetes were corrected for in the statistical analysis. However, in individual cases, interaction with other medication (e.g. with antiplatelet drugs or thyroid preparations) of the participants cannot be excluded.

To increase the statistical power the associations we reported in this paper were calculated for drug classes. It is not possible to trace the observed metabolic changes back to single pharmaceuticals. However, for each drug class specific associations could be clearly observed.

In summary (Fig. 1), for beta-blockers we observed associations possibly related to side-effects, including a decreased concentration of free fatty acids and a possible relation to depression via a decreased serotonin level. The intake of ACE inhibitors and statins associated with metabolites that reflect the direct action of the agent on its target. For fibrates we observed a possible breakdown product. Taken together, these results now provide a starting point for further functional research on the action and side-effects of these drugs.

Fig. 1
figure 1

Summary of the main drug-metabolite associations. Beta-blockers associate with an impaired lipid metabolism; the action of ACE inhibitors can be seen from associations with changed levels of substrate and products of ACE; statins associate with metabolites involved in the biosynthesis and degradation of cholesterol; fibrates associate with a hypothesized breakdown product of fenofibrate; pyroglutamine shows contrariwise associations with fibrates and beta-blockers