ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 30 April 2025

Sec. Biomechanics

Volume 13 - 2025 | https://doi.org/10.3389/fbioe.2025.1541240

Beyond clinical scales: an observational study on instrumental gait analysis and biomechanical patterns in patients with Parkinson’s disease

Paolo De Pasquale&#x;Paolo De Pasquale1Mirjam Bonanno
&#x;Mirjam Bonanno1*Cristiano De MarchisCristiano De Marchis2Luca PergolizziLuca Pergolizzi1Antonino Lombardo FaccialeAntonino Lombardo Facciale1Giuseppe PaladinaGiuseppe Paladina1Maria Grazia MaggioMaria Grazia Maggio1Federica ImpellizzeriFederica Impellizzeri1Irene CiancarelliIrene Ciancarelli3Angelo Quartarone,Angelo Quartarone1,4Rocco Salvatore CalabrRocco Salvatore Calabrò1
  • 1IRCCS Centro Neurolesi Bonino-Pulejo, Messina, Italy
  • 2Department of Engineering, University of Messina, Messina, Italy
  • 3Department of Life, Health and Environmental Sciences, University of L'Aquila, L’Aquila, Italy
  • 4Department of Biomedical, Dental Sciences and Morphological and Functional Images, University of Messina, Messina, Italy

Introduction: Parkinson’s disease (PD), a common neurodegenerative disorder affecting motor functions, is associated with abnormal gait patterns characterized by altered kinematic, kinetic, and electrophysiological parameters. This observational study aims to instrumentally identify and quantify these gait dysfunctions in PD patients compared to normal values from healthy subjects.

Methods: Sixty-nine PD patients underwent clinical and instrumental evaluations to assess gait. Demographic and clinical data were collected before motor assessment. Clinical scales evaluated the level of impairment, gait, balance, risk of falls and ability to complete activities of daily living. Instrumental evaluations were conducted using optoelectronic, force plates and electromyographic (EMG) systems in a motion analysis laboratory. Statistical analysis involved a non-parametric test to compare pathological and normal data, clustering methods to identify groups based on clinical evaluations, and a combination of non-parametric analysis and linear models to assess dependencies on clinical scales.

Results: The results showed that PD patients had significant gait kinematic differences compared to normal values, with increased temporal and shortened spatial parameters. In addition, PD patients were grouped into four clusters based on clinical scales. While some gait features were influenced by clinical scales reflecting impairment, gait and balance, and independence, others were more affected by the perceived fear of falling (FoF).

Discussion: In conclusion, the study identified specific biomechanical gait dysfunctions in kinematic, kinetic, and electrophysiological parameters in PD patients, undetectable by standard clinical scales. Additionally, higher FoF was associated with dysfunctional biomechanical patterns, independent of impairment severity, gait and balance dysfunction, or overall independence.

1 Introduction

Parkinson’s disease (PD) is one of the most common neurodegenerative disorders that affects the substantia nigra of the midbrain. PD is associated with motor, cognitive, and autonomic alterations (Hallett, 2014; Jellinger, 2014; Giustiniani et al., 2022). From the early stages of PD, patients may present dysfunctional gait patterns, mainly characterized by a shortened stride length, increased stride variability, reduced walking speed and festinating gait (Bonanno et al., 2023). The most detectable gait modifications in PD encompass various alterations, including decreased walking speed and step length, asymmetrical gait patterns, reduced arm swing, compromised dissociation between arm and trunk movements during ambulation, and postural instability (Zanardi et al., 2021). Pistacchi et al. (Pistacchi et al., 2017) found that reduced ankle flexion during the swing phase as well as in the stance phase, and mild knee flexion reduction can be detected even in early-stage PD patients. Some gait alterations in these patients are specific and get worse throughout the disease (Zanardi et al., 2021). Gait impairment in PD evolves as the disease progresses, showing distinct patterns at early, mild-to-moderate, and advanced stages. In the early phase, specific changes include reduced arm swing, less fluid movement, and increased asymmetry between limbs. Impaired muscle contraction, rigidity, and postural instability also contribute to reduced limb propulsion, negatively impacting gait parameters like speed and step length. Notably, reduced step length is a characteristic feature of Parkinsonian gait. As PD reaches the mild-to-moderate stage, symptoms become more symmetric, reducing earlier asymmetries. Gait difficulties intensify, with patients often exhibiting shuffling steps, prolonged double-limb support, and increased cadence. At this stage, motor automaticity further reduces, leading to fragmented movements, such as turning with a “block-like” motion (turning in block) and challenges with initiating steps. This pattern reflects both a loss of automaticity and impaired forward propulsion, hallmark characteristics of Parkinsonian gait. Together, these results further confirm that PD gait disturbances are characterized by slower, shorter, and less coordinated movements, with changes in both timing and spatial dimensions of movement that contribute to decreased mobility and stability. This reduced stability is consequently perceived in patients with PD as an enhanced fear of falls (FoF). According to Bryant et al., people with PD with higher levels of FoF tend to exhibit specific biomechanical gait features, such as slow gait speed, short stride length and balance deficits may predispose to falls (Bryant et al., 2014). In particular, gait speed is an important factor that should be considered in the rehabilitation training of Parkinsonian patients. Indeed, some authors (Paker et al., 2015) found that gait speed is strongly correlated with not only demographic factors (e.g., age, gender, height) but also with FoF, balance and mobility. In this sense, early identification and quantification of biomechanical gait parameters are key factors in establishing an effective and customized rehabilitation therapy (Biase et al., 2020; Zanardi et al., 2021). The quantification of disease progression in PD is mainly measured through the Hoehn and Yahr scale (H&Y) (Zhao et al., 2010) which is one of the most widely applied indices of disease severity. However, this clinical assessment tool is focused on gross motor functions, analysing roughly the patients’ quality of movements.

On the other hand, clinical scales or tests to assess gait functions are also often used in clinical practice to evaluate patients’ motor abilities, and “how” they walk, considering also the activities of daily living, and the level of fatigue (Biase et al., 2020). Specifically, this kind of assessment is commonly administered by physiotherapists to investigate the level of functional status of patients, before and after the rehabilitation path (Bonanno et al., 2023). Among the most used rating scales, in the PD population, are the Tinetti scale (TS) for gait and balance (Tinetti, 1986), and Berg Balance Scale (BBS) to assess static and dynamic balance (Berg et al., 1995), and the Falls Efficacy Scale-International (FES-I) to assess the perceived fear of falling (Caronni et al., 2022). However, clinical scales do not always reflect the “objective” stage of the disease, and they cannot detect specific movement changes that are not visible to the clinician’s eye, as some instrumental assessments do. Current approaches to evaluating objectively and quantitatively gait impairments include marker-based motion capture (MoCap) systems (Romero-Flores et al., 2024). In particular, optoelectronic tools are considered as the “gold standard”. These can use either passive systems consisting in a set of multiple cameras to track 3D trajectories (i.e., Vicon, Oxford UK), or active systems that can identify automatic movements through infrared light-emitting diode (LED) (i.e., Optotrack, NDI, Waterloo, ON, Canada), placed on different body landmarks. In this way, these systems allow a comprehensive recording of kinematics and kinetics. Additionally, they can be equipped with force platforms to detect ground reaction forces as well as electromyographic (EMG) sensors to record muscle activity. In this way, the conjunct use of both instrumented gait analysis (GA) and clinical scales became fundamental in the rehabilitation context, especially in PD patients, who manifest gait alterations at the early stages of the disease, but also to monitor the progression of the disease. Indeed, clinical scales, while valuable for assessing general symptoms and overall disease severity in PD, usually focus on broad aspects of motor symptoms (e.g., tremor, rigidity, bradykinesia) rather than specific gait parameters like spatial-temporal parameters and joint angles. These specific parameters are essential to understand the true nature of gait impairment in PD, but they are not well represented in most standard scales. In addition, gait alterations, such as slight reductions in step length or small increases in gait variability like in early stages of the disease, may go unnoticed in a clinical setting, where scoring often lacks the granularity needed to capture these early signs of impairment. For these reasons clinical scales cannot be the only assessment tool for PD motor evaluation. Another aspect that should not be neglected is that instrumental GA parameters could be more effective in detecting changes in biomechanical patterns of gait in PD patients after a rehabilitation program. In this context, Peppe and colleagues highlighted the efficacy of a rehabilitation intervention, which aimed at improving temporal parameters of gait, measured with instrumental GA (Peppe et al., 2007). These results could be valuable in creating a more objective assessment of motor rehabilitation programs and in demonstrating their impact on neurodegenerative conditions like PD. The findings on kinematic data of the lower limbs align with those documented in existing literature. In this sense, a complete assessment of gait, which comprises both instrumental and clinical measures, can be useful to personalize rehabilitation treatment according to patients’ needs. In addition, this combined approach can be also useful to monitor patients over time, during the progression of neurodegenerative disease.

