Skip to main content

ORIGINAL RESEARCH article

Front. Comput. Neurosci., 05 October 2023

Analysis of abnormal posture in patients with Parkinson's disease using a computational model considering muscle tones

  • 1Department of Precision Engineering, School of Engineering, The University of Tokyo, Tokyo, Japan
  • 2Department of Integrated Neuroanatomy and Neuroimaging, Kyoto University Graduate School of Medicine and Faculty of Medicine, Kyoto, Japan
  • 3Department of Advanced Neuroimaging, Integrative Brain Imaging Center, National Center of Neurology and Psychiatry, Tokyo, Japan
  • 4Research into Artifacts, Center for Engineering, School of Engineering, The University of Tokyo, Tokyo, Japan
  • 5Division on Neuroscience, Department of Physiology, Asahikawa Medical University, Asahikawa, Japan
  • 6Department of Neurology, National Center Hospital, National Center of Neurology and Psychiatry, Tokyo, Japan

Patients with Parkinson's disease (PD) exhibit distinct abnormal postures, including neck-down, stooped postures, and Pisa syndrome, collectively termed “abnormal posture” henceforth. In the previous study, when assuming an upright stance, patients with PD exhibit heightened instability in contrast to healthy individuals with disturbance, implying that abnormal postures serve as compensatory mechanisms to mitigate sway during static standing. However, limited studies have explored the relationship between abnormal posture and sway in the context of static standing. Increased muscle tone (i.e., constant muscle activity against the gravity) has been proposed as an underlying reason for abnormal postures. Therefore, this study aimed to investigate the following hypothesis: abnormal posture with increased muscle tone leads to a smaller sway compared with that in other postures, including normal upright standing, under the sway minimization criterion. To investigate the hypothesis, we assessed the sway in multiple postures, which is determined by joint angles, including cases with bended hip joints. Our approach involved conducting forward dynamics simulations using a computational model comprising a musculoskeletal model and a neural controller model. The neural controller model proposed integrates two types of control mechanisms: feedforward control (representing muscle tone as a vector) and feedback control using proprioceptive and vestibular sensory information. An optimization was performed to determine the posture of the musculoskeletal model and the accompanied parameters of the neural controller model for each of the given muscle tone vector to minimize sway. The optimized postures to minimize sway for the optimal muscle tone vector of patients with PD were compared to the actual postures observed in these patients. The results revealed that on average, the joint-angle differences between these postures was <4°, which was less than one-tenth of the typical joint range of motion. These results suggest that patients with PD exhibit less sway in the abnormal posture than in other postures. Thus, adopting an abnormal posture with increased muscle tone can potentially serve as a valid strategy for minimizing sway in patients with PD.

1. Introduction

Parkinson's disease (PD) is a neurodegenerative disorder that primarily affects patients aged >60 years. The incidence of PD increases alongside the aging population (Van Den Eeden et al., 2003; Poewe et al., 2017). The motor symptoms of PD encompass rigidity and postural instability (Balestrino and Schapira, 2020), with the former being linked to increased muscle tone, which is continuous muscle activity against gravitational forces (Doherty et al., 2011; Ivanenko and Gurfinkel, 2018).

Patients with PD exhibit a distinct abnormal posture during standing, typified by conditions such as neck-down and stooped postures and Pisa syndrome (Doherty et al., 2011). These postures often cause back pain and compromise lung capacity, exerting a substantial influence on overall quality of life of patients (Bloch et al., 2006; Doherty et al., 2011). Moreover, these abnormal postures vary among patients with PD, presenting a range of postural configurations collectively categorized as abnormal posture. The pathology underlying these abnormal postures remains ambiguous, prompting the formulation of multiple hypotheses. One prevailing hypothesis points toward rigidity (Doherty et al., 2011). Studies investigating the interplay between the truncal muscles and motor symptoms have reported that hypertonia within the truncal muscles could potentially contribute to postural deficits (Wright et al., 2007). Additionally, studies have shown that lidocaine administration to the external oblique muscle ameliorates abnormal postures, thereby implying that dystonia within the truncal flexor muscles plays a role in the emergence of these abnormal postures (Furusawa et al., 2015). These results underscore a conceivable connection between abnormal posture and increased muscle tone.

Notably, a relationship between abnormal postures and postural sway during standing has also been reported. Bloem et al. proposed that an abnormal posture serves as a compensatory mechanism to avert backward falls by stabilizing the standing posture against foot rotation, such as toe lifting (Bloem et al., 1999). Dietz et al. examined the sway of patients with PD standing on a treadmill who were subjected to anteroposterior sinusoidal movement and reported that the patients exhibited increased sway in upright postures compared with the control subjects, and the upper body of the former failed to adequately adapt to the imposed vibration (Dietz et al., 1993). Conversely, Jacobs et al. observed a reduction in the stability margin [i.e., the difference between the center of mass (CoM) and center of pressure (CoP)] when healthy individuals assumed an abnormal posture, suggesting that individuals who adopt an abnormal posture are more prone to falls when exposed to disturbances (Jacobs et al., 2005). This is inconsistent with the idea that abnormal posture compensates for fall prevention (Jacobs et al., 2005). However, in light of the aforementioned experiments conducted by Dietz et al., it was determined that under perturbed conditions, patients with PD in an upright posture exhibit greater sway than healthy individuals standing in an upright posture (Dietz et al., 1993). Thus, patients with PD may be able to stand with less sway in an abnormal posture than in an upright posture.

Minimizing sway is one of the objectives of human postural control during static standing. To achieve this, one approach involves adjusting the positions of the CoP and CoM to remain within the base of support (Massion et al., 2004; Horak, 2006). However, even when the position of the CoM is within the base of support, a high CoM velocity can lead to instability and potential falling. Therefore, CoM velocity is equally important to keeping the CoM within the base of support to maintain balance while standing (Pai and Patton, 1997). Based on these findings, in the case of patients with PD exhibiting increased muscle tone, sway in abnormal postures (particularly CoM velocity) may be less pronounced than in other postures. However, studies regarding abnormal postures and sway during static standing are limited, and the aforementioned studies only focused on abnormal posture and sway under disturbance.

