Pronounced species divergence in corticospinal tract reorganization and functional recovery after lateralized spinal cord injury favors primates

Generating level-dependent models of cervical and thoracic spinal cord injury: Exploring the interplay of neuroanatomy, physiology, and function

The majority of spinal cord injuries (SCI) occur at the cervical level, which results in significant impairment. Neurologic level and severity of injury are primary endpoints in clinical trials; however, how level-specific damages relate to behavioural performance in cervical injury is incompletely understood. We hypothesized that ascending level of injury leads to worsening forelimb performance, and correlates with loss of neural tissue and muscle-specific neuron pools. A direct comparison of multiple models was made with injury realized at the C5, C6, C7 and T7 vertebral levels using clip compression with sham-operated controls. Animals were assessed for 10weeks post-injury with numerous (40) outcome measures, including: classic behavioural tests, CatWalk, non-invasive MRI, electrophysiology, histologic lesion morphometry, neuron counts, and motor compartment quantification, and multivariate statistics on the total dataset. Histologic staining and T1-weighted MR imaging revealed similar structural changes and distinct tissue loss with cystic cavitation across all injuries. Forelimb tests, including grip strength, F-WARP motor scale, Inclined Plane, and forelimb ladder walk, exhibited stratification between all groups and marked impairment with C5 and C6 injuries. Classic hindlimb tests including BBB, hindlimb ladder walk, bladder recovery, and mortality were not different between cervical and thoracic injuries. CatWalk multivariate gait analysis showed reciprocal and progressive changes forelimb and hindlimb function with ascending level of injury. Electrophysiology revealed poor forelimb axonal conduction in cervical C5 and C6 groups alone. The cervical enlargement (C5-T2) showed progressive ventral horn atrophy and loss of specific motor neuron populations with ascending injury. Multivariate statistics revealed a robust dataset, rank-order contribution of outcomes, and allowed prediction of injury level with single-level discrimination using forelimb performance and neuron counts. Level-dependent models were generated using clip-compression SCI, with marked and reliable differences in forelimb performance and specific neuron pool loss.

Developing a data sharing community for spinal cord injury research

The rapid growth in data sharing presents new opportunities across the spectrum of biomedical research. Global efforts are underway to develop practical guidance for implementation of data sharing and open data resources. These include the recent recommendation of 'FAIR Data Principles', which assert that if data is to have broad scientific value, then digital representations of that data should be Findable, Accessible, Interoperable and Reusable (FAIR). The spinal cord injury (SCI) research field has a long history of collaborative initiatives that include sharing of preclinical research models and outcome measures. In addition, new tools and resources are being developed by the SCI research community to enhance opportunities for data sharing and access. With this in mind, the National Institute of Neurological Disorders and Stroke (NINDS) at the National Institutes of Health (NIH) hosted a workshop on October 5-6, 2016 in Bethesda, MD, in collaboration with the Open Data Commons for Spinal Cord Injury (ODC-SCI) titled "Preclinical SCI Data: Creating a FAIR Share Community". Workshop invitees were nominated by the workshop steering committee (co-chairs: ARF and VPL; members: AC, KDA, MSB, KF, LBJ, PGP, JMS), to bring together junior and senior level experts including preclinical and basic SCI researchers from academia and industry, data science and bioinformatics experts, investigators with expertise in other neurological disease fields, clinical researchers, members of the SCI community, and program staff representing federal and private funding agencies. The workshop and ODC-SCI efforts were sponsored by the International Spinal Research Trust (ISRT), the Rick Hansen Institute, Wings for Life, the Craig H. Neilsen Foundation and NINDS. The number of attendees was limited to ensure active participation and feedback in small groups. The goals were to examine the current landscape for data sharing in SCI research and provide a path to its future. Below are highlights from the workshop, including perspectives on the value of data sharing in SCI research, workshop participant perspectives and concerns, descriptions of existing resources and actionable directions for further engaging the SCI research community in a model that may be applicable to many other areas of neuroscience. This manuscript is intended to share these initial findings with the broader research community, and to provide talking points for continued feedback from the SCI field, as it continues to move forward in the age of data sharing.

RNA-Seq Analysis of Microglia Reveals Time-Dependent Activation of Specific Genetic Programs following Spinal Cord Injury

