ORIGINAL RESEARCH article

Front. Pharmacol., 28 August 2026

Sec. Pharmacoepidemiology

Volume 17 - 2026 | https://doi.org/10.3389/fphar.2026.1825631

Clinical, laboratory, and medication-related factors associated with in-hospital mortality among oncology patients at a tertiary hospital in Saudi Arabia

  • 1. Department of Pharmacology and Toxicology, College of Pharmacy, University of Hail, Hail, Saudi Arabia

  • 2. College of Pharmacy, University of Hail, Hail, Saudi Arabia

  • 3. Institute of Pharmaceutical Sciences (IPS), University of Veterinary and Animal Sciences (UVAS), Lahore, Pakistan

  • 4. Department of Clinical Pharmacy, Faculty of Pharmacy, Al-Baha University, Al-Baha, Saudi Arabia

  • 5. School of Dentistry and Medical Sciences, Charles Sturt University, Orange, NSW, Australia

Abstract

Background:

Oncology inpatients face physiological decline alongside complex pharmacology, yet clinical mortality predictors and drug-related risks remain understudied in Middle Eastern populations. This study identified independent predictors of in-hospital mortality and characterized potential drug–drug interactions (pDDIs) among hospitalized cancer patients in Saudi Arabia.

Methods:

This retrospective cohort study evaluated 401 adult oncology inpatients admitted to a tertiary care center in Hail, Saudi Arabia (2020–2022). Demographics, clinical parameters, and medication profiles were extracted from medical records. Bivariate correlations were assessed using point-biserial () and phi () coefficients. Independent mortality predictors were identified via multivariable binary logistic regression, and biomarker discrimination was evaluated using receiver operating characteristic (ROC) curves. Potential DDIs were screened using Drugs.com.

Results:

The cohort was predominantly female (67.6%) and Saudi (85.0%), with solid malignancies accounting for 94.3% of cases. In-hospital mortality was 25.2% (). Bivariate correlations demonstrated significant associations with mortality across clinical parameters, including mechanical ventilation requirements (). Multivariable regression identified four independent physiological predictors of mortality: renal failure (), baseline albumin , serum creatinine (), and oxygen saturation (). Serum albumin showed high individual discrimination (), while the combined multivariable model achieved an of 0.953. Pharmacological screening revealed near-universal interaction prevalence: 98.6% of 2,148 prescriptions contained a pDDI flag (54.9% severe), and 96.5% of patients had severe pDDI flag.

Conclusion:

In-hospital mortality is strongly driven by acute organ dysfunction and physiological decline rather than database-flagged interactions alone. However, the high prevalence of severe pDDI flags highlights a critical need to integrate structured clinical pharmacy reviews into oncology admission protocols. Targeted medication safety interventions should focus on avoiding nephrotoxic combinations and establishing pre-therapeutic pharmacogenetic screening to mitigate potential iatrogenic risks.

1 Introduction

Cancer remains a leading cause of premature death worldwide, posing one of the greatest continuing public health challenges of our time. According to one survey carried out by the World Health Organization (WHO), people in most countries can expect cancer to be ranked among the top three causes of death before reaching 70 years of age (). In 2020, nearly 10 million people died from cancer globally, and this figure is projected to exceed 16 million per year by 2040, driven primarily by population ageing, tobacco exposure, obesity, and metabolic disturbances (; ). In many developing and middle-income nations, tobacco consumption, rising obesity rates, and metabolic comorbidities have substantially altered cancer epidemiological patterns over recent decades ().

Cancer patients constitute one of the most pharmacologically complex populations in contemporary clinical practice. Beyond cytotoxic chemotherapy, they typically receive corticosteroids, antiemetics, proton pump inhibitors, anticoagulants, antibiotics, antifungals, analgesics, and various supportive care agents concurrently. This polypharmaceutical environment creates a substantial risk of clinically meaningful drug - drug interactions (DDIs). Most antineoplastic agents - taxanes and platinum-based regimens in particular - are metabolised through cytochrome P450 pathways, rendering them vulnerable to pharmacokinetic interactions with co-prescribed drugs (). Additive toxicity through pharmacodynamic mechanisms, including nephrotoxicity, myelosuppression, and QT prolongation, compounds this risk further and is difficult to anticipate without systematic medication review. Hospitalised cancer patients carry a higher burden of adverse drug events than non-oncology populations, a pattern attributable in part to polypharmacy and cumulative DDI exposure ().

Despite progress in prevention and treatment, clinical outcomes at discharge vary considerably across populations, attributed to differences in age distribution, comorbidity burden, and healthcare access (). Available evidence suggests that age, organ dysfunction, and the presence of significant comorbidities are among the key determinants of cancer prognosis (). Specifically, conditions such as diabetes mellitus, renal failure, and hypertension accelerate disease progression and frequently limit the delivery of full therapeutic doses, thereby impacting overall discharge status (; ). Additionally, smoking remains a well-established independent predictor of in-hospital mortality across multiple cancer types (). Over the past 2 decades, both incidence and the mortality proportion from cancer have risen steadily in Saudi Arabia (). Colorectal, breast, and hematologic malignancies account for the greatest clinical burden, yet vital status at discharge within each category depends considerably on the individual patient profile and associated risk factors (; ). Prior work has concentrated largely on site-specific cancers, leaving a gap in the evidence base for studies that assess these short-term clinical outcomes under the combined influence of demographic and clinical variables ().

While numerous studies have evaluated mortality predictors in critically ill patients, fewer have focused specifically on hospitalised oncology patients outside intensive care settings. Mortality in this population is multifactorial: infection, organ dysfunction, treatment-related toxicity, and disease progression each contribute in varying degrees. Identifying modifiable predictors, such as laboratory derangements or medication-related risk factors, may offer opportunities for earlier intervention (). Evaluating both clinical and pharmacological determinants of in-hospital mortality therefore provides a more complete basis for risk stratification in this patient group.

Although previous research has increasingly recognized the clinical burden of potential drug-drug interactions (pDDIs) in oncology, literature characterizing these events remains skewed toward specific care environments. Historically, investigations into inpatient drug safety have primarily focused on critical care or intensive care settings, where complex medication changes are standard (). In contrast, studies targeting non-intensive care unit (non-ICU) hospitalized oncology patients are relatively sparse, spanning isolated international and regional cohorts. Prior regional assessments have evaluated the prevalence and distinct predictors of pDDIs in general inpatient oncology wards, identifying advanced age, prolonged length of hospital stay, the presence of hematologic malignancies, and underlying polypharmacy as prominent drivers of exposure (; ). Furthermore, observed clinical outcomes in these populations range from prolonged hospitalization to severe toxicities, such as dangerous QT-interval prolongations resulting from antiemetic and systemic chemotherapy co-prescription ().

However, these existing studies suffer from collective limitations; they are frequently restricted by single-center designs, small sample sizes, or a narrow focus on exclusive drug classes—such as interactions driven solely by antineoplastic regimens or supportive care agents. Crucially, there is a distinct lack of comprehensive literature integrating full profile pDDI screening data that captures the overlapping use of chemotherapy, supportive medications, and drugs managed for concurrent, non-oncological comorbidities within general oncology floors. This omission creates a significant knowledge gap, as it leaves the full scope of real-world medication risk unquantified for hospitalized cancer patients outside the ICU.

To address this gap, the present study retrospectively analyzed 401 adult oncology inpatients treated at a tertiary referral center in Hail, Saudi Arabia. This study is an exploratory observational analysis designed to identify clinical, laboratory, and medication-related factors (specifically the prevalence and severity of potential DDIs) associated with in-hospital mortality among non-ICU hospitalized oncology patients, and to evaluate whether these pDDIs are independently associated with mortality.

2 Materials and methods

2.1 Study design and setting