Given that patients with PD often display abnormal postures even in the absence of external perturbation while static standing, it is important to uncover the relationship between such abnormal postures and postural sway during static standing to understand abnormal postures in depth and apply corrective measures. Hence, this study aimed to scrutinize the hypotheses concerning the association between abnormal postures and postural sway during static standing. As previously noted, existing research has hinted at a connection between abnormal posture, increased muscle tone, and postural sway (Dietz et al., 1993; Bloem et al., 1999; Jacobs et al., 2005; Doherty et al., 2011). Therefore, this study examined the following hypothesis: abnormal posture with increased muscle tone leads to smaller sway compared with that in other postures, including normal upright standing, under the sway minimization criterion.

2. Methods

2.1. Overview

To assess the abovementioned hypothesis, it is essential to evaluate the extent of postural sway across diverse standing postures determined via the combinations of joint angles in the presence of increased muscle tone. Herein, we performed forward dynamics simulations involving a computational model encompassing a musculoskeletal model representing the human body and a neural controller model responsible for regulating the musculoskeletal model while also representing muscle tone. This model enables the determination of postural sway magnitude for various given sets of parameters, including overall body muscle tone, posture of the musculoskeletal model (joint angle), and parameters of the neural controller model.

While muscle tone might be elevated in patients with PD, the precise measurement of comprehensive whole-body muscle tone in such patients is challenging. To tackle this concern, the concept of muscle tone vector is introduced, which represents the muscle tone across the entire body, to accomplish the study objective through the following approaches:

(1) We calculated the muscle tone vector for each patient with PD so that the simulated posture and sway could reproduce the experimental posture and sway.

(2) We compared the simulated postures at the calculated muscle tone vectors with the experimental postures.

2.2. Problem setting

2.2.1. Tasks

Herein, we focused on the abnormal postures and postural sways observed in patients with PD during static standing. Consequently, the objective assigned to the computational model is to maintain a static standing posture for 5 s. Scaling and forward dynamics simulations were performed on the musculoskeletal simulators, OpenSim (Delp et al., 2007; Seth et al., 2018) and SCONE (Hyfydy) (Geijtenbeek, 2019, 2021).

2.2.2. Attempts

We made the following assumptions about human postural control.

• One of the objectives of human postural control is minimizing of the CoM velocity during static standing.

The aims of human postural control encompass the manipulation of both the positions of the CoM and CoP to ensure that they are within the base of support (Massion et al., 2004; Horak, 2006). Studies have highlighted the significance of the position and velocity of CoM as an escalation in the velocity poses a challenge to maintaining equilibrium (Pai and Patton, 1997). Furthermore, it has been proposed that among the CoM position, velocity, and acceleration utilized for stabilizing human postures during static standing, CoM velocity is the most precise indicator (Jeka et al., 2004). This perspective underscores the vital role of CoM velocity within postural control in static standing. Accordingly, the objectives of human postural control include maintaining the CoM position within the base of support and minimizing the CoM velocity.

2.3. Computational model

This section delineates the computational model utilized in the present study. An overview of the model is depicted in Figure 1. Further details concerning the computational model can be found in our previous publication (Omura et al., 2022).

FIGURE 1
www.frontiersin.org

Figure 1. Computational model (Omura et al., 2022). It comprises musculoskeletal and neural controller models involving feedforward control, feedback control, and time delay. uff(t), ufb(t), and u(t) are vectors representing the signals of the FF and FB controls and their sum for all the muscles at t seconds, respectively. Kp, Kd, and Kves are the FB gains. LMT(t) and L˙MT(t) are the muscle length and stretch rate at t seconds, respectively, with subscript 0 representing the initial value. eves is a vector representing vestibular sensory information. sig is 1 for an extensor muscle and −1 for a flexor muscle, providing an excitatory signal to the extensor and an inhibitory signal to the flexor. a(t) is the muscle activity at t seconds.

2.3.1. Musculoskeletal model

We utilized a musculoskeletal model encompassing 21 degrees of freedom (DOF) and 94 muscles, as proposed in our study, to represent the human body (Omura et al., 2022). Typical abnormal postures among patients with PD encompass hip and knee flexion, while severely abnormal posture also includes neck flexion (Doherty et al., 2011). In conjunction with these DOFs, the musculoskeletal model incorporates joint DOF for the back and ankles. The musculoskeletal model with these joint DOF can reproduce the abnormal postures of patients with PD. The musculoskeletal model was scaled to the motion capture data described below.

2.3.2. Neural controller model

The neural controller model proposed in our study was used to control the musculoskeletal model (Omura et al., 2022). In the previous study, this neural controller model was utilized to replicate human postural sway, and the reproduced sway was valid as that of humans (Omura et al., 2022). The neural controller model comprises an FF control mechanism to model muscle tones and an FB control mechanism employing proprioceptive and vestibular information and accounts for time delays as well. These control mechanisms were implemented based on the reticulospinal and vestibulospinal tracts, which are important for human postural control (Takakusaki, 2017; Omura et al., 2022). The FB controller operates with a time delay owing to multiple sources of delay, including synaptic and muscle-activation dynamics as well as signal transmission along sensory fibers (τfb). The summation of time delay is up to 120 ms (Zajac, 1989; Winters, 1995; Masani et al., 2006). The output of the neural controller model directed to the ith muscle is expressed as follows (Jiang et al., 2016; Omura et al., 2022).

ui=uff,i+ufb,i    (1)
ufb,i=ufb,prop,i+ufb,ves,i    (2)
ufb,prop=Kp,iLiMT(t-τfb)-Li,0MTLi,0MT+Kd,iL˙iMT(t-τfb)-L˙i,0MTVmax    (3)
ufb,ves,i=sigKves,ieves    (4)

uff,i and ufb,i are the FF and FB outputs to the ith muscle, respectively. ufb,prop,i and ufb,ves,i are the FB outputs using proprioceptive and vestibular information to the ith muscle, respectively. Kp,i and Kd,i are the FB gains of the ith muscle, respectively. Kves,i is the vector of the FB gain vector of the ith muscle, respectively. LiMT and L˙iMT are the length and lengthening speed of the ith muscle. Li,0MT and L˙i,0MT are the initial muscle length and lengthening speed of the ith muscle. τfb is the time delay due to the signal transmission along sensory fibers (τfb = 40ms). The sig outputs an excitatory signal to the extensors and an inhibitory signal to the flexors by being 1 for extensors and −1 for flexors. eves is the vector that represents the deviations of vestibular information from the target values.

