Quantification of patient-specific coronary material properties and their correlations with plaque morphological characteristics: An in vivo IVUS study.

BACKGROUND: A method using in vivo Cine IVUS and VH-IVUS data has been proposed to quantify material properties of coronary plaques. However, correlations between plaque morphological characteristics and mechanical properties have not been studied in vivo. METHOD: In vivo Cine IVUS and VH-IVUS data were acquired at 32 plaque cross-sections from 19 patients. Six morphological factors were extracted for each plaque. These samples were categorized into healthy vessel, fibrous plaque, lipid-rich plaque and calcified plaque for comparisons. Three-dimensional thin-slice models were constructed using VH-IVUS data to quantify in vivo plaque material properties following a finite element updating approach by matching Cine IVUS data. Effective Young's moduli were calculated to represent plaque stiffness for easy comparison. Spearman's rank correlation analysis was performed to identify correlations between plaque stiffness and morphological factor. Kruskal-Wallis test with Bonferroni correction was used to determine whether significant differences in plaque stiffness exist among four plaque groups. RESULT: Our results show that lumen circumference change has a significantly negative correlation with plaque stiffness (r = -0.7807, p = 0.0001). Plaque burden and calcification percent also had significant positive correlations with plaque stiffness (r = 0.5105, p < 0.0272 and r = 0.5312, p < 0.0193) respectively. Among the four categorized groups, calcified plaques had highest stiffness while healthy segments had the lowest. CONCLUSION: There is a close link between plaque morphological characteristics and mechanical properties in vivo. Plaque stiffness tends to be higher as coronary atherosclerosis advances, indicating the potential to assess plaque mechanical properties in vivo based on plaque compositions.

Optimization of Left Ventricle Pace Maker Location Using Echo-Based Fluid-Structure Interaction Models.

INTRODUCTION: Cardiac pacing has been an effective treatment in the management of patients with bradyarrhythmia and tachyarrhythmia. Different pacemaker location has different responses, and pacemaker effectiveness to each individual can also be different. A novel image-based ventricle animal modeling approach was proposed to optimize ventricular pacemaker site for better cardiac outcome. METHOD: One health female adult pig (weight 42.5 kg) was used to make a pacing animal model with different ventricle pacing locations. Ventricle surface electric signal, blood pressure and echo image were acquired 15 min after the pacemaker was implanted. Echo-based left ventricle fluid-structure interaction models were constructed to perform ventricle function analysis and investigate impact of pacemaker location on cardiac outcome. With the measured electric signal map from the pig associated with the actual pacemaker site, electric potential conduction of myocardium was modeled by material stiffening and softening in our model, with stiffening simulating contraction and softening simulating relaxation. Ventricle model without pacemaker (NP model) and three ventricle models with the following pacemaker locations were simulated: right ventricular apex (RVA model), posterior interventricular septum (PIVS model) and right ventricular outflow tract (RVOT model). Since higher peak flow velocity, flow shear stress (FSS), ventricle stress and strain are linked to better cardiac function, those data were collected for model comparisons. RESULTS: At the peak of filling, velocity magnitude, FSS, stress and strain for RVOT and PIVS models were 13%, 45%, 18%, 13% and 5%, 30%, 10%, 5% higher than NP model, respectively. At the peak of ejection, velocity magnitude, FSS, stress and strain for RVOT and PIVS models were 50%, 44%, 54%, 59% and 23%, 36%, 39%, 53% higher than NP model, respectively. RVA model had lower velocity, FSS, stress and strain than NP model. RVOT model had higher peak flow velocity and stress/strain than PIVS model. It indicated RVOT pacemaker site may be the best location. CONCLUSION: This preliminary study indicated that RVOT model had the best performance among the four models compared. This modeling approach could be used as "virtual surgery" to try various pacemaker locations and avoid risky and dangerous surgical experiments on real patients.

Design, Synthesis, and Characterization of I-BET567, a Pan-Bromodomain and Extra Terminal (BET) Bromodomain Oral Candidate.