In this observational study, our primary aim is to identify and quantify kinematic, kinetic, and EMG gait parameters obtained by a GA laboratory, which integrates optoelectronic MoCap system, force platforms and EMG sensors, in a group of PD patients, who are naive to gait training, compared to normal data.

Then, our secondary aim is to investigate biomechanical differences among PD patients by clustering them based on clinical rating scales (BBS, TS, Barthel index, H&Y and FES-I) and analysing variations in instrumental biomechanical data across these clusters.

2 Materials and methods

2.1 Study design and population

Sixty-nine PD patients, with a mean age of 66 ± 9 (see Table 1), were consecutively recruited from the Outpatient Movement Disorders Clinic of the IRCCS Centro Neurolesi “Bonino-Pulejo” (Messina, Italy) between June 2018 and March 2024.

Table 1
www.frontiersin.org

Table 1. Socio-demographic and clinical data of the study sample.

Patients were included if they: (1) had a diagnosis of PD according to the Movement Disorder Society Clinical Diagnostic Criteria for Parkinson’s Disease; (2) were naive to gait training; (3) were in “on” phase at the time of the evaluation and (4) were able to walk independently. Patients were excluded if they had: (1) cognitive, visual, or auditory deficits that could impair the comprehension and/or execution of the proposed motor tasks for the evaluation; (2) presence of comorbidities that affected upright posture and walking (e.g., hypotension); (3) refused consent or were unable to provide informed consent.

All experiments were conducted according to the ethical policies and procedures approved by the local ethics committee (IRCCS-ME-23/2022). All participants gave their written informed consent.

2.2 Procedures

All patients were assessed during their “on” phase while taking their own treatment dose, as per clinical indications. Evaluations took place in the morning, approximately 2 hours after the medication intake. Both the instrumented GA and clinical motor assessments were conducted on the same day as the routine neurological consultation by a multiprofessional team composed of a neurologist, a physiotherapist, a bioengineer and a motor science technician. This latter prepared PD patients for instrumented GA, the physiotherapist administered the clinical scales/tests, while the neurologist carried out the neurological examination. The GA and the subsequent analyses were performed by the bioengineer.

2.3 Clinical assessment

A skilled physiotherapist evaluated motor performances of PD patients using the following clinical scales/tests: i) BBS assesses static and dynamic balance using 14 tasks. The score ranges from 0 to 56, with scores below 40 indicating a moderate to high risk of falling; ii) TS is a 16-item measure (7 items for gait and nine for balance), where a total score of 18 or lower indicates a high risk of falls, and a score between 19 and 24 signifies a moderate risk of falls; iii) FES-I is a 16-items measure of perceived FoF, it ranges from a minimum of 16 (no fear/concerns of falling) to a maximum of 64 (strong concern about falling); iv) Barthel index (BI) to measure an individual’s ability to perform basic activities of daily living (ADLs) independently. Each activity is scored based on the level of assistance required, with a total possible score ranging from 0 to 100. Higher scores indicate greater independence, while lower scores suggest higher levels of dependency. Moreover, the neurologist administered the H&Y scale to classify PD patients according to their disability level.

2.4 Instrumental gait analysis

A skilled physiotherapist and a biomedical engineer assessed patients’ gait strides using the BTS Gaitlab (BTS Bioengineering, Milan, Italy) (Materia, 2024). During the instrumental GA assessment, patients were instructed to walk at a self-selected/comfortable speed to obtain the most accurate and representative data on their walking capabilities. This advanced GA system has fully integrated tools that enable objective, quantitative assessments for clinical use. Through instrumental analysis, clinicians and therapists gain an objective and quantitative view of changes in posture and gait, load imbalances, and muscle deficits often undetectable with standard clinical scales or tests.

The BTS Gaitlab system comprises.

- 8 infrared cameras (BTS SMART-DX)

- 4 force plates (BTS P-6000)

- 8 wireless EMG probes (BTS FREEEMG 1000)

The BTS SMART-Clinic software (Figure 1) includes libraries with scientifically validated protocols (Kadaba et al., 1989; Davis et al., 1991). In our study, we applied the “DAVIS Heel: multifactorial gait analysis” protocol, which provides quantitative, objective data on kinematics, kinetics, and associated muscle activity to assess gait functionality. This protocol is based on the “Newington marker set” from the Davis protocol (Davis et al., 1991) and requires measuring the participant’s anthropometric parameters, including weight, height, tibia length, femoral condyle distance, knee and ankle diameters, iliac crest distance, and pelvis thickness. The optoelectronic system captures the kinematic data by tracking marker positions on the patient’s body, calculating hip, knee, pelvis, trunk, and ankle joint angles in flexion-extension, abduction-adduction, and external-internal rotation. Events in the gait stride, such as initial ground contact and toe-off, are identified automatically through marker and force plate data, with manual verification by therapists.

Figure 1
www.frontiersin.org

Figure 1. BTS SMART-Clinic software. The figure shows the data processing procedure for one of the participants. A 3D human body model, reconstructed using the optoelectronic system, is displayed along with the detected events and the synchronized video recordings from the gait trial.

The initial position (standing phase) is recorded by asking the patient to stand neutrally for 5 s. The patient then walks at a self-selected pace. After an average of six acquisitions (mean ± SD = 2.8 ± 0.9 strides per acquisition), baseline evaluation is complete. During analysis, EMG activity is also recorded through the eight wireless EMG probes (four on each limb), connected to the Smart Analyzer system (Version 1.10.469.0; BTS, Milan, Italy). In this study, EMG signals from the gastrocnemius lateralis (GL), tibialis anterior (TA), rectus femoris (RF), and semitendinosus (ST) muscles were evaluated following European guidelines for surface EMG (SENIAM) (Hermens et al., 2000), with skin prepared (cleaned and dried) before positioning bipolar surface electrodes aligned with muscle fiber orientation.

Normative data from healthy adults (40 participants: 28 males, 12 females, ages 18–40) (Kadaba et al., 1990) were used to compare the patients’ kinematic and kinetic, helping to identify gait dysfunction. Furthermore, data from BTS Gaitlab processed during offline analysis included:

2.4.1 Kinematics

Spatio-temporal parameters, hip, knee, pelvis, trunk, and ankle joint angles in flexion-extension, abduction-adduction, and external-internal rotation measured during walking trials.

Spatio-temporal parameters.

