A greedy classifier optimization strategy to assess ion channel blocking activity and pro-arrhythmia in hiPSC-cardiomyocytes.

Novel studies conducting cardiac safety assessment using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are promising but might be limited by their specificity and predictivity. It is often challenging to correctly classify ion channel blockers or to sufficiently predict the risk for Torsade de Pointes (TdP). In this study, we developed a method combining in vitro and in silico experiments to improve machine learning approaches in delivering fast and reliable prediction of drug-induced ion-channel blockade and proarrhythmic behaviour. The algorithm is based on the construction of a dictionary and a greedy optimization, leading to the definition of optimal classifiers. Finally, we present a numerical tool that can accurately predict compound-induced pro-arrhythmic risk and involvement of sodium, calcium and potassium channels, based on hiPSC-CM field potential data.

Reducing Hypermuscularization of the Transitional Segment Between Arterioles and Capillaries Protects Against Spontaneous Intracerebral Hemorrhage.

BACKGROUND: Spontaneous deep intracerebral hemorrhage (ICH) is a devastating subtype of stroke without specific treatments. It has been thought that smooth muscle cell (SMC) degeneration at the site of arteriolar wall rupture may be sufficient to cause hemorrhage. However, deep ICHs are rare in some aggressive small vessel diseases that are characterized by significant arteriolar SMC degeneration. Here we hypothesized that a second cellular defect may be required for the occurrence of ICH. METHODS: We studied a genetic model of spontaneous deep ICH using Col4a1+/G498V and Col4a1+/G1064D mouse lines that are mutated for the α1 chain of collagen type IV. We analyzed cerebroretinal microvessels, performed genetic rescue experiments, vascular reactivity analysis, and computational modeling. We examined postmortem brain tissues from patients with sporadic deep ICH. RESULTS: We identified in the normal cerebroretinal vasculature a novel segment between arterioles and capillaries, herein called the transitional segment (TS), which is covered by mural cells distinct from SMCs and pericytes. In Col4a1 mutant mice, this TS was hypermuscularized, with a hyperplasia of mural cells expressing more contractile proteins, whereas the upstream arteriole exhibited a loss of SMCs. TSs mechanistically showed a transient increase in proliferation of mural cells during postnatal maturation. Mutant brain microvessels, unlike mutant arteries, displayed a significant upregulation of SM genes and Notch3 target genes, and genetic reduction of Notch3 in Col4a1+/G498V mice protected against ICH. Retina analysis showed that hypermuscularization of the TS was attenuated, but arteriolar SMC loss was unchanged in Col4a1+/G498V, Notch3+/- mice. Moreover, hypermuscularization of the retinal TS increased its contractility and tone and raised the intravascular pressure in the upstream feeding arteriole. We similarly found hypermuscularization of the TS and focal arteriolar SMC loss in brain tissues from patients with sporadic deep ICH. CONCLUSIONS: Our results suggest that hypermuscularization of the TS, through increased Notch3 activity, is involved in the occurrence of ICH in Col4a1 mutant mice, by raising the intravascular pressure in the upstream feeding arteriole and promoting its rupture at the site of SMC loss. Our human data indicate that these 2 mutually reinforcing vascular defects may represent a general mechanism of deep ICH.

Augmented resistive immersed surfaces valve model for the simulation of cardiac hemodynamics with isovolumetric phases.

In order to reduce the complexity of heart hemodynamics simulations, uncoupling approaches are often considered for the modeling of the immersed valves as an alternative to complex fluid-structure interaction (FSI) models. A possible shortcoming of these simplified approaches is the difficulty to correctly capture the pressure dynamics during the isovolumetric phases. In this work, we propose an enhanced resistive immersed surfaces (RIS) model of cardiac valves, which overcomes this issue. The benefits of the model are investigated and tested in blood flow simulations of the left heart where the physiological behavior of the intracavity pressure during the isovolumetric phases is recovered without using fully coupled fluid-structure models and without important alteration of the associated velocity field.

How to choose biomarkers in view of parameter estimation.

In numerous applications in biophysics, physiology and medicine, the system of interest is studied by monitoring quantities, called biomarkers, extracted from measurements. These biomarkers convey some information about relevant hidden quantities, which can be seen as parameters of an underlying model. In this paper we propose a strategy to automatically design biomarkers to estimate a given parameter. Such biomarkers are chosen as the solution of a sparse optimization problem given a user-supplied dictionary of candidate features. The method is in particular illustrated with two realistic applications, one in electrophysiology and the other in hemodynamics. In both cases, our algorithm provides composite biomarkers which improve the parameter estimation problem.

