Giulia Pederzani,Andrii Grytsan,Alfons G. Hoekstra,Anne M. Robertson,Paul N. Watton
Abstract:Cerebral vasospasm, a prolonged constriction of cerebral arteries, is the first cause of morbidity and mortality for patients who survive hospitalisation after aneurysmal subarachnoid haemorrhage. The recent finding that stent-retrievers can successfully treat the disease has challenged the viewpoint that damage to the extracellular matrix is necessary. We apply a 3D finite element rate-based constrained mixture model (rb-CMM) to simulate vasospasm, remodelling and treatment with stents. The artery is modelled as a thick-walled fibre-reinforced constrained mixture subject to physiological pressure and axial stretch. The model accounts for distributions of collagen fibre homeostatic stretches, VSMC active response, remodelling and damage. After simulating vasospasm and subsequent remodelling of the artery to a new homeostatic state, we simulate treatment with commonly available stent-retrievers. We perform a parameter study to examine how arterial diameter and thickness affect the success of stent treatment. The model predictions on the pressure required to mechanically resolve the constriction are consistent with stent-retrievers. In agreement with clinical observations, our model predicts that stent-retrievers tend to be effective in arteries of up to 3mm diameter, but fail in larger ones. Variations in arterial wall thickness significantly affect stent pressure requirements. We have developed a novel rb-CMM that accounts for VSMC active response, remodelling and damage. Consistently with clinical observations, simulations predict that stent-retrievers can mechanically resolve vasospasm. Moreover, accounting for a patient's arterial properties is important for predicting likelihood of stent success. This in silico tool has the potential to support clinical decision-making and guide the development and evaluation of dedicated stents for personalised treatment of vasospasm.
What problem does this paper attempt to address?
The problems that this paper attempts to solve are: **How to understand and simulate the occurrence, remodeling, and treatment mechanisms of cerebral vasospasm through computational modeling, especially the effectiveness of mechanical treatment using stent - retrievers**.
Specifically, the paper aims to:
1. **Explain the pathophysiological mechanisms of cerebral vasospasm**: In particular, the role of the active response and remodeling of vascular smooth muscle cells (VSMCs) in cerebral vasospasm after subarachnoid hemorrhage (SAH).
2. **Evaluate the treatment effect of stent - retrievers**: Explore how stent - retrievers can successfully treat cerebral vasospasm without damaging the blood vessel wall.
3. **Predict the influence of different arterial diameters and thicknesses on the treatment effect of stents**: Study how the geometric characteristics of arteries affect the success rate of stent treatment.
4. **Develop a new computational model**: Based on the rate - based constrained mixture model (rb - CMM), this model can simulate the multi - layer structure of the blood vessel wall and its changes in healthy, spastic, and post - mechanical treatment states.
### Background of the Paper
Cerebral vasospasm is a serious phenomenon of continuous constriction of cerebral arteries, usually caused by subarachnoid hemorrhage. It can lead to a decrease in cerebral blood flow, which in turn causes brain tissue hypoxia, and is the main cause of death and disability in SAH patients after hospitalization. Traditional treatment methods include drug treatment and mechanical treatment (such as balloon angioplasty). However, in recent years, it has been found that low - pressure stent - retrievers can also effectively treat cerebral vasospasm, which challenges the previous view that the extracellular matrix (ECM) must be destroyed to relieve the constriction.
### Research Methods
The researchers adopted a three - dimensional finite - element rate - based constrained mixture model (3D finite element rate - based constrained mixture model, rb - CMM) to simulate the multi - layer structure of the arterial wall, including elastic fibers, vascular smooth muscle cells, and collagen fibers. The model takes the following factors into account:
- **Physiological pressure and axial stretch**: To simulate in - vivo conditions.
- **Active response and remodeling of VSMCs**: To simulate the changes of vascular smooth muscle cells at different stages.
- **Distribution and damage of collagen fibers**: Consider the initial waviness of collagen fibers and their changes at different stages.
- **Pressure application in mechanical treatment**: To simulate the pressure applied by stent - retrievers.
### Main Results
The results predicted by the model are consistent with clinical observations, indicating that:
- Stent - retrievers can successfully relieve cerebral vasospasm without destroying the ECM.
- For arteries with a diameter less than 3 millimeters, stent - retrievers are usually effective; while for larger arteries, they usually fail.
- Changes in arterial wall thickness can affect the required stent pressure, and a higher pressure may be required for successful treatment.
### Conclusions
This study has developed a new rb - CMM model that can simulate the active response, remodeling, and damage of VSMCs. The model prediction results are consistent with clinical observations, indicating that stent - retrievers can relieve cerebral vasospasm mechanically. In addition, considering the arterial characteristics of patients is crucial for predicting the success rate of stent treatment. This computational tool is expected to support clinical decision - making and guide the design and evaluation of personalized stents.
### Formula Summary
The formulas involved in the paper are as follows:
1. **Decomposition of the deformation gradient**:
\[
F = F_e\cdot F_g
\]
where \( F_e \) represents elastic deformation and \( F_g \) represents the growth component.
2. **Strain - energy - density function (SEDF)**:
\[
\Psi=\Psi(F)-p(\det F_e - 1)
\]
where \( p \) is the Lagrange multiplier, representing hydrostatic pressure.
3. **Second Piola - Kirchhoff stress**:
\[
S=-pC^{-1}+2\left[\frac{\partial\Psi}{\partial I_1}I+\frac{\partial\Psi}{\partial