Validation of the Reduced Unified Continuum Formulation Against In Vitro 4D-Flow MRI.

We previously introduced and verified the reduced unified continuum formulation for vascular fluid-structure interaction (FSI) against Womersley's deformable wall theory. Our present work seeks to investigate its performance in a patient-specific aortic setting in which assumptions of idealized geometries and velocity profiles are invalid. Specifically, we leveraged 2D magnetic resonance imaging (MRI) and 4D-flow MRI to extract high-resolution anatomical and hemodynamic information from an in vitro flow circuit embedding a compliant 3D-printed aortic phantom. To accurately reflect experimental conditions, we numerically implemented viscoelastic external tissue support, vascular tissue prestressing, and skew boundary conditions enabling in-plane vascular motion at each inlet and outlet. Validation of our formulation is achieved through close quantitative agreement in pressures, lumen area changes, pulse wave velocity, and early systolic velocities, as well as qualitative agreement in late systolic flow structures. Our validated suite of FSI techniques offers a computationally efficient approach for numerical simulation of vascular hemodynamics. This study is among the first to validate a cardiovascular FSI formulation against an in vitro flow circuit involving a compliant vascular phantom of complex patient-specific anatomy.

A matched-pair case control study identifying hemodynamic predictors of cerebral aneurysm growth using computational fluid dynamics.

Introduction: Initiation and progression of cerebral aneurysms is known to be driven by complex interactions between biological and hemodynamic factors, but the hemodynamic mechanism which drives aneurysm growth is unclear. We employed robust modeling and computational methods, including temporal and spatial convergence studies, to study hemodynamic characteristics of cerebral aneurysms and identify differences in these characteristics between growing and stable aneurysms. Methods: Eleven pairs of growing and non-growing cerebral aneurysms, matched in both size and location, were modeled from MRA and CTA images, then simulated using computational fluid dynamics (CFD). Key hemodynamic characteristics, including wall shear stress (WSS), oscillatory shear index (OSI), and portion of the aneurysm under low shear, were evaluated. Statistical analysis was then performed using paired Wilcoxon rank sum tests. Results: The portion of the aneurysm dome under 70% of the parent artery mean wall shear stress was higher in growing aneurysms than in stable aneurysms and had the highest significance among the tested metrics (p = 0.08). Other metrics of area under low shear had similar levels of significance. Discussion: These results align with previously observed hemodynamic trends in cerebral aneurysms, indicating a promising direction for future study of low shear area and aneurysm growth. We also found that mesh resolution significantly affected simulated WSS in cerebral aneurysms. This establishes that robust computational modeling methods are necessary for high fidelity results. Together, this work demonstrates that complex hemodynamics are at play within cerebral aneurysms, and robust modeling and simulation methods are needed to further study this topic.

Virtual Transcatheter Interventions for Peripheral Pulmonary Artery Stenosis in Williams and Alagille Syndromes.

Background Despite favorable outcomes of surgical pulmonary artery (PA) reconstruction, isolated proximal stenting of the central PAs is common clinical practice for patients with peripheral PA stenosis in association with Williams and Alagille syndromes. Given the technical challenges of PA reconstruction and the morbidities associated with transcatheter interventions, the hemodynamic consequences of all treatment strategies must be rigorously assessed. Our study aims to model, assess, and predict hemodynamic outcomes of transcatheter interventions in these patients. Methods and Results Isolated proximal and "extensive" interventions (stenting and/or balloon angioplasty of proximal and lobar vessels) were performed in silico on 6 patient-specific PA models. Autoregulatory adaptation of the cardiac output and downstream arterial resistance was modeled in response to intervention-induced hemodynamic perturbations. Postintervention computational fluid dynamics predictions were validated in 2 stented patients and quantitatively assessed in 4 surgical patients. Our computational methods accurately predicted postinterventional PA pressures, the primary indicators of success for treatment of peripheral PA stenosis. Proximal and extensive treatment achieved median reductions of 14% and 40% in main PA systolic pressure, 27% and 56% in pulmonary vascular resistance, and 10% and 45% in right ventricular stroke work, respectively. Conclusions In patients with Williams and Alagille syndromes, extensive transcatheter intervention is required to sufficiently reduce PA pressures and right ventricular stroke work. Transcatheter therapy was shown to be ineffective for long-segment stenosis and pales hemodynamically in comparison with published outcomes of surgical reconstruction. Regardless of the chosen strategy, a virtual treatment planning platform could identify lesions most critical for optimizing right ventricular afterload.