In silico assessment of the effects of various compounds in MEA/hiPSC-CM assays: Modeling and numerical simulations.

We propose a mathematical approach for the analysis of drugs effects on the electrical activity of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) based on multi-electrode array (MEA) experiments. Our goal is to produce an in silico tool able to simulate drugs action in MEA/hiPSC-CM assays. The mathematical model takes into account the geometry of the MEA and the electrodes' properties. The electrical activity of the stem cells at the ion-channel level is governed by a system of ordinary differential equations (ODEs). The ODEs are coupled to the bidomain equations, describing the propagation of the electrical wave in the stem cells preparation. The field potential (FP) measured by the MEA is modeled by the extracellular potential of the bidomain equations. First, we propose a strategy allowing us to generate a field potential in good agreement with the experimental data. We show that we are able to reproduce realistic field potentials by introducing different scenarios of heterogeneity in the action potential. This heterogeneity reflects the differentiation atria/ventricles and the age of the cells. Second, we introduce a drug/ion channels interaction based on a pore block model. We conduct different simulations for five drugs (mexiletine, dofetilide, bepridil, ivabradine and BayK). We compare the simulation results with the field potential collected from experimental measurements. Different biomarkers computed on the FP are considered, including depolarization amplitude, repolarization delay, repolarization amplitude and depolarization-repolarization segment. The simulation results show that the model reflect properly the main effects of these drugs on the FP.

Identification of Ion Currents Components Generating Field Potential Recorded in MEA From hiPSC-CM.

OBJECTIVE: Multi electrodes arrays (MEAs) combined with cardiomyocytes derived from human-induced pluripotent stem cells (hiPSC-CMs) can enable high- or medium-throughput drug screening in safety pharmacology. This technology has recently attracted a lot of attention, in particular from an international initiative named CiPA. But it is currently limited by the difficulty to analyze the measured signals. We propose a strategy to analyze the signals acquired by the MEA and to automatically deduce the channels affected by the drug. METHODS: Our method is based on the bidomain equations, a model for the MEA electrodes, and an inverse problem strategy. RESULTS: in silico MEA signals are obtained for two commercial devices and an example of early after depolarization is presented. Then, by processing real signals obtained for four different compounds, our algorithm was able to provide dose-response curves for potassium, sodium, and calcium channels. For ivabradine and moxifloxacin, the IC50 and dose-response curves are in very good agreement with known values. SIGNIFICANCE: The proposed strategy offers a possible answer to a major question raised by the community of safety pharmacology. By allowing a more automated analysis of the signals, our approach could contribute to promote the technology based on MEA and hiPSC-CMs and, therefore, improve reliability and efficiency of drug screening.

Modelling variability in cardiac electrophysiology: a moment-matching approach.

The variability observed in action potential (AP) cardiomyocyte measurements is the consequence of many different sources of randomness. Often ignored, this variability may be studied to gain insight into the cell ionic properties. In this paper, we focus on the study of ionic channel conductances and describe a methodology to estimate their probability density function (PDF) from AP recordings. The method relies on the matching of observable statistical moments and on the maximum entropy principle. We present four case studies using synthetic and sets of experimental AP measurements from human and canine cardiomyocytes. In each case, the proposed methodology is applied to infer the PDF of key conductances from the exhibited variability. The estimated PDFs are discussed and, when possible, compared to the true distributions. We conclude that it is possible to extract relevant information from the variability in AP measurements and discuss the limitations and possible implications of the proposed approach.

Composite Biomarkers Derived from Micro-Electrode Array Measurements and Computer Simulations Improve the Classification of Drug-Induced Channel Block.

The Micro-Electrode Array (MEA) device enables high-throughput electrophysiology measurements that are less labor-intensive than patch-clamp based techniques. Combined with human-induced pluripotent stem cells cardiomyocytes (hiPSC-CM), it represents a new and promising paradigm for automated and accurate in vitro drug safety evaluation. In this article, the following question is addressed: which features of the MEA signals should be measured to better classify the effects of drugs? A framework for the classification of drugs using MEA measurements is proposed. The classification is based on the ion channels blockades induced by the drugs. It relies on an in silico electrophysiology model of the MEA, a feature selection algorithm and automatic classification tools. An in silico model of the MEA is developed and is used to generate synthetic measurements. An algorithm that extracts MEA measurements features designed to perform well in a classification context is described. These features are called composite biomarkers. A state-of-the-art machine learning program is used to carry out the classification of drugs using experimental MEA measurements. The experiments are carried out using five different drugs: mexiletine, flecainide, diltiazem, moxifloxacin, and dofetilide. We show that the composite biomarkers outperform the classical ones in different classification scenarios. We show that using both synthetic and experimental MEA measurements improves the robustness of the composite biomarkers and that the classification scores are increased.

