A dynamical system that describes vein graft adaptation and failure

Near-wall hemodynamic parameters of finger arteries altered by hand-transmitted vibration

BACKGROUND

Sustained exposure to high-level hand-transmitted vibrations may result in angioneurotic disorders, which partly originate from vibration-altered hemodynamics in the finger arteries when repeating these disturbances throughout working life. Hence, the aim of this study is to assess the most relevant hemodynamic descriptors in the digital arteries, determine the relationship between the latter and vibration features, and gain better understanding of the physiological mechanisms involved.

METHODS

An experimental setup, mainly comprised of an ultra-high frequency ultrasound scanner and a vibration shaker, was used to image the digital proper volar arteries of the forefinger. Raw ultrasound data were post-processed by custom-made numerical routines to supply a pulsatile fluid mechanics model for computing the hemodynamic descriptors. Twenty-four healthy volunteers participated in the measurement campaign. Classical statistical methods were then applied to the dataset and also the wavelet transform for calculating the signal power in the frequency bands matching cardiac, respiratory, myogenic and neurogenic activities.

RESULTS

The artery diameter, the wall shear stress - WSS - and the WSS temporal gradient - WSSTG - were found to be the most relevant descriptors. Vibration-induced WSS was divided by three compared to its basal value whatever the vibration frequency and it was proportional to log2 of the acceleration level. Marked increases in WSSTG when stopping vibration might also lead to adverse health effects. Vibration caused a drop in WSS power for the frequency band associated with the neurogenic activity of the sympathetic nervous system.

CONCLUSION

This study may pave the way for a new framework to prevent vibration-induced vascular risk.

Distinct Temporal Pattern of the Prediction of Lumen Remodeling of Lower Extremity Vein Bypass Grafts by Initial Local Hemodynamics

We predicted human lower extremity vein bypass graft remodeling by hemodynamics. Computed tomography and duplex ultrasound scans of 55 patients were performed at 1 week and 1, 6, and 12 months post-implantation to obtain wall shear stress (WSS) and oscillatory shear index (OSI) at 1-mm intervals via computational fluid dynamics simulations. Graft remodeling was quantified by computed tomography-measured lumen diameter changes in the early (1 week-1 month), intermediate (1-6 months), and late (6-12 months) periods. Linear mixed-effect models were constructed to examine the overall relationship between remodeling and initial hemodynamics using the average data of all cross sections within the same graft. A significant association of graft remodeling with WSS (p < 0.001) and time (p = 0.001) was found; however, the effect size decreased with time (every 2.7 dyne/cm² increase of WSS was associated with a 0.39, 0.35, 0.002 mm diameter increase in the three periods, respectively). The association of remodeling with OSI was significant only in the intermediate period (every 0.1 increase of OSI was associated with a 0.25 mm lumen diameter decrease, p = 0.004). Therefore, the association of graft lumen remodeling with local hemodynamics has a distinct temporal pattern; WSS and OSI are predictive of remodeling only in certain postoperative periods.

Assessing mechanical vibration-altered wall shear stress in digital arteries

The aim of this study is to implement and validate a method for assessing acute vibration-altered Wall Shear Stress (WSS) in the proper volar digital artery of the non-exposed left forefinger when subjecting the right hand to mechanical vibration. These changes of WSS may be involved in Vibration White Finger. Hence, an experimental device was set-up to link a vibration shaker and an ultra-high frequency ultrasound scanner. The Womersley-based WSS was computed by picking up the maximum velocity from pulse Wave Doppler measurements and extracting the artery diameter from B-mode images through an in-house image processing technique. The parameters of the former method were optimised on numerical ultrasound phantoms of cylindrical and lifelike arteries. These phantoms were computed with the FIELD II and FOCUS platforms which mimicked our true ultrasound device. The Womersley-based WSS were compared to full Fluid Structure Interaction (FSI) and rigid wall models built from resonance magnetic images of a volunteer-specific forefinger artery. Our FSI model took into account the artery's surrounding tissues. The diameter computing procedure led to a bias of 4%. The Womersley-based WSS resulted in misestimating the FSI model by roughly 10% to 20%. No difference was found between the rigid wall computational model and FSI simulations. Regarding the WSS measured on a group of 20 volunteers, the group-averaged basal value was 3 Pa, while the vibration-altered WSS was reduced to 1 Pa, possibly triggering intimal hyperplasia mechanisms and leading to the arterial stenoses encountered in patients suffering from vibration-induced Raynaud's syndrome.

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.

