Host-virus evolutionary dynamics with specialist and generalist infection strategies: Bifurcations, bistability, and chaos.

In this work, we have investigated the evolutionary dynamics of a generalist pathogen, e.g., a virus population, that evolves toward specialization in an environment with multiple host types. We have particularly explored under which conditions generalist viral strains may rise in frequency and coexist with specialist strains or even dominate the population. By means of a nonlinear mathematical model and bifurcation analysis, we have determined the theoretical conditions for stability of nine identified equilibria and provided biological interpretation in terms of the infection rates for the viral specialist and generalist strains. By means of a stability diagram, we identified stable fixed points and stable periodic orbits, as well as regions of bistability. For arbitrary biologically feasible initial population sizes, the probability of evolving toward stable solutions is obtained for each point of the analyzed parameter space. This probability map shows combinations of infection rates of the generalist and specialist strains that might lead to equal chances for each type becoming the dominant strategy. Furthermore, we have identified infection rates for which the model predicts the onset of chaotic dynamics. Several degenerate Bogdanov-Takens and zero-Hopf bifurcations are detected along with generalized Hopf and zero-Hopf bifurcations. This manuscript provides additional insights into the dynamical complexity of host-pathogen evolution toward different infection strategies.

Monomer diffusion into static and evolving polymer networks during frontal photopolymerisation.

Frontal photopolymerisation (FPP) is a directional solidification process that converts monomer-rich liquid into crosslinked polymer solid by light exposure and finds applications ranging from lithography to 3D printing. Inherent to this process is the creation of an evolving polymer network that is exposed to a monomer bath. A combined theoretical and experimental investigation is performed to determine the conditions under which monomer from this bath can diffuse into the propagating polymer network and cause it to swell. First, the growth and swelling processes are decoupled by immersing pre-made polymer networks into monomer baths held at various temperatures. The experimental measurements of the network thickness are found to be in good agreement with theoretical predictions obtained from a nonlinear poroelastic model. FPP propagation experiments are then carried out under conditions that lead to swelling. Unexpectedly, for a fixed exposure time, swelling is found to increase with incident light intensity. The experimental data is well described by a novel FPP model accounting for mass transport and the mechanical response of the polymer network, providing key insights into how monomer diffusion affects the conversion profile of the polymer solid and the stresses that are generated during its growth. The predictive capability of the model will enable the fabrication of gradient materials with tuned mechanical properties and controlled stress development.

Surface waves on a soft viscoelastic layer produced by an oscillating microbubble.

Ultrasound-driven bubbles can cause significant deformation of soft viscoelastic layers, for instance in surface cleaning and biomedical applications. The effect of the viscoelastic properties of a boundary on the bubble-boundary interaction has been explored only qualitatively, and remains poorly understood. We investigate the dynamic deformation of a viscoelastic layer induced by the volumetric oscillations of an ultrasound-driven microbubble. High-speed video microscopy is used to observe the deformation produced by a bubble oscillating at 17-20 kHz in contact with the surface of a hydrogel. The localised oscillating pressure applied by the bubble generates surface elastic (Rayleigh) waves on the gel, characterised by elliptical particle trajectories. The tilt angle of the elliptical trajectories varies with increasing distance from the bubble. Unexpectedly, the direction of rotation of the surface elements on the elliptical trajectories shifts from prograde to retrograde at a distance from the bubble that depends on the viscoelastic properties of the gel. To explain these behaviours, we develop a simple three-dimensional model for the deformation of a viscoelastic solid by a localised oscillating force. By using as input for the model the values of the shear modulus obtained from the propagation velocity of the Rayleigh waves, we find good qualitative agreement with the experimental observations.

Exploiting in silico modelling to enhance translation of liver cell therapies from bench to bedside.

Cell therapies are emerging as promising treatments for a range of liver diseases but translational bottlenecks still remain including: securing and assessing the safe and effective delivery of cells to the disease site; ensuring successful cell engraftment and function; and preventing immunogenic responses. Here we highlight three therapies, each utilising a different cell type, at different stages in their clinical translation journey: transplantation of multipotent mesenchymal stromal/signalling cells, hepatocytes and macrophages. To overcome bottlenecks impeding clinical progression, we advocate for wider use of mechanistic in silico modelling approaches. We discuss how in silico approaches, alongside complementary experimental approaches, can enhance our understanding of the mechanisms underlying successful cell delivery and engraftment. Furthermore, such combined theoretical-experimental approaches can be exploited to develop novel therapies, address safety and efficacy challenges, bridge the gap between in vitro and in vivo model systems, and compensate for the inherent differences between animal model systems and humans. We also highlight how in silico model development can result in fewer and more targeted in vivo experiments, thereby reducing preclinical costs and experimental animal numbers and potentially accelerating translation to the clinic. The development of biologically-accurate in silico models that capture the mechanisms underpinning the behaviour of these complex systems must be reinforced by quantitative methods to assess cell survival post-transplant, and we argue that non-invasive in vivo imaging strategies should be routinely integrated into transplant studies.

Drying-induced stresses in poroelastic drops on rigid substrates.

We develop a theory for drying-induced stresses in sessile, poroelastic drops undergoing evaporation on rigid surfaces. Using a lubrication-like approximation, the governing equations of three-dimensional nonlinear poroelasticity are reduced to a single thin-film equation for the drop thickness. We find that thin drops experience compressive elastic stresses but the total in-plane stresses are tensile. The mechanical response of the drop is dictated by the initial profile of the solid skeleton, which controls the in-plane deformation, the dominant components of elastic stress, and sets a limit on the depth of delamination that can potentially occur. Our theory suggests that the alignment of desiccation fractures in colloidal drops is selected by the shape of the drop at the point of gelation. We propose that the emergence of three distinct fracture patterns in dried blood drops is a consequence of a nonmonotonic drop profile at gelation. We also show that depletion fronts, which separate wet and dry solid, can invade the drop from the contact line and localize the generation of mechanical stress during drying. Finally, the finite element method is used to explore the stress profiles in drops with large contact angles.