Neurons have inherent competence to regrow following injury, although not spontaneously. Spinal cord injury (SCI) induces a pronounced neuroinflammation driven by resident microglia and infiltrating peripheral macrophages. Microglia are the first reactive glial population after SCI and participate in recruitment of monocyte-derived macrophages to the lesion site. Both positive and negative influence of microglia and macrophages on axonal regeneration had been reported after SCI, raising the issue whether their response depends on time post-lesion or different lesion severity. We analyzed molecular alterations in microglia at several time-points after different SCI severities using RNA-sequencing. We demonstrate that activation of microglia is time-dependent post-injury but is independent of lesion severity. Early transcriptomic response of microglia after SCI involves proliferation and neuroprotection, which is then switched to neuroinflammation at later stages. Moreover, SCI induces an autologous microglial expression of astrocytic markers with over 6% of microglia expressing glial fibrillary acidic protein and vimentin from as early as 72 h post-lesion and up to 6 weeks after injury. We also identified the potential involvement of DNA damage and in particular tumor suppressor gene breast cancer susceptibility gene 1 (Brca1) in microglia after SCI. Finally, we established that BRCA1 protein is specifically expressed in non-human primate spinal microglia and is upregulated after SCI. Our data provide the first transcriptomic analysis of microglia at multiple stages after different SCI severities. Injury-induced microglia expression of astrocytic markers at RNA and protein levels demonstrates novel insights into microglia plasticity. Finally, increased microglia expression of BRCA1 in rodents and non-human primate model of SCI, suggests the involvement of oncogenic proteins after CNS lesion.

A data-driven approach for evaluating multi-modal therapy in traumatic brain injury

Combination therapies targeting multiple recovery mechanisms have the potential for additive or synergistic effects, but experimental design and analyses of multimodal therapeutic trials are challenging. To address this problem, we developed a data-driven approach to integrate and analyze raw source data from separate pre-clinical studies and evaluated interactions between four treatments following traumatic brain injury. Histologic and behavioral outcomes were measured in 202 rats treated with combinations of an anti-inflammatory agent (minocycline), a neurotrophic agent (LM11A-31), and physical therapy consisting of assisted exercise with or without botulinum toxin-induced limb constraint. Data was curated and analyzed in a linked workflow involving non-linear principal component analysis followed by hypothesis testing with a linear mixed model. Results revealed significant benefits of the neurotrophic agent LM11A-31 on learning and memory outcomes after traumatic brain injury. In addition, modulations of LM11A-31 effects by co-administration of minocycline and by the type of physical therapy applied reached statistical significance. These results suggest a combinatorial effect of drug and physical therapy interventions that was not evident by univariate analysis. The study designs and analytic techniques applied here form a structured, unbiased, internally validated workflow that may be applied to other combinatorial studies, both in animals and humans.

Biodegradable scaffolds promote tissue remodeling and functional improvement in non-human primates with acute spinal cord injury

Tissue loss significantly reduces the potential for functional recovery after spinal cord injury. We previously showed that implantation of porous scaffolds composed of a biodegradable and biocompatible block copolymer of Poly-lactic-co-glycolic acid and Poly-l-lysine improves functional recovery and reduces spinal cord tissue injury after spinal cord hemisection injury in rats. Here, we evaluated the safety and efficacy of porous scaffolds in non-human Old-World primates (Chlorocebus sabaeus) after a partial and complete lateral hemisection of the thoracic spinal cord. Detailed analyses of kinematics and muscle activity revealed that by twelve weeks after injury fully hemisected monkeys implanted with scaffolds exhibited significantly improved recovery of locomotion compared to non-implanted control animals. Twelve weeks after injury, histological analysis demonstrated that the spinal cords of monkeys with a hemisection injury implanted with scaffolds underwent appositional healing characterized by a significant increase in remodeled tissue in the region of the hemisection compared to non-implanted controls. The number of glial fibrillary acidic protein immunopositive astrocytes was diminished within the inner regions of the remodeled tissue layer in treated animals. Activated macrophage and microglia were present diffusely throughout the remodeled tissue and concentrated at the interface between the preserved spinal cord tissue and the remodeled tissue layer. Numerous unphosphorylated neurofilament H and neuronal growth associated protein positive fibers and myelin basic protein positive cells may indicate neural sprouting inside the remodeled tissue layer of treated monkeys. These results support the safety and efficacy of polymer scaffolds in a primate model of acute spinal cord injury. A device substantially similar to the device described here is the subject of an ongoing human clinical trial.

