Physiological control of an in-series connected pulsatile VAD: numerical simulation study

This paper investigates ventricular assist device (VAD)-assisted cardiovascular dynamics under proportion–integration–differentiation (PID) feedback control. Previously, we have studied the cardiovascular responses under the support of an in-series connected reciprocating-valve VAD through numerical simulation, and no feedback control was applied in the VAD. In this research, we explore the contribution of the VAD control on the circulatory dynamics assisted by the reciprocating-valve VAD, in response to the changing physiological conditions. The classical PID control algorithm is implemented to regulate the VAD stroke beat-to-beat, based on the error signal between the expected and the realistic mean aortic pressures. Simulation results show that under the PID VAD control, physiological variables such as left atrial, ventricular and systemic arterial pressures, cardiac output and ventricular volumes are satisfactorily maintained in the physiological ranges. With the online PID feedback control, operation of the reciprocating-valve VAD can be satisfactorily regulated to accommodate metabolic requirements under various physiological conditions including normal resting and exercise situations.     

Numerical simulation of cardiovascular dynamics with different types of VAD assistance

A variety of methods by which mechanical circulatory support (MCS) can be provided have been described. However, the haemodynamic benefits of the different methods have not been adequately quantified. The aim of this paper is to compare the haemodynamic effects of six forms of MCS by numerical simulation. Three types of ventricular assist device (VAD) are studied: positive displacement; impeller and a novel reciprocating-valve design. Similarly, three pumping modes are modelled: constant flow; counterpulsation and copulsation. The cardiovascular system is modelled using an approach developed previously, using the concentrated parameter method by considering flow resistance, vessel elasticity and inertial effects of blood in individual conduit segments. The dynamic modelling of displacement and impeller pumps is represented by VAD inlet/outlet flow-rate changes. The dynamics of the reciprocating-valve pump is modelled with a specified displacement profile. Results show that in each simulation, the physiological variables of mean arterial pressure and systemic flow are adequately maintained. Modulation of the impeller pump flow profile produces a small (5 mmHg) oscillatory component to arterial pressure, whereas the displacement and reciprocating-valve pumps generate substantial arterial pressure and flow pulsatility. The impeller pump requires the least power input, the reciprocating valve pump slightly more, and the displacement pump the most. The in parallel configuration of the impeller and displacement pump designs with respect to the left ventricle provides near complete unloading and can cause the aortic valve to remain closed throughout the entire cardiac cycle with the attendant risk of aortic valve leaflet fusion following prolonged support. The in series configuration of the reciprocating-valve pump avoids this shortcoming but activation must be carefully synchronized to the cardiac cycle to allow adequate coronary perfusion. The reciprocating-valve pump is associated with haemodynamic advantages and a favourable power consumption.     

Vitamin A Deficiency Impairs Adaptive B and T Cell Responses to a Prototype Monovalent Attenuated Human Rotavirus Vaccine and Virulent Human Rotavirus Challenge in a Gnotobiotic Piglet Model

