Explainable machine learning model of coronary artery disease combined with diabetes: development and validation study

Abstract

Background:

Coronary artery disease (CAD) demonstrates a strong bidirectional association with diabetes mellitus, which not only elevates cardiovascular disease risk but also correlates with poorer clinical prognosis. Prognostication in patients with comorbid CAD and diabetes remains a critical clinical challenge, significantly influencing therapeutic decision-making. Leveraging readily available clinical parameters for predicting adverse outcomes in this population offers substantial clinical value. This investigation employs machine learning algorithms to develop predictive models for prognostic assessment in CAD patients with diabetes comorbidity.

Method:

We conducted a retrospective cohort study of 389 patients with comorbid coronary artery disease (CAD) and diabetes mellitus. The cohort was randomly allocated into a training set (n = 273) and an internal validation set (n = 116). Feature selection utilized LASSO regression followed by backward stepwise Cox regression analysis. A nomogram incorporating independent predictors was developed for clinical application. Model performance was assessed through discrimination metrics, calibration plots, and decision curve analysis (DCA). Random survival forest analysis validated the clinical significance of selected variables.

Result:

Our modeling approach employed a systematic methodology: LASSO regression for initial feature selection followed by backward stepwise Cox regression analysis, yielding eight independent predictors.The final model incorporated hemoglobin, INR, albumin, NT-proBNP, age, fibrinogen, diuretic use, and digitalis therapy. The integrated model demonstrated strong discriminative performance for mortality prediction across both training (AUC = 0.846, 0.838, 0.82) and validation cohorts (AUC = 0.824, 0.813, 0.798) at 3-, 5-, and 8-year intervals. Calibration plots and decision curve analysis confirmed model reliability and clinical utility over time. A nomogram was developed to facilitate individualized risk stratification. Kaplan–Meier analysis showed significant survival stratification by individual predictors, and restricted cubic spline analysis identified non-linear associations between continuous variables and mortality. Random survival forest analysis prioritized five key predictors (hemoglobin, INR, albumin, NT-proBNP, age). Comparative evaluation against the 9-variable model confirmed superior performance of the comprehensive model across all timepoints.

Conclusion:

Our multimodal prognostic model demonstrated robust performance in predicting all-cause mortality among patients with coronary artery disease and diabetes comorbidity. The nomogram's capacity for personalized risk estimation offers potential utility in clinical decision-making and patient stratification.

Introduction

The primary cause of death worldwide and a significant contributor to disability is coronary artery disease (CAD) (1). Diabetic patients are more likely to experience severe acute coronary syndromes and have worse outcomes than non-diabetic patients (2, 3). Although microvascular complications substantially impact diabetes-related morbidity, mortality primarily stems from macrovascular manifestations including CAD and cerebrovascular events (4). This clinical landscape underscores the critical need for accurate risk stratification tools to inform therapeutic strategies and healthcare resource utilization.

Conventional prognostic models have predominantly utilized isolated parameters including natriuretic peptide levels, glycated hemoglobin A1c (HbA1c), and New York Heart Association (NYHA) functional classification (5, 6). Established risk assessment tools include the Framingham Risk Score (FRS) in North America and the Systematic Coronary Risk Evaluation (SCORE) model in Europe (7, 8). However, these models demonstrate limited discriminative capacity for CAD mortality prediction and insufficient integration of modern prognostic factors for acute coronary syndrome (9, 10). Furthermore, FRS applicability shows significant variability across ethnic and geographic populations. Most current prognostic models focus on identifying high-risk subgroups rather than providing comprehensive risk stratification. Emerging biomarkers like the triglyceride-glucose (TyG) index show promise for predicting major adverse cardiovascular events (MACE) in post-CABG patients with CAD and diabetes, potentially enabling targeted monitoring and preventive interventions (11). Notably, validated clinical scoring systems for long-term prognosis assessment in comorbid CAD and diabetes populations remain largely unavailable.

The rapid advancement of artificial intelligence has enabled increasing machine learning applications in clinical prediction. Supervised learning algorithms enhance prognostic accuracy by iteratively optimizing model parameters against clinical outcome data (12). This technological evolution creates new opportunities for developing precision prognostic algorithms in diabetes-associated CAD management.

Our multimodal approach comprehensively integrates demographic profiles, clinical histories, pharmacological data, and echocardiographic parameters to enhance risk assessment in comorbid CAD-diabetes patients. We hypothesize this multidimensional strategy will surpass conventional single-parameter models in prognostic precision, ultimately improving clinical risk categorization and therapeutic decision-making.

Methods

Study design and participants

Out of 834 patients with CAD and diabetes, 389 met the inclusion and exclusion criteria and were included in the final analysis. This study analyzed clinical data from 389 CAD and diabetes patients treated between December 2014 and December 2024. The primary outcome was cardiovascular death. The study included 389 participants who satisfied the following inclusion criteria and had both diabetes and CAD: (1) aged over 18 years old; (2) diagnosed with CAD according to the European Society of Cardiology (ESC) diagnosis and treatment guidelines (13): a pathological process characterized by atherosclerotic plaque accumulation in the epicardial arteries, whether obstructive or non-obstructive; (3) fasting plasma glucose (FPG) ≥ 7 mmol/L (126 mg/dL) is introduced as the criterion for diagnosis of diabetes (14). Exclusion criteria consisted of (1) Acute or chronic infections are present concurrently; (2) the patient has a tumor, autoimmune illness, or hematological disease; (3) and the patient's information is lacking. Survival outcomes were verified through structured telephone interviews or clinical follow-up documentation.This study adhered to the Declaration of Helsinki and was approved by the Ethics Committee of the Second Hospital of Lanzhou University (2024A-1270), with a waiver of informed consent.

