Model-based analysis of influenza A virus replication in genetically engineered cell lines elucidates the impact of host cell factors on key kinetic parameters of virus growth

Mechanistic modeling of viral particle production

Viral systems such as wild-type viruses, viral vectors, and virus-like particles are essential components of modern biotechnology and medicine. Despite their importance, the commercial-scale production of viral systems remains highly inefficient for multiple reasons. Computational strategies are a promising avenue for improving process development, optimization, and control, but require a mathematical description of the system. This article reviews mechanistic modeling strategies for the production of viral particles, both at the cellular and bioreactor scales. In many cases, techniques and models from adjacent fields such as epidemiology and wild-type viral infection kinetics can be adapted to construct a suitable process model. These process models can then be employed for various purposes such as in-silico testing of novel process operating strategies and/or advanced process control.

Impact of Influenza A Virus Infection on Growth and Metabolism of Suspension MDCK Cells Using a Dynamic Model

Cell cultured-based influenza virus production is a viable option for vaccine manufacturing. In order to achieve a high concentration of viable cells, is requirement to have not only optimal process conditions, but also an active metabolism capable of intracellular synthesis of viral components. Experimental metabolic data collected in such processes are complex and difficult to interpret, for which mathematical models are an appropriate way to simulate and analyze the complex and dynamic interaction between the virus and its host cell. A dynamic model with 35 states was developed in this study to describe growth, metabolism, and influenza A virus production in shake flask cultivations of suspension Madin-Darby Canine Kidney (MDCK) cells. It considers cell growth (concentration of viable cells, mean cell diameters, volume of viable cells), concentrations of key metabolites both at the intracellular and extracellular level and virus titers. Using one set of parameters, the model accurately simulates the dynamics of mock-infected cells and correctly predicts the overall dynamics of virus-infected cells for up to 60 h post infection (hpi). The model clearly suggests that most changes observed after infection are related to cessation of cell growth and the subsequent transition to apoptosis and cell death. However, predictions do not cover late phases of infection, particularly for the extracellular concentrations of glutamate and ammonium after about 12 hpi. Results obtained from additional in silico studies performed indicated that amino acid degradation by extracellular enzymes resulting from cell lysis during late infection stages may contribute to this observed discrepancy.

Absolute quantification of viral proteins during single-round replication of MDCK suspension cells

Madin-Darby canine kidney (MDCK) cells are widely used in basic research and for the propagation of influenza A viruses (IAV) for vaccine production. To identify targets for antiviral therapies and to optimize vaccine manufacturing, a detailed understanding of the viral life cycle is important. This includes the characterization of virus entry, the synthesis of the various viral RNAs and proteins, the transfer of viral compounds in the cell and virus budding. In case quantitative information is available, the analysis can be complemented by mathematical modelling approaches. While comprehensive studies focusing on IAV entry as well as viral mRNA, vRNA and cRNA accumulation in the nucleus of cells have been performed, quantitative data regarding IAV protein synthesis and accumulation was mostly lacking. In this study, we present a mass spectrometry (MS)-based method to evaluate whether an absolute quantification of viral proteins is possible for single-round replication in suspension MDCK cells. Using influenza A/PR/8/34 (H1N1, RKI) as a model strain at a multiplicity of infection of ten, defined amounts of isotopically labelled peptides of synthetic origin of four IAV proteins (hemagglutinin, neuraminidase, nucleoprotein, matrix protein 1) were added as an internal standard before tryptic digestion of samples for absolute quantification (AQUA). The first intracellular protein detected was NP at 1 h post infection (hpi). A maximum extracellular concentration of 7.7E+12 copies/mL was achieved. This was followed by hemagglutinin (3 hpi, maximum 4.1E+12 copies/mL at 13 hpi), matrix protein 1 (5 hpi, maximum 2.2E+12 copies/mL at 13 hpi) and neuraminidase (5 hpi, 6.0E+11 copies/mL at 13 hpi). In sum, for the first time absolute IAV protein copy numbers were quantified by a MS-based method for infected MDCK cells providing important insights into viral protein dynamics during single-round virus replication. SIGNIFICANCE: Influenza A virus is a significant human pathogen worldwide. To improve therapies against influenza and overcome bottlenecks in vaccine production in cell culture, it is critical to gain a detailed understanding of the viral life cycle. In addition to qPCR-based models, this study will examine the dynamics of influenza virus proteins during infection of producer cells to gain initial insights into changes in absolute copy numbers.