Rotaviruses (RV) are a major cause of gastroenteritis in children. Widespread vitamin A deficiency is associated with reduced efficacy of vaccines and higher incidence of diarrheal infections in children in developing countries. We established a vitamin A deficient (VAD) gnotobiotic piglet model that mimics subclinical vitamin A deficiency in children to study its effects on an oral human rotavirus (HRV) vaccine and virulent HRV challenge. Piglets derived from VAD and vitamin A sufficient (VAS) sows were orally vaccinated with attenuated HRV or mock, with/without supplemental vitamin A and challenged with virulent HRV. Unvaccinated VAD control piglets had significantly lower hepatic vitamin A, higher severity and duration of diarrhea and HRV fecal shedding post-challenge as compared to VAS control pigs. Reduced protection coincided with significantly higher innate (IFNα) cytokine and CD8 T cell frequencies in the blood and intestinal tissues, higher pro-inflammatory (IL12) and 2-3 fold lower anti-inflammatory (IL10) cytokines, in VAD compared to VAS control pigs. Vaccinated VAD pigs had higher diarrhea severity scores compared to vaccinated VAS pigs, which coincided with lower serum IgA HRV antibody titers and significantly lower intestinal IgA antibody secreting cells post-challenge in the former groups suggesting lower anamnestic responses. A trend for higher serum HRV IgG antibodies was observed in VAD vs VAS vaccinated groups post-challenge. The vaccinated VAD (non-vitamin A supplemented) pigs had significantly higher serum IL12 (PID2) and IFNγ (PID6) compared to vaccinated VAS groups suggesting higher Th1 responses in VAD conditions. Furthermore, regulatory T-cell responses were compromised in VAD pigs. Supplemental vitamin A in VAD pigs did not fully restore the dysregulated immune responses to AttHRV vaccine or moderate virulent HRV diarrhea. Our findings suggest that that VAD in children in developing countries may partially contribute to more severe rotavirus infection and lower HRV vaccine efficacy.     

Transforming growth factor β2 is negatively regulated by endogenous retinoic acid during early heart morphogenesis

Vitamin A‐deficient (VAD) quail embryos lack the vitamin A‐active form, retinoic acid (RA) and are characterized by a phenotype that includes a grossly abnormal cardiovascular system that can be rescued by RA. Here we report that the transforming growth factor, TGFβ2 is involved in RA‐regulated cardiovascular development. In VAD embryos TGFβ2 mRNA and protein expression are greatly elevated. The expression of TGFβ receptor II is also elevated in VAD embryos but is normalized by treatment with TGFβ2‐specific antisense oligonucleotides (AS). Administration of this AS or an antibody specific for TGFβ2 to VAD embryos normalizes posterior heart development and vascularization, while the administration of exogenous active TGFβ2 protein to normal quail embryos mimics the excessive TGFβ2 status of VAD embryos and induces VAD cardiovascular phenotype. In VAD embryos pSmad2/3 and pErk1 are not activated, while pErk2 and pcRaf are elevated and pSmad1/5/8 is diminished. We conclude that in the early avian embryo TGFβ2 has a major role in the retinoic acid‐regulated posterior heart morphogenesis for which it does not use Smad2/3 pathways, but may use other signaling pathways. Importantly, we conclude that retinoic acid is a critical negative physiological regulator of the magnitude of TGFβ2 signals during vertebrate heart formation.     

Retinoic acid is a negative physiological regulator of N‐cadherin during early avian heart morphogenesis

The vitamin A‐deficient (VAD) early avian embryo has a grossly abnormal cardiovascular system that is rescued by treating the embryo with the vitamin A‐active form, retinoic acid (RA). Here we examine the role of N‐cadherin (N‐cad) in RA‐regulated early cardiovascular morphogenesis. N‐cad mRNA and protein are expressed globally in the presomite through HH14 normal and VAD quail embryos. The expression in VAD embryos prior to HH10 is significantly higher than that in normal embryos. Functional analyses of the N‐cad overproducing VAD embryos reveal N‐cad involvement in the RA‐regulated cardiovascular development and suggest that N‐cad expression may be mediated by Msx1. We provide evidence that in the early avian embryo, endogenous RA is a negative physiological regulator of N‐cad. We hypothesize that a critical endogenous level of N‐cad is needed for normal early cardiovascular morphogenesis to occur and that this level is ensured by stage‐specific, developmentally regulated RA signaling.     

Analysis of “integrate-to-threshold” neural coding schemes

Methods of analysis for some deterministic and stochastic variants of the integrate-to-threshold neural coding scheme are presented. Adaptation phenomena are modeled by means of feedforward and feedback adaptive threshold control. Simulations of sinusoidal and step responses reproduce satisfactorily the qualitative characteristics of adaptation as compared with physiological data. It is postulated that such adaptive threshold control may be accomplished by the release, or conformation change, of molecules involved in the control of excitable-channel dynamics.     