Through regulation of the epigenome, the bromodomain and extra terminal (BET) family of proteins represent important therapeutic targets for the treatment of human disease. Through mimicking the endogenous N-acetyl-lysine group and disrupting the protein-protein interaction between histone tails and the bromodomain, several small molecule pan-BET inhibitors have progressed to oncology clinical trials. This work describes the medicinal chemistry strategy and execution to deliver an orally bioavailable tetrahydroquinoline (THQ) pan-BET candidate. Critical to the success of this endeavor was a potency agnostic analysis of a data set of 1999 THQ BET inhibitors within the GSK collection which enabled identification of appropriate lipophilicity space to deliver compounds with a higher probability of desired oral candidate quality properties. SAR knowledge was leveraged via Free-Wilson analysis within this design space to identify a small group of targets which ultimately delivered I-BET567 (27), a pan-BET candidate inhibitor that demonstrated efficacy in mouse models of oncology and inflammation.

A Novel Pulmonary Valve Replacement Surgery Strategy Using Contracting Band for Patients With Repaired Tetralogy of Fallot: An MRI-Based Multipatient Modeling Study.

Patients with repaired Tetralogy of Fallot (ToF), a congenital heart defect which includes a ventricular septal defect and severe right ventricular outflow obstruction, account for the majority of cases with late-onset right ventricle (RV) failure. Current surgery procedures, including pulmonary valve replacement (PVR) with right ventricle remodeling, yield mixed results. PVR with active band insertion was hypothesized to be of clinical usage on improving RV function measured by ejection fraction (EF). In lieu of risky open-heart surgeries and experiments on animal and human, computational biomechanical models were adapted to study the impact of PVR with five band insertion options. Cardiac magnetic resonance (CMR) images were acquired from seven TOF patients before PVR surgery for model construction. For each patient, five different surgery plans combined with passive and active contraction band with contraction ratio of 20, 15, and 10% were studied. Those five plans include three single-band plans with different band locations; one plan with two bands, and one plan with three bands. Including the seven no-band models, 147 computational bi-ventricle models were constructed to simulate RV cardiac functions and identify optimal band plans. Patient variations with different band plans were investigated. Surgery plan with three active contraction bands and band active contraction ratio of 20% had the best performance on improving RV function. The mean ± SD RV ejection fraction value from the seven patients was 42.90 ± 5.68%, presenting a 4.19% absolute improvement or a 10.82% relative improvement, when compared with the baseline models (38.71 ± 5.73%, p = 0.016). The EF improvements from the seven patients varied from 2.87 to 6.01%. Surgical procedures using active contraction bands have great potential to improve RV function measured by ejection fraction for patients with repaired ToF. It is possible to have higher right ventricle ejection fraction improvement with more bands and higher band active contraction ratio. Our findings with computational models need to be further validated by animal experiments before clinical trial could become possible.

Comparisons of simulation results between passive and active fluid structure interaction models for left ventricle in hypertrophic obstructive cardiomyopathy.

BACKGROUND: Patient-specific active fluid-structure interactions (FSI) model is a useful approach to non-invasively investigate the hemodynamics in the heart. However, it takes a lot of effort to obtain the proper external force boundary conditions for active models, which heavily restrained the time-sensitive clinical applications of active computational models. METHODS: The simulation results of 12 passive FSI models based on 6 patients' pre-operative and post-operative CT images were compared with corresponding active models to investigate the differences in hemodynamics and cardiac mechanics between these models. RESULTS: In comparing the passive and active models, it was found that there was no significant difference in pressure difference and shear stress on mitral valve leaflet (MVL) at the pre-SAM time point, but a significant difference was found in wall stress on the inner boundary of left ventricle (endocardium). It was also found that pressure difference on the coapted MVL and the shear stress on MVL were significantly decreased after successful surgery in both active and passive models. CONCLUSION: Our results suggested that the passive models may provide good approximated hemodynamic results at 5% RR interval, which is crucial for analyzing the initiation of systolic anterior motion (SAM). Comparing to active models, the passive models decrease the complexity of the modeling construction and the difficulty of convergence significantly. These findings suggest that, with proper boundary conditions and sufficient clinical data, the passive computational model may be a good substitution model for the active model to perform hemodynamic analysis of the initiation of SAM.

Multi-Band Surgery for Repaired Tetralogy of Fallot Patients With Reduced Right Ventricle Ejection Fraction: A Pilot Study.