• Spatio-temporal parameters were recorded during the walking trial. Normative data were obtained from BTS Gaitlab. Temporal aspects, including gait stride and stance duration (measured in seconds) as well as stance, swing, single and double support phases normalized to gait stride (%), were included alongside spatial parameters such as stride and step length, step width (m), and spatio-temporal parameters such as gait speed in (m/s). Spatio-temporal parameters were compared against the corresponding average gait feature for health people to identify gait impairments. Moreover, the gait profile score (Speciali et al., 2014), calculated as the Euclidean distance between the patient’s kinematic features and the corresponding normative features over the entire gait cycle, and the gait deviation index (Schwartz and Rozumalski, 2008), which indicates whether the subject’s gait features are statistically indistinguishable from those of the control group, were also evaluated. Additionally, the RCIRC shape symmetry index (Equation 1) (Bonanno et al., 2024) was used to compare PD patients and healthy subjects in terms of joint rotation angles, reported with respect to the percentage of gait stride (%), quantifying gait dysfunction by estimating differences between pathological and normal gait patterns.

RCIRC=Cxyn=1101xn2n=1101yn2(1)

In which x is the waveform related to the PD subjects, while y corresponds to the average waveform of the healthy population. Cxy is the circular cross-correlation function at lag 0. RCIRC ranges from −1 to 1 (i.e., identical amplitude profiles shape).

2.4.2 Kinetics

Joint moments and powers at the hip, knee, and ankle joints:

• Kinetic results, including joint moments and powers normalized to the subject’s weight (Newton*meter/kg and Watt/kg) as well as ground reaction forces expressed as a percentage of body weight, were averaged across all gait strides for each participant. Additionally, the RCIRC shape symmetry index (Equation 1) was calculated to compare joint moments, powers, and ground reaction forces between pathological and healthy groups.

2.4.3 EMG

Muscle activation and co-activation patterns:

• Muscle activation signals were recorded with surface electrodes, capturing raw EMG signals (mV), which were filtered with a 20–450 Hz bandpass filter and rectified. The amplitude was normalized to the maximum of any channel and the time to the gait stride duration (% of gait stride). Signal amplitude, which is proportional to the force generated by the muscle, was analysed by calculating the co-contraction of agonist and antagonist muscles during gait enhances support, balance, propulsion, and movement efficiency. Estimating muscle co-contraction provides valuable insights into how a disorder may affect muscle coordination strategies. For this purpose, we used a co-contraction or co-activation index (CoAct) method based on (Equation 2), (Bonanno et al., 2024), where the normalized EMG of the antagonist muscle (norm EMG antago (t)) was the lower value and the normalized EMG of the agonist muscle (norm EMG ago (t)) was the higher value.

CoActt=2×normEMGtantagonormEMGtantago+normEMGtago×100(2)

• In which the agonist-antagonist muscle pairs used for this analysis were TA-GL and RF-ST.

2.5 Statistical analysis

Kinematic and kinetic instrumental outcomes were compared to the normal values (Kadaba et al., 1990). The normal distribution of the samples for each parameter—including kinematic, kinetic and EMG—was investigated through Shapiro-Wilk test (MATLAB function swtest). According to the non-normality of most parameters, we chose a non-parametric analysis. The Wilcoxon rank sum test (MATLAB function ranksum) was used to compare overall pathological subjects’ parameters with healthy normative values. The same test was also applied to assess left-right sides differences. Data were analyzed both as averages between body sides and separately for each side, considering all gait strides of each participant.

We performed a clustering to partition the dataset into groups. To determine the optimal number of clusters, we performed clustering using the k-means algorithm (MATLAB function kmeans). We set a maximum of n = 20 clusters, iterating over possible cluster numbers from 2 to 20. We executed the k-means function on clinical scale evaluations, with 50 replicates to ensure stable convergence of the cluster centers. This generated the within-cluster sum of squared distances. To identify the ideal number of clusters, we evaluated the cumulative percentage of the explained variance. We selected the minimum number of clusters at which the explained variance reached or exceeded 90%, ensuring a balance between data representation accuracy and model simplicity.

Instrumental and clinical scales outcomes were analyzed with a nonparametric one-way ANOVA (Kruskal Wallis, MATLAB function kruskalwallis) and a LM (MATLAB function fitlm) because of the non-normal distribution of the kinematic data. The dependency of the response variable from the experimental factor clusters was tested with the Kruskal Wallis method, while the dependency from the experimental factor FES-I was tested with the LM in a post hoc analysis. The robust fitting type used for the LM was the “bisquare” weight function with the default tuning constant. All the analyses were implemented in Matlab (MATLAB (R2022a), Natick, Massachusetts: The MathWorks Inc.; 2022).

To assess whether our sample size was adequate for the analyses performed, we conducted two post hoc power analyses using G*Power (version 3.1.9.7). For the comparison between Parkinson’s patients and normative values, we used a t-test (Means: difference from constant, one sample case), with α = 0.05 and a total sample size of 69. For the comparison across clusters, we used an F test (ANOVA: fixed effects, omnibus, one-way), with α = 0.05, a total sample size of 69, and four groups. In both cases, the effect size was estimated from the data collected in our study, based on the variables compared. Specifically, for each analysis, we computed the effect size for each parameter separately using group means and standard deviations and then used the average value as input for the power calculation.

3 Results

Comparisons between participants’ instrumental and normative data revealed statistically significant differences in temporal and spatial kinematic parameters. Specifically, stride duration, stance duration, and swing duration average data were significantly higher than normative values (see Figure 2) showing a statistically significant difference (p < 0.001). The stance phase and double support phase also showed higher average values, whereas the swing phase and single support phase showed lower values compared with normative data showing a statistical difference (p < 0.001).

Figure 2
www.frontiersin.org

Figure 2. Kinematic gait parameters. The figure shows some of the temporal and spatial kinematic parameters, reported as mean ± standard deviation (A) Compares temporal parameters (gait stride, stance and swing duration) between participants with PD (red) and healthy subjects (grey) (B) Shows the same comparison for spatial parameters (stride length, step length and step width). Statistical difference significance, evaluated using Wilcoxon rank sum test, is reported as *** < 0.001, ** < 0.01 and * < 0.05.

Both speed and cadence differed significantly from normative data (p < 0.001), with PD participants exhibiting lower values. Additionally, a reduced stride length and step length and an increased step width were observed, with a statistical difference (p < 0.001). The gait profile score was notably higher than normal, whereas the gait deviation index was lower, showing a statistical difference (p < 0.001).

The analysis of side-to-side symmetry revealed a statistically significant asymmetry in RCIRC between left and right sides for both trunk obliquity (p < 0.001) and trunk rotation (p = 0.016) (see Supplementary Table S1).

The minimum number of clusters explaining at least 90% of the variance was four clusters, capturing 92% of the variance. Clusters were created based on clinical scales, presenting a gradient from the highest (Cluster 1) to the lowest motor functional state (Cluster 4) (See Table 2).

Table 2
www.frontiersin.org

Table 2. Socio-demographic and clinical data of the clusters. Statistical difference between clusters was evaluated using nonparametric one-way ANOVA (Kruskal Wallis).

This descending trend was observed across all the assessed scales (BBS, TS, BI, and H&Y), except for the FES-I scale, in which clusters 1 and 3 showed lower perceived FoF. Clusters two and 4 showed higher levels of FoF. Effectively, clusters categorized patients from least to most severe motor impairment, although they displayed distinct perceptions of FoF (see Figure 3).

Figure 3
www.frontiersin.org