Let uff, a vector with uff,i in Equation (1) as each element, be a muscle tone vector. The muscle tone vector can have various values even for the same posture. uff2, the square of the norm of the muscle tone vector, represents the index of the height of the muscle tone of the whole body, which was obtained from the previous study for comparison (Jiang et al., 2016). A large value of this index indicates a high muscle tone throughout the body. ufb,prop,i is the output of the FB control using the muscle length and lengthening speed of the musculoskeletal model as the FB information to the ith muscle. ufb,prop,i is the difference from the initial value of each FB information (Li,0MT and L˙i,0MT, which are the initial muscle length and lengthening speed of the ith muscle, respectively) multiplied by the FB gains, as shown in Equation (3). ufb,ves,i is the output of the FB control using the position information (position, angle, velocity, angular velocity, acceleration, and angular acceleration) of the head of the musculoskeletal model as the FB information. Moreover, ufb,ves,i is the difference from the target value of each FB information multiplied by the FB gains, as expressed in Equation (3). See our previous papers for details (Omura et al., 2022).

2.4. Posture calculation

This section outlines the approach employed to identify postures using a computational model. Specifically, our approach encompasses the following steps: (i) calculating candidates for muscle tone vectors with varying norms corresponding to the postures of patients with PD and (ii) for each candidate of muscle tone vector, determining the posture of the musculoskeletal model and parameters of the neural controller model. These parameters are calculated to ensure a static standing with minimal sway, which was accomplished through an optimization method. The methodological overview is visually depicted in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2. Flowchart of the entire method. The left part of the figure shows the outline of the entire method of the data of one participant. (A) Outline of the muscle tone vector estimation. (B) Outline of the parameter adjustment procedure. The blue part shows the procedure for the feedback gain adjustment, and the red part shows the procedure for the posture (target posture) adjustment.

The optimization method employed was the covariance matrix adaptation evolution strategy (CMA-ES), according to the previous study (Jiang et al., 2016). This method can be used to find postures with less noticeable swaying at different muscle tone vectors.

2.4.1. Measurement of the posture of patients with PD

We utilized marker data derived from the standing posture of patients with PD measured via a motion capture system.

The analysis included eight patients with PD (4 men and 4 women, aged 66.9 ± 8.4 years). Participant information is present in Table 1. We attached 31 or 59 markers to the entire body of the participants. The measurements were conducted during static standing for 5 s using a motion capture system (6-camera, Oqus 3+/Oqus 5+, Qualisys, Gothenburg, Sweden). One trial involved the patient standing static for 5 s. Each participant performed this trial twice or thrice. The studies involving human participants were reviewed and approved by the Ethics Committee at the National Center of Neurology and Psychiatry (A2019-126). The patients provided their written informed consent to participate in this study. We analyzed the first trial with no measurement defects among the two or three trials. Additionally, in the analysis of static standing, the measurement data of the initial 4 s of the entire 5 s were considered owing to some deficiencies in the last 1 s, such as movement exhibited by several participants. In the case of Sub3, the marker position measurements were incomplete during the first 0.25 s, thereby necessitating the exclusion of this segment from our analysis.

TABLE 1
www.frontiersin.org

Table 1. Data obtained from patients with PD.

To scale the musculoskeletal model, a set of 23 or 27 marker data points was employed. These markers encompassed the anterior and posterior superior iliac spines, greater trochanter, medial and lateral femoral epicondyles, medial and lateral malleoli, calcanei, 1st and 5th metatarsal heads, acromion processes, head, C7 vertebrae, TH12, episternum, and xiphoid process of the sternum. To scale the musculoskeletal model for each patient, the postures of the patients were replicated using a 2 s duration of marker data. This choice of a 2 s for scaling was motivated by the potential for error escalation if the posture lacked stability and encompassed swaying. For scaling and reproduction, we used the scaling and inverse kinematics (IK) tools of OpenSim (Delp et al., 2007; Seth et al., 2018). Scaling was executed by adapting the distances between markers on the musculoskeletal model according to the corresponding markers on patients with PD. The subsequent computation of body weight relied on the ideal body weight calculations (BMI: 22) (Tokunaga et al., 1991) based on the height derived from the scaling outcome. Employing IK alongside marker data, the joint angles of the musculoskeletal model were computed. This ensured congruence between every marker on the model and their corresponding counterparts in patients with PD. By adopting this methodology, the sway exhibited by each patient with PD was reproduced on the musculoskeletal model.

The posture replicated through IK involved issues related to ground contact due to disparities between the marker positions within the musculoskeletal model and those observed in patients with PD. This discrepancy led to instances where the model's feet appeared to hover or were embedded within the floor. These challenges had adverse implications; for instance, they precipitated occurrences where the musculoskeletal model either collapsed or instantly emerged at the onset of forward dynamics simulations. Consequently, the joint angles of the lower limbs were adjusted via an optimization method to rectify the ground contact issue between the feet and the floor of the posture reproduced by IK. The following evaluation function sought to minimize alterations in the CoM position of the feet to obtain an appropriate posture. The optimization method employed for this purpose was the CMA-ES, akin to its application in the adjustment of the posture and FB gains within the neural controller model, which has been elaborated upon subsequently. CMA-ES is an effective evolution strategy for non-linear and non-convex functions that do not require derivatives (Hansen, 2016). Therefore, CMA-ES was employed to adjust the posture and FB gains, which are non-linear and have derivatives that are difficult to calculate.

Jpos=wforceJforce+wCOMposJCOMpos+wfeetvelJfeetvel    (5)
Jforce=00.1(|fr(t)-fr(0)|+|fl(t)-fl(0)|)dt    (6)
JCOMpos=|xCOM,l-xCOM,exp,l|+|xCOM,r-xCOM,exp,r|    (7)
Jfeetvel=(00.1|vankle,r(t)|+|vankle,l(t)|+|vsubtalar,r(t)|            +|vsubtalar,l(t)|)dt    (8)

