Guidelines for mechanistic modeling and analysis in cardiovascular research.

Computational, or in silico, models are an effective, noninvasive tool for investigating cardiovascular function. These models can be used in the analysis of experimental and clinical data to identify possible mechanisms of (ab)normal cardiovascular physiology. Recent advances in computing power and data management have led to innovative and complex modeling frameworks that simulate cardiovascular function across multiple scales. While commonly used in multiple disciplines, there is a lack of concise guidelines for the implementation of computer models in cardiovascular research. In line with recent calls for more reproducible research, it is imperative that scientists adhere to credible practices when developing and applying computational models to their research. The goal of this manuscript is to provide a consensus document that identifies best practices for in silico computational modeling in cardiovascular research. These guidelines provide the necessary methods for mechanistic model development, model analysis, and formal model calibration using fundamentals from statistics. We outline rigorous practices for computational, mechanistic modeling in cardiovascular research and discuss its synergistic value to experimental and clinical data.

Application and reduction of a nonlinear hyperelastic wall model capturing ex vivo relationships between fluid pressure, area, and wall thickness in normal and hypertensive murine left pulmonary arteries.

Pulmonary hypertension is a cardiovascular disorder manifested by elevated mean arterial blood pressure (>20 mmHg) together with vessel wall stiffening and thickening due to alterations in collagen, elastin, and smooth muscle cells. Hypoxia-induced (type 3) pulmonary hypertension can be studied in animals exposed to a low oxygen environment for prolonged time periods leading to biomechanical alterations in vessel wall structure. This study introduces a novel approach to formulating a reduced order nonlinear elastic structural wall model for a large pulmonary artery. The model relating blood pressure and area is calibrated using ex vivo measurements of vessel diameter and wall thickness changes, under controlled pressure conditions, in left pulmonary arteries isolated from control and hypertensive mice. A two-layer, hyperelastic, and anisotropic model incorporating residual stresses is formulated using the Holzapfel-Gasser-Ogden model. Complex relations predicting vessel area and wall thickness with increasing blood pressure are derived and calibrated using the data. Sensitivity analysis, parameter estimation, subset selection, and physical plausibility arguments are used to systematically reduce the 16-parameter model to one in which a much smaller subset of identifiable parameters is estimated via solution of an inverse problem. Our final reduced one layer model includes a single set of three elastic moduli. Estimated ranges of these parameters demonstrate that nonlinear stiffening is dominated by elastin in the control animals and by collagen in the hypertensive animals. The pressure-area relation developed in this novel manner has potential impact on one-dimensional fluids network models of vessel wall remodeling in the presence of cardiovascular disease.

Estimating pulmonary arterial remodeling via an animal-specific computational model of pulmonary artery stenosis.

Pulmonary artery stenosis (PAS) often presents in children with congenital heart disease, altering blood flow and pressure during critical periods of growth and development. Variability in stenosis onset, duration, and severity result in variable growth and remodeling of the pulmonary vasculature. Computational fluid dynamics (CFD) models enable investigation into the hemodynamic impact and altered mechanics associated with PAS. In this study, a one-dimensional (1D) fluid dynamics model was used to simulate hemodynamics throughout the pulmonary arteries of individual animals. The geometry of the large pulmonary arteries was prescribed by animal-specific imaging, whereas the distal vasculature was simulated by a three-element Windkessel model at each terminal vessel outlet. Remodeling of the pulmonary vasculature, which cannot be measured in vivo, was estimated via model-fitted parameters. The large artery stiffness was significantly higher on the left side of the vasculature in the left pulmonary artery (LPA) stenosis group, but neither side differed from the sham group. The sham group exhibited a balanced distribution of total distal vascular resistance, whereas the left side was generally larger in the LPA stenosis group, with no significant differences between groups. In contrast, the peripheral compliance on the right side of the LPA stenosis group was significantly greater than the corresponding side of the sham group. Further analysis indicated the underperfused distal vasculature likely moderately decreased in radius with little change in stiffness given the increase in thickness observed with histology. Ultimately, our model enables greater understanding of pulmonary arterial adaptation due to LPA stenosis and has potential for use as a tool to noninvasively estimate remodeling of the pulmonary vasculature.

Right ventricular outflow tract diameter change with exercise: a prospective exercise echocardiography and invasive CPET study.

While cardiac output reserve with exercise predicts outcomes in cardiac and pulmonary vascular disease, precise quantification of exercise cardiac output requires invasive cardiopulmonary testing (iCPET). To improve the accuracy of cardiac output reserve estimation with transthoracic echocardiography (TTE), this prospective study aims to define changes in right ventricular outflow tract diameter (RVOTd) with exercise and its relationship with invasively measured haemodynamics. Twenty subjects underwent simultaneous TTE and iCPET, with data collected at rest, leg-raise, 25 W, 50 W (n = 16), 75 W (n = 14), and 100 W (n = 6). This was followed by a second exercise study with real-time RV pressure-volume loops at similar stages (except leg-raise). The overall cohort included heart failure with preserved ejection fraction (n = 12), pulmonary arterial hypertension (n = 5), and non-cardiac dyspnoea (n = 3). RVOTd was reverse engineered from the TTE-derived RVOT velocity time integral (VTI) and iCPET-derived stroke volume, using the formula: Fick stroke volume = RVOT VTI × RVOT area (wherein RVOT area = π × [RVOTd/2]2). RVOTd increased by nearly 3-4% at every 25 W increment. Using linear regression models, where each subject is treated as a categorical variable and adjusting for subject intercept, RVOTd was correlated with haemodynamic variables (cardiac output, heart rate, pulmonary artery and RV pressures). Of all the predictor haemodynamic variables, cardiac output had the highest r2 model fit (adjusted r2 = 0.68), with a unit increase in cardiac output associated with a 0.0678 increase in RVOTd (P < 0.001). Our findings indicate that RVOTd increases by 3-4% with every 25 W increment, predominantly correlated with cardiac output augmentation. These results can improve the accuracy of cardiac output reserve estimation by adjusting for RVOTd with graded exercise during non-invasive CPET and echocardiogram. However, future studies are needed to define these relationships for left ventricular outflow tract diameter.

Rescuing lung development through embryonic inhibition of histone acetylation.

A major barrier to the impact of genomic diagnosis in patients with congenital malformations is the lack of understanding regarding how sequence variants contribute to disease pathogenesis and whether this information could be used to generate patient-specific therapies. Congenital diaphragmatic hernia (CDH) is among the most common and severe of all structural malformations; however, its underlying mechanisms are unclear. We identified loss-of-function sequence variants in the epigenomic regulator gene SIN3A in two patients with complex CDH. Tissue-specific deletion of Sin3a in mice resulted in defects in diaphragm development, lung hypoplasia, and pulmonary hypertension, the cardinal features of CDH and major causes of CDH-associated mortality. Loss of SIN3A in the lung mesenchyme resulted in reduced cellular differentiation, impaired cell proliferation, and increased DNA damage. Treatment of embryonic Sin3a mutant mice with anacardic acid, an inhibitor of histone acetyltransferase, reduced DNA damage, increased cell proliferation and differentiation, improved lung and pulmonary vascular development, and reduced pulmonary hypertension. These findings demonstrate that restoring the balance of histone acetylation can improve lung development in the Sin3a mouse model of CDH.