Rate Of Virus Production Research Articles

Of vascular defense, hemostasis, cancer, and platelet biology: an evolutionary perspective.

We have established considerable expertise in studying the role of platelets in cancer biology. From this expertise, we were keen to recognize the numerous venous-, arterial-, microvascular-, and macrovascular thrombotic events and immunologic disorders are caused by severe, acute-respiratory-syndrome coronavirus 2 (SARS-CoV-2) infections. With this offering, we explore the evolutionary connections that place platelets at the center of hemostasis, immunity, and adaptive phylogeny. Coevolutionary changes have also occurred in vertebrate viruses and their vertebrate hosts that reflect their respective evolutionary interactions. As mammals adapted from aquatic to terrestrial life and the heavy blood loss associated with placentalization-based live birth, platelets evolved phylogenetically from thrombocytes toward higher megakaryocyte-blebbing-based production rates and the lack of nuclei. With no nuclei and robust RNA synthesis, this adaptation may have influenced viral replication to become less efficient after virus particles are engulfed. Human platelets express numerous receptors that bind viral particles, which developed from archetypal origins to initiate aggregation and exocytic-release of thrombo-, immuno-, angiogenic-, growth-, and repair-stimulatory granule contents. Whether by direct, evolutionary, selective pressure, or not, these responses may help to contain virus spread, attract immune cells for eradication, and stimulate angiogenesis, growth, and wound repair after viral damage. Because mammalian and marsupial platelets became smaller and more plate-like their biophysical properties improved in function, which facilitated distribution near vessel walls in fluid-shear fields. This adaptation increased the probability that platelets could then interact with and engulf shedding virus particles. Platelets also generate circulating microvesicles that increase membrane surface-area encounters and mark viral targets. In order to match virus-production rates, billions of platelets are generated and turned over per day to continually provide active defenses and adaptation to suppress the spectrum of evolving threats like SARS-CoV-2.

Effects of UV Radiation on the Chlorophyte Micromonas polaris Host-Virus Interactions and MpoV-45T Virus Infectivity.

Polar seas are under threat of enhanced UV-radiation as well as increasing shipping activities. Considering the ecological importance of marine viruses, it is timely to study the impact of UV-AB on Arctic phytoplankton host–virus interactions and also test the efficacy of ballast water (BW) UV-C treatment on virus infectivity. This study examined the effects of: (i) ecologically relevant doses of UV-AB radiation on Micromonas polaris RCC2258 and its virus MpoV-45T, and (ii) UV-C radiation (doses 25–800 mJ cm−2) on MpoV-45T and other temperate algal viruses. Total UV-AB exposure was 6, 12, 28 and 48 h (during the light periods, over 72 h total). Strongest reduction in algal growth and photosynthetic efficiency occurred for 28 and 48 h UV-AB treatments, and consequently the virus production rates and burst sizes were reduced by more than half (compared with PAR-only controls). For the shorter UV-AB exposed cultures, negative effects by UV (especially Fv/Fm) were overcome without impacting virus proliferation. To obtain the BW desired log−4 reduction in virus infectivity, a UV-C dose of at least 400 mJ cm−2 was needed for MpoV-45T and the temperate algal viruses. This is higher than the commonly used dose of 300 mJ cm−2 in BW treatment.