Data acquisition

At baseline, all participants underwent comprehensive demographic evaluation including sex, age, and clinical histories (hypertension, diabetes mellitus, atrial fibrillation, coronary heart disease, anemia, smoking status, and alcohol use). Medication profiles were documented, encompassing antiplatelet agents, anticoagulants, statins, nitrates, diuretics, and cardiac remodeling therapies. Standardized laboratory assessments included fasting glucose, hemoglobin, lipid profile (total cholesterol, triglycerides, HDL, LDL), renal function (serum creatinine), hepatic enzymes (AST, ALT), albumin, INR (International Normalized Ratio), CRP, and NT-proBNP (N-terminal pro B-type Natriuretic Peptide). Supplementary cardiovascular evaluations comprised 12-lead ECG and echocardiographic parameters.

In the feature selection process, we used LASSO regression to identify significant predictors of patient outcomes. We included time-dependent covariates and interaction terms in the LASSO regression. Conditions like diabetes and hypertension develop over time due to long-term factors. To account for this, we incorporated time-dependent variables to reduce confounding from the chronic nature of these diseases. Additionally, for medication use, we classified patients as “users” if they had been regularly prescribed medications for coronary artery disease and diabetes for six months or more, defining “non-users” as those who did not meet this threshold. This approach allowed us to effectively capture the temporal relationships and reduce potential biases in our analysis. Although time-dependent covariates were considered during variable selection, the final prognostic model was a baseline Cox proportional hazards model. Medication exposure variables (e.g., diuretics and digoxin) were constructed using prescription records from the first six months after the index date; patients who continuously received these medications for at least six months were classified as “users”, whereas those who did not meet this threshold were classified as “non-users”. Consequently, a formal time-dependent Cox model was not fitted and only baseline covariates were used in the final analysis.

The final cohort comprised 389 patients with confirmed CAD and diabetes. The training set was defined as 116 patients, and their early morning fasting peripheral blood samples were retained for the test of laboratory examination results. The training set consisted of 273 patients, as shown in Figure 1.

Figure 1

DStatistical analysis and model development

All statistical analyses and model construction were done using R software (version 4.3.1). The survival and rms packages were used for Cox proportional hazards modelling, random ForestSRC was used to construct the random survival forest, time-ROC and risk Regression were used to estimate time-dependent ROC curves and AUCs, rms and pec were used to generate calibration curves, and dcurves and rmda were used for decision-curve analysis. A summary of functions and packages is provided in Supplementary Table S7. Group comparisons employed χ² tests for categorical variables and Mann–Whitney U-tests for nonparametric continuous data. Statistical significance was determined using two-tailed α = 0.05. for the handling of missing data and continuous predictors, continuous variables with a missing rate of less than 5% were imputed using the random forest method. Feature selection followed a two-step process: LASSO regression with 10-fold cross-validation was first applied to the training cohort to identify potential predictors, followed by backward stepwise Cox regression. The variance inflation factor (VIF) was used to assess multicollinearity among selected variable. Several methods were used to assess the model's performance. R software (version 4.3.1; R Foundation for Statistical Computing, Vienna, Austria) was used for all statistical analyses. Categorical variables were compared using χ² tests, and continuous variables with non-normal distribution were compared using the Mann–Whitney U-test. A two-tailed alpha of 0.05 was considered statistically significant. Feature selection followed a two-step process: first, LASSO regression with 10-fold cross-validation was applied to identify predictors. These variables chosen by LASSO were then subjected to backward stepwise Cox regression analysis. The variance inflation factor (VIF) was used to evaluate multicollinearity among the final variables that were chosen. Several methods were used to assess the model's performance. Time-dependent receiver operating characteristic (ROC) curves for 3-, 5-, and 8-year mortality forecasts in both training and validation sets were used to evaluate discrimination. These time horizons were selected to reflect short−, intermediate− and long−term follow−up. The median follow−up (survival) time in our cohort was 4.68 years (interquartile range 0.55–8.71 years), so 5 years approximated the mid-point of follow−up and 8 years corresponded to the upper quartile. As shown in Supplementary Figure S6, follow−up durations ranged from less than 1 year to more than 10 years, with a small peak of patients completing long−term follow−up around 8 years. At 5 years, 184 participants were still under observation, and 111 participants remained at risk at 8 years (Supplementary Figure S6). Evaluating the model at 3, 5 and 8 years therefore captures short−, medium− and longer−term risk while maintaining adequate numbers at risk for reliable estimates. To evaluate the degree of agreement between expected and observed results, calibration curves were plotted. Decision curve analysis (DCA) was used to assess the model's clinical utility. Using Kaplan–Meier survival analyses, the prognostic value of categorical factors was evaluated. The “survminer” package was used to find the best cutoff values for continuous variables. Using (restricted cubic splines) RCS with three knots positioned at the 10th, 50th, and 90th percentiles to evaluate possible non-linear connections, the relationships between continuous variables and mortality were investigated. To offer a useful tool for predicting personal risk, a nomogram was created. The top 8 most significant factors from our 9 predictors were determined by additional analysis using a random survival forest with 1,000 trees. Iterative testing from 1 to 75 was used to establish the ideal node size. We next performed comparative studies between models built utilizing these top 8 variables and our 9-variable model. Time-dependent ROC curves, calibration plots, and decision curve analyses at 3, 5, and 8-year time points in both training and validation cohorts were used to compare the models. Internal validation was performed via bootstrap resampling: we estimated optimism−corrected concordance index (C−index) values and computed 95% bias−corrected and accelerated (BCa) confidence intervals from the bootstrap resamples.