一套研究机械电反馈的心室压力钳系统

在心脏机械电反馈的研究中准确控制机械刺激是非常重要的。本研究室构建了一套适用于离体家兔心脏的心室压力钳系统。该系统通过计算机控制压力钳,不仅能模拟正常生理条件下左心室的压力波形,还能在心室活动周期的特定时相、以适当波形对心室施加机械刺激。该系统集心脏灌流与起搏、表面心电图记录、单相动作电位记录、心室压力钳制与测定等多种功能于一体,特别适用于器官水平上观察机械电反馈现象并探讨其机制。     

Is severe pulmonary hypertension a contraindication for orthotopic heart transplantation? Not any more

Pulmonary hypertension (PH) unresponsive to pharmacological intervention is considered a contraindication for orthotopic heart transplantation (OHTX) due to risk of postoperative right-heart failure. In this prospective study, we describe our experience with a treatment strategy of improving severe PH in heart transplant candidates by means of ventricular assist device (VAD) implantation and subsequent OHTX. In 11 heart transplantation candidates with severe PH unresponsive to pharmacological intervention we implanted VAD with the aim of achieving PH to values acceptable for OHTX. In all patients we observed significant drop in pulmonary pressures, PVR and TPG (p < 0.001 for all) 3 months after VAD implantation to values sufficient to allow OHTX. Seven patients underwent transplantation (mean duration of support 216 days) while none of patients suffered right-side heart failure in postoperative period. Two patients died after transplantation and five patients are living in very good condition with a mean duration of 286 days after OHTX. In our opinion, severe PH is not a contraindication for orthotopic heart transplantation any more.     

An approximate stochastic optimal control framework to simulate nonlinear neuro-musculoskeletal models in the presence of noise

Optimal control simulations have shown that both musculoskeletal dynamics and physiological noise are important determinants of movement. However, due to the limited efficiency of available computational tools, deterministic simulations of movement focus on accurately modelling the musculoskeletal system while neglecting physiological noise, and stochastic simulations account for noise while simplifying the dynamics. We took advantage of recent approaches where stochastic optimal control problems are approximated using deterministic optimal control problems, which can be solved efficiently using direct collocation. We were thus able to extend predictions of stochastic optimal control as a theory of motor coordination to include muscle coordination and movement patterns emerging from non-linear musculoskeletal dynamics. In stochastic optimal control simulations of human standing balance, we demonstrated that the inclusion of muscle dynamics can predict muscle co-contraction as minimal effort strategy that complements sensorimotor feedback control in the presence of sensory noise. In simulations of reaching, we demonstrated that nonlinear multi-segment musculoskeletal dynamics enables complex perturbed and unperturbed reach trajectories under a variety of task conditions to be predicted. In both behaviors, we demonstrated how interactions between task constraint, sensory noise, and the intrinsic properties of muscle influence optimal muscle coordination patterns, including muscle co-contraction, and the resulting movement trajectories. Our approach enables a true minimum effort solution to be identified as task constraints, such as movement accuracy, can be explicitly imposed, rather than being approximated using penalty terms in the cost function. Our approximate stochastic optimal control framework predicts complex features, not captured by previous simulation approaches, providing a generalizable and valuable tool to study how musculoskeletal dynamics and physiological noise may alter neural control of movement in both healthy and pathological movements.     

The intraaortic balloon pump: a nonlinear digital computer model

A computer based model for in-series cardiac assistance by intraaortic balloon pumping was developed in this study. The model, obtained from the Navier-Stokes and Continuity equations, was capable of computing pressures, volumetric flow rates and radii through the arterial system. The model was used to study the effects of a wide range of assist device timing adjustments on the benefits of ventricular assistance under conditions corresponding to those measured during animal experiments. The model was also used to study the relationship between device timing adjustments and the benefits of ventricular assistance under constant cardiovascular state conditions. Such studies are important in isolating the response of the system to assist device phasing from the response associated with system state. The results obtained in this study demonstrate that the hemodynamic response of the cardiovascular system to intraaortic balloon pumping is a sensitive function of both the state of the cardiovascular system and phasing of the assist device.     