Figure 3. Distribution of clinical scale scores across identified clusters. The figure shows clustering results based on clinical scales, with a colour gradient representing severity levels on each scale. The gradient ranges from less severe (green) to more severe (red). This allows for an accurate visualization of severity patterns specific to each cluster and clinical scale.

Statistical analysis indicated a significant difference between clusters and a significant effect of FES-I on the clinical scales: BBS (p < 0.001; R2 = 0.29, p < 0.001, see Figure 4A), TS (p < 0.001; R2 = 0.28, p < 0.001), BI (p < 0.001; R2 = 0.16, p = 0.001), H&Y (p < 0.001; R2 = 0.31, p < 0.001) and FES-I (p < 0.001; R2 = 1, p < 0.001, see Figure 4B).

Figure 4
www.frontiersin.org

Figure 4. Motor clinical scales and instrumental biomechanical gait parameters. The left column shows instrumental parameters (C,E,G) with trends similar to BBS (A) while the right column includes those (D,F,H,I) following trends similar to FES-I (B). Dependency of instrumental parameters from both clinical scales (A,B) was statistically confirmed. Each panel shows the box plots for the four identified clusters (represented by different colors), with red crosses indicating outlier values.

The ANOVA revealed significant differences among clusters across several instrumental parameters. In terms of temporal kinematics, we observed significant differences for stride duration (p = 0.049, see Figure 4C), stance duration (p = 0.013), stance phase percentage (p = 0.023, see Figure 4E), and single support phase percentage (p = 0.025). The spatial kinematic variable step width (p = 0.024) also showed significant cluster-based variation. Regarding results obtained from circular cross-correlation, the spatial kinematic hip abduction (p = 0.024, see Figure 4G) and the kinetics such as hip moment (p = 0.001), knee moment (p = 0.003), ankle moment (p = 0.004), hip power (p = 0.001), ankle power (p = 0.008), and medio-lateral force (p = 0.002) were significantly different across clusters. In EMG measurements, significant cluster differences were noted for the CoAct of the TA-GL muscles (p = 0.003).

For parameters influenced by both cluster and FES-I, we observed significant effects for temporal kinematics parameters including swing phase (p = 0.025; R2 = 0.06, p = 0.049, see Figure 4F) and double support phase (p = 0.003; R2 = 0.07, p = 0.046). Similarly, spatial kinematic parameters, such as stride length (p < 0.001; R2 = 0.16, p = 0.001, see Figure 4D) and step length (p < 0.001; R2 = 0.15, p = 0.001) showed significant effects. The spatiotemporal kinematic variable, velocity, was significantly influenced by both cluster and FES-I (p < 0.001; R2 = 0.13, p = 0.002, see Figure 4I). Joint angles as well as kinetic measurements also demonstrated dual significance for the circular cross-correlation of pelvic rotation (p = 0.032; R2 = 0.11, p = 0.006, see Figure 4H), knee flexion-extension (p = 0.007; R2 = 0.09, p = 0.015), anterior-posterior force (p < 0.001; R2 = 0.13, p = 0.004), and vertical force (p < 0.001; R2 = 0.14, p = 0.002).

In line with our findings, the dependency of some instrumental biomechanical gait parameters on motor clinical scales was confirmed, as illustrated in Figure 4 (left column), in which the boxplots of the four clusters follow the BBS’s trend. On the other hand, the parameters in the right column followed the FES-I’s trend among clusters.

Lastly, parameters affected solely by the FES-I score included EMG CoAct of the RF-ST muscles (R2 = 0.10, p = 0.01).

To sum up, we provided a detailed table reporting mean and SD values for each cluster across all the evaluated parameters and the relative statistical comparison results (see Table 3). We also reported the mean values for each instrumental parameter (columns) and cluster (rows) in the matrix shown in Supplementary Figure S1. The figure shows the instrumental results for each cluster and parameter, with a color gradient representing the relative performance on each scale. The gradient ranges from better performance (green) to worse performance (red), allowing for a clear visualization of how each cluster performs across the different instrumental parameters.

Table 3
www.frontiersin.org

Table 3. Instrumental data across clusters. The table reports instrumental variables (mean ± SD) for the four identified clusters. Statistical differences between clusters were assessed using nonparametric one-way ANOVA (Kruskal Wallis test). The influence of the FES-I score was evaluated through a LM.

The results from the post hoc power analysis indicated that our sample size was sufficient to support the main comparisons. We conducted two post hoc power analyses to assess this formally. For the comparison between PD patients and normative values, the power analysis using the observed average effect size (Cohen’s d = 0.5) indicated a high statistical power (1–β = 0.99). For the comparison across the four Parkinson clusters, the observed average effect size (Cohen’s f = 0.4) yielded a power of 0.78, indicating sufficient sensitivity for detecting medium-to-large differences among groups.

4 Discussion

In this study, we highlight the importance of complementing clinical assessments with instrumental GA. While comparing gait biomechanics between PD patients and healthy individuals is not novel in itself, our study enhances and strengthens existing literature by providing a comprehensive biomechanical assessment, including kinematics, kinetics, and EMG analysis, as detailed in (see Supplementary Table S1). Notably, we also report joint angular measurements, which are underrepresented in the literature, as noted by Zanardi et al. in their systematic review with meta-analysis.

Additionally, the clustering analysis, which classified patients based on their clinical scores, revealed that certain instrumental GA parameters, such as stride length, swing phase, pelvic rotation, and gait speed, do not follow the same pattern as clinical scales, highlighting the importance of objective gait assessment. Similarly, Russo et al. conducted a cluster analysis on spatiotemporal parameters extracted from an optoelectronic device in a sample of people with PD (Russo et al., 2023). However, our study goes further by incorporating kinematic angular, kinetic, and EMG parameters alongside clinical scales, adding a novel dimension to our findings. Furthermore, we identified specific dysfunctional gait patterns in individuals with PD who experience higher levels FoF, providing valuable insights into the prevention of potential falls in this patient population.

Consequently, these findings can provide new perspectives to better guide rehabilitation strategies for people with PD, according to their specific needs.

4.1 Instrumental and biomechanical gait parameters

In this observational study, we first aimed to characterize the biomechanical pattern of gait in a cohort of PD patients compared to normal values. Our results suggest that PD patients had significant differences in the kinematic parameters of gait. Indeed, PD patients showed higher temporal (stride, stance and swing duration) parameters than normative values. Conversely, spatial parameters such as stride length and step length were notably reduced in PD patients, compared to normative values. These findings align with previous research (Pistacchi et al., 2017; Zanardi et al., 2021; Russo et al., 2025) indicating that PD-related gait alterations commonly involve prolonged time spent in each phase of the gait cycle, accompanied by shorter step and stride lengths. According to Peppe et al., patients with PD exhibit an increased flexion in the hip, knee, and ankle, both in standing and walking positions, compared to control subjects. Moreover, these authors (Peppe et al., 2007) found that following the motor rehabilitation program, the joint angle values in PD patients, except for the knees, were closer to those observed in the control group. Regarding kinematic results, we also found that people with PD exhibited higher gait profile score compared to normative values. Gait profile score is a single score that represents the overall deviation of gait kinematics from normative data, providing a quantitative assessment of the gait profile. This score is derived from the analysis of nine kinematic variables, evaluated on both sides of the body (Speciali et al., 2014). This result is in line with findings of other studies showing an altered gait profile, possibly explained by the variation found in some of the gait profile score domains (including knee flexion/extension and pelvic obliquity). These angular changes in gait are commonly found in PD and are justified by the reduction in range of motion of lower limbs and difficulty in regulating and coordinating movements during gait (Speciali et al., 2014; Pistacchi et al., 2017; Callais Franco do Nascimento et al., 2021). In addition, we found that the gait deviation index was lower in people with PD compared to normative values. This index is a multivariate measure of overall gait pathology based on 15 gait features derived from 3D kinematic data, offering a clear, comprehensive, and clinically relevant measure of overall gait function (Schwartz and Rozumalski, 2008). In the field of PD, only two studies have explored the role of the Gait deviation index in Parkinsonian gait. Specifically, Galli et al., suggested gait deviation index was more useful for assessing the impact of levodopa treatment on gait rather than for identifying the severity of gait pathology in PD patients (Galli et al., 2012). Other findings were reported by Speciali et al., who stated that this index effectively captures gait alterations during dual-task exercises. Additionally, we found a significant difference between the left and right side in the trunk obliquity and trunk rotation RCIRC, suggesting a potential musculoskeletal asymmetry between the two sides. According to Cano-de-la-Cuerda et al., trunk alterations may result from axial muscle stiffness. Patients may exhibit a reduced joint range of motion due to excessively high muscle tone, which restricts full articulation and leads to asymmetries between the two sides.

