No Arabic abstract
Influenza viruses enter a cell via endocytosis after binding to the surface. During the endosomal journey, acidification triggers a conformational change of the virus spike protein hemagglutinin (HA) that results in escape of the viral genome from the endosome to the cytoplasm. A quantitative understanding of the processes involved in HA mediated fusion with the endosome is still missing. We develop here a stochastic model to estimate the change of conformation of HAs inside the endosome nanodomain. Using a Markov-jump process to model the arrival of protons to HA binding sites, we compute the kinetics of their accumulation and the mean first time for HAs to be activated. This analysis reveals that HA proton binding sites possess a high chemical barrier, ensuring a stability of the spike protein at sub-acidic pH. Finally, we predict that activating more than 3 adjacent HAs is necessary to prevent a premature fusion.
The evolutionary dynamics of human Influenza A virus presents a challenging theoretical problem. An extremely high mutation rate allows the virus to escape, at each epidemic season, the host immune protection elicited by previous infections. At the same time, at each given epidemic season a single quasi-species, that is a set of closely related strains, is observed. A non-trivial relation between the genetic (i.e., at the sequence level) and the antigenic (i.e., related to the host immune response) distances can shed light into this puzzle. In this paper we introduce a model in which, in accordance with experimental observations, a simple interaction rule based on spatial correlations among point mutations dynamically defines an immunity space in the space of sequences. We investigate the static and dynamic structure of this space and we discuss how it affects the dynamics of the virus-host interaction. Interestingly we observe a staggered time structure in the virus evolution as in the real Influenza evolutionary dynamics.
In this paper, we carry out a computational study using the spectral decomposition of the fluctuations of a two-pathogen epidemic model around its deterministic attractor, i.e., steady state or limit cycle, to examine the role of partial vaccination and between-host pathogen interaction on early pathogen replacement during seasonal epidemics of influenza and respiratory syncytial virus.
We study the dynamics of membrane vesicle motor transport into dendritic spines, which are bulbous intracellular compartments in neurons that play a key role in transmitting signals between neurons. We consider the stochastic analog of the vesicle transport model in [Park and Fai, The Dynamics of Vesicles Driven Into Closed Constrictions by Molecular Motors. Bull. Math. Biol. 82, 141 (2020)]. The stochastic version, which may be considered as an agent-based model, relies mostly on the action of individual myosin motors to produce vesicle motion. To aid in our analysis, we coarse-grain this agent-based model using a master equation combined with a partial differential equation describing the probability of local motor positions. We confirm through convergence studies that the coarse-graining captures the essential features of bistability in velocity (observed in experiments) and waiting-time distributions to switch between steady-state velocities. Interestingly, these results allow us to reformulate the translocation problem in terms of the mean first passage time for a run-and-tumble particle moving on a finite domain with absorbing boundaries at the two ends. We conclude by presenting numerical and analytical calculations of vesicle translocation.
We consider self-assembly of proteins into a virus capsid by the methods of molecular dynamics. The capsid corresponds either to SPMV or CCMV and is studied with and without the RNA molecule inside. The proteins are flexible and described by the structure-based coarse-grained model augmented by electrostatic interactions. Previous studies of the capsid self-assembly involved solid objects of a supramolecular scale, e.g. corresponding to capsomeres, with engineered couplings and stochastic movements. In our approach, a single capsid is dissociated by an application of a high temperature for a variable period and then the system is cooled down to allow for self-assembly. The restoration of the capsid proceeds to various extent, depending on the nature of the dissociated state, but is rarely complete because some proteins depart too far unless the process takes place in a confined space.
This paper describes a mathematical model for the spread of a virus through an isolated population of a given size. The model uses three, color-coded components, called molecules (red for infected and still contagious; green for infected, but no longer contagious; and blue for uninfected). In retrospect, the model turns out to be a digital analogue for the well-known SIR model of Kermac and McKendrick (1927). In our RGB model, the number of accumulated infections goes through three phases, beginning at a very low level, then changing to a transition ramp of rapid growth, and ending in a plateau of final values. Consequently, the differential change or growth rate begins at 0, rises to a peak corresponding to the maximum slope of the transition ramp, and then falls back to 0. The properties of these time variations, including the slope, duration, and height of the transition ramp, and the width and height of the infection rate, depend on a single parameter - the time that a red molecule is contagious divided by the average time between collisions of the molecules. Various temporal milestones, including the starting time of the transition ramp, the time that the accumulating number of infections obtains its maximum slope, and the location of the peak of the infection rate depend on the size of the population in addition to the contagious lifetime ratio. Explicit formulas for these quantities are derived and summarized. Finally, Appendix E has been added to describe the effect of vaccinations.