A computational bio-chemo-mechanical model of in vivo tissue-engineered vascular graft development

Stenosis is the primary complication of current tissue-engineered vascular grafts used in pediatric congenital cardiac surgery. Murine models provide considerable insight into the possible mechanisms underlying this situation, but they are not efficient for identifying optimal changes in scaffold design or therapeutic strategies to prevent narrowing. In contrast, computational modeling promises to enable time- and cost-efficient examinations of factors leading to narrowing. Whereas past models have been limited by their phenomenological basis, we present a new mechanistic model that integrates molecular- and cellular-driven immuno- and mechano-mediated contributions to in vivo neotissue development within implanted polymeric scaffolds. Model parameters are inferred directly from in vivo measurements for an inferior vena cava interposition graft model in the mouse that are augmented by data from the literature. By complementing Bayesian estimation with identifiability analysis and simplex optimization, we found optimal parameter values that match model outputs with experimental targets and quantify variability due to measurement uncertainty. Utility is illustrated by parametrically exploring possible graft narrowing as a function of scaffold pore size, macrophage activity, and the immunomodulatory cytokine transforming growth factor beta 1 (TGF-β1). The model captures salient temporal profiles of infiltrating immune and synthetic cells and associated secretion of cytokines, proteases, and matrix constituents throughout neovessel evolution, and parametric studies suggest that modulating scaffold immunogenicity with early immunomodulatory therapies may reduce graft narrowing without compromising compliance.

Extracellular signal-regulated kinase inhibition prevents venous adaptive remodeling via regulation of Eph-B4

OBJECTIVES

Vein graft adaptation (VGA) is a process that vein as a vascular graft conduits in arterial reconstructive surgery; VGA can lead to postoperative vein graft stenosis (VGS) and complications after coronary artery bypass graft and other peripheral artery bypass surgeries. VGA is characterized by vein graft loss the venous features without exhibiting arterial features; furthermore, the activation of ERK inhibited the maintenance of venous properties of the vein graft. We hypothesized that ERK inhibition can affect vein VGS through regulating the expression of EphB4.

METHODS

Rat vein transplantation model was established using wild-type and EphB4^+/- Sprague-Dawley rats. Hematoxylin-eosin, Masson, Verhoeff, actin staining, and immunohistochemistry were applied to observe the structure of the vein grafts. Vascular smooth muscle cells (VSMCs) were isolated from the vein and vein grafts. Western blotting was used to determine the expression of p-ERK1/2 and EphB4, and immunofluorescence was applied to detect the expression and location of EphB4. Cell wound scratch assay and CCK8 assay were used to determine the migration and proliferation of VSMCs. Real-time polymerase chain reaction was used to determine the mRNA expression of EphB4.

RESULTS

Western blotting in vein sample and vein graft sample detected p-ERK1/2 and ERK1/2 expression in both EphB4^+/+ and EphB4^+/- rats. The expression of p-ERK was increased in vein graft compared to vein. Immunofluorescence in VSMCs form EphB4^+/+ and EphB4^+/- rats detected EphB4 expression in both cells, and the expression of EphB4 was increased in VSMCs form EphB4^+/+ rats. SCH772984 reduces the proliferation and migration of VSMCs. Inhibition of ERK suppressed the increase of vein graft wall thickness, and the expression of collagen fibers, elastic fibers, and α-actin was decreased. Vein graft from EphB4^+/- rats reduces the expression of EphB4, and SCH772984 suppressed the decrease of EphB4 in vivo. Vein graft from EphB4^+/- rats increased the expression of EphB4, and SCH772984 suppressed the increase of EphB4 in vivo.

CONCLUSIONS