In this vein, it is noteworthy that the use of instrumental GA can provide objective outcome measures for a rehabilitation program. It can also contribute to give additional information on specific dysfunctional gait patterns that could aid in understanding the complex pathophysiology of PD.

4.2 Biomechanical gait and clinical differences among clusters

Secondly, we aimed to investigate biomechanical gait differences among PD patients by clustering them into four groups based on clinical rating scales (H&Y, BBS, TS, BI, and FES-I) and analysing variations in instrumental biomechanical data across these clusters. We found statistically significant differences among clusters in kinetic gait parameters, such as hip, knee, and ankle moments, as well as hip and knee power. These differences could be related to disease progression, as patients with more severe impairment (such as those in Cluster 4) may have developed compensatory mechanisms, such as a hip-driven strategy to compensate for reduced ankle and knee function. This hypothesis is further supported by the observed CoAct of the RF-ST muscles, which may indicate increased recruitment of proximal muscles rather than distal ones, such as TA-GL (Islam et al., 2020). Notably, patients with higher levels of disability may need to reduce their knee power and moment to minimize lower limb instability, which could increase the risk of falls. Additionally, the typical muscle rigidity and bradykinesia in individuals with PD may lead to a reduced range of motion, consequently diminishing lower limb power. This reduction appears to be closely related to the severity of impairment within each cluster. Furthermore, we observed a statistical difference among clusters regarding medio-lateral force. In particular, Cluster 4 showed a reduction in medio-lateral force, which may indicate a more cautious gait pattern aimed at minimizing balance perturbations. This reduced medio-lateral force could reflect compensatory strategy to minimize balance perturbations, given the increased postural instability in patients with severe PD. This aligns with previous research (Kawami et al., 2022) indicating that PD patients with medio-lateral balance impairments demonstrate reduced variability in medio-lateral center of mass movements, potentially as a response to increased postural sway during gait. While this strategy may initially help stabilize gait, it may also lead to a narrower movement area, reducing the ability to make necessary postural adjustments and ultimately increasing fall risk. In this vein, rehabilitation interventions should aim to enhance medio-lateral stability and optimize stance width to improve balance control and prevent falls in patients with advanced PD (Kawami et al., 2022).

Interestingly, according to our cluster analysis, we found that PD patients can be grouped into two main categories, those with high levels of FoF and those with low levels of FoF independent of their motor impairment. FoF is a common symptom in patients with PD, and it is considered one of the most stressful physical symptoms (Jonasson et al., 2018). In addition, FoF was found to be a predictor of future falls and near falls already in mild PD and is negatively associated with participation and health related as well as overall quality of life. In our study, we found that PD patients with high levels of FoF had specific dysfunctional gait patterns that were different from those with low levels of FoF. In line with our results, some authors suggested that gait speed and stride lengths were poorer in people with a high level of FoF (Bryant et al., 2014), while Uhlig et al., (Uhlig and Prell, 2023a), found that only gait speed was associated with FoF.

This fear, along with decreased confidence in one’s ability to perform daily activities, can contribute to or worsen cautious walking. In PD, reduced gait speed—linked to a more cautious gait and FoF—can be strongly associated with previous fall incidents (Mak and Pang, 2009; Kataoka et al., 2011; Liu et al., 2022; Uhlig and Prell, 2023a). However, we did not analyse this aspect, but in future studies this aspect should be further investigated, correlating previous fall incidents with the FoF and biomechanical parameters.

A cautious gait, which is also commonly seen in healthy older adults (Uhlig and Prell, 2023b), involves slower gait speed, shorter step length, and lower toe clearance, allowing the person to stay “as close as possible” to the floor for further stability. Compared to a previous study (Bryant et al., 2014) in literature, our findings indicate that in individuals with high levels of FoF, not only gait speed is reduced, but we also observed an increase in cadence and a decrease in stride length. These results suggest that PD patients who experience higher FoF tend to adopt shorter, more controlled steps, characterized by reduced speed but an increased step frequency. This altered gait pattern could potentially create a vicious cycle, leading to greater energy expenditure at the muscular level. A recent meta-analysis suggested that PD patients compared to healthy controls manifest a reduced muscle activity of gastrocnemius medial and higher activity of TA, accompanied by a higher co-contraction of these ankle muscles during gait. These alterations can influence adequate transfer weight in preparation for stepping and it can reflect in a higher metabolic cost of walking (Zanardi et al., 2021). However, we found that patients with higher FoF manifested an increased CoAct of RF-ST that was significantly correlated with FES-I score. From a biomechanical point of view, EMG co-contraction patterns could have a compensatory role in increased proximal muscle activity (e.g., RF-ST) in response to reduced distal muscle function. Enhanced quadriceps activation during the stance phase promotes greater knee extension, improving joint stability in single stance and potentially compensating for reduced ankle stability (Islam et al., 2020). Similarly, increased hamstring activity during the swing phase enhances hip extension and knee flexion, partially substituting for the foot placement and initial loading functions typically performed by distal ankle muscles. However, this increased muscle activity also raises metabolic demand, which may contribute to reduced walking speed and mobility. Moreover, it is worth noting that muscle rigidity is a cardinal symptom of PD, as defined by the International Parkinson and Movement Disorder Society. Rigidity in PD is characterized by a velocity-independent increase in muscle tone, assessed through passive muscle stretch at specific joints. Importantly, rigidity and postural instability, two hallmark motor features of PD, are likely to shape EMG patterns during gait. Baradaran et al. demonstrated that rigidity is linked to altered cortical/subcortical connectivity, including changes in the supplementary motor area and putamen, along with increased motor cortex excitability (Baradaran et al., 2013). A potential consequence of this disrupted neural control is heightened agonist/antagonist co-contraction and reduced selective muscle recruitment, which may contribute to gait dynamic instability. This instability, defined as difficulty transitioning between gait phases, can manifest as prolonged double-support time and increased stride time variability. Other elements that could have an influence of biomechanical gait parameters are cognitive and psychological factors, like anxiety (Zhang et al., 2024). While assessing cognitive and psychological status was beyond the scope of our study, future research integrating a comprehensive biomechanical assessment with motor and cognitive/psychological evaluations could offer deeper insights into gait characteristics in individuals with PD.