This retrospective analytical cohort study was conducted at King Salman Hospital, a tertiary care academic medical center and regional referral hub in the Hail region, Saudi Arabia. The oncology service at this facility features a dedicated to inpatient ward. Electronic health records were systematically reviewed over a 3-year calendar period, from 1 January 2022, to 31 December 2024, to identify baseline clinical, laboratory, and pharmacological predictors (including potential drug-drug interactions [pDDIs]) of in-hospital mortality among adult cancer patients. Eligible patients were anchored at hospital admission (“Time Zero”), and baseline predictors were defined as the first clinical or laboratory measurements obtained within 24 h of admission. Post-admission longitudinal data, such as subsequent laboratory fluctuations or mechanical ventilation events, were isolated from the baseline framework and analyzed as markers of hospital-course severity through to hospital discharge. The institutional ethics committee reviewed and approved the study protocol, which was conducted in strict accordance with the Declaration of Helsinki and national guidelines for biomedical research ().

2.2 Participant selection and eligibility

A total of 499 hospitalized oncology patients were initially screened for study eligibility. Of these, 98 patients were excluded for the following reasons: incomplete potential drug-drug interaction (pDDI) screening data (), age under 18 years (), a short hospital stay of less than 48 h (), ICU admission during the study period (), absence of systemic antineoplastic therapy (), or declined participation and other reasons (). Consequently, 401 patients met all eligibility criteria and were enrolled, establishing a final analytical cohort of 401 patients with integrated, comprehensive pDDI data (Figure 1). The cohort included both solid and hematological malignancies and comprised Saudi and non-Saudi nationals, reflecting the mixed patient population of the study center.

FIGURE 1

2.3 Operational definitions

To ensure reproducibility, the following criteria were applied during cohort selection:

Proven Malignancy: Histopathologically or cytologically confirmed solid tumors or hematological malignancies. Supportive/Palliative Care Admission: Non-curative hospitalizations intended strictly for symptom management, pain control, or end-of-life care, rather than active tumor-directed therapeutic interventions. Complete Data: Having full electronic health record documentation for baseline vital signs, initial laboratory panels, and active medication reconciliation sheets at the index admission.

2.4 Data collection and variables

Data were extracted from the hospital electronic health record system and the institutional cancer registry using a standardized data-extraction form. Demographic variables collected included age, sex, and nationality. Comorbid conditions recorded at baseline were diabetes mellitus, hypertension, renal disease (defined as a documented diagnosis of nephropathy or a baseline serum creatinine above the institutional upper reference limit), pulmonary disease, cardiac disease, and hypothyroidism. Smoking status was recorded as a lifestyle variable. Clinical variables included cancer type, disease extent, anticancer therapy profiles, and length of hospital stay. To prevent outcome leakage, all baseline laboratory parameters—including aspartate aminotransferase (AST), alanine aminotransferase (ALT), serum albumin (hypoalbuminemia defined as g/L), total protein, serum creatinine, hemoglobin, red blood cell count, white blood cell count, international normalized ratio (INR), blood pressure measurements, and oxygen saturation—were strictly restricted to the first measurement obtained within 24 h of hospital admission (“Time Zero”). Polypharmacy was defined as five or more medications prescribed concurrently during the index hospitalization. Post-admission events, such as requirements for mechanical ventilation and in-hospital discharge status, were tracked through to discharge and isolated as markers of hospital-course severity ().

2.5 Data validation and quality control

Data were extracted using a structured abstraction form applied systematically across all medical records. Baseline laboratory values were taken from results obtained on admission; where multiple values were available within the first 24 h, the first documented result was used. Medication lists were cross-checked against physician prescribing orders and nursing administration charts to minimize transcription error. To ensure data integrity, a PharmD investigator conducted the data extraction, and a random 10% sample () of all charts was independently double reviewed by the Principal Investigators. Inter-rater reliability was formally assessed using Cohen’s kappa, yielding an excellent agreement score (). Remaining discrepancies were documented and resolved through consensus review involving a third senior investigator

2.6 Drug interaction analysis

Potential drug-drug interactions (pDDIs) were systematically evaluated using the web-based Drugs.com Drug Interaction Checker (Drugsite Trust, Auckland, New Zealand; accessed 2024–2026). This platform functions as a dynamic interface that integrates clinical interaction data from three independent reference engines: IBM Micromedex (Truven Health Analytics, Greenwood Village, CO, United States), the Cerner Multum Drug Information Suite (Cerner Corporation, Kansas City, MO, United States), and the American Society of Health-System Pharmacists (ASHP) Drug Interactions Database (ASHP, Bethesda, MD, United States). The utility of this tool as a screening instrument in retrospective oncology research is supported by a published evaluation of nine DDI screening tools across 145 oral oncology drug pairs, in which Drugs.com achieved a performance score of 0.83, with no statistically significant difference compared to the top-performing subscription-based tool, Lexicomp® (0.89) (). For each drug prescribed during the index admission, the following screening outputs were recorded in a structured database: The potential interaction type, the name of the interacting drug partner (where identified), the potential clinical effect, the screening severity score, the evidence level. Screening severity was classified on a four-point scale: 0 = no significant interaction, 1 = moderate, 2 = severe, and 3 = most severe. Potential interaction types were categorized under one of the following mechanisms or risks: metabolic interaction, toxicity overlap, nephrotoxicity, cardiotoxicity, neuropathy, hepatotoxicity, immune-related toxicity, thrombosis risk, or central nervous system (CNS) depression. Organ-specific interaction flags were recorded at the drug prescription level for potential nephrotoxicity, hepatotoxicity, neuropathy, cardiotoxicity, and thrombosis risk. These were then aggregated to the patient level to determine whether a patient was prescribed at least one medication carrying that specific potential interaction class. Finally, the prevalence of each organ-specific potential interaction category was compared between survivors and non-survivors using the chi-square test, or Fisher’s exact test where expected cell counts fell below five, to evaluate statistical associations with mortality.

2.7 Handling of confounding variables

Potential confounding factors were taken into account, which were dependent on clinical plausibility and previous literature. Variables entering the multivariable model were selected on the basis of statistical significance in univariate analysis (p < 0.05) or known biological plausibility from existing literature. Multicollinearity was assessed using Pearson correlation matrices and variance inflation factors. Highly correlated variables were not entered simultaneously into the regression model.

2.8 Receiver operating characteristics analysis

The discriminatory ability of individual laboratory biomarkers and the combined multivariable logistic regression model to predict in-hospital mortality was assessed using receiver operating characteristic (ROC) curve analysis. The area under the ROC curve (AUC) was calculated for serum albumin, serum creatinine, INR, and oxygen saturation as individual predictors. For continuous variables where higher values indicate lower risk, the variable sign was reversed before ROC analysis to ensure the AUC reflected predictive performance in the direction of mortality. The optimal cut-off value for each biomarker was determined using the Youden index (J = sensitivity + specificity − 1). The combined-model ROC curve was derived from the predicted probabilities of the multivariable logistic regression model, fitted in complete-case analysis to all patients with available data on the five predictor variables (n = 305 of 401 patients). AUC values of 0.70–0.80 represent acceptable discrimination, 0.80–0.90 excellent discrimination, and above 0.90 outstanding discrimination. ROC analysis was performed using Python (scikit-learn version 1.3). All AUC values are reported alongside the Youden-optimal cut-off, sensitivity, and specificity.

2.9 Statistical analysis

Continuous variables were assessed for normality of distribution using the Shapiro–Wilk test alongside visual inspection of normal Q–Q plots. Normally distributed continuous data are presented as mean standard deviation (SD) and compared between groups using the independent-samples -test. Non-normally distributed continuous variables are expressed as median [interquartile range: Q1, Q3] and evaluated using the non-parametric Mann–Whitney test. Categorical variables are expressed as frequencies and percentages () and compared using the Chi-square () test, or Fisher’s exact test if any expected contingency cell count fell below 5.

Continuous and categorical clinical parameters were systematically evaluated for completeness prior to statistical modeling. Missing data were observed in select clinical parameters, most notably oxygen saturation ( missingness) and laboratory indices, concentrated among patients admitted in acute distress. To prevent sample bias and maintain model stability, multivariable binary logistic regression and ROC curve analyses were conducted using a complete-case analysis approach, resulting in an analytical sample of patients with complete records across all model predictors. Variable-level missingness counts and percentages are detailed in Supplementary Table S1

Bivariate associations with the binary outcome of in-hospital mortality (; ) were evaluated using point-biserial correlation () for continuous variables and the phi coefficient () for dichotomous categorical parameters.