The inhibition of ERK1/2 suppressed the process of VGS by decreasing the proliferation of VSMCs. The ERK-inhibitor SCH772984 suppressed the level of VGS by extending the time of EphB4 expression during the process of VGA, thus maintaining the venousization of vein graft. The mechanism may be that the inhibitor SCH772984 suppresses the level of VGS by extending the time of EphB4 expression during the process of VGA. Therefore, our research provides a new target of VGS treatment by inhibiting the expression of ERK1/2 through the process of VGA.

A Twofold Usage of an Agent-Based Model of Vascular Adaptation to Design Clinical Experiments

Several computational models of Vein Graft Bypass (VGB) adaptation have been developed in order to improve the surgical outcome and they all share a common property: their accuracy relies on a winning choice of their driving coefficients which are best to be retrieved from experimental data. Since experiments are time-consuming and resources-demanding, the golden standard is to know in advance which measures need to be retrieved on the experimental table and out of how many samples. Accordingly, our goal is to build a computational framework able to pre-design an effective experimental structure to optimize the computational models setup. Our hypothesis is that an Agent-Based Model (ABM) developed by our group is comparable enough to a true set of experiments to be used to generate reliable virtual experimental data. Thanks to a twofold usage of our ABM, we created a filter to be posed before the real experiment in order to drive its optimal design. This work is the natural continuation of a previous study from our group [1], where the attention was posed on simple single-cellular events models. With this new version we focused on more complex models with the purpose of verifying that the complexity of the experimental setup grows proportionally with the accuracy of the model itself.

A versatile hybrid agent-based, particle and partial differential equations method to analyze vascular adaptation

Peripheral arterial occlusive disease is a chronic pathology affecting at least 8–12 million people in the USA, typically treated with a vein graft bypass or through the deployment of a stent in order to restore the physiological circulation. Failure of peripheral endovascular interventions occurs at the intersection of vascular biology, biomechanics, and clinical decision making. It is our hypothesis that the majority of endovascular treatment approaches share the same driving mechanisms and that a deep understanding of the adaptation process is pivotal in order to improve the current outcome of the procedure. The postsurgical adaptation of vein graft bypasses offers the perfect example of how the balance between intimal hyperplasia and wall remodeling determines the failure or the success of the intervention. Accordingly, this work presents a versatile computational model able to capture the feedback loop that describes the interaction between events at cellular/tissue level and mechano-environmental conditions. The work here presented is a generalization and an improvement of a previous work by our group of investigators, where an agent-based model uses a cellular automata principle on a fixed hexagonal grid to reproduce the leading events of the graft’s restenosis. The new hybrid model here presented allows a more realistic simulation both of the biological laws that drive the cellular behavior and of the active role of the membranes that separate the various layers of the vein. The novel feature is to use an immersed boundary implementation of a highly viscous flow to represent SMC motility and matrix reorganization in response to graft adaptation. Our implementation is modular, and this makes us able to choose the right compromise between closeness to the physiological reality and complexity of the model. The focus of this paper is to offer a new modular implementation that combines the best features of an agent-based model, continuum mechanics, and particle-tracking methods to cope with the multiscale nature of the adaptation phenomena. This hybrid method allows us to quickly test various hypotheses with a particular attention to cellular motility, a process that we demonstrated should be driven by mechanical homeostasis in order to maintain the right balance between cells and extracellular matrix in order to reproduce a distribution similar to histological experimental data from vein grafts.

A two-phase model of early fibrous cap formation in atherosclerosis

