- 1Department of Aeronautics and Astronautics, Graduate School of Engineering, Kyoto University, Kyoto Daigaku-Katsura, Kyoto, Japan
- 2Cyber-Physical Engineering Information Research Core (Cypher), Okayama University, Okayama, Japan
- 3Department of Mathematics, Graduate School of Science, Kyoto University, Kyoto, Japan
- 4Department of Mechanical Science and Bioengineering, Graduate School of Engineering Science, Osaka University, Osaka, Japan
Stride intervals in human walking fluctuate from one stride to the next, exhibiting statistical persistence. This statistical property is changed by aging, neural disorders, and experimental interventions. It has been hypothesized that the central nervous system is responsible for the statistical persistence. Human walking is a complex phenomenon generated through the dynamic interactions between the central nervous system and the biomechanical system. It has also been hypothesized that the statistical persistence emerges through the dynamic interactions during walking. In particular, a previous study integrated a biomechanical model composed of seven rigid links with a central pattern generator (CPG) model, which incorporated a phase resetting mechanism as sensory feedback as well as feedforward, trajectory tracking, and intermittent feedback controllers, and suggested that phase resetting contributes to the statistical persistence in stride intervals. However, the essential mechanisms remain largely unclear due to the complexity of the neuromechanical model. In this study, we reproduced the statistical persistence in stride intervals using a simplified neuromechanical model composed of a simple compass-type biomechanical model and a simple CPG model that incorporates only phase resetting and a feedforward controller. A lack of phase resetting induced a loss of statistical persistence, as observed for aging, neural disorders, and experimental interventions. These mechanisms were clarified based on the phase response characteristics of our model. These findings provide useful insight into the mechanisms responsible for the statistical persistence of stride intervals in human walking.
1. Introduction
Human walking is not perfectly periodic. The stride interval fluctuates from one stride to the next, exhibiting statistical persistence (Hausdorff et al., 1995; West and Griffin, 1998, 1999; Dingwell and Cusumano, 2010), which indicates that deviations in a time series are statistically more likely to be followed by subsequent deviations in the same direction. Although the stride interval fluctuations change depending on the gait speed and during development from childhood to adulthood, the statistical persistence remains unchanged (Hausdorff et al., 1996, 1999). However, the stride interval fluctuations for elderly subjects (Hausdorff et al., 1997) and patients with Huntington's disease (Hausdorff et al., 1997) or Parkinson's disease (Frenkel-Toledo et al., 2005) become uncorrelated. Experimental interventions for walking, such as the use of a metronome, also make the stride interval fluctuations uncorrelated (Hausdorff et al., 1996). It is largely unclear why statistical persistence appears in stride intervals in human walking and why this statistical property is changed by aging, neural disorders, and experimental interventions.
It has been hypothesized that the central nervous system has an underlying persistence and is responsible for the statistical persistence in stride intervals. This is supported by the finding that statistical persistence remains in patients with significant peripheral nerve degeneration (Gates and Dingwell, 2007). Various neural system models have been developed to reproduce the statistical persistence and investigate the associated mechanisms. Hausdorff et al. (1995) developed a model of the central pattern generators (CPGs) in the spinal cord and introduced “memory” into the CPG model by allowing transitions from frequency to frequency. Ashkenazy et al. (2002) extended this model by introducing a random walk for the signal transmission of neural circuits. West and Scafetta (2003) developed a “Super CPG” model that introduces external interventions via a forced van der Pol oscillator.
Human walking is a complex phenomenon generated through dynamic interactions between the central nervous system and the biomechanical system. It has also been hypothesized that the statistical persistence in stride intervals emerges through complex interactions during walking. Fu et al. (2020) integrated a biomechanical model composed of seven rigid links with a CPG model, which incorporated a phase resetting mechanism as sensory feedback as well as feedforward, trajectory tracking, and intermittent feedback controllers, to reproduce statistical persistence. They showed that a lack of phase resetting induces a loss of statistical persistence. However, it is difficult to fully understand the essential mechanisms responsible for generating and changing this statistical property because of the complexity of the neural and biomechanical models.
In human walking, the stance leg, which is almost straight, rotates around the foot contact point like an inverted pendulum. To investigate the essential mechanisms responsible for generating human walking from a dynamic viewpoint, simple compass-type mechanical models have been used (Kuo, 2001; Donelan et al., 2002; Kuo et al., 2005; Bruijn et al., 2011; Okamoto et al., 2020). Gates et al. (2007) and Ahn and Hogan (2013) reproduced the statistical persistence in stride intervals using simple compass-type models with sensory feedback controllers. However, they did not investigate the contribution of the feedback controllers to changes in the statistical persistence; thus, the essential mechanisms remain unclear.
The aim of this study is to clarify the contribution of phase resetting to the generation and change in the statistical persistence using a simple model. Specifically, we used a simplified neuromechanical model composed of a simple compass-type biomechanical model and a simple CPG model that incorporates phase resetting and a feedforward controller. Our model reproduced the statistical persistence in stride intervals. A lack of phase resetting induced a loss of statistical persistence, as observed in Fu et al. (2020). Furthermore, we clarified the mechanisms responsible for changes in this statistical property caused by phase resetting based on the phase response characteristics. Our findings provide important insights into the mechanisms underlying the generation and change of the statistical persistence in the stride intervals in human walking.
2. Methods
2.1. Mechanical Model
We used a simple compass-type model (Figure 1). This model has two legs (swing and stance legs), the lengths of which are both l, connected by a frictionless hip joint. The masses are located at the hip and on the legs at a distance b from the hip joint; M is the hip mass and m is the leg mass. θ1 is the angle of the stance leg with respect to the vertical, and θ2 is the relative angle between the stance and swing legs. The tip of the stance leg, which corresponds to the ankle, is fixed on the ground. The stance leg rotates freely without friction. This model walks on level ground via joint torques u1 (at the ankle) and u2 (at the hip). g is the acceleration due to gravity. We used the following model parameters based on Winter (2004): M = 50 kg, m = 11 kg, l = 1 m, b = 0.4 m, and g = 9.8 m/s2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/836121/fncir-16-836121-HTML-r1/image_m/fncir-16-836121-g001.jpg)
Figure 1. Neuromechanical model of human walking composed of CPG model with phase resetting and compass-type biomechanical model.
When the tip of the swing leg is in the air, the equations of motion for our model are
The tip of the swing leg touches the ground (touchdown) when the following conditions are satisfied:
We used condition (3) so that touchdown occurs only in front of the model to move forward, and condition (4) to ignore the scuffing of the leg tip on the ground when the swing leg is swung forward. We assumed that touchdown is a fully inelastic collision (no slip, no bounce) and that the stance leg lifts off the ground just after touchdown. Because the roles of the swing and stance legs are reversed just after touchdown, we obtain
where ∗− and ∗+ are the state ∗ just before and after touchdown, respectively. Due to this collision, the angular velocities discontinuously change. We assumed that when the stance leg leaves the ground, it does not interact with the ground and the work of the joint torques can be neglected. These assumptions yield
where
2.2. CPG Model
The CPGs in the spinal cord are largely responsible for rhythmic leg movements, such as during locomotion (Grillner, 1975; Shik and Orlovsky, 1976; Orlovsky et al., 1999). They can produce oscillatory behavior even in the absence of rhythmic input and sensory feedback. However, sensory feedback is crucial for producing adaptive locomotor behavior. To investigate the contribution of CPGs to adaptive locomotion in humans, various oscillator models, such as the van der Pol oscillator (Dutra et al., 2003; West and Scafetta, 2003), Matsuoka oscillator (Matsuoka, 1987; Taga et al., 1991; Taga, 1995a,b; Ogihara and Yamazaki, 2001; Hase et al., 2003; Kim et al., 2011), and phase oscillator (Yamasaki et al., 2003; Aoi et al., 2010, 2019; Dzeladini et al., 2014; Aoi and Funato, 2016; Fu et al., 2020; Tamura et al., 2020; Owaki et al., 2021), have been developed.
In this study, we used a phase oscillator, whose phase is ϕ (0 ≤ ϕ < 2π), to generate the motor commands for our model. The oscillator phase follows the dynamics expressed by
where ω is the basic frequency. We determined the joint torques u1 and u2 as
where A1 and A2 are the amplitudes, σ1 and σ2 are noise terms, and Δ is the phase difference between u1 and u2.
It has been reported that locomotion rhythm and phase are regulated by the production of a phase shift and rhythm resetting (phase resetting) for periodic motor commands in response to sensory feedback (Lafreniere-Roula and McCrea, 2005; Rybak et al., 2006). Cutaneous feedback has been observed to contribute to phase shift and rhythm resetting behavior (Duysens, 1977; Schomburg et al., 1998). Phase resetting has thus been modeled so that the oscillator phase is reset based on foot contact information (Yamasaki et al., 2003; Aoi et al., 2010; Aoi and Funato, 2016; Fu et al., 2020; Tamura et al., 2020). In this study, we used the following relationship at touchdown:
where ϕ0 is a constant. When phase resetting is not applied, ϕ is not regulated at touchdown. However, because the roles of the swing and stance legs are reversed just after touchdown so that (i = 1, 2), we used the following relationship at touchdown:
so that (i = 1, 2) when the noise terms σ1 and σ2 are neglected. We designated ϕ0 as the value to which ϕ+ converged during steady walking (limit cycle) for the model without phase resetting and noise. Therefore, steady walking is identical between the models with and without phase resetting in the absence of noise. This allows us to clearly investigate the difference in the response to torque noise between cases with and without phase resetting.
This CPG model has four parameters, namely ω, A1, A2, and Δ. We used ω = 4.8 rad/s based on Hausdorff et al. (1996). Without noise (σ1 = σ2 = 0), we first investigated the dependence of gait speed during steady walking on A1, A2, and Δ, and then calculated the energy cost for one step cycle for A1, A2, and Δ. We determined the parameter set (A1, A2, Δ) required to minimize ε for each gait speed. When phase resetting was used, we determined ϕ0 for each gait speed using the obtained parameter set.
2.3. Torque Noise
To simulate the stochastic fluctuation of the gait, we used two independent series of white Gaussian noise for torque noise terms σ1 and σ2 in (9) and (10), respectively, as follows:
where ξ is the amplitude of the noise, and U1 and U2 are independent white Gaussian noise with standard deviation 1. This torque noise never induces consecutive touchdowns at extremely short intervals because of discontinuous and large changes in the state variables (5)–(7) at touchdown. We numerically solved the governing equations using the Euler-Maruyama method (Higham, 2001) with a time step of 10−5 s.
To be consistent with previous experiments on humans (Hausdorff et al., 1995, 1996, 1997), a stride was defined as two consecutive steps. Stride intervals were calculated based on the time difference between every other touchdown (strides did not overlap). Each simulation trial required the model to walk 1300 steps (650 strides). The first 150 strides were omitted from the analysis to remove transient behavior due to initial conditions.
2.4. Detrended Fluctuation Analysis
We used detrended fluctuation analysis (DFA) to determine the statistical persistence in the time series of stride intervals for each trial of the computer simulation. This method decreases the effect of noise and removes local trends, making it less affected by non-stationarities. The details of the method can be found elsewhere (e.g., Peng et al., 1993, 1994a,b; Hausdorff et al., 1995; Hardstone et al., 2012; Ihlen, 2012). Briefly, the feature amount F(n) constructed from segments of length n of the time series exhibits a power-law relationship, indicating the presence of scaling as F(n) ~ nα. We investigate the scaling exponent α to determine the statistical persistence for the time series data.
In this study, we first formed the following accumulated sum using the sequence of stride intervals x(i) for i = 1, 2, …, N, where N is the total number of strides (N = 500):
where is the mean stride interval from x(1) to x(N). We then divided the integrated series y(i) into segments of length n (n < N), yj(s) (j = 1, 2, …, N/n, s = 1, 2, …, n), so that each segment is equal in length and non-overlapping. We next detrended each segment yj(s) by subtracting a least squares linear regression line ŷj(s) fit to yj(s), and averaged the squares of the detrended data (i.e., the residuals). We thus obtained the standard deviation F(n) as
We used a set of n distributed equally on a logarithmic scale between 4 and N/4 (Jordan et al., 2006), specifically, n = 4, 5, 6, …, 87, 104, and 125 (sample size is 20).
In general, F(n) increases with increasing n and a graph of log F(n) vs. log n exhibits a power-law relationship, indicating the presence of scaling as F(n) ~ nα. We fit log F(n) vs. log n plots with a linear function using a standard least squares regression approach, and obtained the scaling exponent α from the slope of this line. In particular, α = 0.5 indicates that the stride intervals are completely uncorrelated (i.e., white noise). That is, DFA will still produce α = 0.5 even if the time series is rearranged in any manner (through surrogate data analysis). In contrast, α < 0.5 indicates statistical anti-persistence in stride intervals and 0.5 < α ≤ 1.0 indicates statistical persistence. When α > 1.0, the time series is brown noise (i.e., integrated white noise) (Hausdorff et al., 1995).
3. Results
3.1. Determination of Parameters for Each Gait Speed
Without noise (ξ = 0), our model achieved stable walking with a gait speed v of 0.25 to 0.6 m/s depending on the parameters A1, A2, and Δ. Figure 2A shows the contour of the evaluation criterion ε for A1, A2, and Δ, which generated v = 0.3, 0.4, and 0.5 m/s. Figure 2B shows the parameter sets (A1, A2, Δ), each of which minimized ε for a given gait speed v. The use of phase resetting did not affect these results. We use the parameter set A1 = A1(v), A2 = A2(v), and Δ = Δ(v) in the following sections.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/836121/fncir-16-836121-HTML-r1/image_m/fncir-16-836121-g002.jpg)
Figure 2. Dependence of gait performance on parameters A1, A2, and Δ without noise. (A) Contour of evaluation criterion ε for parameters that generate gait speed v = 0.3, 0.4, and 0.5 m/s. Data point indicates the parameter set that minimizes ε. (B) Parameter sets that minimize ε for each gait speed v.
3.2. Stride Interval Fluctuations
Figure 3 compares the simulation results between the models with and without phase resetting at a walking speed of 0.4 m/s (A1 = 4.9, A2 = 10, Δ = 0.47) using the noise amplitude ξ = 1. Figures 3A,B show the angles θ1 and θ2 and the stride intervals, respectively, during 500 strides. Although ξ is identical between the models, the model without phase resetting has larger stride interval fluctuations than those for the model with phase resetting. Figure 3C shows a plot of log F(n) for log n and the scaling exponent α obtained from the slope of the fitted line. The model with phase resetting exhibits statistical persistence in stride intervals (0.5 < α ≤ 1.0), which is consistent with observations of healthy adults (Hausdorff et al., 1995). Furthermore, the standard deviation of stride interval fluctuations of the model with phase resetting is 0.03, which is also consistent with observations of healthy adults (Hausdorff et al., 1995). In contrast, the model without phase resetting exhibits statistical anti-persistence in stride intervals (α < 0.5). Figure 4 shows the dependence of α on ξ. The models with and without phase resetting, both of which kept walking when ξ ≤ 1, exhibited statistical persistence and anti-persistence, respectively, regardless of ξ.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/836121/fncir-16-836121-HTML-r1/image_m/fncir-16-836121-g003.jpg)
Figure 3. Comparison of gait fluctuations between models with and without phase resetting at gait speed v = 0.4 m/s using noise amplitude ξ = 1 (see Supplementary Movie). (A) Angles θ1 and θ2. Black lines and colored areas indicate the average and standard deviation, respectively. (B) Stride intervals. (C) Plot of log F(n) for log n and scaling exponent α obtained from slope of fitted line.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/836121/fncir-16-836121-HTML-r1/image_m/fncir-16-836121-g004.jpg)
Figure 4. Comparison of scaling exponent α for noise amplitude ξ between models with and without phase resetting at gait speed v = 0.4 m/s. Data points and error bars correspond to the means and standard deviations, respectively, of the results of 10 simulations.
Figure 5 compares the simulation results for the models with and without phase resetting for various values of gait speed v using ξ = 10−2. Figures 5A,B show the stride intervals and log F(n) plot, respectively, for v = 0.3 m/s (A1 = 1.3, A2 = 6.1, Δ = 0.57), 0.4 m/s (A1 = 4.9, A2 = 10, Δ = 0.47), and 0.5 m/s (A1 = 14, A2 = 15, Δ = 0.37). Figure 5C shows the dependence of α on v. The model with phase resetting exhibits statistical persistence regardless of v, which is consistent with observations of healthy adults (Hausdorff et al., 1996). In contrast, the model without phase resetting exhibits statistical anti-persistence regardless of v.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/836121/fncir-16-836121-HTML-r1/image_m/fncir-16-836121-g005.jpg)
Figure 5. Comparison of stride interval fluctuations for various values of gait speed v between models with and without phase resetting using noise amplitude ξ = 10−2. (A) Stride intervals and (B) plot of log F(n) for log n for gait speed v = 0.3, 0.4, and 0.5 m/s. (C) Scaling exponent α vs. gait speed v. Data points and error bars correspond to the means and standard deviations, respectively, of the results of 10 simulations.
4. Discussion
4.1. Mechanisms for Statistical Persistence and Anti-persistence of Stride Intervals
In this study, the model with phase resetting exhibited statistical persistence in stride intervals (0.5 < α ≤ 1.0), whereas the model without phase resetting exhibited statistical anti-persistence (α < 0.5) (Figure 3), as observed in a previous modeling study (Fu et al., 2020). Statistical anti-persistence is characterized by the alternation of large and small values. Fu et al. (2020) performed a linearized stability analysis on a model without phase resetting and noise, and showed that the dominant mode (least stable mode) characterized by Floquet multipliers was a pair of complex conjugates whose amplitude was less than but close to unity and whose argument was greater than π/2. This suggests that the fluctuation ξn of the stride number n can be approximately written as , where r ~ 1 (r < 1) and ξ1 is an initial deviation, corresponding to a slowly damped period-2 oscillation. They explained that this period-2 oscillation induced the alternation of long and short stride intervals and statistical anti-persistence. Although we performed the same stability analysis for our model, the dominant mode of our model without phase resetting and noise was positive real, whose amplitude is less than 1, indicating that the initial deviation monotonically decreases. In addition, our model with phase resetting had almost the same dominant mode as that for our model without phase resetting and it is difficult to conclude that these stability characteristics explain the difference in the statistical properties in stride intervals between the models with and without phase resetting. Furthermore, the amplitude of our dominant mode was 0.65 and the damping was relatively fast.
Next, we directly consider the difference in the response of the stride interval to disturbances. Specifically, we focus on the phase response curve in phase reduction theory (Winfree, 1980; Kuramoto, 1984), which explains how the phase of a limit cycle oscillator shifts by a perturbation at an arbitrary phase (Figure 6). The model with phase resetting shows a shift of the locomotion phase after the recovery due to phase resetting in (11) at foot contact, whereas the model without phase resetting shows no phase shift (Tamura et al., 2020). Furthermore, the phase shift for the model with phase resetting varies depending on the timing of the disturbance. Therefore, the accumulated sum y of stride intervals in (14) tends to move to the cumulative sum of the amount of phase shifts induced by input noise in the model with phase resetting, which results in a relatively smooth signal with large low-frequency components, as shown in Figure 7. In contrast, y tends to converge to 0 in the model without phase resetting, which results in a rough signal with large high-frequency components. Because the scaling exponent α increases with the degree of smoothness (Eke et al., 2000), this difference induces the difference in the scaling exponent α and statistical properties between the models with and without phase resetting.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/836121/fncir-16-836121-HTML-r1/image_m/fncir-16-836121-g006.jpg)
Figure 6. Phase shift caused by disturbance to limit cycle of walking. After recovery, locomotion phase is shifted (x1 + ⋯ + xN > Nτ).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/836121/fncir-16-836121-HTML-r1/image_m/fncir-16-836121-g007.jpg)
Figure 7. Comparison of accumulated sum y of stride intervals between models with and without phase resetting at gait speed v = 0.4 m/s and using noise amplitude ξ = 1 in Figure 3.
4.2. Biological Relevance of Our Findings
The scaling exponent α greatly decreases during walking to a metronome in humans (Hausdorff et al., 1996), where the stride interval is constrained by an external cadence (i.e., metronome). This corresponds to the walking of the model without phase resetting, where the stride interval is constrained by the frequency ω in (8) of the phase oscillator. Therefore, the locomotion phase remains almost unchanged during walking to a metronome, and α decreases as in the model without phase resetting (Figure 3), as discussed in Section 4.1. It has been reported that α also greatly decreases for the stride interval fluctuations of elderly subjects (Hausdorff et al., 1997) and patients with Huntington's disease (Hausdorff et al., 1997) or Parkinson's disease (Frenkel-Toledo et al., 2005). Although the phase response characteristics have been clarified during walking for healthy adults (Yamasaki et al., 2003; Funato et al., 2016; Nessler et al., 2016), those during walking for elderly subjects and patients with neural disorders remain unclear. Investigating them would help clarify the mechanisms responsible for changes in the statistical persistence caused by aging and neural disorders.
Although stride interval fluctuations change depending on gait speed in humans, the statistical persistence remains unchanged (Hausdorff et al., 1996). Our model with phase resetting also exhibited statistical persistence regardless of the gait speed (Figure 5). The constraint on gait rhythm seems more crucial for the statistical persistence than the constraint on gait speed, as observed for walking to a metronome (Hausdorff et al., 1996).
The standard deviation of stride interval fluctuations is about 0.04 s in human walking, which is 3% of the mean stride interval (Hausdorff et al., 1995). It was difficult for previous studies (Gates et al., 2007; Fu et al., 2020) using biomechanical models to reproduce a magnitude of stride interval fluctuations similar to that for humans. Although Gates et al. (2007) reproduced statistical persistence in stride intervals (0.5 < α ≤ 1.0) using a simple biomechanical model as in this study, their model was not robust and the noise amplitude was limited. Therefore, their stride interval fluctuations were much smaller than those in humans. Furthermore, the scaling exponent α was sensitive to the noise amplitude, and the fluctuations exhibited brown noise at high noise levels (α > 1.0). In contrast, phase resetting made our model robust, which allowed a magnitude of stride interval fluctuations similar to that for healthy adults (Figure 3). Furthermore, α was 0.5 to 1.0, which is consistent with observations of healthy adults, and was not sensitive to the noise amplitude (Figure 4), but sensitive to the controller (i.e., whether phase resetting was used).
Previous studies (Yamasaki et al., 2003; Aoi et al., 2010; Fujiki et al., 2018; Tamura et al., 2020) have shown that phase resetting contributes to adaptive walking. In this study, we found that it also contributes to the statistical persistence of gait. In addition to the fact that statistical persistence is impaired by aging (Hausdorff et al., 1997), central nervous system diseases, such as Parkinson's disease (Frenkel-Toledo et al., 2005) and Huntington's disease (Hausdorff et al., 1997), and experimental intervention for walking (Hausdorff et al., 1996), it has been suggested that statistical persistence is linked to important characteristics of gait. Bohnsack-McLagan et al. (2016) suggested that fluctuation persistence leads to redundancies in gait and helps predict and prevent fall risk. Ahn and Hogan (2013) and Fu et al. (2020) showed that fluctuation persistence appears in gait with low gait stability. Gates et al. (2007) showed that a decrease in the ability to perform finely controlled movements leads to an increase in motor output noise and impairs the persistence of fluctuations.
Many studies have reported long-range correlations in stride intervals in human walking based on the results of DFA (Hausdorff et al., 1995, 1996, 1997; Ashkenazy et al., 2002), which indicates that stride-to-stride correlations decay in a scale-free (fractal-like) power-law fashion and suggests that each stride depends explicitly on many previous strides. However, DFA is highly sensitive to yielding false positive results (Maraun et al., 2004; Höll and Kantz, 2015), and it is difficult to conclude the presence of long-range correlations from DFA alone. Instead, DFA provides a valid indicator of statistical persistence and anti-persistence in a time series (Maraun et al., 2004). In this study, we used statistical persistence instead of long-range correlations to interpret the results of DFA, as discussed in Dingwell and Cusumano (2010).
4.3. Limitations of Our Model and Future Work
Based on the hypothesis that the statistical persistence in stride intervals emerge through dynamic interactions between the neural and biomechanical systems, we integrated a simple neural model and a simple biomechanical model to reproduce statistical persistence in stride intervals and change in this statistical property. However, our model is very simple and has limitations with regard to replicating many aspects of human walking. In particular, because the feedforward torques (9) and (10) were simply composed of a sinusoidal wave, the gait speeds of our model were slower than those of healthy adults (Figure 2). In addition, although statistical persistence could be associated with low gait stability (low convergence speed to the limit cycle) (Ahn and Hogan, 2013; Fu et al., 2020), our model had higher stability than that of complicated models due to its simplicity. The high stability of our model with phase resetting might have caused the scaling exponent α to be ~ 0.6, which is smaller than that (~ 1) in healthy adults (Hausdorff et al., 1995). Furthermore, stochastic noise is ubiquitous in the central nervous system and peripheral sensory-motor systems (Jones et al., 2002; van Beers et al., 2004; Churchland et al., 2006). However, our model used only torque noise, which may result in the difference between the statistical anti-persistence in the model without phase resetting and the white noise in walking to a metronome in humans (Hausdorff et al., 1996; Bohnsack-McLagan et al., 2016).
Based on the findings in this study, it is important to verify the essential mechanisms responsible for changes in the statistical persistence by using biologically detailed neuromusculoskeletal models. In a previous study (Tamura et al., 2020), we integrated a musculoskeletal model composed of seven rigid links and 18 muscles with a CPG model with a muscle synergy-based controller to investigate the contribution of phase resetting to the phase response characteristics during walking. In another previous study (Fujiki et al., 2019), we used a half-center type CPG model composed of a rhythm generator network, which was modeled using neuron populations of flexor and extensor centers based on Danner et al. (2016, 2017) and Rybak et al. (2006), to clarify the mechanisms responsible for the CPG responses to afferent stimulation using dynamic systems theory based on nullclines. We plan to incorporate these biologically detailed models to further investigate the mechanisms responsible for changes in the statistical persistence.
5. Conclusion
In this study, we clarified the contribution of phase resetting to the generation and change of statistical persistence using a simple neuromechanical model. Specifically, our model reproduced the statistical persistence in stride intervals. A lack of phase resetting induced a loss of statistical persistence. Furthermore, we clarified the mechanisms responsible for changes in statistical persistence caused by phase resetting based on the phase response characteristics. Our findings provide important insight into the mechanisms underlying the generation and change of the statistical persistence in the stride intervals in human walking.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
SA developed the study design. KO performed simulation experiments and analyzed the data in consultation with SA, IO, HK, KS, and KT. KO and SA wrote the manuscript. All authors reviewed and approved it.
Funding
This study was supported in part by JSPS KAKENHI Grant Numbers JP21J23164 and JP20H00229; and JST FOREST Program Grant Number JPMJFR2021.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer TN declared a shared affiliation with the author SA to the handling editor at time of review.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncir.2022.836121/full#supplementary-material
References
Ahn, J., and Hogan, N. (2013). Long-range correlations in stride intervals may emerge from non-chaotic walking dynamics. PLoS ONE 8, e73239. doi: 10.1371/journal.pone.0073239
Aoi, S., and Funato, T. (2016). Neuromusculoskeletal models based on the muscle synergy hypothesis for the investigation of adaptive motor control in locomotion via sensory-motor coordination. Neurosci. Res. 104, 88–95. doi: 10.1016/j.neures.2015.11.005
Aoi, S., Ogihara, N., Funato, T., Sugimoto, Y., and Tsuchiya, K. (2010). Evaluating functional roles of phase resetting in generation of adaptive human bipedal walking with a physiologically based model of the spinal pattern generator. Biol. Cybern. 102, 373–387. doi: 10.1007/s00422-010-0373-y
Aoi, S., Ohashi, T., Bamba, R., Fujiki, S., Tamura, D., Funato, T., et al. (2019). Neuromusculoskeletal model that walks and runs across a speed range with a few motor control parameter changes based on the muscle synergy hypothesis. Sci. Rep. 9, 369. doi: 10.1038/s41598-018-37460-3
Ashkenazy, Y., Hausdorff, J. M., Ivanov, P. C., and Stanley, H. E. (2002). A stochastic model of human gait dynamics. Phys. A 316, 662–670. doi: 10.1016/S0378-4371(02)01453-X
Bohnsack-McLagan, N. K., Cusumano, J. P., and Dingwell, J. B. (2016). Adaptability of stride-to-stride control of stepping movements in human walking. J. Biomech. 49, 229–237. doi: 10.1016/j.jbiomech.2015.12.010
Bruijn, S. M., Bregman, D. J., Meijer, O. G., Beek, P. J., and van Dieën, J. H. (2011). The validity of stability measures: a modelling approach. J. Biomech. 44, 2401–2408. doi: 10.1016/j.jbiomech.2011.06.031
Churchland, M. M., Afshar, A., and Shenoy, K. V. (2006). A central source of movement variability. Neuron 52, 1085–1096. doi: 10.1016/j.neuron.2006.10.034
Danner, S. M., Shevtsova, N. A., Frigon, A., and Rybak, I. A. (2017). Computational modeling of spinal circuits controlling limb coordination and gaits in quadrupeds. eLife 6, e31050. doi: 10.7554/eLife.31050
Danner, S. M., Wilshin, S. D., Shevtsova, N. A., and Rybak, I. A. (2016). Central control of interlimb coordination and speed-dependent gait expression in quadrupeds. J. Physiol. 594, 6947–6967. doi: 10.1113/JP272787
Dingwell, J. B., and Cusumano, J. P. (2010). Re-interpreting detrended fluctuation analyses of stride-to-stride variability in human walking. Gait Post. 32, 348–353. doi: 10.1016/j.gaitpost.2010.06.004
Donelan, J. M., Kram, R., and Kuo, A. D. (2002). Mechanical work for step-to-step transitions is a major determinant of the metabolic cost of human walking. J. Exp. Biol. 205, 3717–3727. doi: 10.1242/jeb.205.23.3717
Dutra, M. S., de Pina Filho, A. C., and Romano, V. F. (2003). Modeling of a bipedal locomotor using coupled nonlinear oscillators of Van der Pol. Biol. Cybern. 88, 286–292. doi: 10.1007/s00422-002-0380-8
Duysens, J. (1977). Fluctuations in sensitivity to rhythm resetting effects during the cat's step cycle. Brain Res. 133, 190–195. doi: 10.1016/0006-8993(77)90063-4
Dzeladini, F., van den Kieboom, J., and Ijspeert, A. (2014). The contribution of a central pattern generator in a reflex-based neuromuscular model. Front. Hum. Neurosci. 8, 371. doi: 10.3389/fnhum.2014.00371
Eke, A., Hermán, P., Bassingthwaighte, J. B., Raymond, G. M., Percival, D. B., Cannon, M., et al. (2000). Physiological time series: distinguishing fractal noises from motions. Pflügers Arch. Eur. J. Physiol. 439, 403–415. doi: 10.1007/s004249900135
Frenkel-Toledo, S., Giladi, N., Peretz, C., Herman, T., Gruendlinger, L., and Hausdorff, J. M. (2005). Treadmill walking as an external pacemaker to improve gait rhythm and stability in Parkinson's disease. Mov. Disord. 20, 1109–1114. doi: 10.1002/mds.20507
Fu, C., Suzuki, Y., Morasso, P., and Nomura, T. (2020). Phase resetting and intermittent control at the edge of stability in a simple biped model generates 1/f-like gait cycle variability. Biol. Cybern. 114, 95–111. doi: 10.1007/s00422-020-00816-y
Fujiki, S., Aoi, S., Funato, T., Sato, Y., Tsuchiya, and Yanagihara, D. (2018). Adaptive hindlimb split-belt treadmill walking in rats by controlling basic muscle activation patterns via phase resetting. Sci. Rep. 8, 17341. doi: 10.1038/s41598-018-35714-8
Fujiki, S., Aoi, S., Tsuchiya, K., Danner, S. M., Rybak, I. A., and Yanagihara, D. (2019). Phase-dependent response to afferent stimulation during fictive locomotion: a computational modeling study. Front. Neurosci. 13, 1288. doi: 10.3389/fnins.2019.01288
Funato, T., Yamamoto, Y., Aoi, S., Imai, T., Aoyagi, T., Tomita, N., et al. (2016). Evaluation of the phase-dependent rhythm control of human walking using phase response curves. PLoS Comput. Biol. 12, e1004950. doi: 10.1371/journal.pcbi.1004950
Gates, D. H., and Dingwell, J. B. (2007). Peripheral neuropathy does not alter the fractal dynamics of stride intervals of gait. J. Appl. Physiol. 102, 965–971. doi: 10.1152/japplphysiol.00413.2006
Gates, D. H., Su, J. L., and Dingwell, J. B. (2007). Possible biomechanical origins of the long-range correlations in stride intervals of walking. Phys. A 380, 259–270. doi: 10.1016/j.physa.2007.02.061
Grillner, S. (1975). Locomotion in vertebrates: central mechanisms and reflex interaction. Physiol. Rev. 55, 247–304. doi: 10.1152/physrev.1975.55.2.247
Hardstone, R., Poil, S.-S., Schiavone, G., Jansen, R., Nikulin, V. V., Mansvelder, H. D., et al. (2012). Detrended fluctuation analysis: a scale-free view on neuronal oscillations. Front. Physiol. 3, 450. doi: 10.3389/fphys.2012.00450
Hase, K., Miyashita, K., Ok, S., and Arakawa, Y. (2003). Human gait simulation with a neuromusculoskeletal model and evolutionary computation. J. Vis. Comput. Animat. 14, 73–92. doi: 10.1002/vis.306
Hausdorff, J. M., Mitchell, S. L., Firtion, R., Peng, C. K., Cudkowicz, M. E., Wei, J. Y., et al. (1997). Altered fractal dynamics of gait: reduced stride-interval correlations with aging and Huntington's disease. J. Appl. Physiol. 82, 262–269. doi: 10.1152/jappl.1997.82.1.262
Hausdorff, J. M., Peng, C. K., Ladin, Z., Wei, J. Y., and Goldberger, A. L. (1995). Is walking a random walk? Evidence for long-range correlations in stride interval of human gait. J. Appl. Physiol. 78, 349–358. doi: 10.1152/jappl.1995.78.1.349
Hausdorff, J. M., Purdon, P. L., Peng, C. K., Ladin, Z., Wei, J. Y., and Goldberger, A. L. (1996). Fractal dynamics of human gait: stability of long-range correlations in stride interval fluctuations. J. Appl. Physiol. 80, 1448–1457. doi: 10.1152/jappl.1996.80.5.1448
Hausdorff, J. M., Zemany, L., Peng, C.-K., and Goldberger, A. L. (1999). Maturation of gait dynamics: stride-to-stride variability and its temporal organization in children. J. Appl. Physiol. 86, 1040–1047. doi: 10.1152/jappl.1999.86.3.1040
Higham, D. J. (2001). An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Rev. 43, 525–546. doi: 10.1137/S0036144500378302
Höll, M., and Kantz, H. (2015). The fluctuation function of the detrended fluctuation analysis - investigation on the AR(1) process. Eur. Phys. J. B 88, 126. doi: 10.1140/epjb/e2015-60143-1
Ihlen, E. A. F. (2012). Introduction to multifractal detrended fluctuation analysis in matlab. Front. Physiol. 3, 141. doi: 10.3389/fphys.2012.00141
Jones, K. E., Hamilton, A. F. C., and Wolpert, D. M. (2002). Sources of signal-dependent noise during isometric force production. J. Neurophysiol. 88, 1533–1544. doi: 10.1152/jn.2002.88.3.1533
Jordan, K., Challis, J. H., and Newell, K. M. (2006). Long range correlations in the stride interval of running. Gait Post. 24, 120–125. doi: 10.1016/j.gaitpost.2005.08.003
Kim, Y., Tagawa, Y., Obinata, G., and Hase, K. (2011). Robust control of CPG-based 3D neuromusculoskeletal walking model. Biol. Cybern. 105, 269–282. doi: 10.1007/s00422-011-0464-4
Kuo, A. D. (2001). A simple model of bipedal walking predicts the preferred speed-step length relationship. J. Biomech. Eng. 123, 264–269. doi: 10.1115/1.1372322
Kuo, A. D., Donelan, J. M., and Ruina, A. (2005). Energetic consequences of walking like an inverted pendulum: step-to-step transitions. Exerc. Sport Sci. Rev. 33, 88–97. doi: 10.1097/00003677-200504000-00006
Kuramoto, Y. (1984). Chemical Oscillations, Waves, and Turbulence. Berlin: Springer-Verlag. doi: 10.1007/978-3-642-69689-3
Lafreniere-Roula, M., and McCrea, D. A. (2005). Deletions of rhythmic motoneuron activity during fictive locomotion and scratch provide clues to the organization of the mammalian central pattern generator. J. Neurophysiol. 94, 1120–1132. doi: 10.1152/jn.00216.2005
Maraun, D., Rust, H. W., and Timmer, J. (2004). Tempting long-memory - on the interpretation of DFA results. Neuromechanical model. Nonlinear Process. Geophys. 11, 495–503. doi: 10.5194/npg-11-495-2004
Matsuoka, K. (1987). Mechanisms of frequency and pattern control in the neural rhythm generators. Biol. Cybern. 56, 345–353. doi: 10.1007/BF00319514
Nessler, J. A., Spargo, T., Craig-Jones, A., and Milton, J. G. (2016). Phase resetting behavior in human gait is influenced by treadmill walking speed. Gait Post. 43, 187–191. doi: 10.1016/j.gaitpost.2015.09.021
Ogihara, N., and Yamazaki, N. (2001). Generation of human bipedal locomotion by a bio-mimetic neuro-musculo-skeletal model. Biol. Cybern. 84, 1–11. doi: 10.1007/PL00007977
Okamoto, K., Aoi, S., Obayashi, I., Kokubu, H., Senda, K., and Tsuchiya, K. (2020). Fractal mechanism of basin of attraction in passive dynamic walking. Bioinspir. Biomim. 15, 055002. doi: 10.1088/1748-3190/ab9283
Orlovsky, G., Deliagina, T. G., and Grillner, S. (1999). Neuronal Control of Locomotion from Mollusc to Man. New York, NY: Oxford University Press. doi: 10.1093/acprof:oso/9780198524052.001.0001
Owaki, D., Horikiri, S., Nishii, J., and Ishiguro, A. (2021). Tegotae-based control produces adaptive inter- and intra-limb coordination in bipedal walking. Front. Neurorobot. 15, 629595. doi: 10.3389/fnbot.2021.629595
Peng, C.-K., Buldyrev, S. V., Goldberger, A. L., Havlin, S., Simons, M., and Stanley, H. E. (1993). Finite-size effects on long-range correlations: implications for analyzing DNA sequences. Phys. Rev. E 47, 3730–3733. doi: 10.1103/PhysRevE.47.3730
Peng, C.-K., Buldyrev, S. V., Hausdorff, J. M., Havlin, S., Mietus, J. E., Simons, M., et al. (1994a). Non-equilibrium dynamics as an indispensable characteristic of a healthy biological system. Integr. Physiol. Behav. Sci. 29, 283–293. doi: 10.1007/BF02691332
Peng, C.-K., Buldyrev, S. V., Havlin, S., Simons, M., Stanley, H. E., and Goldberger, A. L. (1994b). Mosaic organization of DNA nucleotides. Phys. Rev. E 49, 1685–1689. doi: 10.1103/PhysRevE.49.1685
Rybak, I. A., Shevtsova, N. A., Lafreniere-Roula, M., and McCrea, D. A. (2006). Modelling spinal circuitry involved in locomotor pattern generation: insights from deletions during fictive locomotion. J. Physiol. 577, 617–639. doi: 10.1113/jphysiol.2006.118703
Schomburg, E. D., Petersen, N., Barajon, I., and Hultborn, H. (1998). Flexor reflex afferents reset the step cycle during fictive locomotion in the cat. Exp. Brain Res. 122, 339–350. doi: 10.1007/s002210050522
Shik, M. L., and Orlovsky, G. N. (1976). Neurophysiology of locomotor automatism. Physiol. Rev. 56, 465–501. doi: 10.1152/physrev.1976.56.3.465
Taga, G. (1995a). A model of the neuro-musculo-skeletal system for human locomotion I. Emergence of basic gait. Biol. Cybern. 73, 97–111. doi: 10.1007/BF00204048
Taga, G. (1995b). A model of the neuro-musculo-skeletal system for human locomotion II. Real-time adaptability under various constraints. Biol. Cybern. 73, 113–121. doi: 10.1007/BF00204049
Taga, G., Yamaguchi, Y., and Shimizu, H. (1991). Self-organized control of bipedal locomotion by neural oscillators in unpredictable environment. Biol. Cybern. 65, 147–159. doi: 10.1007/BF00198086
Tamura, D., Aoi, S., Funato, T., Fujiki, S., Senda, K., and Tsuchiya, K. (2020). Contribution of phase resetting to adaptive rhythm control in human walking based on the phase response curves of a neuromusculoskeletal model. Front. Neurosci. 14, 17. doi: 10.3389/fnins.2020.00017
van Beers, R. J., Haggard, P., and Wolpert, D. M. (2004). The role of execution noise in movement variability. J. Neurophysiol. 91, 1050–1063. doi: 10.1152/jn.00652.2003
West, B. J., and Griffin, L. (1998). Allometric control of human gait. Fractals 6, 101–108. doi: 10.1142/S0218348X98000122
West, B. J., and Griffin, L. (1999). Allometric control, inverse power laws and human gait. Chaos Soliton. Fract. 10, 1519–1527. doi: 10.1016/S0960-0779(98)00149-0
West, B. J., and Scafetta, N. (2003). Nonlinear dynamical model of human gait. Phys. Rev. E 67, 051917. doi: 10.1103/PhysRevE.67.051917
Winfree, A. T. (1980). The Geometry of Biological Time. New York, NY: Springer-Verlag. doi: 10.1007/978-3-662-22492-2
Winter, D. A. (2004). Biomechanics and Motor Control of Human Movement, 3rd Edn. New York, NY: Wiley.
Keywords: human walking, stride interval fluctuation, neuromechanical model, central pattern generator, phase resetting, statistical persistence
Citation: Okamoto K, Obayashi I, Kokubu H, Senda K, Tsuchiya K and Aoi S (2022) Contribution of Phase Resetting to Statistical Persistence in Stride Intervals: A Modeling Study. Front. Neural Circuits 16:836121. doi: 10.3389/fncir.2022.836121
Received: 15 December 2021; Accepted: 19 May 2022;
Published: 22 June 2022.
Edited by:
Monica A. Daley, University of California, Irvine, United StatesReviewed by:
Auke Ijspeert, Swiss Federal Institute of Technology Lausanne, SwitzerlandTaishin Nomura, Osaka University, Japan
Jonathan B. Dingwell, The Pennsylvania State University (PSU), United States
Copyright © 2022 Okamoto, Obayashi, Kokubu, Senda, Tsuchiya and Aoi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shinya Aoi, aoi.shinya.es@osaka-u.ac.jp