1D simulation of blood flow characteristics in the circle of Willis using THINkS

One-dimensional (1D) simulation of the complete vascular network, so called THINkS (Total Human Intravascular Network Simulation) is developed to investigate changes of blood flow characteristics caused by the variation of CoW. THINkS contains 158 major veins, 85 major arteries, and 77 venous and 43 arterial junctions. THINkS is validated with available in vivo blood flow waveform data. The overall trends of flow rates in variations of the CoW, such as the missing anterior cerebral artery (missing-A1) or missing posterior cerebral artery (missing-P1), are confirmed by in vivo experimental data. It is demonstrated that the CoW has the ability to shunt blood flow to different areas in the brain. Flow rates in efferent arteries remain unaffected under the variation of CoW, while the flow rates in afferent vessels can be subject to substantial changes. The redistribution of blood flow can cause particular vessels to undergo extra flow rate and hemodynamic stresses.     

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were \(C_1=3000\pm 1700\,\hbox {Pa}\) and \(\alpha =45\pm 25\) (\(b_f=25\pm 14\), \(b_{ft}=11\pm 6\), \(b_{t}=9\pm 5\)) for the HF cases and \(C_1=1700\,\hbox {Pa}\) and \(\alpha =15\) (\(b_f=8\), \(b_{ft}=4\), \(b_{t}=3\)) for the healthy case.     

Numerical simulations of the blood flow in the patient-specific arterial cerebral circle region

The Cerebral Circle Region, also known as the Circle of Willis (CoW), is a loop of arteries that form arterial connections between supply arteries to distribute blood throughout the cerebral mass. Among the population, only 25% to 50% have a complete system of arteries forming the CoW. 3D time-varying simulations for three different patient-specific artery anatomies of CoW were performed in order to gain a better insight into the phenomena existing in the cerebral blood flow. The models reconstructed on the basis of computer tomography images start from the aorta and include the largest arteries that supply the CoW and the arteries of CoW. Velocity values measured during the ultrasound examination have been compared with the results of simulations. It is shown that the flow in the right anterior artery in some cases may be supplied from the left internal carotid artery via the anterior communicating artery. The investigations conducted show that the computational fluid dynamic tool, which provides high resolution in both time and space domains, can be used to support physicians in diagnosing patients of different ages and various anatomical arterial structures.     

3D models of blood flow in the cerebral vasculature

The circle of Willis (CoW) is a ring-like arterial structure located in the base of the brain and is responsible for the distribution of oxygenated blood throughout the cerebral mass. To investigate the effects of the complex 3D geometry and anatomical variability of the CoW on the cerebral hemodynamics, a technique for generating physiologically accurate models of the CoW has been created using a combination of magnetic resonance data and computer-aided design software. A mathematical model of the body's cerebral autoregulation mechanism has been developed and numerous computational fluid dynamics simulations performed to model the hemodynamics in response to changes in afferent blood pressure. Three pathological conditions were explored, namely a complete CoW, a fetal P1 and a missing A1. The methodology of the cerebral hemodynamic modelling is proposed with the potential for future clinical application in mind, as a diagnostic tool.     

Structural alteration of the endothelial glycocalyx: contribution of the actin cytoskeleton

The endothelial glycocalyx is a carbohydrate–protein layer that lines the luminal surface of the endothelium. It anchors to the cell membrane via its core proteins that share extended link to the actin cytoskeleton. It is widely accepted that those protein domains and the attached carbohydrates are susceptible to pathological changes. It is unclear, however, to what extent the actin cytoskeleton contributes to the glycocalyx stability. In this study, we investigate the role of the actin cytoskeleton in the maintenance of the glycocalyx under static and laminar flow conditions in vitro. Our results show that in the static culture medium neither rapid actin depolymerisation nor prolonged actin disturbance leads to glycocalyx disruption from the apical surface of human umbilical vein endothelial cells. However, when endothelial cells are exposed to laminar flow for 24 h, the glycocalyx is seen to shift to the downstream peripheral region of the cell surface. The mean fluorescence intensity decreases to \(91.9 \pm 2.5\%\) of the control. When actin depolymerisation is introduced, the intensity decreases significantly to \(54.7 \pm 1.3\%\), indicating a severe disruption of the glycocalyx. Similar changes are observed in human aortic endothelial cells, where the intensity of the glycocalyx is reduced to \(72.8 \pm 1.6\%\) of the control. Collectively, we demonstrate that the actin cytoskeleton contributes to structural stability of the glycocalyx under shear stress. Our results can be used to develop new strategies to prevent shedding of the glycocalyx in cardiovascular diseases.     