Atherosclerotic plaque growth is characterised by chronic, non-resolving inflammation that promotes the accumulation of cellular debris and extracellular fat in the inner artery wall. This material is highly thrombogenic, and plaque rupture can lead to the formation of blood clots that occlude major arteries and cause myocardial infarction or stroke. In advanced plaques, vascular smooth muscle cells (SMCs) are recruited from deeper in the artery wall to synthesise a cap of fibrous tissue that stabilises the plaque and sequesters the thrombogenic plaque content from the bloodstream. The fibrous cap provides crucial protection against the clinical consequences of atherosclerosis, but the mechanisms of cap formation are poorly understood. In particular, it is unclear why certain plaques become stable and robust while others become fragile and dangerously vulnerable to rupture. We develop a multiphase model with non-standard boundary conditions to investigate early fibrous cap formation in the atherosclerotic plaque. The model is parameterised using data from a range of in vitro and in vivo studies, and includes highly nonlinear mechanisms of SMC proliferation and migration in response to an endothelium-derived chemical signal. We demonstrate that the model SMC population naturally evolves towards a steady-state, and predict a rate of cap formation and a final plaque SMC content consistent with experimental observations in mice. Parameter sensitivity simulations show that SMC proliferation makes a limited contribution to cap formation, and demonstrate that stable cap formation relies primarily on a critical balance between the rates of SMC recruitment to the plaque, chemotactic SMC migration within the plaque and SMC loss by apoptosis or phenotype change. This model represents the first detailed in silico study of fibrous cap formation in atherosclerosis, and establishes a multiphase modelling framework that can be readily extended to investigate many other aspects of plaque development.

A Computational Model-Based Framework to Plan Clinical Experiments - an Application to Vascular Adaptation Biology

Several computational models have been developed in order to improve the outcome of Vein Graft Bypasses in response to arterial occlusions and they all share a common property: their accuracy relies on a winning choice of the coefficients' value related to biological functions that drive them. Our goal is to optimize the retrieval of these unknown coefficients on the base of experimental data and accordingly, as biological experiments are noisy in terms of statistical analysis and the models are typically stochastic and complex, this work wants first to elucidate which experimental measurements might be sufficient to retrieve the targeted coefficients and second how many specimens would constitute a good dataset to guarantee a sufficient level of accuracy. Since experiments are often costly and time consuming, the planning stage is critical to the success of the operation and, on the base of this consideration, the present work shows how, thanks to an ad hoc use of a computational model of vascular adaptation, it is possible to estimate in advance the entity and the quantity of resources needed in order to efficiently reproduce the experimental reality.

Linking gene dynamics to vascular hyperplasia - Toward a predictive model of vein graft adaptation

Reductionist approaches, where individual pieces of a process are examined in isolation, have been the mainstay of biomedical research. While these methods are effective in highly compartmentalized systems, they fail to account for the inherent plasticity and non-linearity within the signaling structure. In the current manuscript, we present the computational architecture for tracking an acute perturbation in a biologic system through a multiscale model that links gene dynamics to cell kinetics, with the overall goal of predicting tissue adaptation. Given the complexity of the genome, the problem is made tractable by clustering temporal changes in gene expression into unique patterns. These cluster elements form the core of an integrated network that serves as the driving force for the response of the biologic system. This modeling approach is illustrated using the clinical scenario of vein bypass graft adaptation. Vein segments placed in the arterial circulation for treatment of advanced occlusive disease can develop an aggressive hyperplastic response that narrows the lumen, reduces blood flow, and induces in situ thrombosis. Reducing this hyperplastic response has been a long-standing but unrealized goal of biologic researchers in the field. With repeated failures of single target therapies, the redundant response pathways are thought to be a fundamental issue preventing progress towards a solution. Using the current framework, we demonstrate how theoretical genomic manipulations can be introduced into the system to shift the adaptation to a more beneficial phenotype, where the hyperplastic response is mitigated and the risk of thrombosis reduced. Utilizing our previously published rabbit vein graft genomic data, where grafts were harvested at time points ranging from 2 hours to 28 days and under differential flow conditions, and a customized clustering algorithm, five gene clusters that differentiated the low flow (i.e., pro-hyperplastic) from high flow (i.e., anti-hyperplastic) response were identified. The current analysis advances these general associations to create a model that identifies those genes sets most likely to be of therapeutic benefit. Using this approach, we examine the range of potential opportunities for intervention via gene cluster over-expression or inhibition, delivered in isolation or combination, at the time of vein graft implantation.

Limiting Injury During Saphenous Vein Graft Preparation For Coronary Arterial Bypass Prevents Metabolic Decompensation

