Nonlinear control in the nematode C. elegans

Megan Morrison,Charles Fieseler,J. Nathan Kutz

DOI: https://doi.org/10.48550/arXiv.2001.08332

2020-02-04

Abstract:Recent whole-brain calcium imaging recordings of the nematode C. elegans have demonstrated that neural activity is dominated by dynamics on a low-dimensional manifold that can be clustered according to behavioral states. Despite progress in modeling the dynamics with linear or locally linear models, it remains unclear how a single network of neurons can produce the observed features. In particular, there are multiple clusters, or fixed points, observed in the data which cannot be characterized by a single linear model. We propose a nonlinear control model which is global and parameterized by only four free parameters that match the features displayed by the low-dimensional C. elegans neural activity. In addition to reproducing the average probability distribution of the data, long and short time-scale changes in transition statistics can be characterized via changes in a single parameter. Some of these macro-scale transitions have experimental correlates to single neuro-modulators that seem to act as biological controls, allowing this model to generate testable hypotheses about the effect of these neuro-modulators on the global dynamics. The theory provides an elegant characterization of the neuron population dynamics in C. elegans. Moreover, the mathematical structure of the nonlinear control framework provides a paradigm that can be generalized to more complex systems with an arbitrary number of behavioral states.

Neurons and Cognition,Dynamical Systems,Biological Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to construct a global nonlinear control model that can explain the relationship between the neural activity and behavioral states of Caenorhabditis elegans (C. elegans). Specifically, the existing linear or locally linear models cannot fully describe the multiple fixed points (i.e., different behavioral states) observed in the neural activity of C. elegans. Therefore, a more flexible and global nonlinear model is required to capture these features.

### Core of the problem

1. **Multi - fixed - point problem**: Experimental data show that the neural activity of C. elegans exhibits multiple fixed points on a low - dimensional manifold, and these fixed points correspond to different behavioral states (such as forward movement, backward movement, etc.). Traditional linear models can only describe a single fixed point and thus cannot explain this multi - fixed - point phenomenon.

2. **Dynamics of behavioral transitions**: It is necessary not only to describe the existence of different behavioral states but also to explain the transition mechanism from one behavioral state to another. In particular, the time scales of these transitions can vary greatly, from short - time rapid transitions to long - time changes in behavioral patterns.

3. **Influence of neuromodulators**: Certain neuromodulators have a significant impact on the behavior of C. elegans, for example, changing the reversal frequency. The model needs to be able to explain how these neuromodulators affect the global dynamics.

### Solution

The paper proposes a model based on nonlinear control theory, which has the following characteristics:

- **Globality and flexibility**: By introducing nonlinear terms, the model can describe multiple fixed points globally and can simulate the transitions between different behavioral states by adjusting a small number of parameters.

- **Simple parameterization**: The model only needs four free parameters to match the neural activity features observed in experiments, which makes the model both concise and highly explanatory.

- **Biological relevance**: Some parameters in the model may correspond to specific neuromodulators, thus providing a theoretical basis for understanding the impact of these substances on behavior.

### Specific form of the model

The form of this nonlinear control model is as follows:

\[
x' = y \quad (1a)
\]
\[
y' = f(x, \beta)+\gamma y + u(t) \quad (1b)
\]

where the nonlinear dynamic term \( f(x, \beta) \) is defined as:

\[
f(x, \beta)=-(x + 1)(x-\beta)(x - 1) \quad (2)
\]

This cubic polynomial constructs two stable fixed points \( x=\pm1 \) and an unstable fixed point determined by the parameter \( \beta \). In addition, the model also includes a damping parameter \( \gamma \), a control input \( u(t) \), and a random term \( \sigma dW_t \) to describe the noise fluctuations observed in the experiment.

### Conclusion

Through the above - mentioned nonlinear control model, the paper successfully explains the existence of multiple fixed points in the neural activity of C. elegans and their transition mechanisms, and can reproduce behavioral changes on different time scales by adjusting a single parameter. This model not only provides a new perspective for understanding the behavior of C. elegans but also lays the foundation for further research on neural dynamics in other complex biological systems.