Next Article in Journal
Immunohistochemical Analysis of Single-Stranded DNA Binding Protein 2 in Non-Melanoma Skin Cancers
Previous Article in Journal
Altered Auditory and Visual Evoked Potentials following Single and Repeated Low-Velocity Head Rotations in 4-Week-Old Swine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Aortic Remodeling Kinetics in Response to Coarctation-Induced Mechanical Perturbations

1
Joint Department of Biomedical Engineering, Medical College of Wisconsin, Marquette University, Milwaukee, WI 53226, USA
2
Section of Pediatric Cardiology, Department of Pediatrics, Medical College of Wisconsin, Milwaukee, WI 53226, USA
3
Herma Heart Institute, Children’s Wisconsin, Milwaukee, WI 53226, USA
4
Department of Mathematical and Statistical Sciences, Marquette University, Milwaukee, WI 53233, USA
5
Division of Pediatric Cardiothoracic Surgery, Department of Surgery, Medical College of Wisconsin, Herma Heart Institute, Children’s Wisconsin, Milwaukee, WI 53226, USA
6
PoliToBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, 10129 Turin, Italy
7
Department of Radiology, Medical College of Wisconsin, Milwaukee, WI 53226, USA
8
Department of Physiology, Medical College of Wisconsin, Milwaukee, WI 53226, USA
*
Author to whom correspondence should be addressed.
Biomedicines 2023, 11(7), 1817; https://doi.org/10.3390/biomedicines11071817
Submission received: 24 May 2023 / Revised: 16 June 2023 / Accepted: 21 June 2023 / Published: 25 June 2023
(This article belongs to the Section Biomedical Engineering and Materials)

Abstract

:
Background: Coarctation of the aorta (CoA; constriction of the proximal descending thoracic aorta) is among the most common congenital cardiovascular defects. Coarctation-induced mechanical perturbations trigger a cycle of mechano-transduction events leading to irreversible precursors of hypertension including arterial thickening, stiffening, and vasoactive dysfunction in proximal conduit arteries. This study sought to identify kinetics of the stress-mediated compensatory response leading to these alterations using a preclinical rabbit model of CoA. Methods: A prior growth and remodeling (G&R) framework was reformulated and fit to empirical measurements from CoA rabbits classified into one control and nine CoA groups of various severities and durations (n = 63, 5–11/group). Empirical measurements included Doppler ultrasound imaging, uniaxial extension testing, catheter-based blood pressure, and wire myography, yielding the time evolution of arterial thickening, stiffening, and vasoactive dysfunction required to fit G&R constitutive parameters. Results: Excellent agreement was observed between model predictions and observed patterns of arterial thickening, stiffening, and dysfunction among all CoA groups. For example, predicted vascular impairment was not significantly different from empirical observations via wire myography (p-value > 0.13). Specifically, 48% and 45% impairment was observed in smooth muscle contraction and endothelial-dependent relaxation, respectively, which were accurately predicted using the G&R model. Conclusions: The resulting G&R model, for the first time, allows for prediction of hypertension precursors at neonatal ages that is currently challenging to examine in preclinical models. These findings provide a validated computational tool for prediction of persistent arterial dysfunction and identification of revised severity–duration thresholds that may ultimately avoid hypertension from CoA.

1. Introduction

Coarctation of the aorta (CoA) is a constriction of the proximal descending thoracic aorta and is one of the most common congenital cardiovascular defects [1]. Surgical repair remains the standard approach for treatment in most young infants, while catheter intervention with balloon angioplasty and/or stenting is the standard of care in older children and adults [2,3]. Guidelines for intervention include a transcatheter peak-to-peak blood pressure (BP) gradient ≥20 mmHg, and many published reports regard this criterion as a hemodynamically significant CoA [3]. Unfortunately, hypertension and left-ventricular hypertrophy are common despite successful CoA repair based upon this threshold [4,5].
Mechanistic understanding of hypertension in CoA is hindered by the paucity of data to investigate changes in the kinetics of growth and remodeling (G&R) associated with coarctation-induced mechanical stimuli [6,7]. Hence, animal models are used to investigate underling mechanisms and precursors of hypertension in response to CoA [8]. In contrast to the current putative threshold, recent studies using a rabbit model of CoA suggest that precursors of hypertension are present among treated rabbits with CoA gradients <20 mmHg [8,9,10]. This finding raises doubts about the efficacy of the current guidelines in preventing refractory hypertension after treatment. Interestingly, the duration of coarctation-induced mechanical stimuli also correlates with hypertension precursors [7,9], further supporting the importance of earlier treatment that has been shown to reduce hypertension prevalence among treated CoA patients [11].
In general, the arterial wall is assumed to adapt to mechanical stimuli (perturbation from normal state) and recover a homeostatic state via smooth muscle (SM) cell synthesis and matrix turnover [12,13,14]. For example, in CoA, the presence of the narrowing can cause elevated BP proximally and a turbulent velocity jet distally. This perturbation in blood pressure-induced intramural wall stress (IWS: σ ) proximally and blood flow-induced wall shear stress (WSS: ϕ ) distally changes the vascular tone toward recovering IWS and WSS back to preferred ranges that are assumed to be ‘known’ to the local vasculature. Consequently, fibrillar proteins such as collagen and SM are deposited to remodel the extracellular matrix and modify thickness, arterial caliber, and material properties. The kinematics of this compensatory thickening/stiffening response is stress-driven and has been previously quantified as a linear function of deviations from homeostatic values [14,15,16].
The aim of this work is to quantify G&R parameters applicable to CoA using a computational approach tuned to replicate arterial thickening, stiffening, and dysfunction measured in rabbit aortas exposed to mechanical stimuli within ranges seen clinically. Completion of the study objective will provide a validated computational tool for prediction of persistent arterial dysfunction and allow for the much-needed identification of severity–duration thresholds that prevent arterial dysfunctions contributing to refractory hypertension among treated CoA patients.

2. Material and Methods

2.1. Stress-Mediated G&R Formulation

