Rule-based model of vein graft remodeling

A predictive multiscale model of in-stent restenosis in femoral arteries: linking haemodynamics and gene expression with an agent-based model of cellular dynamics

In-stent restenosis (ISR) is a maladaptive inflammatory-driven response of femoral arteries to percutaneous transluminal angioplasty and stent deployment, leading to lumen re-narrowing as consequence of excessive cellular proliferative and synthetic activities. A thorough understanding of the underlying mechanobiological factors contributing to ISR is still lacking. Computational multiscale models integrating both continuous- and agent-based approaches have been identified as promising tools to capture key aspects of the complex network of events encompassing molecular, cellular and tissue response to the intervention. In this regard, this work presents a multiscale framework integrating the effects of local haemodynamics and monocyte gene expression data on cellular dynamics to simulate ISR mechanobiological processes in a patient-specific model of stented superficial femoral artery. The framework is based on the coupling of computational fluid dynamics simulations (haemodynamics module) with an agent-based model (ABM) of cellular activities (tissue remodelling module). Sensitivity analysis and surrogate modelling combined with genetic algorithm optimization were adopted to explore the model behaviour and calibrate the ABM parameters. The proposed framework successfully described the patient lumen area reduction from baseline to one-month follow-up, demonstrating the potential capabilities of this approach in predicting the short-term arterial response to the endovascular procedure.

Multiscale Modeling of Vascular Remodeling Induced by Wall Shear Stress

Objective

Hemodynamics-induced low wall shear stress (WSS) is one of the critical reasons leading to vascular remodeling. However, the coupling effects of WSS and cellular kinetics have not been clearly modeled. The aim of this study was to establish a multiscale modeling approach to reveal the vascular remodeling behavior under the interaction between the macroscale of WSS loading and the microscale of cell evolution.

Methods

Computational fluid dynamics (CFD) method and agent-based model (ABM), which have significantly different characteristics in temporal and spatial scales, were adopted to establish the multiscale model. The CFD method is for the second/organ scale, and the ABM is for the month/cell scale. The CFD method was used to simulate blood flow in a vessel and obtain the WSS in a vessel cross-section. The simulations of the smooth muscle cell (SMC) proliferation/apoptosis and extracellular matrix (ECM) generation/degradation in a vessel cross-section were performed by using ABM. During the simulation of the vascular remodeling procedure, the damage index of the SMC and ECM was defined as deviation from the obtained WSS. The damage index decreased gradually to mimic the recovery of WSS-induced vessel damage.

Results

(1) The significant wall thickening region was consistent with the low WSS region. (2) There was no evident change of wall thickness in the normal WSS region. (3) When the damage index approached to 0, the amount and distribution of SMCs and ECM achieved a stable state, and the vessel reached vascular homeostasis.

Conclusion

The established multiscale model can be used to simulate the vascular remodeling behavior over time under various WSS conditions.

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.

Monitoring the Impact of Aggressive Glycemic Intervention during Critical Care after Cardiac Surgery with a Glycemic Expert System for Nurse-Implemented Euglycemia: The MAGIC GENIE Project

A novel, multi-dimensional protocol named GENIE has been in use for intensive insulin therapy (IIT, target glucose <140 mg/dL) in the surgical intensive care unit (SICU) after open heart surgery (OHS) at VA Pittsburgh since 2005. Despite concerns over increased mortality from IIT after the publication of the NICE-SUGAR Trial, it remains in use, with ongoing monitoring under the MAGIC GENIE Project showing that GENIE performance over 12 years (2005-2016) aligns with the current consensus that IIT with target blood glucose (BG) <140 mg/dL is advisable only if it does not provoke severe hypoglycemia (SH). Two studies have been conducted to monitor glucometrics and outcomes during GENIE use in the SICU. One compares GENIE (n = 382) with a traditional IIT protocol (FORMULA, n = 289) during four years of contemporaneous use (2005-2008). The other compares GENIE's impact overall (n = 1404) with a cohort of patients who maintained euglycemia after OHS (euglycemic no-insulin [ENo-I], n = 111) extending across 12 years (2005-2016). GENIE performed significantly better than FORMULA during contemporaneous use, maintaining lower time-averaged glucose, provoking less frequent, severe, prolonged, or repetitive hypoglycemia, and achieving 50% lower one-year mortality, with no deaths from mediastinitis (0 of 8 cases vs 4 of 9 on FORMULA). Those benefits were sustained over the subsequent eight years of exclusive use in OHS patients, with an overall one-year mortality rate (4.2%) equivalent to the ENo-I cohort (4.5%). The results of the MAGIC GENIE Project show that GENIE can maintain tight glycemic control without provoking SH in patients undergoing OHS, and may be associated with a durable survival benefit. The results, however, await confirmation in a randomized control trial.

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.

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.

