Hana M. Dobrovolny

The restitution portrait: A new method for investigating rate-dependent restitution

Introduction: Electrical restitution, relating action potential duration (APD) to diastolic interval (DI), was believed to determine the stability of heart rhythm. However, recent studies demonstrate that stability also depends on long‐term APD changes caused by memory. This study presents a new method for investigation of rate‐ and memory‐dependent aspects of restitution and for assessment of mapping models of APD.

Assessing Mathematical Models of Influenza Infections Using Features of the Immune Response

The role of the host immune response in determining the severity and duration of an influenza infection is still unclear. In order to identify severity factors and more accurately predict the course of an influenza infection within a human host, an understanding of the impact of host factors on the infection process is required. Despite the lack of sufficiently diverse experimental data describing the time course of the various immune response components, published mathematical models were constructed from limited human or animal data using various strategies and simplifying assumptions. To assess the validity of these models, we assemble previously published experimental data of the dynamics and role of cytotoxic T lymphocytes, antibodies, and interferon and determined qualitative key features of their effect that should be captured by mathematical models. We test these existing models by confronting them with experimental data and find that no single model agrees completely with the variety of influenza viral kinetics responses observed experimentally when various immune response components are suppressed. Our analysis highlights the strong and weak points of each mathematical model and highlights areas where additional experimental data could elucidate specific mechanisms, constrain model design, and complete our understanding of the immune response to influenza.

Neuraminidase inhibitors for treatment of human and avian strain influenza: A comparative modeling study.

Treatment of seasonal influenza viral infections using antivirals such as neuraminidase inhibitors (NAIs) has been proven effective if administered within 48h post-infection. However, there is growing evidence that antiviral treatment of infections with avian-derived strains even as late as 6 days post-infection (dpi) can significantly reduce infection severity and duration. Using a mathematical model of in-host influenza viral infections which can capture the kinetics of both a short-lived, typical, seasonal infection and a severe infection exhibiting sustained viral titer, we explore differences in the effects of NAI treatment on both types of influenza viral infections. Comparison of our models behavior against experimental data from patients naturally infected with avian strains yields estimates for the times at which patients were infected that are consistent with those reported by the patients, and estimates of drug efficacies that are lower for patients who died than for those who recovered. In addition, our model suggests that the sustained, high, viral titers often seen in more severe influenza virus infections are the reason why antiviral treatment delayed by as much as 6 dpi will still lead to reduced viral titers and shortened illness. We conclude that NAIs may be an effective and beneficial treatment strategy against more severe strains of influenza virus characterized by high, sustained, viral titers. We believe that our mathematical model will be an effective tool in guiding treatment of severe influenza viral infections with antivirals.

Exploring Cell Tropism as a Possible Contributor to Influenza Infection Severity

Several mechanisms have been proposed to account for the marked increase in severity of human infections with avian compared to human influenza strains, including increased cytokine expression, poor immune response, and differences in target cell receptor affinity. Here, the potential effect of target cell tropism on disease severity is studied using a mathematical model for in-host influenza viral infection in a cell population consisting of two different cell types. The two cell types differ only in their susceptibility to infection and rate of virus production. We show the existence of a parameter regime which is characterized by high viral loads sustained long after the onset of infection. This finding suggests that differences in cell tropism between influenza strains could be sufficient to cause significant differences in viral titer profiles, similar to those observed in infections with certain strains of influenza A virus. The two target cell mathematical model offers good agreement with experimental data from severe influenza infections, as does the usual, single target cell model albeit with biologically unrealistic parameters. Both models predict that while neuraminidase inhibitors and adamantanes are only effective when administered early to treat an uncomplicated seasonal infection, they can be effective against more severe influenza infections even when administered late.

Coinfections of the Respiratory Tract: Viral Competition for Resources

Studies have shown that simultaneous infection of the respiratory tract with at least two viruses is common in hospitalized patients, although it is not clear whether these infections are more or less severe than single virus infections. We use a mathematical model to study the dynamics of viral coinfection of the respiratory tract in an effort to understand the kinetics of these infections. Specifically, we use our model to investigate coinfections of influenza, respiratory syncytial virus, rhinovirus, parainfluenza virus, and human metapneumovirus. Our study shows that during coinfections, one virus can block another simply by being the first to infect the available host cells; there is no need for viral interference through immune response interactions. We use the model to calculate the duration of detectable coinfection and examine how it varies as initial viral dose and time of infection are varied. We find that rhinovirus, the fastest-growing virus, reduces replication of the remaining viruses during a coinfection, while parainfluenza virus, the slowest-growing virus is suppressed in the presence of other viruses.