For a section of the aorta, at any time t , the current configuration can be represented simply as a pressurized cylindrical thin-wall tube. Deformation can, therefore, be described by mapping the reference point ( R , Θ , Z ) onto a new position ( r , θ , z ) in the same coordinate system when deformed:
r = r R ,     θ = Θ ,     z = κ Z ,
for some function r R and parameter κ . That is, deformations are restricted so that the principal axes of deformation coincide with the coordinate directions of the principal stretches.
λ r = r R ,     λ θ = r R ,     λ z = κ .
Assuming incompressibility, i.e., λ θ λ z λ r = 1 , one can reduce the number of independent stretch parameters to λ θ = λ and λ z
λ r = λ λ z 1 = 1 κ λ .
Equilibrium in the radial direction gives the following [17]:
σ r t r + σ r t σ θ t r t = 0 ,
with σ representing the IWS in the corresponding directions. Thickness integration with the σ r   r i = p as a boundary condition on the luminal surface ( r = r i ), gives
p t = r i t r o t σ θ t σ r t d r r t .
This equation can be rewritten in reference configuration by observing that r = λ R and d r = λ κ 1 d R from Equations (2) and (3), with the constitutive relation [17,18] σ θ σ r = λ   W / λ :
p t = R i t R o t 1 κ λ   λ k W k t + W a t d R R t .
Neglecting through thickness variation of the integrand with the thin-wall assumption, i.e., h t < < r i t , yields
p t = 1 κ λ ln R o R i λ W t + W a t ,
where W t is the total strain energy per unit area at time t representing passive forces in the vasculature. On the other hand, W a t is the SM work density representing active forces acting on the arterial wall.
Strain energy per unit area can be described by modeling the artery’s ability to adapt to changes of mechanical environment to maintain a homeostatic state via vascular G&R [16]. This process occurs through removal of old vascular constituents and incorporation of new constituents, which can be described in the following time evolution form:
M k t = M k 0 Q k t + 0 t m k τ q k t τ d τ ,
where M k t describes the current mass per unit area for the vascular constituent k . The first term on the right side of the equation is the natural turnover of the initial mass. Q k t is the current remaining fraction for constituent k produced at time 0. The second term describes the natural turnover of the newly produced constituents, with q k t τ being the remaining fraction of the constituent k produced at time τ and mass production rate m k τ . Here, it is assumed that Q k t and q k t follow an exponential decay, i.e., exp ν k t with rate constant ν k capturing turnover of constituent k , i.e., 1 / 80 days for collagen and SM [19]. Elastin is assumed to neither degrade nor grow [20], i.e., Q e t = 1 and q e t = 0 .
For a constrained mixture, based on the mass-average principle, the total strain energy per unit area at time t is
W t = k M k 0 Q k t W k t , 0 + 0 t m k τ q k t τ W k t , τ d τ ,
where W k t , 0 and W k t , τ denote the strain energy density that constituent k   c , m , e produces at times 0 and τ , respectively. The vessel is assumed a constrained mixture of collagen, SM, and elastin, i.e., c , m , e .
Active stress in the circumferential direction σ a t can be defined from a constitutive relation assuming that vascular SM is oriented nearly circumferentially [21] with zero active stress in the radial and longitudinal directions.
σ a t = λ W a t λ ,   σ z a t = σ r a t = 0 .
The active contribution of circumferential stress is modulated as a linear function of generated stress per unit relaxed SM area, i.e., σ a t = λ m S . Here, G&R takes place on a relatively long timescale where the SM response is more importantly characterized by the phenotypic dedifferentiation between a contractile and a synthetic state [9,22,23] as compared to short-term calcium-driven cyclic interactions. This phenotypic shift has been reported as increased non-muscle myosin expression with decreased SM myosin [9] in response to coarctation-induced mechanical stimuli, including those derived from increased systolic, mean, and pulse BP in regions proximal to the coarctation. Although these arterial alterations are in the direction of increased medial thickness through sustained remodeling, the active contractile force is somewhat impaired possibly due to a phenotypic shift from the contractile to synthetic state [9]. Moreover, arterial thickening is at the expense of increased structural stiffness leading to augmented forward and backward pressure waves that further increase pulse pressure and BP [24,25]. This creates a positive feedback loop that may become irreversible leading to refractory hypertension even after removal of the CoA.
We assume herein that phenotypic dedifferentiation is a continuous process in which the SM cell gradually gains or loses its contractile response according to the function S = S   C , where C 0 is the phenotypic modulation stimulus. It represents the accumulation of stimuli descending from alterations in the stretch field. The generated active stress S is, therefore, associated with the slowly varying tone under normal physiological conditions [26]. Since not all infinitesimal areas of SM experience the same stretch, the expression for active stress is adapted to a situation with variable λ m , so that the total active stress difference σ a t σ r a t is computed as a sum of active stress difference contributions from remaining infinitesimal area fractions produced at different times. To derive an equation for active stress differences, the derivative of Equation (9) is taken with respect to λ for SM constituent only. Multiplying both sides by λ then gives
M m t λ W m λ = M m 0 Q m t λ t , 0 d W m t , 0 + 0 t m m τ q m t τ d W m t , τ d τ ,
where the notation d f x = d f / d x . Substituting Equation (10) gives
M m t λ W a t λ = M m 0 Q m t λ m t , 0 S C + 0 t m m τ q m t τ λ m t , τ S C d τ ,
where λ m t , 0 and λ m t , τ are the stretch ratio for SM constituents generated at times 0 and τ , respectively, and can be described in terms of stretch developed during growth, i.e., λ t / λ τ , and pre-stretch of SM G ^ m that develops at the time of production, which is assumed to be a constant equal to 1.3 according to Beak et al. [27,28]:
λ m t , τ = λ t λ τ G ^ m .
The phenotypic modulation of contractile activity is taken to be
S C = S M 1 exp C 2 ,
where S M is the maximum stress capacity of SM in its fully contractile state and assumed to be a constant equal to 100 kPa [29]. For a target lumen radius r ^ i , the stretch is λ ^ = r ^ i / R i and we refer to λ ^ as the target homeostatic stretch. From the Hagen–Poiseuille law, the λ λ ^ = r i   r ^ i =   ϕ ^ ϕ 1 / 3 is available to G&R through WSS signaling from endothelial cells. Hence, the time window average phenotypic modulation stimulus, C , of the vascular SM depends only on this known ratio λ / λ ^ , giving
M m t λ W a t λ = M m 0 Q m t λ λ 0 G ^ m S C λ λ ^ + 0 t m m τ q m t τ λ λ τ G ^ m S C λ λ ^ d τ .
Dividing both sides by λ gives
W a t λ = G ^ m λ m S C λ λ ^ ,
where the effect of the growth history is captured by a SM area-weighted harmonic mean of the stretch:
λ m = M m t × M m 0 λ 0 Q m t + 0 t m m τ λ τ q t τ d τ 1 .
The SM work density can then be formally defined as
W a λ = G ^ m λ m λ ^ λ S C λ λ ^ d λ ,
where the lower limit of the integration is arbitrary owing to the constant term of the primitive function. In a homeostatic state, taking the limit t for Equation (17) gives λ m λ ^ , so that the SM work density in the homeostatic state becomes
W ^ a λ = G ^ m λ ^ λ ^ λ S C λ λ ^ d λ .
The derivatives of the SM work density (i.e., W ^ a ), according to the constitutive relation, give the active force in pathological conditions, which can then be normalized to its homeostatic value obtained from a normal hemodynamics. Here, it was assumed that impaired active response is linearly correlated to the phenotypic shift of SM from the contractile to synthesis state. Hence, impairment was quantified in the G&R model through the ratio of CoA to homeostatic active stress, named the impairment index, and compared to corresponding values obtained empirically using wire myography, i.e., percentage impaired contractile or relaxation behavior of CoA rabbits relative to the control group.
It is worth mentioning that, since the hemodynamic modulation stimulus ( C λ λ ^ ) is an unknown function, it was represented using a truncated Taylor expansion around the homeostatic condition, λ λ ^ = 1 .
C λ λ ^ = C ^ + C λ λ ^ 1 ,
where C ^ and C are nondimensional constants representing the phenotypic modulation stimulus of SM under homeostatic and pathological conditions, respectively. Under homeostasis, S ^ = S(C[1]) = S( C ^ ) , which, according to Equation (14), gives C ^ = ln 1 S ^ / S M = 0.8326 for our choice of parameters. Therefore, we assumed C ^ as constant and tuned C in the G&R model to replicate phenotypic modulation and associated active vascular impairment obtained empirically using an impairment index defined in myography results as detailed below.

2.2. Strain Energy Density Function for Vessel Constituents

Currently, we define a general form of a strain energy density function to describe the passive action of vessel constituents. In general, if it is assumed only one family of collagen and SM is circumferentially oriented, the strain energy density function for both can be defined similar to elastin using a hybrid neo-Hookean and Fung-type material [30]:
W k λ , κ = c 1 k 2 λ 2 + κ λ 2 + κ 2 3 + c 2 k 4 c 3 k exp c 3 k λ 2 1 2 1 ,
where c i k ,   i 1 , 2 , 3 are constitutive parameters characterizing material properties of constituent k c , m , e . Now, we have all strain energy and SM work densities in terms of λ , which is a function of current and generation time, i.e., λ t , τ . Hence, the governing equation of G&R can be derived and discretized for the numerical simulation described in the next section.

2.3. Simulation of Stress Mediated G&R

For a pressurized thin-walled incompressible elastic tube with radius r and thickness h , Hoop stress can be derived by force equilibrium in the radial direction:
σ θ = p r i h ,
where h r i and denote wall thickness and internal radius, respectively.
For laminar flow of a Newtonian fluid with viscosity, η , and velocity, u , the Hagen–Poiseuille equation describes WSS as
ϕ w = 4 η u π r i 3 .
Considering that, for a thin-walled tube, area can be described as A 2 π r i h , rewriting Hoop stress and Hagen–Poiseuille in terms of A and dividing by the hemostatic state denoted by hat sign (   ^ ), one can derive the current to homeostatic area ratio in terms of mechanical stresses acting on the wall:
A ^ A = σ θ σ ^ θ   ϕ ϕ ^   2 / 3 .
This equation demonstrates how two known local stimuli, namely, the WSS and IWS, mediate to reproduce the ratio A ^ / A and, thus, the appropriate direction of growth: increased rate of production for A ^ / A > 1 and decreased rate of production for A ^ / A < 1 . Taylor series expansion of Equation (24) around the homeostatic state A ^ / A = 1 gives
A ^ A 1 + Δ σ θ σ ^ θ + 2 3 Δ ϕ ϕ ^   .
Stress-mediated remodeling theory identifies K σ = σ ^ θ 1 and K ϕ = 2 3 ϕ ^ 1 , and the basal production rate analogous to the definition of mass production rate as a linear function of mechanical stimulus is as follows [31,32]:
m k t = m 0 k 1 + K σ k Δ σ θ k + K ϕ k Δ ϕ ,
with m 0 k being the basal mass production rate of constituent k c , m , e .
In summary, Equation (26) governs the evolution of the vessel constituents through Equation (1) where the mass per unit area for each vessel constituent is defined. WSS can be derived either from the Hagen–Poiseuille law or from computational fluid dynamics (CFD) using a reduced-order 1D model (ROM). The intramural stress is quantified from constitutive relations and strain energy, and SM work densities are defined according to the growth history in the previous section, 0 τ t . Since the right-hand side of Equation (26) depends on the growth history, Equation (1) is the integral form of a delay differential equation (DDE) that is solved using an explicit numerical scheme described below.
The distribution of remaining mass per unit area is defined as
M k t , τ = m k τ q k t τ ,   t 0 ,   0 τ t ,
such that M k t , τ d τ represents the amount of constituent k that remains at time t , and produced within the time interval τ , τ + d τ . It follows that
M k t + Δ t , τ = q k t + Δ t τ q k t τ ,
for any constant time-step Δ t , and that
M k t , t = m k t q k 0   M k t , t = m k t .
Time is discretized using time-step Δ t , with the following discretization points:
t n = n Δ t ,       n = 1 ,   2 , , N .  
Then, M k becomes discretized as
M n , j k = M k n Δ t , j Δ t ,               j = 0 ,   , n .
From Equation (28), we also have
M n + 1 , j k = b n , j k M n , j k     ,       b n , j k = q k n j + 1 Δ t q k n j Δ t .
An explicit numerical scheme is used to integrate the evolution of the vessel constituents using Equation (1) with Equation (26). Slowly varying quantities, such as m k t and λ m , are treated explicitly, while λ , which changes instantaneously due to changes in the controlling parameters λ ^ t and p t , is treated implicitly. Therefore, integrals are approximated using a trapezoidal rule throughout. For a compact notation, we define a summation symbol used for the trapezoidal rule:
i = m n * a i = 0 m n 1 2 a m + i = m + 1 n 1 a i + 1 2 a n o t h e r w i s e .
The equilibrium Equations (7) and (16) give
k λ W k   n + G ^ n m λ n m S C λ n + 1 λ ^ n + 1 + ln R i R o κ λ n + 1 p n + 1 = 0 .
This equation is solved numerically for λ n + 1 with λ ^ n + 1 as an initial guess using the MATLAB numerical solver (MathWorks, Natick, MA, USA). If the solution is not unique, the closest solution to λ ^ n + 1 is chosen [26].
The stress measure σ n + 1 is computed using stress difference terms quantified through differentiating strain energy or SM work density with respect to λ and then multiplying by λ :
σ n + 1 = k λ n + 1 λ W k   n + λ n + 1 G ^ n m λ n m S C λ n + 1 λ ^ n + 1 .
For each constituent k c , m , e , the newly produced materials are computed according to Equations (26) and (29):
M n + 1 , n + 1 k = m 0 k 1 + K σ k σ n + 1 σ ^ n + 1 + K ϕ k ϕ n + 1 ϕ ^ n + 1 .
Moreover, the remaining materials M n + 1 , j k ,   j = 0 ,   1 , , n from production at earlier times is computed using Equation (32). The mass per unit area of each constituent is calculated from Equation (1) using the trapezoidal rule:
M n + 1 k = M k 0 Q k n Δ t + Δ t + Δ t j = 0 n + 1 * M n + 1 , j k .
The weighted harmonic mean of the stretch is also expressed through Equation (17) and the trapezoidal rule:
λ n + 1 m = M n + 1 m M m 0 λ 0 Q m n + 1 Δ t + Δ t j = 0 n + 1 * M n + 1 , j k 1 λ j 1 .
Using the trapezoidal rule and Equation (9) for one constituent gives
λ M k W k   n + 1 = M 0 k Q k n + 1 Δ t G ^ k λ 0 d W k λ λ 0 G ^ k + Δ t j = 0 n + 1 * M n + 1 , j k   G ^ k λ j d W k λ λ j G ^ k ,
which is a function of λ to be used for implicit calculation of λ n + 2 in the subsequent time-step. With these operations, all variables are obtained for the timestep n + 1 Δ t , and the new functions λ M k W k are defined for use in the following timestep. Here, G ^ k represents the pre-stretch for constituent k and assumed to be constant and equal to 1.08 and 1.4 for collagen and elastin, respectively [28].