Furthermore, we found that pelvic rotation and knee flexion extension were reduced when compared to normative values, especially in PD patients with high levels of FoF. This finding could be explained by the fact that patients with high FoF tend to adopt a more cautious gait to avoid falling, limiting the range of motion in their knees and pelvis. What is more, FoF could increase muscle stiffness, as we found within the CoAct of lower limb muscles, further reducing the motion flexibility of the knee and pelvis. In addition, the reduced knee flexion-extension and pelvic rotation may indicate a diminished ability to generate joint power in the hip and knee, as we found in Cluster 4. This limitation in movement could be associated with decreased gait propulsion and coordination, which are further exacerbated by muscle rigidity and difficulty executing broad and rapid movements (Paul et al., 2012). Another aspect we investigated involves ground reaction forces, specifically the anterior-posterior and vertical forces exerted by the ground on the body during movement (Chockalingam et al., 2016; Alam et al., 2017). We observed that PD patients with high levels of FoF tend to exhibit ground reaction forces in the anterior-posterior and vertical directions that deviate more significantly from normal values during gait, which is in line with the current literature (Russo et al., 2025). This finding could suggest that they may exert dysfunctional forces in these directions, possibly because high FoF leads them to move more cautiously to avoid falling. In this context, Oh et al. demonstrated that people with PD exhibited a lower peak in both the vertical and anterior-posterior ground reaction forces during propulsion (Oh et al., 2020), which aligns with our observation of altered ground reaction forces patterns in PD patients with high levels of FoF.

4.3 Strengths, limitations, and future perspectives

In this study, we highlighted the fundamental role of instrumental GA and clinical motor assessment, suggesting that the combined use of both these tools can be useful to assess quantitatively and objectively motor and biomechanical features of gait in PD patients. However, one of our strengths is that we analyzed several biomechanical parameters, including not only kinematics, but also kinetics and electrophysiological data, which are often neglected in clinical contexts.

Another strength of our study is the use of an optoelectronic MoCap in addition to force platforms and EMG probes to analyze gait biomechanics. In this sense, only a few authors have used this comprehensive GA laboratory (Pistacchi et al., 2017; Russo et al., 2023), which is considered the gold standard for the analysis of biomechanics of human movement (Ceseracciu et al., 2014). Despite this important advantage of this type of MoCap, it does not allow real-world movements detection, as wearable devices can do. For instance, wearable devices can be used outside the laboratory setting giving a relatively realistic representation of daily life relevant mobility aspects (Scott et al., 2022), providing kinematic gait data such as step length, stride width (Rossanigo et al., 2023). However, the amount of data that can be gathered in a laboratory is generally greater than that obtained from wearable sensors, unless multiple systems are used synchronously (e.g., EMG, a suit with multiple IMUs, insoles for plantar pressure recording). Nevertheless, we chose to use a controlled recording environment to ensure both the standardization of the acquisition protocol and safety reasons. In comparison with other previous studies (Pistacchi et al., 2017; Schniepp et al., 2021; Zanardi et al., 2021), we used methods to classify PD patients according to clinical assessment in order to evaluate the impact of the level of impairment on the instrumental gait parameters, by using clustering that aimed to discover naturally occurring patterns or groupings within the data. This type of analysis was also performed by Russo et al., in which they extracted two clusters based on spatial-temporal gait parameters from a sample of people with PD. However, they did not include kinetic factors and EMG signals. In this sense, our work expanded and enriched the literature in this field, suggesting that a comprehensive instrumental and clinical assessment is fundamental to address personalized patients’ needs for rehabilitation.

In addition, to analyze gait patterns across gait stride, we used methods, like the circular cross-correlation. On the other hand, to describe the energetic costs related to dysfunctional of gait, we performed methods to quantify the co-contraction between antagonist and agonist muscles.

Our study has some limitations that need to be acknowledged. For instance, a generalization of our results is impossible due to the small number of PD patients involved in our sample. This is why, in some biomechanical gait parameters we observed only a trend but not a statistical significance. However, the sample size aligns with similar observational studies in the field (Zanardi et al., 2021). Moreover, post hoc power analyses confirmed that the collected sample size was sufficient to detect medium to large effects in both the comparison with normative values and the cluster-based analysis, supporting the reliability of the reported findings. Nevertheless, in the cluster-based comparison, smaller differences between groups may have gone undetected due to slightly reduced statistical power, and future studies with larger samples may help clarify subtler distinctions.

Another limitation is the lack of additional clinical data, such as history of previous falls, cognitive functioning, and/or pharmacological dosages, which could have strengthened our results. Moreover, a limitation of this study is the absence of a formal priori power analysis, as the dataset was derived from real-life clinical evaluations conducted during routine practice. Although we have not collected normative data in our own gait analysis laboratory, which could lead to potential biases related to the recording procedures, the normative datasets used in the literature (Kadaba et al., 1990) use instrumentation with similar accuracy and comparable walking distances. In addition, a potential limitation of this study is the lack of an explicit assessment of multicollinearity among the analyzed gait parameters. Given the high correlation that may exist between spatiotemporal, kinematic, and kinetic variables, future studies should incorporate specific statistical approaches to address this issue and minimize the risk of redundancy or misinterpretation of results.

In future studies, the level of FoF should be addressed during the rehabilitation interventions. In fact, it is important to note that the FoF is not merely a perceptual or psychological factor, but it also has a clear impact on walking performance, including kinematics, kinetics and electrophysiological aspects. For example, cognitive and motor training based on VR technologies that add multisensorial stimulation, and feedback can be useful to reduce FoF, since they can simulate different scenarios in a safe and controlled environment (Bonanno et al., 2024; Impellizzeri et al., 2024).

5 Conclusion

Our results highlight that PD affects a wide range of biomechanical gait parameters, including the kinematic, kinetic, and electrophysiological aspects. These gait abnormalities cannot be detected by clinical scales especially in the early stages, as they primarily assess the gross motor functions. Furthermore, we found that patients with higher levels of FoF exhibited distinct biomechanical alterations compared to those with lower levels of FoF, independent of impairment severity, gait and balance dysfunction, and overall independence. However, the question remains open as to whether PD patients tend to manifest FoF due to the motor alterations caused by the disease, or if it is the FoF following the onset of the disease that leads to these alterations. Further larger sample longitudinal studies are needed to better understand gait abnormalities in PD patients and then tailor personalized rehabilitation plans.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the local Ethics Committee approved the study (IRCCS-ME-23/2022). 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. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

PDP: Conceptualization, Data curation, Formal Analysis, Software, Writing – original draft. MB: Conceptualization, Investigation, Methodology, Writing – original draft. CDM: Conceptualization, Data curation, Methodology, Software, Supervision, Writing – original draft. LP: Data curation, Investigation, Methodology, Writing – original draft. ALF: Data curation, Investigation, Methodology, Writing – original draft. GP: Investigation, Methodology, Validation, Visualization, Writing – original draft. MGM: Investigation, Methodology, Validation, Visualization, Writing – original draft. FI: Investigation, Methodology, Validation, Visualization, Writing – original draft. IC: Data curation, Supervision, Validation, Visualization, Writing – original draft. AQ: Funding acquisition, Resources, Supervision, Writing – original draft. RSC: Conceptualization, Methodology, Supervision, Validation, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by Current Research Funds, 2024, Ministry of Health, Italy.

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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/fbioe.2025.1541240/full#supplementary-material

References

Alam, M. N., Garg, A., Munia, T. T. K., Fazel-Rezai, R., and Tavakolian, K. (2017). Vertical ground reaction force marker for Parkinson’s disease. PLoS ONE 12, e0175951. doi:10.1371/journal.pone.0175951

PubMed Abstract | CrossRef Full Text | Google Scholar

Baradaran, N., Tan, S. N., Liu, A., Ashoori, A., Palmer, S. J., Wang, Z. J., et al. (2013). Parkinson’s disease rigidity: relation to brain connectivity and motor performance. Front. Neurol. 4, 67. doi:10.3389/fneur.2013.00067