The coronary bed and its role in the cardiovascular system: a review and an introductory single-branch model

To investigate cardiovascular haemodynamics under normal and pathological conditions, a closed-loop model of the cardiovascular system already presented in the literature¹, has been complemented by a model of the coronary bed. Oxygen available to the myocardium is strictly related to the coronary blood flow; we have developed threshold criteria which correlate cardiac output with the coronary flow. The system utilizes control systems related to the cardiac contractility and frequency, and imitates feedback mechanisms peculiar to the heart. The work exemplifies the autoregulation of events that occur when the equilibrium of the system is disturbed. It is suggested that the heart plays an active role in trying to restore the haemodynamic parameters to their physiological values.     

Long-term vitamin A deficiency induces alteration of adult mouse spermatogenesis and spermatogonial differentiation: direct effect on spermatogonial gene expression and indirect effects via somatic cells

The objective of this study was to further understand the genetic mechanisms of vitamin A deficiency (VAD) induced arrest of spermatogonial stem-cell differentiation.Vitamin A and its derivatives (the retinoids) participate in many physiological processes including vision, cellular differentiation and reproduction. VAD affects spermatogenesis, the subject of our present study. Spermatogenesis is a highly regulated process of differentiation and complex morphologic alterations that leads to the formation of sperm in the seminiferous epithelium. VAD causes early cessation of spermatogenesis, characterized by degeneration of meiotic germ cells, leading to seminiferous tubules containing mostly type A spermatogonia and Sertoli cells. These observations led us to the hypothesis that VAD affects not only germ cells but also somatic cells.To investigate the effects of VAD on spermatogenesis in mice we used adult Balb/C mice fed with Control or VAD diet for an extended period of time (6–28 weeks). We first observed the chronology, then the extent of the effects of VAD on the testes. Using microarray analysis of isolated pure populations of spermatogonia, Leydig and Sertoli cells from control and VAD 18- and 25-week mice, we examined the effects of VAD on gene expression and identified target genes involved in the arrest of spermatogonial differentiation and spermatogenesis.Our results provide a more precise definition of the chronology and magnitude of the consequences of VAD on mouse testes than the previously available literature and highlight direct and indirect (via somatic cells) effects of VAD on germ cell differentiation.     

Stabilizing PID controllers for a single-link biomechanical model with position, velocity, and force feedback

In this paper we address the problem of PID stabilization of a single-link inverted pendulum-based biomechanical model with force feedback, two levels of position and velocity feedback, and with delays in all the feedback loops. The novelty of the proposed model lies in its physiological relevance, whereby both small and medium latency sensory feedbacks from muscle spindle (MS), and force feedback from Golgi tendon organ (GTO) are included in the formulation. The biomechanical model also includes active and passive viscoelastic feedback from Hill-type muscle model and a second-order low-pass function for muscle activation. The central nervous system (CNS) regulation of postural movement is represented by a proportional-integral-derivative (PID) controller. Padé approximation of delay terms is employed to arrive at an overall rational transfer function of the biomechanical model. The Hermite-Biehler theorem is then used to derive stability results, leading to the existence of stabilizing PID controllers. An algorithm for selection of stabilizing feedback gains is developed using the linear matrix inequality (LMI) approach.     

Nano-imaging of the beating mouse heart in vivo: Importance of sarcomere dynamics,as opposed to sarcomere length per se,in the regulation of cardiac function