2.4. G&R Model Fitting

Constitutive parameters of the derived G&R model were tuned to simulate the kinetics of growth and remodeling in response to CoA-induced mechanical stimuli. More specifically, results of the simulations for arterial thickening, stiffening, and dysfunction were compared to measurements of these precursors of hypertension quantified among nine CoA groups (details of empirical measurements are provided in the next section). Constitutive parameters were tuned to minimize the following weighted objective function:
χ 2 = k ω k S E E k ,
where S E E k is the sum of square error of the model vs. empirically measured hypertension precursors with k   ϵ   s t i f f n e s s ,   t h i c k n e s s ,   d y s f u c n t i o n , and w k represents the weight constant for the k-th hypertension precursor. A genetic algorithm was used in MATLAB [33] to solve the nonlinear optimization problem by setting the population size of 1000 and function tolerance of 1 × 10−7. Collectively, the resulting fit represented a mathematical model tuned to predict thickness evolution (from longitudinal Doppler ultrasound images) and terminal stiffness and dysfunction (from invasive final material testing and myograph analysis). Table 1 lists the G&R parameters in the model with references and/or values, where applicable.

2.5. Empirical Quantification of Hypertension Precursors

Under pathological conditions, G&R can be described through SM synthesis and collagen turnover [16] that are both assumed to be stress-driven [31] and can, therefore, be correlated to the CoA-induced mechanical stimuli. Identifying associated G&R constitutive parameters, thus, requires determination of (1) mechanical stimuli, and (2) arterial dysfunction for an array of cases with different severity and durations of CoA. For this purpose, a rabbit model of CoA accommodating various severities and durations of the mechanical stimuli was implemented [8]. Hypertension precursors including arterial thickening, stiffening, and dysfunction were quantified using an array of noninvasive and invasive experimental protocols [8].

2.5.1. Experimental Protocol

Hypertension precursors included (1) arterial thickening, (2) stiffening, and (3) dysfunction, which were quantified following previously published protocols for a rabbit model mimicking human CoA. Figure 1 summarizes the workflow for the experimental protocol and quantities measured using each modality. After IACUC approval, CoA was surgically induced under isoflurane anesthesia in New Zealand white rabbits (total 63; n = 5–11/group) by tying suture to varying diameters around the aorta at an age of ~10 weeks [8]. Tying the suture at different diameters, i.e., 2.06, 1.63, and 1.40 mm, resulted in various clinically important severities named mild (peak-to-peak CoA gradient ≤13 mmHg), intermediate (peak-to-peak CoA gradient 13–20 mmHg), and severe (peak-to-peak CoA gradient ≥20 mmHg). Importantly, the putative threshold suggestive of treatment for CoA in humans is ≥20 mmHg [37].
To investigate the effect of duration of the mechanical stimuli caused by the CoA, sutures with different dissolving properties, i.e., rapid dissolvable, dissolvable, and permanent, were used to present stenoses for different durations of ~2, 5, and 22 weeks named short, long, and prolonged CoA, respectively. Importantly, sutures were tied around the proximal descending thoracic aorta (distal to the left subclavian branch) where CoA most often presents. Together, this combination resulted in nine study groups (three severities and three durations) plus a control group.

2.5.2. Temporal Monitoring of Morphology and Hemodynamics

Using body weight at arrival and according to the growth rate for New Zealand White rabbits [38], age was estimated and then scaled to that of humans according to the duration of different phases of life [39] (Figure 1). Temporal hemodynamic and morphological evaluation was performed as a function of age using Doppler ultrasound imaging following protocols similar to human transthoracic echocardiography. Specifically, Doppler color flow and spectral Doppler images (Figure 2A) were used to identify pre- and post-stenotic flow and associated peak velocity for Doppler-based hemodynamic assessment, while b-mode imaging (Figure 2B) was used to quantify morphological properties, such as aorta diameter, wall thickness, and percentage area obstruction, using the protocols detailed elsewhere [8,9]. Moreover, phase-contrast magnetic resonance imaging (PC-MRI) was performed to characterize morphological properties, as well as quantify cardiac output and flow distribution to the main aortic branches [8,9]. Lastly, high-fidelity BP measurements were performed at the end of the protocol via invasive catheterization before harvesting of aortic tissue for further invasive characterization, e.g., material properties through uniaxial extension testing and vascular dysfunction through myograph analysis (Figure 1). Together, these empirical measurements were quantified as surrogates (i.e., precursors) of hypertension including arterial thickening, stiffening, and dysfunction, as detailed in the next subsections.

2.5.3. Vascular Stiffening

To investigate coarctation-induced arterial stiffening, material properties were quantified in the proximal descending thoracic aorta (region between the left subclavian artery and the coarctation). This region was selected due to its exposure to adversely elevated coarctation-induced BP. Material properties were characterized through uniaxial extension testing (MTS Criterion Load Frame, MTS, Minneapolis, MN, USA) at 37 °C in an environmental chamber (MTS Bionix EnviroBath, Minneapolis, MN, USA). Tissues were dissected in the circumferential direction with a length-to-width ratio of ~2.6 and preconditioned by stretching to 10% of the gauge length. Extension testing was performed at 10 mm/min until hyperelastic behavior appeared. Results were interpreted in stress–strain curves to be compared to G&R model predictions quantified via Equation (34) at the experimental endpoint.

2.5.4. Vascular Thickening

Doppler b-mode images were used to quantify aortic wall thickness in the proximal descending thoracic aorta (Figure 2B). All thickness measurements were made in triplicate with the mean value reported for each date. Body weight was also measured weekly using a scientific scale. Lastly, the thickness evolution was normalized to interpolated body weights at each timepoint assuming scaling of the body weight according to a sigmoid function (Figure 3). The sigmoid function was selected specifically to allow for a relatively low initial growth rate during the gestation period (~32 days) assuming an average of 50 g for newborn rabbit weight [39], as well as stable body weight, upon reaching adulthood.

2.5.5. Vascular Dysfunction

Vascular specimens (3–4 mm rings) were carefully sectioned from the proximal descending thoracic aorta and assessed for vascular dysfunction using a protocol described previously [8]. Briefly, experimental assessment of active vascular dysfunction was conducted to observe SM contraction via phenylephrine (PE) in a half-log increasing dose responses from 10−9 to 10−5. To assess endothelium-dependent relaxation via nitric oxide (NO), arteries were precontracted with PE to the EC50 concentration. Once plateau was achieved, cumulative addition of the endothelium-dependent muscarinic receptor agonist acetylcholine (ACh; 10−9 to 10−5 M) was initiated, and relaxation response curves were quantified as a percentage of precontracted active tension. Area under the dose response curves (AUC) were then quantified as an aggregate measure of vasoactive dysfunction and normalized to that of the control group in a normalized quantity we refer to as the impairment index. All quantifications were performed in duplicate for paired channels with mean values reported for each rabbit. The impairment index was used as a quantitative nondimensional measure of active vascular impairment due to shift in SM phenotype to the dedifferentiated state and underlying endothelial dysfunction. This index was also quantified computationally using the CoA-to-control SM active stress ratio (Equation (18)) and used as a tool to tune the phenotypic modulation stimuli constant C in the G&R model.

2.5.6. WSS and IWS Evolution

WSS ( ϕ ) and IWS ( σ ) are among the most important stress components deriving the rate of stress-mediated growth and remodeling as defined by Equation (26). WSS evolution was calculated on the basis of empirical data from the proximal descending thoracic aorta by implementing peak Doppler velocity readings in the Hagen–Poiseuille law assuming a parabolic velocity profile (Equation (23)). Additionally, IWS in the circumferential direction was quantified from the force equilibrium (Equation (22)). The radius and thickness in this equation were quantified from Doppler b-mode images as described earlier. As mentioned above, mean arterial pressure, i.e., MAP ( p ), was quantified through high-fidelity BP measurements using catheterization and described as a linear function of percentage area obstruction (Figure 4) determined from Doppler ultrasound images.

