Advanced methods of treatment are needed to fight the threats of virus-transmitted deseases 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 non-structural 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 (sufPDE). Knodel et al, Viruses 2018, 10, 28 estimated the diffusion coefficient of the NS5a protein by minimizing the discrepancy of 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 and its contour length. The quality of this estimate is explored by 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.

Despite being small and simple structured in comparison to their victims, virus particles have the potential to harm severly and even kill highly developed species such as humans. To face upcoming virus pandemics, detailed quantitative biophysical understanding of intracellular virus replication mechanisms is inevitable. Unveiling the relationship of form and function might allow to determine putative attack points relevant for the systematic development of direct antiviral agents (DAA) and potent vaccines. Biophysical investigations of spatio-temporal dynamics of intracellular virus replication so far are rare. We are developing a framework to allow for fully spatio-temporal resolved virus replication dynamics simulations based on partial differential equations models (PDE) and evaluated with advanced numerical mathematical methods on leading supercomputers. This study presents an advanced highly nonlinear model of the genome replication cycle of a specific RNA virus, the Hepatitis C virus (HCV). The diffusion-reaction model mimics the interplay of the major components of the viral RNA (vRNA) cycle, namely non structural viral proteins (NSP), vRNA and a generic host factor. Technically, we couple surface PDEs (sufPDEs) on the 3D embedded 2D Endoplasmatic Reticulum (ER) manifold with PDEs in the 3D membranous web (MW) and cytosol volume. (The MWs are the replication factories growing on the ER induced by NSPs.) The sufPDE/PDE model is evaluated at realistic reconstructed cell geometries which are based on experimental data. The simulations couple the effects of NSPs which are restricted to the ER surface with effects appearing in the volume. The volume effects include the host factor supply from the cytosol and the MW dynamics. Special emphasis is put to the exchange of components between ER surface, MWs and cytosol volume. As the vRNA spatial properties are not fully understood so far in experiment, the model allows for vRNA both restricted to the ER and moving in the cytosol. The visualization of the simulation resembles a look into some sort of fully 3D resolved “in silico microscopes” to mirror and complement in vitro / in vivo experiments for the intracellular vRNA cycle dynamics. The output data are quantitatively consistent with experimental data and provoke advanced experimental studies to validate the model.

We consider a class of density-driven flow problems. We are particularly interested in the problem of the salinization of coastal aquifers. We consider the Henry saltwater intrusion problem with uncertain porosity, permeability, and recharge parameters as a test case. The reason for the presence of uncertainties is the lack of knowledge, inaccurate measurements, and inability to measure parameters at each spatial or time location. This problem is nonlinear and time-dependent. The solution is the salt mass fraction, which is uncertain and changes in time. Uncertainties in porosity, permeability, recharge, and mass fraction are modeled using random fields. This work investigates the applicability of the well-known multilevel Monte Carlo (MLMC) method for such problems. The MLMC method can reduce the total computational and storage costs. Moreover, the MLMC method runs multiple scenarios on different spatial and time meshes and then estimates the mean value of the mass fraction. The parallelization is performed in both the physical space and stochastic space. To solve every deterministic scenario, we run the parallel multigrid solver ug4 in a black-box fashion. We use the solution obtained from the quasi-Monte Carlo method as a reference solution.