Proposed research model · Version 2.0

PRĀṆA–PD

Prāṇa-Regulated Adaptive Network Architecture for Parkinson's Disease. An integrative mathematical model, inspired by the Praśna Upaniṣad and Ayurvedic principles, for studying multisystem dysfunction and functional recovery in Parkinson's disease (PD).

Run the simulator Read the research paper Clinical protocol Sākṣī companion app Python & code

What this site contains

  1. The framework (this page): corrected five-Prāṇa model, equations, interventions, dashboard, research programme.
  2. An executable model: the master equation as an interactive simulator and a validated Python implementation, including a coupling-identification experiment whose limitations are reported openly.
  3. A journal-style research paper and a clinical research protocol for neurologists and Ayurvedic physicians.
  4. Sākṣī, your uploaded companion app, integrated unchanged in function (two small fixes, listed in Code).

1. One fundamental correction: Prāṇa does not "approach zero"

The original proposal states

\[ \lim_{\text{Disease Progression}}\mathcal P(t)=0 \]

This is philosophically interesting, but it cannot be established as a biological description of Parkinson's disease. Praśna Upaniṣad III.4–III.7 describes the principal Prāṇa as assigning functions to the other Prāṇas, like a king appointing officials. We keep that coordination principle as the philosophical foundation and replace the untestable limit with a measurable quantity: the Prāṇa-inspired physiological coordination reserve

\[ 0\le \rho(t)\le 1 \]

PD does not necessarily degrade all systems together. Some people have prominent motor difficulty, others prominent gastrointestinal, autonomic, sleep or speech problems. So the model describes selective disruption of individual systems and of their connections, rather than assuming one global factor decays to zero.

2. The five-Prāṇa physiological architecture

The five-dimensional state vector is retained, but each component is an estimated functional state with its own measurable indicators, never the literal quantity of a subtle energy.

\[ \mathbf{\Psi}(t)=\big[\psi_1,\dots,\psi_5\big]^{\!\top}=\big[P_{\text{Prāṇa}},P_{\text{Apāna}},P_{\text{Samāna}},P_{\text{Vyāna}},P_{\text{Udāna}}\big]^{\!\top} \]

Prāṇa

Sensory & respiratory regulation

Tradition: eyes, ears, mouth, nose.

  • Respiratory rate and pattern
  • Pulmonary function when indicated
  • Olfactory function
  • Breath support; coordination during speech

Smell loss is a recognised non-motor feature of PD, but smell testing alone is not a diagnostic test.

Apāna

Elimination & pelvic autonomic

Tradition: downward, elimination-related functions.

  • Bowel frequency; constipation severity
  • Bristol Stool Form Scale
  • Urinary urgency, frequency, emptying
  • Pelvic-floor function

Samāna

Digestion & nutritional homeostasis

Tradition: digestion, assimilation, distribution.

  • Gastric emptying where indicated
  • Appetite, intake, weight change
  • Nutritional screening; GI symptoms
  • Medication-related digestive tolerability

Vyāna

Circulation & distributed motor coordination

Tradition: pervasive distribution and movement.

  • Gait speed, balance, walking capacity
  • Upper- and lower-limb movement
  • Orthostatic blood pressure
  • Heart-rate variability (exploratory)

Vyāna ↔ cardiovascular/motor physiology is a hypothesis, not an established equivalence.

Udāna

Speech, swallowing & expiratory function

Tradition: upward-moving functions.

  • Voice intensity; speech clarity
  • Swallowing function (EAT-10, clinical assessment)
  • Cough effectiveness; aspiration risk
  • Respiratory coordination for speech/swallow
Modelling rule against double counting. Prāṇa and Udāna overlap in respiration. We separate respiratory-sensory regulation (Prāṇa) from speech, swallowing and airway protection (Udāna), and never let one measurement feed two domains.

3. The five-system mathematical model

Let \(\mathbf x(t)\in[0,1]^5\) be normalised functional capacity: 1 is a selected healthy reference, 0 is the lower endpoint of the defined assessment scale. These are research scores, not physiological quantities of Prāṇa.

3.1 The central equation

The original \(\dot{\mathbf\Psi}=\mathcal P\mathbf A\mathbf\Psi+\mathbf E(t)\) is expanded to:

\[ \frac{d\mathbf x}{dt}=\Big[\mathbf A_L+\rho(t)\,\mathbf A_C-\mathbf D(t)\Big]\mathbf x+\mathbf B\mathbf u(t)+\mathbf G\mathbf E(t)+\mathbf w(t) \]
TermInterpretation
\(\mathbf A_L\)Individual subsystem dynamics
\(\rho(t)\)Prāṇa-inspired coordination reserve
\(\mathbf A_C\)Interconnection and feedback between systems
\(\mathbf D(t)\)Domain-specific functional impairment
\(\mathbf B\), \(\mathbf u(t)\)Influence of, and inputs from, supervised therapeutic interventions
\(\mathbf E(t)\), \(\mathbf G\)Environmental and pathological disturbances, and sensitivity to them
\(\mathbf w(t)\)Unexplained system variation

Three distinct phenomena are modelled: (1) individual system function, (2) coordination between systems, (3) response to therapeutic and environmental influences.

Three implementation choices made while making it executable

Written exactly, \(\dot{\mathbf x}=[\mathbf A_L+\rho\mathbf A_C-\mathbf D]\mathbf x\) with no input has \(\mathbf x=\mathbf 0\) as its only equilibrium, so the healthy state would not be stable even with \(\mathbf D=0\). To make the equation computable we added the following, and they are explicit modelling decisions, open to review:

  1. Tonic homeostatic drive \(\mathbf h=\mathbf\Lambda\mathbf 1\) with \(\mathbf A_L=-\mathbf\Lambda\), so the healthy state \(\mathbf x=\mathbf 1\) is an equilibrium when \(\mathbf D=0\).
  2. \(\mathbf A_C\) is a graph Laplacian (zero row sums): coordination redistributes function between systems; it neither creates nor destroys it.
  3. Impairment calibration: \(d_i=\lambda_i s_i/(1-s_i)\), so that an uncoupled, untreated domain settles at \(x_i^\ast=1-s_i\), where \(s_i\in[0,1)\) is the clinical impairment score.

3.2 The 5 × 5 coupling matrix

Instead of assuming independent Prāṇas, we investigate their relationships. The off-diagonal entries are hypothesised and must be identified from longitudinal data, never assigned arbitrarily.

\[ \mathbf A(t)=\mathbf A_L+\mathbf C(t)\odot\mathbf A_C-\mathbf D(t) \]

\(\mathbf C(t)\) is a dynamic coupling matrix and \(\odot\) is element-wise multiplication, so one connection can deteriorate while another stays preserved.

RelationshipScientific research question
Prāṇa ↔ UdānaDoes breath support influence speech and swallowing function?
Apāna ↔ SamānaAre bowel function and gastric transit related?
Samāna ↔ VyānaAre nutritional status and physical capacity associated?
Vyāna ↔ UdānaAre physical and respiratory-muscular capacities associated?
Prāṇa ↔ VyānaIs respiratory function associated with exercise capacity?
Greater coupling is not always better. Excessive synchronisation can itself be pathological. The aim is adaptive, appropriately regulated coordination, not maximising every coefficient. In the simulator's steady-state analysis, raising \(\rho\) lifts the weakest domain but pulls the average down, because coupling redistributes rather than creates function (see the trade-off chart on the Simulator).

4. The central nervous system as a regulatory network

A sixth concept is added so that the five-Prāṇa framework connects to modern neurology without claiming Prāṇa is a neurological structure.

Central nervous system
Basal ganglia · Brainstem · Thalamocortical networks · Autonomic regulation
▼ bidirectional regulatory relationships — research hypothesis ▲

Prāṇa

Sensory and breathing

Apāna

Elimination and pelvic

Samāna

Digestive and metabolic

Vyāna

Circulatory and motor

Udāna

Speech and swallowing
▼
Patient functional outcomes
Movement · Digestion · Communication · Independence · Quality of life

Central research question: can measurable changes in neurological network function help explain differences in coordination between the five physiological domains? This can be investigated with neurological examination, autonomic testing, gait analysis, physiological recordings and, in suitable research settings, neuroimaging.

5. The disease-disruption model

\[ \mathbf D(t)=\operatorname{diag}\big(d_1(t),d_2(t),d_3(t),d_4(t),d_5(t)\big) \]

