Introduction

Spontaneous intracerebral hemorrhage (sICH) is the most common subtype of hemorrhagic stroke and is a critical disease, with only 12–39% of survivors achieving long-term functional independence1. Recent studies have demonstrated that the basal ganglia region exhibits the highest incidence of sICH. Its ___location often correlates with decreased levels of consciousness and motor deficits in patients with contralateral hemiplegia2. To reduce secondary injury caused by the physiological response of the hematoma, such as perilesional edema and inflammatory reactions from blood product toxicity3, surgical intervention for hematoma evacuation is often necessary to eliminate neurotoxic metabolites and reduce the mass effect of the hematoma4,5,6,7.

Unfortunately, several large randomized controlled trials (RCTs) have failed to demonstrate significant improvements in overall functional outcomes for surgical treatment compared to best medical care in basal ganglia intracerebral hemorrhage (ICH)4,8,9,10 This highlights the complexity of managing these cases and the need for careful patient selection. Recent clinical trials, however, have explored minimally invasive techniques for hematoma evacuation, offering some promising results. The MISTIE series demonstrated that minimally invasive surgery plus recombinant tissue plasminogen activator (rtPA) showed greater clot resolution than traditional medical management. The ICES trial described intraoperative CT-guided endoscopic surgery as safe and effective for removing acute intracerebral hematomas, with the potential to enhance neurological recovery. Most recently, the ENRICH trial presented that minimally invasive hematoma evacuation resulted in better functional outcomes at 180 days than guideline-based medical management in patients with lobar or anterior basal ganglia ICHs11.

However, the details of ___location-tailored surgical techniques, such as safe entry points, target trajectories, appropriate aspiration pressures, catheter diameters, and postoperative alteplase protocols, have not yet been fully established. In clinical settings, the Kocher’s point (KP) approach is widely used for inserting a catheter in patients with basal ganglia sICH owing to its established safety, feasibility, and accuracy, as supported by previous retrospective studies4,8,9,10,12,13,14 The KP approach was originally devised for ventriculostomy, which targets the foramen Monro of the lateral ventricle starting from the middle frontal gyrus15. Owing to the anatomical proximity between the lateral ventricle and basal ganglia, the KP approach has become the preferred method for performing stereotactic hematoma aspiration procedures for sICH located in the lentiform nucleus of the globus pallidus or putamen. Recently, supraorbital keyhole (SOK) surgeries using the anterior skull base approach have been widely used in the treatment of tumors or neurovascular diseases. Ventriculostomy via the basal frontal lobe has been attempted using the SOK approach, which demonstrated acceptable safety and feasibility16,17,18

Anatomically, the lentiform nucleus of the basal ganglia has an average length (anterior-to-posterior (AP)) of approximately 41–47 mm, which is far longer than its height (superior-to-inferior (SI), approximately 20 mm)19. We hypothesized that positioning a catheter with a longer diameter in the AP direction (SOK approach) would enable it to come in contact with the widest area of a typical hematoma, potentially resulting in a higher aspiration rate than that resulting by placing a catheter in the S-I direction (KP approach) (Fig. 1). Since other types of hematoma aspiration surgeries (tubular or endoscopic hematoma aspirations) frequently accompany the parenchyma injury along the trajectories, this long-axis concept was understood as an important tactic in these kinds of surgeries20. To verify this hypothesis, we conducted: (1) A retrospective analysis of clinical data using propensity score matching (PSM) to compare outcomes between the SOK and KP approaches; (2) Computational simulations to verify the theoretical advantages of the SOK route in hematoma aspiration.

Fig. 1
figure 1

Two trajectories for accessing the basal ganglia sICH; (A) Burr hole sites for two entries and consequent trajectories; (B) Axial CT scans for SOK approach (upper) and coronal CT scans for KP approach (lower). Blue: Kocher’s point (KP); orange: supraorbital keyhole (SOK).

Materials and methods

Patients