2.6. Statistical Analysis

Descriptive statistics are presented for continuous variables as the mean ± standard error of the mean (SEM). Unbalanced one-way analysis of variance (ANOVA) was used to assess significant differences across groups using Dunnett’s multiple comparison post hoc analysis with 5% level of significance in GraphPad Prism version 9.5.1, GraphPad Software, San Diego, CA, USA, www.graphpad.com (accessed on 12 October 2022). Pearson’s correlation and linear regression analysis examined the relationships and r-squared for goodness-of-fit. G&R model prediction was validated through impairment index values obtained from rabbits not used for fitting and in terms of percentage error and R-squared. D’Agostino and Pearson, Anderson–Darling, Shapiro–Wilk, and Kolmogorov–Smirnov tests were all used to assess the normality of prediction errors at the 0.05 significance level.

3. Results

Following the use of permanent suture in the prolonged CoA groups, the stenosis was present for ~22 weeks, resulting in adversely elevated BP. In contrast, degradation of the suture in the short and long CoA groups restored aortic diameter and, consequently, also restored associated elevation in BP. Area obstruction levels vary as a function of experimental group as shown in Figure 5 due to rabbit growth and resorption of dissolvable sutures.
Arterial thickening. Figure 6 shows the normalized thickness evolution for the CoA vs. control groups over the course of the protocol. Normalized thickness initially decreased in all rabbit groups and eventually plateaued with body weight after ~150 days of age. Overall, arterial thickening was observed in intermediate to severe CoA (peak-to-peak CoA gradient ≥13 mmHg) when presented for long (~5 weeks) or prolonged durations (~22 weeks). The pattern of arterial thickening was similar to the control group for short presentation of the CoA (i.e., all severities for ~2 weeks). In the mild CoA group (peak-to-peak CoA gradient <13 mmHg); however, thickness for the prolonged CoA group was slightly larger than that of the control at the plateau. After parameter calibration, the G&R model successfully replicated thickness evolutions. Table 2 reports the parameters of the G&R model obtained after the optimization procedure. The adjusted R-squared values for the resulting thickness evolution fits were in the range of 0.98 to 1.00 for all groups.
Arterial stiffening. Figure 7 represents the empirical stress–strain curves quantified in the proximal descending thoracic aorta exposed to adversely elevated BP in CoA groups. Results from each CoA group are compared to the control group in each plot. Overall, the pattern of arterial stiffening was observed by increased stress at the same stretch ratio for groups exposed to a prolonged CoA (~22 weeks). This was also observed in the long CoA groups (~5 weeks), but only for the intermediate and severe CoA (i.e., peak-to-peak CoA gradient ≥13 mmHg). To investigate the rate of stress-mediated remodeling and its effect on the vessel stiffness, constitutive parameters for collagen, SM, and elastin were tuned to fit the empirically measured stress–stretch curves at the end of the protocol (Figure 7, red lines). The constitutive material parameters obtained for c1, c2, and c3 are listed in Table 2.
Vascular dysfunction. Overall, CoA rabbits showed impaired vascular response to half-log increasing dose of PE and ACh as compared to the control group. The extent of vascular impairment showed patterns correlating with both severity and duration of the CoA-induced mechanical stimuli. Figure 8 shows the normalized effective force response and associated AUC (Figure 8D,H,L–O) for CoA vs. control groups. The normalized effective force showed an intact logistic contraction response in the control rabbits peaking at 1.24 ± 0.07, whereas CoA groups peaked to smaller normalized tension values ranging from 0.90 ± 0.06 to 1.09 ± 0.13. This diminished contractile capacity was also observed through AUC among CoA rabbits (Figure 8D,H,L–O). For example, when CoA was at the most severe level, i.e., peak-to-peak CoA gradient ≥20 mmHg, the AUC was significantly smaller than in control rabbits, regardless of duration (Figure 8L). Similarly, when CoA was present for the prolonged period (~22 weeks), significant impairment was observed in AUC, regardless of severity (Figure 8O). A nonsignificant trend toward impaired PE contraction was observed in all other groups.
Figure 9 shows arterial relaxation curves in response to the endothelium-dependent agonist ACh for the descending thoracic aorta proximal to the CoA (exposed to adversely elevated coarctation-induced BP). Overall, control rabbits (black, n = 11) showed intact endothelial function with relaxation to 53% ± 4% of the precontraction force at 10−7 M ACh concentration. However, this vasodilatory action of ACh was impaired in CoA rabbits (gray, n = 5–10). The extent of impaired ACh relaxation correlated with the severity of CoA-induced mechanical stimuli. More specifically, significant decreases in AUC (Figure 9D,H,L–O) were observed in the severe CoA group (peak-to-peak CoA gradient ≥20 mmHg) regardless of the duration of the CoA (Figure 9L). Consequently, the quantified impairment index (CoA-to-control AUC ratio) was 0.52 ± 0.27, 0.55 ± 0.13, and 0.52 ± 0.16 for short, long, and prolonged duration groups, respectively. Rabbits with intermediate CoAs also developed impaired ACh relaxation response that was significant when exposed for >2 weeks, i.e., long and prolonged duration groups (Figure 9N). Associated impairment indices were 0.67 ± 0.19 and 0.85 ± 0.08, respectively. A trend toward impairment was observed in the mild CoA groups, but differences did not reach significance (Figure 9M).
Phenotypic modulation was characterized computationally through its stimulus term ( C λ / λ ^   ) for SM as a function of the current to homeostasis (i.e., control) stretch ratio ( λ / λ ^   ). As indicated in Equation (20), C ^ and C are the two nondimensional constants that determine the phenotypic modulation in SM proliferation and impaired active response. Figure 10 shows active stress per unit relaxed SM normalized to that of the control group (impairment index) and compared to the corresponding empirical measurements of PE and ACh myograph results. Overall, PE contraction and ACh relaxation were both impaired with indices ranging from 0.55 to 0.94 (PE), and from 0.52 to 0.93 (Ach), as shown in Figure 10. Interestingly, this index correlated with the severity and duration of the mechanical stimuli caused by the CoA with values as low as 0.52 observed in the most severe and prolonged CoA groups. Phenotypic modulation stimulus constitutive parameters ( C ) tuned in the G&R model (n = 2) to replicate impairment quantified from wire myography are summarized for each group in Table 3.

4. Discussion

A G&R framework was implemented on the basis of previously developed models [15,16,19,26,40,41] to describe the evolution of arterial wall constituents in response to coarctation-induced mechanical stimuli. Parameters were fit to empirically quantified hypertension precursors among CoA rabbits. Underlying assumptions included thin-walled tube theory, constrained mixture theory, finite elasticity, Hagen–Poiseuille’s law assumptions, and long timescale phenotypic shift of SM from the contractile to synthetic state. Empirical measurements used to tune the model parameters included Doppler ultrasound imaging, uniaxial extension testing, high-fidelity invasive BP measurement, and wire myograph that together quantified the time evolution of arterial thickening, stiffening, and vasoactive dysfunction in response to mechanically induced CoA within a range of severities and durations seen clinically.
In general, good agreement was observed between G&R model prediction and empirical hypertension precursors among nine CoA study groups. Specifically, adjusted R-squared values >0.98 were obtained for thickness evolution (Figure 6), as well as >0.95 for arterial stiffening (Figure 7). There was also a good agreement between G&R predictions for an SM impairment index and empirical SM contraction and relaxation impairment indices obtained from wire myography (Figure 10).
Pathophysiological indices of arterial remodeling such as stiffening, thickening, and vasoactive dysfunction are strong predictors of adverse cardiovascular events and target organ damage of the heart, kidney, and brain [42]. Excessive increases in BP and pulsatility are believed to be the main hemodynamic factors triggering a cycle of events due to damaging effects on vascular function, and contributing to pathological arterial remodeling. Adversely elevated BP and pulsatility cause elevated IWS, causing increased stretch acting on collagen and the extracellular matrix. Additionally, increased cyclic stretch was shown and confirmed in the current study to stimulate a phenotypic shift in SM from the contractile to synthetic state. Moreover, these pathological conditions also affect endothelial function through adversely altered WSS [36,43,44]. Hence, they can be used as surrogates to modify currently debated pressure-based intervention guidelines toward preventing refractory hypertension.
On one hand, identifying the kinetics of pathophysiological responses is challenging in CoA patients mainly due to antihypertensive agents administered, patient heterogeneities making it difficult to differentiate the severity of hemodynamic perturbations experienced, and time once irreversible vascular remodeling occurred. On the other hand, animal models allow for mechanistic investigation of vascular stiffening, thickening, and dysfunction under experimentally controlled conditions to characterize the kinetics of growth and remodeling. Importantly, the tuned model presented here can now be used as a computational tool to identify the severity and duration of the mechanical stimuli leading to hypertension precursors at different phases of rabbit life analogous to that of human. Particularly at younger ages that are difficult to study experimentally, the newborn and weaning phases can be computationally investigated using the G&R model presented. For example, Figure 11A shows the impairment index in the clinically important range of severity, as well as critical age phases. The lines of constant impairment index shown can be used as a map for CoA assessment considering both severity of the narrowing and rabbit age. Particularly, an active response <50%, i.e., impairment index <0.50, indicates nonviable tissue as reported by Lopez et al. [45]. Using this criterion, our results indicate that a peak-to-peak CoA gradient of >35 mmHg, even when exposed shortly, can permanently diminish the active arterial response (Figure 11A). The current model predictions show that such functional impairment is expected for a variety of CoA severity and durations including a short exposure to >30 mmHg peak-to-peak CoA gradient, as well as prolonged exposure to a milder CoA gradient (~10 mmHg). Interestingly, in ages younger than the rabbits studied, i.e., log10(days) < 1.75, i.e., ~70 rabbit days of age (or ~1 human year), model extrapolation showed arterial impairment being strongly correlated to severity of the CoA gradient, thus emphasizing the importance of early intervention to prevent impaired vascular function. This finding further emphasizes the application of such predictive models in clinical decision making, where both severity and duration of the coarctation-induced mechanical stimuli contribute to hypertension precursors developing in the aorta and major branches.
Figure 11B shows active response prediction errors for impairment indices representing NO active response in SM contraction and relaxation. Variability in prediction error was observed which was mainly considered to be inherited from biological variability of the impairment index reported in Figure 10. Statistical assessment of errors showed a normal distribution of across all study groups for both contraction and relaxation impairment indices using four normality tests mentioned in the methods section with p-values > 0.1.
The current results should be interpreted relative to several potential limitations. In particular, high-fidelity hemodynamic simulation is challenging to perform longitudinally as 3D reconstruction of the aorta over the course of the disease requires multiple MRI acquisitions and often computationally heavy [46,47]. We, therefore, quantified WSS and IWS from the Hagen–Poiseuille law and force equilibrium using empirical measurements of local hemodynamics over the course of disease development. Additionally, the large number of unknowns in the formulation derived for G&R requires a large sample size, and that is challenging to achieve due to the time and cost of experimental procedures needed for longitudinal monitoring of hypertension precursors using our preclinical model. Therefore, many of the constitutive parameters were assumed on the basis of the literature followed by tuning of the stress-mediating G&R parameters (Table 3) to replicate empirical results. It is also worth mentioning that the pathophysiological response of the aorta is also affected by concomitant morphological anomalies including the bicuspid aortic valve, aortic arch/isthmus hypoplasia, and long-segment coarctation. The associated effect on local hemodynamic and aortic remodeling remains to be investigated in future work.