Mathematical modelling of SARS-CoV-2 infection of human and animal host cells reveals differences in the infection rates and delays in viral particle production by infected cells

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV -2), a causative agent of COVID-19 disease, poses a significant threat to public health. Since its outbreak in December 2019, Wuhan, China, extensive collection of diverse data from cell culture and animal infections as well as population level data from an ongoing pandemic, has been vital in assessing strategies to battle its spread. Mathematical modelling plays a key role in quantifying determinants that drive virus infection dynamics, especially those relevant for epidemiological investigations and predictions as well as for proposing efficient mitigation strategies. We utilized a simple mathematical model to describe and explain experimental results on viral replication cycle kinetics during SARS-CoV-2 infection of animal and human derived cell lines, green monkey kidney cells, Vero-E6, and human lung epithelium cells, A549-ACE2, respectively. We conducted cell infections using two distinct initial viral concentrations and quantified viral loads over time. We then fitted the model to our experimental data and quantified the viral parameters. We showed that such cellular tropism generates significant differences in the infection rates and incubation times of SARS-CoV-2, that is, the times to the first release of newly synthesised viral progeny by SARS-CoV-2-infected cells. Specifically, the rate at which A549-ACE2 cells were infected by SARS-CoV-2 was 15 times lower than that in the case of Vero-E6 cell infection and the duration of latent phase of A549-ACE2 cells was 1.6 times longer than that of Vero-E6 cells. On the other hand, we found no statistically significant differences in other viral parameters, such as viral production rate or infected cell death rate. Since in vitro infection assays represent the first stage in the development of antiviral treatments against SARS-CoV-2, discrepancies in the viral parameter values across different cell hosts have to be identified and quantified to better target vaccine and antiviral research.

Experimental Warming Effects on Prokaryotic Growth and Viral Production in Coastal Waters of the Northwest Pacific during the Cold Season

Climate warming can directly affect biological processes in marine environments. Here, we investigated if warming (+2 °C) can change dynamics in viral and prokaryotic populations in the cold seasons in natural seawaters. We monitored the changes in viral production and prokaryotic growth rate. The prokaryotic average gross growth rates were 0.08 and 0.34 h−1 in November and 0.06 and 0.41 h−1 in December in the in situ and warming experiments, respectively. We found that warming water temperature resulted in a significant increase in prokaryotic growth rates. In warming experiments, the overall viral production rate was about 0.77–14.4 × 105 viruses mL−1 h−1, and a rough estimate of prokaryotic mortality was about 5.6–6.8 × 104 cells mL−1 h−1. Based on our estimation, burst sizes of about 21 and 14 viruses prokaryotes−1 were measured under the experimental warming period. Moreover, the results found that an increased water temperature in the subtropical western Pacific coastal waters increases prokaryotic growth rates, enhances viral production, and changes the carbon fluxes in the trophic interactions of microbes.

Modeling Within-Host Dynamics of SARS-CoV-2 Infection: A Case Study in Ferrets.

The pre-clinical development of antiviral agents involves experimental trials in animals and ferrets as an animal model for the study of SARS-CoV-2. Here, we used mathematical models and experimental data to characterize the within-host infection dynamics of SARS-CoV-2 in ferrets. We also performed a global sensitivity analysis of model parameters impacting the characteristics of the viral infection. We provide estimates of the viral dynamic parameters in ferrets, such as the infection rate, the virus production rate, the infectious virus proportion, the infected cell death rate, the virus clearance rate, as well as other related characteristics, including the basic reproduction number, pre-peak infectious viral growth rate, post-peak infectious viral decay rate, pre-peak infectious viral doubling time, post-peak infectious virus half-life, and the target cell loss in the respiratory tract. These parameters and indices are not significantly different between animals infected with viral strains isolated from the environment and isolated from human hosts, indicating a potential for transmission from fomites. While the infection period in ferrets is relatively short, the similarity observed between our results and previous results in humans supports that ferrets can be an appropriate animal model for SARS-CoV-2 dynamics-related studies, and our estimates provide helpful information for such studies.

Quantifying dose-, strain-, and tissue-specific kinetics of parainfluenza virus infection.