Results

Patient characteristics and outcomes

The baseline clinical characteristics, treatment of the total cohort (N = 389) are detailed in Supplementary Table S1. Of the 389 patients with CAD combined with diabetes, the median age was 86.00 [82.00;90.00] years, including 299 (76.86%) males and 90 (23.14%) females. The median body mass index (BMI) was 24.46 [22.32;26.75] kg/m². The most common comorbidities were hypertension (341, 87.66%), anemia (172, 44.22%), chronic heart failure (127, 32.65%), and chronic kidney disease (121, 31.11%). In terms of laboratory tests, the median NT-proBNP level was 0.5 × 10³ pg/mL (IQR: 0.18–1.91), and the median HDL level was 1.01 mmol/L (0.84; 1.19). Baseline characteristics are shown in Table 1. The cohort was randomly allocated into a training set (n = 273) and an internal validation set (n = 116). No substantial differences were observed between the training group and the validation group across the majority of features, it indicates that the baseline characteristics of the training cohort and the validation cohort are generally comparable. Supplementary Figure S6 illustrates the distribution of follow-up durations and the number of participants at risk: follow-up durations ranged from less than 1 year to more than 10 years, with a median follow-up time of 4.68 years; 184 participants remained under follow-up at 5 years, and 111 participants remained at risk at 8 years.

Prediction model built based on lasso-cox regression

Parameters were screened using Lasso regression, and Figure 2A displayed the variation features of the variables' coefficients. After applying the 10-fold cross-validation approach to the iterative analysis, a model with a low number of variables and outstanding performance was achieved when λ was 0.063 (Log λ = −1.20) (Figure 2B). Age, CRP, Hb, ALB, INR, NT-proBNP, fibrinogen, digitalis, and diuretics were among the factors that were screened. Lasso regression-screened parameters were used to further develop the Cox regression model (Table 2). Detailed coefficients of all selected variables are presented in Supplementary Table S2. Variance inflation factor analysis confirmed the absence of significant multicollinearity among these predictors (Supplementary Table S3).

Figure 2

The predictive accuracy of the new prognostic model

Time-dependent receiver operating characteristic (ROC) curves for 3-, 5-, and 8-year death forecasts were used to assess the model's discriminative ability (Figure 3). These results were confirmed in the validation cohort, with AUCs of 0.824 (95% CI: 0.734–0.915), 0.813 (95% CI: 0.73–0.896), and 0.798 (95% CI: 0.72–0.877) for 3-, 5-, and 8-year mortality (Figure 3B). The training set's C-index was 0.753 (95% CI: 0.722–0.784), and the validation set's was 0.74 (95% CI: 0.685–0.795). These C-index estimates and 95% confidence intervals are optimism-corrected values derived using BCa bootstrap resampling; the validation cohort therefore represents an internal validation set.

Figure 3

Model calibration assessment

Calibration plots comparing expected and observed survival probabilities at 3-, 5-, and 8-year time intervals were used to assess the model's calibration (Figure 4). All three time points’ calibration curves nearly matched the 45-degree diagonal line, showing that the anticipated probability and the actual observations agreed well.

Figure 4

Clinical decision curve analysis

Within a reasonable range of threshold probabilities (roughly 0.1 to 0.6), the model outperformed both “treat all” and “treat none” strategies in the training cohort, demonstrating consistent net benefit across a range of threshold probabilities at 3-year (Figure 5A), 5-year (Figure 5C), and 8-year (Figure 5E) follow-up. At the 3-year (Figure 5B), 5-year (Figure 5D), and 8-year (Figure 5F) time points, the validation cohort showed comparable patterns of net benefit, indicating its possible therapeutic utility in directing treatment choices.

Figure 5

Construction of a prediction nomogram in training cohort

A nomogram for predicting 3-, 5-year and 8-years overall survival of coronary artery disease combined with diabetes patients that integrated the above 8 independent prognostic factors is shown in Figure 6. The sum of the points allotted to each predictor—which ranged from 0 to 400 points—was used to determine the total points. The nomogram score of age ranged from 60 to 100, making the largest contribution to the prognosis of coronary artery disease combined with diabetes.

Figure 6

Kaplan–Meier analysis

In the Kaplan–Meier analysis, the model displayed good capability in coronary artery disease combined with diabetes patients into high-risk and low-risk groups (Log-rank P < 0.0001) in the training set (Figure 7).

Figure 7

Random survival forest analysis