\(d_1\ldots d_5\) are the Prāṇa, Apāna, Samāna, Vyāna and Udāna domain impairments, each estimated from clinical measurement. Example: a patient with marked constipation and urinary urgency, preserved mobility and no swallowing difficulty may show \(d_2>d_4\) and \(d_5\approx 0\). That means the Apāna-related domain has greater measured impairment than the Udāna-related domain. It does not mean one Prāṇa is biologically depleted.

5.1 Individualised patient profiles (illustrative, invented to demonstrate the method)

Higher score = greater impairment. These are not clinical assessments. They are the same profiles preloaded in the Sākṣī clinician view and the simulator.

ProfilePrāṇaApānaSamānaVyānaUdāna
A — predominantly motor2015207525
B — autonomic & gastrointestinal1575703015
C — speech & swallowing3510252580

A single treatment plan or intervention score would not represent these differences, which is why a personalised functional model is worth investigating.

6. The recovery model: from disease management to functional restoration

Three goals must be kept separate: (1) improving symptoms; (2) restoring functional ability and independence; (3) preventing or reversing the underlying neurodegeneration. The first two can be pursued through established treatment and rehabilitation. The third remains an unsolved research problem, and there is currently no established cure for PD. We therefore define recovery as:

Regaining, maintaining or improving physiological function and quality of life, even when the underlying disease remains present.

6.1 Target vector and intervention vector

\[ \mathbf e(t)=\mathbf x^{\ast}-\mathbf x(t),\qquad \mathbf u(t)=\big[u_{\text{medical}},u_{\text{movement}},u_{\text{nutrition}},u_{\text{respiratory}},u_{\text{speech}},u_{\text{mind-body}},u_{\text{Ayurvedic adjunct}}\big]^{\!\top} \]

Each component is a clinician-approved intervention category, not a dose or prescription.

6.2 Adaptive intervention equation

\[ \mathbf u(t)=\operatorname{sat}\Big[\mathbf K\big(\mathbf x^{\ast}-\hat{\mathbf x}(t)\big)\Big] \]

\(\mathbf K\): intervention-response relationship. \(\hat{\mathbf x}\): estimated state. \(\operatorname{sat}\): safety constraint on intensity. In practice: measure, identify the key difficulties, build a multidisciplinary plan, measure the response, adjust with the treating professionals. This is a research control-system design. It must not independently change PD medication, exercise intensity or herbal treatments. In the executable model the medical channel is fixed by the clinician and the controller is never allowed to move it.

6.3 Recovery objective

\[ J=\sum_{t=1}^{T}\Big[(\mathbf x^{\ast}-\mathbf x_t)^{\!\top}\mathbf Q(\mathbf x^{\ast}-\mathbf x_t)+\mathbf u_t^{\!\top}\mathbf R\,\mathbf u_t+\lambda S_t\Big] \]

\(\mathbf Q\) prioritises patient-important outcomes, \(\mathbf R\) penalises intervention burden, \(S\) is safety risk, \(\lambda\) is the weight on safety. In plain language: maximise meaningful functional improvement while minimising treatment burden and safety risk. The weights must be elicited from clinicians and patients, not chosen arbitrarily.

7. Interventions for each Prāṇa domain

Established PD care is kept separate from traditional practices that still need evaluation.

DomainEvidence-informed clinical approachAdditional Vedic–Ayurvedic researchKey research questions
PrāṇaRespiratory assessment, exercise, sensory supportGentle, adapted breathing practice and meditationCan supervised gentle breathing improve respiratory comfort? Is respiratory function associated with speech performance? Can meditation reduce perceived stress?
ApānaBowel and bladder management, hydration, nutritionIndividualised dietary assessment and routineCan individualised nutrition reduce constipation? Are meal timing and fluid intake related to bowel function?
SamānaNutritional support, gastrointestinal assessmentĀhāra principles and individualised dietary researchHow does GI function relate to medication response? Is gastric emptying associated with symptom fluctuation?
VyānaPhysiotherapy, walking, balance, strength trainingAdapted yoga and movement practicesDo gait and coordination improve with regular exercise? Does physical capacity track the coordination index?
UdānaSpeech-language therapy, swallowing rehabilitationCarefully adapted vocal and breathing practicesCan structured voice therapy improve intelligibility? Can swallowing therapy improve safety? Can respiratory-muscle training improve cough in suitable patients?