Imaging Flow Cytometry and Confocal Immunofluorescence Microscopy of Virus-Host Cell Interactions

Viruses are diverse pathogens that use host factors to enter cells and cause disease. Imaging the entry and replication phases of viruses and their interactions with host factors is key to fully understanding viral infections. This review will discuss how confocal microscopy and imaging flow cytometry are used to investigate virus entry and replication mechanisms in fixed and live cells. Quantification of viral images and the use of cryo-electron microscopy to gather structural information of viruses is also explored. Using imaging to understand how viruses replicate and interact with host factors, we gain insight into cellular processes and identify novel targets to develop antiviral therapeutics and vaccines.

Impact of Oseltamivir Treatment on Influenza A and B Virus Dynamics in Human Volunteers

Influenza viruses infect millions of humans every year causing an estimated 400,000 deaths globally. Due to continuous virus evolution current vaccines provide only limited protection against the flu. Several antiviral drugs are available to treat influenza infection, and one of the most commonly used drugs is oseltamivir (Tamiflu). While the mechanism of action of oseltamivir as a neuraminidase inhibitor is well-understood, the impact of oseltamivir on influenza virus dynamics in humans has been controversial. Many clinical trials with oseltamivir have been done by pharmaceutical companies such as Roche but the results of these trials until recently have been provided as summary reports or papers. Typically, such reports included median virus shedding curves for placebo and drug-treated influenza virus infected volunteers often indicating high efficacy of the early treatment. However, median shedding curves may be not accurately representing drug impact in individual volunteers. Importantly, due to public pressure clinical trials data testing oseltamivir efficacy has been recently released in the form of redacted PDF documents. We digitized and re-analyzed experimental data on influenza virus shedding in human volunteers from three previously published trials: on influenza A (1 trial) or B viruses (2 trials). Given that not all volunteers exposed to influenza viruses actually start virus shedding we found that impact of oseltamivir on the virus shedding dynamics was dependent on (i) selection of volunteers that were infected with the virus, and (ii) the detection limit in the measurement assay; both of these details were not well-articulated in the published studies. By assuming that any non-zero viral measurement is above the limit of detection we could match previously published data on median influenza A virus (flu A study) shedding but not on influenza B virus shedding (flu B study B) in human volunteers. Additional analyses confirmed that oseltamivir had an impact on the duration of shedding and overall shedding (defined as area under the curve) but this result varied by the trial. Interestingly, treatment had no impact on the rates at which shedding increased or declined with time in individual volunteers. Additional analyses showed that oseltamivir impacted the kinetics of the end of viral shedding, and in about 20-40% of volunteers that shed the virus treatment had no impact on viral shedding duration. Our results suggest an unusual impact of oseltamivir on influenza viruses shedding kinetics and caution about the use of published median data or data from a few individuals for inferences. Furthermore, we call for the need to publish raw data from critical clinical trials that can be independently analyzed.

Eosinophil Responses at the Airway Epithelial Barrier during the Early Phase of Influenza A Virus Infection in C57BL/6 Mice

Eosinophils, previously considered terminally differentiated effector cells, have multifaceted functions in tissues. We previously found that allergic mice with eosinophil-rich inflammation were protected from severe influenza and discovered specialized antiviral effector functions for eosinophils including promoting cellular immunity during influenza. In this study, we hypothesized that eosinophil responses during the early phase of influenza contribute to host protection. Using in vitro and in vivo models, we found that eosinophils were rapidly and dynamically regulated upon influenza A virus (IAV) exposure to gain migratory capabilities to traffic to lymphoid organs after pulmonary infection. Eosinophils were capable of neutralizing virus upon contact and combinations of eosinophil granule proteins reduced virus infectivity through hemagglutinin inactivation. Bi-directional crosstalk between IAV-exposed epithelial cells and eosinophils occurred after IAV infection and cross-regulation promoted barrier responses to improve antiviral defenses in airway epithelial cells. Direct interactions between eosinophils and airway epithelial cells after IAV infection prevented virus-induced cytopathology in airway epithelial cells in vitro, and eosinophil recipient IAV-infected mice also maintained normal airway epithelial cell morphology. Our data suggest that eosinophils are important in the early phase of IAV infection providing immediate protection to the epithelial barrier until adaptive immune responses are deployed during influenza.

