Lung cancer is still a leading cause of death worldwide. In recent years, knowledge has been obtained of the mechanisms modulating ion channel kinetics and thus of cell bioelectric properties, which is promising for oncological biomarkers and targets. The complex interplay of channel expression and its consequences on malignant processes, however, is still insufficiently understood. We here introduce the first approach of an in-silico whole-cell ion current model of a cancer cell, in particular of the A549 human lung adenocarcinoma, including the main functionally expressed ion channels in the plasma membrane as so far known. This hidden Markov-based model represents the electrophysiology behind proliferation of the A549 cell, describing its rhythmic oscillation of the membrane potential able to trigger the transition between cell cycle phases, and it predicts membrane potential changes over the cell cycle provoked by targeted ion channel modulation. This first A549 in-silico cell model opens up a deeper insight and understanding of possible ion channel interactions in tumor development and progression, and is a valuable tool for simulating altered ion channel function in lung cancer electrophysiology.

Metabolic biomarkers may play an important role in the diagnosis, prognostication and assessment of response to pharmacological therapy in complex diseases. The process of discovering new metabolic biomarkers is a non-trivial task which involves a number of bioanalytical processing steps coupled with a computational approach for the search, prioritization and verification of new biomarker candidates. Kinetic analysis provides an additional dimension of complexity in time-series data, allowing for a more precise interpretation of biomarker dynamics in terms of molecular interaction and pathway modulation. A novel network-based computational strategy for the discovery of putative dynamic biomarker candidates is presented, enabling the identification and verification of unexpected metabolic signatures in complex diseases such as myocardial infarction. The novelty of the proposed method lies in combining metabolic time-series data into a superimposed graph representation, highlighting the strength of the underlying kinetic interaction of preselected analytes. Using this approach, we were able to confirm known metabolic signatures and also identify new candidates such as carnosine and glycocholic acid, and pathways that have been previously associated with cardiovascular or related diseases. This computational strategy may serve as a complementary tool for the discovery of dynamic metabolic or proteomic biomarkers in the field of clinical medicine.

Excessive expression of subunit 1 of GIRK1 in ER+ breast tumors is associated with reduced survival times and increased lymph node metastasis in patients. To investigate possible tumor-initiating properties, benign MCF10A and malign MCF7 mammary epithelial cells were engineered to overexpress GIRK1 neoplasia associated vital parameters and resting potentials were measured and compared to controls. The presence of GIRK1 resulted in resting potentials negative to the controls. Upon GIRK1 overexpression, several cellular pathways were regulated towards pro-tumorigenic action as revealed by comparison of transcriptomes of MCF10AGIRK1 with the control (MCF10AeGFP). According to transcriptome analysis, cellular migration was promoted while wound healing and extracellular matrix interactions were impaired. Vital parameters in MCF7 cells were affected akin the benign MCF10A lines, but to a lesser extent. Thus, GIRK1 regulated cellular pathways in mammary epithelial cells are likely to contribute to the development and progression of breast cancer.

The mathematical modeling of ion channel kinetics is an important tool for studying the electrophysiological mechanisms of the nerves, heart, or cancer, from a single cell to an organ. Common approaches use either a Hodgkin–Huxley (HH) or a hidden Markov model (HMM) description, depending on the level of detail of the functionality and structural changes of the underlying channel gating, and taking into account the computational effort for model simulations. Here, we introduce for the first time a novel system theory-based approach for ion channel modeling based on the concept of transfer function characterization, without a priori knowledge of the biological system, using patch clamp measurements. Using the shaker-related voltage-gated potassium channel Kv1.1 (KCNA1) as an example, we compare the established approaches, HH and HMM, with the system theory-based concept in terms of model accuracy, computational effort, the degree of electrophysiological interpretability, and methodological limitations. This highly data-driven modeling concept offers a new opportunity for the phenomenological kinetic modeling of ion channels, exhibiting exceptional accuracy and computational efficiency compared to the conventional methods. The method has a high potential to further improve the quality and computational performance of complex cell and organ model simulations, and could provide a valuable new tool in the field of next-generation in silico electrophysiology.

No abstract