The dynamics of vein graft remodeling induced by hemodynamic forces: a mathematical model

Importance of Flow Waveform and Flow Volume as Prognostic Indicators for the Patency of Infra-Inguinal Autologous Vein Bypass

OBJECTIVE

To investigate the association of the intra-operative flow waveform and the flow volume with graft prognosis of the infra-inguinal vein bypass.

METHODS

This was a retrospective study of intra-operative flowmetry performed for infra-inguinal autologous vein bypass between 2011 and 2020. Flow waveforms were classified as type 0 - IV according to the Kyushu University (KU) classification. The patients (n = 340) were divided into three groups based on the flow waveform predicting the graft patency: type 0/I (long patency), type II (no early occlusion but late occlusion possible), and type III/IV (early occlusion). The graft occlusion rates of popliteal artery bypass (PAB) and infrapopliteal artery bypass (IPAB) within 30 days of surgery were compared between type 0/I + II and type III/IV groups, while the midterm graft patency rates were compared between type 0/I and type II groups. Additionally, a multivariate analysis was performed to identify independent risk factors for early and late graft occlusion.

RESULTS

The early graft occlusion rates of type 0/I + II and type III/IV groups were 3.9% and 0%, respectively, (p = 1.0) for PAB, and 5.3% and 46.2%, respectively, (p < .001) for IPAB. The two year primary patency rates of type 0/I and type II groups were 91% and 75%, respectively, (p = .030) for PAB, and 58% and 63%, respectively, (p = .72) for IPAB. Independent risk factors for early occlusion were none in PAB and flow waveform (type IV) in IPAB. Independent risk factors for patency loss in PAB were flow waveform (type II), end stage renal disease, and dual antiplatelet use, and those in IPAB were older age, women, lower flow volume, and iterative bypass.

CONCLUSION

Intra-operative flowmetry is useful for predicting the graft prognosis in infra-inguinal vein bypass and this is dependent on the distal target artery.

Multiscale Computational Modeling of Vascular Adaptation: A Systems Biology Approach Using Agent-Based Models

The widespread incidence of cardiovascular diseases and associated mortality and morbidity, along with the advent of powerful computational resources, have fostered an extensive research in computational modeling of vascular pathophysiology field and promoted in-silico models as a support for biomedical research. Given the multiscale nature of biological systems, the integration of phenomena at different spatial and temporal scales has emerged to be essential in capturing mechanobiological mechanisms underlying vascular adaptation processes. In this regard, agent-based models have demonstrated to successfully embed the systems biology principles and capture the emergent behavior of cellular systems under different pathophysiological conditions. Furthermore, through their modular structure, agent-based models are suitable to be integrated with continuum-based models within a multiscale framework that can link the molecular pathways to the cell and tissue levels. This can allow improving existing therapies and/or developing new therapeutic strategies. The present review examines the multiscale computational frameworks of vascular adaptation with an emphasis on the integration of agent-based approaches with continuum models to describe vascular pathophysiology in a systems biology perspective. The state-of-the-art highlights the current gaps and limitations in the field, thus shedding light on new areas to be explored that may become the future research focus. The inclusion of molecular intracellular pathways (e.g., genomics or proteomics) within the multiscale agent-based modeling frameworks will certainly provide a great contribution to the promising personalized medicine. Efforts will be also needed to address the challenges encountered for the verification, uncertainty quantification, calibration and validation of these multiscale frameworks.

Plasticity of the Maternal Vasculature During Pregnancy

Maternal cardiovascular changes during pregnancy include an expansion of plasma volume, increased cardiac output, decreased peripheral resistance, and increased uteroplacental blood flow. These adaptations facilitate the progressive increase in uteroplacental perfusion that is required for normal fetal growth and development, prevent the development of hypertension, and provide a reserve of blood in anticipation of the significant blood loss associated with parturition. Each woman's genotype and phenotype determine her ability to adapt in response to molecular signals that emanate from the fetoplacental unit. Here, we provide an overview of the major hemodynamic and cardiac changes and then consider regional changes in the splanchnic, renal, cerebral, and uterine circulations in terms of endothelial and vascular smooth muscle cell plasticity. Although consideration of gestational disease is beyond the scope of this review, aberrant signaling and/or maternal responsiveness contribute to the etiology of several common gestational diseases such as preeclampsia, intrauterine growth restriction, and gestational diabetes.

Venoarterial communication mediates arterial wall shear stress-induced maternal uterine vascular remodeling during pregnancy