We use a training cohort to construct a random survival forest model in order to evaluate the significance of the variables and compare them with the Cox model. Hemoglobin and INR were the best predictors, according to variable importance analysis, followed by albumin, NT-proBNP, patient age, fibrinogen, diuretics and digitalis (Figure 8). After roughly 400 trees, the error rate tends to settle, suggesting that the model has sufficiently converged. In our analysis, we employed Random Survival Forest (RSF) to assess the contribution of each feature to predicting patient outcomes. Additionally, C-index values were compared for both the training and test datasets for both models (Supplementary Figure S3A). The results revealed comparable discriminatory performance, with RSF showing similar or slightly improved performance in both training and test sets. To further evaluate model performance, we compared the Brier scores of Random Survival Forest (RSF) and Cox proportional hazards (CoxPH) models, which are shown in Supplementary Figure S3B. The Brier score comparison demonstrated that RSF had slightly better calibration than CoxPH. To ensure the transparency and interpretability of the model, we utilized SHAP (Shapley Additive Explanations) values, which allow us to explain how each feature influences the model's output. Specifically, we used SHAP summary bar plots (Supplementary Figure S4A) and beeswarm plots (Supplementary Figure S4B) to visualize feature importance. Additionally, we generated SHAP partial dependence plots (Supplementary Figure S5) to explore the non-linear relationships between continuous features and the model's predictions. These plots provide insight into how variations in feature values affect model output, capturing the complex interactions that are often overlooked in traditional linear models. Hyperparameter settings of the random forest are provided in Supplementary Table S5.

Figure 8

Model comparison across different time horizons

We conducted a comprehensive comparison of mortality prediction models across 3-, 5-, and 8-year intervals (Figure 9). In the training cohort, Model A (incorporating all 8 variables) demonstrated superior discriminative performance for 3-year mortality prediction compared to single-variable models (AUC = 0.846, 95% CI: 0.798–0.894). This advantage persisted for 5-year (AUC = 0.838, 95% CI: 0.792–0.884) and 8-year predictions (AUC = 0.82, 95% CI: 0.768–0.872). Model A maintained consistent performance in the validation cohort across all time points: 3-year AUC = 0.824 (95% CI: 0.734–0.915), 5-year AUC = 0.813 (95% CI: 0.73–0.896), and 8-year AUC = 0.798 (95% CI: 0.72–0.877). (Supplementary Figures S1, S2) Calibration plots demonstrated strong agreement between predicted and observed probabilities at all time intervals. Decision curve analysis consistently showed that Model A provided superior net benefit across a wider range of threshold probabilities compared to both simplified models and default strategies.

Figure 9

Survival analysis based on key predictors

Kaplan–Meier analysis was used to evaluate survival differences, with stratification according to key predictors identified in the nomogram analysis (Figure 10). For continuous variables, optimal cutoff values were determined using statistical methods to maximize separation between risk groups (Supplementary Table S4). NT-proBNP (cutoff: 487.2 pg/mL), INR (cutoff: 1.12), age (cutoff: 85 years), fibrinogen (cutoff: 3.72 g/L), and hemoglobin (cutoff: 116 g/L) all demonstrated significant survival outcome stratification (P < 0.0001 for all). For categorical variables, the use of diuretics and digoxin also showed strong prognostic stratification (P < 0.0001).

Figure 10

Nonlinear association analysis

To investigate the nonlinear relationships between continuous predictors and mortality risk, we conducted adjusted analyses of key predictors identified in multivariable Cox regression models, sequentially excluding each variable under investigation (Figure 11). In adjusted analyses, NT-proBNP, albumin, and INR demonstrated significant nonlinear associations with mortality (global P < 0.001, nonlinearity P < 0.001). Age and hemoglobin exhibited strong linear associations with mortality (global P < 0.001) with marginally significant nonlinear components (nonlinearity P = 0.541 and 0.182, respectively). These patterns remained consistent across both unadjusted and adjusted analyses, confirming the robustness of these relationships independent of other clinically significant predictors.

Figure 11

Discussion

In this comprehensive investigation of outcomes in patients with coronary artery disease and diabetes comorbidity, we developed and validated a multimodal prediction model. The model demonstrated robust predictive performance, with notable mortality discrimination accuracy and consistent calibration across various time intervals. The complementary nature of these parameters suggests that multimodal approaches may more effectively capture the complex pathophysiology underlying disease progression in this comorbid population. To facilitate clinical decision-making, we developed a nomogram based on our model, which allows clinicians to estimate a patient's survival probabilities at 3, 5, and 8 years. The nomogram incorporates key predictors such as age, hemoglobin levels, albumin, INR, and NT-proBNP, which are factors that clinicians commonly encounter in daily practice. By assigning a score to each of these variables, the nomogram calculates a total points score that can be translated into predicted survival probabilities at different time intervals. This tool can be integrated into clinical workflows as a decision-support system, assisting clinicians in risk-adjusted treatment planning tailored to individual patients. For instance, when faced with a patient exhibiting specific risk factors (e.g., advanced age or low albumin levels), a clinician can use the nomogram to estimate the patient's long-term survival prognosis, which may influence decisions regarding aggressive interventions, frequency of follow-up visits, or the need for specialized care.

Our study identified anemia as a significant prognostic predictor in patients with comorbid coronary artery disease and diabetes. Anemia, a known contributor to ischemic conditions, reduces blood oxygen capacity, causing myocardial hypoxia. Compensatory mechanisms such as increased heart rate and cardiac output exacerbate pre-existing myocardial ischemia (15). Furthermore, anemia frequently coexists with chronic kidney disease (CKD). Impaired renal erythropoietin (EPO) production in CKD patients significantly elevates anemia risk (16). This creates a pathophysiological cycle where anemia worsens renal hypoxia, accelerating kidney function deterioration. Importantly, renal insufficiency promotes toxin accumulation, intensifying systemic inflammatory responses and oxidative stress that compound cardiovascular damage (17). Current evidence demonstrates no conclusive benefit of anemia correction for improving clinical outcomes in cardiac patients (18). Robust clinical data indicate that erythropoiesis-stimulating agents (ESAs) provide neither symptomatic relief nor outcome improvement in those with mild-to-moderate anemia, while potentially increasing serious adverse risks (19).