Computational simulation-derived hemodynamic and biomechanical properties of the pulmonary arterial tree early in the course of ventricular septal defects.

Untreated ventricular septal defects (VSDs) can lead to pulmonary arterial hypertension (PAH) characterized by elevated pulmonary artery (PA) pressure and vascular remodeling, known as PAH associated with congenital heart disease (PAH-CHD). Though previous studies have investigated hemodynamic effects on vascular mechanobiology in late-stage PAH, hemodynamics leading to PAH-CHD initiation have not been fully quantified. We hypothesize that abnormal hemodynamics from left-to-right shunting in early stage VSDs affects PA biomechanical properties leading to PAH initiation. To model PA hemodynamics in healthy, small, moderate, and large VSD conditions prior to the onset of vascular remodeling, computational fluid dynamics simulations were performed using a 3D finite element model of a healthy 1-year-old's proximal PAs and a body-surface-area-scaled 0D distal PA tree. VSD conditions were modeled with increased pulmonary blood flow to represent degrees of left-to-right shunting. In the proximal PAs, pressure, flow, strain, and wall shear stress (WSS) increased with increasing VSD size; oscillatory shear index decreased with increasing VSD size in the larger PA vessels. WSS was higher in smaller diameter vessels and increased with VSD size, with the large VSD condition exhibiting WSS >100 dyn/cm[Formula: see text], well above values typically used to study dysfunctional mechanotransduction pathways in PAH. This study is the first to estimate hemodynamic and biomechanical metrics in the entire pediatric PA tree with VSD severity at the stage leading to PAH initiation and has implications for future studies assessing effects of abnormal mechanical stimuli on endothelial cells and vascular wall mechanics that occur during PAH-CHD initiation and progression.

Transcriptomic and Functional Analyses of Mitochondrial Dysfunction in Pressure Overload-Induced Right Ventricular Failure.

Background In complex congenital heart disease patients such as those with tetralogy of Fallot, the right ventricle (RV) is subject to pressure overload, leading to RV hypertrophy and eventually RV failure. The mechanisms that promote the transition from stable RV hypertrophy to RV failure are unknown. We evaluated the role of mitochondrial bioenergetics in the development of RV failure. Methods and Results We created a murine model of RV pressure overload by pulmonary artery banding and compared with sham-operated controls. Gene expression by RNA-sequencing, oxidative stress, mitochondrial respiration, dynamics, and structure were assessed in pressure overload-induced RV failure. RV failure was characterized by decreased expression of electron transport chain genes and mitochondrial antioxidant genes (aldehyde dehydrogenase 2 and superoxide dismutase 2) and increased expression of oxidant stress markers (heme oxygenase, 4-hydroxynonenal). The activities of all electron transport chain complexes decreased with RV hypertrophy and further with RV failure (oxidative phosphorylation: sham 552.3±43.07 versus RV hypertrophy 334.3±30.65 versus RV failure 165.4±36.72 pmol/(s×mL), P<0.0001). Mitochondrial fission protein DRP1 (dynamin 1-like) trended toward an increase, while MFF (mitochondrial fission factor) decreased and fusion protein OPA1 (mitochondrial dynamin like GTPase) decreased. In contrast, transcription of electron transport chain genes increased in the left ventricle of RV failure. Conclusions Pressure overload-induced RV failure is characterized by decreased transcription and activity of electron transport chain complexes and increased oxidative stress which are associated with decreased energy generation. An improved understanding of the complex processes of energy generation could aid in developing novel therapies to mitigate mitochondrial dysfunction and delay the onset of RV failure.