Human parainfluenza viruses (HPIVs) are a leading cause of acute respiratory infection hospitalization in children, yet little is known about how dose, strain, tissue tropism, and individual heterogeneity affects the processes driving growth and clearance kinetics. Longitudinal measurements are possible by using reporter Sendai viruses, the murine counterpart of HPIV 1, that express luciferase, where the insertion location yields a wild-type (rSeV-luc(M-F*)) or attenuated (rSeV-luc(P-M)) phenotype. Bioluminescence from individual animals suggests that there is a rapid increase in expression followed by a peak, biphasic clearance, and resolution. However, these kinetics vary between individuals and with dose, strain, and whether the infection was initiated in the upper and/or lower respiratory tract. To quantify the differences, we translated the bioluminescence measurements from the nasopharynx, trachea, and lung into viral loads and used a mathematical model together a nonlinear mixed effects approach to define the mechanisms distinguishing each scenario. The results confirmed a higher rate of virus production with the rSeV-luc(M-F*) virus compared to its attenuated counterpart, and suggested that low doses result in disproportionately fewer infected cells. The analyses indicated faster infectivity and infected cell clearance rates in the lung and that higher viral doses, and concomitantly higher infected cell numbers, resulted in more rapid clearance. This parameter was also highly variable amongst individuals, which was particularly evident during infection in the lung. These critical differences provide important insight into distinct HPIV dynamics, and show how bioluminescence data can be combined with quantitative analyses to dissect host-, virus-, and dose-dependent effects.

Thermodynamic effect of RNA virus infection on the human cardiovascular system.

In this work, analysis of cardiovascular system under the influence of RNA virus infection has been performed from a thermodynamic perspective. An exergetic efficiency of the system has been defined for this purpose. Results show that except for asymptomatic case, the exergetic efficiency reduces as the viral load goes up. Dynamics of viral growth along with change in efficiency is examined under different parameters such as virus production rate, infectivity rate and cell death rate. Results show that the drop in the exergetic efficiency of cardiovascular system under viral infection can be up to about 20%. Under infection, the exergy requirement of the lungs increases significantly as the work rate required by lungs increase by up to 240%.

Stability in distribution for age‐structured HIV model with delay and driven by Ornstein–Uhlenbeck process

AbstractIn this paper, we study stationarity of an age‐structured human immunodeficiency virus (HIV) model with a delay and driven by the Ornstein–Uhlenbeck process. The HIV model involves cell‐to‐virus infection and cell‐to‐cell transmission, and the two transmission modes have a short intracellular latency described by the delay. As the death rate and virus production rate of infected cells have an explicit age‐dependent feature and may be affected by environmental noises, we formulate a stochastic age‐structured HIV model to investigate the impacts of the two transmission modes on viral spread. Utilizing suitable Lyapunov functions, we first establish the existence and uniqueness of the global solution for the HIV model, and then study theth moment boundedness of the solution. Finally, the unique stationary distribution of the solution to the model is obtained. Numerical simulations are presented to demonstrate validity of our theoretical results.

First report of lytic and lysogenic viral production rates in the vicinity of fish farms (Mediterranean Sea)

The vital role of viruses in organic matter recirculation has been rarely studied in the coastal zone and never before in aquaculture areas. In this study, we examine for the first time viral lytic and lysogenic production in the vicinity of fish farms in the Mediterranean Sea. Viral production rates were determined based on the viral reduction approach at two impacted stations (including two downstream distances from cages, i.e. 50 and 100 m, and various depths in the water column, i.e. from the surface down to 50 m) and one control site. We highlight that lytic production rates were significantly higher at 50 and 100 m distances from the cages compared to the control. Maximum levels (up to 1.1 × 107 released virus-like particles mL−1 h−1) were seen at the bottom waters of the impacted sites. Interestingly, lytic production was significantly correlated with the concentrations of chlorophyll a, phaeophytin and particulate organic phosphorus, confirming the role of primary producers in fuelling prokaryotic, thus viral activity too. On the other hand, lysogenic production was not correlated with any environmental variable examined here, and did not differ between stations. We propose that viruses, acting in parallel to grazers and dilution processes, may lead to extensive “disappearance” of suspended particulate matter that is produced due to nutrient additions near fish cages. Further, we provide insights into the possibly-involved mechanisms of this viral activity.