Inflammation serves as a key driver in the initiation and progression of atherosclerosis (20). In diabetic patients, inflammatory responses exacerbate vascular endothelial damage, accelerating lipid deposition and plaque formation. Furthermore, these responses activate the coagulation system, elevating clotting factor levels and platelet activity, thereby increasing thrombotic risk. This reciprocal interaction between coagulation and inflammation creates a self-reinforcing cycle that accelerates atherosclerotic progression (21, 22). INR levels reflect underlying coagulation dysfunction and correlate with clinical outcomes. Elevated INR levels demonstrate strong prognostic value, showing consistent associations with adverse outcomes in cardiovascular disease patients (23–25). The conversion of fibrinogen to insoluble fibrin and subsequent formation of stable clots represents the terminal phase of the coagulation cascade. Impaired fibrin clot architecture is widely recognized as both a prognostic marker in CAD patients and a potential indicator of thromboembolic risk (26, 27). Furthermore, hyperfibrinogenemia may constitute a key mechanism underlying elevated cardiovascular risk in T2DM populations (28, 29).

Albumin, the most abundant circulating protein in humans, exhibits critical physiological functions including anti-inflammatory effects, antioxidant activity, and inhibition of platelet aggregation (30). As the predominant antioxidant in whole blood, serum albumin derives its antioxidant capacity primarily from free cysteine residues. The prognostic significance of serum albumin in cardiovascular diseases stems principally from its associations with nutritional status and systemic inflammation (31). Albumin additionally modulates inflammatory and immune responses while maintaining vascular integrity (32). Novel albumin-based inflammatory indices such as the C-reactive protein-albumin-lymphocyte (CALLY) index (33) and C-reactive protein/albumin ratio (CAR) (34) enhance prognostic assessment by quantifying systemic inflammatory burden. Reduced serum albumin levels demonstrate associations with ischemic heart disease, heart failure, atrial fibrillation, stroke, and venous thromboembolism (35). In diabetic patients, albumin's antioxidant properties become compromised, acquiring pro-oxidant characteristics (36). Elena The work of Succurro et al. (37) reveals that reduced albumin in diabetes may upregulate Nox2 activity—a mediator of pro-inflammatory responses and vasoconstriction—potentially contributing to cardiovascular events.

Cardiac contraction primarily relies on energy derived from mitochondrial fatty acid oxidation, a process supported by the abundance of mitochondria in myocardial tissue. NT-proBNP serves as a key prognostic biomarker for risk stratification in heart failure patients (38). Pathological cardiac hypertrophy is characterized by reduced fatty acid oxidation alongside enhanced glucose uptake and glycolytic activity. In T2DM, insulin resistance restricts myocardial glucose utilization for energy production, while increased fatty acid oxidation, lipid accumulation, and oxygen consumption collectively impair cardiac function (39). The Malachias et al. study identified NT-proBNP as a significant predictor of mortality and cardiovascular events in diabetic populations (40). Incorporating both BNP and NT-proBNP into optimized risk prediction models enhances discriminative capacity beyond established cardiovascular outcome predictors (41). Serial NT-proBNP measurements demonstrate superior predictive value over single assessments for forecasting cardiovascular mortality or heart failure hospitalization in patients with coronary artery disease and T2DM (42).

Our study revealed that diuretic and digoxin use in patients with comorbid coronary artery disease and diabetes was associated with poorer clinical outcomes, likely attributable to their more complex comorbidities and disease progression toward heart failure. However, prior research has demonstrated that loop diuretic use in coronary artery disease patients without systolic heart failure or renal impairment correlates with increased cardiovascular and non-cardiovascular mortality (43). These findings emphasize the necessity for judicious prescription of these medications in patients lacking clear clinical indications.

Given the multifactorial, multisystem complexity of coronary artery disease (CAD) with comorbid diabetes, there is an urgent clinical need for precise prognostic models to guide therapeutic decision-making. Machine learning demonstrates significant potential in risk prediction, early diagnosis, and personalized management of CAD with diabetes. By analyzing extensive clinical datasets, these algorithms can identify complex patterns potentially overlooked by conventional statistical methods, offering novel insights for early intervention (44–46). Our investigation employed a comprehensive multimodal approach with rigorous validation protocols and clinically relevant time frames. The model demonstrated robust discriminative performance and calibration accuracy across temporal intervals, suggesting practical utility for both short- and long-term risk stratification. However, limitations include the single-center study design with a restricted sample size (n = 389), which may limit generalizability, as well as the lack of external validation using multi-ethnic data. Future validation in larger, multicenter cohorts is required to confirm clinical applicability. Given the relatively small sample size of 389 patients and the inclusion of multiple predictors, the risk of overfitting is a concern, especially in a longitudinal study with potential patient loss to follow-up. To mitigate this, we employed Random Survival Forest (RSF), which reduces overfitting through ensemble learning, enhances model generalizability, and allows for the interpretability of predictor contributions. Another limitation is that we did not implement a formal time-dependent Cox model; medication exposures were aggregated over the first six months and treated as baseline covariates, so residual immortal time bias may remain. In addition to the methods used in this study, we acknowledge the potential benefits of exploring auxiliary model-based optimization methods and hybrid nonlinear identification algorithms, such as those discussed by Ali et al. (47, 48), which could offer further improvements in predictive accuracy. These advanced techniques are particularly beneficial for handling complex, non-linear relationships in clinical data, and could enhance the model's robustness and adaptability to diverse patient populations (49, 50). Specifically, the Runge-Kutta optimization algorithm for fractional input nonlinear output error systems, as well as Chameleon Swarm Optimization and Mountain Gazelle optimization algorithms, offer biologically inspired, auxiliary model-driven frameworks that can improve both interpretability and dynamic response estimation of medical survival models, enabling more precise and reliable predictions in clinical settings. We are actively considering these approaches for future iterations of our model to improve its performance in more dynamic clinical settings.

