Optimal phase control of biological oscillators using augmented phase reduction

• Published in: Biological cybernetics Vol. 113; no. 1-2; pp. 161 - 178
• Main Authors: Monga, Bharat, Moehlis, Jeff
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01-04-2019; Springer Nature B.V
• Subjects: Algorithms; Basal ganglia; Bioinformatics; Biomedical and Life Sciences; Biomedicine; Central nervous system diseases; Circadian rhythms; Complex Systems; Computer Appl. in Life Sciences; Control algorithms; Control theory; Gene expression; Gene regulation; Heart; Hopf bifurcation; Jet lag; Movement disorders; Neurobiology; Neurosciences; Optimal control; Original Article; Oscillators; Phase control; Phase transitions; Reduction; Suprachiasmatic nucleus; Thalamus; Tremor; Circadian rhythms; Phase reduction; Alternans; Spike timing control; Optimal control
• ISSN: 0340-1200; 1432-0770
• DOI: 10.1007/s00422-018-0764-z

We develop a novel optimal control algorithm to change the phase of an oscillator using a minimum energy input, which also minimizes the oscillator’s transversal distance to the uncontrolled periodic orbit. Our algorithm uses a two-dimensional reduction technique based on both isochrons and isostables. We develop a novel method to eliminate cardiac alternans by connecting our control algorithm with the underlying physiological problem. We also describe how the devised algorithm can be used for spike timing control which can potentially help with motor symptoms of essential and parkinsonian tremor, and aid in treating jet lag. To demonstrate the advantages of this algorithm, we compare it with a previously proposed optimal control algorithm based on standard phase reduction for the Hopf bifurcation normal form, and models for cardiac pacemaker cells, thalamic neurons, and circadian gene regulation cycle in the suprachiasmatic nucleus. We show that our control algorithm is effective even when a large phase change is required or when the nontrivial Floquet multiplier is close to unity; in such cases, the previously proposed control algorithm fails.