An agent-based model of leukocyte transendothelial migration during atherogenesis

A vast amount of work has been dedicated to the effects of hemodynamics and cytokines on leukocyte adhesion and trans-endothelial migration (TEM) and subsequent accumulation of leukocyte-derived foam cells in the artery wall. However, a comprehensive mechanobiological model to capture these spatiotemporal events and predict the growth and remodeling of an atherosclerotic artery is still lacking. Here, we present a multiscale model of leukocyte TEM and plaque evolution in the left anterior descending (LAD) coronary artery. The approach integrates cellular behaviors via agent-based modeling (ABM) and hemodynamic effects via computational fluid dynamics (CFD). In this computational framework, the ABM implements the diffusion kinetics of key biological proteins, namely Low Density Lipoprotein (LDL), Tissue Necrosis Factor alpha (TNF-α), Interlukin-10 (IL-10) and Interlukin-1 beta (IL-1β), to predict chemotactic driven leukocyte migration into and within the artery wall. The ABM also considers wall shear stress (WSS) dependent leukocyte TEM and compensatory arterial remodeling obeying Glagov's phenomenon. Interestingly, using fully developed steady blood flow does not result in a representative number of leukocyte TEM as compared to pulsatile flow, whereas passing WSS at peak systole of the pulsatile flow waveform does. Moreover, using the model, we have found leukocyte TEM increases monotonically with decreases in luminal volume. At critical plaque shapes the WSS changes rapidly resulting in sudden increases in leukocyte TEM suggesting lumen volumes that will give rise to rapid plaque growth rates if left untreated. Overall this multi-scale and multi-physics approach appropriately captures and integrates the spatiotemporal events occurring at the cellular level in order to predict leukocyte transmigration and plaque evolution.

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.

Vein graft failure: from pathophysiology to clinical outcomes

Occlusive arterial disease is a leading cause of morbidity and mortality worldwide. Aside from balloon angioplasty, bypass graft surgery is the most commonly performed revascularization technique for occlusive arterial disease. Coronary artery bypass graft surgery is performed in patients with left main coronary artery disease and three-vessel coronary disease, whereas peripheral artery bypass graft surgery is used to treat patients with late-stage peripheral artery occlusive disease. The great saphenous veins are commonly used conduits for surgical revascularization; however, they are associated with a high failure rate. Therefore, preservation of vein graft patency is essential for long-term surgical success. With the exception of 'no-touch' techniques and lipid-lowering and antiplatelet (aspirin) therapy, no intervention has hitherto unequivocally proven to be clinically effective in preventing vein graft failure. In this Review, we describe both preclinical and clinical studies evaluating the pathophysiology underlying vein graft failure, and the latest therapeutic options to improve patency for both coronary and peripheral grafts.

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.

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.

Biomechanics and mechanobiology in functional tissue engineering

The field of tissue engineering continues to expand and mature, and several products are now in clinical use, with numerous other preclinical and clinical studies underway. However, specific challenges still remain in the repair or regeneration of tissues that serve a predominantly biomechanical function. Furthermore, it is now clear that mechanobiological interactions between cells and scaffolds can critically influence cell behavior, even in tissues and organs that do not serve an overt biomechanical role. Over the past decade, the field of "functional tissue engineering" has grown as a subfield of tissue engineering to address the challenges and questions on the role of biomechanics and mechanobiology in tissue engineering. Originally posed as a set of principles and guidelines for engineering of load-bearing tissues, functional tissue engineering has grown to encompass several related areas that have proven to have important implications for tissue repair and regeneration. These topics include measurement and modeling of the in vivo biomechanical environment; quantitative analysis of the mechanical properties of native tissues, scaffolds, and repair tissues; development of rationale criteria for the design and assessment of engineered tissues; investigation of the effects biomechanical factors on native and repair tissues, in vivo and in vitro; and development and application of computational models of tissue growth and remodeling. Here we further expand this paradigm and provide examples of the numerous advances in the field over the past decade. Consideration of these principles in the design process will hopefully improve the safety, efficacy, and overall success of engineered tissue replacements.

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.