Standard harvest and preparation of human saphenous vein (HSV) for autologous coronary and peripheral arterial bypass procedures is associated with injury and increased oxidative stress that negatively affect graft performance. In this study we investigated the global metabolomic profiles of HSV before (unprepared; UP) and after standard vein graft preparation (AP). AP-HSV showed impaired vasomotor function that was associated with increased oxidative stress, phospholipid hydrolysis and energy depletion that are characteristic of mechanical and chemical injury. A porcine model (PSV) was utilized to validate these metabolomic changes in HSV and to determine the efficacy of an improved preparation technique (OP) using pressure-regulated distension, a non-toxic vein marker, and graft storage in buffered PlasmaLyte solution in limiting metabolic decompensation due to graft preparation. Deficits in vasomotor function and metabolic signature observed in AP-PSV could be largely mitigated with the OP procedure. These findings suggest that simple strategies aimed at reducing injury during graft harvest and preparation represents a straightforward and viable strategy to preserve conduit function and possibly improve graft patency.

Vascular Adaptation: Pattern Formation and Cross Validation between an Agent Based Model and a Dynamical System

Myocardial infarction is the global leading cause of mortality (Go et al., 2014). Coronary artery occlusion is its main etiology and it is commonly treated by Coronary Artery Bypass Graft (CABG) surgery (Wilson et al, 2007). The long-term outcome remains unsatisfactory (Benedetto, 2016) as the graft faces the phenomenon of restenosis during the post-surgery, which consists of re-occlusion of the lumen and usually requires secondary intervention even within one year after the initial surgery (Harskamp, 2013). In this work, we propose an extensive study of the restenosis phenomenon by implementing two mathematical models previously developed by our group: a heuristic Dynamical System (DS) (Garbey and Berceli, 2013), and a stochastic Agent Based Model (ABM) (Garbey et al., 2015). With an extensive use of the ABM, we retrieved the pattern formations of the cellular events that mainly lead the restenosis, especially focusing on mitosis in intima, caused by alteration in shear stress, and mitosis in media, fostered by alteration in wall tension. A deep understanding of the elements at the base of the restenosis is indeed crucial in order to improve the final outcome of vein graft bypass. We also turned the ABM closer to the physiological reality by abating its original assumption of circumferential symmetry. This allowed us to finely replicate the trigger event of the restenosis, i.e. the loss of the endothelium in the early stage of the post-surgical follow up (Roubos et al., 1995) and to simulate the encroachment of the lumen in a fashion aligned with histological evidences (Owens et al., 2015). Finally, we cross-validated the two models by creating an accurate matching procedure. In this way we added the degree of accuracy given by the ABM to a simplified model (DS) that can serve as powerful predictive tool for the clinic.

Resolution of vascular injury: Specialized lipid mediators and their evolving therapeutic implications

Acute vascular injury occurs in a number of important clinical contexts, including spontaneous disease-related events (e.g. plaque rupture, thrombosis) and therapeutic interventions such as angioplasty, stenting, or bypass surgery. Endothelial cell (EC) disruption exposes the underlying matrix, leading to a rapid deposition of platelets, coagulation proteins, and leukocytes. A thrombo-inflammatory response ensues characterized by leukocyte recruitment, vascular smooth muscle cell (VSMC) activation, and the elaboration of cytokines, reactive oxygen species and growth factors within the vessel wall. A resolution phase of vascular injury may be described in which leukocyte efflux, clearance of debris, and re-endothelialization occurs. VSMC migration and proliferation leads to the development of a thickened neointima that may lead to lumen compromise. Subsequent remodeling involves matrix protein deposition, and return of EC and VSMC to quiescence. Recent studies suggest that specialized pro-resolving lipid mediators (SPM) modulate key aspects of this response, and may constitute an endogenous homeostatic pathway in the vasculature. SPM exert direct effects on vascular cells that counteract inflammatory signals, reduce leukocyte adhesion, and inhibit VSMC migration and proliferation. These effects appear to be largely G-protein coupled receptor-dependent. Across a range of animal models of vascular injury, including balloon angioplasty, bypass grafting, and experimental aneurysm formation, SPM accelerate repair and reduce lesion formation. With bioactivity in the pM-nM range, a lack of discernible cytotoxicity, and a spectrum of vasculo-protective properties, SPM represent a novel class of vascular therapeutics. This review summarizes current research in this field, including a consideration of critical next steps and challenges in translation.

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.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

