Towards real-time simulation of cardiac electrophysiology in a human heart at high resolution.

We have developed the capability to rapidly simulate cardiac electrophysiological phenomena in a human heart discretised at a resolution comparable with the length of a cardiac myocyte. Previous scientific investigation has generally invoked simplified geometries or coarse-resolution hearts, with simulation duration limited to 10s of heartbeats. Using state-of-the-art high-performance computing techniques coupled with one of the most powerful computers available (the 20 PFlop/s IBM BlueGene/Q at Lawrence Livermore National Laboratory), high-resolution simulation of the human heart can now be carried out over 1200 times faster compared with published results in the field.

In silico cardiac risk assessment in patients with long QT syndrome: type 1: clinical predictability of cardiac models.

Objectives: The study was designed to assess the ability of computer-simulated electrocardiography parameters to predict clinical outcomes and to risk-stratify patients with long QT syndrome type 1 (LQT1).

Background: Although attempts have been made to correlate mutation-specific ion channel dysfunction with patient phenotype in long QT syndrome, these have been largely unsuccessful. Systems-level computational models can be used to predict consequences of complex changes in channel function to the overall heart rhythm.

Cardiac electromechanical models: from cell to organ.

The heart is a multiphysics and multiscale system that has driven the development of the most sophisticated mathematical models at the frontiers of computational physiology and medicine. This review focuses on electromechanical (EM) models of the heart from the molecular level of myofilaments to anatomical models of the organ. Because of the coupling in terms of function and emergent behaviors at each level of biological hierarchy, separation of behaviors at a given scale is difficult.

Integrative modeling of the cardiac ventricular myocyte.

Cardiac electrophysiology is a discipline with a rich 50-year history of experimental research coupled with integrative modeling which has enabled us to achieve a quantitative understanding of the relationships between molecular function and the integrated behavior of the cardiac myocyte in health and disease. In this paper, we review the development of integrative computational models of the cardiac myocyte. We begin with a historical overview of key cardiac cell models that helped shape the field.

Computational modeling of cardiac disease: potential for personalized medicine.

Cardiovascular diseases are leading causes of death, reduce life quality and consume almost half a trillion dollars in healthcare expenses in the USA alone. Cardiac modeling and simulation technologies hold promise as important tools to improve cardiac care and are already in use to elucidate the fundamental mechanisms of cardiac physiology and pathophysiology. However, the emphasis has been on simulating average or exemplar cases.

Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations.

We develop a point model of the cardiac myofilament (MF) to simulate a wide variety of experimental muscle characterizations including Force-Ca relations and twitches under isometric, isosarcometric, isotonic, and auxotonic conditions. Complex MF behaviors are difficult to model because spatial interactions cannot be directly implemented as ordinary differential equations. We therefore allow phenomenological approximations with careful consideration to the relationships with the underlying biophysical mechanisms.

Spatially-compressed cardiac myofilament models generate hysteresis that is not found in real muscle.

In the field of cardiac modeling, calcium- (Ca-) based activation is often described by sets of ordinary differential equations that do not explicitly represent spatial interactions of regulatory proteins or crossbridge attachment. These spatially compressed models are most often mean-field representations as opposed to methods that explicitly compute the surrounding field (or equivalently, the surrounding environment) of individual regulatory units and crossbridges. Instead, a mean value is used to represent the whole population.

A plausible model for the digital response of p53 to DNA damage.

Recent observations show that the single-cell response of p53 to ionizing radiation (IR) is "digital" in that it is the number of oscillations rather than the amplitude of p53 that shows dependence on the radiation dose. We present a model of this phenomenon. In our model, double-strand break (DSB) sites induced by IR interact with a limiting pool of DNA repair proteins, forming DSB-protein complexes at DNA damage foci.

Reconstructing biological networks using conditional correlation analysis.

Motivation: One of the present challenges in biological research is the organization of the data originating from high-throughput technologies. One way in which this information can be organized is in the form of networks of influences, physical or statistical, between cellular components. We propose an experimental method for probing biological networks, analyzing the resulting data and reconstructing the network architecture.

Approaches to modeling crossbridges and calcium-dependent activation in cardiac muscle.

While the primary function of the heart is a pump, ironically, the development of myofilament models that predict developed force have generally lagged behind the modeling of the electrophysiological and Ca2+-handling aspects of heart cells. A major impediment is that the basic events in force generating actin-myosin interactions are still not well understood and quantified despite advanced techniques that can probe molecular levels events and identify numerous energetic states. As a result, the modeler must decide how to best abstract the many identified states into useful models with an essential tradeoff in the level of complexity.

Ising model of cardiac thin filament activation with nearest-neighbor cooperative interactions.

We have developed a model of cardiac thin filament activation using an Ising model approach from equilibrium statistical physics. This model explicitly represents nearest-neighbor interactions between 26 troponin/tropomyosin units along a one-dimensional array that represents the cardiac thin filament. With transition rates chosen to match experimental data, the results show that the resulting force-pCa (F-pCa) relations are similar to Hill functions with asymmetries, as seen in experimental data.