Unital quantum operators on the Bloch ball and Bloch region

For one qubit systems, we present a short, elementary argument characterizing unital quantum operators in terms of their action on Bloch vectors. We then show how our approach generalizes to multi-qubit systems, obtaining inequalities that govern when a ``diagonal'' superoperator on the Bloch region is a quantum operator. These inequalities are the n-qubit analogue of the Algoet-Fujiwara conditions. Our work is facilitated by an analysis of operator-sum decompositions in which negative summands are allowed.Comment: Revised and corrected, to appear in Physical Review

Mathematical Biology at an Undergraduate Liberal Arts College

Since 2002 we have offered an undergraduate major in Mathematical Biology at Harvey Mudd College. The major was developed and is administered jointly by the mathematics and biology faculty. In this paper we describe the major, courses, and faculty and student research and discuss some of the challenges and opportunities we have experienced

Characterising a universal cloning machine by maximum-likelihood estimation

We apply a general method for the estimation of completely positive maps to the 1-to-2 universal covariant cloning machine. The method is based on the maximum-likelihood principle, and makes use of random input states, along with random projective measurements on the output clones. The downhill simplex algorithm is applied for the maximisation of the likelihood functional.Comment: 5 pages, 2 figure

Predicting Outcomes of Prostate Cancer Immunotherapy by Personalized Mathematical Models

Therapeutic vaccination against disseminated prostate cancer (PCa) is partially effective in some PCa patients. We hypothesized that the efficacy of treatment will be enhanced by individualized vaccination regimens tailored by simple mathematical models.We developed a general mathematical model encompassing the basic interactions of a vaccine, immune system and PCa cells, and validated it by the results of a clinical trial testing an allogeneic PCa whole-cell vaccine. For model validation in the absence of any other pertinent marker, we used the clinically measured changes in prostate-specific antigen (PSA) levels as a correlate of tumor burden. Up to 26 PSA levels measured per patient were divided into each patient's training set and his validation set. The training set, used for model personalization, contained the patient's initial sequence of PSA levels; the validation set contained his subsequent PSA data points. Personalized models were simulated to predict changes in tumor burden and PSA levels and predictions were compared to the validation set. The model accurately predicted PSA levels over the entire measured period in 12 of the 15 vaccination-responsive patients (the coefficient of determination between the predicted and observed PSA values was R(2) = 0.972). The model could not account for the inconsistent changes in PSA levels in 3 of the 15 responsive patients at the end of treatment. Each validated personalized model was simulated under many hypothetical immunotherapy protocols to suggest alternative vaccination regimens. Personalized regimens predicted to enhance the effects of therapy differed among the patients.Using a few initial measurements, we constructed robust patient-specific models of PCa immunotherapy, which were retrospectively validated by clinical trial results. Our results emphasize the potential value and feasibility of individualized model-suggested immunotherapy protocols

Optimizing Combination Therapies with Existing and Future CML Drugs

Small-molecule inhibitors imatinib, dasatinib and nilotinib have been developed to treat Chromic Myeloid Leukemia (CML). The existence of a triple-cross-resistant mutation, T315I, has been a challenging problem, which can be overcome by finding new inhibitors. Many new compounds active against T315I mutants are now at different stages of development. In this paper we develop an algorithm which can weigh different combination treatment protocols according to their cross-resistance properties, and find the protocols with the highest probability of treatment success. This algorithm also takes into account drug toxicity by minimizing the number of drugs used, and their concentration. Although our methodology is based on a stochastic model of CML microevolution, the algorithm itself does not require measurements of any parameters (such as mutation rates, or division/death rates of cells), and can be used by medical professionals without a mathematical background. For illustration, we apply this algorithm to the mutation data obtained in [1], [2]

Deterministic mathematical modelling for cancer chronotherapeutics: cell population dynamics and treatment optimisation

Chronotherapeutics has been designed and used for more than twenty years as an effective treatment against cancer by a few teams around the world, among whom one of the first is Francis Lévi's at Paul-Brousse hospital (Villejuif, France), in application of circadian clock physiology to determine best infusion times within the 24-hour span for anticancer drug delivery. Mathematical models have been called in the last ten years to give a rational basis to such optimised treatments, for use in the laboratory and ultimately in the clinic. While actual clinical applications of the theoretical optimisation principles found have remained elusive so far to improve chronotherapeutic treatments in use, mathematical models provide proofs of concepts and tracks to be explored experimentally, to progress from theory to bedside. Starting from a simple ordinary differential equation model that allowed setting and numerically solving a drug delivery optimisation problem with toxicity constraints, this modelling enterprise has been extended to represent the division cycle in proliferating cell populations with different molecular targets, to allow for the representation of anticancer drug combinations that are used in clinical oncology. The main point to be made precise in such a therapeutic optimisation problem is to establish, here in the frame of circadian chronobiology, physiologically based differences between healthy and cancer cell populations in their responses to drugs. To this aim, clear biological evidence at the molecular level is still lacking, so that, starting from indirect observations at the experimental and clinical levels and from theoretical considerations on the model, speculations have been made, that will be exposed in this review of cancer chronotherapeutics models with the corresponding optimisation problems and their numerical solutions, to represent these differences between the two cell populations, with regard to circadian clock control

Mathematical models for immunology:current state of the art and future research directions

The advances in genetics and biochemistry that have taken place over the last 10 years led to significant advances in experimental and clinical immunology. In turn, this has led to the development of new mathematical models to investigate qualitatively and quantitatively various open questions in immunology. In this study we present a review of some research areas in mathematical immunology that evolved over the last 10 years. To this end, we take a step-by-step approach in discussing a range of models derived to study the dynamics of both the innate and immune responses at the molecular, cellular and tissue scales. To emphasise the use of mathematics in modelling in this area, we also review some of the mathematical tools used to investigate these models. Finally, we discuss some future trends in both experimental immunology and mathematical immunology for the upcoming years

Comment on: A Validated Mathematical Model of Cell-Mediated Immune Response to Tumor Growth

Gowal et al. erroneously surmise that in ref. 1, we assume “unrealistic” initial natural killer (NK) cell levels of zero in order to generate the simulations of Fig. 4, and they conclude from their own simulations that NK efficacy is stronger than we indicate. We do not draw conclusions about the efficiency of NK cells from our model. Rather, we use the model equations as a possible description of the different cytolytic mechanisms of NK cells versus CD8+ T cells. These mathematical descriptions of the immune cell dynamics are validated when we show that the simulations can reproduce experimentally observed differences in the efficacy of innate and the tumor-specific responses