In silico Identification of Key Factors Driving the Response of Muscle Sensory Neurons to Noxious Stimuli

Nociceptive nerve endings embedded in muscle tissue transduce peripheral noxious stimuli into an electrical signal [i.e., an action potential (AP)] to initiate pain sensations. A major contributor to nociception from the muscles is mechanosensation. However, due to the heterogeneity in the expression of proteins, such as ion channels, pumps, and exchangers, on muscle nociceptors, we currently do not know the relative contributions of different proteins and signaling molecules to the neuronal response due to mechanical stimuli. In this study, we employed an integrated approach combining a customized experimental study in mice with a computational model to identify key proteins that regulate mechanical nociception in muscles. First, using newly collected data from somatosensory recordings in mouse hindpaw muscles, we developed and then validated a computational model of a mechanosensitive mouse muscle nociceptor. Next, by performing global sensitivity analyses that simulated thousands of nociceptors, we identified three ion channels (among the 17 modeled transmembrane proteins and four endoplasmic reticulum proteins) as potential regulators of the nociceptor response to mechanical forces in both the innocuous and noxious range. Moreover, we found that simulating single knockouts of any of the three ion channels, delayed rectifier voltage-gated K(+) channel (Kv1.1) or mechanosensitive channels Piezo2 or TRPA1, considerably altered the excitability of the nociceptor (i.e., each knockout increased or decreased the number of triggered APs compared to when all channels were present). These results suggest that altering expression of the gene encoding Kv1.1, Piezo2, or TRPA1 might regulate the response of mechanosensitive muscle nociceptors.