BeatProfiler: Multimodal In Vitro Analysis of Cardiac Function Enables Machine Learning Classification of Diseases and Drugs

Goal: Contractile response and calcium handling are central to understanding cardiac function and physiology, yet existing methods of analysis to quantify these metrics are often time-consuming, prone to mistakes, or require specialized equipment/license. We developed BeatProfiler, a suite of cardiac analysis tools designed to quantify contractile function, calcium handling, and force generation for multiple in vitro cardiac models and apply downstream machine learning methods for deep phenotyping and classification. Methods: We first validate BeatProfiler's accuracy, robustness, and speed by benchmarking against existing tools with a fixed dataset. We further confirm its ability to robustly characterize disease and dose-dependent drug response. We then demonstrate that the data acquired by our automatic acquisition pipeline can be further harnessed for machine learning (ML) analysis to phenotype a disease model of restrictive cardiomyopathy and profile cardioactive drug functional response. To accurately classify between these biological signals, we apply feature-based ML and deep learning models (temporal convolutional-bidirectional long short-term memory model or TCN-BiLSTM). Results: Benchmarking against existing tools revealed that BeatProfiler detected and analyzed contraction and calcium signals better than existing tools through improved sensitivity in low signal data, reduction in false positives, and analysis speed increase by 7 to 50-fold. Of signals accurately detected by published methods (PMs), BeatProfiler's extracted features showed high correlations to PMs, confirming that it is reliable and consistent with PMs. The features extracted by BeatProfiler classified restrictive cardiomyopathy cardiomyocytes from isogenic healthy controls with 98% accuracy and identified relax90 as a top distinguishing feature in congruence with previous findings. We also show that our TCN-BiLSTM model was able to classify drug-free control and 4 cardiac drugs with different mechanisms of action at 96% accuracy. We further apply Grad-CAM on our convolution-based models to identify signature regions of perturbations by these drugs in calcium signals. Conclusions: We anticipate that the capabilities of BeatProfiler will help advance in vitro studies in cardiac biology through rapid phenotyping, revealing mechanisms underlying cardiac health and disease, and enabling objective classification of cardiac disease and responses to drugs.

Engineered cardiac tissue model of restrictive cardiomyopathy for drug discovery

Restrictive cardiomyopathy (RCM) is defined as increased myocardial stiffness and impaired diastolic relaxation leading to elevated ventricular filling pressures. Human variants in filamin C (FLNC) are linked to a variety of cardiomyopathies, and in this study, we investigate an in-frame deletion (c.7416_7418delGAA, p.Glu2472_Asn2473delinAsp) in a patient with RCM. Induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) with this variant display impaired relaxation and reduced calcium kinetics in 2D culture when compared with a CRISPR-Cas9-corrected isogenic control line. Similarly, mutant engineered cardiac tissues (ECTs) demonstrate increased passive tension and impaired relaxation velocity compared with isogenic controls. High-throughput small-molecule screening identifies phosphodiesterase 3 (PDE3) inhibition by trequinsin as a potential therapy to improve cardiomyocyte relaxation in this genotype. Together, these data demonstrate an engineered cardiac tissue model of RCM and establish the translational potential of this precision medicine approach to identify therapeutics targeting myocardial relaxation.

Alterations in the kallikrein-kinin system predict death after heart transplant

Heart transplantation remains the definitive treatment for end stage heart failure. Because availability is limited, risk stratification of candidates is crucial for optimizing both organ allocations and transplant outcomes. Here we utilize proteomics prior to transplant to identify new biomarkers that predict post-transplant survival in a multi-institutional cohort. Microvesicles were isolated from serum samples and underwent proteomic analysis using mass spectrometry. Monte Carlo cross-validation (MCCV) was used to predict survival after transplant incorporating select recipient pre-transplant clinical characteristics and serum microvesicle proteomic data. We identified six protein markers with prediction performance above AUROC of 0.6, including Prothrombin (F2), anti-plasmin (SERPINF2), Factor IX, carboxypeptidase 2 (CPB2), HGF activator (HGFAC) and low molecular weight kininogen (LK). No clinical characteristics demonstrated an AUROC > 0.6. Putative biological functions and pathways were assessed using gene set enrichment analysis (GSEA). Differential expression analysis identified enriched pathways prior to transplant that were associated with post-transplant survival including activation of platelets and the coagulation pathway prior to transplant. Specifically, upregulation of coagulation cascade components of the kallikrein-kinin system (KKS) and downregulation of kininogen prior to transplant were associated with survival after transplant. Further prospective studies are warranted to determine if alterations in the KKS contributes to overall post-transplant survival.

Plasma kallikrein predicts primary graft dysfunction after heart transplant

BACKGROUND

Primary graft dysfunction (PGD) is the leading cause of early mortality after heart transplant. Pre-transplant predictors of PGD remain elusive and its etiology remains unclear.

METHODS

Microvesicles were isolated from 88 pre-transplant serum samples and underwent proteomic evaluation using TMT mass spectrometry. Monte Carlo cross validation (MCCV) was used to predict the occurrence of severe PGD after transplant using recipient pre-transplant clinical characteristics and serum microvesicle proteomic data. Putative biological functions and pathways were assessed using gene set enrichment analysis (GSEA) within the MCCV prediction methodology.

RESULTS

Using our MCCV prediction methodology, decreased levels of plasma kallikrein (KLKB1), a critical regulator of the kinin-kallikrein system, was the most predictive factor identified for PGD (AUROC 0.6444 [0.6293, 0.6655]; odds 0.1959 [0.0592, 0.3663]. Furthermore, a predictive panel combining KLKB1 with inotrope therapy achieved peak performance (AUROC 0.7181 [0.7020, 0.7372]) across and within (AUROCs of 0.66-0.78) each cohort. A classifier utilizing KLKB1 and inotrope therapy outperforms existing composite scores by more than 50 percent. The diagnostic utility of the classifier was validated on 65 consecutive transplant patients, resulting in an AUROC of 0.71 and a negative predictive value of 0.92-0.96. Differential expression analysis revealed a enrichment in inflammatory and immune pathways prior to PGD.

CONCLUSIONS

Pre-transplant level of KLKB1 is a robust predictor of post-transplant PGD. The combination with pre-transplant inotrope therapy enhances the prediction of PGD compared to pre-transplant KLKB1 levels alone and the resulting classifier equation validates within a prospective validation cohort. Inflammation and immune pathway enrichment characterize the pre-transplant proteomic signature predictive of PGD.