Sarcomeric contraction in cardiomyocytes serves as the basis for the heart’s pump functions in mammals. Although it plays a critical role in the circulatory system, myocardial sarcomere length (SL) change has not been directly measured in vivo under physiological conditions because of technical difficulties. In this study, we developed a high speed (100–frames per second), high resolution (20-nm) imaging system for myocardial sarcomeres in living mice. Using this system, we conducted three-dimensional analysis of sarcomere dynamics in left ventricular myocytes during the cardiac cycle, simultaneously with electrocardiogram and left ventricular pressure measurements. We found that (a) the working range of SL was on the shorter end of the resting distribution, and (b) the left ventricular–developed pressure was positively correlated with the SL change between diastole and systole. The present findings provide the first direct evidence for the tight coupling of sarcomere dynamics and ventricular pump functions in the physiology of the heart.     

MRI/MRS assessment of in vivo murine cardiac metabolism, morphology, and function at physiological heart rates

Transgenic mice are increasingly used to probe genetic aspects of cardiovascular pathophysiology. However, the small size and rapid rates of murine hearts make noninvasive, physiological in vivo studies of cardiac bioenergetics and contractility difficult. The aim of this report was to develop an integrated, noninvasive means of studying in vivo murine cardiac metabolism, morphology, and function under physiological conditions by adapting and modifying noninvasive cardiac magnetic resonance imaging (MRI) with image-guided (31)P magnetic resonance spectroscopy techniques used in humans to mice. Using spatially localized, noninvasive (31)P nuclear magnetic resonance spectroscopy and MRI at 4.7 T, we observe mean murine in vivo myocardial phosphocreatine-to-ATP ratios of 2.0 +/- 0.2 and left ventricular ejection fractions of 65 +/- 7% at physiological heart rates ( approximately 600 beats/min). These values in the smallest species studied to date are similar to those reported in normal humans. Although these observations do not confirm a degree of metabolic scaling with body size proposed by prior predictions, they do suggest that mice can serve, at least at this level, as a model for human cardiovascular physiology. Thus it is now possible to noninvasively study in vivo myocardial bioenergetics, morphology, and contractile function in mice under physiological conditions.     

A finite element study relating to the rapid filling phase of the human ventricles

During the rapid diastolic filling phase at rest, the ventricles of the human heart double approximately in volume. In order to investigate whether the ventricular filling pressures measured under physiological conditions can give rise to such an extensive augmentation in ventricular volumes, a finite element model of the human right and left ventricles has been developed, taking into account the nonlinear mechanical behavior and effective compressibility of the myocardial tissue. The results were compared with the filling phase of the human left ventricle as extrapolated from measurements documented in the literature. We arrived at the conclusion that the ventricular pressures measured during the rapid filling phase cannot be the sole cause of the rise of the observed ventricular volumes. We rather advocate the assumption that further dilating mechanisms might be part of ventricular activity thus heralding a multiple function of the ventricular muscle body. A further result indicates that under normal conditions the influence of the viscoelasticity of the tissue should not be disregarded in ventricular mechanics.     

Low pulse pressure with high pulsatile external left ventricular power: influence of aortic waves

Elevated pulse pressure (pp) is considered to be a risk factor for adverse cardiovascular events since it is directly related to an elevated myocardial workload. Information about both pressure and flow wave must be provided to assess hemodynamic complexity and true level of external left ventricular power (ELVP). pp value as a single feature of aortic waves cannot identify true level of ELVP. However, it is generally presumed that ELVP (and consequently LV workload) is positively correlated with pp. This study examined this positive correlation. The aim of this study was to test the hypothesis that aortic wave dynamics can create destructive hemodynamic conditions that increase the ELVP even though pp appears to be normal. To test this hypothesis, a computational model of the aorta with physiological properties was used. A Finite Element Method with fluid-structure interaction was employed to solve the equations of the solid and fluid. The aortic wall was assumed to be elastic and isotropic. The blood was assumed to be an incompressible Newtonian fluid. Simulations were performed for various heart rates (HR) and different aortic compliances while keeping the shape of the inlet flow and peripheral resistance constant. As expected, in most of the cases studied here, higher pp was associated with higher LV power demand. However, for a given cardiac output, mean pressure, and location of total reflection site, we have found cases where the above-mentioned trend does not hold. Our results suggest that using pp as a single index can result in an underestimation of the LV power demand under certain conditions related to the altered wave dynamics. Hence, in hypertensive patients, a full analysis of aortic wave dynamics is essential for the prevention and management of left ventricular hypertrophy (LVH) and congestive heart failure.     