Although expansive remodeling of the maternal uterine circulation during pregnancy is essential for maintaining uteroplacental perfusion and normal fetal growth, the underlying physiological mechanisms are not well understood. Using a rat model, surgical approaches were used to alter uterine hemodynamics and wall shear stress (WSS) to evaluate the effects of WSS and venoarterial communication (e.g., transfer of placentally derived growth signals from postplacental veins to preplacental arteries) on gestational uterine vascular remodeling. Changes in WSS secondary to ligation of the cervical but not the ovarian end of the main uterine artery and vein provoked significant expansive remodeling at the opposite end of both vessels, but only in pregnant animals. The ≈50% increase in lumen diameter (relative to the contralateral horn) was associated with an upregulation of total endothelial nitric oxide (NO) synthase expression and was abolished by in vivo NO synthase inhibition with N-nitro-l-arginine methyl ester. Complete removal of a venous segment adjacent to the uterine artery to eliminate local venous influences significantly attenuated the WSS-induced remodeling by about one-half ( P < 0.05). These findings indicate that, during pregnancy, 1) increased WSS stimulates uterine artery growth via NO signaling and 2) the presence of an adjacent vein is required for arterial remodeling to fully occur. NEW & NOTEWORTHY This study provides the first in vivo evidence for the importance of venous influences on arterial growth within the uteroplacental circulation.

Patient-Specific, Multi-Scale Modeling of Neointimal Hyperplasia in Vein Grafts

Neointimal hyperplasia is amongst the major causes of failure of bypass grafts. The disease progression varies from patient to patient due to a range of different factors. In this paper, a mathematical model will be used to understand neointimal hyperplasia in individual patients, combining information from biological experiments and patient-specific data to analyze some aspects of the disease, particularly with regard to mechanical stimuli due to shear stresses on the vessel wall. By combining a biochemical model of cell growth and a patient-specific computational fluid dynamics analysis of blood flow in the lumen, remodeling of the blood vessel is studied by means of a novel computational framework. The framework was used to analyze two vein graft bypasses from one patient: a femoro-popliteal and a femoro-distal bypass. The remodeling of the vessel wall and analysis of the flow for each case was then compared to clinical data and discussed as a potential tool for a better understanding of the disease. Simulation results from this first computational approach showed an overall agreement on the locations of hyperplasia in these patients and demonstrated the potential of using new integrative modeling tools to understand disease progression.

Hemodynamic Influence on Smooth Muscle Cell Kinetics and Phenotype During Early Vein Graft Adaptation

Pathologic vascular adaptation following local injury is the primary driver for accelerated intimal hyperplasia and an occlusive phenotype. Smooth muscle cell (SMC) proliferation within the wall, and migration into the developing intima, is a major component of this remodeling response. The primary objective in the current study was to investigate the effect of the local biomechanical forces on early vein graft adaptation, specifically focusing on the spatial and temporal response of SMC proliferation and conversion from a contractile to synthetic architecture. Taking advantage of the differential adaptation that occurs during exposure to divergent flow environments, vein grafts were implanted in rabbits to create two distinct flow environments and harvested at times ranging from 2 h to 28 days. Using an algorithm for the virtual reconstruction of unfixed, histologic specimens, immunohistochemical tracking of DNA synthesis, and high-throughput transcriptional analysis, the spatial and temporal changes in graft morphology, cell proliferation, and SMC phenotype were catalogued. Notable findings include a burst of cell proliferation at 7 days post-implantation, which was significantly augmented by exposure to a reduced flow environment. Compared to the adjacent media, proliferation rates were 3-fold greater in the intima, and a specific spatial distribution of these proliferating cells was identified, with a major peak in the sub-endothelial region and a second peak centering on the internal elastic lamina. Genomic markers of a contractile SMC phenotype were reduced as early as 2 h post-implantation and reached a nadir at 7 days. Network analysis of upstream regulatory pathways identified GATA6 and KLF5 as important transcription factors that regulate this shift in SMC phenotype.

SAPHENOUS VEIN GRAFT PATENCY AFTER GEOMETRY REMODELING

Tortuous saphenous vein graft (SVG) hemodynamics was investigated using computational fluid dynamics (CFD) techniques. Computed tomography (CT) technology is used for noninvasive bypass graft assessment seven days after surgery. CT investigation shows two regions with severe shape remodeling, one is an angle type contortion and the other one is a sharp curvature with tortuous area reduction. The numerical analysis carefully examines the effect of an SVG geometry remodeling through flow separation, particle deposition, and wall shear stress (WSS). During the cardiac cycle, overall pressure drop increases from 2.6[Formula: see text]mmHg to 4.4[Formula: see text]mmHg. In the accelerating part of the systolic phase, particles released in the inlet section move downstream toward the first narrowed part (elbow type contortion) with a helical motion. WSS range along the cardiac cycle varies from 2[Formula: see text]Pa to 42[Formula: see text]Pa, enough to damage the endothelial cells. Vessel torsion induced helical flow can reduce the flow disturbance and separation. Additionally, in the distal end of the graft, the high particle concentrations can promote the inflammatory processes in the vessels.

A mechanical argument for the differential performance of coronary artery grafts

