 Research article
 Open Access
 Published:
A mathematical model to estimate cholesterylester transfer protein (CETP) triglycerides flux in human plasma
BMC Systems Biology volume 13, Article number: 12 (2019)
Abstract
Background
Cholesterylester transfer protein (CETP) modulates the composition of various lipoproteins associated with cardiovascular disease. Despite its central role in lipoprotein metabolism, its mode of action is still not fully understood. Here we present a simple way to estimate CETPmediated lipid fluxes between different lipoprotein fractions.
Results
The model derived adequately describes the observed findings, especially regarding low and high dense lipoproteins (LDL and HDL), delivering correlation coefficients of R^{2} = 0.567 (p < 0.001) and R^{2} = 0.466 (p < 0.001), respectively. These estimated fluxes correlate best among all other measured concentrations and ‘lipid per lipoprotein’ ratios to the observed fluxes.
Conclusion
Our model approach is independent of CETPaction’s exact mechanistic mode. It is simple and easy to apply, and may be a useful tool in revealing CETP’s ambiguous role in lipid metabolism. The model mirrors a diffusionlike exchange of triglycerides between lipoproteins. Cholesteryl ester and triglyceride concentrations measured in HDL, LDL and VLDL are sufficient to apply the model on a plasma sample.
Background
During their retention in human plasma the lipid load of lipoproteins is altered by several enzymes and receptors. Important enzymes are Cholesterylester transfer protein (CETP), lipases, lecithincholesterol acyl transferase (LCAT), and phospholipid transferprotein (PLTP). CETP mediates the transport of lipids, namely phospholipids (PLs), triglycerides (TG) and cholesteryl esters (CE) from one lipoprotein to another. TG and CE are strictly hydrophobic, are confined to the lipoprotein core, and cannot leave a lipoprotein by diffusion via the aqueous phase. However, the dynamics of CETPmediated TG and CE fluxes resemble the dynamics of nonmediated diffusion [1]. Concentration gradients lead to a net flux of TG from TGrich lipoproteins, namely very low density lipoproteins (VLDL) and intermediate density lipoproteins (IDL) to the more CErich lipoproteins low density lipoproteins (LDL) and high density lipoproteins (HDL). Correspondingly, CE is transported from LDL/HDL to VLDL/IDL. Consequently, in a drawn blood sample in vitro, the CE to TG ratio in all lipoprotein cores will ultimately reach equilibrium after within a certain period [2]. The CETPmediated flux of CE from HDL to VLDL contributes to reverse cholesterol transport [3]. Although CETP is a key player in human lipoprotein metabolism, its impact on atherosclerosis and hypertension is not well understood. CETP inhibition has been a pharmaceutical approach to increase HDL cholesterol while decreasing LDL cholesterol and expecting a reduced risk of cardiovascular diseases. However, the clinical results of CETPinhibition have been far from promising and deeper understanding of CETP action and the CETPmediated dynamics between lipoproteins is urgently needed [4].
The exact molecular mechanism of lipid exchanges among lipoproteins via CETP has not been clarified. Two hypothesis are being debated: a) the ternary complex [5, 6], or b) a shuttle model [7].
In the ternary complex model, CETP binds to a lipoprotein followed by binding to another lipoprotein to a different side. In this situation CETP, forms a hydrophobic tunnel, enabling the exchange of neutral lipids between both lipoproteins (channel mechanism).
In the shuttle model, CETP docks to a single lipoprotein to receive one or two hydrophobic lipids, separates, and docks to another lipoprotein.
There are several methods of measuring CETP activity and mass concentration [8, 9]. As Lagrost [10] described, measurements of CETP activity may be classified into 2 categories: First, by measuring the net mass transfer of CE between lipoprotein fractions like VLDL, LDL and HDL, and second by using radiolabeled CE and observing its distribution over time among HDL and Apolipoprotein B100 (ApoB) containing lipoproteins. Transport rates depend on various factors like the CETP mass concentration, number and composition of lipoproteins in plasma.
Aim of this investigation is to present a mathematical model that is independent of detailed molecular mechanisms, establishes a causal relationship of CETPmediated fluxes between lipoproteins, and finally is able to estimate them.
Results
Subjects
Plasma samples in experiment 1 were divided into 4 subgroups: ‘normal’ (n = 43), ‘high TG’ (n = 18), ‘high LDL’ (n = 8) and ‘low HDL’ (n = 22) according to the methods described. Table 1 summarises the characteristics of all subjects.
Measurement imprecision
The data used for model validation are based on small changes due to the 1 h storage at 37 °C in lipid mass of lipoprotein fractions (Fig. 1). Therefore, measurement imprecision was considered. It is less than 4% for all lipoprotein fractions studied. In most cases the expected changes after one hour at 37 °C are relatively small compared to baseline. Comparing the relative mean changes in TG mass to the corresponding coefficients of variation, the change in VLDL is relatively low, while it is high in LDL and HDL.
We used one hour as incubation duration, as longer durations may lead to nonlinear behavior of lipoprotein associated reactions happening in vitro and hence to data, which is less comparable to the in vivo situation.
Changes in CE and TG between plasma stored at 4 °C and 37 °C for one hour, [μmol/L]. Median (1st, 3rd quartile) of differences between ‘37 °C value’‘4 °C value’. Experiment 1: n = 91 samples of various metabolic states. Experiment 2: n = 11 plasma samples with no inhibition (normal), with inhibition of LCAT (LCAT) and with inhibition of CETP (CETP). Wilcoxon signed rank test * < 0.05, ** < 0.001.
Considering measurement imprecision, it is obvious that the data of TGredistribution from a single plasma sample at 37 °C may be too noisy to yield acceptable inferences of real TG fluxes mediated by CETP. However, sufficiently large collectives (n ≥ 10) may enable us to test our model.
CETP action
Figure 1 summarises changes in TG and CE concentrations in VLDL, IDL, LDL and HDL in experiments 1 and 2. The loss of TG in VLDL, as well as the increase in TG in LDL and HDL and the increase in CE in VLDL are significant in experiment 1. If CETP is inhibited (experiment 2), these changes are much smaller or even zero and lose their significance. If LCAT is inhibited the changes in TG in VLDL, LDL and HDL remain significant.
In experiment 1, there is a small but significant loss of 1% in VLDL ApoB (Wilcoxon: p = 0.001), but no significant corresponding changes in IDL or LDL. We assume that this ApoB is recovered in IDL or LDL but due to measurement imprecision we cannot statistically capture this fact.
Correlations
Figure 2 illustrates our model’s capability to estimate TG net flux via CETP in LDL, HDL and VLDL, as well as CE net flux in VLDL (assuming an equimolar exchange of TG and CE): Especially for HDL and LDL, the estimation correlates quite well to the measured changes. For VLDL, ΔCE correlates better than the corresponding ΔTG.
N = 91, subdivided into 4 metabolic states: n = 43 ‘normal’ (plus), n = 8 ‘high LDL’ (red triangle), n = 18 ‘high TG’ (blue diamond) and n = 22 ‘low HDL’ (green cross). Observed change: in TG (or CE) [μmol/L/h] in plasma stored at 37 °C. A: TG change HDL, B: TG change LDL, C: TG change VLDL, D: CE change VLDL (assuming equimolar exchange). Pearson’s correlation coefficient R^{2}.
Table 2 lists the correlations between the measured fluxes of TG (and CE) to VLDL, IDL, LDL and HDL to the estimated fluxes. Additionally, the best correlation out of all our other measured parameters (including lipid and apolipoprotein masses, and the resulting lipid per apolipoprotein ratios) are listed. For HDL and LDL, model predictions are the best correlations among all tested correlations.
Note that in experiment 1 considering ΔTG, there is no significant correlation between HDL and LDL (R = 0.174, p = 0.098), but rather between VLDL and LDL (R = − 0.248, p = 0.018) as well as between VLDL and HDL (R = − 0.478, p < 0.0001), respectively.
In the VLDL fraction, there is no significant difference between LCAT and nonLCATinhibited plasma in ΔCE. This suggests that LCAT activity in VLDL in experiment 2 is negligible.
Assuming an equimolar exchange between TG and CE in VLDL and no significant LCAT action on VLDL, the correlation between ΔTG and ΔCE in VLDL appears to be relatively weak (R = − 0.160, p = 0.130).
In contrast to the three LDL subfractions, the TG flux in the three measured HDL subfractions resulted in a much poorer fit compared to total HDL. This is not surprising, as the observed data are corrupted by the significant redistribution of ApoA1 mass among those fractions during the 1 h storage at 37 °C [11].
CETP_{TG} and r_{cetp}
Not surprisingly, correlations between the total calculated surface (as described in the method section) and corresponding PL concentration in HDL, IDL, LDL and VLDL are strong (R = 0.937, 0.994, 0.989 and 0.996, respectively). However, using calculated surfaces leads to slightly better results than using PL mass. Distribution of surfaces of lipoprotein (sub)fractions is displayed in Fig. 3.
Calculated surfaces median (1st, 3rd quartile) of lipoprotein fractions in the different subgroups (‘normal’: white; ‘high LDL’: light grey; ‘high TG’: dark grey; ‘low HDL’: black).
Figure 4 illustrates quartiles of the predicted and observed TG net fluxes for VLDL, IDL, LDL and HDL. Note that there is not only good correlation between the model and the observed data in each lipoprotein fraction, but the model is also capable of predicting the proportion of TG mass distributed among the lipoprotein fractions by CETP. The corresponding brokendown fluxes of TG into and out of the fractions are displayed in Fig. 5.
Predicted TG net flux (black) and observed flux (white) of TG; Median (1st, 3rd quartile), 4 groups, n = 43/8/18/22, respectively.
Predicted median (1st, 3rd quartile) flux of TG entering via heteroexchange (white), leaving via heteroexchange (grey) and the resulting net flux (black) in ‘normal’ (n = 43).
Figure 6 displays the v/v TG/(TG + CE) threshold at which –according to the model there is no CETP mediated TGnet flux, a ratio that is here defined as “isotransfer point” (ITP). The here called ITP equals the parameter CETP_{TG}. Considering all four metabolic situations investigated, the ITPs are always located between the corresponding means of the IDL and LDL fraction.
In experiment 1 the ITP shifts among the four subgroups: ‘high LDL’ 0.17 (0.16, 0.19), ‘normal’ 0.18 (0.15, 0.21), ‘low HDL’ 0.26 (0.21, 0.29), and ‘high TG’ 0.30 (0.26, 0.34). The correlation between CETP_{TG} and the ratio TG/(TG + CE) in plasma appears to be remarkably strong (R = 0.968, p = 3.7*E^{− 55}).
The ITP equals a functional characteristic of lipoproteins in a defined metabolic state and shifts from LDL in normolipidemia to IDL in hyperlipidemia.
As described in the model the rate for CETP interacting with a lipoprotein in plasma per hour per lipoprotein surface in\( \raisebox{1ex}{${m}^2$}\!\left/ \!\raisebox{1ex}{$l$}\right. \) is estimated as r_{cetp} = 10^{− 5}*4.116.
Discussion
Here we present a relative simple model for TG redistribution by CETP among lipoproteins in human plasma. Our model is based on the model of Potter et al. [1], which itself is based on the work of Morton [12, 13]. We relied on the data on TG composition of lipoproteins in plasma stored for one hour at 37 °C of n = 91 samples compared to baseline to calculate the CETPmediated TGredistribution among lipoproteins.
A slightly different version of the model presented here was already introduced in our group’s published article [11]. The main difference from that model is that in the present work, we are relaxing the assumption that exchanges are equimolar. In our earlier model, only the means of TG in LDL (n = 27, normolipidemic) were used for model calibration. However, in this study, we not only acquired more data from both healthy and pathologic lipoprotein profiles to test the model, we also expanded our model to include VLDL and HDL. Moreover, we do not only compare measured and estimated medians but also considered corresponding correlations.
The model’s structure
Potter’s model describes an equimolar shuttle exchange of TG and CE by CETP between lipoproteins on the molecular level. Hence, on a formal basis, CETP is discriminated into a free and a lipoproteinbound fraction. Correspondingly, lipoproteins are discriminated into free and CETPbound. The model presented here simplifies and generalises this model. Our model is independent of the exact biochemical mode of CETP action (shuttle or ternary complex)  an obvious advantage, since the mode of action has still not been elucidated. Our generalised model is compatible with both the shuttle and the ternary complex models. Although differing on the molecular level, both approaches result in the same dynamics on the macroscopic level (by the law of large numbers).
As we point out later, our data do not support equimolar change in all fractions, thus we model the TG change only. However, at least in VLDL, our model seems to predict the CE flux to a strong degree. Furthermore, we are not modelling free and bound lipoproteins and CETPmolecules explicitly. The ‘CETP binding to lipoprotein’reaction is disregarded, as the dynamics of those reactions are only relevant on a mechanistic level.
Note that there are several mathematical model approaches in the field of human lipid metabolism dealing with for example topics like oxidized lipoproteins [14], LDL endocytosis [15] or the cholesterol biosynthesis pathway [16].
Comparison to other models
There are other models besides Potter’s describing in vivo CETP dynamics. Figure 7 summarises two thereof. The model presented by Lu et al. [17] describes the flux of CE from HDL to VLDL as well as CE’s flux from LDL to HDL and viceversa. In Lu’s model the fluxes are simple firstgrade reactions\( {R}_i:C{E}_{donor}\overset{k_i}{\to }C{E}_{acceptor} \) (i = 1, 2, 3) with three different corresponding reaction rates. Applying our own data to both models, the observed ΔCE correlates much weaker to Lu’s predicted ΔCE in VLDL than to our predicted ΔCE: R = 0.315, p = 0.002 vs. R = 0.705, p = 6.4E^{− 15}. As we are not assuming equimolar exchange and just model TG, we cannot compare Lu’s model directly to ours. However, considering the LCATinhibited data from experiment 2, we detected no correlation between the CE enrichment in HDL and CE mass in LDL (or VLDL).
Hübner’s model [18] considers TG and cholesterol (CH), which corresponds to FC + CE. A shuttle transport is modeled in which CETP is either free or TG or CHloaded. Let us replace CH by CE to better compare their model to ours. Here, five rates are used for five reactions (Fig. 7). Their model neglects CETPmediated transport of TG out of HDL, and CE out of ApoBcontaining particles. However, their model somehow makes use of the ratio of TG to (CE + TG) associated to CETP. Comparing the strength of correlation between Hübner’s model and our model, ours appears to be superior with regard to HDLΔTG, LDLΔTG, and slightly better for VLDLΔTG (R = 0.515, R = 0.563, R = 0.544 vs our model’s: R = 0.683, R = 0.753, R = 0.562, respectively).
Comparison between the CETP model of Lu et al. [17], Hübner et al. [18] and ours presented. Lu’s model uses firstorder reactions; Hübner’s model uses 2 secondorder and 3 firstorder reactions. The model presented here assumes a steady state and uses one reaction rate for CETP. Due to CETP’s diffusionlike property, the amount of TG entering and leaving a lipoprotein depends on the surface and the TG and CE load of the lipoprotein and all other lipoproteins.
Comparing the Hübner and Lu models to ours presented here, note that ours is less complex, as it only requires one reaction rate instead of 5 or 3, respectively. Hence, unlike their models, which describe each flux from one fraction to another with their own reaction, our newly derived model represents a rather holistic model.
Concept of isotransfer point (ITP)
Figure 6a displays the medians of the ITP in the ‘normal’ and the ‘high TG’ group. Intuitively speaking, the point is lower in the ‘normal’group. Based on our model, the ITP ‘CETP_{TG}’ is the mean of all fraction’s TG/(TG + CE) ratios weighted by the corresponding surface. As clarified in Fig. 6 the greater the distance between the TG/(TG + CE) ratio of a lipoprotein to the ITP, the higher is the net flux of TG into/out of this fraction (if the ratio is smaller/greater than the ITP). Thus we expect LDL and HDL to have a higher TG turnover in the ‘high TG’ situation.
The ITP may be changed due to lipidlowering medication or during the postprandial state. This change may also exert a significant influence on HDL metabolism. The TG/(CE + TG) ratio in IDL in the ‘high TG’ group is surprisingly low. This might be caused by a prolonged retention time in plasma, which may be characteristic in the hypertriglyceridemic state. The clinical significance of the shift of the ITP from LDL to IDL in hyperlipidemia remains unclear. However, the atherogenity of certain lipoproteins (e.g. IDL) may be altered if its metabolic function (CETPmediated TG transfer) changes.
Data
Our in vitro data mirror the corresponding in vivo redistribution of lipids to a fairly good degree, as data from Liu et al. [19] as well as our group’s experiments have revealed evidence that the net change in TG in lipoprotein fractions by CETP action during plasma storage at 37 °C is linear for at least four hours.
Further data from experiment 2 suggest that CETP is the only factor responsible for TG redistribution among lipoprotein subfractions. CETP inhibition suggests that the esterification rate on nonHDL is very low or even nonexistent. However, the data may be too noisy to make precise claims.
In contrast to other methods measuring CETPactivity, the data we used to estimate the CETPmediated TGchange enables the estimation of TG in all fractions, namely VLDL, IDL, LDL and HDL, simultaneously. The data cover a broad physiological range of TG changes mediated by CETP, as individuals with healthy and various pathologic lipid profiles have been considered.
Very probably CETP activity differs among individuals. However, as the first step, our model estimates the CETPmediated relative redistribution of TG among the lipoproteins, which is independent of the CETP reaction rate, given only the TG, CE and ApoB (or ApoA1) concentration of VLDL, IDL, LDL and HDL. Based on this estimation, absolute fluxes may be inferred assuming a fixed reaction rate in a second step, if no additional information is given.
Model fit
Despite the noise of TG change data in lipoprotein fractions, the predicted TG fluxes correlate well with the observed corresponding changes (Fig. 2). While HDL and LDL fit well, VLDL and especially IDL appear more problematic. This could be mainly caused by relatively greater measurement imprecision, but also by a shift of ApoB from IDL to VLDL or viceversa during the one hour storage at 37 °C – as our group already observed [11]. Considering the R^{2} values between modelled and measured TG changes in VLDL, LDL and HDL (Fig. 2), mind that on the one hand our rather imprecise method of measurement may strongly lower the ‘true’ R^{2} value and on the other hand that our prediction is superior compared to the other models mentioned above, as well as to possible linear correlations mentioned in the results part (Table 2).
Equimolarity
Assuming an equimolar lipidexchange it is expected that the amount of TG entering/leaving a lipoprotein fraction equals the amount of CE leaving/entering it if LCAT inhibition is given. Considering ΔCE and ΔTG in experiment 2 (Fig. 1) we noted a significant increase in the TGmedian of 11.3 μmol/L and 19.2 μmol/L in LDL and HDL under LCAT inhibition. Surprisingly no corresponding loss of CE (0 and − 5.2 μmol/L, respectively) was observed. However, due to the small sample size and measurement imprecision, we cannot make a solid statement on this point.
Taken together, we identified some indicators suggesting that the CETP exchange may not necessarily be equimolar. However, due to experiment 2’s low case numbers, our data do not suffice to clarify this issue.
Lipid transfer inhibitor protein (LTIP)
There is strong evidence of a protein inhibiting CETP, namely lipid transfer inhibitor protein (LTIP) or Apolipoprotein F, which mainly interacts with LDL particles [20]. We investigated whether the resulting lower availability of LDL would have a positive effect on our model’s fit by reducing s(P_{LDL}), the surface of LDL. Potter uses a factor of 0.73 for LDL inhibition by LTIP. Such a reduction in our case leads in fact to a better correlation between predicted and observed changes in LDL and especially HDL: If s(P_{LDL})*0.4 is used instead of s(P_{LDL}), the correlation is optimised. If s(P_{LDL}) is reduced, the correlation coefficients R in HDL and LDL increase slightly, while the corresponding coefficient in VLDL decreases simultaneously. Further the estimated and observed proportion of TG enrichment between LDL and HDL tends to impair as s(P_{LDL}) is reduced. A possible inhibitory effect of ApoF is thus not implemented in our model.
CETP mass
Estimating TG fluxes in experiment 1, we assumed all r_{cetp} to be equal. This of course did not hold true. The parameter r_{cetp} rather depends on CETP’s mass concentration and the ability of the particular CETPphenotype to mediate lipid exchange.
In the literature, CETP’s plasma concentration does not differ that much depending on gender or pathologies [8, 21,22,23,24,25]. We are thus assuming that the mass is equal in experiment 1’s different subgroups. However, our model may not account for changes in r_{cetp} due to genetic variants of CETP.
Properties of CETP exchange
Considering our model and applying it to several lipidprofiles, two issues deserve attention: First, the driving force for CETP exchange is the mass and composition of triglyceriderich lipoproteins such as VLDL and chylomicrons, as they contain the most TG in plasma. Figure 4 illustrates this fact, as in the ‘high TG’ case, LDL and HDL both revealed the highest increase in TG mass. This issue might be an important factor, when assessing TG influence on lipid metabolism (for example the postprandial state).
Second, TG from VLDL is distributed to LDL and HDL to a similar degree in all four metabolic states investigated (Fig. 4). However, depending on their surface and lipid composition, the proportion of TG entering HDL compared to TG entering LDL may differ strongly, as both lipoproteinspecies compete against each other (consider for example the ‘high LDL’ case in Fig. 4).
Let us assume that TG’s high flux to LDL or HDL is atherogenic, as (in combination with hepatic lipase) it can lead to small, dense LDL and small HDL. Following the model of Lu et al. [17], lipolysis by hepatic lipase is the main pathway causing HDL particles to shrink. Small HDL are more likely to exit the plasma. Hence, via CETP hypertriglyceridemic plasma may lead to an atherogenic lipoprotein phenotype characterized by low HDL particle concentrations and small HDL and LDL particles. Our model might thus help to quantify the causal relationship between hypertriglyceridemia and those atherogenic characteristics. Furthermore, following our model’s dynamics, a high HDL particle concentration may prevent small dense LDL particles from forming by reducing the net flux of TG to LDL in such a hypertriglyceridemic state. However, the focus of this paper making is not claims about HDL functionality by relying on its lipidcomposition.
Note that only 43 normolipidemic persons out of experiment 1 were used for model calibration. Hence, the results for ‘high TG’, high LDL’ and ‘low HDL’ presented in Fig. 4 are extrapolated and demonstrate our model’s ability to predict CETP redistribution of TG based on lipid and apolipoprotein concentrations only.
When studying lipoprotein metabolism, the concepts considered are usually delipidation of ApoB containing lipoproteins and the reversed cholesterol transport mediated by HDL. However, both metabolic pathways are intimately connected by the CETPmediated flux of TG and CE. This paper is intended to give a new perspective on lipoprotein metabolism concerning CETPmediated fluxes of TG and CE.
Conclusion/outlook
Our newly devised model may be employed to estimate CETPflux data among different lipoprotein fractions using only CE, TG, ApoA1, and ApoB concentrations in lipoprotein fractions. Our model’s structure relies on the known biochemical characteristics of CETP, but is independent of its molecular transport mechanism. The means we have presented to estimate CETP action may help to deepen our understanding of CETP’s role in several physiological and pathological scenarios. For example, lipoprotein metabolism in the postprandial state remains difficult to access, as there is strong TG redistribution from chylomicrons to other lipoproteins. Hence, delayed catabolism of chylomicrons may exert strong effects on HDL and LDL metabolism. Thus, applying our approach to address CETP in such conditions it may now be possible to also make more accurate predictions even in the postprandial state.
Methods
Subjects & experimental procedure
For experiment 1 we used n = 91 plasma samples from individuals with a broad range of hyper, normo and hypolipidemic states. Plasma samples are broken down into 4 disjoint subgroups:

‘normal’: TG < 1.65 mmol/L, LDLCholesterol< 4.03 mmol/L, HDLCholesterol> 1.04 mmol/L for males and > 1.17 mmol/L for women;

‘high TG’: Total TG ≥ 1.65 mmol/L;

‘high LDL’: TG < 1.65 mmol/L, LDLCholesterol≥4.03 mmol/L;

‘low HDL’: TG < 1.65 mmol/L, LDLCholesterol< 4.03 mmol/L, HDLChol≤1.04 mmol/L for males and ≤ 1.17 mmol/L for women
For experiment 2 we used plasma samples of n = 11 volunteers (n = 9 normal, n = 2 nonnormal). The limits of subgroup ‘normal’ are based on guidelines from the German medical association ‘Deutsche Gesellschaft zur Bekämpfung von Fettstoffwechselstörungen und ihren Folgeerkrankungen DGFF (LipidLiga) e. V.’ (https://www.lipidliga.de/fueraertze/empfehlungen). We followed a method described by Jansen and colleagues [11]. In short, EDTAplasma was separated into two samples for both experiments. While the first was stored at 4 °C immediately, the second was exposed to 37 °C for one hour (water bath) and then stored at 4 °C as well. Subsequently, lipoprotein fractions were isolated via ultracentrifugation as described before [26]. Lipoproteins were separated into VLDL, IDL, LDL, HDL and lipid deficient serum (LDS) and further into three LDL subfractions LDLI, LDLII and LDLIII and three HDL subfractions HDL2b, HDL2a and HDL3. TG, CE, free cholesterol (FC), PL as well as ApoB and Apolipoprotein A1 (ApoA1) were measured in plasma and each separated lipoprotein fraction. Further the ratios ‘lipid to ApoB’ and ‘lipid to ApoA1’ were calculated to estimate the lipoprotein composition. These values were measured with storage at 37 °C or 4 °C, the corresponding differences in lipid and apolipoprotein mass and lipoprotein composition were calculated as the differences of baseline (4 °C) minus the 1 h 37 °C value and denoted with Δ.In experiment 2, each plasma sample was additionally divided into 4 groups in which only CETP, only LCAT, or none of both enzymes were inhibited. CETP was inhibited via Torcetrapib (1 μM, resolved in dimethyl sulfoxide), LCAT was inhibited by sodium iodacetate (5 mM resolved in Tris buffer [27]).
The coefficient of variation (CV: SD/mean*100) was determined by taking 20 repeated measurements of lipoprotein fractions.
The model
This work is based on the model by Potter et al. [1], which was modified by our group recently [11]. Briefly, metabolism of FC, CE and TG in LDL was modelled on the molecular level. In contrast to the recently published model, which assumed a onetoone exchange of TG and CE, this restriction for equimolarity is now relaxed, as only TG are considered for flux calculations. If CETP mediates an exchange between two lipoproteins (ternary model) or between a CETP loaded with CE or TG (shuttle model), neither CE nor TG is preferred. Further, in line with Morton [12, 13] we assume that the probability of the association of CETP to a lipoprotein is proportional to its surface. To simplify our model, the lipoprotein core consists of only TG and CE.
Let us consider the situation in which a CETP interacts with a lipoprotein particle. We are assuming that two independent events will occur:

Event 1: a TG molecule leaves the lipoprotein or no TG molecule leaves the lipoprotein

Event 2: a TG molecule enters the lipoprotein or no TG molecule enters the lipoprotein
Given those 2 events, a CETP interaction can lead to 4 possible outcomes:

A TG leaves the lipoprotein (heteroexchange)

A TG enters the lipoprotein (heteroexchange)

A TG enters and a TG leaves the lipoprotein (homoexchange)

No TG enters or leaves the lipoprotein
Considering event 1, the probability that a TG molecule leaves the particle equals the particle’s \( \frac{TG}{CE+ TG} \) v/v ratio and the complementary probability (no TG molecule leaves the particle) equals the particle’s \( \frac{CE}{CE+ TG} \) v/v ratio.
Considering event 2, the probability that a TG molecule enters the particle equals a parameter called CETP_{TG} while (1 CETP_{TG}) is the complementary probability (no TG molecule enters the particle).
The parameter CETP_{TG} is derived out of lipoprotein surfaces and TG and CE concentrations: CETP_{TG} mirrors the \( \frac{TG}{CE+ TG} \) v/v ratio of all TG and CE available from all lipoproteinparticles (weighted by their surface) in plasma to CETP. Let P_{i} denote the ith lipoprotein fraction (characterised by its mean TG and CE composition), where i ∈ I = {VLDL, IDL, LDL − I, LDL − II, LDL − III, HDL2b, HDL2a, HDL3}. For each fraction P_{i} respective its particlesurface S(P_{i}), its concentration of particles k(P_{i}), its \( \frac{TG}{CE+ TG} \) v/v ratio q(P_{i}) and its \( \frac{CE}{CE+ TG} \) v/v ratio q^{′}(P_{i}) = 1 − q(P_{i}) can be allocated. Consequently, the surface (corresponding to the given particle concentration) of all lipoproteins equals\( \sum \limits_{i\in I}s\left({P}_i\right)k\left({P}_i\right) \). It holds for CETP_{TG}, the \( \frac{TG}{CE+ TG} \) v/v ratio of all lipoproteinparticles weighted by their surface:
The methods applied for calculating a particle’s surface and the concentration of particles in the HDL subfractions are described previously [11].
Taken together, given an interaction of CETP with a lipoprotein of fraction P_{i}, the lipoprotein acquires a TG molecule in total, if no TG molecule leaves the particle and one TG molecule enters the particle. The corresponding joint probability is q^{′}(P_{i}) ∙ CETP_{TG}. The probability of a TG molecule leaving the lipoprotein in total can be derived analogously as q(P_{i}) ∙ (1 − CETP_{TG}). These exchanges are socalled heteroexchanges, as a TG molecule is replaced by a nondefined number (due to the assumption of nonequimolary exchange) of CE molecules. The corresponding probability of a TGhomoexchange (a TG leaves and a TG enters the lipoprotein) is q(P_{i}) ∙ CETP_{TG}.
Consequently, the fluxes of TG to the quadruple [VLDL, IDL, LDL, HDL] are proportional to
where r_{cetp} is the reaction rate of CETP. The corresponding fluxes of TG out of the quadruple [VLDL, IDL, LDL, HDL] are proportional to
The net flux of TG is proportional to the sum of term 1 and term 2:
With k, s and q calculated by measured data, only the reaction rate r_{cetp} needs to be inferred to calibrate our model. The proportion of TG mass entering and leaving the single fractions is already determined by the model’s structure (it depends on TG, CE and particle concentrations in the single fractions). We used the subgroup ‘normal’ (n = 43) to infer r_{cetp} by minimising the difference between measured and predicted TG net flux in VLDL, IDL, LDL and HDL (weighted by the corresponding surface). All other non‘normal’ cases were used as a validation group.
The Wilcoxon signed rank test was performed to test for differences between the change in CE and TG due to plasma storage at 37 °C for one hour. The Pearson’s correlation coefficient R was used to describe relationships between variables.
Statistics were analysed with IBM SPSS version 21.0 (IBM SPSS Statistics, IBM Corporation, Chicago, IL), all other computations were done using Scilab 5.4.0 (Scilab enterprises, Le Chesnay). Graphs were depicted by GraphPad Prism Software V7.00.
Abbreviations
 ApoA1:

Apolipoprotein A1
 ApoB:

Apolipoprotein B100
 CE:

Cholesteryl ester
 CETP:

Cholesterylester transfer protein
 FC:

Free cholesterol
 HL:

Hepatic lipase
 ITP:

Isotransfer point
 LPL:

Lipoprotein lipase
 PL:

Phospholipid
References
 1.
Potter LK, Sprecher DL, Walker MC, Tobin FL. Mechanism of inhibition defines CETP activity: a mathematical model for CETP in vitro. J Lipid Res. 2009;50(11):2222–34.
 2.
Barter P, Rye KA. Cholesteryl ester transfer protein: its role in plasma lipid transport. Clin Exp Pharmacol Physiol. 1994;21(9):663–72.
 3.
McPherson R, Marcel Y. Role of cholesteryl ester transfer protein in reverse cholesterol transport. Clin Cardiol. 1991;14(2 Suppl 1):I31–4.
 4.
Kosmas CE, DeJesus E, Rosario D, Vittorio TJ. CETP inhibition: past failures and future hopes. Clin Med Insights Cardiology. 2016;10:37–42.
 5.
Lauer ME, GraffMeyer A, Rufer AC, Maugeais C, von der Mark E, Matile H, D'Arcy B, Magg C, Ringler P, Muller SA, et al. Cholesteryl ester transfer between lipoproteins does not require a ternary tunnel complex with CETP. J Struct Biol. 2016;194(2):191–8.
 6.
Charles MA, Kane JP. New molecular insights into CETP structure and function: a review. J Lipid Res. 2012;53(8):1451–8.
 7.
Qiu X, Mistry A, Ammirati MJ, Chrunyk BA, Clark RW, Cong Y, Culp JS, Danley DE, Freeman TB, Geoghegan KF, et al. Crystal structure of cholesteryl ester transfer protein reveals a long tunnel and four bound lipid molecules. Nat Struct Mol Biol. 2007;14(2):106–13.
 8.
Marcel YL, McPherson R, Hogue M, Czarnecka H, Zawadzki Z, Weech PK, Whitlock ME, Tall AR, Milne RW. Distribution and concentration of cholesteryl ester transfer protein in plasma of normolipemic subjects. J Clin Invest. 1990;85(1):10–7.
 9.
Ritsch A, Auer B, Foger B, Schwarz S, Patsch JR. Polyclonal antibodybased immunoradiometric assay for quantification of cholesteryl ester transfer protein. J Lipid Res. 1993;34(4):673–9.
 10.
Lagrost L. Determination of the mass concentration and the activity of the plasma cholesteryl ester transfer protein (CETP). Methods Mol Biol. 1998;110:231–41.
 11.
Jansen M, Pfaffelhuber P, Hoffmann MM, Puetz G, Winkler K. In silico modeling of the dynamics of low density lipoprotein composition via a single plasma sample. J Lipid Res. 2016;57(5):882–93.
 12.
Morton RE. Binding of plasmaderived lipid transfer protein to lipoprotein substrates. The role of binding in the lipid transfer process. J Biol Chem. 1985;260(23):12593–9.
 13.
Morton RE, Greene DJ. The surface cholesteryl ester content of donor and acceptor particles regulates CETP: a liposomebased approach to assess the substrate properties of lipoproteins. J Lipid Res. 2003;44(7):1364–72.
 14.
Cobbold CA, Sherratt JA, Maxwell SR. Lipoprotein oxidation and its significance for atherosclerosis: a mathematical approach. Bull Math Biol. 2002;64(1):65–95.
 15.
Wattis JA, O'Malley B, Blackburn H, Pickersgill L, Panovska J, Byrne HM, Jackson KG. Mathematical model for low density lipoprotein (LDL) endocytosis by hepatocytes. Bull Math Biol. 2008;70(8):2303–33.
 16.
Watterson S, Guerriero ML, Blanc M, Mazein A, Loewe L, Robertson KA, Gibbs H, Shui G, Wenk MR, Hillston J, et al. A model of flux regulation in the cholesterol biosynthesis pathway: immune mediated graduated flux reduction versus statinlike led stepped flux reduction. Biochimie. 2013;95(3):613–21.
 17.
Lu J, Hubner K, Nanjee MN, Brinton EA, Mazer NA. An insilico model of lipoprotein metabolism and kinetics for the evaluation of targets and biomarkers in the reverse cholesterol transport pathway. PLoS Comput Biol. 2014;10(3):e1003509.
 18.
Hubner K, Schwager T, Winkler K, Reich JG, Holzhutter HG. Computational lipidology: predicting lipoprotein density profiles in human blood plasma. PLoS Comput Biol. 2008;4(5):e1000079.
 19.
Liu XQ, Bagdade JD. Neutral lipid mass transfer among lipoproteins in plasma from normolipidemic subjects is not an equimolar heteroexchange. J Lipid Res. 1995;36(12):2574–9.
 20.
Morton RE, Greene DJ. Regulation of lipid transfer between lipoproteins by an endogenous plasma protein: selective inhibition among lipoprotein classes. J Lipid Res. 1994;35(5):836–47.
 21.
Hogue JC, Lamarche B, Gaudet D, Lariviere M, Tremblay AJ, Bergeron J, Lemieux I, Despres JP, Gagne C, Couture P. Relationship between cholesteryl ester transfer protein and LDL heterogeneity in familial hypercholesterolemia. J Lipid Res. 2004;45(6):1077–83.
 22.
Foger B, Ritsch A, Doblinger A, Wessels H, Patsch JR. Relationship of plasma cholesteryl ester transfer protein to HDL cholesterol. Studies in normotriglyceridemia and moderate hypertriglyceridemia. Arterioscler Thromb Vasc Biol. 1996;16(12):1430–6.
 23.
Sandhofer A, Kaser S, Ritsch A, Laimer M, Engl J, Paulweber B, Patsch JR, Ebenbichler CF. Cholesteryl ester transfer protein in metabolic syndrome. Obesity (Silver Spring). 2006;14(5):812–8.
 24.
Caslake MJ, Stewart G, Day SP, Daly E, McTaggart F, Chapman MJ, Durrington P, Laggner P, Mackness M, Pears J, et al. Phenotypedependent and independent actions of rosuvastatin on atherogenic lipoprotein subfractions in hyperlipidaemia. Atherosclerosis. 2003;171(2):245–53.
 25.
Guerin M, Bruckert E, Dolphin PJ, Turpin G, Chapman MJ. Fenofibrate reduces plasma cholesteryl ester transfer from HDL to VLDL and normalizes the atherogenic, dense LDL profile in combined hyperlipidemia. Arterioscler Thromb Vasc Biol. 1996;16(6):763–72.
 26.
Baumstark MW, Kreutz W, Berg A, Frey I, Keul J. Structure of human lowdensity lipoprotein subfractions, determined by Xray smallangle scattering. Biochim Biophys Acta. 1990;1037(1):48–57.
 27.
Kostner GM, Krempler F, Dieplinger H, Zechner R, Teubl I, Sandhofer F. Altered metabolism of low density lipoprotein in humans after prolonged incubation in plasma. Clin Sci (Lond). 1985;68(4):411–8.
Acknowledgements
The authors are grateful to Brigitte König, Bärbel Krumrey, and Armin Buchwald for blood drawing and Nülüfer Ödünc for technical assistance.
Funding
Not applicable.
Availability of data and materials
All data generated or analysed during this study are included in this published article [and its Additional file 1].
Author information
Affiliations
Contributions
MJ: analysed the data, performed the modelling, wrote the manuscript with support from the other authors. GP & MH: discussed the results and contributed to the final manuscript. KW: planned and supervised the experiments, discussed the results and contributed to the final manuscript. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Martin Jansen.
Ethics declarations
Ethics approval and consent to participate
Data were taken from two not yet terminated studies and from volunteers and were approved by the Freiburg Ethics Committee. Written consent was given by each participant.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file
Additional file 1:
Lipoprotein data of all subjects. (CSV 35 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Jansen, M., Puetz, G., Hoffmann, M.M. et al. A mathematical model to estimate cholesterylester transfer protein (CETP) triglycerides flux in human plasma. BMC Syst Biol 13, 12 (2019) doi:10.1186/s129180190679x
Received
Accepted
Published
DOI
Keywords
 Enzymology/enzyme mechanism
 Lipid transfer proteins
 Lipoproteins/kinetics
 Lipid transport
 Lipoprotein metabolism
 Mathematical model