Phase separation in swelling and deswelling hydrogels with a free boundary.

We present a full kinetic model of a hydrogel that undergoes phase separation during swelling and deswelling. The model accounts for the interfacial energy of coexisting phases, finite strain of the polymer network, and solvent transport across free boundaries. For the geometry of an initially dry layer bonded to a rigid substrate, the model predicts that forcing solvent into the gel at a fixed rate can induce a volume phase transition, which gives rise to coexisting phases with different degrees of swelling, in systems where this cannot occur in the free-swelling case. While a nonzero shear modulus assists in the propagation of the transition front separating these phases in the driven-swelling case, increasing it beyond a critical threshold suppresses its formation. Quenching a swollen hydrogel induces spinodal decomposition, which produces several highly localized, highly swollen phases which coarsen and are then ejected from free boundary. The wealth of dynamic scenarios of this system is discussed using phase-plane analysis and numerical solutions in a one-dimensional setting.

Theoretical conditions for the coexistence of viral strains with differences in phenotypic traits: a bifurcation analysis.

We investigate the dynamics of a wild-type viral strain which generates mutant strains differing in phenotypic properties for infectivity, virulence and mutation rates. We study, by means of a mathematical model and bifurcation analysis, conditions under which the wild-type and mutant viruses, which compete for the same host cells, can coexist. The coexistence conditions are formulated in terms of the basic reproductive numbers of the strains, a maximum value of the mutation rate and the virulence of the pathogens. The analysis reveals that parameter space can be divided into five regions, each with distinct dynamics, that are organized around degenerate Bogdanov-Takens and zero-Hopf bifurcations, the latter of which gives rise to a curve of transcritical bifurcations of periodic orbits. These results provide new insights into the conditions by which viral populations may contain multiple coexisting strains in a stable manner.

A minimal model for solvent evaporation and absorption in thin films.

We present a minimal model of solvent evaporation and absorption in thin films consisting of a volatile solvent and non-volatile solutes. An asymptotic analysis yields expressions that facilitate the extraction of physically significant model parameters from experimental data, namely the mass transfer coefficient and composition-dependent diffusivity. The model can be used to predict the dynamics of drying and film formation, as well as sorption/desorption, over a wide range of experimental conditions. A state diagram is used to understand the experimental conditions that lead to the formation of a solute-rich layer, or "skin", at the evaporating surface during drying. In the case of solvent absorption, the model captures the existence of a saturation front that propagates from the film surface towards the substrate. The theoretical results are found to be in excellent agreement with data produced from dynamic vapour sorption experiments of ternary mixtures comprising an aluminium salt, glycerol, and water. Moreover, the model should be generally applicable to a variety of practical contexts, from paints and coatings, to personal care, packaging, and electronics.

Role of heat generation and thermal diffusion during frontal photopolymerization.

Frontal photopolymerization (FPP) is a rapid and versatile solidification process that can be used to fabricate complex three-dimensional structures by selectively exposing a photosensitive monomer-rich bath to light. A characteristic feature of FPP is the appearance of a sharp polymerization front that propagates into the bath as a planar traveling wave. In this paper, we introduce a theoretical model to determine how heat generation during photopolymerization influences the kinetics of wave propagation as well as the monomer-to-polymer conversion profile, both of which are relevant for FPP applications and experimentally measurable. When thermal diffusion is sufficiently fast relative to the rate of polymerization, the system evolves as if it were isothermal. However, when thermal diffusion is slow, a thermal wavefront develops and propagates at the same rate as the polymerization front. This leads to an accumulation of heat behind the polymerization front which can result in a significant sharpening of the conversion profile and acceleration of the growth of the solid. Our results also suggest that a novel way to tailor the dynamics of FPP is by imposing a temperature gradient along the growth direction.

Controlling frontal photopolymerization with optical attenuation and mass diffusion.

Frontal photopolymerization (FPP) is a versatile directional solidification process that can be used to rapidly fabricate polymer network materials by selectively exposing a photosensitive monomer bath to light. A characteristic feature of FPP is that the monomer-to-polymer conversion profiles take on the form of traveling waves that propagate into the unpolymerized bulk from the illuminated surface. Practical implementations of FPP require detailed knowledge about the conversion profile and speed of these traveling waves. The purpose of this theoretical study is to (i) determine the conditions under which FPP occurs and (ii) explore how optical attenuation and mass transport can be used to finely tune the conversion profile and propagation kinetics. Our findings quantify the strong optical attenuation and slow mass transport relative to the rate of polymerization required for FPP. The shape of the traveling wave is primarily controlled by the magnitude of the optical attenuation coefficients of the neat and polymerized material. Unexpectedly, we find that mass diffusion can increase the net extent of polymerization and accelerate the growth of the solid network. The theoretical predictions are found to be in excellent agreement with experimental data acquired for representative systems.

A Unified Approach for Patterning via Frontal Photopolymerization.

A unified patterning strategy via frontal photopolymerization (FPP) that is robust to a wide range of radical photopolymerizing systems, including thiol-ene and acrylic monomers is reported. The factors governing the spatiotemporal solidification process, including front position, profile shape, and thermal effects, are investigated and modeled theoretically, resulting in the predictive FPP patterning of polymer networks with prescribed dimensions.