5. Conclusions

Results showed that the implemented G&R model could accurately represent hypertension precursors among the studied groups of CoA rabbits. Therefore, the G&R model allows for computational identification of the adverse vascular response as a function of the severity and duration of coarctation-induced stimuli. Clinical translation of this predictive computational model can help in modifying current intervention thresholds and preventing irreversible hypertension precursors that seem to persist after CoA treatment using current clinical guidance.

Author Contributions

Conceptualization, A.G., J.F.L. and R.K.W.; methodology, A.G.; software, A.G.; validation, A.G., J.F.L., M.M. and S.H.A.; formal analysis, A.G.; investigation, A.G.; resources, J.F.L.; data curation, A.G. and M.M.; writing—original draft preparation, A.G. and J.F.L.; writing—review and editing, A.G., M.M., E.S.H.I., C.C., B.J.T., R.K.W., S.H.A. and J.F.L.; visualization, A.G.; supervision, M.M, E.S.H.I., C.C., B.J.T., R.K.W., S.H.A. and J.F.L.; project administration, J.F.L.; funding acquisition, J.F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NIH grant number [R01-HL142955] and [R15-HL096096].

Institutional Review Board Statement

The animal study protocol was approved by the Institutional Review Board of Medical College of Wisconsin (protocol code AUA00001175 approved 18 April 2009 through 16 August 2023).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors gratefully acknowledge the technical support from Leanne Harman, Lindsey Kalvin, Jamasp Azarnoosh, Thomas Eddinger, and Chalani Ellepola.

Conflicts of Interest

The terms of graduate student funding for Arash Ghorbannia, the corresponding author of the current work, have been reviewed and approved by the Medical College of Wisconsin and Marquette University in accordance with their policies on objectivity in research. The authors have no other financial disclosures or conflict of interest to report.

Abbreviations

AAAscending aorta
BPBlood pressure
CFPGContinuous-flow pressure gradient
CoACoarctation of the aorta
IWSIntramural wall stress
MBEModified Bernoulli equation
PC-MRIPhase-contrast magnetic resonance image
ROMReduced-order model
SBESimplified Bernoulli equation
SBE-RPSBE minus recovered pressure term
WSSWall shear stress

