Optimal resource allocation in micro-organisms under periodic nutrient fluctuations

Although microorganisms often live in dynamic environments, most studies, both experimental and theoretical, are carried out under static conditions. In this work, we investigate the issue of optimal resource allocation in bacteria growing in periodic environments. We consider a dynamic model describing the microbial metabolism under varying conditions, involving a control variable quantifying the protein precursors allocation. Our objective is to determine the optimal strategies maximizing the long-term growth of cells under a piecewise-constant periodic environment. Firstly, we perform a theoretical analysis of the resulting optimal control problem (OCP), based on the application the Pontryagin's Maximum Principle (PMP). We determine that the structure of the optimal control must be bang-bang, with possibly some singular arcs corresponding to optimal equilibria of the system. If the control presents singular arcs, then these can only be reached and left through chattering arcs. We also use a direct optimization method, implemented in the BOCOP software, to solve the studied OCP. Our study reveals that the optimal solution over a large time horizon is related to the one over a single period of the varying environment with periodic constraints. Moreover, we observe that the maximal average growth rate attainable under periodic conditions can be higher than the one under a constant environment. We further extend our analysis to conduct a qualitative comparison between the predictions from our model and some recent biological experiments on E. coli. This analysis particularly highlights the mechanisms of action of the ppGpp signaling molecule, thus providing relevant explanations of the experimental observations. In conclusion, our study corroborates previous research indicating that this molecule plays a crucial role in the regulation of resource allocation of protein precursors in E. coli.

Nonlinear optimal control of a mean-field model of neural population dynamics

We apply the framework of nonlinear optimal control to a biophysically realistic neural mass model, which consists of two mutually coupled populations of deterministic excitatory and inhibitory neurons. External control signals are realized by time-dependent inputs to both populations. Optimality is defined by two alternative cost functions that trade the deviation of the controlled variable from its target value against the “strength” of the control, which is quantified by the integrated 1- and 2-norms of the control signal. We focus on a bistable region in state space where one low- (“down state”) and one high-activity (“up state”) stable fixed points coexist. With methods of nonlinear optimal control, we search for the most cost-efficient control function to switch between both activity states. For a broad range of parameters, we find that cost-efficient control strategies consist of a pulse of finite duration to push the state variables only minimally into the basin of attraction of the target state. This strategy only breaks down once we impose time constraints that force the system to switch on a time scale comparable to the duration of the control pulse. Penalizing control strength via the integrated 1-norm (2-norm) yields control inputs targeting one or both populations. However, whether control inputs to the excitatory or the inhibitory population dominate, depends on the location in state space relative to the bifurcation lines. Our study highlights the applicability of nonlinear optimal control to understand neuronal processing under constraints better.

Dynamic optimization reveals alveolar epithelial cells as key mediators of host defense in invasive aspergillosis

Aspergillus fumigatus is an important human fungal pathogen and its conidia are constantly inhaled by humans. In immunocompromised individuals, conidia can grow out as hyphae that damage lung epithelium. The resulting invasive aspergillosis is associated with devastating mortality rates. Since infection is a race between the innate immune system and the outgrowth of A. fumigatus conidia, we use dynamic optimization to obtain insight into the recruitment and depletion of alveolar macrophages and neutrophils. Using this model, we obtain key insights into major determinants of infection outcome on host and pathogen side. On the pathogen side, we predict in silico and confirm in vitro that germination speed is an important virulence trait of fungal pathogens due to the vulnerability of conidia against host defense. On the host side, we found that epithelial cells, which have been underappreciated, play a role in fungal clearance and are potent mediators of cytokine release. Both predictions were confirmed by in vitro experiments on established cell lines as well as primary lung cells. Further, our model affirms the importance of neutrophils in invasive aspergillosis and underlines that the role of macrophages remains elusive. We expect that our model will contribute to improvement of treatment protocols by focusing on the critical components of immune response to fungi but also fungal virulence traits.

A Review of Using Mathematical Modeling to Improve Our Understanding of Bacteriophage, Bacteria, and Eukaryotic Interactions

Phage therapy, the therapeutic usage of viruses to treat bacterial infections, has many theoretical benefits in the ‘post antibiotic era.’ Nevertheless, there are currently no approved mainstream phage therapies. One reason for this is a lack of understanding of the complex interactions between bacteriophage, bacteria and eukaryotic hosts. These three-component interactions are complex, with non-linear or synergistic relationships, anatomical barriers and genetic or phenotypic heterogeneity all leading to disparity between performance and efficacy in in vivo versus in vitro environments. Realistic computer or mathematical models of these complex environments are a potential route to improve the predictive power of in vitro studies for the in vivo environment, and to streamline lab work. Here, we introduce and review the current status of mathematical modeling and highlight that data on genetic heterogeneity and mutational stochasticity, time delays and population densities could be critical in the development of realistic phage therapy models in the future. With this in mind, we aim to inform and encourage the collaboration and sharing of knowledge and expertise between microbiologists and theoretical modelers, synergising skills and smoothing the road to regulatory approval and widespread use of phage therapy.