Non-optimal codon usage preferences of coronaviruses determine their promiscuity for infecting multiple hosts.

Circulating animal coronaviruses occasionally infect humans. The SARS-CoV-2 is responsible for the current worldwide outbreak of COVID-19 that has resulted in 2112844 deaths as of late January 2021. We compared genetic code preferences in 496 viruses, including 34 coronaviruses and 242 corresponding hosts, to uncover patterns that distinguish single- and 'promiscuous' multiple-host-infecting viruses. Based on a codon usage preference score, promiscuous viruses were shown to significantly employ nonoptimal codons, namely codons that involve 'wobble' binding to anticodons, as compared to single-host viruses. The codon adaptation index (CAI) and the effective number of codons (ENC) were calculated for all viruses and hosts. Promiscuous viruses were less adapted hosts vs single-host viruses (P-value=4.392e-11). All coronaviruses exploit nonoptimal codons to infect multiple hosts. We found that nonoptimal codon preferences at the beginning of viral coding sequences enhance the translational efficiency of viral proteins within the host. Finally, coronaviruses lack endogenous RNA degradation motifs to a significant degree, thereby increasing viral mRNA burden and infection load. To conclude, we found that promiscuously infecting coronaviruses prefer nonoptimal codon usage to remove degradation motifs from their RNAs and to dramatically increase their viral RNA production rates.

Theoretical Analysis Reveals the Cost and Benefit of Proofreading in Coronavirus Genome Replication.

Severe acute respiratory syndrome coronaviruses have unusually large RNA genomes replicated by a multiprotein complex containing an RNA-dependent RNA polymerase (RdRp). Exonuclease activity enables the RdRp complex to remove wrongly incorporated bases via proofreading, a process not utilized by other RNA viruses. However, it is unclear why the RdRp complex needs proofreading and what the associated trade-offs are. Here we investigate the interplay among the accuracy, speed, and energetic cost of proofreading in the RdRp complex using a kinetic model and bioinformatics analysis. We find that proofreading nearly optimizes the rate of functional virus production. However, we find that further optimization would lead to a significant increase in the proofreading cost. Unexpected importance of the cost minimization is further supported by other global analyses. We speculate that cost optimization could help avoid cell defense responses. Thus, proofreading is essential for the production of functional viruses, but its rate is limited by energy costs.

Seasonal variations in viral distribution, dynamics, and viral-mediated host mortality in the Arabian Sea

Viruses are key players in the marine ecosystem. It is critical to study specific viral processes and their inter-relationship with various biotic and abiotic variables to quantify their impact on the marine environment. This study investigates the influence of seasonality on viral distribution and their mediated processes in the coastal region of the southeastern Arabian Sea (India) for two consecutive years (2014 and 2015). Water samples were collected from four sampling stations, S1, S2, S3 and S4, in the southeastern Arabian Sea on a monthly basis. Samples were analyzed for the variations in viral abundance, viral production, and viral-induced host mortality covering three seasons namely, pre-monsoon, monsoon and post-monsoon. Viruses and bacteria were enumerated by epifluorescence microscopy following staining with SYBR Green I and viral production was estimated using viral dilution method. Seasonal variations in viral-mediated mortality and microzooplankton grazing mortality of phytoplankton hosts were estimated by modified dilution method. The results revealed a pronounced seasonal pattern of viral abundance (0.04–4.03 × 107 viruses mL−1) and viral production (0.71 × 1010–4.94 × 1010 L−1d−1). The observed high viral-induced mortality of prokaryotes (25.85%) and phytoplankton (7.98–28.9%) during the pre-monsoon season was eventually linked to high viral abundance and production rates. Viral-mediated processes were essentially linked to host density and water temperature. Transmission electron microscopy analysis revealed that myoviruses and non-tailed viruses were dominant in the study region with rod-shaped bacteria being more susceptible to viral infection. The results indicated the prevalence of seasonally induced active viral processes in this tropical marine ecosystem which may significantly contribute to rapid recycling of nutrients in this region.

