Search

Ventricular suction is a common adverse event in ventricular assist device (VAD) patients and can be due to multiple underlying causes. The aim of this study is to analyze the potential of different therapeutic interventions to mitigate suction events induced by different pathophysiological conditions. To do so, a suction module was embedded in a cardiovascular hybrid (hydraulic-computational) simulator reproducing the entire cardiovascular system. An HVAD system (Medtronic) was connected between a compliant ventricular apex and a simulated aorta. Starting from a patient profile with severe dilated cardiomyopathy, four different pathophysiological conditions leading to suction were simulated: hypovolemia (blood volume: -900 ml), right ventricular failure (contractility -70%), hypotension (systemic vascular resistance: 8.3 Wood Units), and tachycardia (heart rate:185 bpm). Different therapeutic interventions such as volume infusion, ventricular contractility increase, vasoconstriction, heart rate increase, and pump speed reduction were simulated. Their effects were compared in terms of general hemodynamics and suction mitigation. Each intervention elicited a different effect on the hemodynamics for every pathophysiological condition. Pump speed reduction mitigated suction but did not ameliorate the hemodynamics. Administering volume and inducing a systemic vasoconstriction were the most efficient interventions in both improving the hemodynamics and mitigating suction. When simulating volume infusion, the cardiac powers increased, respectively, by 38%, 25%, 42%, and 43% in the case of hypovolemia, right ventricular failure, hypotension, and tachycardia. Finally, a management algorithm is proposed to identify a therapeutic intervention suited for the underlying physiologic condition causing suction.