INTRODUCTION: Right ventricle (RV) failure is one of the most common symptoms among patients with repaired tetralogy of Fallot (TOF). The current surgery treatment approach including pulmonary valve replacement (PVR) showed mixed post-surgery outcomes. A novel PVR surgical strategy using active contracting bands is proposed to improve the post-PVR outcome. In lieu of testing the risky surgical procedures on real patients, computational simulations (virtual surgery) using biomechanical ventricle models based on patient-specific cardiac magnetic resonance (CMR) data were performed to test the feasibility of the PVR procedures with active contracting bands. Different band combination and insertion options were tested to identify optimal surgery designs. METHOD: Cardiac magnetic resonance data were obtained from one TOF patient (male, age 23) whose informed consent was obtained. A total of 21 finite element models were constructed and solved following our established procedures to investigate the outcomes of the band insertion surgery. The non-linear anisotropic Mooney-Rivlin model was used as the material model. Five different band insertion plans were simulated (three single band models with different band locations, one model with two bands, and one model with three bands). Three band contraction ratios (10, 15, and 20%) and passive bands (0% contraction ratio) were tested. RV ejection fraction was used as the measure for cardiac function. RESULTS: The RV ejection fraction from the three-band model with 20% contraction increased to 41.58% from the baseline of 37.38%, a 4.20% absolute improvement. The RV ejection fractions from the other four band models with 20% contraction rate were 39.70, 39.45, and 40.70% (two-band) and 39.17%, respectively. The mean RV stress and strain values from all of the 21 models showed only modest differences (5-11%). CONCLUSION: This pilot study demonstrated that the three-band model with 20% band contraction ratio led to 4.20% absolute improvement in the RV ejection fraction, which is considered as clinically significant. The passive elastic bands led to the reduction of the RV ejection fractions. The modeling results and surgical strategy need to be further developed and validated by a multi-patient study and animal experiments before clinical trial could become possible. Tissue regeneration techniques are needed to produce materials for the contracting bands.

Patient-specific in vivo right ventricle material parameter estimation for patients with tetralogy of Fallot using MRI-based models with different zero-load diastole and systole morphologies.

Patient-specific in vivo ventricle material parameter determination is important for cardiovascular investigations. A new cardiac magnetic image (CMR)-based modeling approach with different zero-load diastole and systole geometries was adopted to estimate right ventricle material parameter values for healthy and patients with Tetralogy of Fallot (TOF) and seeking potential clinical applications. CMR data were obtained from 6 healthy volunteers and 16 TOF patients with consent obtained. CMR-based RV/LV models were constructed using two zero-load geometries (diastole and systole, 2G model). Material parameter values for begin-filling (BF), end-filling (EF), begin-ejection (BE), and end-ejection (EE) were recorded for analyses. Effective Young's moduli (YM) for fiber direction stress-strain curves were calculated for easy comparisons. The mean EE YM value of TOF patients was 78.6% higher than that of the healthy group (HG). The mean end-ejection YM value from worse-outcome TOF group (WG) post pulmonary valve replacement (PVR) surgery was 59.5% higher than that from the better-outcome TOF group (BG). Using begin-filling YM and end-ejection YM as predictors and the classic logistic regression model to different better-outcome group patients from worse-outcome group patients, the areas under Receiver Operating Characteristic (ROC) curves were found to be 0.797 and 0.883 for begin-filling YM and end-ejection YM, respectively. The sensitivity and specificity 0.761 and 0.755 using end-ejection YM as the predictor. This preliminary study suggests that ventricle material stiffness could be a potential parameter to be used to differentiate BG patients from WG patients with further effort and large-scale patient data validations.

Numerical simulation study on systolic anterior motion of the mitral valve in hypertrophic obstructive cardiomyopathy.

BACKGROUND: The hydrodynamic mechanisms of systolic anterior motion (SAM) of the mitral valve in hypertrophic obstructive cardiomyopathy (HOCM) remain unclear. METHODS: Based on computed tomography (CT) images and clinical data, pre- and post-operative computational models of the left ventricle were constructed for 6 HOCM patients receiving septal myectomy. SAM was abolished in 5 patients and persisted in one after septal myectomy surgery. The obtained simulation results including flow field of the left ventricle and mechanical behaviors of the mitral valve (MV) between pre- and post-operative FSI models were compared. RESULTS: The pressure difference and shear stress on the mitral valve leaflets (MVL) were relatively high pre-operatively, and decreased significantly after satisfactory surgery, but remained high following failed surgery. The significant increase in coaptation-to-septal distance was found when SAM was abolished. CONCLUSIONS: Our results indicated that high pressure difference and shear stress on the MVL might directly initiate SAM in HOCM. Successful septal myectomy enlarged the coaptation-to-septal distance sufficiently to keep the MVL away from the ejection flow, thereby eliminating SAM.