Stochastic Approach for Witnessing the Incubation Period of a Patient

The spread of HIV remains a huge investigation in this present environment. A Mathematical or Statistical model must be developed for estimating parameters related to the epidemic, the death rate of affected cells or the infectious viral production rate. Inability to carry out people evaluates their HIV status has led to widespread lack of correct and comprehensive data on HIV infection, while an individual first involved. Stochastic model measures the predicted point of threshold through discrete and continuous distribution attained by many researchers in last two decades. This paper develops a stochastic model for the time of HIV epidemic in a homosexual population. Expected time of incubation period derived through shock model approach. The fitting of information sets generated through simulation methods that the Alpha statistical distribution ought to be assumed because the epidemic distribution planned the time of stochastic model to search out HIV epidemics. To check the validity of analytical arguments and to explore the dynamics of disease above the epidemic threshold, this study concludes, the possible significance of the result is that transmit HIV in incubation stage is quicker as the intensity of the immune system is lower.

Rich dynamics of a simple delay host-pathogen model of cell-to-cell infection for plant virus

Viral dynamics within plant hosts can be important for understanding plant disease prevalence and impacts. However, few mathematical modeling efforts aim to characterize within-plant viral dynamics. In this paper, we derive a simple system of delay differential equations that describes the spread of infection throughout the plant by barley and cereal yellow dwarf viruses via the cell-to-cell mechanism. By incorporating ratio-dependent incidence function and logistic growth of the healthy cells, the model can capture a wide range of biologically relevant phenomena via the disease-free, endemic, mutual extinction steady states, and a stable periodic orbit. We show that when the basic reproduction number is less than \begin{document}$ 1 $\end{document} ( \begin{document}$ R_0 ), the disease-free steady state is asymptotically stable. When \begin{document}$ R_0>1 $\end{document} , the dynamics either converge to the endemic equilibrium or enter a periodic orbit. Using a ratio-dependent transformation, we show that if the infection rate is very high relative to the growth rate of healthy cells, then the system collapses to the mutual extinction steady state. Numerical and bifurcation simulations are provided to demonstrate our theoretical results. Finally, we carry out parameter estimation using experimental data to characterize the effects of varying nutrients on the dynamics of the system. Our parameter estimates suggest that varying the nutrient supply of nitrogen and phosphorous can alter the dynamics of the infection in plants, specifically reducing the rate of viral production and the rate of infection in certain cases.

A critical virus production rate for efficiency of oncolytic virotherapy