The SKI complex is a broad-spectrum, host-directed antiviral drug target for coronaviruses, influenza, and filoviruses

The SARS-CoV-2 pandemic has made it clear that we have a desperate need for antivirals. We present work that the mammalian SKI complex is a broad-spectrum, host-directed, antiviral drug target. Yeast suppressor screening was utilized to find a functional genetic interaction between proteins from influenza A virus (IAV) and Middle East respiratory syndrome coronavirus (MERS-CoV) with eukaryotic proteins that may be potential host factors involved in replication. This screening identified the SKI complex as a potential host factor for both viruses. In mammalian systems siRNA-mediated knockdown of SKI genes inhibited replication of IAV and MERS-CoV. In silico modeling and database screening identified a binding pocket on the SKI complex and compounds predicted to bind. Experimental assays of those compounds identified three chemical structures that were antiviral against IAV and MERS-CoV along with the filoviruses Ebola and Marburg and two further coronaviruses, SARS-CoV and SARS-CoV-2. The mechanism of antiviral activity is through inhibition of viral RNA production. This work defines the mammalian SKI complex as a broad-spectrum antiviral drug target and identifies lead compounds for further development.

Model-based approach for predicting the impact of genetic modifications on product yield in biopharmaceutical manufacturing-Application to influenza vaccine production

A large group of biopharmaceuticals is produced in cell lines. The yield of such products can be increased by genetic engineering of the corresponding cell lines. The prediction of promising genetic modifications by mathematical modeling is a valuable tool to facilitate experimental screening. Besides information on the intracellular kinetics and genetic modifications the mathematical model has to account for ubiquitous cell-to-cell variability. In this contribution, we establish a novel model-based methodology for influenza vaccine production in cell lines with overexpressed genes. The manipulation of the expression level of genes coding for host cell factors relevant for virus replication is achieved by lentiviral transduction. Since lentiviral transduction causes increased cell-to-cell variability due to different copy numbers and integration sites of the gene constructs we use a population balance modeling approach to account for this heterogeneity in terms of intracellular viral components and distributed kinetic parameters. The latter are estimated from experimental data of intracellular viral RNA levels and virus titers of infection experiments using cells overexpressing a single host cell gene. For experiments with cells overexpressing multiple host cell genes, only final virus titers were measured and thus, no direct estimation of the parameter distributions was possible. Instead, we evaluate four different computational strategies to infer these from single gene parameter sets. Finally, the best computational strategy is used to predict the most promising candidates for future modifications that show the highest potential for an increased virus yield in a combinatorial study. As expected, there is a trend to higher yields the more modifications are included.

In the present work, we use a sophisticated simulation-based methodology to account for the impact of genetic modifications in producer cell lines on the yield of biomanufacturing processes. Furthermore, our approach opens the possibility to predict the most promising genetic modifications instead of identifying them in costly and time-consuming screening experiments. As an example, we apply our methodology to cell culture-based influenza vaccine production, a process that is of tremendous importance for the maintenance of public health. Here, we consider cell lines in which genes coding for one or more cellular factors are up-regulated by genetic engineering to increase the virus yield. However, the gene editing procedure increases the heterogeneity in the producer cell population because genetic modifications do not occur equally in each cell. This cell-to-cell variability is taken into account in a population balance modeling framework, thus providing a more accurate prediction of the virus yield in the heterogeneous population. Finally, we use our approach and a concise experimental data set from cell lines with one gene modification to predict the virus yield of cell lines with multiple genetic modifications. Thereby, we facilitate the experimental screening of potential candidates. We suggest that this methodology is transferable to a wide range of biomanufacturing processes and constitutes a valuable contribution to experimental design.