Mathematical model of acetylcholine kinetics in neuroeffector junctions

Area of Science:

• Cardiovascular Physiology
• Neuropharmacology
• Mathematical Modeling

Background:

• Acetylcholine (ACh) kinetics in cardiac neuroeffector junctions (NEJ) are crucial for vagal heart rate regulation.
• Previous research suggested ACh concentration ([ACh]) in NEJ follows first-order linear kinetics.

Purpose of the Study:

• To mathematically investigate the underlying reasons for the observed linear kinetics of ACh in NEJ.
• To analyze the influence of diffusion, degradation, and receptor binding on ACh kinetics.

Main Methods:

• Mathematical modeling of ACh diffusion, degradation, and muscarinic receptor binding within NEJ.
• Identification and analysis of conditions influencing ACh kinetics.
• Comparison of ACh kinetics at pacemaker cells versus other NEJ sites.

Main Results:

• Seven specific conditions were identified that influence ACh kinetics in NEJ.
• Nonlinearity of ACh binding to muscarinic receptors minimally impacts [ACh] at pacemaker cells.
• [ACh] at pacemaker cells is proportional to vagal activity frequency and ACh release, divided by degradation rate.
• ACh kinetics are nearly first-order linear at pacemaker cells but not at other NEJ locations like nerve endings.

Conclusions:

• The apparent first-order linear kinetics of ACh observed in prior studies are specific to pacemaker cells in the sinus node NEJ.
• ACh kinetics differ significantly at other sites within the NEJ, such as nerve endings.
• Mathematical analysis clarifies the localized nature of linear ACh kinetics in cardiac autonomic control.