Animal models of spinal cord injury: a systematic review

STUDY DESIGN

PRISMA-guided systematic review.

OBJECTIVES

To provide a comprehensive framework of the current animal models for investigating spinal cord injury (SCI) and categorize them based on the aims, patterns and levels of injury, and outcome measurements as well as animal species.

SETTING

Sina Trauma and Surgery Research Center, Tehran University of Medical Sciences, Tehran, Iran.

METHODS

An electronic search of the Medline database for literature describing animal models of SCI was performed on 1 January 2016 using the following keywords: 'spinal cord injuries' and 'animal models'. The search retrieved 2870 articles. Reviews and non-original articles were excluded. Data extraction was independently performed by two reviewers.

RESULTS

Among the 2209 included studies, testing the effects of drug's or growth factor's interventions was the most common aim (36.6%) followed by surveying pathophysiologic changes (30.2%). The most common spinal region involved was thoracic (81%). Contusion was the most common pattern of injury (41%) followed by transection (32.5%) and compression (19.4%). The most common species involved in animal models of SCI was the rat (72.4%). Two or more types of outcome assessments were used in the majority of the studies, and the most common assessment method was biological plus behavioral (50.8%).

CONCLUSIONS

Prior to choosing an animal model, the objectives of the proposed study must precisely be defined. Contusion and compression models better simulate the biomechanics and neuropathology of human injury, whereas transection models are valuable to study anatomic regeneration. Rodents are the most common and probably best-suited species for preliminary SCI studies.

Contribution of propriospinal neurons to recovery of hand dexterity after corticospinal tract lesions in monkeys

The direct cortico-motoneuronal connection is believed to be essential for the control of dexterous hand movements, such as precision grip in primates. It was reported, however, that even after lesion of the corticospinal tract (CST) at the C4-C5 segment, precision grip largely recovered within 1-3 mo, suggesting that the recovery depends on transmission through intercalated neurons rostral to the lesion, such as the propriospinal neurons (PNs) in the midcervical segments. To obtain direct evidence for the contribution of PNs to recovery after CST lesion, we applied a pathway-selective and reversible blocking method using double viral vectors to the PNs in six monkeys after CST lesions at C4-C5. In four monkeys that showed nearly full or partial recovery, transient blockade of PN transmission after recovery caused partial impairment of precision grip. In the other two monkeys, CST lesions were made under continuous blockade of PN transmission that outlasted the entire period of postoperative observation (3-4.5 mo). In these monkeys, precision grip recovery was not achieved. These results provide evidence for causal contribution of the PNs to recovery of hand dexterity after CST lesions; PN transmission is necessary for promoting the initial stage recovery; however, their contribution is only partial once the recovery is achieved.

A brain-spine interface alleviating gait deficits after spinal cord injury in primates

Spinal cord injury disrupts the communication between the brain and the spinal circuits that orchestrate movement. To bypass the lesion, brain–computer interfaces1–3 have directly linked cortical activity to electrical stimulation of muscles, which have restored grasping abilities after hand paralysis1,4. Theoretically, this strategy could also restore control over leg muscle activity for walking5. However, replicating the complex sequence of individual muscle activation patterns underlying natural and adaptive locomotor movements poses formidable conceptual and technological challenges6,7. Recently, we showed in rats that epidural electrical stimulation of the lumbar spinal cord can reproduce the natural activation of synergistic muscle groups producing locomotion8–10. Here, we interfaced leg motor cortex activity with epidural electrical stimulation protocols to establish a brain–spinal interface that alleviated gait deficits after a spinal cord injury in nonhuman primates. Rhesus monkeys were implanted with an intracortical microelectrode array into the leg area of motor cortex; and a spinal cord stimulation system composed of a spatially selective epidural implant and a pulse generator with real-time triggering capabilities. We designed and implemented wireless control systems that linked online neural decoding of extension and flexion motor states with stimulation protocols promoting these movements. These systems allowed the monkeys to behave freely without any restrictions or constraining tethered electronics. After validation of the brain–spinal interface in intact monkeys, we performed a unilateral corticospinal tract lesion at the thoracic level. As early as six days post-injury and without prior training of the monkeys, the brain–spinal interface restored weight-bearing locomotion of the paralyzed leg on a treadmill and overground. The implantable components integrated in the brain–spinal interface have all been approved for investigational applications in similar human research, suggesting a practical translational pathway for proof-of-concept studies in people with spinal cord injury.