4HNE Impairs Myocardial Bioenergetics in Congenital Heart Disease-Induced Right Ventricular Failure.

BACKGROUND: In patients with complex congenital heart disease, such as those with tetralogy of Fallot, the right ventricle (RV) is subject to pressure overload stress, leading to RV hypertrophy and eventually RV failure. The role of lipid peroxidation, a potent form of oxidative stress, in mediating RV hypertrophy and failure in congenital heart disease is unknown. METHODS: Lipid peroxidation and mitochondrial function and structure were assessed in right ventricle (RV) myocardium collected from patients with RV hypertrophy with normal RV systolic function (RV fractional area change, 47.3±3.8%) and in patients with RV failure showing decreased RV systolic function (RV fractional area change, 26.6±3.1%). The mechanism of the effect of lipid peroxidation, mediated by 4-hydroxynonenal ([4HNE] a byproduct of lipid peroxidation) on mitochondrial function and structure was assessed in HL1 murine cardiomyocytes and human induced pluripotent stem cell-derived cardiomyocytes. RESULTS: RV failure was characterized by an increase in 4HNE adduction of metabolic and mitochondrial proteins (16 of 27 identified proteins), in particular electron transport chain proteins. Sarcomeric (myosin) and cytoskeletal proteins (desmin, tubulin) also underwent 4HNE adduction. RV failure showed lower oxidative phosphorylation (moderate RV hypertrophy, 287.6±19.75 versus RV failure, 137.8±11.57 pmol/[sec×mL]; P=0.0004), and mitochondrial structural damage. Using a cell model, we show that 4HNE decreases cell number and oxidative phosphorylation (control, 388.1±23.54 versus 4HNE, 143.7±11.64 pmol/[sec×mL]; P<0.0001). Carvedilol, a known antioxidant did not decrease 4HNE adduction of metabolic and mitochondrial proteins and did not improve oxidative phosphorylation. CONCLUSIONS: Metabolic, mitochondrial, sarcomeric, and cytoskeletal proteins are susceptible to 4HNE-adduction in patients with RV failure. 4HNE decreases mitochondrial oxygen consumption by inhibiting electron transport chain complexes. Carvedilol did not improve the 4HNE-mediated decrease in oxygen consumption. Strategies to decrease lipid peroxidation could improve mitochondrial energy generation and cardiomyocyte survival and improve RV failure in patients with congenital heart disease.

Fluid-structure interaction modeling of blood flow in the pulmonary arteries using the unified continuum and variational multiscale formulation.

In this work, we present a computational fluid-structure interaction (FSI) study for a healthy patient-specific pulmonary arterial tree using the unified continuum and variational multiscale (VMS) formulation we previously developed. The unified framework is particularly well-suited for FSI, as the fluid and solid sub-problems are addressed in essentially the same manner and can thus be uniformly integrated in time with the generalized-α method. In addition, the VMS formulation provides a mechanism for large-eddy simulation in the fluid sub-problem and pressure stabilization in the solid sub-problem. The FSI problem is solved in a quasi-direct approach, in which the pressure and velocity in the unified continuum body are first solved, and the solid displacement is then obtained via a segregated algorithm and prescribed as a boundary condition for the mesh motion. Results of the pulmonary arterial FSI simulation are presented and compared against those of a rigid wall simulation.

Changes in event-related potential functional networks predict traumatic brain injury in piglets.