4D subject-specific inverse modeling of the chick embryonic heart outflow tract hemodynamics

Blood flow plays a critical role in regulating embryonic cardiac growth and development, with altered flow leading to congenital heart disease. Progress in the field, however, is hindered by a lack of quantification of hemodynamic conditions in the developing heart. In this study, we present a methodology to quantify blood flow dynamics in the embryonic heart using subject-specific computational fluid dynamics (CFD) models. While the methodology is general, we focused on a model of the chick embryonic heart outflow tract (OFT), which distally connects the heart to the arterial system, and is the region of origin of many congenital cardiac defects. Using structural and Doppler velocity data collected from optical coherence tomography, we generated 4D (\(\hbox {3D}\,+\,\hbox {time}\)) embryo-specific CFD models of the heart OFT. To replicate the blood flow dynamics over time during the cardiac cycle, we developed an iterative inverse-method optimization algorithm, which determines the CFD model boundary conditions such that differences between computed velocities and measured velocities at one point within the OFT lumen are minimized. Results from our developed CFD model agree with previously measured hemodynamics in the OFT. Further, computed velocities and measured velocities differ by \(<\)15 % at locations that were not used in the optimization, validating the model. The presented methodology can be used in quantifications of embryonic cardiac hemodynamics under normal and altered blood flow conditions, enabling an in-depth quantitative study of how blood flow influences cardiac development.     

Experimental study of hemodynamics in the circle of willis

Background

The Circle of Willis (CoW) is an important collateral pathway of the cerebral blood flow. An experimental study of the cerebral blood flow (CBF) distribution in different anatomical variations may help to a better understanding of the collateral mechanism of the CoW.

Methods

An in-vitro test rig was developed to simulate the physiological cerebral blood flow in the CoW. Ten anatomical variations were considered in this study, include a set of different degrees of stenosis in L-ICA and L-ICA occlusion coexist with common anatomical variations. Volume flow rates of efferent arteries and pressure signals at the end of communicating arteries of each case were recorded. Physiological pressure waveforms were applied as inlet boundary condition.

Results

In the development of L-ICA stenosis, the total CBF decreases with the increase of stenosis degree. The blood supply of ipsilateral middle cerebral artery (MCA) was affected most by the stenosis of L-ICA. Anterior communicating artery (ACoA) and ipsilateral posterior communicating artery (PCoA) function as important collateral pathways of cerebral collateral circulation when unilateral stenosis occurred. The blood supply of anterior cerebral circulation was compensated by the posterior cerebral circulation through ipsilateral PCoA when L-ICA stenosis degree is greater than 40% and the affected side was compensated immediately by the unaffected side through ACoA. Blood flow of the anterior circulation and the total CBF reached the minimum among all cases studied when L-ICA occlusion coexist with the absence of PCoA.

Conclusion

The results demonstrated the flow distribution patterns of the CoW under anatomical variations and clarified the collateral mechanism of the CoW. The flow ACoA is the most sensitive indexes to the morphology change of ipsilateral ICA. The relative independence of the circulation in anterior and posterior sections of the CoW is not broken and the function of ipsilateral PCoA is not activated until a severe stenosis of unilateral ICA occurs. PCoA is the most important collateral pathway of the collateral circulation and the missing of PCoA has the highest risk of stroke when the ipsilateral ICA has severe stenosis. These findings may provide the basis for future therapeutic and diagnosis applications.

     