Rehabilitation Strategies after Spinal Cord Injury: Inquiry into the Mechanisms of Success and Failure

Body-weight supported locomotor training (BWST) promotes recovery of load-bearing stepping in lower mammals, but its efficacy in individuals with a spinal cord injury (SCI) is limited and highly dependent on injury severity. While animal models with complete spinal transections recover stepping with step-training, motor complete SCI individuals do not, despite similarly intensive training. In this review, we examine the significant differences between humans and animal models that may explain this discrepancy in the results obtained with BWST. We also summarize the known effects of SCI and locomotor training on the muscular, motoneuronal, interneuronal, and supraspinal systems in human and non-human models of SCI and address the potential causes for failure to translate to the clinic. The evidence points to a deficiency in neuronal activation as the mechanism of failure, rather than muscular insufficiency. While motoneuronal and interneuronal systems cannot be directly probed in humans, the changes brought upon by step-training in SCI animal models suggest a beneficial re-organization of the systems' responsiveness to descending and afferent feedback that support locomotor recovery. The literature on partial lesions in humans and animal models clearly demonstrate a greater dependency on supraspinal input to the lumbar cord in humans than in non-human mammals for locomotion. Recent results with epidural stimulation that activates the lumbar interneuronal networks and/or increases the overall excitability of the locomotor centers suggest that these centers are much more dependent on the supraspinal tonic drive in humans. Sensory feedback shapes the locomotor output in animal models but does not appear to be sufficient to drive it in humans.

Advancing research in regeneration and repair of the motor circuitry: non-human primate models and imaging scales as the missing links for successfully translating injectable therapeutics to the clinic

Regeneration and repair is the ultimate goal of therapeutics in trauma of the central nervous system (CNS). Stroke and spinal cord injury (SCI) are two highly prevalent CNS disorders that remain incurable, despite numerous research studies and the clinical need for effective treatments. Neural engineering is a diverse biomedical field, that addresses these diseases using new approaches. Research in the field involves principally rodent models and biologically active, biodegradable hydrogels. Promising results have been reported in preclinical studies of CNS repair, demonstrating the great potential for the development of new treatments for the brain, spinal cord and peripheral nerve injury. Several obstacles stand in the way of clinical translation of neuroregeneration research. There seems to be a key gap in the translation of research from rodent models to human applications, namely non-human primate models, which constitute a critical bridging step. Applying injectable therapeutics and multimodal neuroimaging in stroke lesions using experimental rhesus monkey models is an avenue that a few research groups have begun to embark on. Understanding and assessing the changes that the injured brain or spinal cord undergoes after an intervention with biodegradable hydrogels in non-human primates seem to represent critical preclinical research steps. Existing innovative models in non-human primates allow us to evaluate the potential of neural engineering and injectable hydrogels. The results of these preliminary studies will pave the way for translating this research into much needed clinical therapeutic approaches. Cutting edge imaging technology using Connectome scanners represents a tremendous advancement, enabling the in vivo, detailed, high-resolution evaluation of these therapeutic interventions in experimental animals. Most importantly, they also allow quantifiable and clinically meaningful correlations with humans, increasing the translatability of these innovations to the bedside.

Intrinsic Axonal Growth and the Drive for Regeneration

Following damage to the adult nervous system in conditions like stroke, spinal cord injury, or traumatic brain injury, many neurons die and most of the remaining spared neurons fail to regenerate. Injured neurons fail to regrow both because of the inhibitory milieu in which they reside as well as a loss of the intrinsic growth capacity of the neurons. If we are to develop effective therapeutic interventions that promote functional recovery for the devastating injuries described above, we must not only better understand the molecular mechanisms of developmental axonal growth in hopes of re-activating these pathways in the adult, but at the same time be aware that re-activation of adult axonal growth may proceed via distinct mechanisms. With this knowledge in hand, promoting adult regeneration of central nervous system neurons can become a more tractable and realistic therapeutic endeavor.

Plasticity and Recovery After Dorsal Column Spinal Cord Injury in Nonhuman Primates

