Cancer is a worldwide burden in every aspect: physically, emotionally, and financially. A need for innovation in cancer research has led to a vast interdisciplinary effort to search for the next breakthrough. Mathematical modeling allows for a unique look into the underlying cellular dynamics and allows for testing treatment strategies without the need for clinical trials. This dissertation explores several iterations of a dendritic cell (DC) therapy model and correspondingly investigates what each iteration teaches about response to treatment.

In Chapter 2, motivated by the work of de Pillis et al. (2013), a mathematical model employing six ordinary differential (ODEs) and delay differential equations (DDEs) is formulated to understand the effectiveness of DC vaccines, accounting for cell trafficking with a blood and tumor compartment. A preliminary analysis is performed, with numerical simulations used to show the existence of oscillatory behavior. The model is then reduced to a system of four ODEs. Both models are validated using experimental data from melanoma-induced mice. Conditions under which the model admits rich dynamics observed in a clinical setting, such as periodic solutions and bistability, are established. Mathematical analysis proves the existence of a backward bifurcation and establishes thresholds for R0 that ensure tumor elimination or existence. A sensitivity analysis determines which parameters most significantly impact the reproduction number R0. Identifiability analysis reveals parameters of interest for estimation. Results are framed in terms of treatment implications, including effective combination and monotherapy strategies.

In Chapter 3, a study of whether the observed complexity can be represented with a simplified model is conducted. The DC model of Chapter 2 is reduced to a non-dimensional system of two DDEs. Mathematical and numerical analysis explore the impact of immune response time on the stability and eradication of the tumor, including an analytical proof of conditions necessary for the existence of a Hopf bifurcation. In a limiting case, conditions for global stability of the tumor-free equilibrium are outlined.

Lastly, Chapter 4 discusses future directions to explore. There still remain open questions to investigate and much work to be done, particularly involving uncertainty analysis. An outline of these steps is provided for future undertakings.

One explanation for membrane accommodation in response to a slowly rising current, and the phenomenon underlying the dynamics of elliptic bursting in nerves, is the mathematical problem of dynamic Hopf bifurcation. This problem has been studied extensively for linear (deterministic and stochastic) current ramps, nonlinear ramps, and elliptic bursting. These studies primarily investigated dynamic Hopf bifurcation in space-clamped excitable cells. In this study we introduce a new phenomenon associated with dynamic Hopf bifurcation. We show that for excitable spiny cables injected at one end with a slow current ramp, the generation of oscillations may occur an order one distance away from the current injection site. The phenomenon is significant since in the model the geometric and electrical parameters, as well as the ion channels, are uniformly distributed. In addition to demonstrating the phenomenon computationally, we analyze the problem using a singular perturbation method that provides a way to predict when and where the onset will occur in response to the input stimulus. We do not see this phenomenon for excitable cables in which the ion channels are embedded in the cable membrane itself, suggesting that it is essential for the channels to be isolated in the spines.

Dopamine (DA) is a neurotransmitter involved in attention, goal oriented behavior, movement, reward learning, and short term and working memory. For the past four decades, mathematical and computational modeling approaches have been useful in DA research, and although every modeling approach has limitations, a model is an efficient way to generate and explore hypotheses. This work develops a model of DA dynamics in a representative, single DA neuron by integrating previous experimental, theoretical and computational research. The model consists of three compartments: the cytosol, the vesicles, and the extracellular space and forms the basis of a new mathematical paradigm for examining the dynamics of DA synthesis, storage, release and reuptake. The model can be driven by action potentials generated by any model of excitable membrane potential or even from experimentally induced depolarization voltage recordings. Here the model is forced by a previously published model of the excitable membrane of a mesencephalic DA neuron in order to study the biochemical processes involved in extracellular DA production. After demonstrating that the model exhibits realistic dynamics resembling those observed experimentally, the model is used to examine the functional changes in presynaptic mechanisms due to application of cocaine. Sensitivity analysis and numerical studies that focus on various possible mechanisms for the inhibition of DAT by cocaine provide insight for the complex interactions involved in DA dynamics. In particular, comparing numerical results for a mixed inhibition mechanism to those for competitive, non-competitive and uncompetitive inhibition mechanisms reveals many behavioral similarities for these different types of inhibition that depend on inhibition parameters and levels of cocaine. Placing experimental results within this context of mixed inhibition provides a possible explanation for the conflicting views of uptake inhibition mechanisms found in experimental neuroscience literature.

In vertebrate outer retina, changes in the membrane potential of horizontal cells affect the calcium influx and glutamate release of cone photoreceptors via a negative feedback. This feedback has a number of important physiological consequences. One is called background-induced flicker enhancement (BIFE) in which the onset of dim background enhances the center flicker response of horizontal cells. The underlying mechanism for the feedback is still unclear but competing hypotheses have been proposed. One is the GABA hypothesis, which states that the feedback is mediated by gamma-aminobutyric acid (GABA), an inhibitory neurotransmitter released from horizontal cells. Another is the ephaptic hypothesis, which contends that the feedback is non-GABAergic and is achieved through the modulation of electrical potential in the intersynaptic cleft between cones and horizontal cells. In this study, a continuum spine model of the cone-horizontal cell synaptic circuitry is formulated. This model, a partial differential equation system, incorporates both the GABA and ephaptic feedback mechanisms. Simulation results, in comparison with experiments, indicate that the ephaptic mechanism is necessary in order for the model to capture the major spatial and temporal dynamics of the BIFE effect. In addition, simulations indicate that the GABA mechanism may play some minor modulation role.