Coronary artery bypass grafting (CABG) acutely disturbs the homeostatic state of the transplanted vessel making retention of graft patency dependent on chronic remodeling processes. The time course and extent to which remodeling restores vessel homeostasis will depend, in part, on the nature and magnitude of the mechanical disturbances induced upon transplantation. In this investigation, biaxial mechanical testing and histology were performed on the porcine left anterior descending artery (LAD) and analogs of common autografts, including the internal thoracic artery (ITA), radial artery (RA), great saphenous vein (GSV) and lateral saphenous vein (LSV). Experimental data were used to quantify the parameters of a structure-based constitutive model enabling prediction of the acute vessel mechanical response pre-transplantation and under coronary loading conditions. A novel metric Ξ was developed to quantify mechanical differences between each graft vessel in situ and the LAD in situ, while a second metric Ω compares the graft vessels in situ to their state under coronary loading. The relative values of these metrics among candidate autograft sources are consistent with vessel-specific variations in CABG clinical success rates with the ITA as the superior and GSV the inferior graft choices based on mechanical performance. This approach can be used to evaluate other candidate tissues for grafting or to aid in the development of synthetic and tissue engineered alternatives.

Computational simulation of the adaptive capacity of vein grafts in response to increased pressure

Vein maladaptation, leading to poor long-term patency, is a serious clinical problem in patients receiving coronary artery bypass grafts (CABGs) or undergoing related clinical procedures that subject veins to elevated blood flow and pressure. We propose a computational model of venous adaptation to altered pressure based on a constrained mixture theory of growth and remodeling (G&R). We identify constitutive parameters that optimally match biaxial data from a mouse vena cava, then numerically subject the vein to altered pressure conditions and quantify the extent of adaptation for a biologically reasonable set of bounds for G&R parameters. We identify conditions under which a vein graft can adapt optimally and explore physiological constraints that lead to maladaptation. Finally, we test the hypothesis that a gradual, rather than a step, change in pressure will reduce maladaptation. Optimization is used to accelerate parameter identification and numerically evaluate hypotheses of vein remodeling.

Numerical analysis of non-Newtonian blood flow and wall shear stress in realistic single, double and triple aorto-coronary bypasses

Considering the fact that hemodynamics plays an important role in the patency and overall performance of implanted bypass grafts, this work presents a numerical investigation of pulsatile non-Newtonian blood flow in three different patient-specific aorto-coronary bypasses. The three bypass models are distinguished from each other by the number of distal side-to-side and end-to-side anastomoses and denoted as single, double and triple bypasses. The mathematical model in the form of time-dependent nonlinear system of incompressible Navier-Stokes equations is coupled with the Carreau-Yasuda model describing the shear-thinning property of human blood and numerically solved using the principle of the SIMPLE algorithm and cell-centred finite volume method formulated for hybrid unstructured tetrahedral grids. The numerical results computed for non-Newtonian and Newtonian blood flow in the three aorto-coronary bypasses are compared and analysed with emphasis placed on the distribution of cycle-averaged wall shear stress and oscillatory shear index. As shown in this study, the non-Newtonian blood flow in all of the considered bypass models does not significantly differ from the Newtonian one. Our observations further suggest that, especially in the case of sequential grafts, the resulting flow field and shear stimulation are strongly influenced by the diameter of the vessels involved in the bypassing.

A dynamical system that describes vein graft adaptation and failure

Adaptation of vein bypass grafts to the mechanical stresses imposed by the arterial circulation is thought to be the primary determinant for lesion development, yet an understanding of how the various forces dictate local wall remodeling is lacking. We develop a dynamical system that summarizes the complex interplay between the mechanical environment and cell/matrix kinetics, ultimately dictating changes in the vein graft architecture. Based on a systematic mapping of the parameter space, three general remodeling response patterns are observed: (1) shear stabilized intimal thickening, (2) tension induced wall thinning and lumen expansion, and (3) tension stabilized wall thickening. Notable is our observation that the integration of multiple feedback mechanisms leads to a variety of non-linear responses that would be unanticipated by an analysis of each system component independently. This dynamic analysis supports the clinical observation that the majority of vein grafts proceed along an adaptive trajectory, where grafts dilate and mildly thicken in response to the increased tension and shear, but a small portion of the grafts demonstrate a maladaptive phenotype, where progressive inward remodeling and accentuated wall thickening lead to graft failure.

Rule-based model of vein graft remodeling

When vein segments are implanted into the arterial system for use in arterial bypass grafting, adaptation to the higher pressure and flow of the arterial system is accomplished thorough wall thickening and expansion. These early remodeling events have been found to be closely coupled to the local hemodynamic forces, such as shear stress and wall tension, and are believed to be the foundation for later vein graft failure. To further our mechanistic understanding of the cellular and extracellular interactions that lead to global changes in tissue architecture, a rule-based modeling method is developed through the application of basic rules of behaviors for these molecular and cellular activities. In the current method, smooth muscle cell (SMC), extracellular matrix (ECM), and monocytes are selected as the three components that occupy the elements of a grid system that comprise the developing vein graft intima. The probabilities of the cellular behaviors are developed based on data extracted from in vivo experiments. At each time step, the various probabilities are computed and applied to the SMC and ECM elements to determine their next physical state and behavior. One- and two-dimensional models are developed to test and validate the computational approach. The importance of monocyte infiltration, and the associated effect in augmenting extracellular matrix deposition, was evaluated and found to be an important component in model development. Final model validation is performed using an independent set of experiments, where model predictions of intimal growth are evaluated against experimental data obtained from the complex geometry and shear stress patterns offered by a mid-graft focal stenosis, where simulation results show good agreements with the experimental data.