To identify independent predictors of in-hospital mortality, a multivariable binary logistic regression model was constructed. Candidate variables were selected a priori based on a clinical conceptual framework and known pathophysiological risks (including organ-specific potential drug–drug interaction (pDDI) exposures) rather than relying solely on automated univariable screening thresholds to ensure model stability. The underlying assumptions of the logistic regression model were rigorously evaluated: logit linearity for continuous independent variables was verified using Box–Tidwell transformations, and multicollinearity was ruled out via Variance Inflation Factors (). Extreme or influential observations were screened using Cook’s distance, confirming that no individual case exceeded the operational threshold of 1.0.

Multivariable regression results are expressed as adjusted odds ratios () with corresponding confidence intervals (). Overall model calibration and fit were assessed using the Nagelkerke pseudo- and the Hosmer–Lemeshow goodness-of-fit test (). The discriminatory performance of individual clinical predictors and the final combined model was evaluated using receiver operating characteristic (ROC) curve analysis to determine the area under the curve (/C-statistic), with optimal cut-off values determined by the Youden index (). All statistical analyses were conducted using IBM SPSS Statistics (Version 27.0; IBM Corp., Armonk, NY, United States). A two-sided was considered statistically significant across all evaluations ().

3 Results

3.1 Baseline demographic and clinical characteristics

The cohort comprised 401 patients, of whom 300 (74.8%) were discharged alive and 101 (25.2%) experienced in-hospital mortality. The population was predominantly female (67.6%) and Saudi (85.0%), though male patients exhibited a notably higher mortality rate than females (36.2% vs. 19.9%). Solid malignancies accounted for the vast majority of cases (94.3%), and patients with metastatic (Stage IV) disease faced a starkly elevated mortality rate of 51.3%, compared to 18.9% for those with non-metastatic cancer. Chronic comorbidities were highly prevalent, led by hypertension (40.4%) and diabetes mellitus (33.7%); however, renal disease (13.5% prevalence) and pulmonary disease (20.2% prevalence) presented the highest subgroup mortality risks at 48.1% and 34.6%, respectively. Additionally, while current or ex-smokers constituted only 5.0% of the cohort, they experienced a substantial mortality rate of 45.0% (Table 1).

TABLE 1

CharacteristicFull cohort (N = 401) n (%)Discharged alive (n = 300) n (%)In-hospital mortality (n = 101) n (%)
Sex
Female271 (67.6)217 (80.1)54 (19.9)
Male130 (32.4)83 (63.8)47 (36.2)
Nationality
Saudi341 (85.0)256 (75.1)85 (24.9)
Non-Saudi60 (15.0)44 (73.3)16 (26.7)
Cancer type
Solid malignancy378 (94.3)279 (73.8)99 (26.2)
Hematologic malignancy23 (5.7)21 (91.3)2 (8.7)
Disease extent
Non-metastatic323 (80.5)262 (81.1)61 (18.9)
Metastatic (Stage IV)78 (19.5)38 (48.7)40 (51.3)
Comorbidities
Hypertension162 (40.4)119 (73.5)43 (26.5)
Diabetes mellitus135 (33.7)99 (73.3)36 (26.7)
Pulmonary disease81 (20.2)53 (65.4)28 (34.6)
Cardiac disease76 (19.0)55 (72.4)21 (27.6)
Renal disease54 (13.5)28 (51.9)26 (48.1)
Smoking status
Non-smoker381 (95.0)289 (75.9)92 (24.1)
Current/ex-smoker20 (5.0)11 (55.0)9 (45.0)

Baseline demographic and clinical characteristics of the patient cohort stratified by discharge and mortality status ().

Data are expressed as number (percentage). Total cohort percentages are calculated based on the full sample (). Survival status percentages (Discharged Alive vs. In-Hospital Mortality) are calculated horizontally based on the row total for each specific baseline characteristic.

3.2 Bivariable associations with in-hospital mortality

Bivariable analyses evaluating clinical and demographic factors against vital status at discharge are summarized in Table 2. Bivariate analysis revealed that disease extent, renal disease, biological sex, and pulmonary disease were statistically significant predictors of in-hospital mortality (). The strongest clinical correlations were observed for metastatic disease and baseline renal disease (both ), which demonstrated starkly elevated subgroup mortality rates of 51.3% (vs. 18.9% for non-metastatic) and 48.1% (vs. 21.6% for absent), respectively. Male sex was also highly associated with increased mortality compared to female sex (36.2% vs. 19.9%; ), while baseline pulmonary disease carried a significantly higher mortality risk of 34.6% (vs. 22.8%; ). Conversely, smoking status showed an adverse clinical trend that did not reach statistical significance, likely due to small subgroup size (45.0% vs. 24.1%; ), and cancer type demonstrated no significant correlation with survival outcomes ().

TABLE 2

VariableCategoryTotal (N)Mortality n (%)Test appliedExact p-value
SexFemale
Male
271
130
54 (19.9)
47 (36.2)
Pearson 0.001
Cancer typeSolid
Hematologic
378
23
99 (26.2)
2 (8.7)
Fisher’s exact0.103
Disease extentMetastatic
Non-metastatic
78
323
40 (51.3)
61 (18.9)
Pearson <0.001
Renal diseasePresent
Absent
54
347
26 (48.1)
75 (21.6)
Pearson <0.001
Pulmonary diseasePresent
Absent
81
320
28 (34.6)
73 (22.8)
Pearson 0.042
Smoking statusCurrent/ex-smoker
Non-smoker
20
381
9 (45.0)
92 (24.1)
Fisher’s exact0.067

Association between selected demographic and clinical variables and in-hospital mortality ().

Percentages represent row-wise in-hospital mortality within each variable category. Statistical significance was evaluated at the level. Test Applied: Pearson indicates Pearson’s Chi-squared test; Fisher’s Exact indicates Fisher’s exact test used for subgroups with small, expected cell counts ().

3.3 Continuous clinical parameters and survival outcomes

Table 3 evaluates normally distributed continuous clinical parameters between patients who were discharged alive () and those who experienced in-hospital mortality (). Patients who died in the hospital were significantly older (mean vs. years; ) and presented with profoundly lower mean hemoglobin ( vs. g/dL; ) and red blood cell counts ( vs. /µL; ). Hemodynamic markers were also severely compromised in the mortality group, characterized by substantially reduced left systolic blood pressure ( vs. mmHg; ), right systolic blood pressure ( vs. mmHg; ), and right diastolic blood pressure ( vs. mmHg; ). Methodological rigor is underscored by the inclusion of 95% confidence intervals and Cohen’s effect sizes, which reveal that left systolic BP, hemoglobin levels, and red blood cell counts exhibited the most pronounced magnitude of clinical effect between the two groups (Cohen’s ).

TABLE 3

VariableAlive (mean ± SD) (n = 300)Dead (mean ± SD) (n = 101)Mean diff. (95% CI)Cohen’s dp-value
Age (years)57.5 ± 14.862.2 ± 15.7−4.70 (−8.21 to −1.19)0.31 (small)0.007
Haemoglobin (g/dL)11.84 ± 2.029.46 ± 1.872.38 (1.93–2.83)1.20
(Very large)
<0.001
Red blood cells (/µL)4.16 ± 0.733.31 ± 0.720.85 (0.69–1.01)1.17
(Very large)
<0.001
Left systolic BP (mmHg)124.4 ± 20.689.4 ± 32.435.00 (29.21–40.79)1.45 (huge)<0.001
Right systolic BP (mmHg)124.2 ± 24.3108.8 ± 25.115.40 (9.92–20.88)0.63
(medium)
<0.001
Right diastolic BP (mmHg)71.8 ± 11.163.8 ± 15.88.00 (5.08–10.92)0.64
(medium)
<0.001

Comparison of baseline hemodynamic, hematological, and demographic parameters by patient discharge status.

Data are presented as Mean ± SD. Group sizes are for the “Alive” cohort and for the “Dead” cohort. Mean differences (Alive minus Dead) are presented with their corresponding 95% confidence intervals (CI). Cohen’s was calculated using pooled standard deviations to estimate effect size magnitude (; ; ; ). Two-sample independent -tests were applied for all between-group comparisons.