Effects of Residual Stress, Axial Stretch, and Circumferential Shrinkage on Coronary Plaque Stress and Strain Calculations: A Modeling Study Using IVUS-Based Near-Idealized Geometries.

Accurate stress and strain calculations are important for plaque progression and vulnerability assessment. Models based on in vivo data often need to form geometries with zero-stress/strain conditions. The goal of this paper is to use IVUS-based near-idealized geometries and introduce a three-step model construction process to include residual stress, axial shrinkage, and circumferential shrinkage and investigate their impacts on stress and strain calculations. In Vivo intravascular ultrasound (IVUS) data of human coronary were acquired for model construction. In Vivo IVUS movie data were acquired and used to determine patient-specific material parameter values. A three-step modeling procedure was used to make our model: (a) wrap the zero-stress vessel sector to obtain the residual stress; (b) stretch the vessel axially to its length in vivo; and (c) pressurize the vessel to recover its in vivo geometry. Eight models were constructed for our investigation. Wrapping led to reduced lumen and cap stress and increased out boundary stress. The model with axial stretch, circumferential shrink, but no wrapping overestimated lumen and cap stress by 182% and 448%, respectively. The model with wrapping, circumferential shrink, but no axial stretch predicted average lumen stress and cap stress as 0.76 kPa and -15 kPa. The same model with 10% axial stretch had 42.53 kPa lumen stress and 29.0 kPa cap stress, respectively. Skipping circumferential shrinkage leads to overexpansion of the vessel and incorrect stress/strain calculations. Vessel stiffness increase (100%) leads to 75% lumen stress increase and 102% cap stress increase.

Mechanical stress is associated with right ventricular response to pulmonary valve replacement in patients with repaired tetralogy of Fallot.

OBJECTIVE: Patients with repaired tetralogy of Fallot account for a substantial proportion of cases with late-onset right ventricular failure. The current surgical approach, which includes pulmonary valve replacement/insertion, has yielded mixed results. Therefore, it may be clinically useful to identify parameters that can be used to predict right ventricular function response to pulmonary valve replacement. METHODS: Cardiac magnetic resonance data before and 6 months after pulmonary valve replacement were obtained from 16 patients with repaired tetralogy of Fallot (8 male, 8 female; median age, 42.75 years). Right ventricular ejection fraction change from pre- to postpulmonary valve replacement was used as the outcome. The patients were divided into group 1 (n = 8, better outcome) and group 2 (n = 8, worst outcome). Cardiac magnetic resonance-based patient-specific computational right ventricular/left ventricular models were constructed, and right ventricular mechanical stress and strain, wall thickness, curvature, and volumes were obtained for analysis. RESULTS: Our results indicated that right ventricular wall stress was the best single predictor for postpulmonary valve replacement outcome with an area under the receiver operating characteristic curve of 0.819. Mean values of stress, strain, wall thickness, and longitudinal curvature differed significantly between the 2 groups with right ventricular wall stress showing the largest difference. Mean right ventricular stress in group 2 was 103% higher than in group 1. CONCLUSIONS: Computational modeling and right ventricular stress may be used as tools to identify right ventricular function response to pulmonary valve replacement. Large-scale clinical studies are needed to validate these preliminary findings.

Using contracting band to improve right ventricle ejection fraction for patients with repaired tetralogy of Fallot: a modeling study using patient-specific CMR-based 2-layer anisotropic models of human right and left ventricles.