Conclusion

The advancement of machine learning has established random forest models as a robust methodology for developing medical prediction tools. Integration of refined biomarkers with artificial intelligence-driven models shows significant potential to enhance prognostic assessment in patients with coronary artery disease and comorbid diabetes. This study conducted a comprehensive analysis of clinical parameters associated with mortality risk in coronary artery disease patients with diabetes comorbidity. Our findings highlight the critical predictive value of hemoglobin, INR, albumin, NT-proBNP, age, fibrinogen, diuretic use, and digitalis therapy in risk stratification. The developed model demonstrates reliable short-term prognostic accuracy using a parsimonious set of clinical variables. This advancement holds promise for enabling personalized therapeutic strategies to optimize clinical outcomes in this high-risk population.

References

• 1.

  Roth GA Mensah GA Johnson CO Addolorato G Ammirati E Baddour LM et al Global burden of cardiovascular diseases and risk factors, 1990–2019: update from the GBD 2019 study. J Am Coll Cardiol. (2020) 76(25):2982–3021. 10.1016/j.jacc.2020.11.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2.

  Li Y Liu Y Cui J Zhu M Wang W Chen K et al Oral-gut microbial transmission promotes diabetic coronary heart disease. Cardiovasc Diabetol. (2024) 23(1):123. 10.1186/s12933-024-02217-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  Marx N Federici M Schütt K Müller-Wieland D Ajjan RA Antunes MJ et al 2023 ESC guidelines for the management of cardiovascular disease in patients with diabetes. Eur Heart J. (2023) 44(39):4043–140. 10.1093/eurheartj/ehad192

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  Goodarzi MO Rotter JI . Genetics insights in the relationship between type 2 diabetes and coronary heart disease. Circ Res. (2020) 126(11):1526–48. 10.1161/circresaha.119.316065

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  Kragelund C Grønning B Køber L Hildebrandt P Steffensen R . N-terminal pro-B-type natriuretic peptide and long-term mortality in stable coronary heart disease. N Engl J Med. (2005) 352(7):666–75. 10.1056/NEJMoa042330

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  American Diabetes Association Professional Practice Committee. 2. Diagnosis and classification of diabetes: standards of care in diabetes-2024. Diabetes Care. (2024) 47(1):S20–42. 10.2337/dc24-S002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  D'Agostino RB Sr. Vasan RS Pencina MJ Wolf PA Cobain M Massaro JM et al General cardiovascular risk profile for use in primary care: the framingham heart study. Circulation. (2008) 117(6):743–53. 10.1161/circulationaha.107.699579

  • CrossRef
  • Google Scholar

• 8.

  Pavlović J Bos D Ikram MK Ikram MA Kavousi M Leening MJG . Guideline-directed application of coronary artery calcium scores for primary prevention of atherosclerotic cardiovascular disease. JACC Cardiovasc Imaging. (2025) 18(4):465–75. 10.1016/j.jcmg.2024.12.008

  • CrossRef
  • Google Scholar

• 9.

  Cook NR Paynter NP Eaton CB Manson JE Martin LW Robinson JG et al Comparison of the framingham and reynolds risk scores for global cardiovascular risk prediction in the multiethnic Women’s health initiative. Circulation. (2012) 125(14):1748–56, s1–11. 10.1161/circulationaha.111.075929

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  Selvarajah S Kaur G Haniff J Cheong KC Hiong TG van der Graaf Y et al Comparison of the framingham risk score, SCORE and WHO/ISH cardiovascular risk prediction models in an Asian population. Int J Cardiol. (2014) 176(1):211–8. 10.1016/j.ijcard.2014.07.066

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  Zhang H Chong H Li Z Li K Zhang B Xue Y et al Triglyceride-glucose index in the prediction of major adverse cardiovascular events in patients with type 2 diabetes mellitus after coronary artery bypass surgery: a retrospective cohort study. Front Endocrinol. (2022) 13:1015747. 10.3389/fendo.2022.1015747

  • CrossRef
  • Google Scholar