Twist buckling behavior of arteries

Arteries are often subjected to torsion due to body movement and surgical procedures. While it is essential that arteries remain stable and patent under twisting loads, the stability of arteries under torsion is poorly understood. The goal of this work was to experimentally investigate the buckling behavior of arteries under torsion and to determine the critical buckling torque, the critical buckling twist angle, and the buckling shape. Porcine common carotid arteries were slowly twisted in vitro until buckling occurred while subjected to a constant axial stretch ratio (1.1, 1.3, 1.5 (in vivo level) and 1.7) and lumen pressure (20, 40, 70 and 100 mmHg). Upon buckling, the arteries snapped to form a kink. For a group of six arteries, the axial stretch ratio significantly affected the critical buckling torque (\(p<0.002\)) and the critical buckling twist angle (\(p<0.001\)). Lumen pressure also significantly affected the critical buckling torque (\(p<0.001\)) but had no significant effect on the critical twist angle (\(p=0.067\)). Convex material constants for a Fung strain energy function were determined and fit well with the axial force, lumen pressure, and torque data measured pre-buckling. The material constants are valid for axial stretch ratios, lumen pressures, and rotation angles of 1.3–1.5, 20–100 mmHg, and 0–270\(^\circ \), respectively. The current study elucidates the buckling behavior of arteries under torsion and provides new insight into mechanical instability of blood vessels.     

Model-based cardiovascular disease diagnosis: a preliminary in-silico study

In this study, we developed and examined the feasibility of a model-based system identification approach to cardiovascular disease diagnosis. The basic premise of the approach is that it may be possible to diagnose cardiovascular disease from disease-induced alterations in the arterial mechanical properties manifested in the proximal and distal arterial blood pressure waveforms. It first individualizes the lumped-parameter model of wave propagation and reflection in the artery using the measurement of proximal and distal arterial blood pressure waveforms. Then, it employs a diagnosis logic, in the form of disease-specific patterns in model parameters, referred as \(\alpha , \beta \) and pulse transit time. The longitudinal change in these parameters is used to diagnose the presence of peripheral artery disease and arterial stiffening. We illustrated the feasibility of the proposed approach by testing it in a full-scale in-silico arterial tree simulation. The results showed that the approach exhibited superior sensitivity to ankle-brachial index and convenience to carotid-femoral pulse wave velocity: The model parameters \(\alpha \) and \(\beta \) responded with up to 100 and 40 % changes to peripheral artery disease with up to 50 % arterial blockage whereas the change in ankle-brachial index was \({<}5\,\%\); the same parameters responded with up to 300 and 40 % changes to up to 100 % arterial stiffening while pulse transit time changed by up to 24 %. Together with the development of more convenient techniques for the measurement of arterial blood pressure waveforms, the proposed approach may evolve into a viable alternative to the state-of-the-art techniques for cardiovascular disease diagnosis.     

Beyond the aorta: partial transmission of reflected waves from aortic coarctation into supra-aortic branches modulates cerebral hemodynamics and left ventricular load

Wave reflection from the site of aortic coarctation produces a reflected backward compression wave (BCW) that raises left ventricular (LV) afterload. However, not all reflected wave power will propagate back to the LV. This study investigated the hypothesis that the BCW is partially transmitted into supra-aortic vessels as a forward wave and explored the consequences of this phenomenon for cerebral and LV haemodynamic load. In eight sheep, high fidelity pressure and flow were measured in the aortic trunk (AoT) and brachiocephalic trunk (BCT, the single supra-aortic vessel present in sheep) at baseline and during two levels of proximal descending aortic constriction. Wave power analysis showed that aortic constriction produced not only a BCW in the AoT, but also a second forward compression wave (\(\mathrm{FCW}_{2})\) in the BCT that augmented pressure and flow after the initial forward compression wave (\(\mathrm{FCW}_{1})\). Mathematical analysis and a one-dimensional model of the human systemic arteries and aortic coarctation suggested that the relative transmission of waves into supra-aortic vessels versus the aorta was determined by the relative admittances of these vessels. Reducing supra-aortic admittance (1) increased pressure and flow pulsatility in cerebral arteries, (2) produced carotid and middle cerebral arterial flow waveforms with an older adult phenotype, (3) promoted transmission of reflected wave power towards the LV and (4) substantially increased mid- to late-systolic myocardial stress, which may promote LV hypertrophy. These findings suggest that wave transmission into supra-aortic branches has an important impact on both cerebral hemodynamics and LV load in aortic coarctation.     

Simulation of extracellular matrix remodeling by fibroblast cells in soft three-dimensional bioresorbable scaffolds

To culture functional soft tissues and organs in three-dimensional scaffolds, it is essential to elucidate the optimal scaffold mechanical properties. However, mechanoregulated soft tissue remodeling is not well understood. In this study, we hypothesized that individual cells are capable of remodeling extracellular matrix within a short proximity of themselves in order to match the stiffness of the broader surrounding matrix. This theory was implemented in a three-dimensional finite element model to simulate soft tissue remodeling of human fibroblast cells in two collagen–chitosan scaffolds with different mechanical properties. Simulation results closely matched with previously reported experimental data, showing that soft tissue growth in compliant (Scaf-A, 4.30 kPa) and stiff (Scaf-B, 17.03 kPa) scaffolds led to an almost eightfold difference in the resulting stiffnesses after 10 days (8.40 kPa for Scaf-A, 59.25 kPa for Scaf-B). Furthermore, varying the simulated rate for tissue remodeling by \(\pm \)50 % caused unequal changes in the resulting stiffness (+3.6 and \(-\)23 % for Scaf-A, +5 and \(-\)17 % for Scaf-B), and \(\pm \)50 % changes in the assumed upper limit on tissue stiffness only had larger effects on the stiff scaffold (+42 and \(-\)44 % for Scaf-B). These results reinforce the notion that soft tissue remodeling is governed by the stiffness of the surrounding matrix, until meeting an upper limit on tissue stiffness.     

Backbone resonance assignments for the SET domain of the human methyltransferase NSD2

Aberrant NSD2 methyltransferase activity is implicated as the oncogenic driver in multiple myeloma, suggesting opportunities for novel therapeutic intervention. The methyltransferase activity of NSD2 resides in its catalytic SET domain, which is conserved among most lysine methyltransferases. Here we report the backbone \(\hbox {H}^{\mathrm{N}}\), N, C\(^{\prime }\), \(\hbox {C}^\alpha\) and side-chain \(\hbox {C}^\beta\) assignments of a 25 kDa NSD2 SET domain construct, spanning residues 991–1203. A chemical shift analysis of C\(^{\prime }\), \(\hbox {C}^\alpha\) and \(\hbox {C}^\beta\) resonances predicts a secondary structural pattern that is in agreement with homology models.     

Model of cellular mechanotransduction via actin stress fibers

Mechanical stresses due to blood flow regulate vascular endothelial cell structure and function and play a key role in arterial physiology and pathology. In particular, the development of atherosclerosis has been shown to correlate with regions of disturbed blood flow where endothelial cells are round and have a randomly organized cytoskeleton. Thus, deciphering the relation between the mechanical environment, cell structure, and cell function is a key step toward understanding the early development of atherosclerosis. Recent experiments have demonstrated very rapid (\(\sim \)100 ms) and long-distance (\(\sim \)10 \(\upmu \)m) cellular mechanotransduction in which prestressed actin stress fibers play a critical role. Here, we develop a model of mechanical signal transmission within a cell by describing strains in a network of prestressed viscoelastic stress fibers following the application of a force to the cell surface. We find force transmission dynamics that are consistent with experimental results. We also show that the extent of stress fiber alignment and the direction of the applied force relative to this alignment are key determinants of the efficiency of mechanical signal transmission. These results are consistent with the link observed experimentally between cytoskeletal organization, mechanical stress, and cellular responsiveness to stress. Based on these results, we suggest that mechanical strain of actin stress fibers under force constitutes a key link in the mechanotransduction chain.     

Ultrasound palpation for fast in-situ quantification of articular cartilage stiffness,thickness and relaxation capacity

Most current cartilage testing devices require the preparation of excised samples and therefore do not allow intra-operative application for diagnostic purposes. The gold standard during open or arthroscopic surgery is still the subjective perception of manual palpation. This work presents a new diagnostic method of ultrasound palpation (USP) to acquire applied stress and strain data during manual palpation of articular cartilage. With the proposed method, we obtain cartilage thickness and stiffness. Moreover, repeated palpations allow the quantification of relaxation effects. USP measurements on elastomer phantoms demonstrated very good repeatability for both, stage-guided (97.2%) and handheld (96.0%) applications. The USP measurements were compared with conventional indentation experiments and revealed very good agreement on elastomer phantoms (\(r = 0.98\)) and good agreement on porcine cartilage samples (\(r = 0.76\)). Artificially degenerated cartilage samples showed reduced stiffness, weak capacity to relax after palpation and an increase of stiffness of approximately 50% with each single palpation. Intact cartilage was measured by USP directly at the patella (in situ) and after excision and removal of the subchondral bone (ex situ), leading to stiffness values of \(12.1\pm 5.5\) and \(8.5\pm 5.9\,\hbox {MPa}\) (\(p<0.05\)), respectively. The results demonstrate the potential of the USP system for cartilage testing, its sensitivity to degenerative changes and as a method for quantifying relaxation processes by means of repeated palpations. Furthermore, the differences in the results of in-situ and ex-situ measurements are of general interest, since such comparison has not been reported previously. We point out the limited comparability of ex-situ cartilage with its in-situ biomechanical behavior.     

Characterization of the mechanical behavior of the optic nerve sheath and its role in spaceflight-induced ophthalmic changes

Visual impairment and intracranial pressure (VIIP) syndrome is characterized by a number of permanent ophthalmic changes, including loss of visual function. It occurs in some astronauts during long-duration spaceflight missions. Thus, understanding the pathophysiology of VIIP is currently a major priority in space medicine research. It is hypothesized that maladaptive remodeling of the optic nerve sheath (ONS), in response to microgravity-induced elevations in intracranial pressure (ICP), contributes to VIIP. However, little is known about ONS biomechanics. In this study, we developed a custom mechanical testing system that allowed for unconfined lengthening, twisting, and circumferential distension of the porcine ONS during inflation and axial loading. Data were fit to a four-fiber family constitutive equation to extract material and structural parameters. Inflation testing showed a characteristic “cross-over point” in the pressure–diameter curves under different axial loads in all samples that were tested; the cross-over pressure was \(10.3 \pm 0.95\) mmHg (\(\hbox {mean} \pm \hbox {SEM}\)). Large sample-to-sample variations were observed in the circumferential strain, while only modest variations were observed in the circumferential stress. Multiphoton microscopy revealed that the collagen fibers of the ONS were primarily oriented axially when the tissue was loaded. The existence of this cross-over behavior is expected to be neuroprotective, as it would avoid optic nerve compression during routine changes in gaze angle, so long as ICP was within the normal range. Including these observations into computational models of VIIP will help provide insight into the pathophysiology of VIIP and could help identify risk factors and potential interventions.     

Study of the collateral capacity of the circle of Willis of patients with severe carotid artery stenosis by 3D computational modeling

This numerical study aims to investigate the capacity of the circle of Wills (CoW) to provide collateral blood supply for patients with unilateral carotid arterial stenosis. The basic 3D geometry of the CoW was reconstructed based on a magnetic resonance angiogram of a normal human subject. A total of 52 computational fluid dynamics simulations were performed for four geometry configurations of the CoW with an artificially inserted axisymmetric stenosis of different luminal area reductions in an internal carotid artery (ICA) under a variety of boundary conditions. The CoW geometric configurations included (a) a normal CoW with all communicating arteries; (b) as model (a) but with enlarged communicating arterial diameters; (c) as (a) but with the ipsilateral posterior communicating artery missing, and (d) as (c) but with enlarged communicating arteries. It is found that the blood perfusion pressure drop between the ipsilateral ICA and the middle cerebral artery (MCA) only becomes significant when the degree of stenosis is greater than 86%. The cerebral autoregulation range varied significantly between the different CoW configurations for the severe stenosis cases. Without causing the flow rates to decrease at the efferent arterial ends, the mean perfusion pressure in the ipsilateral ICA can drop from 100 to 73, 67, 92 and 84mmHg for the CoW models (a)-(d) with 96% luminal area reduction stenosis, respectively. The additional pathways are able to raise the ipsilateral MCA pressure significantly without reducing the total flow perfusion. Cerebral autoregulation effects were not directly included in the study. Therefore, the findings in the study should be interpreted with cautions when comes to the biological and clinical significance.     

Cell Volume Regulation in the Proximal Tubule of Rat Kidney

We developed a dynamic model of a rat proximal convoluted tubule cell in order to investigate cell volume regulation mechanisms in this nephron segment. We examined whether regulatory volume decrease (RVD), which follows exposure to a hyposmotic peritubular solution, can be achieved solely via stimulation of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. We also determined whether regulatory volume increase (RVI), which follows exposure to a hyperosmotic peritubular solution under certain conditions, may be accomplished by activating basolateral \(\hbox {Na}^+\)/H\(^+\) exchangers. Model predictions were in good agreement with experimental observations in mouse proximal tubule cells assuming that a 10% increase in cell volume induces a fourfold increase in the expression of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. Our results also suggest that in response to a hyposmotic challenge and subsequent cell swelling, \(\hbox {Na}^+\)–\(\hbox {HCO}^-_3\) cotransporters are more efficient than basolateral K\(^+\) and \(\hbox {Cl}^-\) channels at lowering intracellular osmolality and reducing cell volume. Moreover, both RVD and RVI are predicted to stabilize net transcellular \(\hbox {Na}^+\) reabsorption, that is, to limit the net \(\hbox {Na}^+\) flux decrease during a hyposmotic challenge or the net \(\hbox {Na}^+\) flux increase during a hyperosmotic challenge.     

Understanding the effects of isolation on seed and pollen flow,spatial genetic structure and effective population size of the dioecious tropical tree species <Emphasis Type="Italic">Myracrodruon urundeuva</Emphasis>

This study examines the levels of gene flow, the distance and the patterns of pollen and seed dispersal, the intra-population spatial genetic structure (SGS) and the effective population size of a spatially isolated Myracrodruon urundeuva population using five microsatellite loci. The study was carried out in the Paulo de Faria Ecological Station, São Paulo State, Brazil and included the sampling and mapping of 467 adult-trees and 149 juveniles. Open-pollinated seeds (514) from 29 seed-trees were also sampled and genotyped. Significant SGS was detected in both adult (S _p  = 0.0269) and juveniles trees (S _p  = 0.0246), indicating short-distance seed dispersal. Using maternity analysis, all juveniles had the mother-tree assigned within the stand. A father-tree within the stand was also assigned for 97.3% of the juveniles and 98.4% of offspring. The average pollen dispersal distance measured in juveniles \( \left( {\hat{\delta } = 1 3 8\pm 1 6 9 {\text{ m}},{\text{ mean}} \pm {\text{SD}}} \right) \) and offspring \( \left( {\hat{\delta } = 2 5 2\pm 20 4 {\text{ m}}} \right) \) were higher than the average seed dispersal distance measured in juveniles \( \left( {\hat{\delta } = 1 2 4\pm 1 50{\text{ m}}} \right) \). About 70% of the pollen from juveniles and 51% from offspring traveled less than 200 m and, 72% of the seeds traveled less than 50 m. The effective population size of the studied sample indicates that the 467 adult-trees and 145 juveniles correspond respectively to 335 and 63 individuals that are neither inbred nor relatives. The results are discussed in relation to their impact on seed collection practices and genetic conservation.