Quantitative Analysis of Hepatitis C NS5A Viral Protein Dynamics on the ER Surface

Intracellular "In Silico Microscopes"-Comprehensive 3D Spatio-Temporal Virus Replication Model Simulations

Despite their small and simple structure compared with their hosts, virus particles can cause severe harm and even mortality in highly evolved species such as humans. A comprehensive quantitative biophysical understanding of intracellular virus replication mechanisms could aid in preparing for future virus pandemics. By elucidating the relationship between the form and function of intracellular structures from the host cell and viral components, it is possible to identify possible targets for direct antiviral agents and potent vaccines. Biophysical investigations into the spatio-temporal dynamics of intracellular virus replication have thus far been limited. This study introduces a framework to enable simulations of these dynamics using partial differential equation (PDE) models, which are evaluated using advanced numerical mathematical methods on leading supercomputers. In particular, this study presents a model of the replication cycle of a specific RNA virus, the hepatitis C virus. The diffusion-reaction model mimics the interplay of the major components of the viral replication cycle, including non structural viral proteins, viral genomic RNA, and a generic host factor. Technically, surface partial differential equations (sufPDEs) are coupled on the 3D embedded 2D endoplasmic reticulum manifold with partial differential equations (PDEs) in the 3D membranous web and cytosol volume. The membranous web serves as a viral replication factory and is formed on the endoplasmic reticulum after infection and in the presence of nonstructural proteins. The coupled sufPDE/PDE model was evaluated using realistic cell geometries based on experimental data. The simulations incorporate the effects of non structural viral proteins, which are restricted to the endoplasmic reticulum surface, with effects appearing in the volume, such as host factor supply from the cytosol and membranous web dynamics. Because the spatial diffusion properties of genomic viral RNA are not yet fully understood, the model allows for viral RNA movement on the endoplasmic reticulum as well as within the cytosol. Visualizing the simulated intracellular viral replication dynamics provides insights similar to those obtained by microscopy, complementing data from in vitro/in vivo viral replication experiments. The output data demonstrate quantitative consistence with the experimental findings, prompting further advanced experimental studies to validate the model and refine our quantitative biophysical understanding.

Efficient Estimates of Surface Diffusion Parameters for Spatio-Temporally Resolved Virus Replication Dynamics

Advanced methods of treatment are needed to fight the threats of virus-transmitted diseases and pandemics. Often, they are based on an improved biophysical understanding of virus replication strategies and processes in their host cells. For instance, an essential component of the replication of the hepatitis C virus (HCV) proceeds under the influence of nonstructural HCV proteins (NSPs) that are anchored to the endoplasmatic reticulum (ER), such as the NS5A protein. The diffusion of NSPs has been studied by in vitro fluorescence recovery after photobleaching (FRAP) experiments. The diffusive evolution of the concentration field of NSPs on the ER can be described by means of surface partial differential equations (sufPDEs). Previous work estimated the diffusion coefficient of the NS5A protein by minimizing the discrepancy between an extended set of sufPDE simulations and experimental FRAP time-series data. Here, we provide a scaling analysis of the sufPDEs that describe the diffusive evolution of the concentration field of NSPs on the ER. This analysis provides an estimate of the diffusion coefficient that is based only on the ratio of the membrane surface area in the FRAP region to its contour length. The quality of this estimate is explored by a comparison to numerical solutions of the sufPDE for a flat geometry and for ten different 3D embedded 2D ER grids that are derived from fluorescence z-stack data of the ER. Finally, we apply the new data analysis to the experimental FRAP time-series data analyzed in our previous paper, and we discuss the opportunities of the new approach.

Influence of T-Bar on Calcium Concentration Impacting Release Probability