Here, we review recent work on plasticity and recovery after dorsal column spinal cord injury in nonhuman primates. Plasticity in the adult central nervous system has been established and studied for the past several decades; however, capacities and limits of plasticity are still under investigation. Studies of plasticity include assessing multiple measures before and after injury in animal models. Such studies are particularly important for improving recovery after injury in patients. In summarizing work by our research team and others, we suggest how the findings from plasticity studies in nonhuman primate models may affect therapeutic interventions for conditions involving sensory loss due to spinal cord injury.

Re-Establishment of Cortical Motor Output Maps and Spontaneous Functional Recovery via Spared Dorsolaterally Projecting Corticospinal Neurons after Dorsal Column Spinal Cord Injury in Adult Mice

Motor cortical plasticity contributes to spontaneous recovery after incomplete spinal cord injury (SCI), but the pathways underlying this remain poorly understood. We performed optogenetic mapping of motor cortex in channelrhodopsin-2 expressing mice to assess the capacity of the cortex to re-establish motor output longitudinally after a C3/C4 dorsal column SCI that bilaterally ablated the dorsal corticospinal tract (CST) containing ∼96% of corticospinal fibers but spared ∼3% of CST fibers that project via the dorsolateral funiculus. Optogenetic mapping revealed extensive early deficits, but eventual reestablishment of motor cortical output maps to the limbs at the same latency as preoperatively by 4 weeks after injury. Analysis of skilled locomotion on the horizontal ladder revealed early deficits followed by partial spontaneous recovery by 6 weeks after injury. To dissociate between the contributions of injured dorsal projecting versus spared dorsolateral projecting corticospinal neurons, we established a transient silencing approach to inactivate spared dorsolaterally projecting corticospinal neurons specifically by injecting adeno-associated virus (AAV)-expressing Cre-dependent DREADD (designer receptor exclusively activated by designer drug) receptor hM4Di in sensorimotor cortex and AAV-expressing Cre in C7/C8 dorsolateral funiculus. Transient silencing uninjured dorsolaterally projecting corticospinal neurons via activation of the inhibitory DREADD receptor hM4Di abrogated spontaneous recovery and resulted in a greater change in skilled locomotion than in control uninjured mice using the same silencing approach. These data demonstrate the pivotal role of a minor dorsolateral corticospinal pathway in mediating spontaneous recovery after SCI and support a focus on spared corticospinal neurons as a target for therapy.

SIGNIFICANCE STATEMENT

Spontaneous recovery can occur after incomplete spinal cord injury (SCI), but the pathways underlying this remain poorly understood. We performed optogenetic mapping of motor cortex after a cervical SCI that interrupts most corticospinal transmission but results in partial recovery on a horizontal ladder task of sensorimotor function. We demonstrate that the motor cortex can reestablish output to the limbs longitudinally. To dissociate the roles of injured and uninjured corticospinal neurons in mediating recovery, we transiently silenced the minor dorsolateral corticospinal pathway spared by our injury. This abrogated spontaneous recovery and resulted in a greater change in skilled locomotion than in uninjured mice using the same approach. Therefore, uninjured corticospinal neurons substantiate remarkable motor cortical plasticity and partial recovery after SCI.

Multidimensional Analysis of Magnetic Resonance Imaging Predicts Early Impairment in Thoracic and Thoracolumbar Spinal Cord Injury

Literature examining magnetic resonance imaging (MRI) in acute spinal cord injury (SCI) has focused on cervical SCI. Reproducible systems have been developed for MRI-based grading; however, it is unclear how they apply to thoracic SCI. Our hypothesis is that MRI measures will group as coherent multivariate principal component (PC) ensembles, and that distinct PCs and individual variables will show discriminant validity for predicting early impairment in thoracic SCI. We undertook a retrospective cohort study of 25 patients with acute thoracic SCI who underwent MRI on admission and had American Spinal Injury Association Impairment Scale (AIS) assessment at hospital discharge. Imaging variables of axial grade, sagittal grade, length of injury, thoracolumbar injury classification system (TLICS), maximum canal compromise (MCC), and maximum spinal cord compression (MSCC) were collected. We performed an analytical workflow to detect multivariate PC patterns followed by explicit hypothesis testing to predict AIS at discharge. All imaging variables loaded positively on PC1 (64.3% of variance), which was highly related to AIS at discharge. MCC, MSCC, and TLICS also loaded positively on PC2 (22.7% of variance), while variables concerning cord signal abnormality loaded negatively on PC2. PC2 was highly related to the patient undergoing surgical decompression. Variables of signal abnormality were all negatively correlated with AIS at discharge with the highest level of correlation for axial grade as assessed with the Brain and Spinal Injury Center (BASIC) score. A multiple variable model identified BASIC as the only statistically significant predictor of AIS at discharge, signifying that BASIC best captured the variance in AIS within our study population. Our study provides evidence of convergent validity, construct validity, and clinical predictive validity for the sampled MRI measures of SCI when applied in acute thoracic and thoracolumbar SCI.