OBJECTIVE: Patients with repaired tetralogy of Fallot account for most cases of late-onset right ventricle failure. The current surgical approach, which includes pulmonary valve replacement/insertion, has yielded mixed results. A new surgical option of placing an elastic band in the right ventricle is proposed to improve right ventricular cardiac function as measured by the ejection fraction. METHODS: A total of 20 computational right ventricular/left ventricular/patch/band combination models using cardiac magnetic resonance imaging from a patient with tetralogy of Fallot were constructed to investigate the effect of band material stiffness variations, band length, and active contraction. These models included 4 different band material properties, 3 band length, 3 active contracting band materials, and models with patch and scar replaced by contracting tissue. RESULTS: Our results indicated that the band insertion, combined with active band contraction and tissue regeneration techniques that restore right ventricular myocardium, has the potential to improve right ventricular ejection fraction by 7.5% (41.63% ejection fraction from the best active band model to more than 34.10% ejection fraction from baseline passive band model) and 4.2% (41.63% from the best active band model compared with cardiac magnetic resonance imaging-measured ejection fraction of 37.45%). CONCLUSIONS: The cardiac magnetic resonance imaging-based right ventricular/left ventricular/patch/band model provides a proof of concept for using elastic bands to improve right ventricular cardiac function. Band insertion, combined with myocardium regeneration techniques and right ventricular remodeling surgical procedures, has the potential to improve ventricular function in patients with repaired tetralogy of Fallot and other similar forms of right ventricular dysfunction after surgery. Additional investigations using in vitro experiments, animal models, and, finally, patient studies are warranted.

Image-Based Patient-Specific Ventricle Models with Fluid-Structure Interaction for Cardiac Function Assessment and Surgical Design Optimization.

Recent advances in medical imaging technology and computational modeling techniques are making it possible that patient-specific computational ventricle models be constructed and used to test surgical hypotheses and replace empirical and often risky clinical experimentation to examine the efficiency and suitability of various reconstructive procedures in diseased hearts. In this paper, we provide a brief review on recent development in ventricle modeling and its potential application in surgical planning and management of tetralogy of Fallot (ToF) patients. Aspects of data acquisition, model selection and construction, tissue material properties, ventricle layer structure and tissue fiber orientations, pressure condition, model validation and virtual surgery procedures (changing patient-specific ventricle data and perform computer simulation) were reviewed. Results from a case study using patient-specific cardiac magnetic resonance (CMR) imaging and right/left ventricle and patch (RV/LV/Patch) combination model with fluid-structure interactions (FSI) were reported. The models were used to evaluate and optimize human pulmonary valve replacement/insertion (PVR) surgical procedure and patch design and test a surgical hypothesis that PVR with small patch and aggressive scar tissue trimming in PVR surgery may lead to improved recovery of RV function and reduced stress/strain conditions in the patch area.

Patient-specific MRI-based 3D FSI RV/LV/patch models for pulmonary valve replacement surgery and patch optimization.

A patient-specific right/left ventricle and patch (RV/LV/patch) combination model with fluid-structure interactions (FSIs) was introduced to evaluate and optimize human pulmonary valve replacement/insertion (PVR) surgical procedure and patch design. Cardiac magnetic resonance (CMR) imaging studies were performed to acquire ventricle geometry, flow velocity, and flow rate for healthy volunteers and patients needing RV remodeling and PVR before and after scheduled surgeries. CMR-based RV/LV/patch FSI models were constructed to perform mechanical analysis and assess RV cardiac functions. Both pre- and postoperation CMR data were used to adjust and validate the model so that predicted RV volumes reached good agreement with CMR measurements (error <3%). Two RV/LV/patch models were made based on preoperation data to evaluate and compare two PVR surgical procedures: (i) conventional patch with little or no scar tissue trimming, and (ii) small patch with aggressive scar trimming and RV volume reduction. Our modeling results indicated that (a) patient-specific CMR-based computational modeling can provide accurate assessment of RV cardiac functions, and (b) PVR with a smaller patch and more aggressive scar removal led to reduced stress/strain conditions in the patch area and may lead to improved recovery of RV functions. More patient studies are needed to validate our findings.

Two-layer passive/active anisotropic FSI models with fiber orientation: MRI-based patient-specific modeling of right ventricular response to pulmonary valve insertion surgery.

A single-layer patient specific right/left ventricle patch (RV/LV/Patch) combination model with fluid-structure interactions (FSI) was introduced in our previous papers to evaluate and optimize human pulmonary valve replacement/insertion (PVR) surgical procedure and patch design. In this paper, an active anisotropic model with two-layer structure for ventricle wall and tissue fiber orientation was introduced to improve previous isotropic model for more accurate assessment of RV function and potential application in PVR surgery and patch design. A material-stiffening approach was used to model active heart contraction. The computational models were used to conduct "virtual (computational)" surgeries and test the hypothesis that a PVR surgical design with a smaller patch and more aggressive scar tissue trimming would lead to improved RV cardiac function recovery. Results from our models validated by pre-operation data indicated that the small patch design had 11% improvement in RV function as measured by RV ejection fraction, compared to the conventional patch. Maximum Stress-P1 value from the active anisotropic model was 121.2% higher than that from the passive isotropic model. Computational RV volume predictions agreed well with CMR-measured volume data (error < 2%).