Between January 2010 and December 2021, data of 243 patients diagnosed with sICH who underwent surgery were retrospectively reviewed. The patients were classified as (1) those diagnosed with sICH in the “basal ganglia” and (2) those who underwent frameless stereotactic hematoma aspiration surgery alone using either the SOK or KP approach. Patients with a Glasgow coma scale (GCS) score of ≥ 5, a hematoma volume of ≥ 20 mL with neurological signs and symptoms, and no surgical contraindications including an underlying medical disease were included in the study. Meanwhile, patients with ICH caused by other factors (n = 42), those with ICH in other locations (n = 94), or those who underwent other types of surgical interventions (n = 31) were excluded. Finally, only 76 patients were selected. A flowchart of the study is shown in Fig. 2. The 76 reviewed patients were divided into the “typical” and “atypical” groups according to the shape of the hematoma. The typical type was defined when the length (A-P) of the hematoma exceeded its height (S-I) by at least 20% in a three-dimensional (3D) reconstructed image (Fig. 3). The current study was approved by the Institutional Review Board of the Human Research Center of Korea University Anam Hospital. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Due to retrospective nature of the study, informed consent was waived by the Institutional Review Board of the Human Research Center of Korea University Anam Hospital.

Fig. 2
figure 2

Flow chart of the patient selection process.

Fig. 3
figure 3

Three-dimensional reconstructed images of sICH: Representative patients with typical (upper) and atypical (lower) hematoma shapes. Hematomas with an anterior-posterior diameter of > 20% longer than the superior-inferior diameter were defined as “typical”.

Due to retrospective nature of the study Informed consent was waived by (the name of the IRB).

The patient’s general demographic and clinical data were retrospectively gathered, including age, sex, medical history (hypertension, diabetes, dyslipidemia, and previous stroke), medication history (use of antiplatelet or anticoagulant agents), onset to operation time, GCS score upon admission, the level of consciousness, and pupillary responses. To evaluate the outcomes, the length of intensive care unit stays, total length of hospital stays, Glasgow outcome scale scores at discharge, and modified Rankin Scale (mRS) scores at 3 and 6 months were also reviewed.

Radiologic evaluation