3.4 Comparison of non-normally distributed clinical and laboratory parameters by discharge status

Non-survivors experienced significantly longer hospital stays and were characterized by elevated baseline levels of aspartate aminotransferase (AST), alanine aminotransferase (ALT), creatinine, white blood cell (WBC) count, and international normalized ratio (INR). Conversely, patients who died exhibited significantly lower levels of albumin, total protein, oxygen saturation (), and diastolic blood pressure. No significant differences in body temperature or respiratory rate were observed between the groups (Table 4; Figure 2).

TABLE 4

VariableAlive median
[Q1–Q3]
Dead median
[Q1–Q3]
p value
Hospital stay (days)1 [0–4]5 [1, 1–4]<0.001
AST (U/L)21 [15–30]54 [24–205]<0.001
ALT (U/L)19 [13–31]28 [13–74]<0.001
Total protein (g/L)a68 [63–72]62 [53–69]0.001
Albumin (g/L)40 [36–43]25 [22–30]<0.001
Creatinine (µmol/L)64.5 [53–79]88.5 [49–170]<0.001
WBC (×109/L)6.4 [4.8–8.5]10.5 [6.3–16.1]<0.001
INR1.1 [1.1–1.2]1.6 [1.3–2.0]<0.001
Oxygen saturation (%)b97 [95–98]91 [85–95]<0.001

Comparison of non-normally distributed clinical and laboratory parameters by discharge status.

Data are presented as median [Q1–Q3] to reflect the interquartile range (25th and 75th percentiles). Base sample sizes are for the “Alive” cohort and for the “Dead” cohort, unless otherwise indicated. Between-group comparisons were evaluated using the non-parametric Mann–Whitney U test.

a

Total protein values were available for a subset of 146 alive and 65 dead patients.

b

Oxygen saturation values were available for 252 alive and 59 dead patients (reflecting a 22.4% missingness rate in the full cohort).

FIGURE 2

3.5 Comparison of laboratory parameters and outcomes between breast and colorectal cancer patients

Colorectal cancer patients had significantly higher AST, white blood cell count, serum creatinine, and INR compared with breast cancer patients, and significantly lower serum albumin. ALT did not differ between the two groups, with identical medians observed in both. In-hospital mortality was substantially and significantly higher in colorectal cancer patients than in breast cancer patients. Breast cancer patients were younger and had better haemoglobin and lower creatinine, consistent with a less advanced systemic disease burden at the time of admission. The differences in creatinine and albumin between the two cancer groups are clinically meaningful: colorectal cancer is more likely to involve hepatic metastases and peritoneal disease, both of which impair albumin synthesis and increase renal load, whereas breast cancer patients in this cohort had predominantly non-metastatic disease at lower rates of visceral involvement (Table 5).

TABLE 5

VariableBreast cancer
Median [Q1–Q3]
(n = 134)
Colorectal cancer
Median [Q1–Q3]
(n = 95)
p-value
Age (years)a54.4 ± 11.760.9 ± 14.6<0.001
Haemoglobin (g/dL)12.1 [10.9–13.2]10.9 [9.2–12.8]0.002
WBC (×109/L)6.1 [4.7–8.8]7.5 [5.2–10.9]0.010
RBC (×106/µL)4.2 [3.8–4.7]3.9 [3.4–4.5]0.060
INR1.1 [1.0–1.2]1.2 [1.1–1.5]<0.001
AST (U/L)23.0 [16.0–31.0]26.0 [19.0–49.5]0.028
ALT (U/L)21.0 [14.0–33.2]21.0 [13.0–43.0]0.944
Total protein (g/L)69.0 [65.0–73.0]64.0 [57.1–69.0]<0.001
Serum albumin (g/L)41.0 [37.0–43.0]38.0 [29.0–42.0]0.003

Comparison of laboratory parameters and outcomes between breast and colorectal cancer patients.

Data are presented as median [Q1–Q3] to reflect the interquartile range (25th and 75th percentiles) unless otherwise specified. Continuous non-normally distributed variables were compared using the Mann–Whitney U test.

a

Age data are presented as Mean ± SD and were analyzed using an independent samples -test.

3.6 Drug-specific and cancer-specific comparisons

Electronic medical records were comprehensively reviewed to extract data on all prescribed medications, which were subsequently classified into non-oncological (supportive care), oncological (cancer-directed therapies), and other unclassified therapeutic agents. A total of 2,148 prescriptions were analyzed in the clinical cohort. Non-oncological supportive care medications represented the majority of the therapeutic burden, accounting for 59.03% () of the total prescribed agents. This high proportion reflects the intensive supportive care management—including analgesics, antiemetics, and nutritional supplements—required alongside standard oncology regimens. Cancer-directed oncological therapies constituted 40.74% () of the total prescriptions. Within the oncological drug classification, traditional cytotoxic chemotherapy was the predominant treatment strategy, representing 27.98% () of overall prescriptions and 68.69% of all oncological agents. This was followed by targeted therapies at 6.84% (), hormonal therapies at 3.59% (), immunotherapies at 2.10% (), and radiotherapies at 0.23% (). A minor fraction of the dataset consisted of unmapped or miscellaneous placeholder entries, categorized as others/unspecified, which accounted for 0.23% () of the aggregate prescription data (Table 6).

TABLE 6

Drug classificationTreatment type/categoryTotal prescriptionsOverall proportion (%)
Non-oncologicalSupportive care1,26859.03%
OncologicalChemotherapy60127.98%
Targeted therapy1476.84%
Hormonal therapy773.59%
Immunotherapy452.10%
Radiotherapy50.23%
Subtotal (oncological)All cancer therapies87540.74%
OtherOthers/unspecified50.23%
TotalAll prescriptions2,148100.00%

Distribution and proportional breakdown of prescribed oncological and non-oncological medications across the clinical cohort (N = 2,148).

Data represent the absolute number and percentage distribution of prescriptions written across the entire patient cohort (). Percentages are calculated relative to the grand total of all categories combined.

3.6.1 Drug-drug interaction profile

A total of 2,148 drug prescriptions were evaluated across 401 patients, with a median of five drugs per patient [Q1 = 3, Q3 = 7]. Polypharmacy, defined as the concurrent use of five or more medications, was present in 206 patients (51.4%). Of the total prescription records, potential drug–drug interactions (pDDIs) were flagged in 2,117 instances (98.6%), leaving only 31 records (1.4%) completely devoid of identified interactions (Table 7; Figure 3A).

TABLE 7

DDI severity screening scoreRecord count (n = 2,148)% of total recordsPatient count (n = 401)% of total patients
Score 0 — No potential interaction identified311.40%143.50%
Score 1 — moderate (monitor or adjust dose)93543.50%——
Score ≥2 — severe (avoid or use extreme caution)1,17954.90%38796.50%
Score 3 — most severe (contraindicated)30.10%——
Total records/patients with any flagged pDDI2,11798.60%38796.50%

DDI Severity Classification (per drug prescription record; N = 2,148 total records).

pDDI, potential drug-drug interaction. Percentages may not total exactly 100% due to rounding. Patient count totals omit scores where multiple interactions per patient span different categories, representing baseline distribution.

FIGURE 3

Among these flagged interactions, severe risk profiles predominated: 54.9% () were categorized as severe (severity score = 2) and 0.1% () carried the most severe, contraindicated classification (severity score = 3). Moderate risks (severity score = 1) accounted for 43.5% (). Specific named pDDI pairs, where the screening tool identified a distinct interacting partner by name, were recorded in 529 instances; of these, 503 (95.1%) fell within the severe classification tier. At the patient level, 387 out of 401 individuals (96.5%) were prescribed at least one medication carrying a screening flag for a severe potential interaction (Table 7).