References

  1. Ferencz, C.; Rubin, J.D.; McCarter, R.J.; Brenner, J.I.; Neill, C.A.; Perry, L.W.; Hepner, S.I.; Downing, J.W. Congenital Heart Disease: Prevalence at Livebirth: The Baltimore-Washington Infant Study. Am. J. Epidemiol. 1985, 121, 31–36. [Google Scholar] [CrossRef]
  2. Torok, R.D.; Campbell, M.J.; Fleming, G.A.; Hill, K.D. Coarctation of the Aorta: Management from Infancy to Adulthood. World J. Cardiol. 2015, 7, 765–775. [Google Scholar] [CrossRef]
  3. Feltes, T.F.; Bacha, E.; Beekman, R.H.; Cheatham, J.P.; Feinstein, J.A.; Gomes, A.S.; Hijazi, Z.M.; Ing, F.F.; de Moor, M.; Morrow, W.R.; et al. Indications for Cardiac Catheterization and Intervention in Pediatric Cardiac Disease. Circulation 2011, 123, 2607–2652. [Google Scholar] [CrossRef]
  4. Cohen, M.; Fuster, V.; Steele, P.M.; Driscoll, D.; McGoon, D.C. Coarctation of the Aorta. Long-Term Follow-up and Prediction of Outcome after Surgical Correction. Circulation 1989, 80, 840–845. [Google Scholar] [CrossRef]
  5. Sendzikaite, S.; Sudikiene, R.; Tarutis, V.; Lubaua, I.; Silis, P.; Rybak, A.; Jankauskiene, A.; Litwin, M. Prevalence of Arterial Hypertension, Hemodynamic Phenotypes, and Left Ventricular Hypertrophy in Children after Coarctation Repair: A Multicenter Cross-Sectional Study. Pediatr. Nephrol. 2020, 35, 2147–2155. [Google Scholar] [CrossRef]
  6. LaDisa, J.F., Jr.; Bozdag, S.; Olson, J.; Ramchandran, R.; Kersten, J.R.; Eddinger, T.J. Gene Expression in Experimental Aortic Coarctation and Repair: Candidate Genes for Therapeutic Intervention? PLoS ONE 2015, 10, e0133356. [Google Scholar] [CrossRef]
  7. LaDisa, J.F.J.; Tomita-Mitchell, A.; Stamm, K.; Bazan, K.; Mahnke, D.K.; Goetsch, M.A.; Wegter, B.J.; Gerringer, J.W.; Repp, K.; Palygin, O.; et al. Human Genotyping and an Experimental Model Reveal NPR-C as a Possible Contributor to Morbidity in Coarctation of the Aorta. Physiol. Genom. 2019, 51, 177–185. [Google Scholar] [CrossRef]
  8. Menon, A.; Wendell, D.C.; Wang, H.; Eddinger, T.J.; Toth, J.M.; Dholakia, R.J.; Larsen, P.M.; Jensen, E.S.; LaDisa, J.F., Jr. A Coupled Experimental and Computational Approach to Quantify Deleterious Hemodynamics, Vascular Alterations, and Mechanisms of Long-Term Morbidity in Response to Aortic Coarctation. J. Pharmacol. Toxicol. Methods 2012, 65, 18–28. [Google Scholar] [CrossRef]
  9. Menon, A.; Eddinger, T.J.; Wang, H.; Wendell, D.C.; Toth, J.M.; LaDisa, J.F., Jr. Altered Hemodynamics, Endothelial Function, and Protein Expression Occur with Aortic Coarctation and Persist after Repair. Am. J. Physiol. Heart Circ. Physiol. 2012, 303, H1304-18. [Google Scholar] [CrossRef]
  10. Ghorbannia, A.; Ellepola, C.D.; Woods, R.K.; Ibrahim, E.-S.H.; Maadooliat, M.; Ramirez, H.M.; LaDisa, J.F.J. Clinical, Experimental, and Computational Validation of a New Doppler-Based Index for Coarctation Severity Assessment. J. Am. Soc. Echocardiogr. 2022, 35, 1311–1321. [Google Scholar] [CrossRef]
  11. Panzer, J.; Bové, T.; Vandekerckhove, K.; De Wolf, D. Hypertension after Coarctation Repair-a Systematic Review. Transl. Pediatr. 2022, 11, 270–279. [Google Scholar] [CrossRef]
  12. Valentín, A.; Holzapfel, G.A. Constrained Mixture Models as Tools for Testing Competing Hypotheses in Arterial Biomechanics: A Brief Survey. Mech. Res. Commun. 2012, 42, 126–133. [Google Scholar] [CrossRef]
  13. Valentín, A.; Humphrey, J.D. Evaluation of Fundamental Hypotheses Underlying Constrained Mixture Models of Arterial Growth and Remodelling. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2009, 367, 3585–3606. [Google Scholar] [CrossRef]
  14. Humphrey, J.D. Mechanisms of Arterial Remodeling in Hypertension. Hypertension 2008, 52, 195–200. [Google Scholar] [CrossRef]
  15. Cyron, C.J.; Humphrey, J.D. Vascular Homeostasis and the Concept of Mechanobiological Stability. Int. J. Eng. Sci. 2014, 85, 203–223. [Google Scholar] [CrossRef]
  16. Humphrey, J.D.; Rajagopal, K.R. A Constrained Mixture Model for Growth and Remodeling of Soft Tissues. Math. Model. Methods Appl. Sci. 2002, 12, 407–430. [Google Scholar] [CrossRef]
  17. Teng, Z.; Yuan, J.; Feng, J.; Zhang, Y.; Brown, A.J.; Wang, S.; Lu, Q.; Gillard, J.H. The Influence of Constitutive Law Choice Used to Characterise Atherosclerotic Tissue Material Properties on Computing Stress Values in Human Carotid Plaques. J. Biomech. 2015, 48, 3912–3921. [Google Scholar] [CrossRef]
  18. Ogden, R.W. Non-Linear Elastic Deformations; Dover Civil and Mechanical Engineering; Dover Publications: Mineola, NY, USA, 1997; ISBN 9780486696485. [Google Scholar]
  19. Wilson, J.S.; Baek, S.; Humphrey, J.D. Parametric Study of Effects of Collagen Turnover on the Natural History of Abdominal Aortic Aneurysms. Proc. R. Soc. A Math. Phys. Eng. Sci. 2013, 469, 20120556. [Google Scholar] [CrossRef]
  20. Wu, J.; Shadden, S.C. Coupled Simulation of Hemodynamics and Vascular Growth and Remodeling in a Subject-Specific Geometry. Ann. Biomed. Eng. 2015, 43, 1543–1554. [Google Scholar] [CrossRef]
  21. Rhodin, J.A.G. Architecture of the Vessel Wall. In Comprehensive Physiology; Wiley: New York, NY, USA, 2014; pp. 1–31. ISBN 978-0-470-65071-4. [Google Scholar]
  22. Halayko, A.J.; Solway, J. Invited Review: Molecular Mechanisms of Phenotypic Plasticity in Smooth Muscle Cells. J. Appl. Physiol. 2001, 90, 358–368. [Google Scholar] [CrossRef]
  23. Albinsson, S.; Bhattachariya, A.; Hellstrand, P. Stretch-Dependent Smooth Muscle Differentiation in the Portal Vein-Role of Actin Polymerization, Calcium Signaling, and MicroRNAs. Microcirculation 2014, 21, 230–238. [Google Scholar] [CrossRef]
  24. Westerhof, N.; Stergiopulos, N.; Noble, M.I.M. Snapshots of Hemodynamics an Aid for Clinical Research and Graduate Education; Springer: New York, NY, USA, 2005. [Google Scholar]
  25. Van Varik, B.; Rennenberg, R.; Reutelingsperger, C.; Kroon, A.; de Leeuw, P.; Schurgers, L. Mechanisms of Arterial Remodeling: Lessons from Genetic Diseases. Front. Genet. 2012, 3, 290. [Google Scholar] [CrossRef]
  26. Lindström, S.B.; Stålhand, J.; Klarbring, A. Modeling of the Mechanobiological Adaptation in Muscular Arteries. Eur. J. Mech.-A/Solids 2017, 64, 165–177. [Google Scholar] [CrossRef]
  27. Figueroa, C.A.; Baek, S.; Taylor, C.A.; Humphrey, J.D. A Computational Framework for Fluid-Solid-Growth Modeling in Cardiovascular Simulations. Comput. Methods Appl. Mech. Eng. 2009, 198, 3583–3602. [Google Scholar] [CrossRef]
  28. Baek, S.; Valentín, A.; Humphrey, J.D. Biochemomechanics of Cerebral Vasospasm and Its Resolution: II. Constitutive Relations and Model Simulations. Ann. Biomed. Eng. 2007, 35, 1498–1509. [Google Scholar] [CrossRef]
  29. Rachev, A.; Hayashi, K. Theoretical Study of the Effects of Vascular Smooth Muscle Contraction on Strain and Stress Distributions in Arteries. Ann. Biomed. Eng. 1999, 27, 459–468. [Google Scholar] [CrossRef]
  30. Ogden, R.W. Nonlinear Elasticity, Anisotropy, Material Stability and Residual Stresses in Soft Tissue BT—Biomechanics of Soft Tissue in Cardiovascular Systems; Holzapfel, G.A., Ogden, R.W., Eds.; Springer: Vienna, Austria, 2003; pp. 65–108. ISBN 978-3-7091-2736-0. [Google Scholar]
  31. Taber, L.A. A Model for Aortic Growth Based on Fluid Shear and Fiber Stresses. J. Biomech. Eng. 1998, 120, 348–354. [Google Scholar] [CrossRef]
  32. Rachev, A. A Model of Arterial Adaptation to Alterations in Blood Flow. J. Elast. Phys. Sci. Solids 2000, 61, 83–111. [Google Scholar] [CrossRef]
  33. Conn, A.R.; Gould, N.I.M.; Toint, P. A Globally Convergent Augmented Lagrangian Algorithm for Optimization with General Constraints and Simple Bounds. SIAM J. Numer. Anal. 1991, 28, 545–572. [Google Scholar] [CrossRef]
  34. Holzapfel, G.A.; Gasser, T.C.; Ogden, R.W. A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models. J. Elast. Phys. Sci. Solids 2000, 61, 1–48. [Google Scholar] [CrossRef]
  35. Touyz, R.M.; Alves-Lopes, R.; Rios, F.J.; Camargo, L.L.; Anagnostopoulou, A.; Arner, A.; Montezano, A.C. Vascular Smooth Muscle Contraction in Hypertension. Cardiovasc. Res. 2018, 114, 529–539. [Google Scholar] [CrossRef]
  36. DePaola, N.; Gimbrone, M.A., Jr.; Davies, P.F.; Dewey, C.F. Vascular Endothelium Responds to Fluid Shear Stress Gradients. Arterioscler. Thromb. 1992, 12, 1254–1257. [Google Scholar] [CrossRef]
  37. Stout, K.K.; Daniels, C.J.; Aboulhosn, J.A.; Bozkurt, B.; Broberg, C.S.; Colman, J.M.; Crumb, S.R.; Dearani, J.A.; Fuller, S.; Gurvitz, M.; et al. 2018 AHA/ACC Guideline for the Management of Adults With Congenital Heart Disease: A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines. Circulation 2019, 139, e698–e800. [Google Scholar] [CrossRef]
  38. Masoud, I.; Shapiro, F.; Kent, R.; Moses, A. A Longitudinal Study of the Growth of the New Zealand White Rabbit: Cumulative and Biweekly Incremental Growth Rates for Body Length, Body Weight, Femoral Length, and Tibial Length. J. Orthop. Res. Off. Publ. Orthop. Res. Soc. 1986, 4, 221–231. [Google Scholar] [CrossRef]
  39. Dutta, S.; Sengupta, P. Rabbits and Men: Relating Their Ages. J. Basic Clin. Physiol. 2018, 29, 427–435. [Google Scholar] [CrossRef]
  40. Humphrey, J.D. Towards a Theory of Vascular Growth and Remodeling BT—Mechanics of Biological Tissue; Holzapfel, G.A., Ogden, R.W., Eds.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 3–15. ISBN 978-3-540-31184-3. [Google Scholar]
  41. Latorre, M.; Humphrey, J.D. A Mechanobiologically Equilibrated Constrained Mixture Model for Growth and Remodeling of Soft Tissues. Z. Fur. Angew. Math. Und Mech. 2018, 98, 2048–2071. [Google Scholar] [CrossRef]
  42. Laurent, S.; Boutouyrie, P. Arterial Stiffness and Stroke in Hypertension: Therapeutic Implications for Stroke Prevention. CNS Drugs 2005, 19, 1–11. [Google Scholar] [CrossRef]
  43. Tardy, Y.; Resnick, N.; Nagel, T.; Gimbrone, M.A.; Dewey, C.F. Shear Stress Gradients Remodel Endothelial Monolayers in Vitro via a Cell Proliferation-Migration-Loss Cycle. Arterioscler. Thromb. Vasc. Biol. 1997, 17, 3102–3106. [Google Scholar] [CrossRef]
  44. Gimbrone, M.A.J.; García-Cardeña, G. Vascular Endothelium, Hemodynamics, and the Pathobiology of Atherosclerosis. Cardiovasc. Pathol. Off. J. Soc. Cardiovasc. Pathol. 2013, 22, 9–15. [Google Scholar] [CrossRef]
  45. López, D.; Rodríguez-Sinovas, A.; Agulló, L.; Inserte, J.; Cabestrero, A.; García-Dorado, D. Acidic Reoxygenation Protects against Endothelial Dysfunction in Rat Aortic Rings Submitted to Simulated Ischemia. Am. J. Physiol.-Heart Circ. Physiol. 2008, 295, H2409–H2416. [Google Scholar] [CrossRef]
  46. Ghorbanniahassankiadeh, A.; Marks, D.S.; LaDisa, J.F., Jr. Correlation of Computational Instantaneous Wave-Free Ratio With Fractional Flow Reserve for Intermediate Multivessel Coronary Disease. J. Biomech. Eng. 2021, 143, 051011. [Google Scholar] [CrossRef]
  47. LaDisa, J.F.; Ghorbannia, A.; Marks, D.S.; Mason, P.; Otake, H. Advancements and Opportunities in Characterizing Patient-Specific Wall Shear Stress Imposed by Coronary Artery Stenting. Fluids 2022, 7, 325. [Google Scholar] [CrossRef]