Implementation and acceleration of optimal control for systems biology

Optimal control theory provides insight into complex resource allocation decisions. The forward–backward sweep method (FBSM) is an iterative technique commonly implemented to solve two-point boundary value problems arising from the application of Pontryagin’s maximum principle (PMP) in optimal control. The FBSM is popular in systems biology as it scales well with system size and is straightforward to implement. In this review, we discuss the PMP approach to optimal control and the implementation of the FBSM. By conceptualizing the FBSM as a fixed point iteration process, we leverage and adapt existing acceleration techniques to improve its rate of convergence. We show that convergence improvement is attainable without prohibitively costly tuning of the acceleration techniques. Furthermore, we demonstrate that these methods can induce convergence where the underlying FBSM fails to converge. All code used in this work to implement the FBSM and acceleration techniques is available on GitHub at https://github.com/Jesse-Sharp/Sharp2021.

Designing and Constructing Artificial Small RNAs for Gene Regulation and Carbon Flux Redirection in Photosynthetic Cyanobacteria

Photosynthetic cyanobacteria are not only model organisms for studying photosynthesis and biological cycling of carbon in biosphere but also potential "green microbial factories" to produce renewable fuels and chemicals, due to their capability to utilizing solar energy and CO₂. Therefore, strategies for gene regulation and carbon flux redirection are important for both fundamental research and metabolic engineering of cyanobacteria. To address the challenges, regulatory tools based on artificial small RNAs have been developed with satisfactory effects for single or multiple gene(s) regulation in various cyanobacterial species. When combined with the promoters of varying gradient strength and the inducible switches developed in recent years, it is now feasible to realize precise gene regulation in photosynthetic cyanobacteria for producing fuels and chemicals. Here in this chapter, we provide a detailed introduction of the design principles and constructing methods of the artificial sRNA tools to achieve accurate inducible regulation of cyanobacterial gene(s).

Using optimal control to understand complex metabolic pathways

BACKGROUND

Optimality principles have been used to explain the structure and behavior of living matter at different levels of organization, from basic phenomena at the molecular level, up to complex dynamics in whole populations. Most of these studies have assumed a single-criteria approach. Such optimality principles have been justified from an evolutionary perspective. In the context of the cell, previous studies have shown how dynamics of gene expression in small metabolic models can be explained assuming that cells have developed optimal adaptation strategies. Most of these works have considered rather simplified representations, such as small linear pathways, or reduced networks with a single branching point, and a single objective for the optimality criteria.

RESULTS

Here we consider the extension of this approach to more realistic scenarios, i.e. biochemical pathways of arbitrary size and structure. We first show that exploiting optimality principles for these networks poses great challenges due to the complexity of the associated optimal control problems. Second, in order to surmount such challenges, we present a computational framework which has been designed with scalability and efficiency in mind, including mechanisms to avoid the most common pitfalls. Third, we illustrate its performance with several case studies considering the central carbon metabolism of S. cerevisiae and B. subtilis. In particular, we consider metabolic dynamics during nutrient shift experiments.

CONCLUSIONS

We show how multi-objective optimal control can be used to predict temporal profiles of enzyme activation and metabolite concentrations in complex metabolic pathways. Further, we also show how to consider general cost/benefit trade-offs. In this study we have considered metabolic pathways, but this computational framework can also be applied to analyze the dynamics of other complex pathways, such as signal transduction or gene regulatory networks.

Trends in mathematical modeling of host-pathogen interactions

Pathogenic microorganisms entail enormous problems for humans, livestock, and crop plants. A better understanding of the different infection strategies of the pathogens enables us to derive optimal treatments to mitigate infectious diseases or develop vaccinations preventing the occurrence of infections altogether. In this review, we highlight the current trends in mathematical modeling approaches and related methods used for understanding host-pathogen interactions. Since these interactions can be described on vastly different temporal and spatial scales as well as abstraction levels, a variety of computational and mathematical approaches are presented. Particular emphasis is placed on dynamic optimization, game theory, and spatial modeling, as they are attracting more and more interest in systems biology. Furthermore, these approaches are often combined to illuminate the complexities of the interactions between pathogens and their host. We also discuss the phenomena of molecular mimicry and crypsis as well as the interplay between defense and counter defense. As a conclusion, we provide an overview of method characteristics to assist non-experts in their decision for modeling approaches and interdisciplinary understanding.