1. Introduction
Brain disorders affect 6.75% of the American adult population [
1] and constitute a significant cause of morbidity worldwide, with their incidence steadily increasing [
2]. Prediction of treatment outcomes for brain disorders can be difficult, even in patients with similar clinical conditions. Clinical overlaps and similarities in medical imaging can also lead to misdiagnosis between various brain disorders such as normal pressure hydrocephalus (NPH) and Parkinson’s or Alzheimer’s disease [
3,
4,
5]. The complexities in treatment outcomes may be related to the lack of clear and comprehensive knowledge concerning the mechanisms behind brain disorders. Brain tissue is a complex biomaterial with a highly heterogeneous and anisotropic microstructure, variable biomechanical properties, and dynamic behavior that can change in response to different loadings due to various brain disorders [
6]. The brain’s ability to adapt to volume change is important to clarify the complexities of some brain disorders. Intracranial compliance (ICC) refers to the volume-buffering capacity and the ability of the intracranial system to adjust its volume in order to maintain stable intracranial pressure (ICP). Maintaining stable ICP is crucial for cerebral blood flow system, preventing tissue damage, and ensuring optimal neural activity [
7]. ICC is a critical parameter to understanding the mechanisms underlying some brain disorders and the brain’s responses to various pathological processes [
8]. It can also be useful in predicting clinical outcomes in some brain disorders such as hydrocephalus, cerebral atrophy, intracranial hypertension, traumatic brain injury, and intracranial hematoma. Therefore, gaining a conceptual understanding of how to change ICC in various types of brain disorders and how it can be used in different brain pathologies for its clinical applications is of great significance in the diagnosis and management of brain disorders.
The basis for ICC measurement was first suggested 70 years ago in an animal study [
9]. Subsequent research has consistently emphasized the clinical significance, efficacy, and sensitivity of ICC in diagnosing and evaluating treatment outcomes across a spectrum of brain disorders. However, the practical implementation of ICC as an evaluative tool for neurosurgeons presents intricate challenges. In contrast, ICP stands as a well-established practical metric for studying a wide array of brain disorders. The threshold levels of ICP for both healthy individuals and patients have also been thoroughly defined [
10]. Nevertheless, ICP cannot function as a substitute for ICC due to its inability to account for volume changes and its lower sensitivity in evaluating specific disorders, such as NPH. The present review is dedicated to addressing these concerns and discussing the challenges and difficulties associated with the practical application of ICC for neurosurgeons in evaluating brain disorders.
2. Measurement, calculation, and estimation of ICC
The primary techniques for direct measurement of ICC involve monitoring changes in cerebrospinal fluid (CSF) volume within the craniospinal space. These techniques are generally based on the Marmarou method, and they are invasive, carrying a higher risk of complications. Furthermore, in specific clinical scenarios, even a minor alteration in CSF volume can result in a detrimental change in ICP, leading to unpredictable clinical consequences. Okon et al also revealed that these techniques, which generally involve the injection of excessive fluid, might not be suitable in cases of elevated ICP [
11]. Similarly, Smielewski et al demonstrated that these techniques, including constant infusion, constant pressure infusion, bolus manipulation, and lumbar ventricular perfusion, could lead to potential errors in pressure readings due to possible vasomotor responses [
12]. Hence, there was a need to explore alternatives for assessing ICC without adding or removing fluid from the craniospinal space which might involve manipulating the natural CSF circulation system and non-invasive methods emerged as superior options. Pioneering efforts by Alperin et al marked the first attempt to calculate ICC non-invasively through computer simulations [
13,
14]. We also developed and optimized the non-invasive ICC calculation employing computer simulations and mathematical analyses [
15,
16,
17,
18,
19]. Additionally, we designed and fabricated an in vitro model for experimental ICC measurements [
20]. Several other studies have endeavored to gauge ICC non-invasively, although they were unable to directly measure or calculate ICC values. Instead, they attempted to estimate ICC effects indirectly based on measurable parameters such as optic nerve sheath diameter [
21], transcranial Doppler ultrasound [
22], ICP waveform morphology [
23], and bioelectrical models [
24]. These methods for ICC measurement, calculation, and estimation have been introduced and utilized in some studies to assess ICC. For practical applications of these methods, the choice of method for assessing ICC can depend on the individual case and the preferences of the physician. For instance, estimation methods, although they cannot assess exact ICC values and can only estimate the general trend of ICC changes, can be suggested more viable in neurological emergency situations compared to relatively slower computer simulations or invasive ICC measurements.
3. Viscous component of the brain
The brain is composed of different substructures, including white and gray matter, blood vessels, neurons, and fluids such as extracellular fluid. Previous studies modeled the brain as a poroelastic structure based on brain consolidation theory [
25]. The brain poroelastic model assumed the brain to be a porous medium consisting of solid and fluid phases. This model described the time-dependent behavior of both solid and fluid phases of the brain and the interaction between them. On the other hand, brain substructures have different viscous time-dependent properties that contribute to the brain’s overall damping characteristics. Previous studies also confirmed that a viscous time-dependent component must be considered in the definition of the brain model [
26]. Hence, in addition to the importance of the scale of strain (short or large) in the brain model, the viscous time-dependent component is also necessary to achieve agreement between mathematical biphasic characteristics and experimental material properties of the brain [
27]. Hrapko et al used a large strain viscoelastic framework to define one of the most prominent time-dependent stress-strain models for the brain [
28]. Various studies have also tried to consider a viscous time-dependent component in the brain model under a viscoelastic, hyper-viscoelastic, or poro-hyperviscoelastic model (
Table 1). Gholampour et al and Cheng et al used a poro-viscoelastic model to define the most appropriate constitutive model for the brain [
15,
18,
29]. On the other hand, Elkin et al showed that the best conformity with experimental data is obtained when the viscoelastic component of the brain is fitted to the shear modulus (
) using the Prony series [
30]. Hence, we applied this method to consider the parameter of time (t) in our poro-viscoelastic brain model using equation 1 [
15,
18].
where
,
, and
k in equation 1 are input shear modulus, relaxation modulus, and relaxation time, respectively. It can be deduced that, regarding the proven impact of the viscous component in the mathematical and computational models of the brain, the role of time and the time-dependency of the brain are irrefutable in studying the brain and, consequently, brain disorders.
4. The role of time in ICC definition
The measurement and calculation of ICC in all previous methods were based on the classic equation of ICC, conventionally defined using the change in volume divided by the change in ICP (Eq. 2). The industrial materials, except for specific types of composites, are non-viscous and their behavior is independent of time. Hence, the parameters used in the ICC equation (changes in volume and pressure) have adequate potential to define compliance in these materials. Nevertheless, the brain models mentioned above demonstrated that the physiological and biomechanical functions of the brain are time-dependent (as discussed in the preceding section). Our previous study also showed the importance of the creep behavior (time-dependent deformation) of brain tissue in the treatment process of some brain disorders such as hydrocephalus [
16]. Furthermore, the onset and progression of certain brain disorders occur over time. For instance, NPH, intracranial hypertension and hydrocephalus often manifest gradually, with many patients experiencing elevated ICP over extended periods [
15,
31]. Previous studies have also proved time-dependency and long-term alterations in CSF dynamics and brain morphometric parameters in patients with CSF disorders [
15,
19,
32,
33]. Therefore, in addition to the defined parameters in the ICC equation, time may also play a significant role in the formulation and, consequently, in the measurement/calculation of ICC. Certain previous studies formulated compliance based on pressures and volumes of cerebral blood and CSF [
34,
35,
36]. Some studies defined ICC based on the elastic modulus of brain tissue, ICP, and intracranial volume [
37,
38,
39]. Additionally, others formulated ICC based on the dynamics of injected CSF, ICP, and resistance to CSF outflow in invasive ICC measurement methods [
40]. However, these studies have not directly considered the parameter of time in their ICC equations. We considered that the ICC equation can remain in its general form (Eq. 2) and consider the effect of time in the concept of “change”. “Change” of volume or ICP in the ICC equation (Eq. 2) means V
2-V
1 or ICP
2-ICP
1 [
31]. It can be deduced that we can consider the time elapsed (TE) between the measurements of V
2 and V
1, or ICP
2 and ICP
1 as representative of the parameter of time in the ICC equation.
Table 1.
Different time-dependent brain models.
Table 1.
Different time-dependent brain models.
Brain model |
Authors, year |
Type of brain disorder |
Solving method |
Brain regions |
Source |
Poroelastic |
Yuan Et al, 2022 |
Healthy subjects under drug infusion |
Mathematical analysis based on arbitrary Lagrange-Eulerian equations |
White matter |
[41] |
Lambride Et al, 2020 |
Brain injury |
Finite element method |
Single region |
[42] |
Guo Et al, 2018 |
Alzheimer’s disease |
Finite element method |
White matter |
[43] |
Gholampour et al, 2014 and 2015 |
Non-communicating hydrocephalus |
Fluid-structure interaction |
Single region |
[44,45] |
Viscoelastic |
Li et al, 2021 |
Healthy subject |
Finite element method |
Grey and white matter |
[46] |
Siegkas et al, 2019 |
Brain injury |
Finite element method |
Single region |
[47] |
Gholampour et al, 2017 |
Hydrocephalus |
Fluid-structure interaction |
Single region |
[19,33] |
Harpko et al, 2006 |
Healthy subject |
Mathematical analysis |
White matter |
[28] |
Hyper-viscoelastic |
Menghani et al, 2023 |
Head impact |
Finite element method |
Basal ganglia, cerebral hemispheres, and corpus callosum |
[48] |
Wang et al, 2018 |
Brain injury |
Finite element method |
Grey matter, white matter, and pia mater |
[49] |
Wilkie et al, 2012 |
Hydrocephalus |
Mathematical analysis using fractional Zener model |
Single region |
[50] |
Dutta-Roy, 2011 |
Normal pressure hydrocephalus |
Finite element method |
Single region |
[51] |
Poro-viscoelastic |
Gholampour et al, 2022 and 2023 |
Communicating hydrocephalus |
Fluid-structure interaction |
Single region |
[16,17,31] |
Pavan Et al, 2022 |
Brain injury |
Finite element method |
One region |
[52] |
Gholampour, 2018 |
Non-communicating hydrocephalus |
Fluid-structure interaction |
Single region |
[15] |
Cheng et al, 2010 |
Non-communicating hydrocephalus |
Finite element method |
White matter |
[29] |
Poro-hyperviscoelastic |
Hosseini-Farid et al, 2020 |
Healthy subject |
Finite element method |
Grey and white matter |
[53] |
Forte et al, 2017 |
Healthy subject |
Finite element method |
Grey and white matter |
[54] |
5. Approaches to TE in ICC assessment
Differences in ICC values across different assessment methods, as well as different TE values, are evident in
Table 2. These significant differences, characterized by inconsistency ICC results, may provide insight into why, despite the clinical importance of ICC, its monitoring, unlike ICP monitoring, has not yet become a standard component of clinical care [
55].
Table 2 also indicated that the values of TE in previous non-invasive ICC measurement methods (in vitro, lumped model, and computer simulation) were less than one minute [
15,
20]. The corresponding values for TE in previous invasive ICC measurement methods (bolus injection and lumbar or ventricular infusion) had different values (
Table 2). These differences may be attributed to variations in the clinical conditions of the patients and the different time taken to reach plateau ICP during the test. Hence, a question will be raised whether these TEs used in previous non-invasive and invasive ICC measurements are sufficiently long to demonstrate the complete effect of the brain’s viscous time-dependent component. This component, reflecting the brain’s load-history-dependent behavior, pertains to its ability to dissipate external loading and resist deformation. Consequently, considering the influence of this component is essential for accurate assessment of ICC. Hence, addressing this question is essential to ensuring the effective application of ICC in practical scenarios. Recently, we demonstrated that before two months, the complete effect of the viscous component of the brain had not yet appeared in the hydrocephalic brain [
16,
31]. This means that the appropriate TE for measuring ICC in hydrocephalus patients should not be short and should be at least two months. Evaluating the brains of healthy rats subjected to loading due to electrode implantation also revealed that the effects of loading did not completely dissipate within two months and brain material properties, such as shear and elastic modules, did not stabilize during this period, somehow mirroring ICC behavior [
56]. Similarly, Boulet et al obtained comparable results, confirming significant changes in brain stiffness in rat brains over a month-long period [
57]. Therefore, building on previous studies, we can categorize the evaluation of ICC into two approaches: ICC in a short TE and ICC in a large TE [
31]. Comparing ICC behavior using these two different approaches can be helpful in clarifying certain complexities related to application of ICC for some brain disorders. For example, Eide and Brean showed when the ICC decreases in NPH patients in a short TE, the brain material behaves similarly to a linear elastic [
58]. While, another study showed that the changes in elastance, stiffness, and creep of the brain in a large TE were non-linear and non-monotonic in patients with communicating hydrocephalus [
16]. Therefore, the correct choice of TE (short or large) directly affects the degree to which the impacts of the viscous component of the brain are realized and also affects the arrival of brain material properties to a stable condition, allowing for the correct assessment of ICC.
6. Discussion
The complex and dynamic behavior of the brain can change in response to various loadings associated with different brain disorders, and ICC plays a crucial role in understanding the mechanisms and complexities underlying these changes. As a result, the challenges pertaining to the practical application of ICC are of significant importance in the context of brain disorders, and this review is specifically designed to address these concerns. We emphasized the significance of TE as a representation of the time parameter in the ICC equation, drawing from previous research findings. Additionally, our research highlights that in certain brain disorders, the changes in brain material properties do not stabilize within a short TE, emphasizing the dynamic nature of these conditions across various time scales, including a large TE. Furthermore, our exploration reveals that the selection of TE, be it short or large, directly influences the manifestation of the brain’s viscous component. Consequently, this review highlights that the selection of TE (short or large) directly impacts the accuracy of ICC assessment.
In addition to assessing ICC behavior under two TE approaches, it is crucial to analyze the impact of these approaches (short or large TE) on the changes in the volume-ICP curve (VI curve). Because ICC represents the slope of this curve (Eq. 2), and the curve itself plays a prominent role in understanding the compensatory response of the brain in different brain disorders. Previous studies strongly believed the trend of the VI curve to be monotonic over short TE in patients with different types of brain disorders (
Figure 1a) [
13,
67,
68]. However, Okon et al expressed doubts about the monotonic trend of VI curves in the short TE [
69]. They showed that the trend of this curve in patients with idiopathic intracranial hypertension was not necessarily monotonic. Our evaluations showed that the trend of this curve is non-monotonic in hydrocephalus patients in a large TE (fifteen months after treatment) (
Figure 1b and
Figure 2a) [
16,
31]. We also showed that in a large TE, ICC does not necessarily increase with decreasing ICP level (
Figure 2b), and the ICC oscillatory increases during the treatment process (
Figure 2b) [
16,
31]. It can also be understood from a study by Tisell et al that there is a non-monotonic trend in the VI curve over a large TE (three months after treatment) in hydrocephalus patients [
70]. Hence, the compensatory response of the brain can have non-monotonic and variable changes in a large TE in certain brain disorders (
Figure 2). Therefore, it can be deduced that the trend of the compensatory response of the brain, which plays a prominent role in ICC assessments, may be different over short or large TE in certain brain disorders.
The correct choice of TE (short or large) depends on the strain rate and loading time in that particular brain disorder under assessment for ICC. Strain rate is a measure of how quickly the brain is being deformed under loading caused by a brain disorder. The loading time refers to the duration of the loading that is applied to the brain in a specific brain disorder and it is a measure of how long the brain is being subjected to loading. The strain rate and loading time vary across different types of brain disorders. For instance, these parameters significantly differ in traumatic brain injury and some time-dependent brain disorders such as hydrocephalus. Therefore, it can be inferred that the practical utilization of ICC depends on the strain rate and loading duration specific to each individual brain disorder, taking into account the two approaches used in ICC assessment.
Previous studies showed descents in the VI curve at certain times in a large TE for some disorders like hydrocephalus (
Figure 1b) because of the non-monotonic trend of this curve [
16,
31]. Similar descents were also demonstrable in a study by Tisell et al [
70]. This means in the context of ICC assessment using a large TE approach, ICC values and consequently the compensatory response of brain tissue, may partially recover at specific times. This is worth mentioning when the brain is in an uncompensated status in a short TE, it is traditionally defined as a high-risk condition [
55]. However, the temporary uncompensatory responses of the brain caused by the non-monotonic trend of the VI curve in a large TE can not necessarily be deduced as a high-risk condition. In addition, two separate compensatory reserve zones were shown in the VI curve in a short TE (
Figure 1a: blue and green colors). The ICC values and pulse amplitude of ICP change in these two zones. However, the VI curve in a large TE was not divisible into specific meaningful compensatory reserve zones in hydrocephalus patients (
Figure 1b). Overall, the general trend (mono-exponential) and compensatory reserve zones (two zones) of the VI curve in a short TE were approximately the same in different types of brain disorders such as hydrocephalus [
13,
71], brain edema [
67], and brain injury [
68], and they were different from a large TE approach.
In order to measure or calculate ICC, two statuses, labeled 1 and 2, are required to calculate the ‘change’ in volume and ICP (Eq. 2). In emergency conditions of certain brain disorders, such as hydrocephalus, which involve a large strain rate and loading duration, the ICC in a large TE approach may not be applicable due to the time needed to establish these two statuses for the ICC equation. In such scenarios, ICC estimation methods may prove more useful. However, in cases like NPH, ICC can be exclusively used for diagnosis and evaluating treatment outcomes. This is because the conditions of these patients are not always emergencies; NPH is generally a time-dependent disease, and other available indicators such as ICP, gait analysis, or cognitive examinations may not always be effective to evaluate these diseases [
72,
73]. It’s worth noting that there is no established numerical threshold for ICC values for the practical use of ICC as a clinical indicator for these conditions, even for diseases like NPH or in non-emergency situations. Expanding optimized computer simulation methods [
18] for non-invasive ICC calculations in future studies could be beneficial in addressing this concern by defining a threshold level between ICC values in healthy subjects and those in patients.
The recovery behavior of the brain during the treatment process (unloading condition) is of great importance in the management of brain disorders [
16]. The previous studies have demonstrated that the biomechanical response of the brain in loading (i.e., afflicted with a brain disorder) and unloading (i.e., during some treatment processes) conditions is different [
74]. Elastic hysteresis behavior of the brain is one of the primary reasons for this difference specifically in disorders caused by
cyclic loading, i.e. in hydrocephalus and intracranial hypertension that are caused by
pulsatile elevated ICP. The theory of elastic hysteresis of the brain states that the brain stores energy during loading, however, not all of this energy is released during unloading, as some of it is dissipated due to internal friction. This leads to a residual (permanent) deformation in brain tissue after the treatment process. The differences in residual deformation of the hydrocephalic brain were reported in short and large TEs [
16,
27]. This theory which is provided by Lesniak et al confirmed the direct impact of elastic hysteresis behavior of the brain on ICC changes [
75]. Comparing ICC values in loading and unloading conditions over a large TE in future studies can be useful for developing the theory of elastic hysteresis. This comparison can further clarify how ICC changes in loading and unloading conditions in a large TE, potentially enhancing the practical utility of ICC in time-dependent brain disorders. In addition, previous studies showed a decrease in glymphatic drainage and ICC value in a short TE for patients with NPH, traumatic brain injury, and Alzheimer’s diseases [
76,
77,
78]. Evaluation of the changes in ICC and glymphatic drainage in a large TE in future studies may be useful to support the hypothesis of glymphatic dysfunction as the underlying pathophysiology of time-dependent brain disorders such as NPH and Alzheimer’s diseases