The design of controllers for batch bioreactors

The implementation of control algorithms to batch bioreactors is often complicated by variations in process dynamics that occur during the course of fermentation. Such a wide operating range often renders the performance of fixed gain proportional-integral-differential (PID) controllers unsatisfactory. In this work, detailed studies on the control of batch fermentations are per formed. Two simple controller designs are presented with the intent to compensate for changing process dynamics. One design incorporates the concepts of static feedforward-feedback control. While this technique produces tighter control than feedback alone, it is not as successful as a controller based on gain scheduling. The gain-scheduling controller, a subclass of adaptive controllers, uses the oxygen uptake rate as an auxiliary variable to fine-tune the PID controller parameters. The control of oxygen tension in the bioreactor is used as a vehicle to convey the proposed ideas, analyses, and results. Simulation experiments indicate significant improvement in controller performance can be achieved by both of the proposed approaches even in the presence of measurement noise.     

In vitro validation of finite-element model of AAA hemodynamics incorporating realistic outlet boundary conditions

The purpose of this study is to validate numerical simulations of flow and pressure in an abdominal aortic aneurysm (AAA) using phase-contrast magnetic resonance imaging (PCMRI) and an in vitro phantom under physiological flow and pressure conditions. We constructed a two-outlet physical flow phantom based on patient imaging data of an AAA and developed a physical Windkessel model to use as outlet boundary conditions. We then acquired PCMRI data in the phantom while it operated under conditions mimicking a resting and a light exercise physiological state. Next, we performed in silico numerical simulations and compared experimentally measured velocities, flows, and pressures in the in vitro phantom to those computed in the in silico simulations. There was a high degree of agreement in all of the pressure and flow waveform shapes and magnitudes between the experimental measurements and simulated results. The average pressures and flow split difference between experiment and simulation were all within 2%. Velocity patterns showed good agreement between experimental measurements and simulated results, especially in the case of whole-cycle averaged comparisons. We demonstrated methods to perform in vitro phantom experiments with physiological flows and pressures, showing good agreement between numerically simulated and experimentally measured velocity fields and pressure waveforms in a complex patient-specific AAA geometry.     

Superior performance of continuous over pulsatile flow ventricular assist devices in the single ventricle circulation: A computational study

This study compares the physiological responses of systemic-to-pulmonary shunted single ventricle patients to pulsatile and continuous flow ventricular assist devices (VADs). Performance differences between pulsatile and continuous flow VADs have been clinically observed, but the underlying mechanism remains poorly understood. Six systemic-to-pulmonary shunted single ventricle patients (mean BSA=0.30 m²) were computationally simulated using a lumped-parameter network tuned to match patient specific clinical data. A first set of simulations compared current clinical implementation of VADs in single ventricle patients. A second set modified pulsatile flow VAD settings with the goal to optimize cardiac output (CO). For all patients, the best-case continuous flow VAD CO was at least 0.99 L/min greater than the optimized pulsatile flow VAD CO (p=0.001). The 25 and 50 mL pulsatile flow VADs exhibited incomplete filling at higher heart rates that reduced CO as much as 9.7% and 37.3% below expectations respectively. Optimization of pulsatile flow VAD settings did not achieve statistically significant (p<0.05) improvement to CO. Results corroborate clinical experience that continuous flow VADs produce higher CO and superior ventricular unloading in single ventricle patients. Impaired filling leads to performance degradation of pulsatile flow VADs in the single ventricle circulation.