In vivo MRI-Based 3D FSI RV/LV Models for Human Right Ventricle and Patch Design for Potential Computer-Aided Surgery Optimization.

Right ventricular dysfunction is one of the more common causes of heart failure in patients with congenital heart defects. Use of computer-assisted procedures is becoming more popular in clinical decision making process and computer-aided surgeries. A 3D in vivo MRI-based RV/LV combination model with fluid-structure interaction (FSI), RV-LV interaction, and RV-patch interaction was introduced to perform mechanical analysis for human right ventricle with potential clinical applications. Patient-specific RV/LV morphologies were acquired by using planar tagged MRI. The 3D RV/LV FSI model was solved using a commercial finite element package ADINA. Our results indicated that flow and stress/strain distributions in the right ventricle are closely related to RV morphology, material properties and blood pressure conditions. Patches with material properties better matching RV tissue properties and smaller size lead to better RV function recoveries. Computational RV volumes showed very good agreement with MRI data (error < 3%). More patient studies are needed to establish baseline database so that computational simulations can be used to replace empirical and often risky clinical experimentation to examine the efficiency and suitability of various reconstructive procedures in diseased hearts and optimal design can be found.

Quantifying effects of plaque structure and material properties on stress distributions in human atherosclerotic plaques using 3D FSI models.

BACKGROUND: Atherosclerotic plaques may rupture without warning and cause acute cardiovascular syndromes such as heart attack and stroke. Methods to assess plaque vulnerability noninvasively and predict possible plaque rupture are urgently needed. METHOD: MRI-based three-dimensional unsteady models for human atherosclerotic plaques with multi-component plaque structure and fluid-structure interactions are introduced to perform mechanical analysis for human atherosclerotic plaques. RESULTS: Stress variations on critical sites such as a thin cap in the plaque can be 300% higher than that at other normal sites. Large calcification block considerably changes stress/strain distributions. Stiffness variations of plaque components (50% reduction or 100% increase) may affect maximal stress values by 20-50%. Plaque cap erosion causes almost no change on maximal stress level at the cap, but leads to 50% increase in maximal strain value. CONCLUSIONS: Effects caused by atherosclerotic plaque structure, cap thickness and erosion, material properties, and pulsating pressure conditions on stress/strain distributions in the plaque are quantified by extensive computational case studies and parameter evaluations. Computational mechanical analysis has good potential to improve accuracy of plaque vulnerability assessment.

Pattern formation of vascular smooth muscle cells subject to nonuniform fluid shear stress: role of PDGF-beta receptor and Src.

Blood vessels are subject to fluid shear stress, a hemodynamic factor that inhibits the mitogenic activities of vascular cells. The presence of nonuniform shear stress has been shown to exert graded suppression of cell proliferation and induces the formation of cell density gradients, which in turn regulate the direction of smooth muscle cell (SMC) migration and alignment. Here, we investigated the role of platelet-derived growth factor (PDGF)-beta receptor and Src in the regulation of such processes. In experimental models with vascular polymer implants, SMCs migrated from the vessel media into the neointima of the implant under defined fluid shear stress. In a nonuniform shear model, blood shear stress suppressed the expression of PDGF-beta receptor and the phosphorylation of Src in a shear level-dependent manner, resulting in the formation of mitogen gradients, which were consistent with the gradient of cell density as well as the alignment of SMCs. In contrast, uniform shear stress in a control model elicited an even influence on the activity of mitogenic molecules without modulating the uniformity of cell density and did not significantly influence the direction of SMC alignment. The suppression of the PDGF-beta receptor tyrosine kinase and Src with pharmacological substances diminished the gradients of mitogens and cell density and reduced the influence of nonuniform shear stress on SMC alignment. These observations suggest that PDGF-beta receptor and Src possibly serve as mediating factors in nonuniform shear-induced formation of cell density gradients and alignment of SMCs in the neointima of vascular polymer implants.