Safety boundaries

  • Breathing programmes must be adapted to cardiopulmonary status. Forceful breathing and breath retention should not be prescribed for PD without clinical assessment.
  • Persistent constipation must be medically assessed, not assumed to be "Vata imbalance".
  • Medication and protein timing for levodopa is individual advice from the neurologist and dietitian, not a generic rule.
  • Exercise is individualised to fall risk and disease stage.
  • Speech and swallowing therapy are clinical services, not optional wellness.

Red flags: seek prompt medical assessment

New choking, coughing during meals, a wet-sounding voice, repeated chest infections or unexplained weight loss need prompt medical and swallowing assessment. Breathing difficulty or severe choking needs emergency care.

8. Ayurvedic variables: a second layer

\[ \mathbf{A}(t)=\big[\text{Vata},\text{Pitta},\text{Kapha},\text{Agni},\text{Āma},\text{Dhātu},\text{Ojas}\big]^{\!\top},\qquad \dot{\mathbf x}=\mathbf f(\mathbf x,\mathbf A,\rho,\mathbf u,\mathbf D,\mathbf E) \]

These are traditional assessment constructs that need standardised definitions and reliability testing before they can enter a model. Two levels result: the five Prāṇas (functional systems and their coordination) and the seven Ayurvedic constructs (individualised traditional assessment and therapeutic context). Vata and the Apāna functional score may be associated, but they remain independently measured variables. The scientific test is whether Ayurvedic assessment improves prediction beyond conventional neurological and functional measures.

9. The PRĀṆA–PD measurement dashboard

DomainSuggested measurements
PrāṇaRespiratory assessments, olfactory assessment, relevant sensory measures
ApānaBristol Stool Form Scale, bowel diary, SCOPA-AUT
SamānaNutritional assessment, GI symptom measures, gastric emptying where indicated
VyānaMDS-UPDRS III, 10-metre walk test, Timed Up and Go, blood pressure (incl. orthostatic)
UdānaSpeech intelligibility, EAT-10, clinical swallowing assessment, instrumental testing when indicated
Overall outcomesMDS-UPDRS, PDQ-39, activities of daily living, patient-defined goals

Analysis: correlation and multivariate regression; longitudinal mixed-effects models; dynamic network analysis; state-space modelling; Bayesian inference; machine learning with independent validation. The same patient should be assessed under consistent medication timing and testing conditions so that normal fluctuations are not mistaken for improvement or deterioration. The Sākṣī app is the prototype capture instrument: it pairs measured and reported scores for each domain and surfaces discordance without labelling it "poor awareness".

10. How to prove whether the model works

Indicative study-planning timeline, not a completed or approved clinical trial.

  1. Months 1–3 · Textual and scientific foundations. Review the Praśna Upaniṣad, Ayurvedic concepts and the current PD literature. Define each mapping and its limits.
  2. Months 4–6 · Measurement development. Build the five-domain instrument. Review by Ayurvedic experts, neurologists, physiotherapists and measurement scientists.
  3. Months 7–9 · Pilot observational study. Repeated clinical measurements in people with PD. Inter-domain associations and reliability.
  4. Months 10–12 · Mathematical validation. Estimate the coupling matrix, coordination reserve and impairment parameters. Test reproducibility and prediction.
  5. Months 13–15 · Intervention feasibility. Feasibility, adherence and safety of adding standardised complementary practices to usual care.
  6. Months 16–18 · Independent review and trial design. Prepare a larger randomised controlled study if justified.

Eighteen months can support model and feasibility development. It cannot prove disease reversal or long-term neuroprotection.

CodeHypothesis
H1Five-Prāṇa functional domains can be operationalised with reliable measurements.
H2Some domains are statistically associated with specific PD symptoms.
H3Inter-domain coordination measures add information beyond individual symptom scores.
H4Individualised multidisciplinary rehabilitation improves selected functional outcomes.
H5Some standardised Ayurvedic adjuncts give additional benefit over conventional care alone.
H6The combined model predicts individual functional outcomes better than a conventional baseline model.
On H5, note the evidence. A 2025 systematic review (Hammoud et al.) found only five clinical trials of Mucuna pruriens in PD, 108 participants in total, with quality limitations. Mucuna contains L-DOPA, so it is not an innocuous supplement to add on top of prescribed PD medication. Some Ayurvedic preparations have contained toxic amounts of lead, mercury or arsenic, so product quality and independent safety testing are essential. The protocol therefore excludes L-DOPA-containing adjuncts from the feasibility phase.