• 12.

  Wei S Guo X He S Zhang C Chen Z Chen J et al Application of machine learning for patients with cardiac arrest: systematic review and meta-analysis. J Med Internet Res. (2025) 27:e67871. 10.2196/67871

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  Knuuti J Wijns W Saraste A Capodanno D Barbato E Funck-Brentano C et al 2019 ESC guidelines for the diagnosis and management of chronic coronary syndromes. Eur Heart J. (2020) 41(3):407–77. 10.1093/eurheartj/ehz425

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14.

  Alberti KG Zimmet PZ . Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet Med. (1998) 15(7):539–53. 10.1002/(sici)1096-9136(199807)15:7%3C539::Aid-dia668%3E3.0.Co;2-s

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  da Silveira AD Ribeiro RA Rossini AP Stella SF Ritta HA Stein R et al Association of anemia with clinical outcomes in stable coronary artery disease. Coron Artery Dis. (2008) 19(1):21–6. 10.1097/MCA.0b013e3282f27c0a

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16.

  Zafrir B Leviner DB Saliba W Sharoni E . Prognostic interplay of chronic kidney disease, Anemia, and diabetes in coronary bypass surgery. Ann Thorac Surg. (2021) 111(1):94–101. 10.1016/j.athoracsur.2020.04.124

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  Carracedo J Alique M Vida C Bodega G Ceprián N Morales E et al Mechanisms of cardiovascular disorders in patients with chronic kidney disease: a process related to accelerated senescence. Front Cell Dev Biol. (2020) 8:185. 10.3389/fcell.2020.00185

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18.

  Maraey A Elzanaty A Elsharnoby H Younes M Salem M Khalil M . Restrictive vs liberal blood transfusion strategy in patients with symptomatic coronary artery disease: meta-analysis from randomized controlled trials. Eur Heart J. (2021) 42(1):ehab724.1257. 10.1093/eurheartj/ehab724.1257

  • CrossRef
  • Google Scholar

• 19.

  Kansagara D Dyer E Englander H Fu R Freeman M Kagen D . Treatment of anemia in patients with heart disease: a systematic review. Ann Intern Med. (2013) 159(11):746–57. 10.7326/0003-4819-159-11-201312030-00007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20.

  Soehnlein O Libby P . Targeting inflammation in atherosclerosis—from experimental insights to the clinic. Nat Rev Drug Discovery. (2021) 20(8):589–610. 10.1038/s41573-021-00198-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  Bryk AH Prior SM Plens K Konieczynska M Hohendorff J Malecki MT et al Predictors of neutrophil extracellular traps markers in type 2 diabetes mellitus: associations with a prothrombotic state and hypofibrinolysis. Cardiovasc Diabetol. (2019) 18(1):49. 10.1186/s12933-019-0850-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22.

  Babes EE Bustea C Behl T Abdel-Daim MM Nechifor AC Stoicescu M et al Acute coronary syndromes in diabetic patients, outcome, revascularization, and antithrombotic therapy. Biomed Pharmacother. (2022) 148:112772. 10.1016/j.biopha.2022.112772

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23.

  Tang Y Sun J Yu Z Liang B Peng B Ma J et al Association between prothrombin time-international normalized ratio and prognosis of post-cardiac arrest patients: a retrospective cohort study. Front Public Health. (2023) 11:1112623. 10.3389/fpubh.2023.1112623

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24.

  Delgado GE Zirlik A Gruber R Scheffold T Krämer BK März W et al The association of high-normal international-normalized-ratio (INR) with mortality in patients referred for coronary angiography. PLoS One. (2019) 14(8):e0221112. 10.1371/journal.pone.0221112

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25.

  Liu L Ying M Chen S Li Q Chen G Li H et al The association between prothrombin time-international normalized ratio and long-term mortality in patients with coronary artery disease: a large cohort retrospective study with 44,662 patients. BMC Cardiovasc Disord. (2022) 22(1):297. 10.1186/s12872-022-02619-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26.

  Ząbczyk M Ariëns RAS Undas A . Fibrin clot properties in cardiovascular disease: from basic mechanisms to clinical practice. Cardiovasc Res. (2023) 119(1):94–111. 10.1093/cvr/cvad017

  • CrossRef
  • Google Scholar

• 27.

  Larsen JB Hvas AM . Fibrin clot properties in coronary artery disease: new determinants and prognostic markers. Pol Arch Intern Med. (2021) 131(11):16113. 10.20452/pamw.16113

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28.

  Ravi Kiran M Harshavardhan L . Fibrinogen levels in type 2 diabetes Mellitus and it’s correlation with microvascular microvascular complications. J Assoc Physicians India. (2022) 70(4):11–2. eng. Epub 2022/04/22. Cited in: Pubmed; PMID 35443328.

  • Google Scholar