The correlation between computed tomography and duplex evaluation of autogenous vein bypass grafts and their relationship to failure

OBJECTIVE

Duplex ultrasound (DUS) imaging for vein bypass graft (VBG) surveillance is confounded by technical and physiologic factors that reduce the sensitivity for detecting impending graft failure. In contrast, three-dimensional computed tomography angiography (CTA) offers high-fidelity anatomic characterization of VBGs, but its utility in detecting at risk grafts is unknown. The current study analyzed the correlation between DUS and CTA for detection of vein graft stenosis and evaluated the relationship of the observed abnormalities to VBG failure.

METHODS

Consecutive lower extremity VBG patients underwent surveillance with concurrent DUS imaging and CTA at 1 week and at 1, 6, and 12 months postoperatively. A standardized algorithm was used for CT reconstruction and extraction of the lumen geometries at 1-mm intervals. At each interval, CT-derived cross-sectional areas were coregistered and correlated to DUS peak systolic velocities (PSVs) within six predesignated anatomic zones and then analyzed for outcome association. Vein graft failure was defined as pathologic change within a given anatomic zone resulting in thrombosis, amputation, or reintervention within the 6-month period after the observed time point.

RESULTS

The study recruited 54 patients, and 10 (18%) experienced failure ≤18 months of implantation. The expected inverse relationship between cross-sectional area and PSV was only weakly correlated (Spearman rank coefficient = -0.19). Moderate elevations in the PSV ratio (PSVr; 2-3.5) were frequently transient, with 14 of 18 grafts (78%) demonstrating ratio reduction on subsequent imaging. A PSVr ≥3.5 was associated with a 67% failure rate. CT stenosis <50% was highly correlated with success (0 failures); however, high-grade (>80%) CT stenosis was more likely to succeed than to fail (25%). Significant discordance between CT and DUS was found in 18 patients. Although 14 of these patients had CT stenosis >70% with a PSVr <3.5, subsequent failure occurred in only two. Conversely, graft failure occurred in three of four patients with CT stenosis <70% but PSVr >3.5. Focused analysis of these patients using computational fluid dynamic modeling demonstrated that vein side branches, local tortuosity, regional diameter variations, and venovenostomies were the drivers of these discrepancies.

CONCLUSIONS

This analysis demonstrated that a PSVr ≥3.5 is strongly correlated with VBG failure, whereas the natural history of moderately elevated PSVr (2-3.5) is largely clinically benign. Although minimum stenosis on the CT scan was highly predictive of success, high-grade CT stenosis was infrequently associated with failure. The interaction of anatomic features with the local flow dynamics was identified as the primary confounder for a direct correlation between CT and DUS imaging.

A Multiscale Computational Framework to Understand Vascular Adaptation

The failure rate for vascular interventions (vein bypass grafting, arterial angioplasty/stenting) remains unacceptably high. Over the past two decades, researchers have applied a wide variety of approaches to investigate the primary failure mechanisms, neointimal hyperplasia and aberrant remodeling of the wall, in an effort to identify novel therapeutic strategies. Despite incremental progress, specific cause/effect linkages among the primary drivers of the pathology, (hemodynamic factors, inflammatory biochemical mediators, cellular effectors) and vascular occlusive phenotype remain lacking. We propose a multiscale computational framework of vascular adaptation to develop a bridge between theory and experimental observation and to provide a method for the systematic testing of relevant clinical hypotheses. Cornerstone to our model is a feedback mechanism between environmental conditions and dynamic tissue plasticity described at the cellular level with an agent based model. Our implementation (i) is modular, (ii) starts from basic mechano-biology principle at the cell level and (iii) facilitates the agile development of the model.