11. Five scientific contributions the model could pursue

1 · Multisystem coordination measurement

Can PD functional problems be described by measurable relationships between five domains?

2 · Personalised assessment

A model that recognises individual functional profiles instead of treating PD as one uniform condition.

3 · Functional restoration

Can a coordinated rehabilitation programme improve several areas of functioning at once?

4 · Evidence-based Ayurveda integration

Which traditional practices give measurable added benefit, for whom, and under what safety conditions?

5 · Predictive monitoring

Can longitudinal data flag worsening and signal when professional assessment is needed?

12. Final scientific architecture: three layers

Layer 1 · Vedic–Ayurvedic foundation

Five Prāṇas + Ayurveda: Prāṇa · Apāna · Samāna · Vyāna · Udāna · Doṣas · Agni · Dhātu

Layer 2 · Neurological and mathematical foundation

Coordination + network dynamics: central regulatory networks · dynamic coupling · disease-related impairment · physiological measurements

Layer 3 · Functional restoration

Adaptive rehabilitation: individualised clinical care · exercise · nutrition · speech · swallowing · complementary approaches

Ultimate research outcome: better function, greater independence, improved quality of life, and evidence-based identification of potentially beneficial integrative interventions.

The final master equation

\[ \frac{d\mathbf x}{dt}=\underbrace{\Big[\mathbf A_L+\rho(t)\mathbf A_C-\mathbf D(t)\Big]\mathbf x}_{\text{functional system dynamics}}+\underbrace{\mathbf B\mathbf u(t)}_{\text{intervention}}+\underbrace{\mathbf G\mathbf E(t)}_{\text{external influences}}+\mathbf w(t) \]

subject to \(0\le\rho(t)\le1\), \(\mathbf x(t)\in\mathcal X_{\text{safe}}\), and the individualised objective

\[ \min_{\mathbf u}\sum_{t=1}^{T}\Big[\|\mathbf x^{\ast}-\mathbf x(t)\|_{\mathbf Q}^{2}+\|\mathbf u(t)\|_{\mathbf R}^{2}+\lambda S(t)\Big] \]

13. What has been built, and what the numbers do and do not show

Verified

  • The master equation is implemented twice (TypeScript for the web simulator, Python with NumPy), and the two agree to four decimals on all nine deterministic profile × arm scenarios.
  • With \(\mathbf D=0\) and no input, the healthy state is an equilibrium; with coupling as a Laplacian, coordination conserves total function.
  • A coupling-identification experiment runs end to end and returns confidence intervals.

Not shown, and why it matters

  • The simulated advantage of the adaptive arm over usual care exists because the placeholder \(\mathbf B\) matrix assumes interventions work. It illustrates the control logic. It is not evidence of efficacy.
  • In a synthetic cohort with known coupling, the estimator recovers the sign of all five links, but three of five are biased upward by about 0.04–0.05 (roughly 60–90 %) even at n = 200, while the other two are near-unbiased. The estimator needs recorded disturbances and instrumenting, and its remaining bias is unresolved. Real data would need a pre-registered simulation-based calibration before any coupling value is interpreted.
  • Identification was tested on weekly observations only. Monthly clinic data is untested and likely weaker.

14. The three deliverables, built

Academic paper

Journal-style manuscript with abstract, theory, methods, simulation results, limitations and references. Print to PDF from the browser.

Open paper

Clinical protocol

For neurologists and Ayurvedic physicians: phased design, eligibility, schedule of assessments, safety rules, analysis plan.

Open protocol

Model, simulation, Python

Interactive simulator, downloadable Python implementation, Monte-Carlo and identifiability scripts.

Simulator Code

The aim is to help people living with Parkinson's regain and preserve as much physical, cognitive, social and functional independence as possible. The five-Prāṇa framework is offered as a research hypothesis for studying that aim. Any claim that it cures PD, restores lost neurons or reverses disease progression must wait for convincing clinical evidence.