In a planar smoothly bounded domain$\Omega$, we consider the model for oncolytic virotherapy given by$$\left\{ \begin{array}{l} u_t = \Delta u - \nabla \cdot (u\nabla v) - uz, \\[1mm] v_t = - (u+w)v, \\[1mm] w_t = d_w \Delta w - w + uz, \\[1mm] z_t = d_z \Delta z - z - uz + \beta w, \end{array} \right.$$with positive parameters$ D_w $,$ D_z $and$\beta$. It is firstly shown that whenever$\beta \lt 1$, for any choice of$M \gt 0$, one can find initial data such that the solution of an associated no-flux initial-boundary value problem, well known to exist globally actually for any choice of$\beta \gt 0$, satisfies$$u\ge M \qquad \mbox{in } \Omega\times (0,\infty).$$If$\beta \gt 1$, however, then for arbitrary initial data the corresponding is seen to have the property that$$\liminf_{t\to\infty} \inf_{x\in\Omega} u(x,t)\le \frac{1}{\beta-1}.$$This may be interpreted as indicating that$\beta$plays the role of a critical virus replication rate with regard to efficiency of the considered virotherapy, with corresponding threshold value given by$\beta = 1$.

A mathematical model describing the localization and spread of influenza A virus infection within the human respiratory tract.

Within the human respiratory tract (HRT), virus diffuses through the periciliary fluid (PCF) bathing the epithelium. But virus also undergoes advection: as the mucus layer sitting atop the PCF is pushed along by the ciliated cell's beating cilia, the PCF and its virus content are also pushed along, upwards towards the nose and mouth. While many mathematical models (MMs) have described the course of influenza A virus (IAV) infections in vivo, none have considered the impact of both diffusion and advection on the kinetics and localization of the infection. The MM herein represents the HRT as a one-dimensional track extending from the nose down towards the lower HRT, wherein stationary cells interact with IAV which moves within (diffusion) and along with (advection) the PCF. Diffusion was found to be negligible in the presence of advection which effectively sweeps away IAV, preventing infection from disseminating below the depth at which virus first deposits. Higher virus production rates (10-fold) are required at higher advection speeds (40 μm/s) to maintain equivalent infection severity and timing. Because virus is entrained upwards, upper parts of the HRT see more virus than lower parts. As such, infection peaks and resolves faster in the upper than in the lower HRT, making it appear as though infection progresses from the upper towards the lower HRT, as reported in mice. When the spatial MM is expanded to include cellular regeneration and an immune response, it reproduces tissue damage levels reported in patients. It also captures the kinetics of seasonal and avian IAV infections, via parameter changes consistent with reported differences between these strains, enabling comparison of their treatment with antivirals. This new MM offers a convenient and unique platform from which to study the localization and spread of respiratory viral infections within the HRT.

A critical virus production rate for blow-up suppression in a haptotaxis model for oncolytic virotherapy

In a two-dimensional smoothly bounded domain, we consider the oncolytic virotherapy model given by ut=Δu−∇⋅(u∇v)−ρuz,vt=−(u+w)v,wt=DwΔw−w+uz,zt=DzΔz−z−uz+βw,with parameters Dw>0,Dz>0,β>0 and ρ≥0.According to previous findings, a corresponding Neumann initial–boundary value problem possesses a globally defined classical solution whenever the prescribed initial data are suitably regular and satisfy appropriate positivity assumptions. In the present study it is shown that if β<1, then any such solution is uniformly bounded. This complements a recent result which in the case ρ=0 has asserted the existence of unbounded solutions whenever β>1.

Drivers of Bacterial α- and β-Diversity Patterns and Functioning in Subsurface Hadal Sediments.

Oceanic trenches at hadal (>6,000 m) depths are hot spots of organic matter deposition and mineralization and can host abundant and active bacterial assemblages. However, the factors able to shape their biodiversity and functioning remain largely unexplored, especially in subsurface sediments. Here, we investigated the patterns and drivers of benthic bacterial α- and β-diversity (i.e., OTU richness and turnover diversity) along the vertical profile down to 1.5 m sediment depth in the Izu-Bonin Trench (at ~10,000 m water depth). The protease and glucosidase enzymatic activity rates were also determined, as a proxy of organic matter degradation potential in the different sediment layers. Molecular fingerprinting based on automated ribosomal intergenic spacer analysis (ARISA) indicated that the α-diversity of bacterial assemblages remained high throughout the vertical profile and that the turnover (β-) diversity among sediment horizons reached values up to 90% of dissimilarity. Multivariate distance-based linear modeling (DISTLM) pointed out that the diversity and functioning of the hadal bacterial assemblages were influenced by the variability of environmental conditions (including the availability of organic resources and electron donors/acceptors) and of viral production rates along the sediment vertical profile. Based on our results, we can argue that the heterogeneity of physical-chemical features of the hadal sediments of the Izu-Bonin Trench contribute to increase the niches availability for different bacterial taxa, while viruses contribute to maintain high levels of bacterial turnover diversity and to enhance organic matter cycling in these extremely remote and isolated ecosystems.

Revealing uninfected and infected target cell dynamics from peripheral blood data in highly and less pathogenic simian/human immunodeficiency virus infected Rhesus macaque

Since chimeric simian and human immunodeficiency viruses (SHIVs) used here, that is, SHIV-#64 and -KS661 utilize both CCR5 and CXCR4 chemokine receptors, they have broad target cell properties. A highly pathogenic SHIV strain, SHIV-KS661, causes an infection that systemically depletes the CD4+ T cells of Rhesus macaques, while a less pathogenic strain, SHIV-#64, does not cause severe symptoms in the macaques. In our previous studies, we established in vitro quantification system for virus infection dynamics, and concluded that SHIV-KS661 effectively produces infectious virions compared with SHIV-#64 in the HSC-F cell culture. However, in vivo dynamics of SHIV infection have not been well understood. To quantify SHIV-#64 and -KS661 infection dynamics in Rhesus macaques, we developed a novel approach and analyzed total CD4+ T cells and viral load in peripheral blood, and reproduced the expected dynamics for the uninfected and infected CD4+ T cells in silico. Using our previous cell culture experimental datasets, we revealed that an infection rate constant is different between SHIV-#64 and -KS661, but the viral production rate and the death rate are similar for the both strains. Thus, here, we assumed these relations in our in vivo data and carried out the data fitting. We performed Bayesian estimation for the whole dataset using MCMC sampling, and simultaneously fitted our novel model to total CD4+ T cells and viral load of SHIV-#64 and -KS661 infection. Our analyses explained that the Malthusian parameter (i.e., fitness of virus infection) and the basic reproduction number (i.e., potential of virus infection) for SHIV-KS661 are significantly higher than those of SHIV-#64. In addition, we demonstrated that the number of uninfected CD4+ T cells in SHIV-KS661 infected Rhesus macaques decreases to the significantly lower value than that before the inoculation several days earlier compared with SHIV-#64 infection. Taken together, the differences between SHIV-#64 and -KS661 infection before the peak viral load might determine the subsequent destiny, that is, whether the systemic CD4+ T cell depletion occurs or the host immune response develop.

Within-host infectious disease models accommodating cellular coinfection, with an application to influenza.

Within-host models are useful tools for understanding the processes regulating viral load dynamics. While existing models have considered a wide range of within-host processes, at their core these models have shown remarkable structural similarity. Specifically, the structure of these models generally consider target cells to be either uninfected or infected, with the possibility of accommodating further resolution (e.g. cells that are in an eclipse phase). Recent findings, however, indicate that cellular coinfection is the norm rather than the exception for many viral infectious diseases, and that cells with high multiplicity of infection are present over at least some duration of an infection. The reality of these cellular coinfection dynamics is not accommodated in current within-host models although it may be critical for understanding within-host dynamics. This is particularly the case if multiplicity of infection impacts infected cell phenotypes such as their death rate and their viral production rates. Here, we present a new class of within-host disease models that allow for cellular coinfection in a scalable manner by retaining the low-dimensionality that is a desirable feature of many current within-host models. The models we propose adopt the general structure of epidemiological ‘macroparasite’ models that allow hosts to be variably infected by parasites such as nematodes and host phenotypes to flexibly depend on parasite burden. Specifically, our within-host models consider target cells as ‘hosts’ and viral particles as ‘macroparasites’, and allow viral output and infected cell lifespans, among other phenotypes, to depend on a cell’s multiplicity of infection. We show with an application to influenza that these models can be statistically fit to viral load and other within-host data, and demonstrate using model selection approaches that they have the ability to outperform traditional within-host viral dynamic models. Important in vivo quantities such as the mean multiplicity of cellular infection and time-evolving reassortant frequencies can also be quantified in a straightforward manner once these macroparasite models have been parameterized. The within-host model structure we develop here provides a mathematical way forward to address questions related to the roles of cellular coinfection, collective viral interactions, and viral complementation in within-host viral dynamics and evolution.