Chromatin organization regulates viral egress dynamics.

Various types of DNA viruses are known to elicit the formation of a large nuclear viral replication compartment and marginalization of the cell chromatin. We used three-dimensional soft x-ray tomography, confocal and electron microscopy, combined with numerical modelling of capsid diffusion to analyse the molecular organization of chromatin in herpes simplex virus 1 infection and its effect on the transport of progeny viral capsids to the nuclear envelope. Our data showed that the formation of the viral replication compartment at late infection resulted in the enrichment of heterochromatin in the nuclear periphery accompanied by the compaction of chromatin. Random walk modelling of herpes simplex virus 1-sized particles in a three-dimensional soft x-ray tomography reconstruction of an infected cell nucleus demonstrated that the peripheral, compacted chromatin restricts viral capsid diffusion, but due to interchromatin channels capsids are able to reach the nuclear envelope, the site of their nuclear egress.

Diffusion through thin membranes: Modeling across scales.

From macroscopic to microscopic scales it is demonstrated that diffusion through membranes can be modeled using specific boundary conditions across them. The membranes are here considered thin in comparison to the overall size of the system. In a macroscopic scale the membrane is introduced as a transmission boundary condition, which enables an effective modeling of systems that involve multiple scales. In a mesoscopic scale, a numerical lattice-Boltzmann scheme with a partial-bounceback condition at the membrane is proposed and analyzed. It is shown that this mesoscopic approach provides a consistent approximation of the transmission boundary condition. Furthermore, analysis of the mesoscopic scheme gives rise to an expression for the permeability of a thin membrane as a function of a mesoscopic transmission parameter. In a microscopic model, the mean waiting time for a passage of a particle through the membrane is in accordance with this permeability. Numerical results computed with the mesoscopic scheme are then compared successfully with analytical solutions derived in a macroscopic scale, and the membrane model introduced here is used to simulate diffusive transport between the cell nucleus and cytoplasm through the nuclear envelope in a realistic cell model based on fluorescence microscopy data. By comparing the simulated fluorophore transport to the experimental one, we determine the permeability of the nuclear envelope of HeLa cells to enhanced yellow fluorescent protein.

Investigation of an entropic stabilizer for the lattice-Boltzmann method.

The lattice-Boltzmann (LB) method is commonly used for the simulation of fluid flows at the hydrodynamic level of description. Due to its kinetic theory origins, the standard LB schemes carry more degrees of freedom than strictly needed, e.g., for the approximation of solutions to the Navier-stokes equation. In particular, there is freedom in the details of the so-called collision operator. This aspect was recently utilized when an entropic stabilizer, based on the principle of maximizing local entropy, was proposed for the LB method [I. V. Karlin, F. Bösch, and S. S. Chikatamarla, Phys. Rev. E 90, 031302(R) (2014)]. The proposed stabilizer can be considered as an add-on or extension to basic LB schemes. Here the entropic stabilizer is investigated numerically using the perturbed double periodic shear layer flow as a benchmark case. The investigation is carried out by comparing numerical results obtained with six distinct LB schemes. The main observation is that the unbounded, and not explicitly controllable, relaxation time for the higher-order moments will directly influence the leading-order error terms. As a consequence, the order of accuracy and, in general, the numerical behavior of LB schemes are substantially altered. Hence, in addition to systematic numerical validation, more detailed theoretical analysis of the entropic stabilizer is still required in order to properly understand its properties.

High-accuracy approximation of high-rank derivatives: isotropic finite differences based on lattice-Boltzmann stencils.

We propose isotropic finite differences for high-accuracy approximation of high-rank derivatives. These finite differences are based on direct application of lattice-Boltzmann stencils. The presented finite-difference expressions are valid in any dimension, particularly in two and three dimensions, and any lattice-Boltzmann stencil isotropic enough can be utilized. A theoretical basis for the proposed utilization of lattice-Boltzmann stencils in the approximation of high-rank derivatives is established. In particular, the isotropy and accuracy properties of the proposed approximations are derived directly from this basis. Furthermore, in this formal development, we extend the theory of Hermite polynomial tensors in the case of discrete spaces and present expressions for the discrete inner products between monomials and Hermite polynomial tensors. In addition, we prove an equivalency between two approaches for constructing lattice-Boltzmann stencils. For the numerical verification of the presented finite differences, we introduce 5th-, 6th-, and 8th-order two-dimensional lattice-Boltzmann stencils.