BACKGROUND: Traumatic brain injury is a leading cause of cognitive and behavioral deficits in children in the US each year. None of the current diagnostic tools, such as quantitative cognitive and balance tests, have been validated to identify mild traumatic brain injury in infants, adults and animals. In this preliminary study, we report a novel, quantitative tool that has the potential to quickly and reliably diagnose traumatic brain injury and which can track the state of the brain during recovery across multiple ages and species. METHODS: Using 32 scalp electrodes, we recorded involuntary auditory event-related potentials from 22 awake four-week-old piglets one day before and one, four, and seven days after two different injury types (diffuse and focal) or sham. From these recordings, we generated event-related potential functional networks and assessed whether the patterns of the observed changes in these networks could distinguish brain-injured piglets from non-injured. FINDINGS: Piglet brains exhibited significant changes after injury, as evaluated by five network metrics. The injury prediction algorithm developed from our analysis of the changes in the event-related potentials functional networks ultimately produced a tool with 82% predictive accuracy. INTERPRETATION: This novel approach is the first application of auditory event-related potential functional networks to the prediction of traumatic brain injury. The resulting tool is a robust, objective and predictive method that offers promise for detecting mild traumatic brain injury, in particular because collecting event-related potentials data is noninvasive and inexpensive.

Infant skull fracture risk for low height falls.

Skull fractures are common injuries in young children, typically caused by accidental falls and child abuse. The paucity of detailed biomechanical data from real-world trauma in children has hampered development of biomechanical thresholds for skull fracture in infants. The objectives of this study were to identify biomechanical metrics to predict skull fracture, determine threshold values associated with fracture, and develop skull fracture risk curves for low-height falls in infants. To achieve these objectives, we utilized an integrated approach consisting of case evaluation, anthropomorphic reconstruction, and finite element simulation. Four biomechanical candidates for predicting skull fracture were identified (first principal stress, first principal strain, shear stress, and von Mises stress) and evaluated against well-witnessed falls in infants (0-6 months). Among the predictor candidates, first principal stress and strain correlated best with the occurrence of parietal skull fracture. The principal stress and strain thresholds associated with 50 and 95% probability of parietal skull fracture were 25.229 and 36.015 MPa and 0.0464 and 0.0699, respectively. Risk curves using these predictors determined that infant falls from 0.3 m had a low probability (0-54%) to result in parietal skull fracture, particularly with carpet impact (0-1%). Head-first falls from 0.9 m had a high probability of fracture (86-100%) for concrete impact and a moderate probability (34-81%) for carpet impact. Probabilities of fracture in 0.6 m falls were dependent on impact surface. Occipital impacts from 0.9 m onto the concrete also had the potential (27-90% probability) to generate parietal skull fracture. These data represent a multi-faceted biomechanical assessment of infant skull fracture risk and can assist in the differential diagnosis for head trauma in children.

EGF augments TGFß-induced epithelial-mesenchymal transition by promoting SHP2 binding to GAB1.

In many epithelial cells, epidermal growth factor (EGF) augments the epithelial-mesenchymal transition (EMT) that occurs when cells are treated with transforming growth factor ß (TGFß). We demonstrate that this augmentation requires activation of SH2 domain-containing phosphatase-2 (SHP2; also known as PTPN11), a proto-oncogene. In lung and pancreatic cancer cell lines, reductions in E-cadherin expression, increases in vimentin expression and increases in cell scatter rates were larger when cells were treated with TGFß and EGF versus TGFß or EGF alone. SHP2 knockdown promoted epithelial characteristics basally and antagonized EMT in response to TGFß alone or in combination with EGF. Whereas EGF promoted SHP2 binding to tyrosine phosphorylated GAB1, which promotes SHP2 activity, TGFß did not induce SHP2 association with phosphotyrosine-containing proteins. Knockdown of endogenous SHP2 and reconstitution with an SHP2 mutant with impaired phosphotyrosine binding ability eliminated the EGF-mediated EMT augmentation that was otherwise restored with wild-type SHP2 reconstitution. These results demonstrate roles for basal and ligand-induced SHP2 activity in EMT and further motivate efforts to identify specific ways to inhibit SHP2, given the role of EMT in tumor dissemination and chemoresistance.