PubMed Abstract | CrossRef Full Text | Google Scholar

Berg, K., Wood-Dauphinee, S., and Williams, J. I. (1995). The Balance Scale: reliability assessment with elderly residents and patients with an acute stroke. Scand. J. Rehabil. Med. 27, 27–36. doi:10.2340/1650197719952736

PubMed Abstract | CrossRef Full Text | Google Scholar

Biase, L. di, Santo, A. D., Caminiti, M. L., Liso, A. D., Shah, S. A., Ricci, L., et al. (2020). Gait analysis in Parkinson’s disease: an overview of the most accurate markers for diagnosis and symptoms monitoring. Sensors Basel, Switz. 20, 3529. doi:10.3390/s20123529

CrossRef Full Text | Google Scholar

Bonanno, M., De Nunzio, A. M., Quartarone, A., Militi, A., Petralito, F., and Calabrò, R. S. (2023). Gait analysis in neurorehabilitation: from research to clinical practice. Bioengineering 10, 785. doi:10.3390/bioengineering10070785

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonanno, M., De Pasquale, P., De Marchis, C., Lombardo Facciale, A., Paladina, G., Fonti, B., et al. (2024). Might patients with cerebellar ataxia benefit from the computer assisted rehabilitation ENvironment (CAREN)? A pilot study focusing on gait and balance. Front. Bioeng. Biotechnol. 12, 1385280. doi:10.3389/fbioe.2024.1385280

PubMed Abstract | CrossRef Full Text | Google Scholar

Bryant, M. S., Rintala, D. H., Hou, J.-G., and Protas, E. J. (2014). Influence of fear of falling on gait and balance in Parkinson’s disease. Disabil. Rehabil. 36, 744–748. doi:10.3109/09638288.2013.814722

PubMed Abstract | CrossRef Full Text | Google Scholar

Callais Franco do Nascimento, T., Martins Gervásio, F., Pignolo, A., Augusto Santos Bueno, G., Araújo do Carmo, A., Martins Ribeiro, D., et al. (2021). Assessment of the kinematic adaptations in Parkinson’s disease using the gait profile score: influences of trunk posture, a pilot study. Brain Sci. 11, 1605. doi:10.3390/brainsci11121605

PubMed Abstract | CrossRef Full Text | Google Scholar

Caronni, A., Picardi, M., Redaelli, V., Antoniotti, P., Pintavalle, G., Aristidou, E., et al. (2022). The Falls Efficacy Scale International is a valid measure to assess the concern about falling and its changes induced by treatments. Clin. Rehabil. 36, 558–570. doi:10.1177/02692155211062110

PubMed Abstract | CrossRef Full Text | Google Scholar

Ceseracciu, E., Sawacha, Z., and Cobelli, C. (2014). Comparison of markerless and marker-based motion capture technologies through simultaneous data collection during gait: proof of concept. PLoS ONE 9, e87640. doi:10.1371/journal.pone.0087640

PubMed Abstract | CrossRef Full Text | Google Scholar

Chockalingam, N., Healy, A., and Needham, R. (2016). “Interpreting ground reaction forces in gait,” in Handbook of human motion. Editors B. Müller, S. I. Wolf, G.-P. Brueggemann, Z. Deng, A. McIntosh, and F. Miller (Cham: Springer International Publishing), 1–15. doi:10.1007/978-3-319-30808-1_31-1

CrossRef Full Text | Google Scholar

Davis, R. B., Õunpuu, S., Tyburski, D., and Gage, J. R. (1991). A gait analysis data collection and reduction technique. Hum. Mov. Sci. 10, 575–587. doi:10.1016/0167-9457(91)90046-Z

CrossRef Full Text | Google Scholar

Galli, M., Cimolin, V., De Pandis, M. F., Schwartz, M. H., and Albertini, G. (2012). Use of the Gait Deviation index for the evaluation of patients with Parkinson’s disease. J. Mot. Behav. 44, 161–167. doi:10.1080/00222895.2012.664180

PubMed Abstract | CrossRef Full Text | Google Scholar

Giustiniani, A., Maistrello, L., Danesin, L., Rigon, E., and Burgio, F. (2022). Effects of cognitive rehabilitation in Parkinson disease: a meta-analysis. Neurol. Sci. 43, 2323–2337. doi:10.1007/s10072-021-05772-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hermens,, H. J., Freriks, B., Disselhorst-Klug, C., and Rau, G. (2000). Development of recommendations for SEMG sensors and sensor placement procedures. Journal of Electromyography and Kinesiology 10, 361–374. doi:10.1016/S1050-6411(00)00027-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallett, M. (2014). Tremor: pathophysiology. Park. and Relat. Disord. 20, S118–S122. doi:10.1016/S1353-8020(13)70029-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Impellizzeri, F., Maggio, M. G., De Pasquale, P., Bonanno, M., Bonanno, L., De Luca, R., et al. (2024). Coupling neurologic music therapy with immersive virtual reality to improve executive functions in individuals with Parkinson’s disease: a Quasi-Randomized Clinical Trial. Clin. Park. and Relat. Disord. 11, 100277. doi:10.1016/j.prdoa.2024.100277

PubMed Abstract | CrossRef Full Text | Google Scholar

Islam, A., Alcock, L., Nazarpour, K., Rochester, L., and Pantall, A. (2020). Effect of Parkinson’s disease and two therapeutic interventions on muscle activity during walking: a systematic review. npj Park. Dis. 6, 22–16. doi:10.1038/s41531-020-00119-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Jellinger, K. A. (2014). “Parkinson’s disease,” in Pathobiology of human disease. Editors L. M. McManus, and R. N. Mitchell (San Diego: Academic Press), 2021–2035. doi:10.1016/B978-0-12-386456-7.04605-0

CrossRef Full Text | Google Scholar

Jonasson, S. B., Nilsson, M. H., Lexell, J., and Carlsson, G. (2018). Experiences of fear of falling in persons with Parkinson’s disease – a qualitative study. BMC Geriatr. 18, 44. doi:10.1186/s12877-018-0735-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kadaba, M. P., Ramakrishnan, H. K., and Wootten, M. E. (1990). Measurement of lower extremity kinematics during level walking. J. Orthop. Res. 8, 383–392. doi:10.1002/jor.1100080310

PubMed Abstract | CrossRef Full Text | Google Scholar

Kadaba, M. P., Ramakrishnan, H. K., Wootten, M. E., Gainey, J., Gorton, G., and Cochran, G. V. (1989). Repeatability of kinematic, kinetic, and electromyographic data in normal adult gait. J. Orthop. Res. 7, 849–860. doi:10.1002/jor.1100070611

PubMed Abstract | CrossRef Full Text | Google Scholar

Kataoka, H., Tanaka, N., Eng, M., Saeki, K., Kiriyama, T., Eura, N., et al. (2011). Risk of falling in Parkinson’s disease at the hoehn-yahr stage III. Eur. Neurol. 66, 298–304. doi:10.1159/000331635

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawami, Y., Nikaido, Y., Nose, S., Unekawa, M., Marumoto, K., Kawami, M., et al. (2022). Mediolateral postural control during gait in Parkinson’s disease. Prog. Rehabil. Med. 7, 20220048. doi:10.2490/prm.20220048

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W.-Y., Tung, T.-H., Zhang, C., and Shi, L. (2022). Systematic review for the prevention and management of falls and fear of falling in patients with Parkinson’s disease. Brain Behav. 12, e2690. doi:10.1002/brb3.2690

PubMed Abstract | CrossRef Full Text | Google Scholar