Differences in predictions of ODE models of tumor growth: a cautionary example

BackgroundWhile mathematical models are often used to predict progression of cancer and treatment outcomes, there is still uncertainty over how to best model tumor growth. Seven ordinary differential equation (ODE) models of tumor growth (exponential, Mendelsohn, logistic, linear, surface, Gompertz, and Bertalanffy) have been proposed, but there is no clear guidance on how to choose the most appropriate model for a particular cancer.MethodsWe examined all seven of the previously proposed ODE models in the presence and absence of chemotherapy. We derived equations for the maximum tumor size, doubling time, and the minimum amount of chemotherapy needed to suppress the tumor and used a sample data set to compare how these quantities differ based on choice of growth model.ResultsWe find that there is a 12-fold difference in predicting doubling times and a 6-fold difference in the predicted amount of chemotherapy needed for suppression depending on which growth model was used.ConclusionOur results highlight the need for careful consideration of model assumptions when developing mathematical models for use in cancer treatment planning.

Determining drug efficacy parameters for mathematical models of influenza

Antivirals are the first line of defence against influenza, so drug efficacy should be re-evaluated for each new strain. However, due to the time and expense involved in assessing the efficacy of drug treatments both in vitro and in vivo, treatment regimens are largely not re-evaluated even when strains are found to be resistant to antivirals. Mathematical models of the infection process can help in this assessment, but for accurate model predictions, we need to measure model parameters characterizing the efficacy of antivirals. We use computer simulations to explore whether in vitro experiments can be used to extract drug efficacy parameters for use in viral kinetics models. We find that the efficacy of neuraminidase inhibitors can be determined by measuring viral load during a single cycle assay, while the efficacy of adamantanes can be determined by measuring infected cells during the preparation stage for the single cycle assay.

High-Resolution High-Speed Panoramic Cardiac Imaging System

A panoramic cardiac imaging system consisting of three high-speed CCD cameras has been developed to image the surface electrophysiology of a rabbit heart via fluorescence imaging using a voltage-sensitive fluorescent dye. A robust, unique mechanical system was designed to accommodate the three cameras and to adapt to the requirements of future experiments. A unified computer interface was created for this application-a single workstation controls all three CCD cameras, illumination, stimulation, and a stepping motor that rotates the heart. The geometric reconstruction algorithms were adapted from a previous cardiac imaging system. We demonstrate the system by imaging a polymorphic cardiac tachycardia.

The in vivo efficacy of neuraminidase inhibitors cannot be determined from the decay rates of influenza viral titers observed in treated patients

Antiviral therapy is a first line of defence against new influenza strains. Current pandemic preparations involve stock- piling oseltamivir, an oral neuraminidase inhibitor (NAI), so rapidly determining the effectiveness of NAIs against new viral strains is vital for deciding how to use the stockpile. Previous studies have shown that it is possible to extract the drug efficacy of antivirals from the viral decay rate of chronic infections. In the present work, we use a nonlinear mathematical model representing the course of an influenza infection to explore the possibility of extracting NAI drug efficacy using only the observed viral titer decay rates seen in patients. We first show that the effect of a time-varying antiviral concentration can be accurately approximated by a constant efficacy. We derive a relationship relating the true treatment dose and time elapsed between doses to the constant drug dose required to approximate the time- varying dose. Unfortunately, even with the simplification of a constant drug efficacy, we show that the viral decay rate depends not just on drug efficacy, but also on several viral infection parameters, such as infection and production rate, so that it is not possible to extract drug efficacy from viral decay rate alone.

Assessing Uncertainty in A2 Respiratory Syncytial Virus Viral Dynamics

Respiratory syncytial virus (RSV) is the most common cause of bronchiolitis and pneumonia in children younger than 1 year of age in the United States. Moreover, RSV is being recognized more often as a significant cause of respiratory illness in older adults. Although RSV has been studied both clinically and in vitro, a quantitative understanding of the infection dynamics is still lacking. In this paper, we study the effect of uncertainty in the main parameters of a viral kinetics model of RSV. We first characterize the RSV replication cycle and extract parameter values by fitting the mathematical model to in vivo data from eight human subjects. We then use Monte Carlo numerical simulations to determine how uncertainty in the parameter values will affect model predictions. We find that uncertainty in the infection rate, eclipse phase duration, and infectious lifespan most affect the predicted dynamics of RSV. This study provides the first estimate of in vivo RSV infection parameters, helping to quantify RSV dynamics. Our assessment of the effect of uncertainty will help guide future experimental design to obtain more precise parameter values.