STRENGTHS AND LIMITATIONS OF THIS STUDY
An innovative serious game Virtual Reality (VR) software precisely emulating rehabilitation based on neuroproprioceptive ‘facilitation, inhibition’ principles.
Original approach to physiotherapy fusing VR, and physical therapy methods to facilitate physiological movement of impaired upper extremity function and sit-to-stand transitions.
Methodically near-identical physiotherapy content with VR enables identifying the contribution of the tested software on the final effect on the selected outcome measures.
Outcome measures include advanced functional MRI and tractography evaluating the brain functional connectivity and plasticity triggered by physical therapy intervention.
Unique exploratory assessment with molecular biomarkers of long non-coding RNA, peripheral blood leucocytes and dendritic cell analysis to document the neurobiological impact of the therapy.
Introduction
Upper limb (UL) mobility dysfunction affects 60%–75% of people with multiple sclerosis (PwMS), resulting in reduced participation in activities of daily living. 1 2 It has been documented that physical therapy (PT) brings important physical and psychological benefits. PT plays a key role in coping and optimising functional ability 3 and improvement of activity level. 4
Technologies such as virtual reality (VR) might be of benefit in rehabilitation, 5 but only a few studies focused on ULs in PwMS. Those indicated improvement in manual dexterity, coordination, precision, execution time and gross manual dexterity. 2 6 7
The current use of VR in PT is mainly based on the principles of sensorimotor learning. Game mechanisms focus on concentration and motivation, aiming at increased task performance repetitions. 8 9 Multisensory stimulation entrains the mirror neuron system, 10 and stored motor plans are activated, potentiating motor performance 11 and producing whole-body illusions. 10 VR through multifactorial stimulation influences dopamine centres in the brain. 12 Our team has developed a VR rehabilitation programme based on NFI physiotherapeutic methods. Systematic repetition induces a substantial cortical network reorganisation topographically closely related to the trained movement causing synaptogenesis. 8 13 NFI-based PT enhances synaptic connections among neurons forming functional networks evocating movement. A suitable combination of afferent stimuli repeatedly activates motor programmes at the subcortical level, inducing adaptive and plastic processes. 14 Brain plasticity changes following VR have been documented in poststroke people. 15 , 17 NFI approach affects brain structures involved in motor control, 1418 , 20 inducing structural changes. 19 21 Combining NFI and sensorimotor learning in VR has the potential to further induce adaptive and plastic processes. This is the first study to map the effect of VR PT using functional MRI (fMRI) and magnetic resonance tractography (MRT) in PwMS.
Molecular biomarkers could indicate improvements that support repair or compensatory strategies. Long non-coding RNAs (lncRNAs) have the potential to serve as a biomarker of MS activity degree, its progression phase and the response to disease-modifying drugs. 22 High expression of lncRNA MALAT1 inhibits neuronal apoptosis and lowers the proinflammatory interleukin 6 (IL-6). 23 High expression of lncRNA MEG3 contributes to the sustained proinflammatory state of microglial cells. 24 LncRNA H19 induces microglial apoptosis. 25 LncRNA GAS5 and miRNA miR-137 gene expression changes suggest a key regulatory role in MS pathogenesis. 26 Our group discovered lncRNA PARTICLE that takes part in the regulation of gene imprinting, epigenetic chromatin modification, 27 transcription interference and nuclear export. 28 29 We have identified selected lncRNAs (eg, MALAT1, MEG3, H19, GAS5, PARTICLE, TINA, etc) as suitable blood biomarkers for monitoring the efficacy of the intervention.
At the cellular level, this study analyses the composition of peripheral blood leucocytes Natural Killer (NK) cells, monocytes, T cells, B cells and their subpopulations) and the functional condition of generated monocyte-derived dendritic cells (MoDCs) by determining the expression of relevant surface molecules (see methods). The strongest genetic risk factor in MS is HLA-DRB1*15:01 allele encoding the beta chain of HLA (human leukocyte antigen)-DR molecule. 30 The research on the determination of dendritic cell (DC) populations is ongoing. 31 The conventional DCs play a pivotal role in the pathogenesis of MS 32 33 by cooperating with the autoreactive T cells. 34 The demyelination process is supported by the microglia and MoDCs. 35 The conventional DCs in the relapsing-remitting MS or the secondary progressive MS had a higher proportion of CD40 expression compared with the healthy controls. 36 In MS patients, DC frequency and maturity may correlate with disease stage and clinical presentation, but it has not been studied in the context of PT intervention.
Our team proposes a double-arm randomised controlled trial (RCT) that optimises and enhances VR-based PT PwMS with impaired upper extremity function, trunk coordination or transitions from a sitting position. We hypothesise that the use of VR in PT will have a higher effect on clinical functions than identical neuroproprioceptive ’facilitation, inhibition’ (NFI) PT without VR. Moreover, plastic and adaptive processes will be monitored using imaging methods fMRI, MRT and molecular genetic methods (long non-coding RNAs) and monitoring of immune system (the activation of DCs).
Methods and analysis
Study design
This study is designed as a single-blinded RCT with parallel assignment, registered under ID: NCT04807738 in clinical trials.gov. Participants will be randomly assigned into two groups of outpatient rehabilitation programme focusing on ULs, trunk control, stability in sitting and sit-to-stand ability. Both groups will receive physiotherapy based on the NFI principles. The intervention group engages in the VR with active involvement of the physiotherapist, while the control group participates in the same activities without the use of VR. Examination at baseline (E1) and after completion of the PT programme (E2) will consist of clinical examination assessing UL function and sit-to-stand transitions, fMRI, MRT and blood sample analysis.
Study setting
This study will be conducted at the Department of Rheumatology and Physiotherapy, Third Faculty of Medicine, Charles University and Thomayer University Hospital, Czech Republic and Department of Rehabilitation Medicine, Third Faculty of Medicine, Charles University, Prague, Czechia. Both workplaces are equipped with HTC Vive VR devices and tracking devices required for the VR-based PT. The therapeutic VR software for rehabilitation was developed and is further managed by the computer science researchers at the University of West Bohemia, Czech Republic.
Research and management team
The study team includes a coordinator, physiotherapists, occupational therapists (OTs), software developers, molecular biologists, MRI specialists and a statistician. The coordinator manages participant recruitment, eligibility assessment, allocation, assignment and study flow while providing participants with detailed study information. All personnel involved in assessments are blinded to treatment allocation. A blinded OT assessor conducts clinical measurements and administers questionnaire surveys. Molecular biologists handle the processing of peripheral blood samples, including peripheral blood mononuclear cell (PBMC) extraction and lncRNA expression analysis from in-vitro cultivated DCs. Blinded MRI specialists oversee fMRI examinations and data analysis. Four physiotherapists, specifically trained in the NFI-based PT method and certified in this technique, conduct both arms of the intervention. These physiotherapists, as they have an average of 2 years of clinical experience, work under the supervision of an experienced physiotherapist who developed the methodology (KR). The statistician prepares the dataset, will conduct regular audits, ensure data quality, and determine the appropriateness of parametric or non-parametric testing for significance analyses.
Participants
Sample size
Due to the limited number of PwMS patients, the planned number of participants is 70 (35 per group). This sample size is sufficient to detect a difference in the nine hole peg test (NHPT) and in other quantitative outcome measures (OMs) between the intervention and control group with an effect size of Cohen’s d=0.68 (a medium-to-large effect size), a power of 80% and a significance level of 5%, when using the two-tailed independent samples t-test. Expecting 10% of patients not meeting the inclusion criteria and a drop-out rate of 20%, 98 participants should be recruited to achieve the sample size of 70. The power analysis (pwr package 37 ) in R software was used for sample size calculation.
Patient recruitment
Patient recruitment is in cooperation with MS centres in Prague. The study is also advertised in local patient organisations in Prague (Unie Roska). Recruitment of the patients is continuous. All participants will be referred to the coordinator. The recruitment period is from April 2021 to January 2025.
Eligibility criteria
The inclusion criteria are age >18 and <75 years, diagnosis of multiple sclerosis according to the 2017 revised Mc Donald criteria, 38 an Expanded Disability Status Scale (EDSS) score ≥2 a ≤7 (27) (determined by a neurologist in the latest patient record) and in the past 3 months prior to recruitment no history of relapse, change in disease-modifying treatment nor corticosteroid therapy.
Exclusion criteria include other factors influencing mobility (history of stroke, pregnancy, traumatic injury of limb/s, severe cardiovascular or orthopaedic dysfunction or impaired cognitive functions during examination and/or consecutive therapy).
The following data will be collected on patient population characteristics: age, sex, weight and height, laterality, 39 EDSS, disease duration, MS therapy, use of mobility aids, history of falls and the International Classification of Functioning, Disability and Health (ICF) Multiple Sclerosis Brief Core set. 40
Randomisation, blinding
Patients eligible will be provided with the study details. After obtaining written consent, patients will be randomly assigned using simple randomisation (Random Allocation Software for parallel group randomised trials 41 ) and opaque sealed envelopes concealment. The allocation ratio is 1:1. Unblinding is permissible regarding participant safety, during the analysis, data interpretation or for regulatory compliance.
Retention and withdrawal, participation conditions
To promote adherence, the therapists will use their professional know-how and effective instruction to target each therapy to the current condition to increase adherence. Permitted concomitant care includes pharmacotherapy and other forms of supportive therapy—psychotherapy. Reasons for participant’s drop-out include discontinuation of the therapy for more than 2.5 weeks or five missed PTs, change in medication (by patient’s neurologist). The patient may step out at any time ( Figure 1 ).
Intervention
Participants in both groups will undergo 2 months (8 weeks, 60 min two times per week) of ambulatory NFI physiotherapy, combining the key principles of the Proprioceptive Neuromuscular Facilitation (PNF) and the Motor Programme Activating Therapy (MPAT) methods. In MPAT, suitable afferent somatosensory stimuli are applied, affecting the associative brain cortex and initiating desirable motor programmes. The repetition of these activated programmes through a set of stimuli adjusts posture with anatomical centration of the joints under various conditions. This will lead to better support for postural stabilisation while seated, getting up, stepping forward and standing. This will help participants to automatically apply the acquired motor skills in their daily lives. 42 PNF is a method used to learn effective movement patterns with high biomechanical efficiency based on repetitive stimulation of cooperating alpha-motoneurons and proprioceptors in muscles, tendons and joint capsules, using spiral and diagonal movements. 43
Experimental group: NFI in VR
The VR environment is programmed as one of the facilitation stimuli to attract attention and motivate participants to perform movements in the practised seven therapeutic elements while guiding them into qualitatively correct patterns. Participants will be instructed to focus on the suggested trajectory of movement, visualised in immersive stereoscopic 3D. Each movement is further specified with a location point, rotation, orientation or directional hint arrow. Wrist, arm and trunk sensors continuously log positions and 3D orientations using HTC Vive Tracker 2.0 devices. The patient’s movements are projected in the VR for immediate feedback and displayed to the physiotherapist overseeing the therapy.
Each therapeutic session will consist of seven therapeutic elements: the therapeutic session will begin with the activation of the motor programme for sitting according to MPAT. This activated position, UL movements will be subsequently trained through a ball game, a counterclockwise spiral, and the activation of the qualitatively correct function of the upper extremities based on the PNF principles (diagonal I. II. and II. with flexion and rotation of trunk), followed by practice in standing up and sitting down ( figure 2 ).
Each element will be repeated, with the number of repetitions gradually increasing based on the participant’s health condition—from seven repetitions in the first session to 12 by the end. The intensity is moderate. The physiotherapist actively engaged in the therapy, ensuring proper execution of tasks and providing facilitation, inhibition stimuli or support as needed. Emphasis was placed on movement quality, and if performed correctly, the therapist increased the activity’s difficulty using aids such as a TheraBand. Warm-up and cool-down exercises are not explicitly included in the protocol. Each therapeutic session will commence with feedback where participants will be asked about their postsession experience, including any perceived improvements or discomfort. In cases of musculoskeletal pain, the therapist will apply soft tissue techniques to address the issue. At the conclusion of each session, the therapist and participant will have time to review the recorded data, analyse the participant’s performance, highlight successful aspects and identify areas for correction. A more detailed description of the intervention is available in the online supplemental material 1 .
Active comparator group: NFI
The control group will receive therapy of equal duration, intensity and frequency (eight consecutive weeks, two times per week for 60 min), with the same seven therapeutic elements as outlined above, but without the use of VR. In the control group, the physiotherapist will also actively administer the NFI-based PT. The instructions include progressively increasing the number of repetitions based on individual performance (from seven to 12 per task).
Data collection and management
To promote the consistency of the data, instruments will be regularly calibrated. Both hands will be assessed. Separate Excel-based data sheets will be created for the baseline (E1) and the assessment after 2 months of therapy (E2; follow-up). All OMs will be examined within 1 week of the last therapy. The statistician and coordinator will conduct data vetting.
Outcome measures
( Figure 3 )
Primary OMs
Nine hole peg test (NHPT)
The NHPT is a standardised OM. It is a highly recommended OM. 44 Its psychometric properties include excellent test-retest reliability 45 and inter-rater and intrarater reliability 46 with a minimal clinically important difference of 15,3%. 47 The NHPT requires participants to repeatedly place and then remove nine pegs into nine holes, one at a time, as quickly as possible. Both ULs will be assessed with the NHPT (two trials per limb and averaged).
Box and block test (BBT)
The BBT is among the most used tests for assessing unilateral gross manual dexterity in MS 48 and 49 provides a valid assessment of UL function in MS with higher disability. Normative data are available for healthy adults. 50 The number of cubes transferred with one hand from one compartment of a box to another of equal size within 60 s is counted. Two trials in each hand will be performed and averaged.
Hand grip strength (HGS)
Jamar Hand Dynamometer is a validated tool that measures isometric grip force and strength, measured three times for each position available (9 cm, 12 cm, 14,5 cm, 17 cm and 20 cm) in a standardised position (sitting with neutrally rotated shoulder, elbow flexed at 90°, forearm in a neutral position). Jamar Hydraulic Pinch Gauge measures finger pinch force—key pinch, pinch tripod and pinch tip-tip. In each position, three trials will be measured in each hand.
Tremor
An accelerometer device placed on the index finger records tremor for 1 min while holding one hand raised forward in a standing or sitting position. First, we will record with open eyes, and a second record will be measured with closed eyes. Four trials in total will be recorded, one for each position and one for each arm.
fMAX—the spectral characteristic of postural tremor measured by the 3-axis accelerometer and 3-axis gyroscope chip (motion tracking sensor MPU-6050)—a frequency for which the smoothed power spectral density is maximal.
Pf1-f2—the spectral characteristic of postural tremor, a power of the signal in a band from f1 to f2—lower value, lower tremor.
Five times sit to stand test (5STS)
The 5TSTS test is a highly reliable tool for assessing lower limbs strength, balance control and mobility. 51 The time when the subject sits and stands repeatedly five times will be measured.
Secondary OM information
Multiple sclerosis impact scale (MSIS-29)
MSIS 29 is a highly recommended OM by MS EDGE and MSTF 46 for use in clinical trials and also in an outpatient setting, due to the consistency of the results. 52 A 29-item self-report measure with 20 items associated with a physical scale and nine items with a psychological scale. Summed responses are converted to a 0–100 scale, indicating the impact of the disease on daily function. 53
EQ-5D-5L health questionnaire (EuroQol)
EuroQol is a descriptive system and an established and validated generic instrument for assessing health-related quality of life 54 55 consisting of five dimensions.
Prespecified OMs
Kinematic analysis
The kinematic data from HCT Vive VR deliver suitable data quality for kinematic analysis 56 to track the learning process and the subject’s understanding of the proposed trajectory in examination E1 versus E2, providing feedback on the accuracy (distance from the outlined trajectory, in metres) and duration of the movement. It is possible to measure the spatial distance from the suggested movement trajectory (scaled to fit each patient’s proportions) averaged for the exercise, as well as angular deviation from the suggested orientation of each tracked locus.
fMRI examination and analysis
The fMRI examination substudy will be performed in the Magnetic Resonance (MR) system Siemens Vida 3T, using a 64-channel Radiofrequency (RF) head coil. MR examination (E1 and E2) will include measurements of functional stimulated data related to motoric inputs, functional connectivity measurement (resting state measurement, individually assessed), structural images (three-dimensional SPACE (sampling perfection with application-optimized contrasts using different flip angle evolutions) FLAIR (fluid-attenuated inversion recovery) and three-dimensional magnetisation prepared-RApid gradient echo (3D MP-RAGE)) and structural connectivity, with 106 spatial directions and 3b-factors allowing the reconstruction not only of fractional anisotropy but also the diffusion kurtosis maps. All functional data will be coregistered to structural images acquired with isotropic high spatial resolution.
The protocol includes two stimulation intervals. Imaginary execution of screened exercise (10 s), followed by a resting state of 10 s and a 10 s unilateral hand active movement-driven stimulation interval. The fMRI measurement includes 320 volumes (each covering the brain with 72 slices) acquired using the echo-planar imaging (EPI) sequence with multiband slice excitation (four simultaneous slices, SMS4) and the following basic parameters: TE of 30 ms, TR of 1250 ms (temporal resolution), voxel size of 1.9×1.9×1.9 mm 3 (spatial resolution). The resting state fMRI is measured with higher temporal resolution (TR of 650 ms) but slightly lower spatial resolution of 3×3×3 mm 3 .
The analysis will be performed using SPM 12 and CONN software. The fMRI data will be coregistered to individual structural 3D data with high spatial resolution (MP-RAGE with 1 mm 3 isotropic spatial resolution). Statistical analysis on an individual level uses general linear model statistics and threshold for final statistical maps is p=0.05 with family-wise error (FWE) correction. When testing the fMRI examination protocol, four healthy volunteers repeatedly underwent examination to calculate the reproducibility and repeatability of searched effects (size and strength of brain motor activations).
Diffusion-weighted imaging data
MR scans will be performed according to an optimised protocol for diffusion-weighted imaging using a spin-echo EPI (SE-EPI) sequence; repetition time (TR)=4500 ms, echo time (TE)=92 ms, 69 axial slices, field of view (FOV)=200 mm, 106 diffusion directions and a voxel size of 2×2×2 mm³. A multishell diffusion scheme with b-values of 0 s/mm², 1000 s/mm² and 3000 s/mm² will be employed for baseline and follow-up scans. To compare longitudinal changes in white matter tracts differential tractography 57 will be performed.
Immunological and epigenetic examination
Generation of MoDCs
A blood sample of 20 mL peripheral blood will be acquired in two 10 mL K 3 EDTA test tubes in examinations E1 and E2. PBMCs will be isolated by density gradient centrifugation using Ficoll-Paque media. Monocytic cell fraction will be separated by 2-hour cultivation in a 75 cm 2 plastic culture flask in complete medium (CM) containing Roswell Park Memorial Institute (RPMI) 1640 medium, 10% fetal bovine serum, 1% antibiotic-antimycotic solution at 37°C in a 5% CO 2 atmosphere. Adherent monocytes will be cultivated in CM supplemented with human granulocyte-macrophage colony‐stimulating factor and human IL‐4 for 6 days. MoDC maturation is induced by the addition of lipopolysaccharides to cells for 24 hours. Expression of the following CD cell markers will be analysed on leukocytes and MoDC: CD3, CD19, CD4, CD8, CD25, CD127, CD279, CD122, CD14, CD16, CD33, CD64, CD86, CD11c, HLA-DR and CD63 on a flow cytometer using antihuman monoclonal antibodies.
Analysis of lncRNA
Total RNA will be isolated from the second 10 mL K 3 EDTA test tube using the NucleoSpin RNA Midi blood Mini, according to the manual. Total RNA (1 µg) is converted into cDNA using standard protocol procedures and reagents (Thermo Fischer Scientific, Waltham, Massachusetts, USA).
Real-time quantitative PCR (qPCR)
The expression levels of preselected lncRNA ( MALAT1, MEG3, H19, GAS5, PARTICLE and TINA ) are determined by qPCR. For normalisation purposes and relative gene expression analysis, endogenous control genes are assessed ( GAPDH, B2M ) and measured in a StepOne eal-time PCR system instrument (VWR, Prague). All samples will be measured in triplets. The thermocycling conditions will be as follows: 50°C 2 min and, 95°C 10 min, followed by 40 cycles: 95°C 15 s and 60°C 1 min. Relative gene expression is quantified using the 2 (-ΔΔCt) method.
Statistical analysis and software
The intervention and control groups will be compared in their characteristics and OMs. The normality of quantitative (continuous or discrete) outcomes will be explored using a quantile-quantile plot and Shapiro–Wilk’s test. The homogeneity of variances in the groups will be assessed by Levene’s test. The differences in means between the groups are evaluated using the independent samples t-test and within the groups using the paired t-test. The two-way mixed analysis of variance (ANOVA, between-subjects: group, within-subjects: time) will be performed to summarise and visualise the group and time effects on quantitative outcomes. If the assumptions of the above-mentioned parametric tests are not met, the data will be transformed using a Box–Cox normal transformation or non-parametric versions of tests will be applied (Wilcoxon rank-sum test, Wilcoxon signed-rank test, the two-way mixed ANOVA for aligned rank transformed data). The differences in the distribution of qualitative (categorical) outcomes between the groups will be evaluated using Pearson’s χ 2 test or Fisher’s test, as appropriate, and within the groups using McNemar’s test. When multiple testing is done, procedures such as Bonferroni’s or Benjamini-Hochberg’s will be used to control for the so-called FWE.
We do not expect any significant differences in potential confounders (eg, the level of disability) between the groups at baseline. If they are, the statistical analysis will be adjusted (eg, using regression analysis, stratification or weighting) to control for selection bias. The attained level of significance (p value)<0.05 is considered statistically significant. R software is used for statistical data analysis. 58
Harms
The PT intervention will be carried out by qualified specialists actively approaching all subjects concerning individual disabilities. In VR, the potential harms of immersion include cybersickness and/or spatial disorientation. The physiotherapists will attentively and actively approach to prevent these situations, and most exercises for ULs are carried out while sitting. The programme will be individualised to prevent patient exertion or harm. All adverse events must be reported to the principal investigator and relevant specialist and will be stored.
Strengths and limitations of the study
Specificity, precision and technique of the proposed movement are key features to the therapeutic effect combining flexibility regarding content, number of repetitions and movement performance with a motivating VR setting. The tested VR features seven unilateral UL tasks. To target bilateral engagement of the UL 2 59 a VR software update is planned. Evaluation of VR contribution to the final effect is based on the intentional similarity of physiotherapy content in the real environment with VR PT. Possible sources of bias might occur from the simple randomisation process and include type, disease severity and MS stage, personal history and medication. Cognitive deficits regarding the participant’s ability to engage in the study will be evaluated by a treating neurologist. Participants will not be excluded for spasticity, as the therapeutic intent is to address it. However, its assessment could still be beneficial in future trials.
The duration of the therapeutic programme may be debatable concerning the potential to demonstrate clinical effects and plastic and adaptive changes in the CNS. However, based on our previous experience, our team believes this duration is sufficient and aligns with the lengths used in other studies. 760 , 63
The proposed methods offer a complex evaluation of new VR using novel methods. There is no outright consensus on the fMRI protocol setting, design and interpretation found in the literature. 64 Therefore, the proposed methodology could serve as a steppingstone to future comparisons of related RCTs with fMRI evaluation. The immunological and epigenetic substudy design serves as an exploration of molecular and cellular dynamics in response to PT, and it is a novelty within currently conducted PT studies in the multiple sclerosis patient population.
Results of this study have the potential to significantly contribute to the accessibility of innovative PT in outpatient care in PwMS, having a lifelong need for rehabilitation.
Ethics and dissemination
The multicentric ethics committee Faculty Hospital Kralovske Vinohrady and Medical Faculty of Charles University has approved the study under the code EK-VP/38/0/2021 (2 June 2021) and amendments (EK-VP/38/1/2021 on 4 April 2022). The ethics committee of the Institute for the Clinical and Experimental Medicine and Thomayer University Hospital have approved the study under No. 1983/21+4772/21 (G-21–02) (2 February 2021, and amendments (06958/22 (G21-26), 21 April 2022), based on submitted protocol, informed consent forms and participant education materials in Czech language (see English translation in the appendix). The estimated study completion is in January 2025. Study results will be disseminated to both the public and the scientific or medical community through publications in peer-reviewed journal articles and presentations at conferences, aiming to reach researchers and clinicians in the field. Anonymised data and research results may be shared on request for research purposes, such as meta-analyses or systematic reviews. Additionally, efforts will be made to communicate findings to the general public through collaborations with relevant organisations and the use of accessible briefs, ensuring broad accessibility and engagement.
supplementary material
Footnotes
Funding: The study was supported by the Grantová Agentura, Univerzita Karlova (Charles University Grant Agency) No 202322. The research was supported by Univerzita Karlova v Praze (Charles University), programme Cooperatio (Neuroscience and Medical Diagnostics and Basic Medical Sciences, the field ‘Medical Genetics’) No 260533/SVV/2024. The study was supported by SGS23/169/OHK3/3T/13 of the Czech Technical University in Prague. The study was supported (partly) by the long-term strategic development financing of the Institute of Computer Science (RVO:67985807) by Akademie Věd České Republiky (the Czech Academy of Sciences).
Prepublication history and additional supplemental material for this paper are available online. To view these files, please visit the journal online ( https://doi.org/10.1136/bmjopen-2024-088046 ).
Provenance and peer review: Not commissioned; externally peer reviewed.
Patient consent for publication: Consent obtained directly from patient(s).
Collaborators: Not applicable.
Patient and public involvement: Patients and/or the public were not involved in the design, conduct, reporting or dissemination plans of this research.
Author note: The authors declare that all claims expressed in this article are solely those of the authors. The study sponsors provided financial support for the study's execution, including materials and operational costs. The VR software and tremor analysing tool were developed in cooperation with the principal investigator and other team members. Any other 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 or manufacturer. A data monitoring committee was not established. Postclinical care is not contracted. All personal data are strictly confidential. Only pseudoanonymised data are analysed. Biological specimens are not stored. Access to the final trial dataset and statistical code will be restricted to the principal investigators and authorised research team members and will be made available upon reasonable request. Authorship for results will be granted if the coworker substantially contributes, drafts and revises manuscripts, is accountable and approves the final manuscript. The protocol SPIRIT 2013 and 2022 checklist is attached. Protocol V.V4/2024-10-24.
References
1. Johansson S, Ytterberg C, Claesson IM, et al. High concurrent presence of disability in multiple sclerosis. Associations with perceived health. J Neurol. 2007;254:767–73. doi: 10.1007/s00415-006-0431-5.
2. Bertoni R, Lamers I, Chen CC, et al. Unilateral and bilateral upper limb dysfunction at body functions, activity and participation levels in people with multiple sclerosis. Mult Scler. 2015;21:1566–74. doi: 10.1177/1352458514567553.
3. Buraga I, Davidescu I, Nica S, et al. Physical rehabilitation and quality of life in multiple sclerosis. Ro J Neurol. 2009;8:157–9. doi: 10.37897/RJN.2009.4.1.
4. Spooren AIF, Timmermans AAA, Seelen HAM. Motor training programs of arm and hand in patients with MS according to different levels of the ICF: a systematic review. BMC Neurol. 2012;12:49. doi: 10.1186/1471-2377-12-49.
5. Chen J, Or CK, Chen T. Effectiveness of Using Virtual Reality-Supported Exercise Therapy for Upper Extremity Motor Rehabilitation in Patients With Stroke: Systematic Review and Meta-analysis of Randomized Controlled Trials. J Med Internet Res. 2022;24:e24111. doi: 10.2196/24111.
6. Feys P, Coninx K, Kerkhofs L, et al. Robot-supported upper limb training in a virtual learning environment: a pilot randomized controlled trial in persons with MS. J Neuroeng Rehabil. 2015;12:60. doi: 10.1186/s12984-015-0043-3.
7. Norouzi E, Gerber M, Pühse U, et al. Combined virtual reality and physical training improved the bimanual coordination of women with multiple sclerosis. Neuropsychol Rehabil. 2021;31:552–69. doi: 10.1080/09602011.2020.1715231.
8. Kerr AL, Cheng SY, Jones TA. Experience-dependent neural plasticity in the adult damaged brain. J Commun Disord. 2011;44:538–48. doi: 10.1016/j.jcomdis.2011.04.011.
9. Chen J, Or CK, Li Z, et al. Effectiveness, safety and patients’ perceptions of an immersive virtual reality–based exercise system for poststroke upper limb motor rehabilitation: A proof-of-concept and feasibility randomized controlled trial. Dig Health. 2023;9 doi: 10.1177/20552076231203599.
10. Normand JM, Giannopoulos E, Spanlang B, et al. Multisensory stimulation can induce an illusion of larger belly size in immersive virtual reality. PLoS One. 2011;6:e16128. doi: 10.1371/journal.pone.0016128.
11. Calabrò RS, Naro A, Russo M, et al. The role of virtual reality in improving motor performance as revealed by EEG: a randomized clinical trial. J Neuroeng Rehabil. 2017;14:53. doi: 10.1186/s12984-017-0268-4.
12. Janssen J, Verschuren O, Renger WJ, et al. Gamification in Physical Therapy: More Than Using Games. Pediatr Phys Ther. 2017;29:95–9. doi: 10.1097/PEP.0000000000000326.
13. Straudi S, Severini G, Sabbagh Charabati A, et al. The effects of video game therapy on balance and attention in chronic ambulatory traumatic brain injury: an exploratory study. BMC Neurol. 2017;17:86. doi: 10.1186/s12883-017-0871-9.
14. Rasova K, Prochazkova M, Tintera J, et al. Motor programme activating therapy influences adaptive brain functions in multiple sclerosis: clinical and MRI study. Int J Rehabil Res. 2015;38:49–54. doi: 10.1097/MRR.0000000000000090.
15. Keller J, Štětkářová I, Macri V, et al. Virtual reality-based treatment for regaining upper extremity function induces cortex grey matter changes in persons with acquired brain injury. J Neuroeng Rehabil. 2020;17:127. doi: 10.1186/s12984-020-00754-7.
16. Patel J, Fluet G, Qiu Q, et al. Intensive virtual reality and robotic based upper limb training compared to usual care, and associated cortical reorganization, in the acute and early sub-acute periods post-stroke: a feasibility study. J Neuroeng Rehabil. 2019;16:92. doi: 10.1186/s12984-019-0563-3.
17. You SH, Jang SH, Kim YH, et al. Virtual reality-induced cortical reorganization and associated locomotor recovery in chronic stroke: an experimenter-blind randomized study. Stroke. 2005;36:1166–71. doi: 10.1161/01.STR.0000162715.43417.91.
18. Rasova K, Krasensky J, Havrdova E, et al. Is it possible to actively and purposely make use of plasticity and adaptability in the neurorehabilitation treatment of multiple sclerosis patients? A pilot project. Clin Rehabil. 2005;19:170–81. doi: 10.1191/0269215505cr831oa.
19. Prochazkova M, Tintera J, Spanhelova S, et al. Brain activity changes following neuroproprioceptive 'facilitation, inhibition' physiotherapy in multiple sclerosis: a parallel group randomized comparison of two approaches. Eur J Phys Rehabil Med. 2021;57:356–65. doi: 10.23736/S1973-9087.20.06336-4.
20. Angelova G, Skodova T, Prokopiusova T, et al. Ambulatory Neuroproprioceptive Facilitation and Inhibition Physical Therapy Improves Clinical Outcomes in Multiple Sclerosis and Modulates Serum Level of Neuroactive Steroids: A Two-Arm Parallel-Group Exploratory Trial. Life (Basel) 2020;10:1–16. doi: 10.3390/life10110267.
21. Bučková B, Kopal J, Řasová K, et al. Open Access: The Effect of Neurorehabilitation on Multiple Sclerosis-Unlocking the Resting-State fMRI Data. Front Neurosci. 2021;15 doi: 10.3389/fnins.2021.662784.
22. Nociti V, Santoro M. What do we know about the role of lncRNAs in multiple sclerosis? Neural Regen Res. 2021;16:1715–22. doi: 10.4103/1673-5374.306061.
23. Ma P, Li Y, Zhang W, et al. Long Non-coding RNA MALAT1 Inhibits Neuron Apoptosis and Neuroinflammation While Stimulates Neurite Outgrowth and Its Correlation With MiR-125b Mediates PTGS2, CDK5 and FOXQ1 in Alzheimer’s Disease. Curr Alzheimer Res. 2019;16:596–612. doi: 10.2174/1567205016666190725130134.
24. Meng J, Ding T, Chen Y, et al. LncRNA-Meg3 promotes Nlrp3-mediated microglial inflammation by targeting miR-7a-5p. Int Immunopharmacol. 2021;90:107141. doi: 10.1016/j.intimp.2020.107141.
25. Wan P, Su W, Zhang Y, et al. LncRNA H19 initiates microglial pyroptosis and neuronal death in retinal ischemia/reperfusion injury. Cell Death Differ. 2020;27:176–91. doi: 10.1038/s41418-019-0351-4.
26. Senousy MA, Shaker OG, Sayed NH, et al. LncRNA GAS5 and miR-137 Polymorphisms and Expression are Associated with Multiple Sclerosis Risk: Mechanistic Insights and Potential Clinical Impact. ACS Chem Neurosci. 2020;11:1651–60. doi: 10.1021/acschemneuro.0c00150.
27. O’Leary VB, Hain S, Maugg D, et al. Long non-coding RNA PARTICLE bridges histone and DNA methylation. Sci Rep. 2017;7 doi: 10.1038/s41598-017-01875-1.
28. Zhang X, Wang W, Zhu W, et al. Mechanisms and Functions of Long Non-Coding RNAs at Multiple Regulatory Levels. IJMS. 2019;20 doi: 10.3390/ijms20225573.
29. Khorkova O, Hsiao J, Wahlestedt C. Basic biology and therapeutic implications of lncRNA. Adv Drug Deliv Rev. 2015;87:15–24. doi: 10.1016/j.addr.2015.05.012.
30. Olerup O, Hillert J. HLA class II‐associated genetic susceptibility in multiple sclerosis: A critical evaluation. Tissue Antigens. 1991;38:1–15. doi: 10.1111/j.1399-0039.1991.tb02029.x.
31. Villani AC, Satija R, Reynolds G, et al. Single-cell RNA-seq reveals new types of human blood dendritic cells, monocytes, and progenitors. Science. 2017;356 doi: 10.1126/science.aah4573.
32. Mancini VSBW, Di Pietro AA, Pasquini LA. Microglia depletion as a therapeutic strategy: friend or foe in multiple sclerosis models? Neural Regen Res. 2023;18:267–72. doi: 10.4103/1673-5374.346538.
33. Mundt S, Mrdjen D, Utz SG, et al. Conventional DCs sample and present myelin antigens in the healthy CNS and allow parenchymal T cell entry to initiate neuroinflammation. Sci Immunol. 2019;4 doi: 10.1126/sciimmunol.aau8380.
34. Grigoriadis N, van Pesch V. A basic overview of multiple sclerosis immunopathology. Eur J Neurol. 2015;22:3–13. doi: 10.1111/ene.12798.
35. Zang YCQ, Skinner SM, Robinson RR, et al. Regulation of differentiation and functional properties of monocytes and monocyte-derived dendritic cells by interferon beta in multiple sclerosis. Mult Scler. 2004;10:499–506. doi: 10.1191/1352458504ms1081oa.
36. Karni A, Abraham M, Monsonego A, et al. Innate immunity in multiple sclerosis: myeloid dendritic cells in secondary progressive multiple sclerosis are activated and drive a proinflammatory immune response. J Immunol. 2006;177:4196–202. doi: 10.4049/jimmunol.177.6.4196.
37. Champely S, Claus E, Jeffrey G, et al. Basic Functions for Power Analysis. 2020.
38. Thompson AJ, Banwell BL, Barkhof F, et al. Diagnosis of multiple sclerosis: 2017 revisions of the McDonald criteria. Lancet Neurol. 2018;17:162–73. doi: 10.1016/S1474-4422(17)30470-2.
39. Oldfield RC. The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. 1971;9:97–113. doi: 10.1016/0028-3932(71)90067-4.
40. Coenen M, Cieza A, Freeman J, et al. The development of ICF Core Sets for multiple sclerosis: results of the International Consensus Conference. J Neurol. 2011;258:1477–88. doi: 10.1007/s00415-011-5963-7.
41. Saghaei M. Random allocation software for parallel group randomized trials. BMC Med Res Methodol. 2004;4 doi: 10.1186/1471-2288-4-26.
42. Řasová K, Bučková B, Prokopiusová T, et al. A Three-Arm Parallel-group Exploratory Trial documents balance improvement without much evidence of white matter integrity changes in people with multiple sclerosis following two months ambulatory neuroproprioceptive 'facilitation and inhibition' physical therapy. Eur J Phys Rehabil Med. 2021;57:889–99. doi: 10.23736/S1973-9087.21.06701-0.
43. Adler SS, Beckers D, Buck M. PNF in Practice. Berlin, Heidelberg: 2014. PNF in practice . http://link.springer.com/10.1007/978-3-642-34988-1 Available.
44. Cohen ET, Potter K, Allen DD, et al. Selecting Rehabilitation Outcome Measures for People with Multiple Sclerosis. Int J MS Care. 2015;17:181–9. doi: 10.7224/1537-2073.2014-067.
45. Hervault M, Balto JM, Hubbard EA, et al. Reliability, precision, and clinically important change of the Nine-Hole Peg Test in individuals with multiple sclerosis. Int J Rehabil Res. 2017;40:91–3. doi: 10.1097/MRR.0000000000000209.
46. Rosti-Otajärvi E, Hämäläinen P, Koivisto K, et al. The reliability of the MSFC and its components. Acta Neurol Scand. 2008;117:421–7. doi: 10.1111/j.1600-0404.2007.00972.x.
47. Jensen HB, Mamoei S, Ravnborg M, et al. Distribution-based estimates of minimum clinically important difference in cognition, arm function and lower body function after slow release-fampridine treatment of patients with multiple sclerosis. Mult Scler Relat Disord. 2016;7:58–60. doi: 10.1016/j.msard.2016.03.007.
48. Solaro C, Di Giovanni R, Grange E, et al. Box and block test, hand grip strength and nine‐hole peg test: correlations between three upper limb objective measures in multiple sclerosis. Eur J Neurology. 2020;27:2523–30. doi: 10.1111/ene.14427.
49. Solaro C, Di Giovanni R, Grange E, et al. Correlation between patient-reported manual ability and three objective measures of upper limb function in people with multiple sclerosis. Eur J Neurol. 2023;30:172–8. doi: 10.1111/ene.15560.
50. Mathiowetz V, Volland G, Kashman N, et al. Adult norms for the Box and Block Test of manual dexterity. Am J Occup Ther. 1985;39:386–91. doi: 10.5014/ajot.39.6.386.
51. Muñoz-Bermejo L, Adsuar JC, Mendoza-Muñoz M, et al. Test-Retest Reliability of Five Times Sit to Stand Test (FTSST) in Adults: A Systematic Review and Meta-Analysis. Biol (Basel) 2021;10:510. doi: 10.3390/biology10060510.
52. McGuigan C, Hutchinson M. The multiple sclerosis impact scale (MSIS-29) is a reliable and sensitive measure. J Neurol Neurosurg Psychiatry. 2004;75:266–9.
53. Hobart J, Lamping D, Fitzpatrick R. The Multiple Sclerosis Impact Scale (MSIS-29): A new patient-based outcome measure. Brain (Bacau) 2001;124:962–73. doi: 10.1093/brain/124.5.962.
54. Jongen PJ. Health-Related Quality of Life in Patients with Multiple Sclerosis: Impact of Disease-Modifying Drugs. CNS Drugs. 2017;31:585–602. doi: 10.1007/s40263-017-0444-x.
55. Visser LA, Louapre C, Uyl-de Groot CA, et al. Health-related quality of life of multiple sclerosis patients: a European multi-country study. Arch Public Health. 2021;79:1–12. doi: 10.1186/s13690-021-00561-z.
56. Spitzley KA, Karduna AR. Feasibility of using a fully immersive virtual reality system for kinematic data collection. 2019.
57. Yeh FC, Zaydan IM, Suski VR, et al. Differential tractography as a track-based biomarker for neuronal injury. Neuroimage. 2019;202 doi: 10.1016/j.neuroimage.2019.116131.
58. R Core Team R: a language and environment for statistical computing. 2023.
59. Ovacik U, Tarakci E, Gungor F, et al. The minnesota manual dexterity test as a bimanual performance measure in people with multiple sclerosis. Mult Scler Relat Disord. 2022;64:103943. doi: 10.1016/j.msard.2022.103943.
60. Cuesta-Gómez A, Martín-Díaz P, Sánchez-Herrera Baeza P, et al. Nintendo Switch Joy-Cons’ Infrared Motion Camera Sensor for Training Manual Dexterity in People with Multiple Sclerosis: A Randomized Controlled Trial. J Clin Med. 2022;11 doi: 10.3390/jcm11123261.
61. Doğan M, Ayvat E, Kılınç M. Telerehabilitation versus virtual reality supported task-oriented circuit therapy on upper limbs and trunk functions in patients with multiple sclerosis: A randomized controlled study. Mult Scler Relat Disord. 2023;71:104558. doi: 10.1016/j.msard.2023.104558.
62. Marcos-Antón S, Jardón-Huete A, Oña-Simbaña ED, et al. sEMG-controlled forearm bracelet and serious game-based rehabilitation for training manual dexterity in people with multiple sclerosis: a randomised controlled trial. J Neuroeng Rehabil. 2023;20 doi: 10.1186/s12984-023-01233-5.
63. Feys P, Coninx K, Kerkhofs L, et al. Robot-supported upper limb training in a virtual learning environment: a pilot randomized controlled trial in persons with MS. J Neuroeng Rehabil. 2015;12 doi: 10.1186/s12984-015-0043-3.
64. Tavazzi E, Cazzoli M, Pirastru A, et al. Neuroplasticity and Motor Rehabilitation in Multiple Sclerosis: A Systematic Review on MRI Markers of Functional and Structural Changes. Front Neurosci. 2021;15 doi: 10.3389/fnins.2021.707675.