Mak, M. K. Y., and Pang, M. Y. C. (2009). Fear of falling is independently associated with recurrent falls in patients with Parkinson’s disease: a 1-year prospective study. J. Neurol. 256, 1689–1695. doi:10.1007/s00415-009-5184-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Materia, S. I. (2024). BTS GAITLAB | Complete laboratory for multifactorial motion analysis. BTS. Available online at: https://www.btsbioengineering.com/products/bts-gaitlab/(Accessed February 7, 2024).

Google Scholar

Oh, J., Eltoukhy, M., Kuenze, C., Andersen, M. S., and Signorile, J. F. (2020). Comparison of predicted kinetic variables between Parkinson’s disease patients and healthy age-matched control using a depth sensor-driven full-body musculoskeletal model. Gait and Posture 76, 151–156. doi:10.1016/j.gaitpost.2019.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Paker, N., Bugdayci, D., Goksenoglu, G., Demircioğlu, D. T., Kesiktas, N., and Ince, N. (2015). Gait speed and related factors in Parkinson’s disease. J. Phys. Ther. Sci. 27, 3675–3679. doi:10.1589/jpts.27.3675

PubMed Abstract | CrossRef Full Text | Google Scholar

Paul, S. S., Canning, C. G., Sherrington, C., and Fung, V. S. C. (2012). Reduced muscle strength is the major determinant of reduced leg muscle power in Parkinson’s disease. Park. Relat. Disord. 18, 974–977. doi:10.1016/j.parkreldis.2012.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Peppe, A., Chiavalon, C., Pasqualetti, P., Crovato, D., and Caltagirone, C. (2007). Does gait analysis quantify motor rehabilitation efficacy in Parkinson’s disease patients? Gait and Posture 26, 452–462. doi:10.1016/j.gaitpost.2006.11.207

PubMed Abstract | CrossRef Full Text | Google Scholar

Pistacchi, M., Gioulis, M., Sanson, F., De Giovannini, E., Filippi, G., Rossetto, F., et al. (2017). Gait analysis and clinical correlations in early Parkinson’s disease. Funct. Neurol. 32, 28–34. doi:10.11138/FNeur/2017.32.1.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero-Flores, C. F., Bustamante-Bello, R., Moya Bencomo, M., Martinez-Ríos, E. A., and Montesinos, L. (2024). Optical marker-based motion capture of the human spine: a scoping review of study design and outcomes. Ann. Biomed. Eng. 52, 2373–2387. doi:10.1007/s10439-024-03567-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossanigo, R., Caruso, M., Bertuletti, S., Deriu, F., Knaflitz, M., Della Croce, U., et al. (2023). Base of support, step length and stride width estimation during walking using an inertial and infrared wearable system. Sensors (Basel) 23, 3921. doi:10.3390/s23083921

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, M., Amboni, M., Pisani, N., Volzone, A., Calderone, D., Barone, P., et al. (2025). Biomechanics parameters of gait analysis to characterize Parkinson’s disease: a scoping review. Sensors 25, 338. doi:10.3390/s25020338

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, M., Ricciardi, C., Amboni, M., Volzone, A., Barone, P., Romano, M., et al. (2023). “A cluster analysis for Parkinson’s disease phenotyping with gait parameters,” in 2023 IEEE international conference on metrology for eXtended reality, artificial intelligence and neural engineering (MetroXRAINE), 882–887. doi:10.1109/MetroXRAINE58569.2023.10405572

CrossRef Full Text | Google Scholar

Schniepp, R., Huppert, A., Decker, J., Schenkel, F., Schlick, C., Rasoul, A., et al. (2021). Fall prediction in neurological gait disorders: differential contributions from clinical assessment, gait analysis, and daily-life mobility monitoring. J. Neurol. 268, 3421–3434. doi:10.1007/s00415-021-10504-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartz, M. H., and Rozumalski, A. (2008). The gait deviation index: a new comprehensive index of gait pathology. Gait and Posture 28, 351–357. doi:10.1016/j.gaitpost.2008.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Scott, K., Bonci, T., Salis, F., Alcock, L., Buckley, E., Gazit, E., et al. (2022). Design and validation of a multi-task, multi-context protocol for real-world gait simulation. J. Neuroeng Rehabil. 19, 141. doi:10.1186/s12984-022-01116-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Speciali, D. S., Oliveira, E. M., Cardoso, J. R., Correa, J. C. F., Baker, R., and Lucareli, P. R. G. (2014). Gait profile score and movement analysis profile in patients with Parkinson’s disease during concurrent cognitive load. Braz J. Phys. Ther. 18, 315–322. doi:10.1590/bjpt-rbf.2014.0049

PubMed Abstract | CrossRef Full Text | Google Scholar

Tinetti, M. E. (1986). Performance-oriented assessment of mobility problems in elderly patients. J. Am. Geriatrics Soc. 34, 119–126. doi:10.1111/j.1532-5415.1986.tb05480.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Uhlig, M., and Prell, T. (2023a). Gait characteristics associated with fear of falling in hospitalized people with Parkinson’s disease. Sensors Basel, Switz. 23, 1111. doi:10.3390/s23031111

PubMed Abstract | CrossRef Full Text | Google Scholar

Uhlig, M., and Prell, T. (2023b). Gait characteristics associated with fear of falling in hospitalized people with Parkinson’s disease. Sensors Basel, Switz. 23, 1111. doi:10.3390/s23031111

PubMed Abstract | CrossRef Full Text | Google Scholar

Zanardi, A. P. J., da Silva, E. S., Costa, R. R., Passos-Monteiro, E., dos Santos, I. O., Kruel, L. F. M., et al. (2021). Gait parameters of Parkinson’s disease compared with healthy controls: a systematic review and meta-analysis. Sci. Rep. 11, 752. doi:10.1038/s41598-020-80768-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Jin, Y., Wang, M., Ji, C., Chen, Z., Fan, W., et al. (2024). The impact of anxiety on gait impairments in Parkinson’s disease: insights from sensor-based gait analysis. J. Neuroeng Rehabil. 21, 68. doi:10.1186/s12984-024-01364-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y. J., Wee, H. L., Chan, Y.-H., Seah, S. H., Au, W. L., Lau, P. N., et al. (2010). Progression of Parkinson’s disease as evaluated by Hoehn and Yahr stage transition times. Mov. Disord. 25, 710–716. doi:10.1002/mds.22875

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: gait analysis, biomechanics of gait, optoelectronic motion capture system, Parkinson’s disease, fear of falling, neurorehabilitation gait analysis, neurorehabilitation

Citation: De Pasquale P, Bonanno M, De Marchis C, Pergolizzi L, Lombardo Facciale A, Paladina G, Maggio MG, Impellizzeri F, Ciancarelli I, Quartarone A and Calabrò RS (2025) Beyond clinical scales: an observational study on instrumental gait analysis and biomechanical patterns in patients with Parkinson’s disease. Front. Bioeng. Biotechnol. 13:1541240. doi: 10.3389/fbioe.2025.1541240

Received: 07 December 2024; Accepted: 07 April 2025;
Published: 30 April 2025.

Edited by:

Marco Ghislieri, Polytechnic University of Turin, Italy

Reviewed by:

Marco Caruso, Polytechnic University of Turin, Italy
Michela Russo, University of Naples Federico II, Italy

Copyright © 2025 De Pasquale, Bonanno, De Marchis, Pergolizzi, Lombardo Facciale, Paladina, Maggio, Impellizzeri, Ciancarelli, Quartarone and Calabrò. 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: Mirjam Bonanno, bWlyamFtLmJvbmFubm9AaXJjY3NtZS5pdA==

These authors have contributed equally to this work

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.