Comprehensive in vitro Proarrhythmia Assay (CiPA): Pending issues for successful validation and implementation.

INTRODUCTION: The Comprehensive in vitro Proarrhythmia Assay (CiPA) is a nonclinical Safety Pharmacology paradigm for discovering electrophysiological mechanisms that are likely to confer proarrhythmic liability to drug candidates intended for human use. TOPICS COVERED: Key talks delivered at the 'CiPA on my mind' session, held during the 2015 Annual Meeting of the Safety Pharmacology Society (SPS), are summarized. Issues and potential solutions relating to crucial constituents [e.g., biological materials (ion channels and pluripotent stem cell-derived cardiomyocytes), study platforms, drug solutions, and data analysis] of CiPA core assays are critically examined. DISCUSSION: In order to advance the CiPA paradigm from the current testing and validation stages to a research and regulatory drug development strategy, systematic guidance by CiPA stakeholders is necessary to expedite solutions to pending and newly arising issues. Once a study protocol is proved to yield robust and reproducible results within and across laboratories, it can be implemented as qualified regulatory procedure.

Numerical simulation of electrocardiograms for full cardiac cycles in healthy and pathological conditions.

This work is dedicated to the simulation of full cycles of the electrical activity of the heart and the corresponding body surface potential. The model is based on a realistic torso and heart anatomy, including ventricles and atria. One of the specificities of our approach is to model the atria as a surface, which is the kind of data typically provided by medical imaging for thin volumes. The bidomain equations are considered in their usual formulation in the ventricles, and in a surface formulation on the atria. Two ionic models are used: the Courtemanche-Ramirez-Nattel model on the atria and the 'minimal model for human ventricular action potentials' by Bueno-Orovio, Cherry, and Fenton in the ventricles. The heart is weakly coupled to the torso by a Robin boundary condition based on a resistor-capacitor transmission condition. Various electrocardiograms (ECGs) are simulated in healthy and pathological conditions (left and right bundle branch blocks, Bachmann's bundle block, and Wolff-Parkinson-White syndrome). To assess the numerical ECGs, we use several qualitative and quantitative criteria found in the medical literature. Our simulator can also be used to generate the signals measured by a vest of electrodes. This capability is illustrated at the end of the article. Copyright © 2015 John Wiley & Sons, Ltd.

Identification of artery wall stiffness: in vitro validation and in vivo results of a data assimilation procedure applied to a 3D fluid-structure interaction model.

We consider the problem of estimating the stiffness of an artery wall using a data assimilation method applied to a 3D fluid-structure interaction (FSI) model. Recalling previous works, we briefly present the FSI model, the data assimilation procedure and the segmentation algorithm. We present then two examples of the procedure using real data. First, we estimate the stiffness distribution of a silicon rubber tube from image data. Second, we present the estimation of aortic wall stiffness from real clinical data.

Group-wise construction of reduced models for understanding and characterization of pulmonary blood flows from medical images.

3D computational fluid dynamics (CFD) in patient-specific geometries provides complementary insights to clinical imaging, to better understand how heart disease, and the side effects of treating heart disease, affect and are affected by hemodynamics. This information can be useful in treatment planning for designing artificial devices that are subject to stress and pressure from blood flow. Yet, these simulations remain relatively costly within a clinical context. The aim of this work is to reduce the complexity of patient-specific simulations by combining image analysis, computational fluid dynamics and model order reduction techniques. The proposed method makes use of a reference geometry estimated as an average of the population, within an efficient statistical framework based on the currents representation of shapes. Snapshots of blood flow simulations performed in the reference geometry are used to build a POD (Proper Orthogonal Decomposition) basis, which can then be mapped on new patients to perform reduced order blood flow simulations with patient specific boundary conditions. This approach is applied to a data-set of 17 tetralogy of Fallot patients to simulate blood flow through the pulmonary artery under normal (healthy or synthetic valves with almost no backflow) and pathological (leaky or absent valve with backflow) conditions to better understand the impact of regurgitated blood on pressure and velocity at the outflow tracts. The model reduction approach is further tested by performing patient simulations under exercise and varying degrees of pathophysiological conditions based on reduction of reference solutions (rest and medium backflow conditions respectively).

A robust and efficient valve model based on resistive immersed surfaces.

A procedure for modeling the heart valves is presented. Instead of modeling complete leaflet motion, leaflets are modeled in open and closed configurations. The geometry of each configuration can be defined, for example, from in vivo image data. This method enables significant computational savings compared with complete fluid-structure interaction and contact modeling, while maintaining realistic three-dimensional velocity and pressure distributions near the valve, which is not possible from lumped parameter modeling. Leaflets are modeled as immersed, fixed surfaces over which a resistance to flow is assigned. On the basis of local flow conditions, the resistance values assigned for each configuration are changed to switch the valve between open and closed states. This formulation allows for the pressure to be discontinuous across the valve. To illustrate the versatility of the model, realistic and patient-specific simulations are presented, as well as comparison with complete fluid-structure interaction simulation.

MagnetoHemoDynamics in the aorta and electrocardiograms.

This paper addresses a complex multi-physical phenomenon involving cardiac electrophysiology and hemodynamics. The purpose is to model and simulate a phenomenon that has been observed in magnetic resonance imaging machines: in the presence of a strong magnetic field, the T-wave of the electrocardiogram (ECG) gets bigger, which may perturb ECG-gated imaging. This is due to a magnetohydrodynamic (MHD) effect occurring in the aorta. We reproduce this experimental observation through computer simulations on a realistic anatomy, and with a three-compartment model: inductionless MHD equations in the aorta, bi-domain equations in the heart and electrical diffusion in the rest of the body. These compartments are strongly coupled and solved using finite elements. Several benchmark tests are proposed to assess the numerical solutions and the validity of some modeling assumptions. Then, ECGs are simulated for a wide range of magnetic field intensities (from 0 to 20 T).

Sequential parameter estimation for fluid-structure problems: application to hemodynamics.

We present a robust and computationally efficient parameter estimation strategy for fluid-structure interaction problems. The method is based on a filtering algorithm restricted to the parameter space, known as the reduced-order unscented Kalman filter. It does not require any adjoint or tangent problems. In addition, it can easily be run in parallel, which is of great interest in fluid-structure problems where the computational cost of the forward simulation is already a challenge in itself. We illustrate our methodology with the estimation of the artery wall stiffness from the wall displacement measurements - as they could be extracted from medical imaging - in a three-dimensional idealized abdominal aortic aneurysm. We also show preliminary results about the estimation of the proximal Windkessel resistance, which is an important parameter for setting appropriate fluid boundary conditions.

Computational analysis of an aortic valve jet with Lagrangian coherent structures.

Important progress has been achieved in recent years in simulating the fluid-structure interaction around cardiac valves. An important step in making these computational tools useful to clinical practice is the development of postprocessing techniques to extract clinically relevant information from these simulations. This work focuses on flow through the aortic valve and illustrates how the computation of Lagrangian coherent structures can be used to improve insight into the transport mechanics of the flow downstream of the valve, toward the goal of aiding clinical decision making and the understanding of pathophysiology.

Mathematical modeling of electrocardiograms: a numerical study.

This paper deals with the numerical simulation of electrocardiograms (ECG). Our aim is to devise a mathematical model, based on partial differential equations, which is able to provide realistic 12-lead ECGs. The main ingredients of this model are classical: the bidomain equations coupled to a phenomenological ionic model in the heart, and a generalized Laplace equation in the torso. The obtention of realistic ECGs relies on other important features--including heart-torso transmission conditions, anisotropy, cell heterogeneity and His bundle modeling--that are discussed in detail. The numerical implementation is based on state-of-the-art numerical methods: domain decomposition techniques and second order semi-implicit time marching schemes, offering a good compromise between accuracy, stability and efficiency. The numerical ECGs obtained with this approach show correct amplitudes, shapes and polarities, in all the 12 standard leads. The relevance of every modeling choice is carefully discussed and the numerical ECG sensitivity to the model parameters investigated.

[Numerical simulation of the cardiovascular system]. / Simulation numérique du système cardiovasculaire.

In this article, we aim at giving a non-technical overview of some mathematical models currently used in the numerical simulation of the cardiovascular system. A hierarchy of models for blood flows in large arteries is briefly described as well as an electromechanical model for the heart. We discuss some possible applications of the numerical simulations of such models, for example the optimization of prostheses. We also address the issue of the data assimilation for the calibration of the models.