Mechanical Design and Analysis of a Unilateral Cervical Spinal Cord Contusion Injury Model in Non-Human Primates

Non-human primate (NHP) models of spinal cord injury better reflect human injury and provide a better foundation to evaluate potential treatments and functional outcomes. We combined finite element (FE) and surrogate models with impact data derived from in vivo experiments to define the impact mechanics needed to generate a moderate severity unilateral cervical contusion injury in NHPs (Macaca mulatta). Three independent variables (impactor displacement, alignment, and pre-load) were examined to determine their effects on tissue level stresses and strains. Mechanical measures of peak force, peak displacement, peak energy, and tissue stiffness were analyzed as potential determinants of injury severity. Data generated from FE simulations predicted a lateral shift of the spinal cord at high levels of compression (>64%) during impact. Submillimeter changes in mediolateral impactor position over the midline increased peak impact forces (>50%). Surrogate cords established a 0.5 N pre-load protocol for positioning the impactor tip onto the dural surface to define a consistent dorsoventral baseline position before impact, which corresponded with cerebrospinal fluid displacement and entrapment of the spinal cord against the vertebral canal. Based on our simulations, impactor alignment and pre-load were strong contributors to the variable mechanical and functional outcomes observed in in vivo experiments. Peak displacement of 4 mm after a 0.5N pre-load aligned 0.5-1.0 mm over the midline should result in a moderate severity injury; however, the observed peak force and calculated peak energy and tissue stiffness are required to properly characterize the severity and variability of in vivo NHP contusion injuries.

Induced Pluripotent Stem Cell Therapies for Cervical Spinal Cord Injury

Cervical-level injuries account for the majority of presented spinal cord injuries (SCIs) to date. Despite the increase in survival rates due to emergency medicine improvements, overall quality of life remains poor, with patients facing variable deficits in respiratory and motor function. Therapies aiming to ameliorate symptoms and restore function, even partially, are urgently needed. Current therapeutic avenues in SCI seek to increase regenerative capacities through trophic and immunomodulatory factors, provide scaffolding to bridge the lesion site and promote regeneration of native axons, and to replace SCI-lost neurons and glia via intraspinal transplantation. Induced pluripotent stem cells (iPSCs) are a clinically viable means to accomplish this; they have no major ethical barriers, sources can be patient-matched and collected using non-invasive methods. In addition, the patient’s own cells can be used to establish a starter population capable of producing multiple cell types. To date, there is only a limited pool of research examining iPSC-derived transplants in SCI—even less research that is specific to cervical injury. The purpose of the review herein is to explore both preclinical and clinical recent advances in iPSC therapies with a detailed focus on cervical spinal cord injury.

Mechanisms Underlying the Neuromodulation of Spinal Circuits for Correcting Gait and Balance Deficits after Spinal Cord Injury

Epidural electrical stimulation of lumbar segments facilitates standing and walking in animal models and humans with spinal cord injury. However, the mechanisms through which this neuromodulation therapy engages spinal circuits remain enigmatic. Using computer simulations and behavioral experiments, we provide evidence that epidural electrical stimulation interacts with muscle spindle feedback circuits to modulate muscle activity during locomotion. Hypothesis-driven strategies emerging from simulations steered the design of stimulation protocols that adjust bilateral hindlimb kinematics throughout gait execution. These stimulation strategies corrected subject-specific gait and balance deficits in rats with incomplete and complete spinal cord injury. The conservation of muscle spindle feedback circuits across mammals suggests that the same mechanisms may facilitate motor control in humans. These results provide a conceptual framework to improve stimulation protocols for clinical applications.