Background Many new tools for abdominal aortic aneurysm (AAA) rupture risk evaluation have been developed. These new tools need detailed hemodynamic information in AAA. However, hemodynamic data obtained from in vivo research are lacking. Thus, the objective of this study was to analyze blood flow patterns in an in vivo AAA model to acquire real-time hemodynamic information using AneurysmFlow, a novel flow evaluation system. Methods Digital subtraction angiography images of patients who underwent endovascular aneurysm repair were analyzed using the visualization function of the AneurysmFlow to classify blood flow patterns as laminar or turbulent flow. The presence of boundary layer separation was also evaluated. The time taken for contrast medium to travel from the infrarenal aortic neck to aortic bifurcation was acquired to calculate the flow velocity. Associations between characteristics of aneurysm including lumen occupying ratio of intraluminal thrombus (ILT) and the hemodynamic flow pattern were evaluated. Results A total of 37 AAA patients was enrolled. Their blood flow patterns were evaluated using the AneurysmFlow. Logistic regression analyses with lumen occupying ratio of ILT as an independent variable showed that the larger the lumen occupying ratio of ILT, the more likely the aneurysm was to show a laminar pattern (P = 0.03) and the more likely the boundary layer separation would not exist (P = 0.04). The flow velocity from the infrarenal aortic neck to the aortic bifurcation showed a positive association with the lumen occupying ratio of the ILT in linear regression analysis (P Conclusion Hemodynamic analysis of AAA with the AneurysmFlow using real-time individual patient models showed different flow patterns and flow velocities depending on ILT. This novel analytic approach using AneurysmFlow has potential to play an important role in obtaining clinically meaningful hemodynamic information of AAA.

INTROUCTION Hemodynamic effects on the retrograde visceral reconstruction (RVR) for TAAA treatment by anastomotic angle remains unclear. This study aims to qualitatively and quantitatively investigate the effects of different anastomotic angles on hemodynamics and patency. METHODS Three RVR models with 45-degree, 60-degree and 90-degree anastomotic angles were reconstructed respectively by manipulating apostoperative patient-specific model. The manipulated models of the RVRs were numerically simulated and analyzed in terms of hemodynamics including theinstant and cumulative patency, flow pattern and indicators based on wall shear stress (WSS). RESULTS Although a smaller anastomotic angle may decrease the patency rate of common iliac arteries, it can improve the visceral perfusion during a cardiac cycle. More importantly, RVR with the smallest anastomotic angle experienced a minimal low time-averaged wall shear stress(TAWSS), high oscillatory shear index(OSI) and relative residence time(RRT) in the anastomosis region, whereas the largest anastomotic angle can introduce more unfavorable WSS in the graft trunk. Furthermore, a spiral flow pattern was observed in the proximal graft trunk of all three models, where no high-risk shear distribution was detected in this region. CONCLUSION A smaller anastomotic angle may have more benefits of hemodynamic environment in RVR, especially the wall shear stress (WSS) distribution and flow pattern in the graft trunk. We may also suggest that additional stents or an extended cuff for the graft can be used to induce spiral flow intentionally, which can further improve local hemodynamic environment and long-term prognosis.

Background Previous models describing the fetal-to-neonatal transition often lack oxygen saturation levels, homeostatic control mechanisms, phasic hemodynamic signals, or describe the heart with a time-varying elastance model. Methods We incorporated these elements in the adapted CircAdapt model with the one-fiber model for myocardial contraction, to simulate the hemodynamics of the healthy term human fetal circulation and its transition during the first 24 h after birth. The fetal-to-neonatal model was controlled by a time- and event-based script of changes occurring at birth, such as lung aeration and umbilical cord clamping. Model parameters were based on and validated with human and animal data. Results The fetal circulation showed low pulmonary blood flow, right ventricular dominance, and inverted mitral and tricuspid flow velocity patterns, as well as high mean ductus venosus flow velocity. The neonatal circulation showed oxygen saturation levels to gradually increase to 98% in the first 15 min after birth as well as temporary left ventricular volume overload. Conclusions Hemodynamics of the term fetus and 24-h-old neonate, as well as the events occurring directly after birth and the transition during the first 24 h after birth, were realistically represented, allowing the model to be used for educational purposes and future research. ImpactWith the addition of oxygen saturation levels, homeostatic pressure-flow control mechanisms, and the one-fiber model for myocardial contraction, a new closed-loop cardiovascular model was constructed to give more insight into the healthy term human fetal circulation and its cardiovascular transition during the first 24 h after birth. Extensive validation confirmed that the hemodynamics of the term fetus and the fetal-to-neonatal transition were realistically represented with the model. This well-validated and versatile model can serve as an education as well as a research platform for in silico investigation of fetal-to-neonatal hemodynamic changes under a wide range of physiological and pathophysiological conditions.

Supplemental Digital Content is available in the text., Background: Many heart diseases can result in reduced pumping capacity of the heart muscle. A mismatch between ATP demand and ATP production of cardiomyocytes is one of the possible causes. Assessment of the relation between myocardial ATP production (MVATP) and cardiac workload is important for better understanding disease development and choice of nutritional or pharmacologic treatment strategies. Because there is no method for measuring MVATP in vivo, the use of physiology-based metabolic models in conjunction with protein abundance data is an attractive approach. METHOD: We developed a comprehensive kinetic model of cardiac energy metabolism (CARDIOKIN1) that recapitulates numerous experimental findings on cardiac metabolism obtained with isolated cardiomyocytes, perfused animal hearts, and in vivo studies with humans. We used the model to assess the energy status of the left ventricle of healthy participants and patients with aortic stenosis and mitral valve insufficiency. Maximal enzyme activities were individually scaled by means of protein abundances in left ventricle tissue samples. The energy status of the left ventricle was quantified by the ATP consumption at rest (MVATP[rest]), at maximal workload (MVATP[max]), and by the myocardial ATP production reserve, representing the span between MVATP(rest) and MVATP(max). Results: Compared with controls, in both groups of patients, MVATP(rest) was increased and MVATP(max) was decreased, resulting in a decreased myocardial ATP production reserve, although all patients had preserved ejection fraction. The variance of the energetic status was high, ranging from decreased to normal values. In both patient groups, the energetic status was tightly associated with mechanic energy demand. A decrease of MVATP(max) was associated with a decrease of the cardiac output, indicating that cardiac functionality and energetic performance of the ventricle are closely coupled. Conclusions: Our analysis suggests that the ATP-producing capacity of the left ventricle of patients with valvular dysfunction is generally diminished and correlates positively with mechanical energy demand and cardiac output. However, large differences exist in the energetic state of the myocardium even in patients with similar clinical or image-based markers of hypertrophy and pump function. Registration: URL: https://www.clinicaltrials.gov; Unique identifiers: NCT03172338 and NCT04068740.

PURPOSE Leaflet thrombosis is a significant adverse event after transcatheter aortic valve (TAV) replacement (TAVR). The purpose of our study was to present a semi-empirical, mathematical model that links patient-specific anatomic, valve, and flow parameters to predict likelihood of leaflet thrombosis. METHODS The two main energy sources of neo-sinus (NS) washout after TAVR include the jet flow downstream of the TAV and NS geometric change in volume due to the leaflets opening and closing. Both are highly dependent on patient anatomic and hemodynamic factors. As rotation of blood flow is prevalent in both the sinus of Valsalva and then the NS, we adopted the vorticity flux or circulation (Г) as a metric quantifying overall washout. Leaflet thrombus volumes were segmented based on hypo-attenuating leaflet thickening (HALT) in post-TAVR patient's gated computed tomography. Г was assessed using dimensional scaling as well as computational fluid dynamics (CFD) respectively and correlated to the thrombosis volumes using sensitivity and specificity analysis. RESULTS Г in the NS, that accounted for patient flow and anatomic conditions derived from scaling arguments significantly better predicted the occurrence of leaflet thrombus than CFD derived measures such as stasis volumes or wall shear stress. Given results from the six patient datasets considered herein, a threshold Г value of 28.0 yielded a sensitivity and specificity of 100% where patients with Gamma < 28 developed valve thrombosis. A 10% error in measurements of all variables can bring the sensitivity specificity down to 87%. CONCLUSION A predictive model relating likelihood of valve thrombosis using Г in the NS was developed with promising sensitivity and specificity. With further studies and improvements, this predictive technology may lead to alerting physicians on the risk for thrombus formation following TAVR.

Cardiac voltage-gated sodium channel gain-of-function prolongs repolarization in the long-QT syndrome type 3 (LQT3). Previous studies suggest that narrowing the perinexus within the intercalated disc, leading to rapid sodium depletion, attenuates LQT3-associated action potential duration (APD) prolongation. However, it remains unknown whether extracellular sodium concentration modulates APD prolongation during sodium channel gain-of-function. We hypothesized that elevated extracellular sodium concentration and widened perinexus synergistically prolong APD in LQT3. LQT3 was induced with sea anemone toxin (ATXII) in Langendorff-perfused guinea pig hearts (n = 34). Sodium concentration was increased from 145 to 160 mM. Perinexal expansion was induced with mannitol or the sodium channel β1-subunit adhesion domain antagonist (βadp1). Epicardial ventricular action potentials were optically mapped. Individual and combined effects of varying clefts and sodium concentrations were simulated in a computational model. With ATXII, both mannitol and βadp1 significantly widened the perinexus and prolonged APD, respectively. The elevated sodium concentration alone significantly prolonged APD as well. Importantly, the combination of elevated sodium concentration and perinexal widening synergistically prolonged APD. Computational modeling results were consistent with animal experiments. Concurrently elevating extracellular sodium and increasing intercalated disc edema prolongs repolarization more than the individual interventions alone in LQT3. This synergistic effect suggests an important clinical implication that hypernatremia in the presence of cardiac edema can markedly increase LQT3-associated APD prolongation. Therefore, to our knowledge, this is the first study to provide evidence of a tractable and effective strategy to mitigate LQT3 phenotype by means of managing sodium levels and preventing cardiac edema in patients. NEW & NOTEWORTHY This is the first study to demonstrate that the long-QT syndrome type 3 (LQT3) phenotype can be exacerbated or concealed by regulating extracellular sodium concentrations and/or the intercalated disc separation. The animal experiments and computational modeling in the current study reveal a critically important clinical implication: sodium dysregulation in the presence of edema within the intercalated disc can markedly increase the risk of arrhythmia in LQT3. These findings strongly suggest that maintaining extracellular sodium within normal physiological limits may be an effective and inexpensive therapeutic option for patients with congenital or acquired sodium channel gain-of-function diseases.

The human aorta is a high-risk area for vascular diseases, which are commonly restored by thoracic endovascular aortic repair. In this paper, we report a promising shear-activated targeted nanoparticle drug delivery strategy to assist in the treatment of coarctation of the aorta and aortic aneurysm. Idealized three-dimensional geometric models of coarctation of the aorta and aortic aneurysm are designed, respectively. The unique hemodynamic environment of the diseased aorta is used to improve nanoparticle drug delivery. Micro-carriers with nanoparticle drugs would be targeting activated to release nanoparticle drugs by local abnormal shear stress rate (SSR). Coarctation of the aorta provides a high SSR hemodynamic environment, while the aortic aneurysm is exposed to low SSR. We propose a method to calculate the SSR thresholds for the diseased aorta. Results show that the upstream near-wall area of the diseased location is an ideal injection location for the micro-carriers, which could be activated by the abnormal SSR. Released nanoparticle drugs would be successfully targeted delivered to the aortic diseased wall. Besides, the high diffusivity of the micro-carriers and nanoparticle drugs has a significant impact on the surface drug concentrations of the diseased aortic walls, especially for aortic aneurysms. This study preliminary demonstrates the feasibility of shear-activated targeted nanoparticle drug delivery in the treatment of aortic diseases and provides a theoretical basis for developing the drug delivery system and novel therapy., 22 pages, 7 figures

Myocardial function and exercise reserve are important determinants of outcome in pulmonary arterial hypertension (PAH) but are incompletely understood. For this study, we performed subject-specific computer simulations, based on invasive measurements and cardiac magnetic resonance imaging (CMR), to investigate whole circulation properties in PAH at rest and exercise and determinants of exercise reserve. CMR and right heart catheterization were performed in nine patients with idiopathic PAH, and CMR in 10 healthy controls. CMR during exercise was performed in seven patients with PAH. A full-circulation computer model was developed, and model parameters were optimized at the individual level. Patient-specific simulations were used to analyze the effect of right ventricular (RV) inotropic reserve on exercise performance. Simulations achieved a high consistency with observed data. RV contractile force was increased in patients with PAH (127.1 ± 28.7 kPa vs. 70.5 ± 14.5 kPa

Aims Because evidence regarding risk stratification predicting prognosis of patients with heart failure (HF) decompensation attended in primary care is lacking, we developed and externally validated a model to forecast death/hospitalization during the first 30 days after an episode of decompensation. The predictive model is based on variables easily obtained in primary care settings. Methods and results HEFESTOS is a multinational study consisting of a derivation cohort of HF patients recruited in 14 primary healthcare centres in Barcelona and a validation cohort from primary healthcare in 9 other European countries. The derivation and validation cohorts included 561 and 250 patients, respectively. Percentages of women in the derivation and validation cohorts were 56.3% and 47.6% (P = 0.026), respectively. Mean age was 82.2 years (SD 8.03) in the derivation cohort, and 79.3 years (SD 10.3) in the validation one (P = 0.001). HF with preserved ejection fraction represented 72.1% in the derivation cohort and 58.8% in the validation one (P = 0.004). Mortality/hospitalization during the first 30 days after a decompensation episode was 30.5% and 26% (P = 0.225) for the derivation and validation cohorts, respectively. Multivariable logistic regression models were performed to develop a score of risk. The identified predictors were worsening of dyspnoea [odds ratio (OR): 2.5; P = 0.001], orthopnoea (OR: 2.16; P = 0.01), paroxysmal nocturnal dyspnoea (OR: 2.25; P = 0.01), crackles (OR: 2.35; P = 0.01), New York Heart Association functional class III/IV (OR: 2.11; P = 0.001), oxygen saturation ≤ 90% (OR: 4.98; P 100 b.p.m. (OR: 2.72; P = 0.002), and previous hospitalization due to HF (OR: 2.45; P 20%. Outcome incidence was 2.7% for the low‐risk group, 12.8% for medium risk, and 46.2% for high risk in the derivation cohort, and 9.1%, 12.9%, and 39.6% in the validation one. Conclusions The HEFESTOS score, based on variables easily accessible in a community setting and validated in an external European cohort, properly predicted the risk of death/hospitalization during the first 30 days after an HF decompensation episode.

Objective Aortic incompetence (AI) is observed to be accelerated in the continuous-flow left ventricular assist device (LVAD) population and is related to increased mortality. Using computational fluid dynamics (CFD), we investigated the hemodynamic conditions related to the orientation of the LVAD outflow in these patients. Method We identified 10 patients with new aortic regurgitation, and 20 who did not, after LVAD implantation between 2009 and 2018. Three-dimensional models of patients' aortas were created from their computed tomography scans. The geometry of the LVAD outflow graft in relation to the aorta was quantified using azimuth angles (AA), polar angles (PAs), and distance from aortic root. The models were used to run CFD simulations, which calculated the pressures and wall shear stress (rWSS) exerted on the aortic root. Results The AA and PA were found to be similar. However, for combinations of high values of AA and low values of PA, there were no patients with AI. The distance from aortic root to the outflow graft was also smaller in patients who developed AI (3.39 ± 0.7 vs 4.07 ± 0.77 cm, P = .04). There was no significant difference in aortic root pressures in the 2 groups. The rWSS was greater in AI patients (4.60 ± 5.70 vs 2.37 ± 1.20 dyne/cm2, P Conclusions Using CFD simulations, we demonstrated that patients who developed de novo AI have greater rWSS at the aortic root, and their outflow grafts were placed closer to the aortic roots than those patients without de novo AI.

Bicuspid aortic valve (BAV) is the most common congenital heart disease. Calcific aortic valve disease (CAVD) accounts for the majority of aortic stenosis (AS) cases. Half of the patients diagnosed with AS have a BAV, which has an accelerated progression rate. This study aims to develop a computational modeling approach of both the calcification progression in BAV, and its biomechanical response incorporating fluid-structure interaction (FSI) simulations during the disease progression. The calcification is patient-specifically reconstructed from Micro-CT images of excised calcified BAV leaflets, and processed with a novel reverse calcification technique that predicts prior states of CAVD using a density-based criterion, resulting in a multilayered calcified structure. Four progressive multilayered calcified BAV models were generated: healthy, mild, moderate, and severe, and were modeled by FSI simulations during the full cardiac cycle. A valve apparatus model, composed of the excised calcified BAV leaflets, was tested in an in-vitro pulse duplicator, to validate the severe model. The healthy model was validated against echocardiography scans. Progressive AS was characterized by higher systolic jet flow velocities (2.08, 2.3, 3.37, and 3.85 m s(−1)), which induced intense vortices surrounding the jet, coupled with irregular recirculation backflow patterns that elevated viscous shear stresses on the leaflets. This study shed light on the fluid-structure mechanism that drives CAVD progression in BAV patients.

Cardiovascular manifestations are frequent in COVID-19 infection and are predictive of adverse outcomes. Elevated cardiac biomarkers are common findings in patients with cardiovascular comorbidities and severe COVID-19 infection. Troponin, inflammatory and thrombotic markers may also improve risk prediction in COVID-19. In our comprehensive review, we provide an overview of the incidence, potential mechanisms and outcome of acute cardiac injury in COVID-19. Thereby, we discuss coagulation abnormalities in sepsis and altered immune response as contributing factors favoring myocardial injury. We further highlight the role of endothelial damage in the pathophysiological concepts. Finally, observational studies addressing the incidence of myocardial infarction during COVID-19 pandemic are discussed.

Coronary flow is affected by evolving events such as atherosclerotic plaque formation, rupture, and thrombosis, resulting in myocardial ischemia and infarction. Highly resolved 3D hemodynamic data at the stenosis is essential to model shear-sensitive thrombotic events in coronary artery disease. We developed a hybrid 1D–3D simulation framework to compute patient-specific coronary hemodynamics efficiently. A 1D model of the coronary flow is coupled to an image-based 3D model of the region of interest. This framework affords the advantages of reduced-order modeling, decreasing the global computational cost, without sacrificing the accuracy of the quantities of interest. We validated our 1D–3D model against full 3D coronary simulations in healthy and diseased conditions. Our results showed good agreement between the 3D and the 1D–3D models while reducing the computational cost by 40-fold compared to the 3D simulation. The 1D–3D model predicted left/right coronary flow distribution within 3% and provided an accurate estimation of fractional flow reserve and wall shear stress distribution at the stenosis comparable to the 3D simulation. Savings in computational cost may be significant in situations with changing geometry, such as growing thrombosis. Also, this approach would allow quantifying the time-dependent effect of thrombotic growth and occlusion on the global coronary circulation.

Heart failure is a major health risk, and with limited availability of donor organs, there is an increasing need for developing cardiac assist devices (CADs). Mock circulatory loops (MCL) are an important in-vitro test platform for CAD’s performance assessment and optimisation. The MCL is a lumped parameter model constructed out of hydraulic and mechanical components aiming to simulate the native cardiovascular system (CVS) as closely as possible. Further development merged MCLs and numerical circulatory models to improve flexibility and accuracy of the system; commonly known as hybrid MCLs. A total of 128 MCLs were identified in a literature research until 25 September 2020. It was found that the complexity of the MCLs rose over the years, recent MCLs are not only capable of mimicking the healthy and pathological conditions, but also implemented cerebral, renal and coronary circulations and autoregulatory responses. Moreover, the development of anatomical models made flow visualisation studies possible. Mechanical MCLs showed excellent controllability and repeatability, however, often the CVS was overly simplified or lacked autoregulatory responses. In numerical MCLs the CVS is represented with a higher order of lumped parameters compared to mechanical test rigs, however, complex physiological aspects are often simplified. In hybrid MCLs complex physiological aspects are implemented in the hydraulic part of the system, whilst the numerical model represents parts of the CVS that are too difficult to represent by mechanical components per se. This review aims to describe the advances, limitations and future directions of the three types of MCLs.

Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.

Thoracic endovascular aortic repair (TEVAR) has been accepted as the mainstream treatment for type B aortic dissection, but post-TEVAR biomechanical-related complications are still a major drawback. Unfortunately, the stent-graft (SG) configuration after implantation and biomechanical interactions between the SG and local aorta are usually unknown prior to a TEVAR procedure. The ability to obtain such information via personalised computational simulation would greatly assist clinicians in pre-surgical planning. In this study, a virtual SG deployment simulation framework was developed for the treatment for a complicated aortic dissection case. It incorporates patient-specific anatomical information based on pre-TEVAR CT angiographic images, details of the SG design and the mechanical properties of the stent wire, graft and dissected aorta. Hyperelastic material parameters for the aortic wall were determined based on uniaxial tensile testing performed on aortic tissue samples taken from type B aortic dissection patients. Pre-stress conditions of the aortic wall and the action of blood pressure were also accounted for. The simulated post-TEVAR configuration was compared with follow-up CT scans, demonstrating good agreement with mean deviations of 5.8% in local open area and 4.6 mm in stent strut position. Deployment of the SG increased the maximum principal stress by 24.30 kPa in the narrowed true lumen but reduced the stress by 31.38 kPa in the entry tear region where there was an aneurysmal expansion. Comparisons of simulation results with different levels of model complexity suggested that pre-stress of the aortic wall and blood pressure inside the SG should be included in order to accurately predict the deformation of the deployed SG.

Personalized models of cardiac electrophysiology (EP) that match clinical observation with high fidelity, referred to as cardiac digital twins (CDTs), show promise as a tool for tailoring cardiac precision therapies. Building CDTs of cardiac EP relies on the ability of models to replicate the ventricular activation sequence under a broad range of conditions. Of pivotal importance is the His–Purkinje system (HPS) within the ventricles. Workflows for the generation and incorporation of HPS models are needed for use in cardiac digital twinning pipelines that aim to minimize the misfit between model predictions and clinical data such as the 12 lead electrocardiogram (ECG). We thus develop an automated two stage approach for HPS personalization. A fascicular-based model is first introduced that modulates the endocardial Purkinje network. Only emergent features of sites of earliest activation within the ventricular myocardium and a fast-conducting sub-endocardial layer are accounted for. It is then replaced by a topologically realistic Purkinje-based representation of the HPS. Feasibility of the approach is demonstrated. Equivalence between both HPS model representations is investigated by comparing activation patterns and 12 lead ECGs under both sinus rhythm and right-ventricular apical pacing. Predominant ECG morphology is preserved by both HPS models under sinus conditions, but elucidates differences during pacing. Supplementary Information The online version contains supplementary material available at 10.1007/s10439-021-02825-9.

The intraventricular blood flow changed by blood pump flow dynamics may correlate with thrombosis and ventricular suction. The flow velocity, distribution of streamlines, vorticity, and standard deviation of velocity inside a left ventricle failing to different extents throughout the cardiac cycle when supported by an axial blood pump were measured by particle image velocimetry (PIV) in this study. The results show slower and static flow velocities existed in the central region of the left ventricle near the mitral valve and aortic valve and that were not sensitive to left ventricular (LV) failure degree or LV pressure. Strong vorticity located near the inner LV wall around the LV apex and the blood pump inlet was not sensitive to LV failure degree or LV pressure. Higher standard deviation of the blood velocity at the blood pump inlet decreased with increasing LV failure degree, whereas the standard deviation of the velocity near the atrium increased with increasing intraventricular pressure. The experimental results demonstrated that the risk of thrombosis inside the failing left ventricle is not related to heart failure degree. The “washout” performance of the strong vorticity near the inner LV wall could reduce the thrombotic potential inside the left ventricle and was not related to heart failure degree. The vorticity near the aortic valve was sensitive to LV failure degree but not to LV pressure. We concluded that the risk of blood damage caused by adverse flow inside the left ventricle decreased with increasing LV pressure.

COVID-19 bezeichnet eine der schlimmsten Krisen unserer Generation und stellt (nicht nur) für das Gesundheitssystem eine schwer bewältigbare Herausforderung dar. Mortalität und Morbidität sind im Vergleich zu anderen saisonalen Erkrankungen wie der Influenza deutlich erhöht. COVID-19 bedroht allerdings nicht die gesamte Bevölkerung in gleichem Maße. Hochrisikopatienten sind älter und leiden an kardiovaskulären Erkrankungen wie Bluthochdruck, Diabetes mellitus oder einer koronaren Herzerkrankung. Um das Risiko für einen schweren Erkrankungsverlaufs zu quantifizieren bedarf es einer multimodalen Herangehensweise. Verschiedene Risikostratifizierungssysteme stehen zu Verfügung um ungünstige Verläufe wie Intensivbehandlung oder Gesamtmortalität vorauszusagen. Biomarker wie Troponin-I, D-Dimere und NT pro-BNP kombiniert mit echokardiographischen Parametern wie links- und rechtsventrikulärer Pumpfunktion sowie pulmonalarteriellem Druck können hilfreich sein um Hochrisikopatienten zu identifizieren, die ein intensiviertes Monitoring und eine stringentere Behandlung benötigen. Da kardiovaskuläre Risikofaktoren und Komorbiditäten von großer Bedeutung zur Abschätzung des Verlaufs einer SARS-CoV-2 Infektion sind, könnten alle hospitalisierten COVID-19 Patienten von einer routinemäßigen kardiologischen Betreuung durch ein COVID-19-Heart-Team profitieren. Ein frühzeitiges Erkennen von (kardiovaskulären) Hochrisikopatienten könnte das Management erleichtern sowie die Prognose einer schweren SARS-CoV-2 Infektion verbessern.The COVID-19 pandemic represents one of the largest burdens of our generation with enormous consequences on socioeconomics and healthcare around the globe. COVID-19 is associated with increased morbidity and mortality. However, the course of disease differs among the population. Advanced age and a compromised immune system among others are associated with significantly worse outcomes. Furthermore, cardiovascular comorbidities are of utmost importance and may influence severity of SARS-CoV-2 infection to a considerable degree. This review aims to identify risk factors, biomarkers, and echocardiographic parameters associated with a severe course of SARS-CoV-2 infection and give a brief overview of the management of COVID-19 patients with cardiovascular comorbidities. Furthermore, available models identifying high-risk patients and their performance in prediction of severe outcomes of SARS-CoV-2 are addressed.

Background Heart failure (HF) is a common disease globally. Ventricular assist devices (VADs) are widely used to treat HF. In contrast to the natural heart, different VADs generate different blood flow waves in the aorta. Objective To explore whether the different inflow rate waveforms from the ascending aorta generate far-reaching hemodynamic influences on the human aortic arch. Methods An aortic geometric model was reconstructed based on computed tomography data of a patient with HF. A total of five numerical simulations were conducted, including a case with the inflow rate waveforms from the ascending aorta with normal physiological conditions, two HF, and two with typical VAD support. The hemodynamic parameters, wall shear stress (WSS), oscillatory shear index (OSI), relative residence time (RRT), and the strength of the helical flow, were calculated. Results In contrast to the natural heart, numerical simulation showed HF decreasing WSS and inducing higher OSI and RRT. Moreover, HF weakened helical flow strength. Pulsatile flow VADs will elevated the WSS, inducing some helical flow, while continuous flow VAD could not produce any helical flow. Conclusions HF leads to an adverse hemodynamic environment by decreasing WSS and reducing the helical flow strength. Pulsatile flow VADs are more advantageous than the continuous flow VADs on hemodynamic effects. Thus, pulsatile flow VADs may be a better option for patients with HF.

Background Despite the great evolution of endograft devices for thoracic endovascular aortic repair (TEVAR), threatening related complication such as graft migration and endoleaks still occur during follow up. The Drag Forces (DF), that is the displacement forces that play a role in graft migration and endoleaks caused by the blood flow against the thoracic graft, can be studied by means of Computational Fluid Dynamics (CFD). Method A general review of papers found in current literature was performed. CFD studies available on the topic of thoracic aortic diseases and DF were analyzed. All anatomic, hemodynamics or graft related factors which could have an impact on DF were reported. Results Different factors deeply influence DF magnitude in the different site of the Ishimaru's zones classification: angulation, tortuosity and length of the landing zone, graft diameter, length and deployment position, blood pressure, pulse waveform, blood viscosity and patient heart rate have been related to the magnitude of DF. Moreover, also the three-dimensional orientation of DF is emerging as a fundamental issue from CFD studies. DF can be divided in sideways and upward components. The former, even of higher magnitude in zone 0, maintain always an orthogonal orientation and does not change in any type of aortic arch; the latter result strictly related to the anatomic complexity of the aortic arch with values up to four times higher in zone 3. Conclusion Different DF magnitude and orientation could explain how TEVAR have higher rate of migration and endoleaks when we face with more complex aortic anatomies. All these aspects should be foreseen during the planning of TEVAR procedure. In this field, collaboration between physicians and engineers is crucial, as both parts have a primary role in understanding and describing hidden aspects involved in TEVAR procedures.

Chronic thromboembolic pulmonary hypertension (CTEPH) is caused by recurrent or unresolved pulmonary thromboemboli, leading to perfusion defects and increased arterial wave reflections. CTEPH treatment aims to reduce pulmonary arterial pressure and reestablish adequate lung perfusion, yet patients with distal lesions are inoperable by standard surgical intervention. Instead, these patients undergo balloon pulmonary angioplasty (BPA), a multi-session, minimally invasive surgery that disrupts the thromboembolic material within the vessel lumen using a catheter balloon. However, there still lacks an integrative, holistic tool for identifying optimal target lesions for treatment. To address this insufficiency, we simulate CTEPH hemodynamics and BPA therapy using a multiscale fluid dynamics model. The large pulmonary arterial geometry is derived from a computed tomography (CT) image, whereas a fractal tree represents the small vessels. We model ring- and web-like lesions, common in CTEPH, and simulate normotensive conditions and four CTEPH disease scenarios; the latter includes both large artery lesions and vascular remodeling. BPA therapy is simulated by simultaneously reducing lesion severity in three locations. Our predictions mimic severe CTEPH, manifested by an increase in mean proximal pulmonary arterial pressure above 20 mmHg and prominent wave reflections. Both flow and pressure decrease in vessels distal to the lesions and increase in unobstructed vascular regions. We use the main pulmonary artery (MPA) pressure, a wave reflection index, and a measure of flow heterogeneity to select optimal target lesions for BPA. In summary, this study provides a multiscale, image-to-hemodynamics pipeline for BPA therapy planning for inoperable CTEPH patients., Comment: 41 pages, 9 figures, 4 tables

Introduction Thanks to shared circulation in monochorionic twins, single intrauterine fetal death (IUD) may lead to acute feto-fetal transfusion (aFFTR). The objective of the study was to describe our model of aFFTR simulation after IUD in monochorionic (MC) twins. Methods Prospective study analyzed 99 fresh MC placentas with the physiological course. A specially designed protocol was used for the preparation and analysis of the placentas. A pair of infusion sets fixed together using a mechanical mercury sphygmomanometer cuff was connected to the cannulated umbilical arteries. The tonometer was pressurized up to 30 and 40 mmHg. A positive finding of aFFTR was determined as the amount exceeding 1 ml of dye flowed out of the umbilical cord simulating a dead fetus. The number and types of anastomoses, types, and distances between cords insertions, and the size of the placental areas for each fetus were also statistically analyzed. The placental angioarchitecture with and without proven aFFTR was statistically compared, odds ratio (OR) and multivariable logistic analysis were performed. Results A total of 49/99 (49.5%) cases of aFFTR was proven, and the average transfusion time of 1 ml was 30 s (19–46 s). aFFTR was present in 49/78 (62.8%) of placentas with arterio-arterial (AA) anastomosis. The median diameter of AA anastomoses with the present, and absent aFFTRF was 2.0 mm and 1.0 mm, respectively. The proven interfetal transfusion was 8%, 31%, and 61% in AA anastomoses with a diameter below 0,5 mm, 0,5–1,5 mm, and above 1,5 mm, respectively (p 1.5 mm had OR of 44.2 (95% CI 5.54–352.39). In the case of coexistence of AA anastomosis and umbilical cord distance ≤5th percentile, the aFFTRF occurred in 90.9%. Discussion The potential risk of aFFTR in monochorionic twins is mainly due to the presence and nature of AA anastomoses. The diameter and length of the vessels play a crucial role, which is clinically related to the distance of the umbilical cords insertions.

Artificial organs 46(1), 71-82 (2022). doi:10.1111/aor.14042, Published by Wiley-Blackwell, Oxford [u.a.]

Improper left ventricular assist device (LVAD) inflow cannula (IC) positioning creates areas of stasis and low pulsatility that predispose thromboembolism, but may be mitigated with LVAD speed modulation. A mock loop study was performed to assess the sensitivity of left ventricle (LV) flow architecture to IC position and speed modulation during HeartMate3 support. System pressure, flow, and the time-resolved velocity field were measured within a transparent silicone LV for three IC angles and three IC insertion depths at matched levels of cardiac function and LVAD speed. Inflow cannula angulation towards the septum increased the resistance to LVAD flow as well as increasing the size and energy of the counter-clockwise (CCW) vortex. Apical velocity was reduced compared to IC angulation towards the mitral valve, but regional pulsatility was maintained across all angles and LVAD speeds. Increased IC protrusion decreased LVAD flow resistance, increasing velocity within the IC but reducing flow and pulsatility in the adjacent apical region. Increasing LVAD flow resistance improves aortic valve opening and strengthens the CCW vortex which directs inflow towards the septum, producing higher blood residence time and shear activation potential. Despite this impact on flow architecture, pulsatility reduction with increased LVAD speed was minimal with the HeartMate3 speed modulation feature.

Clinical observation and ex vivo studies have established a strong association between inflammation and postoperative atrial fibrillation (POAF). However, it is unclear whether the inflammatory phenotype is causally linked to this event or is an epiphenomenon, and it is not known which inflammatory meditators may increase susceptibility to POAF. The limitations of available animal models of spontaneous POAF (sPOAF) makes it difficult to select an experimental system. Here, we provide experimental and clinical evidence for mechanistic involvement of interleukin-6 (IL-6) in sPOAF. PHASE I: We established a mouse model of cardiac surgery with nonpaced sPOAF. IL-6 knockout mice were protected from sPOAF compared with wild-type mice. PHASE II: At 48 hours after surgery, the heart was separated into 6 regions and cultured. IL-6 was expressed in all regions, with highest abundance in the left atrium (LA). In PHASE III, we demonstrated that IL-6 in the LA elicited early profibrotic properties in atria via the pSTAT3/STAT3 signaling pathway and contributed to sPOAF. PHASE IV: In a translational prospective clinical study, we demonstrated that humans with POAF had a higher IL-6 concentration in pericardial drainage (PD). This study provides preliminary evidence of a causal relationship between IL-6 and POAF in a novel nonpaced sPOAF mouse model. IL-6 is a crucial prerequisite for eliciting profibrotic properties in cardiac myocytes via the pSTAT3 pathway during the early postoperative period, leading to an increased susceptibility to POAF. Measuring IL-6 in PD could be a new noninvasive biomarker for the clinical prediction of POAF.

Aims: Wall shear stress (WSS) has been considered a major determinant of aortic atherosclerosis. Recently, non-obstructive general angioscopy (NOGA) was developed to visualize various atherosclerotic pathologies, including in vivo ruptured plaque (RP) in the aorta. However, the relationship between aortic RP and WSS distribution within the aortic wall is unclear. This study aimed to investigate the relationship between aortic NOGA-derived RP and the stereographic distribution of WSS by computational fluid dynamics (CFD) modeling using three-dimensional computed tomography (3D-CT) angiography. Methods: We investigated 45 consecutive patients who underwent 3D-CT before coronary angiography and NOGA during coronary angiography. WSS in the aortic arch was measured by CFD analysis based on the finite element method using uniform inlet and outlet flow conditions. Aortic RP was detected by NOGA. Results: Patients with a distinct RP showed a significantly higher maximum WSS value in the aortic arch than those without aortic RP (56.2±30.6 Pa vs 36.2±19.8 Pa, p =0.017), no significant difference was noted in the mean WSS between those with and without aortic RP. In a multivariate logistic regression analysis, the presence of a maximum WSS value more than a specific value was a significant predictor of aortic RP (odds ratio 7.21, 95% confidence interval 1.78-37.1, p =0.005). Conclusions: Aortic RP detected by NOGA was strongly associated with a higher maximum WSS in the aortic arch derived by CFD using 3D-CT. The maximum WSS value may have an important role in the underlying mechanism of not only aortic atherosclerosis, but also aortic RP.

Type A aortic dissection (TAAD) involves the ascending aorta or the arch. Acute TAAD usually requires urgent replacement of the ascending aorta. However, a subset of these patients develops aortic rupture due to further dilatation of the residual dissected aorta. There is currently no reliable means to predict the risk of dilatation following TAAD repair. In this study, we performed a comprehensive morphological and hemodynamic analysis for patients with and without progressive aortic dilatation following surgical replacement of the ascending aorta. Patient-specific models of repaired TAAD were reconstructed from post-surgery computed tomography images for detailed computational fluid dynamic analysis. Geometric and hemodynamic parameters were evaluated and compared between patients with stable aortic diameters (N = 9) and those with aortic dilatation (N = 8). Our results showed that the number of re-entry tears and true/false lumen pressure difference were significantly different between the two groups. Patients with progressive aortic dilatation had higher luminal pressure difference (6.7 [4.6, 10.9] vs. 0.9 [0.5, 2.3] mmHg; P = 0.001) and fewer re-entry tears (1.5 [1, 2.8] vs. 5 [3.3, 7.5]; P = 0.02) compared to patients with stable aortic diameters, suggesting that these factors may serve as potential predictors of aneurysmal dilatation following surgical repair of TAAD.

Dilatació aòrtica; Aneurisma de l'aorta ascendent; Anàlisi espai-temporal Dilatación aórtica; Aneurisma de la aorta ascendente; Análisis espacio-temporal Aortic dilation; Ascending aorta aneurysm; Spatiotemporal analysis Motivated by the evidence that the onset and progression of the aneurysm of the ascending aorta (AAo) is intertwined with an adverse hemodynamic environment, the present study characterized in vivo the hemodynamic spatiotemporal complexity and organization in human aortas, with and without dilated AAo, exploring the relations with clinically relevant hemodynamic and geometric parameters. The Complex Networks (CNs) theory was applied for the first time to 4D flow magnetic resonance imaging (MRI) velocity data of ten patients, five of them presenting with AAo dilation. The time-histories along the cardiac cycle of velocity-based quantities were used to build correlation-based CNs. The CNs approach succeeded in capturing large-scale coherent flow features, delimiting flow separation and recirculation regions. CNs metrics highlighted that an increasing AAo dilation (expressed in terms of the ratio between the maximum AAo and aortic root diameter) disrupts the correlation in forward flow reducing the correlation persistence length, while preserving the spatiotemporal homogeneity of secondary flows. The application of CNs to in vivo 4D MRI data holds promise for a mechanistic understanding of the spatiotemporal complexity and organization of aortic flows, opening possibilities for the integration of in vivo quantitative hemodynamic information into risk stratification and classification criteria. Open access funding provided by Politecnico di Torino within the CRUI-CARE Agreement.

High wall shear stress (WSS) and near-infrared spectroscopy (NIRS) detected lipid-rich plaque (LRP) are both known to be associated with plaque destabilization and future adverse cardiovascular events. However, knowledge of spatial co-localization of LRP and high WSS is lacking. This study investigated the co-localization of LRP based on NIRS and high WSS. Fifty-three patients presenting acute coronary syndrome underwent NIRS-intravascular-ultrasound (NIRS-IVUS) imaging of a non-culprit coronary artery. WSS was obtained using WSS profiling in 3D-reconstructions of the coronary arteries based on fusion of IVUS-segmented lumen and CT-derived 3D-centerline. Thirty-eight vessels were available for final analysis and divided into 0.5 mm/45° sectors. LRP sectors, as identified by NIRS, were more often colocalized with high WSS than sectors without LRP. Moreover, there was a dose-dependent relationship between lipid content and high WSS exposure. This study is a first step in understanding the evolution of LRPs to vulnerable plaques. Graphical Abstract Electronic supplementary material The online version of this article (10.1007/s12265-020-10072-x) contains supplementary material, which is available to authorized users.

Acute Ischemic Stroke (AIS) is defined as the acute condition of occlusion of a cerebral artery and is often caused by a Hypertensive Condition (HC). Due to its sudden occurrence, AIS is not observable the right moment it occurs, thus information about instantaneous changes in hemodynamics is limited. This study aimed to propose an integrated Lumped Parameter (LP) model of the cardiovascular system to simulate an AIS and describe instantaneous changes in hemodynamics. In the integrated LP model of the cardiovascular system, heart chambers have been modelled with elastance systems with controlled pressure inputs; heart valves have been modelled with static open/closed pressure-controlled valves; eventually, the vasculature has been modelled with resistor-inductor-capacitor (RLC) direct circuits and have been linked to the rest of the system through a series connection. After simulating physiological conditions, HC has been simulated by changing pressure inputs and constant RLC parameters. Then, AIS occurring in arteries of different sizes have been simulated by considering time-dependent RLC parameters due to the elimination from the model of the occluding artery; instantaneous changes in hemodynamics have been evaluated by Systemic Arteriolar Flow (Qa) and Systemic Arteriolar Pressure (Pa) drop with respect to those measured in HC. Occlusion of arteries of different sizes leaded to an average Qa drop of 0.38 ml/s per cardiac cycle (with minimum and maximum values of 0.04 ml/s and 1.93 ml/s) and average Pa drop of 0.39 mmHg, (with minimum and maximum values of 0.04 mmHg and 1.98 mmHg). In conclusion, hemodynamic variations due to AIS are very small with respect to HC. A direct relation between the inverse of the length of the artery in which the occlusion occurs and the hemodynamic variations has been highlighted; this may allow to link the severity of AIS to the length of the interested artery.

Background: Besides biological factors, abdominal aortic aneurysm rupture is also caused by mechanical parameters, which are constantly affecting the wall’s tissue due to their abnormal values. The ability to evaluate these parameters could vastly improve the clinical treatment of patients with abdominal aortic aneurysms. The objective of this study was to develop and demonstrate a methodology to analyze the fluid dynamics that cause the wall stress distribution in abdominal aortic aneurysms, using accurate 3D geometry and a realistic, nonlinear, elastic biomechanical model using a computer-aided software. Methods: The geometry of the abdominal aortic aneurysm; was constructed on a 3D scale using computer-aided software SolidWorks (Dassault Systems SolidWorksCorp., Waltham MA). Due to the complex nature of the abdominal aortic aneurysm geometry, the physiological forces and constraints acting on the abdominal aortic aneurysm wall were measured by using a simulation setup using boundary conditions and initial conditions for different studies such as finite element analysis or computational fluid dynamics. Results: The flow pattern showed an increase velocity at the angular neck, followed by a stagnated flow inside the aneurysm sack. Furthermore, the wall shear stress analysis showed to focalized points of higher stress, the top and bottom of the aneurysm sack, where the flow collides against the wall. An increase of the viscosity showed no significant velocity changed but results in a slight increase in overall pressure and wall shear stress. Conclusions: Conducting computational fluid dynamics modeling of the abdominal aortic aneurysm using computer-aided software SolidWorks (Dassault Systems SolidWorksCorp., Waltham MA) proves to be an insightful approach for the clinical setting. The careful consideration of the biomechanics of the abdominal aortic aneurysm may lead to an improved, case-specific prediction of the abdominal aortic aneurysm rupture potential, which could significantly improve the clinical management of these patients.

A central shunt (CS) was an important surgery of systemic-to-pulmonary shunt (SPS) for the treatment of complex congenital heart diseases with decreased pulmonary blood flow (CCHDs-DPBF). There was no clear conclusion on how to deal with unclosed patent ductus arteriosus (PDA) during CS surgery. This study expanded the knowledge base on PDA by exploring the influence of the closing process of the PDA on the hemodynamic parameters for the CS model. The initial three-dimensional (3D) geometry was reconstructed based on the patient’s computed tomography (CT) data. Then, a CS configuration with three typical pulmonary artery (PA) dysplasia structures and different sizes of PDA was established. The three-element windkessel (3WK) multiscale coupling model was used to define boundary conditions for transient simulation through computational fluid dynamics (CFD). The results showed that the larger size of PDA led to a greater systemic-to-pulmonary shunt ratio ( Q S / A ), and the flow ratio of the left pulmonary artery (LPA) to right pulmonary artery (RPA) ( Q L / R ) was more close to 1, while both the proportion of high wall shear stress (WSS) areas and power loss decreased. The case of PDA nonclosure demonstrates that the aortic oxygen saturation (Sao2) increased, while the systemic oxygen delivery (Do2) decreased. In general, for the CS model with three typical PA dysplasia, the changing trends of hemodynamic parameters during the spontaneous closing process of PDA were roughly identical, and nonclosure of PDA had a series of hemodynamic advantages, but a larger PDA may cause excessive PA perfusion and was not conducive to reducing cyanosis symptoms.

Despite significant advances in its detection, understanding and management, atrial fibrillation (AF) remains a highly prevalent cardiac arrhythmia with a major impact on morbidity and mortality of millions of patients. AF results from complex, dynamic interactions between risk factors and comorbidities that induce diverse atrial remodelling processes. Atrial remodelling increases AF vulnerability and persistence, while promoting disease progression. The variability in presentation and wide range of mechanisms involved in initiation, maintenance and progression of AF, as well as its associated adverse outcomes, make the early identification of causal factors modifiable with therapeutic interventions challenging, likely contributing to suboptimal efficacy of current AF management. Computational modelling facilitates the multilevel integration of multiple datasets and offers new opportunities for mechanistic understanding, risk prediction and personalized therapy. Mathematical simulations of cardiac electrophysiology have been around for 60 years and are being increasingly used to improve our understanding of AF mechanisms and guide AF therapy. This narrative review focuses on the emerging and future applications of computational modelling in AF management. We summarize clinical challenges that may benefit from computational modelling, provide an overview of the different in silico approaches that are available together with their notable achievements, and discuss the major limitations that hinder the routine clinical application of these approaches. Finally, future perspectives are addressed. With the rapid progress in electronic technologies including computing, clinical applications of computational modelling are advancing rapidly. We expect that their application will progressively increase in prominence, especially if their added value can be demonstrated in clinical trials., Graphical Abstract

Successful implementation of automated blood sampling (ABS) into a telemetry instrumented canine cardiovascular model provides simultaneous cardiovascular assessment of novel compounds while collecting multiple blood samples for analysis of drug level, cytokines, and biomarkers. Purpose-bred male Beagle dogs (n = 36) were instrumented with a dual-pressure telemetry transmitter and vascular access port. Modifications to acclimation practices, surgical procedures, and housing were required for implementation of ABS in our established cardiovascular canine telemetry colony. These modifications have increased the use and reproducibility of the model by combining early pharmacokinetic and cardiovascular studies, thus achieving both refinement and reduction from a 3R perspective. In addition, the modified model can shorten timelines and reduce the compound requirement in early stages of drug development. This telemetry–ABS model provides an efficient means to quickly identify potential effects on key cardiovascular parameters in a large animal species and to obtain a more complete pharmacokinetic–pharmacodynamic profile for discovery compounds.

This paper aims to investigate detailed mechanical interactions between the pulmonary haemodynamics and left heart function in pathophysiological situations (e.g. atrial fibrillation and acute mitral regurgitation). This is achieved by developing a complex computational framework for a coupled pulmonary circulation, left atrium and mitral valve model. The left atrium and mitral valve are modelled with physiologically realistic three-dimensional geometries, fibre-reinforced hyperelastic materials and fluid–structure interaction, and the pulmonary vessels are modelled as one-dimensional network ended with structured trees, with specified vessel geometries and wall material properties. This new coupled model reveals some interesting results which could be of diagnostic values. For example, the wave propagation through the pulmonary vasculature can lead to different arrival times for the second systolic flow wave (S2 wave) among the pulmonary veins, forming vortex rings inside the left atrium. In the case of acute mitral regurgitation, the left atrium experiences an increased energy dissipation and pressure elevation. The pulmonary veins can experience increased wave intensities, reversal flow during systole and increased early-diastolic flow wave (D wave), which in turn causes an additional flow wave across the mitral valve (L wave), as well as a reversal flow at the left atrial appendage orifice. In the case of atrial fibrillation, we show that the loss of active contraction is associated with a slower flow inside the left atrial appendage and disappearances of the late-diastole atrial reversal wave (AR wave) and the first systolic wave (S1 wave) in pulmonary veins. The haemodynamic changes along the pulmonary vessel trees on different scales from microscopic vessels to the main pulmonary artery can all be captured in this model. The work promises a potential in quantifying disease progression and medical treatments of various pulmonary diseases such as the pulmonary hypertension due to a left heart dysfunction.

The achievement of clinically viable methodologies to simulate the hemodynamics in patient-specific coronary arteries is still a major challenge. Therefore, the novelty of this work is attained by the introduction of the viscoelastic property of blood in the numerical simulations, to study the role of the left coronary artery (LCA) geometry configuration in the atherosusceptibility. Apparently healthy patients were used and four different methodologies were tested. The methodology giving the most accurate results at the same time of having the lowest computational time is the one considering the viscoelastic property of blood and computational fluid dynamics. A Pearson correlation analysis was used to highlight relationships between geometric configuration and hemodynamic descriptors based on the simulated wall shear stress (WSS). The left main stem (LMS) has the greatest atherosusceptibility followed by the left anterior descending artery (LAD) since the relative residence time (RRT) average values are 3.81 and 3.70 Pa

The calf muscle pump is a major determinate of venous return in the legs but has not been studied as a risk factor for venous thromboembolism (VTE). A population-based cohort study of Olmsted County, Minnesota residents was performed using calf pump function (CPF) measurements from venous plethysmography studies from 1998 to 2015. Patients with a history of VTE were excluded. Nursing validated VTE outcomes from the Rochester Epidemiology Project were identified after the index study date, and patients with reduced CPF (rCPF) were compared with patients with normal CPF. A total of 1532 patients with recorded CPF (28% air and 72% strain gauge plethysmography) were included; 591 (38.5%) had normal CPF, 353 (23.0%) had unilateral rCPF, and 588 (38.3%) had bilateral rCPF. Any VTE occurred in 87 patients (5.7%) after a median follow-up of 11.7 years (range, 0-22.0 years). Comparing patients with bilateral reduced to bilateral normal CPF, the unadjusted hazard ratio (HR) for incident VTE was 2.0 (95% confidence interval [CI], 1.2-3.4) and after adjusting for age, BMI, and Charlson Comorbidity Index, the HR was 1.68 (95% CI, 0.98-2.89). The adjusted HR for ipsilateral deep vein thrombosis was evaluated in 3064 legs comparing legs with reduced to normal CPF and was 1.71 (95% CI, 1.03-2.84). Mortality was significantly higher in both the bilateral (P < .001) and unilateral (P < .001) rCPF groups compared with normal CPF. Our results demonstrate that CPF is a risk factor for VTE in an otherwise low-risk ambulatory population and might be a useful component in risk stratification models.

Aims Arrhythmogenic cardiomyopathy (AC) is an inherited cardiac disease, characterized by life-threatening ventricular arrhythmias and progressive cardiac dysfunction. The aim of this study is to use computer simulations to non-invasively estimate the individual patient’s myocardial tissue substrates underlying regional right ventricular (RV) deformation abnormalities in a cohort of AC mutation carriers. Methods and results In 68 AC mutation carriers and 20 control subjects, regional longitudinal deformation patterns of the RV free wall (RVfw), interventricular septum (IVS), and left ventricular free wall (LVfw) were obtained using speckle-tracking echocardiography. We developed and used a patient-specific parameter estimation protocol based on the multi-scale CircAdapt cardiovascular system model to create virtual AC subjects. Using the individual’s deformation data as model input, this protocol automatically estimated regional RVfw and global IVS and LVfw tissue properties. The computational model was able to reproduce clinically measured regional deformation patterns for all subjects, with highly reproducible parameter estimations. Simulations revealed that regional RVfw heterogeneity of both contractile function and compliance were increased in subjects with clinically advanced disease compared to mutation carriers without clinically established disease (17 ± 13% vs. 8 ± 4%, P = 0.01 and 18 ± 11% vs. 10 ± 7%, P Conclusion Regional RV deformation abnormalities in AC mutation carriers were related to reduced regional contractile function and tissue compliance. In clinically advanced disease stages, a characteristic apex-to-base heterogeneity of tissue abnormalities was present in the majority of the subjects, with most pronounced disease in the basal region of the RVfw.

Atrial fibrillation (AF) can alter intra-cardiac flow and cardiac output that subsequently affects aortic flow circulation. These changes may become more significant where they occur concomitantly with ageing. Aortic ageing is accompanied with morphological changes such as dilation, lengthening, and arch unfolding. While the recognition of AF mechanism has been the subject of numerous studies, less focus has been devoted to the aortic circulation during the AF and there is a lack of such investigation at different ages. The current work aims to address the present gap. First, we analyse aortic flow distribution in three configurations, which attribute to young, middle and old people, using geometries constructed via clinical data. We then introduce two transient inlet flow conditions representative of key AF-associated defects. Results demonstrate that both AF and ageing negatively affect flow circulation. The main consequence of concomitant occurrence is enhancement of endothelial cell activation potential (ECAP) throughout the vascular domain, mainly at aortic arch and descending thoracic aorta, which is consistent with some clinical observations. The outcome of the current study suggests that AF exacerbates the vascular defects occurred due to the ageing, which increases the possibility of cardiovascular diseases per se. Supplementary Information The online version contains supplementary material available at (10.1007/s10439-021-02744-9).

Mechanical dyssynchrony (MD) affects left ventricular (LV) mechanics and coronary perfusion. To understand the multifactorial effects of MD, we developed a computational model that bidirectionally couples the systemic circulation with the LV and coronary perfusion with flow regulation. In the model, coronary flow in the left anterior descending (LAD) and left circumflex (LCX) arteries affects the corresponding regional contractility based on a prescribed linear LV contractility-coronary flow relationship. The model is calibrated with experimental measurements of LV pressure and volume, as well as LAD and LCX flow rate waveforms acquired under regulated and fully dilated conditions from a swine under right atrial (RA) pacing. The calibrated model is applied to simulate MD. The model can simultaneously reproduce the reduction in mean LV pressure (39.3%), regulated flow (LAD: 7.9%; LCX: 1.9%), LAD passive flow (21.6%), and increase in LCX passive flow (15.9%). These changes are associated with right ventricular pacing compared with RA pacing measured in the same swine only when LV contractility is affected by flow alterations with a slope of 1.4 mmHg/mL(2) in a contractility-flow relationship. In sensitivity analyses, the model predicts that coronary flow reserve (CFR) decreases and increases in the LAD and LCX with increasing delay in LV free wall contraction. These findings suggest that asynchronous activation associated with MD impacts 1) the loading conditions that further affect the coronary flow, which may explain some of the changes in CFR, and 2) the coronary flow that reduces global contractility, which contributes to the reduction in LV pressure. NEW & NOTEWORTHY A computational model that couples the systemic circulation of the left ventricular (LV) and coronary perfusion with flow regulation is developed to study the effects of mechanical dyssynchrony. The delayed contraction in the LV free wall with respect to the septum has a significant effect on LV function and coronary flow reserve.

OBJECTIVES The actual underlying mechanisms of acute type A aortic dissection (AAAD) are not well understood. The present study aimed to elucidate the mechanism of AAAD using computational fluid dynamics (CFD) analysis. METHODS We performed CFD analysis using patient-specific computed tomography imaging in 3 healthy control cases and 3 patients with AAAD. From computed tomography images, we made a healthy control model or pre-dissection model for CFD analysis. Pulsatile cardiac flow during one cardiac cycle was simulated, and a three-dimensional flow streamline was visualized to evaluate flow velocity, wall shear stress and oscillatory shear index (OSI). RESULTS In healthy controls, the transvalvular aortic flow was parallel to the ascending aorta. There was no spotty high OSI area at the ascending aorta. In pre-dissection patients, accelerated transvalvular aortic flow was towards the posterolateral ascending aorta. The vortex flow was observed on the side of the lesser curvature in mid-systole and expanded throughout the entire ascending aorta during diastole. Systolic wall shear stress was high due to the accelerated aortic blood flow on the side of the greater curvature of the ascending aorta. On the side of the lesser curvature, high OSI areas were observed around the vortex flow. In all pre-dissection cases, a spotty high OSI area was in close proximity to the actual primary entry site of the future AAAD. CONCLUSIONS The pre-onset high OSI area with vortex flow is closely associated with the future primary entry site. Therefore, we can elucidate the mechanism of AAAD with CFD analysis.

Despite numerous experimental observations regarding heart failure with preserved ejection fraction (HFpEF), which is characterized mainly by left ventricular hypertrophy and a left ventricular ejection fraction over 50%, myocardial dynamics under HFpEF have not yet been fully clarified, particularly regarding the relationship between myocardial strain distribution and myocardial work. To address this issue, we numerically investigated radial distribution of myocardial strain during a cardiac cycle with fixed internal volume at the end of the systolic and diastolic phases under different mechanical conditions, such as those involving myocardial thickness and elasticity of myocardial fibers. The myocardium was a modeled as a visco-hyperelastic continuous material. This model was taken into account that active contractile stress along the myocardial fiber direction depends on membrane potential change. Our numerical results showed that both radial and circumferential strains decreased as wall thickness increased, which reflected cardiac hypertrophy, but that myocardial work became larger than that observed with thin ventricular walls. Further, the change in left ventricular diastolic internal pressure caused circumferential strain, while fiber stiffness contributed to radial strain. Since peak circumferential strain was well estimated by the maximum difference between total internal and myocardial volumes, measuring the epicardial contraction rate should be helpful in understanding patients with HFpEF.

Under normal conditions, coronary blood flow (CBF) provides critical blood supply to the myocardium so that it can appropriately meet the metabolic demands of the body. Dogmatically, there exist several known regulators and modulators of CBF that include local metabolites and neurohormonal factors that can influence the function of the coronary circulation. In disease states such as diabetes and myocardial ischemia, these regulators are impaired or shifted such that CBF is reduced. Although functional considerations have been and continued to be well studied, more recent evidence builds upon established studies that collectively suggest that the relative roles of coronary structure, biomechanics, and the influence of cardiac biomechanics via extravascular compression may also play a significant role in dictating CBF. In this mini review, we discuss these regulators of CBF under normal and pathophysiological conditions and their potential influence on the control of CBF.

Hemodynamics play a central role in hemostasis and thrombosis by affecting all aspects linked to platelet functions and coagulation. In vitro flow devices are extensively used in basic research, pharmacological studies, antiplatelet agent screening, and development of diagnostic tools. Because hemodynamic conditions vary tremendously throughout the vascular tree and among different (patho)physiological processes, it is important to use flow conditions based on relevant biorheological reference ranges. Surprisingly, it is particularly difficult to find a concise overview of relevant hemodynamic parameters in various human and mouse vessels. To our knowledge, this is the first time an inventory of flow conditions in healthy, non-diseased, human and mouse vessels has been created. The objective of providing such a repertoire is to aid researchers in the field of hemostasis and thrombosis in choosing rheological conditions relevant in in vitro flow experiments and to promote harmonization of flow-based assays to facilitate comparative evaluations between studies. With reference to the human, we discuss relevant similarities and discrepancies in wall shear rates in the mouse, which are typically one order of magnitude greater in agreement with allometric scaling laws between species. Importantly, we bring the attention of the researchers to the fact that the relevant range of average wall shear rates in human arteries where clinically relevant arterial thrombosis occurs may fall as low as 100 to 200 s(-1), thus significantly overlapping with what are considered "venous" shear rates. The same range for the murine arteries used for arterial thrombosis models may significantly exceed 1000 s(-1) reaching values considered to be "pathological."

The role played by the right ventricular (RV) dysfunction has long been underestimated in clinical practice. Recent findings are progressively confirming that when the RV efficiency deteriorates both the right and the left circulation is (significantly) affected, but studies dedicated to a detailed description of RV hemodynamic role still lack. In response to such a gap in knowledge, this work proposes a numerical model that for the first time evaluates the effect of isolated RV dysfunction on the whole circulation. Lumped parameter modelling was applied to represent the physio-pathological hemodynamics. Different grades of impairment were simulated for three dysfunctions i.e., systolic, diastolic, and combined systolic and diastolic. Hemodynamic alterations (i.e., of blood pressure, flow, global hemodynamic parameters), arising from the dysfunctions, are calculated and analysed. Results well accord with clinical observations, showing that RV dysfunction significantly affects both the pulmonary and systemic hemodynamics. Successful verification against in vivo data proved the clinical potentiality of the model i.e., the capability of identifying the degree of RV impairment for given hemodynamic conditions. This study aims at contributing to the improvement of RV dysfunction recognition and treatment, and to the development of tools for the clinical management of pathologies involving the right heart.

The decompensated univentricular circulation is identified as one of the most challenging conditions and the application of the mechanical circulatory support (MCS) devices is proposed as therapeutic option for Fontan failure. Modelling methodologies are reported to identify the optimized types, extent and duration of required hemodynamic support using MCS. The specific parameters of device-body interaction during support of failing Fontan circulation within the design points of dedicated pediatric ventricular assist devices has not been previously defined. In this work, we introduce a mathematical model developed to simulate the interaction between the Fontan single-ventricular circulation and a constant-flow pediatric ventricular assist device (VAD) Sputnik. The interaction is studied at a pump rotor speed of 5000-9000 rpm. This simulation demonstrates that the pump replacing pulmonary ventricle of the heart creates necessary pressure differential between the systemic veins (7 mmHg) and the pulmonary artery (17.3 mmHg). Moreover, it increases the venous return that, according to the Frank-Starling mechanism, increases the stroke volume up to 32 ml/bpm (26 ml/bpm - without using a pump). For the first time, a simulation for the pediatric VAD Sputnik has been carried out. The simulation results confirm pediatric VAD Sputnik can be a possible tool to normalize the Fontan circulation in pediatric patients.

Object. Use a fluid-structure interaction (FSI) model to do a safety risk assessment of cervical spinal manipulative therapy (cSMT) for atherosclerosis to study the biomechanical changes of atherosclerosis plaque and the carotid artery (CA) lumen during the process of cSMT. Method. The FSI model, based on MRI data of an atherosclerosis patient, was used to simulate the deformations of the plaque and lumen during the process of two kinds of typical cSMT (the high-speed, low-amplitude spinal manipulation and the cervical rotatory manipulation). The biomechanical parameters were recorded, such as the highest wall shear stress (WSS), the maximum plaque wall stress (PWS), the wall tensile stress (Von mises stress, VWTS), and the strain. Result. The max_WSS was 33.77 kPa in the most extensive deformation. The highest WSS region on the plaque surface was also the highest PWS region. The max_PWS in a 12% stretch was 55.11 kPa, which was lower than the rupture threshold. The max_VWTS of the cap in 12% stretch which approached the fracture stress level was 116.75 kPa. Moreover, the vessel’s max_VWTS values in 10% and 12% stretch were 554.21 and 855.19 kPa. They were higher than the fracture threshold, which might cause media fracture. Meanwhile, the 7% stretched strain was 0.29, closed to the smallest experimental green strains at rupture. Conclusion. The carotid arteries’ higher stretch generated the higher stress level of the plaque. Cervical rotatory manipulation might cause plaque at a high risk of rupture in deformation after 12% stretch and more. Lower deformation of the plaque and artery caused by the high-speed, low-amplitude spinal manipulation might be safer.

PURPOSE OF REVIEW Cardiovascular disease (CVD) is the leading cause of death in patients with chronic kidney disease (CKD). However, traditional CVD risk prediction equations do not work well in patients with CKD, and inclusion of kidney disease metrics such as albuminuria and estimated glomerular filtration rate have a modest to no benefit in improving prediction. RECENT FINDINGS As CKD progresses, the strength of traditional CVD risk factors in predicting clinical outcomes weakens. A pooled cohort equation used for CVD risk prediction is a useful tool for guiding clinicians on management of patients with CVD risk, but these equations do not calibrate well in patients with CKD, although a number of studies have developed modifications of the traditional equations to improve risk prediction. The reason for the poor calibration may be related to the fact that as CKD progresses, associations of traditional risk factors such as BMI, lipids and blood pressure with CVD outcomes are attenuated or reverse, and other risk factors may become more important. SUMMARY Large national cohorts such as the US Veteran cohort with many patients with evolving CKD may be useful resources for the developing CVD prediction models; however, additional considerations are needed for the unique composition of patients receiving care in these healthcare systems, including those with multiple comorbidities, as well as mental health issues, homelessness, posttraumatic stress disorders, frailty, malnutrition and polypharmacy. Machine learning over conventional risk prediction models may be better suited to handle the complexity needed for these CVD prediction models.