Stratification by mechanistic classification (Table 8; Figure 3B) revealed that potential metabolic interactions were the most frequent tier (894 records, 41.6%), reflecting predicted cytochrome P450 (CYP)-mediated competition between antineoplastic agents and co-prescribed antimicrobials, corticosteroids, or antiemetics. Potential toxicity overlap was the second most frequent category (395 records, 18.4%), followed by screening flags for potential immune-related toxicity (147 records, 6.8%), potential central nervous system (CNS) depression (118 records, 5.5%), potential nephrotoxicity (115 records, 5.4%), and potential neuropathy (106 records, 4.9%).

TABLE 8

Interaction typen%
Metabolic interaction89441.6%
Toxicity overlap39518.4%
Immune-related toxicity1476.8%
CNS depression1185.5%
Nephrotoxicity1155.4%
Neuropathy1064.9%
Cardiotoxicity874.1%
Hepatotoxicity854.0%
Thrombosis prevention/risk743.4%
Other (neurotoxicity, GI, pulmonary, other)351.6%

Interaction type distribution (2,148 records).

3.6.2 Statistical evaluation of organ-specific toxicity exposures and subgroup prevalence

To assess whether the pharmacological profile of prescribed regimens was statistically associated with clinical outcomes, exposure to medications carrying screening flags for potential organ-specific toxicities was evaluated by survival cohort (Table 9). Potential hepatotoxic medication exposure was highly prevalent but did not vary significantly between groups (82.7% of survivors vs. 73.3% of non-survivors; ). Notably, potential neuropathic medication exposure was significantly more prevalent among survivors than non-survivors (66.7% vs. 50.5%; ). Rather than reflecting a pharmacologically protective effect, this highlights a survivorship pattern: taxane-based regimens, which drive these neuropathic flags, were disproportionately administered to stable patients who completed treatment, whereas critically ill patients were less likely to receive full cytotoxic cycles before clinical deterioration. Exposure to medications with potential thrombosis risks (29.0% vs. 27.7%; ) and nephrotoxic drug exposure (20.3% vs. 23.8%; ) did not differ significantly between groups, while potential cardiotoxic exposure remained rare across the entire cohort (0.5% overall; ).

TABLE 9

Organ categoryn patients (%)Alive n (%)Dead n (%)
Hepatotoxic drug exposure322 (80.3%)248 (82.7%)74 (73.3%) p = 0.056
Neuropathic drug exposure251 (62.6%)200 (66.7%)51 (50.5%) p = 0.005 *
Thrombosis-risk drug exposure115 (28.7%)87 (29.0%)28 (27.7%) p = 0.906
Nephrotoxic drug exposure85 (21.2%)61 (20.3%)24 (23.8%) p = 0.556
Cardiotoxic drug exposure2 (0.5%)1 (0.3%)1 (1.0%) p = 1.000

Organ-Specific Drug Exposure by Discharge Outcome (patient level).

* Indicates statistical significance at p < 0.05.

Collectively, these screening data indicate that pre-existing patient-specific organ dysfunction—such as baseline renal disease ()—rather than database-flagged interaction categories per se, acts as the dominant determinant of mortality risk. Descriptively, patients with potential nephrotoxic medication exposure in the absence of renal disease had a mortality rate of 23.5% (16 of 68), tracking closely with the 21.1% mortality rate (59 of 279) of those with neither risk factor. Conversely, patients presenting with pre-existing renal disease alone experienced a mortality rate of 48.6% (18 of 37), which converged with the 47.1% mortality rate (8 of 17) observed in patients with both risk factors, underscoring baseline renal impairment as the primary risk driver.

Patient-level screening for organ-specific potential toxicities (Figure 3C) confirmed that potential hepatotoxicity flags were the most widespread (80.3%), followed by flags for potential neuropathy (62.6%), potential thrombosis risks (28.7%), potential nephrotoxicity (21.2%), and potential cardiotoxicity (0.5%).

When evaluating severe interactions across specific oncological classifications (Table 10; Figure 3D), near-universal exposure to at least one high-risk flag () was observed across all cohorts, yielding an overall patient-level baseline penetration rate of 96.5% (). The highest risk profile was observed in the hematologic malignancy group, which demonstrated an absolute 100.0% penetration rate (). The severe interaction rates remained profoundly elevated in both the breast cancer and colorectal cancer subgroups at 98.50% () and 97.90% (), respectively. The lowest comparative rates were found in the lung cancer subgroup (92.90%) and the heterogeneous “Other” malignancy classification (93.30%), though both subsets still closely mirrored the pervasive interaction burden observed across the wider study population.

TABLE 10

Cancer groupPatient subgroup (n)Patients with ≥1 potential severe interaction flag (n)aSubgroup severe pDDI rate (%)Reference cohort rate (%)
Breast13513398.50%96.50%
Colorectal959397.90%96.50%
Lung131292.90%96.50%
Hematologic2323100.00%96.50%
Other13512693.30%96.50%
Overall cohort40138796.50%—

Stratification of potential organ-specific toxicities and patient-level high-risk interaction rates by malignancy category.

pDDI = potential drug–drug interaction. a ‐ Potential severe drug–drug interactions were defined as major, contraindicated, or high-risk interactions requiring clinical intervention, intensive monitoring, or alternative therapy, identified using standard interaction screening software.

3.6.3 Comparative profile of paclitaxel and oxaliplatin treatment sub-cohorts

Among patients receiving paclitaxel-based therapy () and those receiving oxaliplatin-based regimens (), no significant differences were found in baseline liver enzymes (AST, ALT), serum albumin, white blood cell count (WBC), or overall length of hospital stay. However, serum creatinine was significantly lower in the paclitaxel cohort compared to the oxaliplatin cohort ( [52–75] vs. [66–91], ), reflecting the known nephrotoxic profile and platinum-class effect associated with oxaliplatin. The international normalized ratio (INR) was also significantly lower in the paclitaxel group ( [1.1–1.2] vs. [1.1–1.6], ), suggesting a greater underlying coagulopathy burden in patients receiving oxaliplatin, which may be secondary to a higher degree of hepatic involvement or tumor burden in their primary colorectal or gastric malignancies. Overall, in-hospital outcomes did not differ significantly between the two therapeutic groups (Table 11).

TABLE 11

VariablePaclitaxel median
[Q1, Q3](n = 85)
Oxaliplatin median
[Q1, Q3](n = 81)
p-value
AST (U/L)25 [17, 44]24 [19, 41]0.874
ALT (U/L)25 [16, 37]20 [14, 27]0.247
Albumin (g/L)39 [34, 43]38 [28, 43]0.383
Creatinine (μmol/L)61 [52, 75]72 [66, 91]0.007
WBC (×109/L)6.4 [4.7, 10.4]7.6 [6.0, 10.5]0.167
INR1.1 [1.1, 1.2]1.2 [1.1, 1.6]0.017
Hospital stay (days)1 [0, 4]1 [0, 5]0.284

Comparison of laboratory parameters and outcomes between paclitaxel and oxaliplatin.

Data are presented as Median [Q1, Q3]. -values were calculated using the Mann–Whitney U test for continuous variables. A -value indicates statistical significance.

3.7 Potential drug–drug interaction profile

Across the final analytical cohort ( patients) with a total of active inpatient prescriptions evaluated, potential drug–drug interactions (pDDIs) were extensively characterized across mechanistic, severity, organ-specific, and oncological domains (Table 12). Mechanistically, Pharmacokinetic Pathways driven by metabolic interactions (e.g., Cytochrome P450 [] pathways) represented the primary driver of interactions, accounting for () of the total identified pDDIs, while Pharmacodynamic Pathways characterized by toxicity overlaps constituted a substantial secondary driver, comprising general toxicity overlap (; ), nephrotoxicity (; ), neuropathy (; ), and cardiotoxicity (, ). When evaluated by clinical severity per individual interaction, the vast majority were classified as severe (score ), comprising () of the profile, followed by moderate interactions (score 1) at () and minor interactions (score 0) at (). Regarding specific antineoplastic protocols, paclitaxel (combined regimens) was the most frequent driver of pDDIs, accounting for () of the cohort’s profile, followed closely by the doxorubicin/AC protocol (, ), docetaxel (, ), oxaliplatin (, ), and irinotecan (, ). At the individual patient level (), organ-specific exposures carrying DDI risks were highly prevalent; a total of of patients () were exposed to at least one hepatotoxic DDI during their stay, while () experienced an exposure carrying a risk of induced neuropathy, () were exposed to a thrombosis-risk interaction, and () faced a nephrotoxic DDI exposure. Finally, stratification of severe pDDIs across primary tumor classifications revealed distinct variations among the clinical subgroups, with the highest prevalence of severe interactions observed within the colorectal cancer subgroup (, of ), followed by the gynecologic cancer group (, of ) and the breast cancer cohort (, of ) (Table 12).

TABLE 12

DomainCategoryn%
Most frequent drugs in DDIsPaclitaxel (combined)266.5
Docetaxel133.2
Doxorubicin/AC protocol256.2
Oxaliplatin71.7
Irinotecan61.5
Mechanistic classification (per interaction)
Pharmacokinetic pathwaysMetabolic interactions (e.g., CYP3A4 pathways)89441.6
Pharmacodynamic pathwaysGeneral toxicity overlap39518.4
Nephrotoxicity1155.4
Cardiotoxicity874.1
Neuropathy1064.9
DDI severity (per interaction)Severe (score ≥2)34485.8
Moderate (score 1)4411.0
Minor (score 0)133.2
Organ-specific DDI exposure (per patient)Any hepatotoxic DDI33583.5
Any nephrotoxic DDI8521.2
Any neuropathy DDI25162.6
Any thrombosis-risk DDI12230.4
Severe DDIs by cancer groupBreast cancer (n = 135)11383.7
Colorectal cancer (n = 95)8791.6
Gynecologic cancer (n = 26)2388.5

Drug–drug interaction profile of the study cohort (N = 401, total drugs in prescriptions N = 2,148).

3.8 Correlation between clinical parameters and in-hospital mortality

Correlation analyses demonstrated statistically significant relationships with in-hospital mortality (all ; Table 13 and Figure 4). Moderate-to-strong negative correlations were found for baseline albumin (), left systolic blood pressure (), and left diastolic blood pressure (), alongside a moderate negative correlation with oxygen saturation (). Conversely, positive correlations tracked with increased mortality risk: moderate relationships were noted for international normalized ratio (INR; ) and serum creatinine (), while weak-to-moderate positive associations were identified for serum aspartate aminotransferase (AST; ), alanine aminotransferase (ALT; ), white blood cell count (), mechanical ventilation requirements (), and radiological disease extent ().

TABLE 13

VariableVariable typeAppropriate statistical methodCorrelation/metric (rpb or ϕ)Interpretationp-value
AlbuminContinuousPoint-biserial (rpb)−0.657Moderate-to-strong negative: higher baseline albumin is strongly associated with survivorship (protective factor)<0.001
Left systolic BPContinuousPoint-biserial (rpb)−0.533Moderate-to-strong negative: higher systolic blood pressure is associated with survivorship (hypotension tracks with mortality)<0.001
Left diastolic BPContinuousPoint-biserial (rpb)−0.527Moderate-to-strong negative: higher diastolic blood pressure is associated with survivorship (hypotension tracks with mortality)<0.001
Oxygen saturationContinuousPoint-biserial (rpb)−0.399Moderate negative: lower baseline SpO2 is significantly associated with increased mortality<0.001
INRContinuousPoint-biserial (rpb)0.418Moderate positive: coagulopathy (higher INR) is moderately associated with increased mortality risk<0.001
CreatinineContinuousPoint-biserial (rpb)0.337Moderate positive: renal impairment (higher serum creatinine) is moderately associated with increased mortality risk<0.001
ASTContinuousPoint-biserial (rpb)0.245Weak-to-moderate positive: hepatic injury indicator (elevated AST) tracks with higher mortality<0.001
ALTContinuousPoint-biserial (rpb​)0.234Weak-to-moderate positive: hepatic injury indicator (elevated ALT) tracks with higher mortality<0.001
WBC countContinuousPoint-biserial (rpb)0.201Weak-to-moderate positive: leukocytosis (systemic inflammation/infection) tracks with higher mortality<0.001
Mechanical ventilationBinary (0/1)Phi coefficient (ϕ)/Chi-square0.118Weak positive: requirement for mechanical ventilation is associated with increased mortality0.018
Disease extentOrdinal/continuousPoint-biserial (rpb) or spearman ρ0.102Weak positive: greater radiological extent of disease weakly tracks with increased mortality0.042

Point-biserial and phi correlation coefficients between clinical parameters and in-hospital mortality.

In-hospital mortality is recorded as a binary outcome (; ). Bivariate correlations between clinical parameters and mortality were calculated using point-biserial correlation () for continuous metrics and the phi coefficient () for dichotomous categorical metrics. Both yield numerical values mathematically equivalent to Pearson’s under binary coding.

FIGURE 4

3.9 Multivariable logistic regression analysis

Multivariable binary logistic regression identified four independent physiological predictors of in-hospital mortality: renal failure, baseline serum albumin, serum creatinine, and oxygen saturation. The multivariable model demonstrated exceptional discrimination and overall classification performance (; Table 14). Overall in-hospital mortality in this cohort reflects a complex interplay between underlying oncological disease severity and acute physiological decompensation. Across non-survivors, hemodynamic instability, anemia, acute renal dysfunction, and systemic inflammatory activation followed a consistent, predictable trend. Taken together, these findings indicate that overall physiological maladaptation and multi-organ strain—rather than isolated laboratory aberrations or database-flagged drug interactions alone—characterize the patients at the highest risk of in-hospital mortality (Table 14; Figures 4, 5).

TABLE 14

VariableBSEWaldAdjusted ORp value
Renal failure3.791.903.9944.290.046
Albumin−0.350.137.140.710.008
Creatinine0.020.014.221.020.040
Oxygen saturation−0.390.166.250.680.012

Multivariable logistic regression analysis for predictors of in-hospital mortality.

Treatment interventions (e.g., mechanical ventilation) and highly collinear redundant markers were excluded from the final multivariable binary logistic regression model to prevent statistical suppression, overadjustment, and severe multicollinearity with primary organ dysfunction variables (e.g., renal failure and low ).

FIGURE 5

3.10 Polypharmacy, interaction effects, and mortality

Polypharmacy was evident in about half of the cohort, and it did not solely correlate with mortality. But drug-drug interaction effects showed a statistically significant association with mortality. Moreover, even though the level of polypharmacy was very high, it was not associated with mortality when it was adjusted with confounders. This observation suggests that the only number of drugs might not be significant but the clinical circumstances surrounding their prescription Table 15.

TABLE 15

Polypharmacy/Variablen (%)Alive nDead nORp value
Polypharmacy status
No195 (48.6)14154——
Yes206 (51.4)160460.7620.244
Interaction-related variables
Interaction effect———0.9940.015
Interaction type———1.0120.597

Polypharmacy status, association with in-hospital mortality, and regression analysis.