The source data of preoperative 2-mm thin-sliced brain computed tomography (CT) and postoperative 5-mm thin-sliced brain CT images were acquired in the Digital Imaging and Communications in Medicine format. Only intraparenchymal hematomas were manually identified pixel by pixel in each slice, after the setting thresholds of 40–100 Hounsfield units. Hematoma segmentation modeling was performed using the free open-source software, 3D Slicer (version 4.10.2, 3D Slicer, https://www.slicer.org). Two physicians independently segmented the hematoma, and the median hematoma volumes were used (Supplementary Fig. 1). Using the preoperative images, the following radiologic findings were evaluated: midline shifting, the side of the involved hemisphere, the presence of intraventricular hemorrhage, the absence of a basal cistern (ipsilateral), hematoma shape, and hematoma volume. Using the postoperative images, catheter-related complications, such as poor catheter positioning, hematoma expansion, and catheter-related bleeding (axis injury), were evaluated.

Operative technique

Prior to surgery, the stereotactic indicator was stitched to the forehead of the patient using a frameless technique. The stereotactic coordinates of the lesion were determined from a 2-mm thin-sliced brain CT scan using a neurosurgical navigation system (Stealth Station, Medtronic, Broomfield, CO, USA). Under general endotracheal anesthesia, the patient was placed in the supine position, with the head gently immobilized. The skin was carefully prepared with an aseptic betadine solution to prevent infection, followed by sterile draping. A skin incision was created either at the KP or supraorbital area. For the SOK approach, a small incision was made at the lateral half of the superior margin of the eyebrow, and a burr hole was created above the supraorbital ridge lateral to the supraorbital notch. A small dural and cortical incision was then made through bipolar cauterization, and a catheter was inserted to evacuate the hematoma following the calculated coordinates under the navigation guide. The catheter was initially inserted 5 mm into the deepest portion of the trajectory and then withdrawn. Manual aspiration with negative pressure was carefully performed during withdrawal using a 10-mL syringe. Hematoma aspiration was conducted under minimal pressure (< 2 mL in a syringe scale). The closed drainage system was connected to a catheter for external drainage and additional thrombolytic therapy. Thrombolytic therapy with recombinant tissue plasminogen activator (rtPA) was administered if the residual hematoma volume was more than 20 cc or 50% of the preoperative hematoma volume. A 1.0 mg/1 mL dose of rtPA with 2 cc of normal saline was injected via the catheter every 12 h. The catheter was clamped for 2 h and then released for natural drainage. After surgery, neurological complications such as supraorbital neuralgia or brow ptosis were evaluated based on clinical findings.

PSM and statistical analysis

After patient enrollment and classification into the KP (n = 28) and SOK (n = 22) groups, logistic regression analysis was conducted to estimate the correlation of each variable with the selected treatment. The coefficient estimates from this regression were subsequently used to retrospectively calculate the predicted probability (propensity score: ranging from 0 to 1) of assigning such treatment to each patient based on their characteristics. Each patient from the SOK group was then matched with another patient (1:1 matching) from the control group based on the closest probability of treatment and size of the available control group. Once each treatment patient was matched, the unused controls were removed, and the outcomes of the SOK group were compared with those of the newly matched control group (18 SOK patients vs. 18 KP patients)21. A flowchart of the patient classification and matched classification after PSM is provided in Fig. 2.

Continuous variables were expressed as the means and standard deviations or median and interquartile ranges, while categorical variables were expressed as numbers and percentages. Independent t-tests or Mann–Whitney U tests were used for analyzing continuous variables according to the normality of distributions, while chi-square tests were used for analyzing categorical variables. If more than 20% of the cells in the chi-square test had an expected frequency of less than 5, the likelihood ratio or Fisher’s exact test was used. A P value of < 0.05 was considered significant. Statistical analyses were performed using the standard software (version 23.0; SPSS, IBM, Chicago, IL, USA).

3D modeling for computational simulation

As depicted in Fig. 2, we selected eight representative patients after study enrollment and gathered their source data: four with typical sICH (two in the KP group and two in the SOK group) and four with atypical sICH (two in the KP group and two in the SOK group). Initially, the region corresponding to the hematoma was isolated from the rest by applying a threshold in the Slicer 3D program (Fig. 4A). Subsequently, the interior and exterior cavities of the hematoma were filled using the “filling holes” function to avoid unnecessary intricate surfaces and reduce simulation time (Fig. 4B). After setting the boundary conditions and removing the separate redundancies, a seamless one-body hematoma model was generated (Fig. 4C). The number of meshes was reduced using the re-mesh function in Autodesk Fusion 360 (Autodesk, Inc., San Francisco, CA, USA) while preserving the original hematoma shape and volume [Fig. 4D]. The simplified hematoma model had the same dimensions and volume as the original model. Regardless of the number of meshes (> 500), the simulation results exhibited a strong correlation.

Fig. 4
figure 4

Schematic of hematoma modeling. (A) Three-dimensional reconstruction of a thrombus from brain computed tomography data. (B) Surface simplification using the “filling holes” function of the simulator. (C) Elimination of superfluous regions. (D) Thrombus shapes with varying quantities of mesh and (E) their respective volumes. (F) Residues following suction with various mesh numbers.

Computational simulations

The hematoma aspiration simulation was performed using the COMSOL Multiphysics software (COMSOL, Inc., Burlington, MA, USA). The model was constructed using a two-phase solution comprising a hematoma (phase 1) and the brain (phase 2). The catheter tip had a hole with a boundary condition corresponding to the actual inner diameter of the catheter (2.5 mm) (Fig. 5). Pressure was applied to the boundary area between the circumference of the brain and catheter tip, and the Deformed Geometry function was used to simulate suction by retracting the catheter. The pressure boundary conditions were set at the boundary of the brain (Table 1). The catheter path was delineated at 5-mm intervals between the SOK or KP and the other sides of each condition. An image of the hematoma shape was generated using the isosurface function of the simulation software, and the remaining volume of the hematoma was calculated using the volume fraction function. Hematoma suction was simulated using the Navier–Stokes equation, wherein flow velocity (V), the density of blood plasma (ρ), pressure (p), and viscosity (µ) were considered. The physical values of these variations were obtained from previous studies (Table 1). 22,23,24

Fig. 5
figure 5

Schematic of computational simulation. (A) Two catheter trajectories (SOK and KP). (B) Simulated hematoma changes aspirated using a catheter (left SOK, right KP). KP Kocher’s point, SOK supraorbital keyhole.

Table 1 Physical values used for the simulation.

Results

Results of PSM

Initially, 50 patients who showed “typical basal ganglia sICH” were selected for PSM analysis, and 18 patients in each group were matched. The detailed results on PSM analysis were presented in Supplementary Fig. 4 and the logit propensity score of the matched group was 0.04 ± 0.0726 with three parameters (Level of consciousness, antiplatelet usage, and initial hematoma volume).

Results of clinical comparison analyses

The baseline demographic data of the 36 matched patients are presented in Supplementary Table 1. The patient’s mean age was 55 years, while the mean hematoma volume was 40.92 ± 22.80 3.

The results of the comparative analyses (KP vs. SOK) of the preoperative factors are shown in Table 2. After PSM, the preoperative factors showed no significant differences between the groups. Specifically, the preoperative hematoma volumes were 39.41 ± 20.11 and 42.42 ± 25.71 cm3, respectively.

Table 2 Results of the comparative analyses between the KP and SOK groups: preoperative factors.

The results of the main comparison analysis of the outcome parameters between the KP and SOK groups are presented in Table 3. In terms of radiologic outcomes, the incidence of catheter-related complications was lower in the SOK group (p = 0.074). In addition, a higher hematoma volume reduction was observed in the SOK group (13.018 ± 14.682 in KP vs. 21.693 ± 12.345 in SOK, p = 0.063). The reduction in hematoma percentage (= aspiration rate) was significantly lower in the SOK group (31.841 ± 40.131 in KP vs. 55.191 ± 25.387 in SOK, p = 0.045). In terms of clinical outcomes, the proportion of patients who achieved good functional outcomes (mRS score, 0–2 at 6 months) was significantly higher in the SOK group (27.8% vs. 61.1%, p = 0.044).

Table 3 Results of comparative analyses between the KP and SOK groups: outcome parameters.

Results of computational simulations

The results of the computational simulation of hematoma aspiration in the KP and SOK groups are depicted in Fig. 5B. Hematoma aspiration occurs in hematoma-mimicking phase 1 when negative pressure is applied to a catheter moving along a defined path (Fig. 6A). To simplify the simulation, a re-meshing operation was implemented to reduce the number of meshes and ultimately reduce computation time (Fig. 6B,C). To determine the maximum mesh size (minimum mesh number) required to sustain the hematoma structure and volume, the residual volume of the hematoma after suction was compared by adjusting the number of meshes (Fig. 6D). When the volumes of hematoma from eight patients were compared with the different numbers of meshes, a significant change was observed in the residual volume of patients with less than 400 meshes (Fig. 4E and F). Simulations were conducted using the actual hematoma data from eight patients divided into the typical (four patients) and atypical (four patients) groups. The simulation results demonstrate the suction value along the SOK path, particularly in patients with typical ICH. On the SOK path, the average residual volume after aspiration was 21.5%, which was significantly less than the average residual volume (43.4%) along the KP path (Fig. 6E,F). Supplementary Figs. 2 and 3 present the simulation outcomes of additional patients not shown in the main figure. A video of each simulation is also supplied in Supplementary Movies 1 and 2.

Fig. 6
figure 6

Results of computational simulation. (A) Schematic of computational hematoma aspiration analysis. (B) The morphology of the hematoma model, both typical and atypical, dependent on the mesh numbers. (C) Simulation time required based on the mesh numbers. (D) Residual volumes of hematoma after aspiration based on the mesh numbers. (E) 3D morphologies of residual hematoma after aspiration along one of the KP or SOK paths. (F) Quantified volumes of the residual hematoma. 3D three-dimensional, KP Kocher’s point, SOK supraorbital keyhole.

Discussion

The current study demonstrated that stereotactic hematoma aspiration via the SOK route is a safe and effective method for patients with basal ganglia ICH, especially in those with typical hematoma shapes. Computational simulation demonstrated a higher aspiration rate using the SOK route, and comparison analysis results revealed improved radiological and clinical outcomes in the SOK group.

Various studies have explored the advantages of open hematoma evacuation compared with those of medical treatment in patients with sICH12,13,14,25; however, serial RCTs revealed that early surgery had no overall benefits or had minimal relative advantages on hematomas in superficial locations26,27. Instead, several studies regarding minimally invasive surgery using frameless stereotactic hematoma aspiration techniques combined with thrombolytic therapy have been performed in patients with sICH and have revealed positive results. The MISTIE I trial revealed that minimally invasive surgery plus rtPA showed greater clot resolution than traditional medical management. Moreover, rtPA does not seem to exacerbate perilesional edema4. The MISTIE II trial showed that minimally invasive surgery plus rtPA appears safe with better functional outcomes at 180 days8. The ICES trial also described intraoperative CT-guided endoscopic surgery as a safe and effective method for removing acute intracerebral hematomas, with the potential to enhance neurological recovery10. A recent meta-analysis suggested the superiority of a minimally invasive strategy over conventional craniotomy, with a lower rate of rebleeding and a higher rate of good recovery28. Although the MISTIE III trial did not show demonstrate significant improvements in the proportion of patients achieving a good long-term response within 365 days after the occurrence of moderate to large ICH9, early and rapid hematoma evacuation “theoretically” improves cerebral perfusion, eliminates blood toxic products, minimizes oxidative stresses, and ultimately leads to better functional outcomes29,30. As most previous RCTs commonly presented the safety of the procedures4,8,9,10, it is now widely performed in most centers in selective patients with sICH.

In cases for minimally invasive surgery, the KP route is widely used for accessing the lateral ventricle or neighboring deep brain regions (basal ganglia or thalamus) owing to its guaranteed safety and convenience (familiarity with surgeons). To the best of our knowledge, most centers use the KP route as a routine trajectory for assessing basal ganglia sICH without considering the natural shape and anatomy of the basal ganglia. With the increasing adoption of minimally invasive surgeries for accessing the anterior skull base through SOK craniotomy, growing evidence supports the safety and feasibility of accessing the lateral ventricle through the basal frontal lobe16,17,18 We hypothesize that the SOK route is a safe and effective alternative for sICH evacuation, considering the AP-elongated shape of the basal ganglia19.

While there have been studies exploring the use of SOK approaches for various intracranial procedures, research specifically focusing on hematoma removal in basal ganglia sICH via this route is limited. It’s important to note that our study is unique in its specific enrollment of patients with deep-seated basal ganglia hematoma. For instance, Oh et al. (2022) investigated supraorbital endoscopic evacuation in patients with frontal traumatic ICH, which primarily involved subcortical hematomas31. This patient population differs significantly from our focus on deep-seated basal ganglia ICH. Similarly, Zheng et al. (2024) examined the utility of the supraorbital eyebrow arch keyhole approach for various skull base procedures, but did not specifically address catheter insertion for basal ganglia hematoma evacuation32. Our study, therefore, addresses a gap in the literature by specifically evaluating the efficacy of the supraorbital keyhole approach for stereotactic aspiration of typical basal ganglia ICH.

In terms of safety, eyebrow incision-related complications such as brow ptosis (associated with facial nerve injury) or frontal neuralgia (associated with supraorbital nerve injury) did not occur. This low incidence of complications is attributed to the small incision (2 cm) created on the lateral supraciliary line during a small burr hole craniotomy, which avoids the dangerous anatomy of the facial or supraorbital nerve (Fig. 1A). Furthermore, no significant differences were observed in the incidence of catheter-related complications between the SOK approach and the KP approach33. Every aspect of the procedure mirrored that of the KP approach, except for a slightly deeper engagement of the catheter by 1–2 cm through the white matter. Importantly, this difference did not affect the significant complication rate. When focusing on the AP-elongated ellipsoid type of hematoma, the total length from the bone margin to the hematoma is relatively the same in both approaches (Fig. 1B). Compared with traditional KP, the SOK approach enables a longer trajectory within the hematoma. This approach enables hematoma evacuation with just a single pass of the catheter, thus minimizing iatrogenic brain injury. This trajectory provides a multifocal aspiration core with only a single tract by reciprocating motion34. Moreover, considering the cerebral superficial venous system, the SOK approach has a lower likelihood of encountering midline superficial cortical veins compared with the KP approach35. Thus, the risk of postoperative venous infarction can be avoided.

In terms of efficacy, the aspiration rate achieved using the SOK approach was notably higher when the cohorts were compared. Specifically, the subgroup analysis revealed significantly improved radiological (aspiration rate) and clinical (favorable mRS score) outcomes in the SOK group, suggesting that patients with typical basal ganglia sICH may be suitable for the supraorbital approach. Our results indicate that 50 of 76 (65.8%) patients showed a typical AP-elongated ellipsoid pattern, allowing for the insertion of a catheter with a longer diameter via the SOK route. This suggests that the majority of patients with basal ganglia sICH can be assessed using the simple eyebrow approach. Although the efficacy of the SOK route was sufficiently verified using clinical data, it has inherent limitations, including its retrospective design; small sample size; and biases such as surgeon preference, patient selection, and differences in baseline characteristics. To address these limitations, we conducted a computational simulation to establish foundational qualifications and provide academic verification.

Computational simulations confirmed that aspiration via the SOK route generated less residual hematoma in patients with typical ICH. Further simulations using an imaginary simple oval hematoma revealed the formation of a large notch when the aspiration catheter moved vertically along the shorter axis (Supplementary Fig. 5). Notches were also observed around the horizontally inserted non-moving catheters. Under negative pressure, the interfacial tension at the boundary surface of the hematoma near the entrance of the catheter caused the formation of notches. When the catheter remained stationary, two notches were formed on the upper and lower planes of the hematoma surface, closer to the catheter entrance. In dynamic simulations that account for catheter movement, the starting point of catheter motion was close to the hematoma surface when the catheter was vertically pulled out. As the catheter moved, a small notch initially formed on the hematoma surface plane parallel to and near the catheter entrance, which subsequently grew in size (Supplementary Movie 3). The developed notch disrupts the aspirated flow of hematoma, thereby increasing the likelihood of flow disruption and residue formation. In the simulation, the advantage of the SOK route in the typical group can be attributed to the (i) reduced movement at the hematoma surface parallel (orange) and close to the catheter entrance at the beginning, (ii) reduced distance between the hematoma surface planes (green) to the catheter entrance during catheter movement, and (iii) increased total stroke along the long axis of the hematoma. The distances from the catheter entrance to the parallel (orange) or normal planes (green) are shown [Supplementary Fig. 5C,D].

Simulations under two-phase flow conditions can predict both the volume and morphology of the residual hematoma after aspiration. To avoid sharp edges between meshes, phase 1 (hematoma) must be attached inside phase 2 (brain) under such large deformation conditions (Supplementary Fig. 6). During a rapid change in hematoma volume, meshes in the single-phase simulation move only under the constraints of neighboring meshes and form unbound acute angles. There was an error regardless of the mesh size. The catheter speed and aspiration pressure can be adjusted by clinicians. In our simulation, the negative pressure generated by a 10-cc syringe was calculated assuming that the piston was fully retracted. In actual clinical treatment, the catheter moves back and forth multiple times; during our simulation, the catheter moved in a single stroke. Using individual CT data and personal computers, the simulation was completed within a few hours. The developed method can calculate the residual hematoma after aspiration along a specific path, along with the size and morphology of each hematoma. This approach enables physicians to estimate the efficacy of moving a catheter along a specific aspiration path during the treatment of typical ICH.

Conclusions

The current study demonstrated that stereotactic hematoma aspiration via the SOK route is safe and effective in patients with basal ganglia ICH. Notably, patients with the typical shape of ICH demonstrated a significantly improved aspiration rate and more favorable outcomes. To overcome the limitations of a retrospective analysis, we compared these two routes using a computational simulation of hematoma aspiration. The findings from our simulation revealed a superior aspiration rate via the SOK route in patients with typical basal ganglia ICH.