Jforce is a term that evaluates the magnitude of the floor reaction force. f(t) denotes the floor reaction force in t seconds. The subscripts r and l denote the right and left foot, respectively. JCOMpos is a term that evaluates the CoM position relative to the foot. xCOM and xCOM,exp denote the CoM position relative to the foot at 0 s and the CoM position relative to the foot as reproduced by IK, respectively. Jfeetvel is a term that evaluates the angular velocity of the joints around the foot. vankle(t) and vsubtalar(t) denote the angular velocities of the ankle and subtalar joints in t seconds, respectively.

2.4.2. Calculation of candidates for muscle tone vector

To determine the candidates of the muscle tone vector associated with the measured posture of each patient with PD, we first aligned the posture of the musculoskeletal model to replicate that of the patient. The candidates for the muscle tone vector were derived using a musculoskeletal model, scaled to the posture of a patient with PD. The procedural framework is depicted in Figure 2A. Muscle tone indicates constant muscle activity and is independent of time delay. Therefore, the experimental posture was set to the posture of the musculoskeletal model and the posture was maintained for 5 s under the condition without time delay and with FB control only, based on previous studies (Jiang et al., 2016). From these results, the average muscle activity of each muscle in the second half was calculated as the muscle tone. The muscle tone was obtained using the following equation (Jiang et al., 2016):

uff,i=t1t2ai(t)dtt2-t1    (9)

uff,i denotes the muscle tone in the ith muscle. ai(t) denotes the muscle activity of the ith muscle in t seconds. t1 and t2 are 3 and 5 s, respectively, based on the previous study (Jiang et al., 2016). The muscle tone vector is given in the following equation:

uff=[uff,1uff,2uff,94]    (10)

The muscle tone vectors that could maintain the standing posture for 5 s for each given posture were not unique and different muscle tone vectors could do the same job depending on the values of FB gains used for the simulations. Furthermore, the muscle tone vectors procured might not necessarily ensure the model's ability to maintain a standing posture under circumstances involving time delays. Moreover, the task of measuring muscle tone throughout the entire body of actual patients with PD presents inherent challenges. Therefore, multiple muscle tone vectors need to be calculated as candidates, from which we can narrow down a list of muscle tone vectors to the “optimal” one.

To calculate candidates for muscle tone vector, the given FB gains were determined via the CMA-ES. It is possible to calculate muscle tone vectors that have different values by incorporating the cosine similarity with the already calculated muscle tone vectors into the evaluation function. To calculate the muscle tone vector with arbitrary uff2, a term evaluating uff2 is added to the evaluation function. In addition, two terms for evaluating falls and the CoM velocity based on assumptions were incorporated. Thus, the evaluation function is expressed as follows:

Jmt=wfallJfall+wcomJcom+wcosJcos+wuffJuff    (11)
Jfall=1Tfail(Tsimu-Tfail)    (12)
Jcom=(vcom, x(t)2+vcom, y(t)2+vcom, z(t)2)dtTfall    (13)
Jcos=1ni=1nai·baib    (14)
Juff=b||2-α    (15)
w uff={1 (|Juff|<0.5)1000 (|Juff|0.5)    (16)

Jfall is a term that evaluates the time to fall. Tfail and Tsimu denote the time to fall and the simulation time (5 s), respectively. Jcom is a term that evaluates the CoM velocity. vcom,x(t), vcom,y(t), and vcom,z(t) denote the CoM velocities in the x, y, and z axes in t seconds, respectively. Jcos is a term that evaluates the cosine similarity. n denotes the number of muscle tone vectors calculated so far. ai denotes the ith calculated muscle tone vector. b denotes the muscle tone vector to be calculated. Juff is a term that evaluates uff2 of the muscle tone vector. α denotes uff2 of the muscle tone vector to be calculated. wfall, wcom, wuff and wcos are 50,000,000, 5,000, 1,000, and 0.1, respectively, based on previous research, to search for solutions that do not fall and reduce the CoM speed preferentially (Kaminishi et al., 2019). wuff was considered according to Equation (16), and the value of the term to evaluate the cosine similarity relatively increased when the muscle tone vector was calculated to ensure that the difference between the specified value and uff2 was < 0.5. Consequently, muscle tone vectors with varying cosine similarities can be preferentially searched when the condition of uff2 is fulfilled.

In a previous study employing a comparable computational model, muscle tone vectors were computed to achieve equidistant values of uff2. The resulting muscle tone vectors were then compared with the experimental results of healthy individuals. This analysis unveiled that the muscle tone vector corresponding to uff2=2.07 yielded muscle activity most akin to that of healthy individuals in the upright stance (Jiang et al., 2016). Herein, muscle tone vectors were calculated for uff2 values exceeding uff2=2.07. Through the adjustment of α in Equation 15, specifically setting it at 2, 4, 6, and 8, the computation encompassed the derivation of muscle tone vectors with uff2 values lying within a range of ±0.5 from these designated magnitudes. To calculate muscle tone vectors with multiple directions, we repeatedly computed the muscle tone vector until n = 6 at Equation (14) empirically (Figure 2A).

2.4.3. FB gains and posture calculation

Using the candidates of muscle tone vectors in the postures of patients with PD that were calculated in the previous section, we searched for a posture in which the musculoskeletal model can maintain a standing posture with minimal sway under conditions with a time delay. To achieve this, the parameters of the computational model were adjusted using the optimization method. The procedure is shown in Figure 2B. The target parameters are the FB gains and posture (target posture for FB control). We utilized the posture data of a patient with Parkinson's Disease as the starting point. From there, we used forward dynamics simulation to recreate sway and evaluated the recreated sway by adjusting the parameters accordingly. The FB gains were adjusted to maintain the standing posture despite time delays, and the posture of the musculoskeletal model was modified to find a posture with less sway. CMA-ES was used as the optimization method. The evaluation function is expressed as the following function that evaluates the CoM velocity based on the assumption.

Jparam=wfallJfall+wcomJcom    (17)

The FB gains encompass a total of 280 parameters, while posture involves 21 parameters in total. Simultaneously adjusting these parameters would lead to an expansive solution space owing to their high count.

To navigate this challenge, these parameters are alternately adjusted by tuning the FB gain and posture. The optimization method seeks to find the optimal values for these parameters. In the FB gain optimization, the optimal value of the evaluation function is denoted as J1, whereas in the optimization of the posture, it is denoted as J2. Each optimization iteration continues until J1 becomes < J2 during FB gain optimization or until J2 becomes < J1 during posture optimization.

We adjusted the parameters for each muscle tone vector that was calculated. To evaluate the posture and sway, we used the muscle tone vector with the smallest evaluation in each uff2. This resulted in obtaining simulation results with optimized parameters for 32 conditions, which included 8 participants and 4 muscle tone vectors each.

2.5. Evaluation

The resulting posture and sway from the simulation using the optimized parameters were evaluated using the following procedure. The procedures correspond to the approach (1) and (2) in Section 2.1..

(1) We selected the muscle tone vector for each patient with PD so that the simulated posture and sway could reproduce the experimental posture and sway.

(2) We compared the simulated postures at the calculated muscle tone vectors with the experimental postures.

To assess and compare the standing posture and sway as described in (1), several metrics were calculated: (a) mean difference across all joint angles, (b) difference in the mean CoM velocity, and (c) standard deviation of all the joint angles. The evaluation of differences in the standing posture was quantitatively performed via (a), while (b) and (c) were employed to evaluate sway. The difference in the mean CoM velocity indicated by metric (b) was used to assess sway, grounded in certain assumptions. However, (b) might be inadequate to gauge the sway of individual body parts such as in cases where two body parts sway in opposing phases. Consequently, (c) was introduced to evaluate the sway of each specific body part.

As these metrics possess distinct units, each value was standardized using the T-score. The aggregate of squared T-scores was employed as an indicator of posture and sway. The muscle tone vector that yields the smallest index value for the combined posture and sway evaluation can be deemed the most optimal. Herein, we considered it as the muscle tone vector for the patients with PD, established through the optimization method. During step (2), we used (a) to compare the standing posture at the selected muscle tone vector with the experimental posture.

3. Results

3.1. Posture adjustment

Figure 3 shows the patient postures reproduced from the measured marker data using IK along with the outcomes of the posture adjustment process described in 2.5. Notably, Subs 1 and 4 exhibited pronounced trunk flexion. For detailed variations in joint angles before and after posture adjustment, please refer to the Appendix in Supplementary material.

FIGURE 3
www.frontiersin.org

Figure 3. The abnormal posture of each participant before and after postural adjustment. Different postures are represented for each participant.

3.2. Parameter adjustment

For each participant, six potential candidates were computed for muscle tone vectors corresponding to α = 2, 4, 6, and 8. In other words, we could calculate the candidates for muscle tone vectors with distinct elements using cosine similarities.

For every muscle tone vector, the parameters were adjusted according to Section 2.4.3. Consequently, across all uff2 values for each participant, the musculoskeletal model was able to maintain the standing posture for a minimum of 5 s employing at least one muscle tone vector. This suggests that the posture and FB gain necessary for maintaining a standing posture with time delays can be determined through parameter adjustment using the optimization method.

3.3. Muscle tones

The index for assessing posture and sway, as detailed in Section 2.5, was computed for each uff2 condition. The indices to evaluate posture and sway under each condition are depicted in Figure 4. Among the participants, the muscle tone vectors yielding the smallest index were uff2=4 for Sub2, uff2=6 for Subs1, 4, 6, 7, and 8, and uff2=8 for Subs3 and 5. These results indicate that the postural sway in these muscle tone vectors closely matches the experimental results. In other words, the posture and sway observed in these muscle tone vectors were the most representative of actual patients with PD. Previous studies utilizing computational models (Jiang et al., 2016) have suggested that a muscle tone vector with uff2 = 2.07 is most fitting for healthy subjects. Notably, uff2 = 4, 6, 8 exceed uff2 = 2.07. Consequently, within the set of calculated muscle tone vectors (uff2 = 2, 4, 6, 8), the muscle tone vector with a larger uff2 than that of healthy individuals yielded the most suitable posture and sway for patients with PD. These outcomes are consistent with the observation that muscle tone is elevated in patients with PD compared with that in healthy individuals (Burleigh et al., 1995).

FIGURE 4
www.frontiersin.org

Figure 4. Index to evaluate the posture and sway under each condition. uff2 is a measure of the muscle tone vector. In each condition, the differences in the mean CoM velocities, standard deviations of the joint angles, and mean T-score of the differences in the joint angles are shown in different colors.

3.4. Comparing postures

Figure 5 demonstrates the posture at the muscle tone vector corresponding to the minimum index for evaluating posture and sway in Figure 4, juxtaposed with the posture from experimental results. This depiction presents postures within the sagittal and frontal planes. The simulation result with optimized parameters represents the average standing posture from 1 to 4 s.

FIGURE 5
www.frontiersin.org

Figure 5. Comparison of the standing posture in the simulation results obtained using optimized parameters and the posture in the experiment results at the muscle tone vector where the index to evaluate the posture and sway in each participant was the smallest. The simulation result shows the average posture in 1–4 s. (A) Sagittal plane posture, and (B) frontal plane posture.

Furthermore, Table 2 provides mean values ± standard deviation (SD) of the joint angle differences for these postures. The comparison of these postures revealed that the mean joint-angle difference was <4° in all the participants except Sub4 where it was 5.2°. Calculating the mean + 2SD corresponding to a 95% confidence interval under the assumption of a normal distribution of joint angle differences revealed a mean + 2SD of 5.5° for Sub1 and 8.8° for Sub4.

TABLE 2
www.frontiersin.org

Table 2. The mean value of all the joint-angle differences between the simulation results with optimized parameters and the experimental results at the muscle tone vector where the index to evaluate the posture and sway was the smallest ± standard deviations.

To summarize, a quantitative assessment of posture differences using joint angles (Table 2) displayed a mean difference of <4° in all the participants except Sub4. Additionally, the mean difference +2SD was <4° for all the subjects except Sub1 and Sub4. Subsequently, a comparison was conducted between the joint angle differences and range of motion. A characteristic abnormal posture of patients with PD often involves flexed knees and hips (Doherty et al., 2011). Consequently, the discrepancy in joint angle and range of motion was compared in the sagittal plane for the knee and lumbar joints, as representative examples. The maximum difference in joint angle was 5.2° for the knee joint and 5.0° for the lumbar joint, except for in Sub4 (Appendix in Supplementary material). The ranges of motion for the knee and lumbar joints were approximately 130° and 90°, respectively (Roach and Miles, 1991; Dvořák et al., 1995). In essence, the joint angle difference stood at merely 4.0% of the range of motion for the knee joint and 5.6% of that for the lumbar joint. Consequently, even in comparison to the range of motion, the joint angle difference between the simulation results and experimental outcomes using optimized parameters remained at < 1/10 of the joint motion range.

4. Discussion

We selected the muscle tone vectors identified in Section 3.3 as the ones responsible for actuating the joints of patients with PD, and proceeded to compare the resulting postures. As detailed in Section 3.4, the mean differences in the joint angles between the postures optimized to minimize sway and the postures of patients with PD were ≤ 1/10 of the range of joint motion. This outcome indicates that the extent of sway is smaller in the abnormal posture than in other postures.

Conversely, the mean + 2SD of the joint angle exceeded 4° only in Sub1 and Sub4. These variations are likely attributable to differences in disease severity among the patients. For instance, Sub1 exhibited the longest disease duration and a Hoehn & Yahr (H&Y) scale score of 3. Similarly, Sub4 exhibited the highest Unified Parkinson's Disease Rating Scale (UPDRS) Part III score compared with those of the other participants, as indicated in Table 1. These observations point to greater overall PD severity in Sub1 and more severe motor symptoms in Sub4 compared with that in the other participants. Furthermore, as illustrated in Figure 3, Subs 1 and 4 demonstrated more pronounced abnormal postures than the other participants. Reportedly, patients with advanced PD might experience alterations in soft tissue, muscles, and even the spinal cord, potentially leading to more pronounced postural deformities (Doherty et al., 2011). Thus, it is conceivable that, for Subs 1 and 4, maintaining a specific posture for an extended period could have contributed to skeletal deformation, thereby hindering their ability to assume different postures. The omission of skeletal deformation as a consideration in our study might have influenced the disparity between the experimental results and the simulation outcomes obtained with optimized parameters.

Herein, we calculated multiple muscle tone vectors based on the observed experimental postures. Additionally, the FB gains and postures of our computational model were adjusted using an optimization method for various muscle tone vectors. This adjustment allowed the musculoskeletal model to sustain the standing posture, and we sought to identify a posture with reduced sway starting from the experimental posture. By comparing the simulation results achieved using the optimized parameters to the experimental outcomes, we found that the muscle tone vectors with uff2 values surpassing those of healthy participants yielded the closest resemblance to the experimental data. In essence, the posture and sway achieved through these larger uff2 muscle tone vectors best replicate the characteristics of patients with PD. This observation aligns with prior research that highlights increased muscle tone among patients with PD compared with their healthy individuals (Burleigh et al., 1995). Furthermore, the disparities in joint angles between the optimized postures (which minimized sway) at these muscle tone vectors and the experimental postures were < 1/10 of the joint range of motion. These findings suggest that patients with PD who exhibit increased muscle tone experience reduced sway in their abnormal postures compared with other stances. Consequently, to effectively mitigate sway in the postural control of patients with PD, it might be crucial to maintain an abnormal posture with their increased muscle tone across the entire body. Therefore, the hypothesis of abnormal posture with an increased muscle tone leading to a smaller sway compared with that in other postures, including normal upright standing, under the sway minimization criterion may be effective for postural control in patients with PD.

Prior research has indicated variations in muscle tone among different patients, with some exhibiting higher muscle tone while others exhibiting lower muscle tone levels (Horak et al., 1992; Yamamoto et al., 2011). The findings of our study revealed that in all the participants, the muscle tone observed was generally higher than that in healthy individuals. Therefore, the results of this study may represent a strategy for patients with higher muscle tone. As mentioned in the Introduction section, prior experiments involving lidocaine administration to the external oblique muscle reported improvements in abnormal postures (Furusawa et al., 2015). In the same experiment, as the tilt of the tilt table approached 90°, abnormal posture became more pronounced alongside an increase in muscle tone (Furusawa et al., 2015). In addition, a model-based categorization of patients with PD revealed that the strategies they employ differ based on ankle flexibility, which is closely associated with muscle tone levels (Suzuki et al., 2020). Therefore, patients with lower muscle tone may not exhibit abnormal postures as prominently.

4.1. Limitation

Herein, we computed muscle tone vectors based on measured abnormal postures and used an optimization method to determine postures that minimized sway compared with other postures. As a result, the postures to minimize sway corresponded to abnormal postures. Therefore, the muscle tone vectors calculated for patients with PD in this study reflect abnormal postures. The precise reasons for the differences in muscle tone vectors between patients with PD and healthy individuals were not addressed within this study. Additional research into alterations in muscle tone vectors is warranted.

As previously mentioned, musculoskeletal anomalies have been associated with abnormal posture, particularly in patients with advanced or late-stage PD (Doherty et al., 2011). Herein, we replicated the postures of patients with PD using a musculoskeletal model. However, considering musculoskeletal deformities in patients who exhibited such postures was not feasible within our study. Gathering comprehensive skeletal system data for each participant to accurately simulate such deformities poses significant challenges. Nevertheless, a more comprehensive exploration of abnormal posture, especially in patients with later-stage of PD, could be achieved by focusing on established musculoskeletal deformities and conducting simulations based on these assumptions.

We focused on the static standing posture and did not delve into dynamic posture control, such as walking or walking initiation. Previous research has suggested that abnormal posture functions as a compensatory mechanism to enhance CoM movement at the beginning of walking (Jacobs et al., 2005). Investigating the influence of these postures on dynamic postural control would be valuable, and simulations involving actions such as walking initiation could provide further insights into abnormal postures.

We formulated the evaluation function in this study based on the assumption that minimizing the CoM velocity is essential during human standing posture control. However, other evaluation functions such as energy and muscle activity minimization are also considered in human posture control, but we did not use them in this study. Nevertheless, considering fall prevention, the reduction of CoM velocity is crucial, thereby making the evaluation function we used suitable for exploring our hypothesis.

Considering that abnormal postures occur because of increased muscle tone, it may be possible that the higher muscle tone leads to more abnormal postures. However, abnormal postures can manifest in diverse ways, encompassing factors such as foot width and body leaning. Consequently, defining the extent of postural abnormalities either qualitatively or quantitatively is challenging.

4.2. Conclusions and future perspectives

Herein, we examined the following hypothesis on abnormal posture during static standing: abnormal posture with increased muscle tone leads to smaller sway compared with that in other postures, including normal upright standing, under the sway minimization criterion. We used a computational model comprising musculoskeletal and neural controller models to test the proposed hypothesis. The norm of the muscle tone vector that could reproduce the most appropriate experimental posture and sway for patients with PD was larger than that for healthy individuals, consistent with the increased muscle tone in patients with PD. The average difference in the joint angles between the experimental posture and the posture to minimize sway was <4° except for one participant. The results of the computational model suggested that patients with PD with increased muscle tone display less sway in abnormal postures than in other postures. These findings suggest that once the muscle tone of patients with PD becomes that at an anterior leaning posture, they need to maintain an abnormal posture because their sway is greater in other postures compared with that in the abnormal posture. Given the importance of sway reduction in human postural control, it may be desirable to maintain abnormal postures more than the other postures and it may be necessary to alter muscle tone to change the standing posture. Therefore, the hypothesis described in this study is considered valid. In future studies, we intend to use the computational model used in this study to explore the interplay between walking initiation, disturbances, and abnormal posture. Additionally, we aim to delve into the connection between dopamine, the underlying cause of PD, and the parameters of the computational model.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by the Ethics Committee at the National Center of Neurology and Psychiatry. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

YO performed the simulations, analyzed the data, and drafted the manuscript. HT, MA, YT, and THan performed the experiments and acquired the data. HT, KK, THas, RC, AY, KT, THan, and JO reviewed the research methods. JO supervised the study. All authors discussed the results and reviewed the article. All authors contributed to the article and approved the submitted version.

Funding

This work was partially supported by JSPS KAKENHI, Grant-in-Aid for Scientific Research on Innovative Areas Hyper-adaptability for overcoming body-brain dysfunction: Integrated empirical and system theoretical approaches (Grant Numbers 19H05722 and 19H05730), Grant-in-Aid for JSPS Fellows (Grant Number 23KJ0522), and the Mohammed bin Salman Center for Future Science and Technology for Saudi–Japan Vision 2030 at The University of Tokyo (MbSC2030).

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.

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/fncom.2023.1218707/full#supplementary-material

References

Balestrino, R., and Schapira, A. H. V. (2020). Parkinson disease. Eur. J. Neurol. 27, 27–42. doi: 10.1111/ene.14108

CrossRef Full Text | Google Scholar

Bloch, F., Houeto, J. L., Tezenas du Montcel, S., Bonneville, F., Etchepare, F., Welter, M. L., et al. (2006). Parkinson's disease with camptocormia. J. Neurol. Neurosurg. Psychiatry 77, 1223–1228. doi: 10.1136/jnnp.2006.087908

CrossRef Full Text | Google Scholar

Bloem, B. R., Beckley, D. J., and Van Dijk, J. G. (1999). Are automatic postural responses in patients with Parkinson's disease abnormal due to their stooped posture? Exp. Brain Res. 124, 481–488.

PubMed Abstract | Google Scholar

Burleigh, A., Horak, F., Nutt, J., and Frank, J. (1995). Levodopa reduces muscle tone and lower extremity tremor in Parkinson's disease. Can. J. Neurol. Sci. 22, 280–285.

PubMed Abstract | Google Scholar

Delp, S. L., Anderson, F. C., Arnold, A. S., Loan, P., Habib, A., John, C. T., et al. (2007). OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans. Biomed. Eng. 54, 1940–1950. doi: 10.1109/TBME.2007.901024

PubMed Abstract | CrossRef Full Text | Google Scholar

Dietz, V., Zijlstra, W., Assaiante, C., Trippel, M., and Berger, W. (1993). Balance control in Parkinson's disease. Gait Post. 1, 77–84.

Google Scholar

Doherty, K. M., van de Warrenburg, B. P., Peralta, M. C., Silveira-Moriyama, L., Azulay, J.-P., Gershanik, O. S., et al. (2011). Postural deformities in Parkinson's disease. Lancet Neurol. 10, 538–549. doi: 10.1016/S1474-4422(11)70067-9

CrossRef Full Text | Google Scholar

Dvořák, J., Vajda, E. G., Grob, D., and Panjabi, M. M. (1995). Normal motion of the lumbar spine as related to age and gender. Eur. Spine J. 4, 18–23.

PubMed Abstract | Google Scholar

Furusawa, Y., Hanakawa, T., Mukai, Y., Aihara, Y., Taminato, T., Iawata, Y., et al. (2015). Mechanism of camptocormia in Parkinson's disease analyzed by tilt table-EMG recording. Parkinson. Relat. Disord. 21, 765–770. doi: 10.1016/j.parkreldis.2015.02.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Geijtenbeek, T. (2019). SCONE: open source software for predictive simulation of biological motion. J. Open Source Softw. 4, 1421. doi: 10.21105/joss.01421

CrossRef Full Text | Google Scholar

Geijtenbeek, T. (2021). The Hyfydy Simulation Software. Available online at: https://hyfydy.com

Google Scholar

Hansen, N. (2016). The CMA evolution strategy: a tutorial. arXiv preprint arXiv:1604.00772. doi: 10.48550/arXiv.1604.00772

CrossRef Full Text | Google Scholar

Horak, F. B. (2006). Postural orientation and equilibrium: what do we need to know about neural control of balance to prevent falls? Age Ageing 35(Suppl 2), ii7–ii11. doi: 10.1093/ageing/afl077

PubMed Abstract | CrossRef Full Text | Google Scholar

Horak, F. B., Nutt, J. G., and Nashner, L. M. (1992). Postural inflexibility in Parkinsonian subjects. J. Neurol. Sci. 111, 46–58.

PubMed Abstract | Google Scholar

Ivanenko, Y., and Gurfinkel, V. S. (2018). Human postural control. Front. Neurosci. 12, 171. doi: 10.3389/fnins.2018.00171

CrossRef Full Text | Google Scholar

Jacobs, J. V., Dimitrova, D. M., Nutt, J. G., and Horak, F. B. (2005). Can stooped posture explain multidirectional postural instability in patients with Parkinson's disease? Exp. Brain Res. 166, 78–88. doi: 10.1007/s00221-005-2346-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeka, J., Kiemel, T., Creath, R., Horak, F., and Peterka, R. (2004). Controlling human upright posture: velocity information is more accurate than position or acceleration. J. Neurophysiol. 92, 2368–2379. doi: 10.1152/jn.00983.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, P., Chiba, R., Takakusaki, K., and Ota, J. (2016). Generation of the human biped stance by a neural controller able to compensate neurological time delay. PLoS ONE 11, e0163212. doi: 10.1371/journal.pone.0163212

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaminishi, K., Jiang, P., Chiba, R., Takakusaki, K., and Ota, J. (2019). Postural control of a musculoskeletal model against multidirectional support surface translations. PLoS ONE 14, e0212613. doi: 10.1371/journal.pone.0212613

PubMed Abstract | CrossRef Full Text | Google Scholar

Masani, K., Vette, A. H., and Popovic, M. R. (2006). Controlling balance during quiet standing: proportional and derivative controller generates preceding motor command to body sway position observed in experiments. Gait Post. 23, 164–172. doi: 10.1016/j.gaitpost.2005.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Massion, J., Alexandrov, A., and Frolov, A. (2004). Why and how are posture and movement coordinated? in Progress in Brain Research, volume 143 (Netherland: Elsevier), 13–27.

Google Scholar

Omura, Y., Kaminishi, K., Chiba, R., Takakusaki, K., and Ota, J. (2022). A neural controller model considering the vestibulospinal tract in human postural control. Front. Comput. Neurosci. 16, 785099. doi: 10.3389/fncom.2022.785099

PubMed Abstract | CrossRef Full Text | Google Scholar

Pai, Y. C., and Patton, J. (1997). Center of mass velocity-position predictions for balance control. J. Biomech. 30, 347–354.

PubMed Abstract | Google Scholar

Poewe, W., Seppi, K., Tanner, C. M., Halliday, G. M., Brundin, P., Volkmann, J., et al. (2017). Parkinson disease. Nat. Rev. Dis. Primers 3, 17013. doi: 10.1038/nrdp.2017.13

CrossRef Full Text | Google Scholar

Roach, K. E., and Miles, T. P. (1991). Normal hip and knee active range of motion: the relationship to age. Phys. Ther. 71, 656–665.

PubMed Abstract | Google Scholar

Seth, A., Hicks, J. L., Uchida, T. K., Habib, A., Dembia, C. L., Dunne, J. J., et al. (2018). OpenSim: simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. PLoS Comput. Biol. 14, e1006223. doi: 10.1371/journal.pcbi.1006223

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, Y., Nakamura, A., Milosevic, M., Nomura, K., Tanahashi, T., Endo, T., et al. (2020). Postural instability via a loss of intermittent control in elderly and patients with Parkinson's disease: a model-based and data-driven approach. Chaos 30, 113140. doi: 10.1063/5.0022319

PubMed Abstract | CrossRef Full Text | Google Scholar

Takakusaki, K. (2017). Functional neuroanatomy for posture and gait control. Mov. Disord. 10, 1–17. doi: 10.14802/jmd.16062

PubMed Abstract | CrossRef Full Text | Google Scholar

Tokunaga, K., Matsuzawa, Y., Kotani, K., Keno, Y., Kobatake, T., Fujioka, S., et al. (1991). Ideal body weight estimated from the body mass index with the lowest morbidity. Int. J. Obes. 15, 1–5.

Google Scholar

Van Den Eeden, S. K., Tanner, C. M., Bernstein, A. L., Fross, R. D., Leimpeter, A., Bloch, D. A., et al. (2003). Incidence of Parkinson's disease: variation by age, gender, and race/ethnicity. Am. J. Epidemiol. 157, 1015–1022. doi: 10.1093/aje/kwg068

PubMed Abstract | CrossRef Full Text | Google Scholar

Winters, J. M. (1995). An improved muscle-reflex actuator for use in large-scale neuromusculoskeletal models. Ann. Biomed. Eng. 23, 359–374.

PubMed Abstract | Google Scholar

Wright, W. G., Gurfinkel, V. S., Nutt, J., Horak, F. B., and Cordo, P. J. (2007). Axial hypertonicity in Parkinson's disease: direct measurements of trunk and hip torque. Exp. Neurol. 208, 38–46. doi: 10.1016/j.expneurol.2007.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamamoto, T., Suzuki, Y., Nomura, K., Nomura, T., Tanahashi, T., Fukada, K., et al. (2011). A classification of postural sway patterns during upright stance in healthy adults and patients with Parkinson's disease. J. Adv. Comput. Intell. Intell. Informat. 15, 997–1010. doi: 10.20965/jaciii.2011.p0997

CrossRef Full Text | Google Scholar

Zajac, F. E. (1989). Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control. Crit. Rev. Biomed. Eng. 17, 359–411.

PubMed Abstract | Google Scholar

Keywords: Parkinson's disease, abnormal posture, computational model, musculoskeletal model, neural controller model

Citation: Omura Y, Togo H, Kaminishi K, Hasegawa T, Chiba R, Yozu A, Takakusaki K, Abe M, Takahashi Y, Hanakawa T and Ota J (2023) Analysis of abnormal posture in patients with Parkinson's disease using a computational model considering muscle tones. Front. Comput. Neurosci. 17:1218707. doi: 10.3389/fncom.2023.1218707

Received: 08 May 2023; Accepted: 18 September 2023;
Published: 05 October 2023.

Edited by:

Petia D. Koprinkova-Hristova, Institute of Information and Communication Technologies (BAS), Bulgaria

Reviewed by:

Tetsuro Funato, The University of Electro-Communications, Japan
Taishin Nomura, Osaka University, Japan

Copyright © 2023 Omura, Togo, Kaminishi, Hasegawa, Chiba, Yozu, Takakusaki, Abe, Takahashi, Hanakawa and Ota. 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: Yuichiro Omura, omura@race.t.u-tokyo.ac.jp

Disclaimer: 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.