Dashes (—) indicate variables or cells for which the measure was not applicable or not estimated in the corresponding analysis (e.g., descriptive or univariate measures not included in regression models.

3.11 Receiver operating characteristics analysis

ROC curve analysis confirmed serum albumin as the strongest single predictor of in-hospital mortality, with an of . At the Youden-optimal cut-off of , albumin yielded a sensitivity of and specificity of , consistent with the established clinical threshold for hypoalbuminemia. INR was the second most discriminating individual variable (): at an optimal cut-off of , sensitivity was and specificity was , reflecting the clinical utility of coagulopathy as a marker of hepatic dysfunction and systemic failure. Oxygen saturation demonstrated fair discrimination (); its Youden-optimal cut-off of yielded a sensitivity of and specificity of . This estimate should be interpreted with caution given the overall missingness in oxygen saturation data, which was concentrated among acutely unwell patients admitted in extremis. Serum creatinine showed more modest individual discrimination (): at the optimal cut-off of , specificity was high at , but sensitivity was limited to . This indicates that while marked creatinine elevation () is highly specific for mortality risk, a significant proportion of non-survivors succumb without reaching extreme renal biomarker thresholds. Individual ROC curves are displayed in Figures 6A,B.

FIGURE 6

The combined multivariable logistic regression model—incorporating renal failure, serum albumin, serum creatinine, and oxygen saturation—achieved an overall of () in patients with complete data across all predictor variables. The Youden-optimal operating point identified a sensitivity of and specificity of , corresponding to a Youden statistic of . This level of discrimination substantially exceeds that of any individual biomarker alone and aligns with published combined risk prediction models in hospitalized oncology populations, where ensemble approaches routinely report values between and (). The combined model AUC reflects apparent rather than cross-validated performance, as the model was derived and evaluated in the same dataset. External validation in a prospective multicenter cohort is required before clinical use.

3.12 Summary of key results

In-hospital mortality within this cohort was primarily driven by acute physiological derangements and multi-organ dysfunction rather than baseline demographics or specific oncological diagnoses. The strongest, most consistent independent predictors of in-hospital mortality included acute renal dysfunction (manifested by elevated serum creatinine and clinical renal failure), severe hypoalbuminemia, and hypoxemia (reduced oxygen saturation), supported by unadjusted bivariate associations with hemodynamic instability and mechanical ventilation requirements. Conversely, baseline demographic factors and primary cancer types alone did not significantly impact discharge status, underscoring that short-term mortality is dictated by acute physiological decline rather than underlying malignant classification.

4 Discussion

The clinical management of patients undergoing oncological treatment involves a delicate balance between administering highly potent antineoplastic regimens and executing robust supportive care to mitigate systemic toxicities and tumor-related symptoms. Analysis of the 2,148 prescriptions within this cohort reveals that non-oncological supportive medications comprised the clear majority of the therapeutic burden (59.03%, ), emphasizing the intensive polypharmacy environment required to sustain oncology patients through active treatment. Conversely, cancer-directed therapies accounted for 40.74% () of the total prescriptions, with traditional cytotoxic chemotherapy serving as the predominant treatment strategy (27.98% of overall prescriptions). This extensive multi-drug landscape underscores an unavoidable exposure to complex drug-drug interactions (DDIs), which demand precise risk stratification to protect patient safety without compromising oncological efficacy. Crucially, the demographic features and cancer type alone were not the main factors contributing to in-hospital deaths in this retrospective study. Instead, acute physiological derangements and multi-organ dysfunction predominated in driving the mortality rate of this multinational oncology cohort treated in a tertiary care centre in Saudi Arabia. The most powerful independent predictors of in-hospital mortality proved to be renal failure, hypoalbuminemia, elevated serum creatinine, and hypoxemia (reduced oxygen saturation), complemented by significant unadjusted associations with hemodynamic instability and the necessity for mechanical ventilation. These findings highlight that systemic clinical decompensation and acute illness severity—rather than baseline cancer characteristics or database-flagged interactions alone—are the focal determinants of short-term outcomes in hospitalized cancer patients (; ).

Among these, renal dysfunction emerged as an exceptionally powerful independent predictor of mortality. Loss of renal function serves as a critical marker of both pre-existing comorbidity and treatment-induced pathology. For example, platinum-based compounds are highly nephrotoxic and frequently induce acute tubular damage, cellular necrosis, and severe electrolyte imbalances. Furthermore, kidney dysfunction fundamentally alters the pharmacokinetics and clearance of antineoplastic agents, consequently exposing the system to a greater free fraction of cytotoxic substances and compounding systemic toxicity. This bidirectional, hazardous interaction between deteriorating renal function and targeted pharmacotherapy underscores its strong correlation with increased in-hospital mortality across multiple malignancies, including lung, breast, and colorectal cancers (; ).

Similarly, hypoalbuminemia demonstrated a significant independent relationship with in-hospital mortality. As a classic surrogate marker of malnutrition, catabolic stress, and chronic systemic inflammation, reduced serum albumin concentrations fundamentally alter drug binding kinetics (; ). A drop in albumin increases the circulating free drug fraction, leading to unpredictable pharmacokinetic shifts that predispose patients to severe, unmanaged toxicities. This underscores the predictive utility of routine nutritional and biochemical profiling, reinforcing the necessity of preemptive nutritional optimization in acute care settings.

Indicators of hemodynamic and physical instability, such as hypotension, hypoxemia, and elevated inflammatory markers, were also tightly coupled with adverse outcomes (). While mechanical ventilation was a dominant predictor of mortality, its clinical value should be interpreted as an indicator of critical illness severity rather than an isolated therapeutic risk factor. Advanced respiratory failure, fulminant sepsis, and multiorgan breakdown characterize the patients requiring ventilatory support in this cohort (). While polypharmacy was highly prevalent due to local prescriptive trends and a heavy baseline burden of infections, it was not an independent predictor of death on its own. Instead, the net presence of severe drug-drug interactions underpinned the statistical correlation with mortality, bringing medication safety to the forefront of oncology care (; ).

Beyond traditional drug-drug interactions, underlying pharmacogenetic variability plays a critical role in modulating systemic drug exposure and idiosyncratic toxicity profiles within this cohort. For instance, polymorphisms in the DPYD gene, which encodes the rate-limiting enzyme dihydropyrimidine dehydrogenase, severely impair the catabolism of fluorouracil (5-FU), predisposing carriers to life-threatening myelosuppression and severe mucositis even under standard dosing regimens (). Similarly, genetic variations in UGT1A1 (such as the UGT1A1*28 allele) reduce the glucuronidation of SN-38, the active metabolite of irinotecan, driving up the risk of severe, dose-limiting neutropenia and debilitating diarrhea ().

This risk landscape is further complicated by the widespread, often unguided, concomitant use of natural products and herbal supplements among oncology patients. Many common botanical supplements function as potent modulators of cytochrome P450 pathways (such as CYP3A4 inhibition or induction) or drug transporters like P-glycoprotein (). When taken alongside narrow-therapeutic-index agents such as paclitaxel or targeted tyrosine kinase inhibitors, these natural products can either precipitously drop systemic drug levels—leading to therapeutic failure—or dangerously elevate them, exacerbating organ-specific toxicities like the hepatotoxicity and nephrotoxicity tracked in this study ().

Collectively, the clinical implications of these overlapping variables emphasize the urgent transition toward personalized medicine in oncology. Retrospective stratification of clinical DDIs alone is no longer sufficient; true optimization of patient care requires a proactive framework. Integrating pre-therapeutic companion diagnostic screening—such as targeted DPYD and UGT1A1 genotyping—alongside comprehensive, pharmacist-led reconciliation of both prescriptive medications and over-the-counter herbal supplements is essential (). By customizing drug selection and preemptively adjusting dosages based on an individual’s unique genomic and polypharmacy profile, clinical teams can significantly mitigate the severe interaction risks (Scores 2) identified in this cohort, ultimately safeguarding therapeutic efficacy while minimizing preventable in-hospital morbidity and mortality.

In unadjusted univariate analyses, advanced age was significantly linked to mortality, a finding consistent with decreased physiological reserve, multimorbidity, and heightened vulnerability to aggressive oncologic regimens. Although age was not retained in the final multivariable model, its prognostic impact is likely mediated through age-related organ frailty and functional decline (; ; ). Similarly, male sex was associated with higher mortality rates in univariate models, mirroring broader cancer epidemiology trends that are often driven by differences in comorbidity loads, delayed healthcare-seeking behaviors, and disease severity at initial presentation rather than immutable biological differences (; ).

Notably, smoking status was significantly related to an increased risk of in-hospital mortality, reinforcing extensive data showing that active smoking impairs immune responses, augments cardiopulmonary complications, and blunts the efficacy of antineoplastic therapies (; ). Conversely, mounting evidence demonstrates that smoking cessation following a cancer diagnosis yields substantially better prognosis and discharge outcomes across various malignancies, particularly lung and head and neck tumors (; ). This strong association highlights the profound clinical significance of implementing structured, compulsory smoking cessation interventions within oncology care, particularly at the time of acute hospital admission when patient engagement with the healthcare infrastructure is maximally enhanced (; ).

The choice of Drugs.com as the primary DDI screening platform warrants direct comment in the context of published performance data. The platform integrates data from Micromedex, Cerner Multum, and the American Society of Health-System Pharmacists (ASHP), efficiently providing a composite output from three independent, validated sources through a single query interface. In a head-to-head comparison of nine DDI screening tools evaluating oral oncology drug pairs, Drugs.com achieved the highest performance accuracy score among freely available tools (0.83), demonstrating no statistically significant difference from the top-performing subscription database, Lexicomp® (0.89) ().

A subsequent evaluation across 229 oral antineoplastic drug pairs confirmed that Drugs.com (accuracy score 352 of 400) and Lexicomp® (355 of 400) both substantially outperformed other major tools, including Micromedex (270 of 400) and Medscape (280 of 400) (). These established benchmarks support its application as a highly credible, accurate screening platform in retrospective oncology research settings where access to cost-prohibitive subscription software may be restricted. For pharmacokinetic interactions—such as metabolic alterations driven by paclitaxel and docetaxel regimens—simple administration-time separation is ineffective due to the prolonged, systemic modulation of hepatic enzyme activity (). Management requires proactive dose modifications of narrow-therapeutic-index supportive care agents or transitioning to non-interacting alternatives (). Conversely, pharmacodynamic pDDIs characterized by cumulative systemic overlaps demand vigilant clinical and laboratory surveillance rather than scheduling adjustments. Specifically, these toxicities necessitate regular complete blood count monitoring to mitigate overlapping myelosuppression, routine serum creatinine tracking to manage the patient-level nephrotoxic exposure rate, and serial electrocardiographic tracking to actively screen for additive corrected QT () interval prolongation ().

Several limitations of the present study must be acknowledged. First, the primary outcome measure was restricted to in-hospital mortality. Valuable intermediate and longitudinal outcomes, such as intensive care unit (ICU) mortality metrics, post-discharge mortality proportions, and exact length of stay parameters, were not evaluated and remain important avenues for future prospective studies. Second, the retrospective, single-center design naturally limits direct causal inferences and immediate geographic generalizability, as findings may partially reflect localized prescription trends, specific supportive care capabilities, or regional infection burdens. Furthermore, the absence of consistently documented functional performance status scores (e.g., ECOG or Karnofsky) and definitive treatment intent (curative vs. palliative) creates a possibility for residual confounding.

Additionally, while the sample size of 401 patients was sufficient for the primary multivariable analysis, it limited the statistical power required for robust subgroup comparisons. The combined model’s excellent area under the receiver operating characteristic curve (AUC = 0.953) represents an apparent rather than cross-validated discrimination metric, as the model was derived and evaluated within the same complete-case dataset ( patients with complete data across all five primary predictors). Finally, interactions were classified using a single DDI database and were not independently cross-referenced against a secondary database at the exact point of data collection. Future prospective multi-center studies utilizing dual-database verification protocols are recommended to further strengthen drug interaction severity classifications and enhance the generalizability of these findings.

5 Conclusion

The findings of this retrospective study demonstrate that short-term, in-hospital mortality among acute oncology inpatients is primarily driven by acute physiological derangements and organ failure rather than baseline oncological characteristics or demographic variables. Specifically, renal failure (), baseline hypoalbuminemia (), elevated serum creatinine (), and reduced oxygen saturation () were identified as robust, independent clinical predictors of in-hospital mortality. Among these, deterioration of the renal axis emerged as the most formidable prognostic indicator. Impaired kidney function not only marks systemic critical illness but directly compromises the clearance of narrow-therapeutic-index antineoplastic agents, exposing patients to prolonged cytotoxic drug exposure and exacerbating systemic treatment-related toxicity. Concurrently, profound hypoalbuminemia serves as a key surrogate for severe catabolic stress and systemic inflammation, fundamentally altering drug–protein binding kinetics and increasing the circulating free fraction of active drugs to dangerous thresholds.

From a pharmacological safety standpoint, this multinational cohort exposed a critical clinical vulnerability: a near-universal exposure to severe-class drug-drug interactions (DDIs), with 97.5% of the active patient population experiencing at least one severe (Score 2) interaction risk. Crucially, multi-variable modeling demonstrated that simple polypharmacy—defined merely by the sheer quantity of concurrent medications—does not independently predict in-hospital death. Instead, it is the qualitative composition of the drug combinations rather than the absolute drug count that constitutes the actual interaction-related clinical risk factor. Cornerstone oncology protocols, such as 5-fluorouracil (5-FU), oxaliplatin, and irinotecan regimens (e.g., FOLFOX, FOLFIRI, XELOX/CapOx), are frequently combined with intensive supportive care agents (analgesics, antiemetics, and proton-pump inhibitors), creating highly intricate metabolic competition and profound toxicity overlaps that challenge hepatic and renal reserves.

Consequently, these outcomes underscore the urgent necessity of shifting from reactive clinical management toward proactive, individualized risk mitigation. Structured, multidisciplinary clinical pharmacy reviews must be systematically integrated directly into oncology admission protocols. These specialized interventions should be double-pronged: Targeted Nephroprotective Surveillance: Pharmacists must actively screen for, flag, and modulate highly nephrotoxic drug combinations—particularly where platinum-based therapies or intense supportive regimens intersect—to protect renal clearance and stabilize baseline creatinine levels. Pre-Therapeutic Companion Genotyping: Given the high utilization of fluoropyrimidines and irinotecan in gastrointestinal and colorectal malignancies within this tertiary care setting, integrating companion pharmacogenetic screening is highly warranted. Preemptively establishing a patient’s DPYD and UGT1A1 polymorphic status allows for precise, genome-guided dose titrations before initial drug exposure. By combining automated DDI screening with genetic profiling and comprehensive over-the-counter herbal supplement reconciliation, clinical teams can safely deliver core antineoplastic regimens, significantly reduce preventable iatrogenic morbidity, and ultimately improve discharge outcomes and lower the in-hospital mortality proportion in hospitalized cancer patients.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by Hail health cluster, Research and Studies department, Institutional review board. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants; legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

AhA: Writing – original draft, Methodology, Writing – review and editing. MA: Writing – review and editing, Writing – original draft, Methodology. SMA: Writing – original draft, Writing – review and editing, Methodology. SFA: Methodology, Writing – original draft, Writing – review and editing. FA: Writing – review and editing, Methodology, Writing – original draft. AmA: Methodology, Writing – review and editing, Writing – original draft. HS: Software, Writing – review and editing, Writing – original draft, Formal Analysis. SaA: Software, Formal Analysis, Writing – review and editing, Writing – original draft. MFS: Writing – original draft, Formal Analysis, Software, Writing – review and editing. SiA: Funding acquisition, Project administration, Writing – original draft, Supervision, Conceptualization, Writing – review and editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This research is funded by Scientific Research Deanship at University of Ha’il- Saudi Arabia through project number RG-24 172.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The authors MFS, SA declared that they were an editorial board member of Frontiers at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declared that generative AI was used in the creation of this manuscript. The authors declare that Generative AI was utilized during the writing process of this manuscript to optimize linguistic clarity and adhere to journal formatting guidelines. The authors have independently verified all data, clinical conclusions, and references, and retain full accountability for the scientific integrity of the final text.

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/fphar.2026.1825631/full#supplementary-material

Abbreviations

INR, International normalized ratio

References

Summary

Keywords

hospitalized cancer patients, hypo-albuminemia, in-hospital mortality, laboratory biomarkers, medication safety, oncology, polypharmacy, potential drug–drug interactions

Citation

Alafnan A, Alorainan MM, Alanazi SM, Algharbi SF, Alshammari FR, Alghaslan A, Saleem H, Alghamdi S, Shaikh MF and Anwar S (2026) Clinical, laboratory, and medication-related factors associated with in-hospital mortality among oncology patients at a tertiary hospital in Saudi Arabia. Front. Pharmacol. 17:1825631. doi: 10.3389/fphar.2026.1825631

Received

08 March 2026

Revised

30 July 2026

Accepted

30 July 2026

Published

28 August 2026

Volume

17 - 2026

Edited by

Mohammed Salahudeen, University of Tasmania, Australia

Reviewed by

Valeria Conti, University of Salerno, Italy

Anna Pradiningsih, Universitas Muhammadiyah Mataram, Indonesia

Updates

Copyright

*Correspondence: Sirajudheen Anwar,

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics