- 1Department of Pharmacy, China-Japan Friendship Hospital, Beijing, China
- 2School of Information, Beijing University of Chemical Technology, Beijing, China
Objective: Physical instability of lipid in parenteral nutrition (PN) poses significant clinical safety risks. As lipid stability is influenced by multiple complex factors and remains incompletely characterized, this study aimed to quantify the relative importance of stability determinants and to develop a machine learning (ML) model for predicting stability in individualized PN prescriptions.
Methods: A retrospective meta-analysis integrated experimental data from multi-laboratory studies. The ML framework employed transfer learning for cross-laboratory data harmonization and Synthetic Minority Over-sampling Technique (SMOTE) for class imbalance mitigation. Model performance was evaluated using the area under the receiver operating characteristic curve (AUC-ROC) and accuracy.
Results: The datasets comprised 17 stability-related features (electrolytes, macronutrients, and storage conditions) extracted from 1,518 samples representing 872 unique PN formulations across 19 studies (2000 and 2024). The XGBoost model achieved exceptional predictive performance (accuracy: 98.2%, AUC 0.968). SHAP-based feature importance analysis identified the concentrations of Amino and phosphate, storage time and lipid composition as key stability determinants.
Conclusion: This study establishes the first interpretable ML framework for predicting lipid emulsions stability in PN, resolving cross-laboratory data heterogeneity. We have provided a high-accuracy prediction tool for assessing lipid emulsion stability in PN, while the methodology demonstrates generalizability for stability studies of complex drugs and nutrients formulations.
Introduction
Parenteral nutrition (PN) is liquid nutrition that is delivered directly to the bloodstream of patients unable to absorb adequate nutrients through the digestive system. PN simultaneously supplies macronutrients (amino acids, dextrose, and lipids), which constitute the caloric and protein supply, and micronutrients (vitamins and trace elements), which complement the diet. As one of the most complex preparations in hospital, PN carries inherent risks of physicochemical interactions among its components. These interactions can compromise emulsion stability, potentially leading to patient safety hazards.
Commercially intravenous lipid emulsion is thermodynamically unstable heterogeneous dispersions two-phase system. Oil droplets, primarily composed of triglycerides, are stabilized by a surfactant layer of phospholipids derived from egg lecithin. The phosphatidic acid outside the lipid droplets is completely ionized to form as zeta potential between −40 and −50 mV, which realizes mutual electrostatic repulsion, thereby preventing aggregation and maintaining a uniform distribution of lipid droplets. The addition of cations (e.g., Na+, K+, Ca2+, Mg2+) can neutralize this negative surface charge, reducing the zeta potential, lead to droplet aggregation and subsequent phase separation within PN. The critical aggregation number (CAN) is a parameter that tries to predict the amount of cationic charge that would disturb lipid emulsion. pH significantly influences zeta potential and stability. A pH range of 6.0–9.0 (1) is optimal for lipid emulsion stability. Consequently, low-pH additives like dextrose solutions can compromise stability (2), while amino acid solutions can exert a stabilizing effect by modulating the PN pH (2). Lipid emulsions become unstable at pH values below 5.0 (2). The American Society for Parenteral and Enteral Nutrition (ASPEN) stated in 2014 that PN containing >4% amino acids, >10% dextrose, and >2% lipid emulsion remained stable at room temperature for up to 30 h or refrigerated (5 °C) for up to 9 days followed by 24 h at room temperature (2), However, the 2024 ASPEN update provides only select concentration examples without comprehensive stability criteria (3). Similarly, other guidelines list validated PN formulations (4), but generally lack specific recommendations on macronutrient and electrolyte concentration limits essential for predicting physicochemical stability.
In summary, lipid emulsion stability in PN is influenced by a multitude of interacting factors, including electrolyte concentrations, macronutrient composition, pH, storage temperature, and duration. Furthermore, the inherent variability of individualized PN prescriptions means that existing stability studies are often limited to specific formulations. Consequently, a generalizable method for evaluating and predicting lipid emulsion stability across diverse clinical PN remains lacking in practice. To address this gap, the present study conducted a systematic overview of current research on lipid emulsion stability in PN. We extracted key features from these studies, detailing the analytical methodologies employed, the stability variables measured, and the principal factors investigated. Ultimately, this work aims to establish a clinically applicable prediction model for lipid emulsion stability in PN utilizing Machine Learning (ML).
Materials and methods
Study design
An ML-powered meta-analysis framework incorporating transfer learning techniques was developed to integrate cross-laboratory evidence. Studies investigating lipid emulsion stability in PN were systematically identified, screened, and subjected to original data. This process established a multidimensional ML database. Domain adaptation techniques were applied to address heterogeneous data distributions across laboratories. Subsequently, ML algorithms were employed to construct a lipid emulsion stability prediction model, with performance validated through cross-validation.
Search strategy
A systematic literature search was conducted across PubMed, Embase, and Web of Science to identify relevant studies on PN stability. The search was restricted to English-language publications from 2000 to 2024. The following search syntax was employed: (“parenteral nutrition” OR “total parenteral nutrition” OR “total nutrition admixture” OR “all in one nutrition” OR “lipid emulsion”) AND (“stability/instability” OR “compatibility/incompatibility” OR “safety/unsafety”).
Inclusion and exclusion criteria
Studies were included if they met the following criteria:
1. Experimental investigations of PN stability relevant to clinical practice.
2. Lipid emulsion stability had been analyzed at least one of Lipid droplet evaluation criteria following USP Chapter 729 <Globule Size Distribution in Lipid Injections>
• PFAT5 (The percentage of weighted volume of fat residing in droplet of diameter >5 microns).
• MDD (Mean droplet diameter).
3. Utilization of validated analytical methods.
• Dynamic Light Scattering (DLS)/Photon Correlation Spectroscopy (PCS) or Laser Diffraction (LD) for MDD.
• Light Extinction/Single Particle Optical Sensing (LE/SPOS) for PFAT5.
Studies were excluded if they:
1. Assessed only non-physical stability aspects (e.g., chemical degradation, precipitation).
2. Evaluated formulations containing non-standard compounding agents.
3. Analyzed lipid emulsion stability outside PN.
4. Focused solely on drug-PN compatibility.
5. Were conference abstracts, reviews, book chapters, or lacked accessible full text/raw data.
Flow diagram shown in Figure 1, summarizes the search strategy following PRISMA (5).
Feature extraction
The features were extracted from the experimental studies reported in the literature including Environmental variables, PN composition characteristics and stability assessment results. Two researchers independently extracted the features data of experimental studies based on the following features definitions and formed the original dataset after cross-verification.
1. Environmental variables included:
• Study year: In which year the study was conducted.
• Temperature storage: Preparation and storage temperature.
• Temperature test: The temperature during the stability analysis.
• Times storage: The storage time of preparation.
• Times test: The time for conducting the stability analysis.
2. PN composition characteristics:
• Lipid composition: Long-Chain Triglycerides-LCTs, Medium Chain Triglycerides-MCTs or MCT/LCT blends in different proportions.
• Addition of micronutrients (within or without) included: Fat-soluble vitamins, Water-soluble Vitamins, Trace elements, Heparin.
• Concentration of electrolyte (mmol/L) included: Potassium; Sodium; Calcium; Magnesium; Phosphate; Chloride.
• CAN: The critical aggregation number is a predictive parameter for cationic charge induced destabilization (CAN = a + 64b + 729c). a, b and c are the sum of the concentrations (mmol/l) of mono-, di- and trivalent cations.
• Concentrations of Macronutrient (g/100 mL) included: Amino, Dextrose, Lipid.
3. Results of lipid emulsion stability assessment: The methods that were used for determination the stability of lipid emulsion was following the United States Pharmacopeia (USP), lipid globule size distribution be controlled within specified limits:
• MDD: Instability threshold is MDD ≥ 0.5 μm (analytical methods: DLS/PCS or LD).
• PFAT5: Instability threshold is PFAT5 ≥ 0.05% PFAT5 < 0.05% (analytical methods: LE/SPOS).
Data preprocessing
1. Missing value handling: An automated protocol classified missingness as: Low (<5%), Moderate (5%–15%), High (15%–30%), and Extreme (>30%). Missing data mechanisms were evaluated using Little’s Missing Completely at Random (MCAR) test (α = 0.05). For variables with 5%–20% missingness, Multiple Imputation by Chained Equations (MICE) was implemented (10 iterations generating 5 imputed datasets).
2. Outlier detection & processing: A tripartite detection framework was employed: Modified Z-score (threshold: |Z| > 3.5); Isolation Forest (contamination = 0.05, n_estimators = 200, max_samples = 256, random_tate = 42); Local Outlier Factor (n_neighbors = 20, contamination = 0.05). A specific constraints rule based on the formulation of PN solution was established, and the violation of the rule was considered as abnormal: CAN >150, Amino acids (>1%), Glucose (>2%), Lipid emulsion (>1%), MDD > 50 nm, PFAT5 < 0.001% or >4%.
3. Dataset partitioning and leakage prevention: A strict random split was applied to partition the dataset: Training set n = 1,062 (70%), and Test set n = 456 (30%). This approach ensured clear separation, effectively preventing leakage of information. Parameters for all preprocessing steps (including SMOTE resampling, missing value imputation and feature scaling) were all determined before random split.
4. Multicollinearity handing: We calculated the Variance Inflation Factor (VIF) for all variables. We employed a stepwise approach: first removing variables with VIF > 10, for any remaining variable pairs with a correlation coefficient >0.8, retaining the variable deemed to have greater physical significance based on domain knowledge.
Statistical analysis
Data were expressed as frequencies (percentages) for categorical variables and as Mean ± SD or median (IQR) based on the test of normality for continuous variables. The Kolmogorov–Smirnov test was applied to test normality. Differences between stability and instability groups were analyzed using Student’s t-test for continuous variables and Fisher’s exact test for categorical variables. A correlation heatmap was used to test potential multicollinearity between features. Statistical analyses and modeling were conducted using IBM SPSS Statistics (v25.0) and scikit-learn (v1.1) within Python 3.7 (Python Software Foundation; accessed January 1, 2025).
Model development
1. Algorithmic diversity: Tree-based models (random forest, decision tree, XGBoost, gradient boosting), Non-tree models (logistic regression, support vector machine, K-Nearest neighbors, Naïve Bayes) and Ensemble/Neural (Multilayer Perceptron and AdaBoost).
2. Class imbalance mitigation: To address significant class imbalance (13.1% unstable emulsions), systematic evaluation of hybrid resampling techniques was implemented including: SMOTE, ADASYN, random under-sampling. Domain-informed optimization of sampling ratios enhanced minority-class (unstable emulsions) recognition performance across validated configurations.
3. Hyperparameter optimization: Model hyperparameters were systematically optimized using Bayesian methods to ensure robustness and stability.
4. Evaluation protocol: Model performance was comprehensively assessed using five core metrics: accuracy, precision, recall, F1-score, AUC-ROC. To ensure clinical reliability and enable transparent assessment of misclassification risks, particular emphasis was placed on a detailed confusion matrix analysis and class-specific performance metrics—including precision, recall, and F1-score—for both stable and unstable lipid emulsion categories, thereby enabling multidimensional validation of the model’s predictive capability.
5. Ensemble strategy: Final model outputs were weighted soft-voting aggregation of probabilistic predictions to enhance cross-laboratory generalization capability.
Model interpretation
We employed SHAP values to quantify feature importance and interpret model predictions. Grounded in cooperative game theory’s Shapley values, SHAP provides consistent, theoretically-grounded attribution of feature contributions to individual predictions. Predictions. Model-specific implementations were utilized: Trees SHAP for tree models and Kernels HAP for non-tree models. Additivity was verified (tolerance <0.01), with ensemble SHAP values computed as weighted averages of base models. Global importance was ranked by mean absolute SHAP values, identifying top features.
Results
Data characteristics
The final analytical dataset comprised 1,518 PN samples representing 872 unique clinical prescriptions, extracted from 19 studies published between 2000 and 2024. These studies involve six different types of emulsions (three kinds of LCTs, two kinds of MCT /LCTs and one kind of SMOF lipid), storage temperature from 4 to 37 °C and storage time from 0 h to 14 days, the concentration of monovalent cation (Na+/K+) from 0 to 230 mmol/L, divalent cations (Ca2+/Mg2+) from 0 to 92.5 mmol/L. The concentration of macronutrient (mean ± SD): Amino (3.4 ± 2.0 g/100 mL), Dextrose (11.2 ± 6.2 g/100 mL), Lipid (2.8 ± 1.0 g/100 mL), Complete study characteristics are detailed in Table 1.
Correlation between features and lipid emulsion stability in PN
Correlational analyses between extracted features and emulsion stability are presented in the correlation heatmap (Supplementary Figure S1). MDD analysis performed on 1,058 bags (69.7%), with instability (MDD ≥ 0.5 μm) detected in 32 samples (3.0% 32/1058). PFAT5 were determined 971 samples (64.0%), instability (PFAT5 ≥ 0.05%) observed in 173 samples (17.8% 173/971). Overall instability defined as violation of either USP threshold (MDD ≥ 0.5 μm or PFAT5 ≥ 0.05%), occurring in 199 samples (13.1% 199/1518). Comparative analysis between stable (n = 1,319) and unstable (n = 199) groups demonstrated significant differences (p < 0.05) in most variables except the concentration of Calcium, Magnesium, Lipid, Lipid composition and Trace elements. Complete statistical comparisons are detailed in Table 2.
Machine learning prediction of PN lipid stability
We compared the performance of ML models to predict PN lipid stability using ROC analysis. The AUC-ROC of the XGBoost model had the highest predictive value (AUC-ROC: 0.968, Accuracy: 0.980), followed by a random forest model (AUC-ROC: 0.962) and Gradient Boosting (AUC-ROC: 0.961). Detailed performance measures for all models are provided in Supplementary Table S1. We demonstrated the ROC curve comparison of 10 prediction models (Figure 2). Imbalance treatment optimization (Supplementary Table S2) confirmed SMOTE significantly enhanced XGBoost performance versus alternative methods. Finally, the XGBoost model achieved a recall of 0.98, a Precision of 0.98 and F1 score of 0.98.
Detailed mechanism of XGBoost: Algorithm Principle, An ensemble learning method based on gradient-boosted decision trees, which builds a strong classifier by sequentially training multiple weak learners; Key Parameters, max_depth = 6 (controls tree depth to balance bias and variance), learning_rate = 0.1, n_estimators = 100 (number of trees); Objective Function, Obj = Σ l (yᵢ, ŷᵢ) + Σ Ω(ƒₖ), where “l” is the loss function and “Ω” is the regularization term; Feature Interaction, Capable of automatically learning non-linear interaction relationships between features without manual feature engineering.
The test set confusion matrix shows: True Positives (correctly identified unstable cases) = 55, False Negatives (unstable cases misclassified as stable) = 4. True Negatives (correctly identified stable cases) = 392, False Positives (stable cases misclassified as unstable) = 5. This results in a Sensitivity/Recall of 93.2% (55/59) and a Specificity of 98.7% (392/397). We demonstrated the confusion matrix of SMOTE+XGBoost model (Supplementary Figure S2).
Feature importance analysis
To provide recommendations for practice with a straightforward understanding of key features related to predicting lipid stability, SHAP analysis of the optimal XGBoost model identified key determinants of lipid emulsion stability with or without Multicollinearity handing (Figure 3). Amino, Phosphate, Dextrose, Time storage and Lipid composition demonstrated predominant predictive influence.
Figure 3. Features SHAP analysis in the XGBoost model. (A) No handling of multicollinearity; (B) Deal with multicollinearity.
Discussion
Core findings and model innovation
This study established an XGBoost-driven multimodal data fusion framework to achieve high-accuracy prediction of lipid emulsion physical stability in PN. The model demonstrated exceptional performance on the test set, significantly outperforming conventional logistic regression and random forest. Key innovations include intelligent multisource data integration through transfer learning and domain adaptation, resolving heterogeneity from environmental conditions and prescription differences by enabling equivalent data transformation across 19 studies. Optimized Feature Engineering, 17 stability relevant features were extracted.
To address the potential impact of multicollinearity among input variables (as shown in Supplementary Figure S1), we systematically evaluated and managed feature correlations during data preprocessing. Variables were carefully reviewed, and those with lower physicochemical relevance (e.g., CAN) were excluded. A comparative sensitivity analysis of SHAP values was conducted before and after this multicollinearity handling. Although the ranking of specific features saw minor shifts, the overall importance trend remained largely consistent. This process refined the contribution weight of individual variables to the model, thereby enhancing the robustness and interpretability of our final model.
Although multiple continuous variables showed statistically significant differences between the stable and unstable PN groups (Table 2), their inherent variability and intercorrelations make it difficult to determine the stability of lipid in individualized TNA. This further highlights the advantage of using a ML approach, which can integrate multidimensional interactions among features—such as amino acids, dextrose, and electrolytes—to produce a robust predictive model. Importantly, the SHAP analysis reaffirmed the model’s ability to synthesize these interrelated influences, with key factors (e.g., amino acids, dextrose, and storage time) consistently emerging as major contributors, thereby validating the model’s interpretability and clinical relevance.
The performance gain achieved through the SMOTE-based resampling strategy can be attributed to the intrinsic characteristics of our dataset. SMOTE-generated synthetic samples effectively populate sparse regions within the feature space. Consequently, SMOTE’s interpolation mechanism effectively populated these sparse regions with plausible synthetic samples, which significantly enhanced the model’s ability to identify the minority class. This approach elevated the AUC from 0.95 to 0.97 and elevated the Accuracy from 0.97 to 0.98, thereby ensuring the predictive safety of the model. We add a sensitivity analysis regarding key SMOTE parameters in Supplementary Table S3.
Analysis of key features and mechanisms
The stability of lipid in TNA is most significantly influenced by the concentration of amino acids, which play a crucial role in regulating the pH of TNA. This study first identifies phosphate concentration as a novel predictor of PN lipid stability through ML modeling—a previously unreported association whose mechanism warrants further investigation, though we hypothesize involvement of phosphate groups in the emulsion’s hydrophilic layer. Following features such as concentration of others electrolyte and dextrose, type of lipid, CAN, or storage time that have been largely evaluated and confirmed as relevant predictors of lipid stability in PN. CAN was excluded due to multicollinearity handling. This was also discovered in all models SHAP analysis (Supplementary Figure S3). Cations can disrupt the electrostatic repulsion between lipid droplets, thereby compromising the stability of lipid emulsions, though modulated by compositional factors (6). The 100-fold higher aqueous solubility of medium-chain caprylic acid (C8) versus palmitic acid (C16) (7), reduces oil-water interfacial tension, explaining the low stability of LCT emulsions (8, 9) observed in our model. While temperature and storage duration influence stability, their limited predictive role may relate to anomalous PFAT5 reduction over time—potentially attributable to compounding-induced air bubbles (10). Our stability assessment aligned with USP and Chinese Pharmacopoeia standards using MDD/PFAT5, as enlarged droplets (>5 μm) pose embolism risks: murine studies demonstrate PFAT5 > 0.05% causes microvascular damage and reticuloendothelial injury (11), consistent with clinical reports of pulmonary fat deposition in neonates (12); critically, PFAT5 > 0.4% induces phase separation visible as surface oil droplets (1), representing immediate life-threatening risks.
Value in clinical practice
As shown above, there has a critical impact of PN stability on patient safety, particularly the complex physicochemical behavior of lipid emulsions influenced by multifactorial interactions. In clinical practice, there are complex decisions in medical order, prescription review, PN preparation and PN safety monitoring. This ML-driven solution empowers pharmacists and physicians to rapidly evaluate PN compatibility by analyzing the influencing factors of specific prescriptions. ML has been increasingly used in the field of nutrition science, considering the diversity and complexity of individuals, ML can be used to predict precision nutrition support through multivariate and high-dimensional data analysis. In 2022–2027, the National Institutes of Health of the United States will provide $170 million to fund the development of algorithms for precision nutrition support (13). Our work provides a distinct contribution. Rather than merely applying XGBoost model, this study establishes a comprehensive, ML-enabled meta-analysis framework tailored specifically to PN stability. Our approach integrates a uniquely standardized clinical dataset, domain-specific feature engineering that combines environmental, electrolyte, and macronutrient variables, and an interpretable predictive tool with direct clinical utility. By adapting ML techniques to this specific challenge, we offer not only a practical system for stability prediction but also a reusable methodological framework for studying complex systems stability in pharmaceutical and nutrition sciences.
Limitations and strengths
As a retrospective investigation, this study necessitates further validation through prospective multi-laboratory trials, particularly to verify phosphate’s mechanistic role in emulsion destabilization. In addition, lipid emulsion stability parameters were insufficient: pH, Zeta potential, peroxide parameter, and free fatty acid were lacking, but we followed scientific physical stability verification criteria. Subsequently, extended multi-task learning was used to predict calcium and phosphorus precipitation, vitamin degradation and other stability risks.
ML approach overcomes the clinical impracticality of experimental verification for individualized PN formulations by modeling multifactorial interactions across heterogeneous datasets. Future research will expand training data to enhance generalizability and integrate physical–chemical stability predictions, ultimately constructing a comprehensive risk assessment paradigm for PN quality assurance.
Conclusion
This study establishes an XGBoost-driven multimodal fusion framework that transforms opaque model predictions into clinically actionable insights via SHAP interpretability. By resolving cross-laboratory data integration barriers in predicting lipid emulsion stability, we deliver a high-accuracy, explainable decision tool for clinical practice. The synergistic effect of key role of Phosphate, Dextrose, Sodium and Lipid composition (MCT/LCT) provide evidence-based reference for individualized PN formulation review. This study delivers both a high-accuracy predictive model for lipid emulsion stability in PN and a pioneering research strategy for assessing stability in complex pharmaceutical formulations.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
SY-g: Methodology, Writing – review & editing, Writing – original draft, Investigation, Data curation, Project administration. WX-l: Data curation, Investigation, Writing – review & editing, Writing – original draft, Formal analysis. CY: Supervision, Visualization, Software, Writing – review & editing, Formal analysis, Project administration, Writing – original draft, Data curation. QW-j: Supervision, Investigation, Methodology, Writing – review & editing, Formal analysis, Funding acquisition, Writing – original draft, Resources. LP-m: Project administration, Writing – review & editing, Conceptualization, Visualization, Supervision. ZL: Funding acquisition, Writing – original draft, Visualization, Resources, Validation, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by National Science and Technology Major Projects of China (2024ZD0522305), Beijing Health Technology Promotion Project (BHTPP2024105).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2025.1668464/full#supplementary-material
Abbreviations
PN, Parenteral nutrition; MDD, Mean droplet diameter; PFAT5, The percentage of weighted volume of fat residing in droplet of diameter >5 microns; CAN, Critical aggregation number; ML, Machine learning; SHAP, Shapley additive explanations; AUC-ROC, Area under the receiver operating characteristic curve; SMOTE, Synthetic Minority Over-sampling Technique.
References
1. Driscoll, DF. Physicochemical assessment of total nutrient admixture stability and safety: quantifying the risk. Nutrition. (1997) 13:166–7. doi: 10.1016/s0899-9007(96)00394-2
2. Boullata, JI, Gilbert, K, Sacks, G, Labossiere, RJ, Crill, C, Goday, P, et al. A.S.P.E.N clinical guidelines parenteral nutrition ordering, order review, compounding, labeling, and dispensing. JPEN J Parenter Enteral Nutr. (2014) 38:334–77. doi: 10.1177/0148607114521833
3. Boullata, JI, Salman, G, Mirtallo, JM, Canada, TW, Monczka, J, Gura, KM, et al. Parenteral nutrition compatibility and stability: practical considerations. Nutr Clin Pract. (2024) 39:1150–63. doi: 10.1002/ncp.11189
4. Bolisetty, S, Osborn, D, Schindler, T, Sinn, J, Deshpande, G, Wong, CS, et al. Standardised neonatal parenteral nutrition formulations - Australasian neonatal parenteral nutrition consensus update 2017. BMC Pediatr. (2020) 20:59. doi: 10.1186/s12887-020-1958-9
5. Page, MJ, and Moher, D. Evaluations of the uptake and impact of the preferred reporting items for systematic reviews and meta-analyses (PRISMA) statement and extensions: a scoping review. Syst Rev. (2017) 6:263. doi: 10.1186/s13643-017-0663-8
6. Lobo, BW, da Veiga, VF, Cabral, LM, Michel, RC, Volpato, NM, and de Sousa, VP. Influence of the relative composition of trace elements and vitamins in physicochemical stability of total parenteral nutrition formulations for neonatal use. Nutr J. (2012) 11:26. doi: 10.1186/1475-2891-11-26
7. Bach, AC, and Babayan, VK. Medium-chain triglycerides: an update. Am J Clin Nutr. (1982) 36:950–62. doi: 10.1093/ajcn/36.5.950
8. Driscoll, DF. Pharmaceutical and clinical aspects of lipid injectable emulsions. JPEN J Parenter Enteral Nutr. (2017) 41:125–34. doi: 10.1177/0148607116673187
9. Driscoll, DF, Parikh, M, Silvestri, AP, Klütsch, K, Bistrian, BR, and Nehne, J. Establishing a stability window for medium- and long-chain-triglyceride lipid-based total nutrient admixtures using USP standards. Am J Health Syst Pharm. (2006) 63:2135–43. doi: 10.2146/ajhp050520
10. Driscoll, DF, Silvestri, AP, Nehne, J, Klütsch, K, Bistrian, BR, and Niemann, W. Physicochemical stability of highly concentrated total nutrient admixtures for fluid-restricted patients. Am J Health Syst Pharm. (2006) 63:79–85. doi: 10.2146/ajhp050122
11. Driscoll, DF, Ling, PR, Andersson, C, and Bistrian, BR. Hepatic indicators of oxidative stress and tissue damage accompanied by systemic inflammation in rats following a 24-hour infusion of an unstable lipid emulsion admixture. JPEN J Parenter Enteral Nutr. (2009) 33:327–35. doi: 10.1177/0148607108327155
12. Klang, MG. PFAT5 and the evolution of lipid admixture stability. JPEN J Parenter Enteral Nutr. (2015) 39:67s–71s. doi: 10.1177/0148607115595976
13. Kirk, D, Kok, E, Tufano, M, Tekinerdogan, B, Feskens, EJ, and Camps, G. Machine learning in nutrition research. Adv Nutr. (2022) 13:2573–89. doi: 10.1016/j.advnut.2023.03.001
14. Driscoll, DF, Bacon, MN, and Bistrian, BR. Physicochemical stability of two types of intravenous lipid emulsion as total nutrient admixtures. JPEN J Parenter Enteral Nutr. (2000) 24:15–22. doi: 10.1177/014860710002400115
15. Driscoll, DF, Giampietro, K, Wichelhaus, DP, Peterss, H, Nehne, J, Niemann, W, et al. Physicochemical stability assessments of lipid emulsions of varying oil composition. Clin Nutr. (2001) 20:151–7. doi: 10.1054/clnu.2001.0375
16. Driscoll, DF, Nehne, J, Peterss, H, Klütsch, K, Bistrian, BR, and Niemann, W. Physicochemical stability of intravenous lipid emulsions as all-in-one admixtures intended for the very young. Clin Nutr. (2003) 22:489–95. doi: 10.1016/s0261-5614(03)00046-3
17. Gonyon, T, Carter, PW, Dahlem, O, Denet, AR, Owen, H, and Trouilly, JL. Container effects on the physicochemical properties of parenteral lipid emulsions. Nutrition. (2008) 24:1182–8. doi: 10.1016/j.nut.2008.06.003
18. Skouroliakou, M, Matthaiou, C, Chiou, A, Panagiotakos, D, Gounaris, A, Nunn, T, et al. Physicochemical stability of parenteral nutrition supplied as all-in-one for neonates. JPEN J Parenter Enteral Nutr. (2008) 32:201–9. doi: 10.1177/0148607108314768
19. Chaieb, DS, Chaumeil, JC, Jebnoun, S, Khrouf, N, Hedhili, A, and Sfar, S. Effect of the intravenous lipid emulsions on the availability of calcium when using organic phosphate in TPN admixtures. Pharm Res. (2008) 25:2545–54. doi: 10.1007/s11095-008-9671-7
20. Driscoll, DF, Silvestri, AP, and Bistrian, BR. Stability of MCT/LCT-based total nutrient admixtures for neonatal use over 30 hours at room temperature: applying pharmacopeial standards. JPEN J Parenter Enteral Nutr. (2010) 34:305–12. doi: 10.1177/0148607109338215
21. Skouroliakou, M, Kountouri, AM, Hatziantoniou, S, Koutri, K, and Chiou, A. Physicochemical stability assessment of all-in-one parenteral emulsion for neonates containing SMOFlipid. Eur J Hosp Pharm. (2012) 19:514–8. doi: 10.1136/ejhpharm-2012-000121
22. Athanasiou, C, Hatziantoniou, S, Skouroliakou, M, and Markantonis‐Kyroudis, S. Assessment of the physicochemical stability of all-in-one parenteral emulsions for neonates according to USP specifications. JPEN J Parenter Enteral Nutr. (2014) 38:867–72. doi: 10.1177/0148607113499589
23. De Cloet, J, Van Biervliet, S, and Van Winckel, M. Physicochemical stable standard all-in-one parenteral nutrition admixtures for infants and children in accordance with the ESPGHAN/ESPEN guidelines. Nutrition. (2018) 49:41–7. doi: 10.1016/j.nut.2017.11.019
24. Sayed, OAE-A, Hassan, SB, Abdelkader, A, Elsabahy, M, Abdelaziz, NHR, and El-Sayed, AM. Stability study and clinical evaluation of lipid injectable emulsion in parenteral nutrition admixtures used for preterm neonates. Nutr Clin Pract. (2021) 36:696–703. doi: 10.1002/ncp.10556
25. Zhao, B, Gao, R, Jiao, L, Zhang, F, Wang, B, and Mei, D. Physical stability of medium-chain triglyceride/long-chain triglyceride emulsion injections from 5 manufacturers in high-concentration electrolyte-based Total nutrient admixtures. JPEN J Parenter Enteral Nutr. (2021) 45:1275–84. doi: 10.1002/jpen.2009
26. Gostynska, A, Starkowska, J, Sobierajska, P, Jelinska, A, and Stawny, M. All-in-one pediatric parenteral nutrition admixtures with an extended shelf life-insight in correlations between composition and physicochemical parameters. Pharmaceutics. (2021) 13:1017. doi: 10.3390/pharmaceutics13071017
27. Gao, S, Fan, J, Jia, X, Li, B, Li, X, Li, Z, et al. PFAT5 stability assessment of Lipovenoes MCT in total nutrient admixtures. Ann Palliat Med. (2021) 10:12244–50. doi: 10.21037/apm-21-3353
28. Otero-Millán, L, Bea-Mascato, B, Legido Soto, JL, Martínez-López-De-Castro, N, and Lago Rivero, N. Electron microscopy for the stability assessment of parenteral nutrition admixtures: focus on precipitation. Nutrients. (2024) 16:1390. doi: 10.3390/nu16091390
Keywords: lipid emulsion stability, parenteral nutrition, XGBoost, SHAP interpretation, machine learning
Citation: Yong-guang S, Xue-lian W, Yong C, Wang-jun Q, Peng-mei L and Lei Z (2025) Machine learning predicts lipid emulsion stability in parenteral nutrition using multi-laboratory literature data. Front. Nutr. 12:1668464. doi: 10.3389/fnut.2025.1668464
Edited by:
Sudeep Tanwar, Nirma University, IndiaReviewed by:
Zhongyuan Zhao, Binghamton University, United StatesVladimir Kondratenko, All-Russian Dairy Research Institute, Russia
Copyright © 2025 Yong-guang, Xue-lian, Yong, Wang-jun, Peng-mei and Lei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhang Lei, emVyZ3psQDE2My5jb20=; Qin Wang-jun, cXdqMjAwNHdhbmdAMTYzLmNvbQ==
Wang Xue-lian1