• 29.

  Bruno G Cavallo-Perin P Bargero G Borra M D'Errico N Pagano G . Association of fibrinogen with glycemic control and albumin excretion rate in patients with non-insulin-dependent diabetes mellitus. Ann Intern Med. (1996) 125(8):653–7. 10.7326/0003-4819-125-8-199610150-00005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30.

  Manolis AA Manolis TA Melita H Mikhailidis DP Manolis AS . Low serum albumin: a neglected predictor in patients with cardiovascular disease. Eur J Intern Med. (2022) 102:24–39. 10.1016/j.ejim.2022.05.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31.

  Don BR Kaysen G . Serum albumin: relationship to inflammation and nutrition. Semin Dial. (2004) 17(6):432–7. 10.1111/j.0894-0959.2004.17603.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32.

  Tufoni M Baldassarre M Zaccherini G Antognoli A Caraceni P . Hemodynamic and systemic effects of albumin in patients with advanced liver disease. Curr Hepatol Rep. (2020) 19(3):147–58. 10.1007/s11901-020-00521-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33.

  Ji H Luo Z Ye L He Y Hao M Yang Y et al Prognostic significance of C-reactive protein-albumin-lymphocyte (CALLY) index after primary percutaneous coronary intervention in patients with ST-segment elevation myocardial infarction. Int Immunopharmacol. (2024) 141:112860. 10.1016/j.intimp.2024.112860

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34.

  Di Rosa M Sabbatinelli J Giuliani A Carella M Magro D Biscetti L et al Inflammation scores based on C-reactive protein and albumin predict mortality in hospitalized older patients independent of the admission diagnosis. Immun Ageing. (2024) 21(1):67. 10.1186/s12979-024-00471-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  Arques S . Human serum albumin in cardiovascular diseases. Eur J Intern Med. (2018) 52:8–12. 10.1016/j.ejim.2018.04.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36.

  Roche M Rondeau P Singh NR Tarnus E Bourdon E . The antioxidant properties of serum albumin. FEBS Lett. (2008) 582(13):1783–7. 10.1016/j.febslet.2008.04.057

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  Succurro E Andreozzi F Carnevale R Sciacqua A Cammisotto V Cassano V et al Nox2 up-regulation and hypoalbuminemia in patients with type 2 diabetes mellitus. Free Radical Biol Med. (2021) 168:1–5. 10.1016/j.freeradbiomed.2021.03.026

  • CrossRef
  • Google Scholar

• 38.

  Jhund PS Anand IS Komajda M Claggett BL McKelvie RS Zile MR et al Changes in N-terminal pro-B-type natriuretic peptide levels and outcomes in heart failure with preserved ejection fraction: an analysis of the I-preserve study. Eur J Heart Fail. (2015) 17(8):809–17. 10.1002/ejhf.274

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39.

  Malachias MVB Wijkman MO Bertoluci MC . NT-proBNP as a predictor of death and cardiovascular events in patients with type 2 diabetes. Diabetol Metab Syndr. (2022) 14(1):64. 10.1186/s13098-022-00837-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40.

  Malachias MVB Jhund PS Claggett BL Wijkman MO Bentley-Lewis R Chaturvedi N et al NT-proBNP by itself predicts death and cardiovascular events in high-risk patients with type 2 diabetes Mellitus. J Am Heart Assoc. (2020) 9(19):e017462. 10.1161/jaha.120.017462

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41.

  Wolsk E Claggett B Pfeffer MA Diaz R Dickstein K Gerstein HC et al Role of B-type natriuretic peptide and N-terminal prohormone BNP as predictors of cardiovascular morbidity and mortality in patients with a recent coronary event and type 2 diabetes Mellitus. J Am Heart Assoc. (2017) 6(6):e004743. 10.1161/jaha.116.004743

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42.

  Wolsk E Claggett B Diaz R Dickstein K Gerstein HC Køber L et al Risk estimates of imminent cardiovascular death and heart failure hospitalization are improved using serial natriuretic peptide measurements in patients with coronary artery disease and type 2 diabetes. J Am Heart Assoc. (2022) 11(8):e021327. 10.1161/jaha.121.021327

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43.

  Schartum-Hansen H Loeland KH Svingen GFT Nordrehaug JE Nygaard O . Use of loop diuretics is strongly associated with increased risk of death in coronary patients without systolic heart failure or renal impairment. Eur Heart J. (2013) 34(1):878. 10.1093/eurheartj/eht308.878

  • CrossRef
  • Google Scholar

• 44.

  Gautam N Mueller J Alqaisi O Gandhi T Malkawi A Tarun T et al Machine learning in cardiovascular risk prediction and precision preventive approaches. Curr Atheroscler Rep. (2023) 25(12):1069–81. 10.1007/s11883-023-01174-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45.

  Oikonomou EK Khera R . Machine learning in precision diabetes care and cardiovascular risk prediction. Cardiovasc Diabetol. (2023) 22(1):259. 10.1186/s12933-023-01985-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46.

  Ji Y Shang H Yi J Zang W Cao W . Machine learning-based models to predict type 2 diabetes combined with coronary heart disease and feature analysis-based on interpretable SHAP. Acta Diabetol. (2025) 62(10):1631–46. 10.1007/s00592-025-02496-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47.

  Samadi ME Guzman-Maldonado J Nikulina K Mirzaieazar H Sharafutdinov K Fritsch SJ et al A hybrid modeling framework for generalizable and interpretable predictions of ICU mortality across multiple hospitals. Sci Rep. (2024) 14(1):5725. 10.1038/s41598-024-55577-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48.

  Ali MA Chaudhary NI Khan TA Mao W-L Lin C-C Raja MAZ . Design of key term separated identification model for fractional input nonlinear output error systems: auxiliary model based Runge Kutta optimization algorithm. Chaos, Solitons Fractals. (2024) 189:115696. 10.1016/j.chaos.2024.115696

  • CrossRef
  • Google Scholar

• 49.

  Ali MA Chaudhary NI Lin C-C Mao W-L . Mountain gazelle optimization algorithm for identification of nonlinear hammer-stein output error systems. J Innov Technol. (2025) 7(1):117–24. 10.29424/JIT.202503_7(1).0010d

  • CrossRef
  • Google Scholar

• 50.

  Ali MA Chaudhary NI Khan TA Mao W-L Lin C-C Khan ZA et al Auxiliary model-based chameleon swarm optimization for robust parameter estimation of fractional order nonlinear hammerstein systems. J Comput Nonlinear Dyn. (2025) 20(9):091005. 10.1115/1.4068718

  • CrossRef
  • Google Scholar