Optimizing cognitive neuroscience experiments for separating event- related fMRI BOLD responses in non-randomized alternating designs

Functional magnetic resonance imaging (fMRI) has revolutionized human brain research. But there exists a fundamental mismatch between the rapid time course of neural events and the sluggish nature of the fMRI blood oxygen level-dependent (BOLD) signal, which presents special challenges for cognitive neuroscience research. This limitation in the temporal resolution of fMRI puts constraints on the information about brain function that can be obtained with fMRI and also presents methodological challenges. Most notably, when using fMRI to measure neural events occurring closely in time, the BOLD signals may temporally overlap one another. This overlap problem may be exacerbated in complex experimental paradigms (stimuli and tasks) that are designed to manipulate and isolate specific cognitive-neural processes involved in perception, cognition, and action. Optimization strategies to deconvolve overlapping BOLD signals have proven effective in providing separate estimates of BOLD signals from temporally overlapping brain activity, but there remains reduced efficacy of such approaches in many cases. For example, when stimulus events necessarily follow a non-random order, like in trial-by-trial cued attention or working memory paradigms. Our goal is to provide guidance to improve the efficiency with which the underlying responses evoked by one event type can be detected, estimated, and distinguished from other events in designs common in cognitive neuroscience research. We pursue this goal using simulations that model the nonlinear and transient properties of fMRI signals, and which use more realistic models of noise. Our simulations manipulated: (i) Inter-Stimulus-Interval (ISI), (ii) proportion of so-called null events, and (iii) nonlinearities in the BOLD signal due to both cognitive and design parameters. We offer a theoretical framework along with a python toolbox called deconvolve to provide guidance on the optimal design parameters that will be of particular utility when using non-random, alternating event sequences in experimental designs. In addition, though, we also highlight the challenges and limitations in simultaneously optimizing both detection and estimation efficiency of BOLD signals in these common, but complex, cognitive neuroscience designs.