Elucidation of molecular kinetic schemes from macroscopic traces using system identification

Overall cellular responses to biologically-relevant stimuli are mediated by networks of simpler lower-level processes. Although information about some of these processes can now be obtained by visualizing and recording events at the molecular level, this is still possible only in especially favorable cases. Therefore the development of methods to extract the dynamics and relationships between the different lower-level (microscopic) processes from the overall (macroscopic) response remains a crucial challenge in the understanding of many aspects of physiology. Here we have devised a hybrid computational-analytical method to accomplish this task, the SYStems-based MOLecular kinetic scheme Extractor (SYSMOLE). SYSMOLE utilizes system-identification input-output analysis to obtain a transfer function between the stimulus and the overall cellular response in the Laplace-transformed domain. It then derives a Markov-chain state molecular kinetic scheme uniquely associated with the transfer function by means of a classification procedure and an analytical step that imposes general biological constraints. We first tested SYSMOLE with synthetic data and evaluated its performance in terms of its rate of convergence to the correct molecular kinetic scheme and its robustness to noise. We then examined its performance on real experimental traces by analyzing macroscopic calcium-current traces elicited by membrane depolarization. SYSMOLE derived the correct, previously known molecular kinetic scheme describing the activation and inactivation of the underlying calcium channels and correctly identified the accepted mechanism of action of nifedipine, a calcium-channel blocker clinically used in patients with cardiovascular disease. Finally, we applied SYSMOLE to study the pharmacology of a new class of glutamate antipsychotic drugs and their crosstalk mechanism through a heteromeric complex of G protein-coupled receptors. Our results indicate that our methodology can be successfully applied to accurately derive molecular kinetic schemes from experimental macroscopic traces, and we anticipate that it may be useful in the study of a wide variety of biological systems.Author summary: Unraveling the lower-level (microscopic) processes underlying the overall (macroscopic) cell response to a given stimulus is a challenging problem in cell physiology. This has been a classic problem in biophysics, where the ability to record the activity of single ion channels that generate a macroscopic ion current has allowed a measure of direct access to the underlying microscopic processes. These classic studies have demonstrated that very different groupings of the microscopic processes can yield extremely similar macroscopic responses. Biologists in fields other than biophysics are frequently confronted with the same macroscopic-to-microscopic problem, usually, however, without any direct access to the microscopic processes. Thus, the development of computational methods to deduce from the available macroscopic measurements the nature of the underlying microscopic processes can be expected to substantially advance the study of many areas of cell physiology. Toward that aim, here we have derived and tested a hybrid computational-analytical method to extract information about the microscopic processes that is hidden in macroscopic experimental traces. Our method is independent of the particular system under study, and thus can be applied to new as well as previously-recorded macroscopic traces obtained in a wide variety of biological systems.