Figure 1. Experimental protocol to characterize hypertension precursors in response to CoA in rabbits. Male New Zealand White rabbits were randomly selected for a control group or designated to undergo descending thoracic CoA (nine study groups: three severities and three durations) via left thoracotomy in the third intercostal space. Weekly Doppler ultrasound imaging characterized local hemodynamic and morphology adaptation to CoA before PC-MRI prior to the final week. At the end of the protocol, invasive characterization was conducted for blood pressure measurement, and the aorta was dissected to quantify arterial stiffening and dysfunction through uniaxial extension testing and myograph analysis, respectively.
Figure 1. Experimental protocol to characterize hypertension precursors in response to CoA in rabbits. Male New Zealand White rabbits were randomly selected for a control group or designated to undergo descending thoracic CoA (nine study groups: three severities and three durations) via left thoracotomy in the third intercostal space. Weekly Doppler ultrasound imaging characterized local hemodynamic and morphology adaptation to CoA before PC-MRI prior to the final week. At the end of the protocol, invasive characterization was conducted for blood pressure measurement, and the aorta was dissected to quantify arterial stiffening and dysfunction through uniaxial extension testing and myograph analysis, respectively.
Biomedicines 11 01817 g001
Figure 2. Quantification of hemodynamics (A) and morphology (B) in a rabbit model of CoA. Doppler ultrasound follow-up was performed weekly to quantify the temporal evolution of hemodynamic changes (A, spectral Doppler) and morphology (B, Doppler b-mode) in response to the CoA.
Figure 2. Quantification of hemodynamics (A) and morphology (B) in a rabbit model of CoA. Doppler ultrasound follow-up was performed weekly to quantify the temporal evolution of hemodynamic changes (A, spectral Doppler) and morphology (B, Doppler b-mode) in response to the CoA.
Biomedicines 11 01817 g002
Figure 3. Representative thickness and body weight evolution over time for a CoA rabbit. Thickness (squares) was quantified as a function of age and through weekly ultrasound imaging at follow-up where 2D sagittal sections of the aorta were identified in b-mode images, and triplicate thickness measurements were performed with mean values shown at each timepoint. Body weight evolution was extrapolated using a sigmoid fit (dashed line) to the measured body weights (circles).
Figure 3. Representative thickness and body weight evolution over time for a CoA rabbit. Thickness (squares) was quantified as a function of age and through weekly ultrasound imaging at follow-up where 2D sagittal sections of the aorta were identified in b-mode images, and triplicate thickness measurements were performed with mean values shown at each timepoint. Body weight evolution was extrapolated using a sigmoid fit (dashed line) to the measured body weights (circles).
Biomedicines 11 01817 g003
Figure 4. High-fidelity BP measurements proximal to CoA via catheterization. Rabbits (n = 55) were anesthetized for measurement of the BP above the coarctation prior to tissue harvest using high-fidelity pressure transducers (Harvard Apparatus, Holliston, MA, USA) attached to 0.86 mm ID noncompliant fluid-filled catheters inserted into the right common carotid artery and advanced to the aortic arch. Systolic blood pressure (SBP: black circles), mean arterial pressure (MAP), and diastolic blood pressure (DBP: open circles) were then quantified from the BP waveforms. Linear regression was then used to mathematically model MAP (dashed line) as a function of percent area obstruction (formula).
Figure 4. High-fidelity BP measurements proximal to CoA via catheterization. Rabbits (n = 55) were anesthetized for measurement of the BP above the coarctation prior to tissue harvest using high-fidelity pressure transducers (Harvard Apparatus, Holliston, MA, USA) attached to 0.86 mm ID noncompliant fluid-filled catheters inserted into the right common carotid artery and advanced to the aortic arch. Systolic blood pressure (SBP: black circles), mean arterial pressure (MAP), and diastolic blood pressure (DBP: open circles) were then quantified from the BP waveforms. Linear regression was then used to mathematically model MAP (dashed line) as a function of percent area obstruction (formula).
Biomedicines 11 01817 g004
Figure 5. Percentage area obstruction vs. days after CoA surgery. CoA morphology was evaluated on a weekly basis noninvasively using ultrasound imaging. Proximal aorta and CoA diameters were quantified in triplicate, and the average value was plotted (scatter) at each date. The percentage area obstruction (%AO = 100 × (Dprox2 − DCoA2)/Dprox2) was, therefore, quantified and plotted vs. days after surgery. For the dissolving sutures (i.e., short and long duration CoA groups), a plateau followed by exponential decay function was fitted to the %AO, while, for the permanent suture (i.e., prolonged CoA group), an exponential function was used to fit the data.
Figure 5. Percentage area obstruction vs. days after CoA surgery. CoA morphology was evaluated on a weekly basis noninvasively using ultrasound imaging. Proximal aorta and CoA diameters were quantified in triplicate, and the average value was plotted (scatter) at each date. The percentage area obstruction (%AO = 100 × (Dprox2 − DCoA2)/Dprox2) was, therefore, quantified and plotted vs. days after surgery. For the dissolving sutures (i.e., short and long duration CoA groups), a plateau followed by exponential decay function was fitted to the %AO, while, for the permanent suture (i.e., prolonged CoA group), an exponential function was used to fit the data.
Biomedicines 11 01817 g005
Figure 6. Thickness evolution for CoA rabbits (scatter) and G&R model predictions (lines). Coarctation-induced arterial thickening (shades of gray) was quantified from Doppler ultrasound b-mode images. To compare with control group, thicknesses were normalized to body weight (black). CoA groups included all combinations of severity (mild, intermediate, and severe) and duration (short, long, and prolonged) studied. G&R model constitutive parameters were tuned to fit empirical measurements (n = 2/group) with adjusted R2 characterizing goodness of fit.
Figure 6. Thickness evolution for CoA rabbits (scatter) and G&R model predictions (lines). Coarctation-induced arterial thickening (shades of gray) was quantified from Doppler ultrasound b-mode images. To compare with control group, thicknesses were normalized to body weight (black). CoA groups included all combinations of severity (mild, intermediate, and severe) and duration (short, long, and prolonged) studied. G&R model constitutive parameters were tuned to fit empirical measurements (n = 2/group) with adjusted R2 characterizing goodness of fit.
Biomedicines 11 01817 g006
Figure 7. Empirical stress–stretch curves for CoA (gray) vs. control (black) groups and G&R model prediction (red). Error bars representing the SEM were quantified and averaged at each strain energy level discretized over the range of 0 to 3 MPa. The constitutive material parameters identified are listed in Table 2.
Figure 7. Empirical stress–stretch curves for CoA (gray) vs. control (black) groups and G&R model prediction (red). Error bars representing the SEM were quantified and averaged at each strain energy level discretized over the range of 0 to 3 MPa. The constitutive material parameters identified are listed in Table 2.
Biomedicines 11 01817 g007
Figure 8. Comparative representation of active contraction response to PE in CoA vs. control. Dose response curves (myography) by increasing half-log dose of phenylephrine (PE: administered every 5 min). The measured contractile force was normalized to maximum K+ contraction achieved for each segment. The nine CoA groups consist of all combinations of severity (mild, intermediate, and severe: different gray levels) and duration (short, long, and prolonged: different line styles) studied, compared to the control group (black) in each plot (AC,EG,IK). Bar plots (D,H,LO) represent the area under each dose–response curve (AUC). Error bars represent the SEM. * Significant difference compared to the control group in a two-tailed Student’s t-test with a 5% significance level, assuming similar variation among the groups.
Figure 8. Comparative representation of active contraction response to PE in CoA vs. control. Dose response curves (myography) by increasing half-log dose of phenylephrine (PE: administered every 5 min). The measured contractile force was normalized to maximum K+ contraction achieved for each segment. The nine CoA groups consist of all combinations of severity (mild, intermediate, and severe: different gray levels) and duration (short, long, and prolonged: different line styles) studied, compared to the control group (black) in each plot (AC,EG,IK). Bar plots (D,H,LO) represent the area under each dose–response curve (AUC). Error bars represent the SEM. * Significant difference compared to the control group in a two-tailed Student’s t-test with a 5% significance level, assuming similar variation among the groups.
Biomedicines 11 01817 g008
Figure 9. Comparative representation of CoA vs. control active relaxation response to ACh. Dose–response curves (myography) by increasing the half-log dose of acetylcholine (ACh: administered every 5 min). Relaxation was measured relative to the initial contraction (EC50 concentration). The nine CoA groups consisted of all combinations of severity (mild, intermediate, and severe: different gray levels) and duration (short, long, and prolonged: different line styles) studied and compared to the control group (black) in each plot (AC,EG,IK). Bar plots (D,H,LO) represent area under each dose–response curve (AUC). Error bars represent SEM. * Significant difference compared to the control group in a two-tailed Student’s t-test with 5% significance level assuming similar variation among the groups.
Figure 9. Comparative representation of CoA vs. control active relaxation response to ACh. Dose–response curves (myography) by increasing the half-log dose of acetylcholine (ACh: administered every 5 min). Relaxation was measured relative to the initial contraction (EC50 concentration). The nine CoA groups consisted of all combinations of severity (mild, intermediate, and severe: different gray levels) and duration (short, long, and prolonged: different line styles) studied and compared to the control group (black) in each plot (AC,EG,IK). Bar plots (D,H,LO) represent area under each dose–response curve (AUC). Error bars represent SEM. * Significant difference compared to the control group in a two-tailed Student’s t-test with 5% significance level assuming similar variation among the groups.
Biomedicines 11 01817 g009
Figure 10. Coarctation-induced active vascular response normalized to the control group. Area under the active dose–response curves (AUC) for each CoA group normalized to that of the control. Dose–response curves were quantified via wire myography and a half-log increase in phenylephrine (PE, gray) and acetylcholine (Ach, pink) to characterize NO active response in SM contraction and relaxation, respectively. A similar index was also quantified by taking the ratio of active stress per unit relaxed SM for CoA relative to control groups using the G&R framework (green) for two rabbits per group. G&R constitutive parameters were tuned to achieve minimum error to empirically quantified active response, i.e., PE and Ach. No significant differences were observed between results from empirical and G&R through two-tailed Student’s t-test with a 5% significance level.
Figure 10. Coarctation-induced active vascular response normalized to the control group. Area under the active dose–response curves (AUC) for each CoA group normalized to that of the control. Dose–response curves were quantified via wire myography and a half-log increase in phenylephrine (PE, gray) and acetylcholine (Ach, pink) to characterize NO active response in SM contraction and relaxation, respectively. A similar index was also quantified by taking the ratio of active stress per unit relaxed SM for CoA relative to control groups using the G&R framework (green) for two rabbits per group. G&R constitutive parameters were tuned to achieve minimum error to empirically quantified active response, i.e., PE and Ach. No significant differences were observed between results from empirical and G&R through two-tailed Student’s t-test with a 5% significance level.
Biomedicines 11 01817 g010
Figure 11. G&R model prediction (A) and prediction error (B) relative to the measured impaired vascular response. (A) Contour plot showing impairment index quantified by taking the ratio of active stress per unit relaxed SM for CoA relative to control using the G&R framework. Contour lines indicate regions of constant impairment indices. CoA severity, shown on the vertical axis, is characterized through peak-to-peak trans-coarctation BP gradient (BPG), whereas duration is represented on a logarithmic scale of age in days. (B) Model error relative to measured impairment index, i.e., area under the active dose–response curves (AUC) for each CoA group normalized to that of the control. Dose–response curves were quantified via wire myography and half-log increase of phenylephrine (PE, gray) and acetylcholine (Ach, pink) to characterize NO active response in SM contraction and relaxation, respectively. Model error was then quantified as a percentage relative to the empirical PE and Ach impairment indices. Rabbits used in the tuning process were excluded for model prediction error.
Figure 11. G&R model prediction (A) and prediction error (B) relative to the measured impaired vascular response. (A) Contour plot showing impairment index quantified by taking the ratio of active stress per unit relaxed SM for CoA relative to control using the G&R framework. Contour lines indicate regions of constant impairment indices. CoA severity, shown on the vertical axis, is characterized through peak-to-peak trans-coarctation BP gradient (BPG), whereas duration is represented on a logarithmic scale of age in days. (B) Model error relative to measured impairment index, i.e., area under the active dose–response curves (AUC) for each CoA group normalized to that of the control. Dose–response curves were quantified via wire myography and half-log increase of phenylephrine (PE, gray) and acetylcholine (Ach, pink) to characterize NO active response in SM contraction and relaxation, respectively. Model error was then quantified as a percentage relative to the empirical PE and Ach impairment indices. Rabbits used in the tuning process were excluded for model prediction error.
Biomedicines 11 01817 g011
Table 1. Constitutive parameters in the derived G&R framework.
Table 1. Constitutive parameters in the derived G&R framework.
ParameterDescriptionReferenceValue
ν c Collagen turnover rateRhobin et al. [21]1/80 [1/days]
ν s SM turnover rateRhobin et al. [21]1/80 [1/days]
Q e Initial elastin decay kineticsHalayko et al. [22]1
q e Remodeled elastin decay kineticsHalayko et al. [22]0
G ^ m SM pre-stretchHolzapfel et al. [34]1.3
S M Maximum stress capacity of SM in the fully contractile stateConn et al. [33]100 [kPa]
C ^ Phenotypic modulation stimulus parameter at normotensive conditionLindstrom et al. [26]0.8326
C Phenotypic modulation stimulus parameter at hypertensive conditionTouyz et al. [35][0.1, 100]
c 1 k Neo-Hookean   material   parameter   for   constituent   k c , s , e Wu et al. [20]
c 2 k Exponential   material   parameter   for   constituent   k c , s , e Wu et al. [20]
c 3 k Exponential   material   parameter   for   constituent   k c , s , e Wu et al. [20]
G ^ c Pre-stretch for collagenRachev et al. [32]1.08
G ^ e Pre-stretch for elastin Rachev et al. [32]1.4
K σ c G&R collagen kinetic parameter for hoop stress DePaola et al. [36] [ 2.0 ,   18 ] × 10 6
K ϕ c G&R collagen kinetic parameter for WSSDePaola et al. [36] [ 0.8 ,   8 ] × 10 6
K σ s G&R SM kinetic parameter for hoop stressDePaola et al. [36] [ 2.0 ,   18 ] × 10 6
K ϕ s G&R SM kinetic parameter for WSSDePaola et al. [36] [ 0.8 ,   8 ] × 10 6
m 0 c Basal net production rate of collagenHalayko et al. [22][0, 10]
m 0 s Basal net production rate of SMHalayko et al. [22][0, 10]
m 0 e Basal net production rate of elastinHalayko et al. [22]0
Table 2. G&R constitutive material parameters for collagen, SM, and elastin obtained from CoA rabbits (n = 2 per group).
Table 2. G&R constitutive material parameters for collagen, SM, and elastin obtained from CoA rabbits (n = 2 per group).
Rabbit GroupConstitutive Parameters
Duration (Weeks)Severity (Peak-to-Peak CoA Gradient)Collagen *Smooth Muscle *Elastin **
c1 ***c2c3c1c2c3c1c2c3
Short
(~2)
Mild (<13 )0.0597 ± 9111.7 ± 4.80.011.5 ± 3.32.3 ± 1.456.2 ± 210.0NA
Intermediate (13–20)0.0646 ± 1016.40 ± 2.10.013.2 ± 4.33.7 ± 1.867.4 ± 180.0NA
Severe (≥20)0.0613 ± 947.70 ± 2.00.013.0 ± 3.91.9 ± 1.058.7 ± 230.0NA
Long
(~5)
Mild (<13)0.0629 ± 9710.5 ± 3.80.010.1 ± 4.23.7 ± 1.959.3 ± 80.0NA
Intermediate (13–20)0.0673 ± 1139.70 ± 4.00.014.2 ± 4.84.8 ± 2.163.8 ± 160.0NA
Severe (≥20)0.0716 ± 11714.9 ± 5.00.018.9 ± 4.37.2 ± 2.872.9 ± 310.0NA
Prolonged
(~22)
Mild (<13)0.0692 ± 8812.7 ± 3.60.013.8 ± 5.111.7 ± 4.865.4 ± 210.0NA
Intermediate (13–20)0.0708 ± 12213.2 ± 4.50.015.1 ± 3.811.7 ± 4.861.4 ± 170.0NA
Severe (≥20)0.0721 ± 4214.0 ± 5.20.015.2 ± 4.211.7 ± 4.868.1 ± 240.0NA
* Collagen and smooth muscle were described by an exponential strain energy density function [30] that requires c1 = 0. ** Elastin was described by a neo-Hookean [17] strain energy density function that requires c2 = 0 and also drops c3 from Equation (21). *** c1 and c2 are in kPa, and c3 is nondimensional. Values represent the average ± SD. Pressures are in units of mmHg.
Table 3. G&R constitutive parameters determining the kinetics of stress-mediated response to CoA.
Table 3. G&R constitutive parameters determining the kinetics of stress-mediated response to CoA.
Rabbit GroupConstitutive Parameters
Duration (Weeks)Severity (Peak-to-Peak CoA Gradient)CollagenSmooth Muscle
K σ (×10−6) K ϕ (×10−6) K σ (×10−6) K ϕ (×10−6)C
Short
(~2)
Mild (<13 )3.10 ± 2.0 *1.80 ± 1.25.30 ± 3.03.90 ± 1.510.1 ± 3.0
Intermediate (13–20)5.20 ± 1.12.90 ± 1.210.2 ± 4.54.10 ± 2.08.80 ± 2.1
Severe (≥20)3.70 ± 0.91.70 ± 0.57.10 ± 2.74.50 ± 1.810.0 ± 2.8
Long
(~5)
Mild (<13)4.30 ± 1.62.30 ± 1.25.60 ± 3.15.60 ± 2.79.60 ± 2.5
Intermediate (13–20)4.10 ± 2.12.00 ± 0.87.60 ± 2.74.60 ± 1.111.6 ± 1.7
Severe (≥20)5.00 ± 1.73.50 ± 1.39.60 ± 1.53.60 ± 1.68.60 ± 2.1
Prolonged
(~22)
Mild (<13)5.40 ± 2.52.80 ± 1.25.60 ± 2.23.60 ± 1.17.60 ± 2.7
Intermediate (13–20)4.60 ± 1.33.10 ± 1.912.60 ± 2.85.60 ± 2.99.60 ± 3.0
Severe (≥20)6.30 ± 2.23.40 ± 1.710.60 ± 3.04.60 ± 2.710.6 ± 1.9
Pressures are represented in mmHg. * Values represent average ± SD.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ghorbannia, A.; Maadooliat, M.; Woods, R.K.; Audi, S.H.; Tefft, B.J.; Chiastra, C.; Ibrahim, E.S.H.; LaDisa, J.F. Aortic Remodeling Kinetics in Response to Coarctation-Induced Mechanical Perturbations. Biomedicines 2023, 11, 1817. https://doi.org/10.3390/biomedicines11071817

AMA Style

Ghorbannia A, Maadooliat M, Woods RK, Audi SH, Tefft BJ, Chiastra C, Ibrahim ESH, LaDisa JF. Aortic Remodeling Kinetics in Response to Coarctation-Induced Mechanical Perturbations. Biomedicines. 2023; 11(7):1817. https://doi.org/10.3390/biomedicines11071817

Chicago/Turabian Style

Ghorbannia, Arash, Mehdi Maadooliat, Ronald K. Woods, Said H. Audi, Brandon J. Tefft, Claudio Chiastra, El Sayed H. Ibrahim, and John F. LaDisa. 2023. "Aortic Remodeling Kinetics in Response to Coarctation-Induced Mechanical Perturbations" Biomedicines 11, no. 7: 1817. https://doi.org/10.3390/biomedicines11071817

APA Style

Ghorbannia, A., Maadooliat, M., Woods, R. K., Audi, S. H., Tefft, B. J., Chiastra, C., Ibrahim, E. S. H., & LaDisa, J. F. (2023). Aortic Remodeling Kinetics in Response to Coarctation-Induced Mechanical Perturbations. Biomedicines, 11(7), 1817. https://doi.org/10.3390/biomedicines11071817

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop