Search

Cardiovascular disease is the leading cause of morbidity and mortality in the Western World. The inability of fully differentiated, load-bearing cardiovascular tissues to in vivo regenerate and the limitations of the current treatment therapies greatly motivate the efforts of cardiovascular tissue engineering to become an effective clinical strategy for injured heart and vessels. For the effective production of organized and functional cardiovascular engineered constructs in vitro, a suitable dynamic environment is essential, and can be achieved and maintained within bioreactors. Bioreactors are technological devices that, while monitoring and controlling the culture environment and stimulating the construct, attempt to mimic the physiological milieu. In this study, a review of the current state of the art of bioreactor solutions for cardiovascular tissue engineering is presented, with emphasis on bioreactors and biophysical stimuli adopted for investigating the mechanisms influencing cardiovascular tissue development, and for eventually generating suitable cardiovascular tissue replacements.

Bone regeneration can be accelerated by localized delivery of appropriate growth factors/biomolecules. Localized delivery can be achieved by a 2-level system: (i) incorporation of biomolecules within biodegradable particulate carriers (nanoparticles), and (ii) inclusion of such particulate carriers (nanoparticles) into suitable porous scaffolds. In this study, freeze-dried porous chitosan-gelatin scaffolds (CH-G: 1:2 ratio by weight) were embedded with various amounts of poly(lactide-co-glycolide) (PLGA) nanoparticles, precisely 16.6%, 33.3% and 66.6% (respect to CH-G weight). Scaffolds loaded with PLGA nanoparticles were subjected to physico-mechanical and biological characterizations including morphological analysis, swelling and dissolution tests, mechanical compression tests and cell viability tests. Results showed that incorporation of PLGA nanoparticles into porous crosslinked CH-G scaffolds: (i) changed the micro-architecture of the scaffolds in terms of mean pore diameter and pore size distribution, (ii) reduced the dissolution degree of the scaffolds, and (iii) increased the compressive modulus. On the other hand, the water uptake behavior of CH-G scaffolds containing PLGA nanoparticles significantly decreased. The incorporation of PLGA nanoparticles did not affect the biocompatibility of CH-G scaffolds., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Here we present a study on the impact of assumptions on image-based hemodynamic simulations of healthy carotid bifurcations. In particular, we evaluate to which extent assumptions on blood rheology influence bulk flow features, driven by the fact that few studies have provided adequate insights into the influence of assumptions to confidently model the 4D hemodynamics within the bifurcation. The final goal is to complement, integrate and extend with a quantitative characterization of the bulk flow the description currently adopted to classify altered hemodynamics, which is based on wall shear stress (WSS). Hemodynamic simulations of two image-based carotid bifurcation geometries were carried out assuming a reference Newtonian viscosity, two non-Newtonian rheology models and Newtonian viscosities based on characteristic shear rates. WSS-based and Lagrangian-based metrics for helical flow quantification and for vorticity dynamics quantification were calculated. Our findings suggest that the assumption of Newtonian rheology: (1) could be reasonable for bulk flow metrics (differences from non-Newtonian behavior are lower than 10%); (2) influences at different levels the WSS-based indicators, depending on the bifurcation model, even if in our study it is lower than the major source of uncertainty as recognized by the literature (i.e., uncertainty on geometry reconstruction)., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

The hemodynamics within the aorta of five healthy humans were investigated to gain insight into the complex helical flow patterns that arise from the existence of asymmetries in the aortic region. The adopted approach is aimed at (1) overcoming the relative paucity of quantitative data regarding helical blood flow dynamics in the human aorta and (2) identifying common characteristics in physiological aortic flow topology, in terms of its helical content. Four-dimensional phase-contrast magnetic resonance imaging (4D PC MRI) was combined with algorithms for the calculation of advanced fluid dynamics in this study. These algorithms allowed us to obtain a 4D representation of intra-aortic flow fields and to quantify the aortic helical flow. For our purposes, helicity was used as a measure of the alignment of the velocity and the vorticity. There were two key findings of our study: (1) intra-individual analysis revealed a statistically significant difference in the helical content at different phases of systole and (2) group analysis suggested that aortic helical blood flow dynamics is an emerging behavior that is common to normal individuals. Our results also suggest that helical flow might be caused by natural optimization of fluid transport processes in the cardiovascular system, aimed at obtaining efficient perfusion. The approach here applied to assess in vivo helical blood flow could be the starting point to elucidate the role played by helicity in the generation and decay of rotating flows in the thoracic aorta.

The actin microfilament (F-actin) is a structural and functional component of the cell cytoskeleton. Notwithstanding the primary role it plays for the mechanics of the cell, the mechanical behaviour of F-actin is still not totally explored. In particular, the relationship between the mechanics of F-actin and its molecular architecture is not completely understood. In this study, the mechanical properties of F-actin were related to the molecular topology of its building monomers (G-actin) by employing a computational multi-level approach. F-actins with lengths up to 500 nm were modelled and characterized, using a combination of equilibrium molecular dynamics (MD) simulations and normal mode analysis (NMA). MD simulations were performed to analyze the molecular rearrangements of G-actin in physiological conditions; NMA was applied to compute the macroscopic properties of F-actin from its vibrational modes of motion. Results from this multi-level approach showed that bending stiffness, bending modulus and persistence length are independent from the length of F-actin. On the contrary, the orientations and motions of selected groups of residues of G-actin play a primary role in determining the filament flexibility. In conclusion, this study (i) demonstrated that a combined computational approach of MD and NMA allows to investigate the biomechanics of F-actin taking into account the molecular topology of the filament (i.e., the molecular conformations of G-actin) and (ii) that this can be done using only crystallographic G-actin, without the need of introducing experimental parameters nor of reducing the number of residues., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

Although flow-driven mechanisms associated with vascular physiopathology also deal with four-dimensional phenomena such as species transport, the majority of the research on the subject focuses primarily on wall shear stress as indicator of disturbed flow. Indeed, the role that bulk flow plays in vascular physiopathology has not been thoroughly investigated, partly because of a lack of descriptors that would be able to reduce the intricacy of arterial hemodynamics. Here, an approach is proposed to investigate, in silico, the bulk flow within the carotid bifurcation. For this purpose, we coupled a three-dimensional model of carotid bifurcation with a lumped model of the downstream vasculature. For the sake of comparison, we also imposed three different fixed flow rate repartitions between the internal and external carotid arteries on the three-dimensional model. The bulk flow was characterized by applying a descriptor of helical motion, the helical flow index (HFI) to the model; the HFI has recently been shown to provide an accurate representation of complex flows. Moreover, a new metric is presented to investigate the vorticity dynamics within the bifurcation. Our results highlight the effectiveness of these metrics in the following contexts: (i) identifying and ranking emerging hemodynamic features and (ii) quantifying the influence of the outflow boundary conditions on the composition of the translational and rotational components of the fluid motion. The metrics applied herein allow for a more comprehensive analysis, which may lead to the development of an instrument to relate the bulk flow to vascular pathophysiological events that involve not only fluid-related forces, but also transport phenomena within blood.

Microtubules are supramolecular structures that make up the cytoskeleton and strongly affect the mechanical properties of the cell. Within the cytoskeleton filaments, the microtubule (MT) exhibits by far the highest bending stiffness. Bending stiffness depends on the mechanical properties and intermolecular interactions of the tubulin dimers (the MT building blocks). Computational molecular modeling has the potential for obtaining quantitative insights into this area. However, to our knowledge, standard molecular modeling techniques, such as molecular dynamics (MD) and normal mode analysis (NMA), are not yet able to simulate large molecular structures like the MTs; in fact, their possibilities are normally limited to much smaller protein complexes. In this work, we developed a multiscale approach by merging the modeling contribution from MD and NMA. In particular, MD simulations were used to refine the molecular conformation and arrangement of the tubulin dimers inside the MT lattice. Subsequently, NMA was used to investigate the vibrational properties of MTs modeled as an elastic network. The coarse-grain model here developed can describe systems of hundreds of interacting tubulin monomers (corresponding to up to 1,000,000 atoms). In particular, we were able to simulate coarse-grain models of entire MTs, with lengths up to 350 nm. A quantitative mechanical investigation was performed; from the bending and stretching modes, we estimated MT macroscopic properties such as bending stiffness, Young modulus, and persistence length, thus allowing a direct comparison with experimental data., (Copyright © 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Motivation: Nonlinear small datasets, which are characterized by low numbers of samples and very high numbers of measures, occur frequently in computational biology, and pose problems in their investigation. Unsupervised hybrid-two-phase (H2P) procedures-specifically dimension reduction (DR), coupled with clustering-provide valuable assistance, not only for unsupervised data classification, but also for visualization of the patterns hidden in high-dimensional feature space., Methods: 'Minimum Curvilinearity' (MC) is a principle that-for small datasets-suggests the approximation of curvilinear sample distances in the feature space by pair-wise distances over their minimum spanning tree (MST), and thus avoids the introduction of any tuning parameter. MC is used to design two novel forms of nonlinear machine learning (NML): Minimum Curvilinear embedding (MCE) for DR, and Minimum Curvilinear affinity propagation (MCAP) for clustering., Results: Compared with several other unsupervised and supervised algorithms, MCE and MCAP, whether individually or combined in H2P, overcome the limits of classical approaches. High performance was attained in the visualization and classification of: (i) pain patients (proteomic measurements) in peripheral neuropathy; (ii) human organ tissues (genomic transcription factor measurements) on the basis of their embryological origin., Conclusion: MC provides a valuable framework to estimate nonlinear distances in small datasets. Its extension to large datasets is prefigured for novel NMLs. Classification of neuropathic pain by proteomic profiles offers new insights for future molecular and systems biology characterization of pain. Improvements in tissue embryological classification refine results obtained in an earlier study, and suggest a possible reinterpretation of skin attribution as mesodermal., Availability: https://sites.google.com/site/carlovittoriocannistraci/home.

Purpose: Since stretching plays a key role in skeletal muscle tissue development in vivo, by making use of an innovative bioreactor and a biodegradable microfibrous scaffold (DegraPol(R)) previously developed by our group, we aimed to investigate the effect of mechanical conditioning on the development of skeletal muscle engineered constructs, obtained by seeding and culturing murine skeletal muscle cells on electrospun membranes., Methods: Following 5 days of static culture, skeletal muscle constructs were transferred into the bioreactor and further cultured for 13 days, while experiencing a stretching pattern adapted from the literature to resemble mouse development and growth. Sample withdrawal occurred at the onset of cyclic stretching and after 7 and 10 days. Myosin heavy chain (MHC) accumulation in stretched constructs (D) was evaluated by Western blot analysis and immunofluorescence staining, using statically cultured samples (S) as controls., Results: Western blot analysis of MHC on dynamically (D) and statically (S) cultured constructs at different time points showed that, at day 10, the applied stretching pattern led to an eight-fold increase in myosin accumulation in cyclically stretched constructs (D) with respect to the corresponding static controls (S). These results were confirmed by immunofluorescence staining of total sarcomeric MHC., Conclusions: Since previous attempts to reproduce skeletal myogenesis in vitro mainly suffered from the difficulty of driving myoblast development into an architecturally organized array of myosin expressing myotubes, the chance of inducing MHC accumulation via mechanical conditioning represents a significant step towards the generation of a functional muscle construct for skeletal muscle tissue engineering applications.

Purpose: Knowledge of the mechanical behavior of myosin and actin monomer is critical for understanding the molecular mechanism of actomyosin-based muscle and non-muscle motility. Different experimental studies concerning actomyosin interaction have been performed in vitro, but studies at the single molecule level have just begun. The aim of this study was to provide a mechanical characterization of myosin II and actin monomer using a numerical approach., Methods: The elastic properties of the two proteins involved in muscle contraction were assessed by performing stretching simulations up to 10% protein elongation using the restraining method. Interaction properties of the actomyosin complex were evaluated at eight intermolecular distances during which the entire system was left free to move., Results: According to our results, the values of the elastic modulus of the myosin motor domain and actin are 0.30 GPa, and 0.08 GPa, respectively. As for the actomyosin complex, the interaction force has a maximum value of 541.15 pN., Conclusions: Mechanical properties of molecular motors are currently being debated. Our results match a number of experimental data, therefore, supporting the idea that molecular mechanics may be a powerful tool to find a way in this complex subject.

Denoising is a fundamental early stage in 2-DE image analysis strongly influencing spot detection or pixel-based methods. A novel nonlinear adaptive spatial filter (median-modified Wiener filter, MMWF), is here compared with five well-established denoising techniques (Median, Wiener, Gaussian, and Polynomial-Savitzky-Golay filters; wavelet denoising) to suggest, by means of fuzzy sets evaluation, the best denoising approach to use in practice. Although median filter and wavelet achieved the best performance in spike and Gaussian denoising respectively, they are unsuitable for contemporary removal of different types of noise, because their best setting is noise-dependent. Vice versa, MMWF that arrived second in each single denoising category, was evaluated as the best filter for global denoising, being its best setting invariant of the type of noise. In addition, median filter eroded the edge of isolated spots and filled the space between close-set spots, whereas MMWF because of a novel filter effect (drop-off-effect) does not suffer from erosion problem, preserves the morphology of close-set spots, and avoids spot and spike fuzzyfication, an aberration encountered for Wiener filter. In our tests, MMWF was assessed as the best choice when the goal is to minimize spot edge aberrations while removing spike and Gaussian noise.

Thromboembolism and the attendant risk of cardioembolic stroke remains an impediment to the development of prosthetic cardiovascular devices. In particular, altered haemodynamics are implicated in the acute blood cell damage that leads to thromboembolic complications, with platelet activation being the underlying mechanism for cardioemboli formation in blood flow past mechanical heart valves (MHVs) and other blood re-circulating devices. In this work, a new modeling paradigm for evaluating the cardioembolic risk of MHVs is described. In silico fluid-structure interaction (FSI) approach is used for providing a realistic representation of the flow through a bileaflet MHV model, and a Lagrangian analysis is adopted for characterizing the mechanism of mechanically induced activation of platelets by means of a mathematical model for platelet activation state prediction. Additionally, the relationship between the thromboembolic potency of the device and the local flow dynamics is quantified by giving a measure of the role played by the local streamwise and spanwise vorticity components. Our methodology indicates that (i) mechanically induced activation of platelets when passing through the valve is dependent on the phase of the cardiac cycle, where the platelet rate of activation is lower at early systole than late systole; (ii) local spanwise vorticity has greater influence on the activation of platelets (R>or=0.94) than streamwise vorticity (R>or=0.78). In conclusion, an integrated Lagrangian description of key flow characteristics could provide a more complete and quantitative picture of blood flow through MHVs and its potential to activate platelets: the proposed "comprehensive scale" approach could represent an efficient and novel assessment tool for MHV performance and may possibly lead to improved valve designs.

A comprehensive computational study modelling the operation of a rotating hollow-fiber bioreactor for artificial liver (BAL) was performed to explore the interactions between the oxygenated culture medium and the cultured hepatocytes. Computational fluid dynamics investigations were carried out using two-dimensional (2D) and 3D time-dependent numerical simulations, integrating calculations of diffusion, convection, and multiphase fluid dynamics. The analysis was aimed at determining the rotational speed value of the chamber to ensure homogenous distribution of the floating microcarrier-attached aggregated cells (microCAACs) and avoid their sedimentation and excessive packing, analyzing oxygen (O(2)) delivery and cellular O(2) consumption as an index of cellular metabolic activity, and analyzing the fluid-induced mechanical stress experienced by cells. According to our results, homogeneous distribution of cells is reached at a rotational speed of 30 rpm; spreading of cellular concentration at around the initial value of 12% was limited (median = 11.97%, 5th percentile = 10.94%, 95th percentile = 13.2%), resulting in uniform suspension of microCAACs, which did not appear to be excessively packed. Mixing within the rotating fluid caused a maximum fluid-induced stress value of 0.05 Pa, which was neither endangering for liver-specific functions of cultured cells, nor causing disruption of the floating aggregates. Moreover, an inlet medium flow rate of 200 mL/m with a partial pressure of oxygen (pO(2)) value of 160 mmHg was found to guarantee an adequate O(2) supply for the hepatocytes (2.7 x 10(8) hepatocytes are simulated); under such conditions, the minimum pO(2) value (23 mmHg) is above the critical threshold value, causing the onset of cellular hypoxia (10 mmHg). We proved that numerical simulation of transport phenomena is a valuable tool for the computer-aided design of BALs, helping overcome the unsolved issues in optimizing the cell-environment conditioning procedure in rotating BALs.

The mechanics of blood flow in arteries plays a key role in the health of individuals. In this framework, the role played by the presence of helical flow in the human aorta is still not clear in its relation to physiology and pathology. We report here a method for quantifying helical flow in vivo employing time-resolved cine phase contrast magnetic resonance imaging to obtain the complete spatio-temporal description of the three-dimensional pulsatile blood flow patterns in aorta. The method is applied to data of one healthy volunteer. Particle traces were calculated from velocity data: to them we applied a Lagrangian-based method for helical flow quantification, the Helical Flow Index, which has been developed and evaluated in silico in order to reveal global organization of blood flow. Our results: (i) put in evidence that the systolic hemodynamics in aorta is characterized by an evolving helical flow (we quantified a 24% difference in terms of the content of helicity in the streaming blood, between mid and early systole); (ii) indicate that in the first part of the systole helicity is ascrivable mainly to the asymmetry of blood flow in the left ventricle, joined with the laterality of the aorta. In conclusion, this study shows that the quantification of helical blood flow in vivo is feasible, and it might allow detection of anomalies in the expected physiological development of helical flow in aorta and accordingly, could be used in a diagnostic/prognostic index for clinical practice.

The assessment of cardiovascular function by means of arterial pulse wave analysis (PWA) is well established in clinical practice. PWA is applied to study risk stratification in hypertension, with emphasis on the measurement of the augmentation index as a measure of aortic pressure wave reflections. Despite the fact that the prognostic power of PWA, in its current form, still remains to be demonstrated in the general population, there is general agreement that analysis and interpretation of the waveform might provide a deeper insight in cardiovascular pathophysiology. We propose here the use of wavelet analysis (WA) as a tool to quantify arterial pressure waveform features, with a twofold aim. First, we discuss a specific use of wavelet transform in the study of pressure waveform morphology, and its potential role in ascertaining the dynamics of temporal properties of arterial pressure waveforms. Second, we apply WA to evaluate a database of carotid artery pressure waveforms of healthy middle-aged women and men. Wavelet analysis has the potential to extract specific features (wavelet details), related to wave reflection and aortic valve closure, from a measured waveform. Analysis showed that the fifth detail, one of the waveform features extracted applying the wavelet decomposition, appeared to be the most appropriate for the analysis of carotid artery pressure waveforms. What remains to be assessed is how the information embedded in this detail can be further processed and transformed into quantitative data, and how it can be rendered useful for automated waveform classification and arterial function parameters with potential clinical applications.

In this work we used molecular simulations to investigate the elastic properties of collagen single chain and triple helix with the aim of understanding its features starting from first principles. We analysed ideal collagen peptides, homotrimeric and heterotrimeric collagen type I and pathological models of collagen. Triple helices were found much more rigid than single chains, thus enlightening the important role of interchain stabilizing forces, like hydrophobic interaction and hydrogen bonds. We obtained Young's moduli close to 4.5GPa for the ideal model of collagen and for the physiological heterotrimer, while the physiological homotrimer presented a Young's modulus of 2.51GPa, that can be related to a mild form of Osteogenesis Imperfecta in which only the homotrimeric form of collagen type I is produced. Otherwise, the pathological model (presenting a glycine to alanine substitution) showed an elastic modulus of 4.32GPa, thus only slightly lower than the ideal model. This suggests that this mutation only slightly affects the mechanical properties of the collagen molecule, but possibly acts on an higher scale, such as the packing of collagen fibrils.

In the current scientific literature, particular attention is dedicated to the study of the mitral valve and to comprehension of the mechanisms that lead to its normal function, as well as those that trigger possible pathological conditions. One of the adopted approaches consists of computational modelling, which allows quantitative analysis of the mechanical behaviour of the valve by means of continuum mechanics theory and numerical techniques. However, none of the currently available models realistically accounts for all of the aspects that characterize the function of the mitral valve. Here, a new computational model of the mitral valve has been developed from in vivo data, as a first step towards the development of patient-specific models for the evaluation of annuloplasty procedures. A structural finite-element model of the mitral valve has been developed to account for all of the main valvular substructures. In particular, it includes the real geometry and the movement of the annulus and papillary muscles, reconstructed from four-dimensional ultrasound data from a healthy human subject, and a realistic description of the complex mechanical properties of mitral tissues. Preliminary simulations allowed mitral valve closure to be realistically mimicked and the role of annulus and papillary muscle dynamics to be quantified.

In order to investigate the reliability of the so called mean velocity/vessel area formula adopted in clinical practice for the estimation of the flow rate using an intravascular Doppler guide wire instrumentation, a multiscale computational model was used to give detailed predictions on flow profiles within Y-shaped coronary artery bypass graft (CABG) models. At this purpose three CABG models were built from clinical patient's data and used to evaluate and compare, in each model, the computed flow rate and the flow rate estimated according to the assumption of parabolic velocity profile. A consistent difference between the exact and the estimated value of the flow rate was found in every branch of all the graft models. In this study we showed that this discrepancy in the flow rate estimation is coherent to the theory of Womersley regarding spatial velocity profiles in unsteady flow conditions. In particular this work put in evidence that the error in flow rate estimation can be reduced by using the estimation formula recently proposed by Ponzini et al. [Ponzini R, Vergara C, Redaelli A, Veneziani A. Reliable CFD-based estimation of flow rate in haemodynamics measures. Ultrasound Med Biol 2006;32(10):1545-55], accounting for the unsteady nature of blood, applicable in the clinical practice without resorting to further measurements.

Osmometry is an essential technique for solution analysis and the investigation of chemical and biological phenomena. Commercially available osmometers rely on the measurements of freezing point, vapor pressure, and osmotic pressure of solutions. Although vapor pressure osmometry (VPO) and freezing point osmometry (FPO) can perform rapid and inexpensive measurements, they are indirect techniques, which rely on thermodynamic assumptions, which limit their applicability. While membrane osmometry (MO) provides a potentially unlimited direct measurement of osmotic pressure and solution osmolality, the conventional technique is often time-consuming and difficult to operate. In the present work, a novel membrane osmometer is presented. The instrument significantly reduces the conventional MO measurement time and is not subject to the limitations of VPO and FPO. For this paper, the osmotic pressure of aqueous sucrose solutions was collected in a molality range 0-5.5, by way of demonstration of the new instrument. When compared with data found in the literature, the experimental data were generally in good agreement. However, differences among results from the three techniques were observed.

Background: Functional mitral regurgitation (FMR) is the inability of mitral leaflets to coapt due to a combination of functional and geometrical factors. Valve competence is commonly restored by undersized annuloplasty, reducing the native annulus anteroposterior dimension. In case of severe FMR, this solution may be inadequate. The use of rings specific for the correction of FMR may lead to better results., Methods: The performance of the Geoform ring, a recently designed FMR-specific prosthesis, was compared with that of a standard Physio annuloplasty ring. Finite element modeling was used to simulate dilated cardiomyopathy-related FMR and compare, at the systolic peak, the valve's pathologic condition with the postoperative scenario corresponding to both devices. Three degrees of the pathology were simulated by progressively displacing papillary muscles apically, up to 5 mm. Three ring sizes were modeled., Results: Regurgitant area, coaptation length, and stresses acting on valve structures were assessed. When the use of the Geoform was modeled, coaptation length was always longer than 7 mm. In the most unfavorable case, the regurgitant area reduction was 74% with respect to baseline, and leaflets stresses were reduced by 20% when undersizing was simulated. When Physio ring implantation was simulated, coaptation length maximum extent was equal to 4.3 mm, the maximum regurgitant area reduction was equal to 60%, and leaflet stress reduction was observed., Conclusions: Disease-specific prostheses may allow for restoration of valve competence even for significant degrees of leaflets tethering and avoid the need for aggressive undersizing, thus leading to more durable results.

Objective: Automatic devices have been recently introduced to make the anastomosis procedure quick and efficient when creating a coronary bypass on the beating heart. However, the implantation of these devices could modify the graft configuration, consistently affecting the hemodynamics usually found in the traditional anastomosis. As local fluid dynamics could play a significant role in the onset of vessel wall pathologies, in this article a computational approach was designed to investigate flow patterns in the presence of the Ventrica magnetic vascular positioner (Ventrica MVP) device., Methods: A model of standard hand-sewn anastomosis and of automated magnetic anastomosis were constructed, and the finite volume method was used to simulate in silico realistic graft hemodynamics. Synthetic analytical descriptors -- i.e., time-averaged wall shear stress (TAWSS), oscillating shear index (OSI) and helical flow index (HFI) -- were calculated and compared for quantitative assessment of the anastomosis geometry hemodynamic performance., Results: In this case study, the same most critical region was identified for the 2 models as the one with the lowest TAWSS and the highest OSI (TAWSS=0.229, OSI=0.255 for the hand-sewn anastomosis; TAWSS=0.297, OSI=0.171 for the Ventrica MVP(R)). However, the shape of the Ventrica MVP does not induce more critical wall shear stresses, oscillating flow and damped helicity in the graft fluid dynamics, as compared with conventional anastomosis., Conclusions: We found that the use of the Ventrica MVP for the case study under investigation was not associated with more critical fluid dynamics than with conventional hand-sewn anastomosis. Thereby, the device could facilitate beating heart and minimally invasive coronary artery bypass grafting without increasing local hemodynamic-related risks of failure.