The relation of form and function, namely the impact of the synaptic anatomy on calcium dynamics in the presynaptic bouton, is a major challenge of present (computational) neuroscience at a cellular level. The Drosophila larval neuromuscular junction (NMJ) is a simple model system, which allows studying basic effects in a rather simple way. This synapse harbors several special structures. In particular, in opposite to standard vertebrate synapses, the presynaptic boutons are rather large, and they have several presynaptic zones. In these zones, different types of anatomical structures are present. Some of the zones bear a so-called T-bar, a particular anatomical structure. The geometric form of the T-bar resembles the shape of the letter “T” or a table with one leg. When an action potential arises, calcium influx is triggered. The probability of vesicle docking and neurotransmitter release is superlinearly proportional to the concentration of calcium close to the vesicular release site. It is tempting to assume that the T-bar causes some sort of calcium accumulation and hence triggers a higher release probability and thus enhances neurotransmitter exocytosis. In order to study this influence in a quantitative manner, we constructed a typical T-bar geometry and compared the calcium concentration close to the active zones (AZs). We compared the case of synapses with and without T-bars. Indeed, we found a substantial influence of the T-bar structure on the presynaptic calcium concentrations close to the AZs, indicating that this anatomical structure increases vesicle release probability. Therefore, our study reveals how the T-bar zone implies a strong relation between form and function. Our study answers the question of experimental studies (namely “Wichmann and Sigrist, Journal of neurogenetics 2010”) concerning the sense of the anatomical structure of the T-bar.

Lyapunov function and global asymptotic stability for a new multiscale viral dynamics model incorporating the immune system response: Implemented upon HCV

In this paper, a new mathematical model is formulated that describes the interaction between uninfected cells, infected cells, viruses, intracellular viral RNA, Cytotoxic T-lymphocytes (CTLs), and antibodies. Hence, the model contains certain biological relations that are thought to be key factors driving this interaction which allow us to obtain precise logical conclusions. Therefore, it improves our perception, that would otherwise not be possible, to comprehend the pathogenesis, to interpret clinical data, to control treatment, and to suggest new relations. This model can be used to study viral dynamics in patients for a wide range of infectious diseases like HIV, HPV, HBV, HCV, and Covid-19. Though, analysis of a new multiscale HCV model incorporating the immune system response is considered in detail, the analysis and results can be applied for all other viruses. The model utilizes a transformed multiscale model in the form of ordinary differential equations (ODE) and incorporates into it the interaction of the immune system. The role of CTLs and the role of antibody responses are investigated. The positivity of the solutions is proven, the basic reproduction number is obtained, and the equilibrium points are specified. The stability at the equilibrium points is analyzed based on the Lyapunov invariance principle. By using appropriate Lyapunov functions, the uninfected equilibrium point is proven to be globally asymptotically stable when the reproduction number is less than one and unstable otherwise. Global stability of the infected equilibrium points is considered, and it has been found that each equilibrium point has a specific domain of stability. Stability regions could be overlapped and a bistable equilibria could be found, which means the coexistence of two stable equilibrium points. Hence, the solution converges to one of them depending on the initial conditions.

Advanced Hepatitis C Virus Replication PDE Models within a Realistic Intracellular Geometric Environment

The hepatitis C virus (HCV) RNA replication cycle is a dynamic intracellular process occurring in three-dimensional space (3D), which is difficult both to capture experimentally and to visualize conceptually. HCV-generated replication factories are housed within virus-induced intracellular structures termed membranous webs (MW), which are derived from the Endoplasmatic Reticulum (ER). Recently, we published 3D spatiotemporal resolved diffusion⁻reaction models of the HCV RNA replication cycle by means of surface partial differential equation (sPDE) descriptions. We distinguished between the basic components of the HCV RNA replication cycle, namely HCV RNA, non-structural viral proteins (NSPs), and a host factor. In particular, we evaluated the sPDE models upon realistic reconstructed intracellular compartments (ER/MW). In this paper, we propose a significant extension of the model based upon two additional parameters: different aggregate states of HCV RNA and NSPs, and population dynamics inspired diffusion and reaction coefficients instead of multilinear ones. The combination of both aspects enables realistic modeling of viral replication at all scales. Specifically, we describe a replication complex state consisting of HCV RNA together with a defined amount of NSPs. As a result of the combination of spatial resolution and different aggregate states, the new model mimics a cis requirement for HCV RNA replication. We used heuristic parameters for our simulations, which were run only on a subsection of the ER. Nevertheless, this was sufficient to allow the fitting of core aspects of virus reproduction, at least qualitatively. Our findings should help stimulate new model approaches and experimental directions for virology.