  We analyze equilibrium behavior and optimal policy within a
Susceptible-Infected-Recovered epidemic model augmented with potentially
undiagnosed agents who infer their health status and a social planner with
imperfect enforcement of social distancing. We define and prove the existence
of a perfect Bayesian Markov competitive equilibrium and contrast it with the
efficient allocation subject to the same informational constraints. We identify
two externalities, static (individual actions affect current risk of infection)
and dynamic (individual actions affect future disease prevalence), and study
how they are affected by limitations on testing and enforcement. We prove that
a planner with imperfect enforcement will always wish to curtail activity, but
that its incentives to do so vanish as testing becomes perfect. When a vaccine
arrives far into the future, the planner with perfect enforcement may encourage
activity before herd immunity. We find that lockdown policies have modest
welfare gains, whereas quarantine policies are effective even with imperfect
testing.


  Protein engineering seeks to identify protein sequences with optimized
properties. When guided by machine learning, protein sequence generation
methods can draw on prior knowledge and experimental efforts to improve this
process. In this review, we highlight recent applications of machine learning
to generate protein sequences, focusing on the emerging field of deep
generative methods.


  Structure-based Deep Fusion models were recently shown to outperform several
physics- and machine learning-based protein-ligand binding affinity prediction
methods. As part of a multi-institutional COVID-19 pandemic response, over 500
million small molecules were computationally screened against four protein
structures from the novel coronavirus (SARS-CoV-2), which causes COVID-19.
Three enhancements to Deep Fusion were made in order to evaluate more than 5
billion docked poses on SARS-CoV-2 protein targets. First, the Deep Fusion
concept was refined by formulating the architecture as one, coherently
backpropagated model (Coherent Fusion) to improve binding-affinity prediction
accuracy. Secondly, the model was trained using a distributed, genetic
hyper-parameter optimization. Finally, a scalable, high-throughput screening
capability was developed to maximize the number of ligands evaluated and
expedite the path to experimental evaluation. In this work, we present both the
methods developed for machine learning-based high-throughput screening and
results from using our computational pipeline to find SARS-CoV-2 inhibitors.


  We show how the standard field theoretical language based on creation and
annihilation operators may be used for a straightforward derivation of an
SIR-type stochastic model for COVID-19 epidemic, from which we obtain the time
evolution of the mean number of infectious (active cases) and deceased
individuals. In order to capture the effects of lockdown and social distancing,
we use a time-dependent infection rate. The results are in good agreement with
the data for three different waves of epidemic activity in South Korea.


  Although reproducibility is a core tenet of the scientific method, it remains
challenging to reproduce many results. Surprisingly, this also holds true for
computational results in domains such as systems biology where there have been
extensive standardization efforts. For example, Tiwari et al. recently found
that they could only repeat 50% of published simulation results in systems
biology. Toward improving the reproducibility of computational systems
research, we identified several resources that investigators can leverage to
make their research more accessible, executable, and comprehensible by others.
In particular, we identified several domain standards and curation services, as
well as powerful approaches pioneered by the software engineering industry that
we believe many investigators could adopt. Together, we believe these
approaches could substantially enhance the reproducibility of systems biology
research. In turn, we believe enhanced reproducibility would accelerate the
development of more sophisticated models that could inform precision medicine
and synthetic biology.


  The response of many governments to the COVID-19 pandemic has involved
measures to control within- and between-household transmission, providing
motivation to improve understanding of the absolute and relative risks in these
contexts. Here, we perform exploratory, residual-based, and
transmission-dynamic household analysis of the Office for National Statistics
(ONS) COVID-19 Infection Survey (CIS) data from 26 April 2020 to 15 July 2021
in England. This provides evidence for: (i) temporally varying rates of
introduction of infection into households broadly following the trajectory of
the overall epidemic and vaccination programme; (ii) Susceptible-Infectious
Transmission Probabilities (SITPs) of within-household transmission in the
15-35% range; (iii) the emergence of the Alpha and Delta variants, with the
former being around 50% more infectious than wildtype and 35% less infectious
than Delta within households; (iv) significantly (in the range 25-300%) more
risk of bringing infection into the household for workers in patient-facing
roles pre-vaccine; (v) increased risk for secondary school-age children of
bringing the infection into the household when schools are open; (vi) increased
risk for primary school-age children of bringing the infection into the
household when schools were open since the emergence of new variants.


  Biomedical networks (or graphs) are universal descriptors for systems of
interacting elements, from molecular interactions and disease co-morbidity to
healthcare systems and scientific knowledge. Advances in artificial
intelligence, specifically deep learning, have enabled us to model, analyze,
and learn with such networked data. In this review, we put forward an
observation that long-standing principles of systems biology and medicine --
while often unspoken in machine learning research -- provide the conceptual
grounding for representation learning on graphs, explain its current successes
and limitations, and even inform future advancements. We synthesize a spectrum
of algorithmic approaches that, at their core, leverage graph topology to embed
networks into compact vector spaces. We also capture the breadth of ways in
which representation learning has dramatically improved the state-of-the-art in
biomedical machine learning. Exemplary domains covered include identifying
variants underlying complex traits, disentangling behaviors of single cells and
their effects on health, assisting in diagnosis and treatment of patients, and
developing safe and effective medicines.


  Data-driven artificial intelligence models fed with published scientific
findings have been used to create powerful prediction engines for scientific
and technological advance, such as the discovery of novel materials with
desired properties and the targeted invention of new therapies and vaccines.
These AI approaches typically ignore the distribution of human prediction
engines -- scientists and inventor -- who continuously alter the landscape of
discovery and invention. As a result, AI hypotheses are designed to substitute
for human experts, failing to complement them for punctuated collective
advance. Here we show that incorporating the distribution of human expertise
into self-supervised models by training on inferences cognitively available to
experts dramatically improves AI prediction of future human discoveries and
inventions. Including expert-awareness into models that propose (a) valuable
energy-relevant materials increases the precision of materials predictions by
~100%, (b) repurposing thousands of drugs to treat new diseases increases
precision by 43%, and (c) COVID-19 vaccine candidates examined in clinical
trials by 260%. These models succeed by predicting human predictions and the
scientists who will make them. By tuning AI to avoid the crowd, however, it
generates scientifically promising "alien" hypotheses unlikely to be imagined
or pursued without intervention, not only accelerating but punctuating
scientific advance. By identifying and correcting for collective human bias,
these models also suggest opportunities to improve human prediction by
reformulating science education for discovery.


  We studied simple random-walk models with asymmetric time delays. Stochastic
simulations were performed for hyperbolic-tangent fitness functions and to
obtain analytical results we approximated them by step functions. A novel
behavior has been observed. Namely, the mean position of a walker depends on
time delays. This is a joint effect of both stochasticity and time delays
present in the system. We also observed that by shifting appropriately fitness
functions we may reverse the effect of time delays - the mean position of the
walker changes the sign.


  Most countries have started vaccinating people against COVID-19. However, due
to limited production capacities and logistical challenges it will take
months/years until herd immunity is achieved. Therefore, vaccination and social
distancing have to be coordinated. In this paper, we provide some insight on
this topic using optimization-based control on an age-differentiated
compartmental model. For real-life decision making, we investigate the impact
of the planning horizon on the optimal vaccination/social distancing strategy.
We find that in order to reduce social distancing in the long run, without
overburdening the healthcare system, it is essential to vaccinate the people
with the highest contact rates first. That is also the case if the objective is
to minimize fatalities provided that the social distancing measures are
sufficiently strict. However, for short-term planning it is optimal to focus on
the high-risk group.


  Non-equilibrium conditions must have been crucial for the assembly of the
first informational polymers of early life, but supporting their formation and
continuous enrichment in a long-lasting environment. Here, we explore how gas
bubbles in water subjected to a thermal gradient, a likely scenario within
crustal mafic rocks on the early Earth, drive a complex, continuous enrichment
of prebiotic molecules. NRA precursors, monomers, active ribozymes,
oligonucleotides and lipids are shown to (1) cycle between dry and wet states,
enabling the central step of RNA phosphorylation, (2) accumulate at the
gas-water interface to drastically increase ribozymatic activity, (3) condense
into hydrogels, (4) form pure crystals and (5) encapsulate into protecting
vesicle aggregates that subsequently undergo fission. These effects occur
within less than 30 min. The findings unite, in one location, the physical
conditions that were crucial for the chemical emergence of biopolymers. They
suggest that heated microbubbles could have hosted the first cycles of
molecular evolution.


  Objective: To evaluate the relationship between coronavirus disease 2019
(COVID-19) diagnosis with SARS-CoV-2 variant B.1.1.7 (also known as Variant of
Concern 202012/01) and the risk of hospitalisation compared to diagnosis with
wildtype SARS-CoV-2 variants.
  Design: Retrospective cohort, analysed using stratified Cox regression.
  Setting: Community-based SARS-CoV-2 testing in England, individually linked
with hospitalisation data.
  Participants: 839,278 laboratory-confirmed COVID-19 patients, of whom 36,233
had been hospitalised within 14 days, tested between 23rd November 2020 and
31st January 2021 and analysed at a laboratory with an available TaqPath assay
that enables assessment of S-gene target failure (SGTF). SGTF is a proxy test
for the B.1.1.7 variant. Patient data were stratified by age, sex, ethnicity,
deprivation, region of residence, and date of positive test.
  Main outcome measures: Hospitalisation between 1 and 14 days after the first
positive SARS-CoV-2 test.
  Results: 27,710 of 592,409 SGTF patients (4.7%) and 8,523 of 246,869 non-SGTF
patients (3.5%) had been hospitalised within 1-14 days. The stratum-adjusted
hazard ratio (HR) of hospitalisation was 1.52 (95% confidence interval [CI]
1.47 to 1.57) for COVID-19 patients infected with SGTF variants, compared to
those infected with non-SGTF variants. The effect was modified by age
(P<0.001), with HRs of 0.93-1.21 for SGTF compared to non-SGTF patients below
age 20 years, 1.29 in those aged 20-29, and 1.45-1.65 in age groups 30 years or
older.
  Conclusions: The results suggest that the risk of hospitalisation is higher
for individuals infected with the B.1.1.7 variant compared to wildtype
SARS-CoV-2, likely reflecting a more severe disease. The higher severity may be
specific to adults above the age of 30.


  NEMF is a novel network-based ecosystem modelling framework. It is a flexible
and easy-to-use tool for modelling ecosystems with low- to intermediate
complexity. It is designed around the idea of visualizing an ecosystem through
a network that implicitly defines a set of differential equations. These
equations are then solved numerically and can also be inverse modelled. The
framework offers the functionality to handle non-equilibrium, non-linear
interactions.


  Many countries are managing COVID-19 epidemic by switching between lighter
and heavier restrictions. While an open-close and a close-open cycle have
comparable socio-economic costs, the former leads to a much heavier burden in
terms of deaths and pressure on the healthcare system. An empirical
demonstration of the toll ensuing from procrastination was recently observed in
Israel, where both cycles were enforced from late August to mid-December 2020,
yielding some 1,600 deaths with open-close compared to 440 with close-open.


  Complex socio-environmental interdependencies drive biological invasions,
causing damages across large spatial scales. For widespread invasions,
targeting of management activities based on optimization approaches may fail
due to computational or data constraints. Here we evaluate an alternative
approach that embraces complexity by representing the invasion as a network and
using network structure to inform management locations. We compare optimal
versus network-guided invasive species management at a landscape-scale,
considering siting of boat decontamination stations targeting 1.6 million
boater movements among 9,182 lakes in Minnesota, USA. Studying performance for
58 counties, we find that when full information is known on invasion status and
boater movements, the best-performing network-guided metric achieves a median
and lower quartile performance of 100% of optimal. We also find that
performance remains relatively high using different network metrics or with
less information (median above 80% and lower quartile above 60% of optimal for
most metrics), but is more variable, particularly at the lower quartile.
Additionally, performance is generally stable across counties with varying lake
counts, suggesting viability for large-scale invasion management.


  Polymer models are a widely used tool to study the prebiotic formation of
metabolism at the origins of life. Counts of the number of reactions in these
models are often crucial in probabilistic arguments concerning the emergence of
autocatalytic networks. In the first part of this paper, we provide the first
exact description of the number of reactions under widely applied model
assumptions. Conclusions from earlier studies rely on either approximations or
asymptotic counting, and we show that the exact counts lead to similar, though
not always identical, asymptotic results. In the second part of the paper, we
investigate a novel model assumption whereby polymers are invariant under
spatial rotation. We outline the biochemical relevance of this condition and
again give exact enumerative and asymptotic formulae for the number of
reactions.


  Geographic ranges of communities of species evolve in response to
environmental, ecological, and evolutionary forces. Understanding the effects
of these forces on species' range dynamics is a major goal of spatial ecology.
Previous mathematical models have jointly captured the dynamic changes in
species' population distributions and the selective evolution of
fitness-related phenotypic traits in the presence of an environmental gradient.
These models inevitably include some unrealistic assumptions, and biologically
reasonable ranges of values for their parameters are not easy to specify. As a
result, simulations of the seminal models of this type can lead to markedly
different conclusions about the behavior of such populations, including the
possibility of maladaptation setting stable range boundaries. Here, we
harmonize such results by developing and simulating a continuum model of range
evolution in a community of species that interact competitively while diffusing
over an environmental gradient. Our model extends existing models by
incorporating both competition and freely changing intraspecific trait
variance. Simulations of this model predict a spatial profile of species' trait
variance that is consistent with experimental measurements available in the
literature. Moreover, they reaffirm interspecific competition as an effective
factor in limiting species' ranges, even when trait variance is not
artificially constrained. These theoretical results can inform the design of,
as yet rare, empirical studies to clarify the evolutionary causes of range
stabilization.


  A multiscale mathematical model is presented to describe the de novo
granulation and the evolution of multispecies granular biofilms within a
continuous reactor. The granule is modelled as a spherical free boundary domain
with radial symmetry. The equation which governs the free boundary is derived
from global mass balance considerations and takes into account the growth of
sessile biomass and the exchange fluxes with the bulk liquid. Starting from a
vanishing initial value, the expansion of the free boundary is initiated by the
attachment process, which depends on the microbial species concentrations
within the bulk liquid and their specific attachment velocity. Nonlinear
hyperbolic PDEs model the growth of the sessile microbial species, while
quasi-linear parabolic PDEs govern the dynamics of substrates and invading
species within the granular biofilm. Nonlinear ODEs govern the evolution of
soluble substrates and planktonic biomass within the bulk liquid. The model is
applied to an anaerobic granular-based system and solved numerically to test
its qualitative behaviour and explore the main aspects of de novo anaerobic
granulation: ecology, biomass distribution, relative abundance, dimensional
evolution of the granules and soluble substrates and planktonic biomass
dynamics within the reactor. The numerical results confirm that the model
accurately describes the ecology and the concentrically-layered structure of
anaerobic granules observed experimentally, and is able to predict the effects
of some significant factors, such as influent wastewater composition,
granulation properties of planktonic biomass, biomass density and hydrodynamic
and shear stress conditions, on the process performance.


  Many biological and social systems show significant levels of collective
action. Several cooperation mechanisms have been proposed, yet they have been
mostly studied independently. Among these, direct reciprocity supports
cooperation on the basis of repeated interactions among individuals. Signals
and quorum dynamics may also drive cooperation. Here, we resort to an
evolutionary game theoretical model to jointly analyse these two mechanisms and
study the conditions in which evolution selects for direct reciprocity,
signalling, or their combination. We show that signalling alone leads to higher
levels of cooperation than when combined with reciprocity, while offering
additional robustness against errors. Specifically, successful strategies in
the realm of direct reciprocity are often not selected in the presence of
signalling, and memory of past interactions is only exploited opportunistically
in case of earlier coordination failure. Differently, signalling always
evolves, even when costly. In the light of these results, it may be easier to
understand why direct reciprocity has been observed only in a limited number of
cases among non-humans, whereas signalling is widespread at all levels of
complexity.


  Data-driven dynamic models of cell biology can be used to predict cell
response to unseen perturbations. Recent work (CellBox) had demonstrated the
derivation of interpretable models with explicit interaction terms, in which
the parameters were optimized using machine learning techniques. While the
previous work was tested only in a single biological setting, this work aims to
extend the range of applicability of this model inference approach to a
diversity of biological systems. Here we adapted CellBox in Julia differential
programming and augmented the method with adjoint algorithms, which has
recently been used in the context of neural ODEs. We trained the models using
simulated data from both abstract and biology-inspired networks, which afford
the ability to evaluate the recovery of the ground truth network structure. The
resulting accuracy of prediction by these models is high both in terms of low
error against data and excellent agreement with the network structure used for
the simulated training data. While there is no analogous ground truth for real
life biological systems, this work demonstrates the ability to construct and
parameterize a considerable diversity of network models with high predictive
ability. The expectation is that this kind of procedure can be used on real
perturbation-response data to derive models applicable to diverse biological
systems.


  We present an approach for automated in-situ monitoring of phytoplankton and
zooplankton communities based on a dual magnification dark-field imaging
microscope/camera. We describe the Dual Scripps Plankton Camera (DSPC) system
and associated image processing, and assess its capabilities in detecting and
characterizing plankton species of different size and taxonomic categories, and
in measuring their abundances in both laboratory and field applications. In the
laboratory, body size and abundance estimates by the DSPC significantly and
robustly scale with the same measurements derived by traditional microscopy. In
the field, a DSPC installed permanently at 3 m depth in Lake Greifensee
(Switzerland), delivered images of plankton individuals, colonies, and
heterospecific aggregates without disrupting natural arrangements of
interacting organisms, their microenvironment or their behavior at hourly
timescales. The DSPC was able to track the dynamics of taxa in the size range
between ~10 $\mu$m to ~ 1 cm, covering virtually all the components of the
planktonic food web (including parasites and potentially toxic cyanobacteria).
Comparing data from the field-deployed DSPC to traditional sampling and
microscopy revealed a general overall agreement in estimates of plankton
diversity and abundances, despite imaging limitations in detecting small
phytoplankton species and rare and large zooplankton taxa (e.g. carnivorous
zooplankton). The most significant disagreements between traditional methods
and the DSPC resided in the measurements of community properties of
zooplankton, organisms that are heterogeneously distributed spatially and
temporally, and whose demography appeared to be better captured by automated
imaging. Time series collected by the DSPC depicted ecological succession
patterns, algal bloom dynamics and circadian fluctuations with a temporal
frequency and morphological [continues...]


  Predictive modeling is the key factor for saving time and resources with
respect to manufacturing processes such as fermentation processes arising e.g.\
in food and chemical manufacturing processes. According to Zhang et al. (2002),
the open-loop dynamics of yeast are highly dependent on the initial cell mass
distribution. This can be modeled via population balance models describing the
single-cell behavior of the yeast cell. There have already been several
population balance models for wine fermentation in the literature. However, the
new model introduced in this paper is much more detailed than the ones studied
previously. This new model for the white wine fermentation process is based on
a combination of components previously introduced in literature. It turns it
into a system of highly nonlinear weakly hyperbolic partial/ordinary
integro-differential equations. This model becomes very challenging from a
theoretical and numerical point of view. Existence and uniqueness of solutions
to a simplified version of the introduced problem is studied based on semigroup
theory. For its numerical solution a numerical methodology based on a finite
volume scheme combined with a time implicit scheme is derived. The impact of
the initial cell mass distribution on the solution is studied and underlined
with numerical results. The detailed model is compared to a simpler model based
on ordinary differential equations. The observed differences for different
initial distributions and the different models turn out to be smaller than
expected. The outcomes of this paper are very interesting and useful for
applied mathematicians, winemakers and process engineers.


  Major advances in public health have resulted from disease prevention.
However, prevention of a new infectious disease by vaccination or
pharmaceuticals is made difficult by the slow process of vaccine and drug
development. We propose an additional intervention that allows rapid control of
emerging infectious diseases, and can also be used to eradicate diseases that
rely almost exclusively on human-to-human transmission. The intervention is
based on (1) testing every individual for the disease, (2) repeatedly, and (3)
isolation of infected individuals. We show here that at a sufficient rate of
testing, the reproduction number is reduced below 1.0 and the epidemic will
rapidly collapse. The approach does not rely on strong or unrealistic
assumptions about test accuracy, isolation compliance, population structure or
epidemiological parameters, and its success can be monitored in real time by
following the test positivity rate. In addition to the compliance rate and
false negatives, the required rate of testing depends on the design of the
testing regime, with concurrent testing outperforming random sampling. Provided
that results are obtained rapidly, the test frequency required to suppress an
epidemic is monotonic and near-linear with respect to R0, the infectious
period, and the fraction of susceptible individuals. The testing regime is
effective against both early phase and established epidemics, and additive to
other interventions (e.g. contact tracing and social distancing). It is also
robust to failure: any rate of testing reduces the number of infections,
improving both public health and economic conditions. These conclusions are
based on rigorous analysis and simulations of appropriate epidemiological
models. A mass-produced, disposable test that could be used at home would be
ideal, due to the optimal performance of concurrent tests that return immediate
results.


  It is largely taken for granted that differential abundance analysis is, by
default, the best first step when analyzing genomic data. We argue that this is
not necessarily the case. In this article, we identify key limitations that are
intrinsic to differential abundance analysis: it is (a) dependent on
unverifiable assumptions, (b) an unreliable construct, and (c) overly
reductionist. We formulate an alternative framework called ratio-based
biomarker analysis which does not suffer from the identified limitations.
Moreover, ratio-based biomarkers are highly flexible. Beyond replacing DAA,
they can also be used for many other bespoke analyses, including dimension
reduction and multi-omics data integration.


  Biocomputing systems based on engineered bacteria can lead to novel tools for
environmental monitoring and detection of metabolic diseases. In this paper, we
propose a Bacterial Molecular Computing on a Chip (BMCoC) using microfluidic
and electrochemical sensing technologies. The computing can be flexibly
integrated into the chip, but we focus on engineered bacterial AND Boolean
logic gate and ON-OFF switch sensors that produces secondary signals to change
the pH and dissolved oxygen concentrations. We present a prototype with
experimental results that shows the electrochemical sensors can detect small pH
and dissolved oxygen concentration changes created by the engineered bacterial
populations' molecular signals. Additionally, we present a theoretical model
analysis of the BMCoC computation reliability when subjected to unwanted
effects, i.e., molecular signal delays and noise, and electrochemical sensors
threshold settings that are based on either standard or blind detectors. Our
numerical analysis found that the variations in the production delay and the
molecular output signal concentration can impact on the computation reliability
for the AND logic gate and ON-OFF switch. The molecular communications of
synthetic engineered cells for logic gates integrated with sensing systems can
lead to a new breed of biochips that can be used for numerous diagnostic
applications.


  The Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2), or
Covid-19, burst into a pandemic in the beginning of 2020. An unprecedented
worldwide effort involving academic institutions, regulatory agencies and
industry is facing the challenges imposed by the rapidly spreading disease.
Emergency use authorization for vaccines were granted in the beginning of
December 2020 in Europe and nine days later in the United States. The urge for
vaccination started a race, forcing governs and health care agencies to take
decisions on the fly regarding the vaccination strategy and logistics. So far,
the vaccination strategies and non-pharmaceutical interventions, such as social
distancing and the use of face masks, are the only efficient actions to stop
the pandemic. In this context, it is of fundamental importance to understand
the dynamical behavior of the Covid-19 spread along with possible vaccination
strategies. In this work a Susceptible - Infected - Removed - Sick with
vaccination (SIRSi-Vaccine) model is proposed. In addtion, the SIRSi-Vaccine
model also accounts for unreported, or asymptomatic, cases and the possibility
of temporary immunity, either after infection or vaccination. Disease free and
endemic equilibrium points existence conditions are determined in the (! ? ?)
vaccine-effort and social distancing parameter space. The model is adjusted to
the data from S\~ao Paulo, Santos and Campinas, three major cities in the State
of S\~ao Paulo, Brazil.


  In this work we look at several mathematical models that have been
constructed during the present pandemic to address dfferent issues of
importance to public health policies about epidemic scenarios and thier causes.
We start by briefly reviewing the most basic properties of the classic
Kermack-McKendrick models and then proceed to look at some generalizations and
their applications to contact structures, cocirculation of viral infections,
growth patterns of epidemic curves, characterization of probability
distributions and passage times necessary to parametrize the compartments that
are added to the basic Kermack-McKendrick model. All of these examples revolve
around the idea that a system of differential equations constructed for an
specifc epidemiological problem, has as central and main theoretical and
conceptual support the epidemiological, medical and biological context that
motivates its construction and analysis. Fitting differential models to data
without a sound perspective on the biological feasibility of the mathematical
model will provide an statistical significant fit but will also fall short of
providing useful information for the management, mitigation and eventual
control of the epidemic of SARS-CoV-2.


  Since the start of the still ongoing COVID-19 pandemic, there have been many
modeling efforts to assess several issues of importance to public health. In
this work, we review the theory behind some important mathematical models that
have been used to answer questions raised by the development of the pandemic.
We start revisiting the basic properties of simple Kermack-McKendrick type
models. Then, we discuss extensions of such models and important
epidemiological quantities applied to investigate the role of heterogeneity in
disease transmission e.g. mixing functions and superspreading events, the
impact of non-pharmaceutical interventions in the control of the pandemic,
vaccine deployment, herd-immunity, viral evolution and the possibility of
vaccine escape. From the perspective of mathematical epidemiology, we highlight
the important properties, findings, and, of course, deficiencies, that all
these models have.


  Background: Gene regulatory networks coordinate the expression of genes
across physiological states and ensure a synchronized expression of genes in
cellular subsystems, critical for the coherent functioning of cells. Here we
address the questions whether it is possible to predict gene synchronization
from network structure alone. We have recently shown that synchronized gene
expression may be predicted from symmetries in the gene regulatory networks
(GRN) and described by the concept of symmetry fibrations. We showed that
symmetry fibrations partition the genes into groups called fibers based on the
symmetries of their 'input trees', the set of paths in the network through
which signals can reach a gene. In idealized dynamic gene expression models,
all genes in a fiber are perfectly synchronized, while less idealized models -
with gene input functions differencing between genes - predict symmetry
breaking and desynchronization.
  Results: To study the functional role of gene fibers and to test whether some
of the fiber-induced coexpression remains in reality, we analyze gene
fibrations for the gene regulatory networks of E. coli and B. subtilis and
confront them with expression data. We find approximate gene coexpression
patterns consistent with symmetry fibrations with idealized gene expression
dynamics. This shows that network structure alone provides useful information
about gene synchronization, and suggest that gene input functions within fibers
may be further streamlined by evolutionary pressures to realize a coexpression
of genes.
  Conclusions: Thus, gene fibrations provides a sound conceptual tool to
describe tunable coexpression induced by network topology and shaped by
mechanistic details of gene expression.


  The endogenous tridecapeptide neurotensin (NT) has emerged as an important
inhibitory modulator of pain transmission, exerting its analgesic action
through the activation of the G protein-coupled receptors, NTS1 and NTS2.
Whereas both NT receptors mediate the analgesic effects of NT, NTS1 activation
also produces hypotension and hypothermia, which may represent obstacles for
the development of new pain medications. In the present study, we implemented
various chemical strategies to improve the metabolic stability of the
biologically active fragment NT(8-13) and assessed their NTS1/NTS2 relative
binding affinities. We then determined their ability to reduce the nociceptive
behaviors in acute, tonic, and chronic pain models and to modulate blood
pressure and body temperature. To this end, we synthesized a series of NT(8-13)
analogs carrying a reduced amide bond at Lys8-Lys9 and harboring site-selective
modifications with unnatural amino acids, such as silaproline (Sip) and
trimethylsilylalanine (TMSAla). Incorporation of Sip and TMSAla respectively in
positions 10 and 13 of NT(8-13) combined with the Lys8-Lys9 reduced amine bond
(JMV5296) greatly prolonged the plasma half-life time over 20 hours. These
modifications also led to a 25-fold peptide selectivity toward NTS2. More
importantly, central delivery of JMV5296 was able to induce a strong
antinociceptive effect in acute (tail-flick), tonic (formalin), and chronic
inflammatory (CFA) pain models without inducing hypothermia. Altogether, these
results demonstrate that the chemically-modified NT(8-13) analog JMV5296
exhibits a better therapeutic profile and may thus represent a promising avenue
to guide the development of new stable NT agonists and improve pain management.


  In recent years, fossil-based fuels have been used to supply the energy needs
of the world. Fossil-based fuels induce accumulation of the atmospheric CO2
which causes global warming. One of CO2 source is flue gas emission from the
power plant. The microalgae have been considered excellent biological materials
for reduction CO2 with ability to photosynthesis. In this study, air, 10% CO2,
15% CO2 and simulated flue gas (containing 15% CO2) feed were used to observe
effect CO2 concentration and flue gas on cell growth and lipid content of
Chlorella protothecoides. The highest dry cell weight (1,5 g/L) and lipid
content (45%) values were obtained with 15% CO2 feed while the highest growth
rate (1,10), biomass productivity (0,125 g/L/day), and lipid weight (0,63 g/g)
were observed in 10% CO2 feed. Cultures fed with flue gas did not inhibited C.
protothecoides growth and showed similar results with those fed with 15% CO2
gas in terms of growth rate, dry cell weight, biomass productivity and lipid
content. These results showed that C. protothecoides has great potential for
reducing CO2 emission from flue gas.


  Can egalitarian norms or conventions survive the presence of dominant
individuals who are ensured of victory in conflicts? We investigate the
interaction of power asymmetry and partner choice in games of conflict over a
contested resource. We introduce three models to study the emergence and
resilience of cooperation among unequals when interaction is random, when
individuals can choose their partners, and where power asymmetries dynamically
depend on accumulated payoffs. We find that the ability to avoid bullies with
higher competitive ability afforded by partner choice mostly restores
cooperative conventions and that the competitive hierarchy never forms. Partner
choice counteracts the hyper dominance of bullies who are isolated in the
network and eliminates the need for others to coordinate in a coalition. When
competitive ability dynamically depends on cumulative payoffs, complex cycles
of coupled network-strategy-rank changes emerge. Effective collaborators gain
popularity (and thus power), adopt aggressive behavior, get isolated, and
ultimately lose power. Neither the network nor behavior converge to a stable
equilibrium. Despite the instability of power dynamics, the cooperative
convention in the population remains stable overall and long-term inequality is
completely eliminated. The interaction between partner choice and dynamic power
asymmetry is crucial for these results: without partner choice, bullies cannot
be isolated, and without dynamic power asymmetry, bullies do not lose their
power even when isolated. We analytically identify a single critical point that
marks a phase transition in all three iterations of our models. This critical
point is where the first individual breaks from the convention and cycles start
to emerge.


  Identifying novel functional protein structures is at the heart of molecular
engineering and molecular biology, requiring an often computationally
exhaustive search. We introduce the use of a Deep Convolutional Generative
Adversarial Network (DCGAN) to classify protein structures based on their
functionality by encoding each sample in a grid object structure using three
features in each object: the generic atom type, the position atom type, and its
occupancy relative to a given atom. We train DCGAN on 3-dimensional (3D) decoy
and native protein structures in order to generate and discriminate 3D protein
structures. At the end of our training, loss converges to a local minimum and
our DCGAN can annotate functional proteins robustly against adversarial protein
samples. In the future we hope to extend the novel structures we found from the
generator in our DCGAN with more samples to explore more granular functionality
with varying functions. We hope that our effort will advance the field of
protein structure prediction.


  As antibiotic resistance grows more frequent for common bacterial infections,
alternative treatment strategies such as phage therapy have become more widely
studied in the medical field. While many studies have explored the efficacy of
antibiotics, phage therapy, or synergistic combinations of phages and
antibiotics, the impact of virus competition on the efficacy of antibiotic
treatment has not yet been considered. Here, we model the synergy between
antibiotics and two viral types, temperate and chronic, in controlling
bacterial infections. We demonstrate that while combinations of antibiotic and
temperate viruses exhibit synergy, competition between temperate and chronic
viruses inhibits bacterial control with antibiotics. In fact, our model reveals
that antibiotic treatment may counterintuitively increase the bacterial load
when a large fraction of the bacteria develop antibiotic-resistance.


  Tropical geometry with the max-plus algebra has been applied to statistical
learning models over tree spaces because geometry with the tropical metric over
tree spaces has some nice properties such as convexity in terms of the tropical
metric. One of the challenges in applications of tropical geometry to tree
spaces is the difficulty interpreting outcomes of statistical models with the
tropical metric. This paper focuses on combinatorics of tree topologies along a
tropical line segment, an intrinsic geodesic with the tropical metric, between
two phylogenetic trees over the tree space and we show some properties of a
tropical line segment between two trees. Specifically we show that a
probability of a tropical line segment of two randomly chosen trees going
through the origin (the star tree) is zero if the number of leave is greater
than four, and we also show that if two given trees differ only one nearest
neighbor interchange (NNI) move, then the tree topology of a tree in the
tropical line segment between them is the same tree topology of one of these
given two trees with possible zero branch lengths.


  We discuss in this survey several network modeling methods and their
applicability to precision medicine. We review several network centrality
methods (degree centrality, closeness centrality, eccentricity centrality,
betweenness centrality, and eigenvector-based prestige) and two systems
controllability methods (minimum dominating sets and network structural
controllability). We demonstrate their applicability to precision medicine on
three multiple myeloma patient disease networks. Each network consists of
protein-protein interactions built around a specific patient's mutated genes,
around the targets of the drugs used in the standard of care in multiple
myeloma, and around multiple myeloma-specific essential genes. For each network
we demonstrate how the network methods we discuss can be used to identify
personalized, targeted drug combinations uniquely suited to that patient.


  Ecosystem services are the direct and indirect contributions of an ecosystem
to human well-being and survival. Ecosystem valuation is a method of assigning
a monetary value to an ecosystem with its goods and services,often referred to
as ecosystem service value (ESV). With the rapid expansion of cities, a
mismatch occurs between urban development and ecological development, and it is
increasingly urgent to establish a valid ecological assessment method. In this
study, we propose an ecological evaluation standard framework by designing an
ESV monitoring workflow based on the establishment of multi-level grids. The
proposed method is able to capture multi-scale features, facilitates
multi-level spatial expression, and can effectively reveal the spatial
heterogeneity of ESV. Taking Haian city in the Jiangsu province as the study
case, we implemented the proposed dynamic multi-level grids-based (DMLG) to
calculate its urban ESV in 2016 and 2019. We found that the ESV of Haian city
showed considerable growth (increased by 24.54 million RMB). Negative ESVs are
concentrated in the central city, which presented a rapid trend of outward
expansion. The results illustrated that the ongoing urban expanse does not
reduce the ecological value in the study area. The proposed unified grid
framework can be applied to other geographical regions and is expected to
benefit future studies in ecosystem service evaluation in terms of capture
multi-level spatial heterogeneity.


  Many studies have proposed a relationship between COVID-19 transmissibility
and ambient pollution levels. However, a major limitation in establishing such
associations is to adequately account for complex disease dynamics, influenced
by e.g. significant differences in control measures and testing policies.
Another difficulty is appropriately controlling the effects of other
potentially important factors, due to both their mutual correlations and a
limited dataset. To overcome these difficulties, we will here use the basic
reproduction number ($R_0$) that we estimate for USA states using non-linear
dynamics methods. To account for a large number of predictors (many of which
are mutually strongly correlated), combined with a limited dataset, we employ
machine-learning methods. Specifically, to reduce dimensionality without
complicating the variable interpretation, we employ Principal Component
Analysis on subsets of mutually related (and correlated) predictors. Methods
that allow feature (predictor) selection, and ranking their importance, are
then used, including both linear regressions with regularization and feature
selection (Lasso and Elastic Net) and non-parametric methods based on ensembles
of weak-learners (Random Forest and Gradient Boost). Through these
substantially different approaches, we robustly obtain that PM$_{2.5}$ is a
major predictor of $R_0$ in USA states, with corrections from factors such as
other pollutants, prosperity measures, population density, chronic disease
levels, and possibly racial composition. As a rough magnitude estimate, we
obtain that a relative change in $R_0$, with variations in pollution levels
observed in the USA, is typically ~30%, which further underscores the
importance of pollution in COVID-19 transmissibility.


  Most experimental sciences now rely on computing, and biological sciences are
no exception. As datasets get bigger, so do the computing costs, making proper
optimization of the codes used by scientists increasingly important. Many of
the codes developed in recent years are based on the Python-based NumPy, due to
its ease of use and good performance characteristics. The composable nature of
NumPy, however, does not generally play well with the multi-tier nature of
modern CPUs, making any non-trivial multi-step algorithm limited by the
external memory access speeds, which are hundreds of times slower than the
CPU's compute capabilities. In order to fully utilize the CPU compute
capabilities, one must keep the working memory footprint small enough to fit in
the CPU caches, which requires splitting the problem into smaller portions and
fusing together as many steps as possible. In this paper, we present changes
based on these principles to two important functions in the scikit-bio library,
principal coordinates analysis and the Mantel test, that resulted in over 100x
speed improvement in these widely used, general-purpose tools.


  In a world being hit by waves of COVID-19, vaccination is a light on the
horizon. However, the roll-out of vaccination strategies and their influence on
the pandemic are still open problems. In order to compare the effect of various
strategies proposed by the World Health Organization and other authorities, a
previously developed SEIRS stochastic model of geographical spreading of the
virus is extended by adding a compartment for vaccinated people. The parameters
of the model were fitted to describe the pandemic evolution in Argentina,
Mexico and Spain to analyze the effect of the proposed vaccination strategies.
The mobility parameters allow to simulate different social behaviors (e.g.
lock-down interventions). Schemes in which vaccines are applied homogeneously
in all the country, or limited to the most densely-populated areas, are
simulated and compared. The second strategy is found to be more effective.
Moreover, under the current global shortage of vaccines, it should be remarked
that immunization is enhanced when mobility is reduced. Additionally,
repetition of vaccination campaigns should be timed considering the immunity
lapse of the vaccinated (and recovered) people. Finally, the model is extended
to include the effect of isolation of detected positive cases, shown to be
important to reduce infections.


  Policymakers often choose a policy bundle that is a combination of different
interventions in different dosages. We develop a new technique -- treatment
variant aggregation (TVA) -- to select a policy from a large factorial design.
TVA pools together policy variants that are not meaningfully different and
prunes those deemed ineffective. This allows us to restrict attention to
aggregated policy variants, consistently estimate their effects on the outcome,
and estimate the best policy effect adjusting for the winner's curse. We apply
TVA to a large randomized controlled trial that tests interventions to
stimulate demand for immunization in Haryana, India. The policies under
consideration include reminders, incentives, and local ambassadors for
community mobilization. Cross-randomizing these interventions, with different
dosages or types of each intervention, yields 75 combinations. The policy with
the largest impact (which combines incentives, ambassadors who are information
hubs, and reminders) increases the number of immunizations by 44% relative to
the status quo. The most cost-effective policy (information hubs, ambassadors,
and SMS reminders but no incentives) increases the number of immunizations per
dollar by 9.1% relative to status quo.


  We studied the COVID-19 pandemic evolution in selected African countries. For
each country considered, we modeled simultaneously the data of the active,
recovered and death cases. In this study, we used a year of data since the
first cases were reported. We estimated the time-dependent basic reproduction
numbers, $R_0$, and the fractions of infected but unaffected populations, to
offer insights into containment and vaccine strategies in African countries. We
found that $R_0\leq 4$ at the start of the pandemic but has since fallen to
$R_0 \sim 1$. The unaffected fractions of the populations studied vary between
$1-10$\% of the recovered cases.


  Microbial communities are ubiquitous in nature and come in a multitude of
forms, ranging from communities dominated by a handful of species to
communities containing a wide variety of metabolically distinct organisms. This
huge range in diversity is not a curiosity - microbial diversity has been
linked to outcomes of substantial ecological and medical importance. However,
the mechanisms underlying microbial diversity are still under debate, as simple
mathematical models only permit as many species to coexist as there are
resources. A plethora of mechanisms have been proposed to explain the origins
of microbial diversity, but many of these analyses omit a key property of real
microbial ecosystems: the propensity of the microbes themselves to change their
growth properties within and across generations. In order to explore the impact
of this key property on microbial diversity, we expand upon a recently
developed model of microbial diversity in fluctuating environments. We
implement changes in growth strategy in two distinct ways. First, we consider
the regulation of a cell's enzyme levels within short, ecological times, and
second we consider evolutionary changes driven by mutations across generations.
Interestingly, we find that these two types of microbial responses to the
environment can have drastically different outcomes. Enzyme regulation may
collapse diversity over long enough times while, conversely,
strategy-randomizing mutations can produce a "rich-get-poorer" effect that
promotes diversity. This work makes explicit, using a simple serial-dilutions
framework, the conflicting ways that microbial adaptation and evolution can
affect community diversity.


  In this work Data Envelopment Analysis (DEA) is employed in thirty-day
windows to quantify temporal evolution of relative pandemic mitigation
inefficiency of Brazilian municipalities. For each thirty-day window the
results of inefficiency scores of over five thousand Brazilian municipalities
are displayed on maps, to address the spatial distribution of the corresponding
values. This phenomenological spatiotemporal approach reveals location of the
hotspots, in terms of relative pandemic inefficiency of the municipalities, at
different points in time along the pandemic. It is expected that the current
approach may subsidize decision making through comparative analysis of previous
practice, and thus contribute to the future pandemic mitigation efforts.


  Evolutionary transitions occur when previously-independent replicating
entities unite to form more complex individuals. Such transitions have
profoundly shaped natural evolutionary history and occur in two forms:
fraternal transitions involve lower-level entities that are kin (e.g.,
transitions to multicellularity or to eusocial colonies), while egalitarian
transitions involve unrelated individuals (e.g., the origins of mitochondria).
The necessary conditions and evolutionary mechanisms for these transitions to
arise continue to be fruitful targets of scientific interest. Here, we examine
a range of fraternal transitions in populations of open-ended self-replicating
computer programs. These digital cells were allowed to form and replicate kin
groups by selectively adjoining or expelling daughter cells. The capability to
recognize kin-group membership enabled preferential communication and
cooperation between cells. We repeatedly observed group-level traits that are
characteristic of a fraternal transition. These included reproductive division
of labor, resource sharing within kin groups, resource investment in offspring
groups, asymmetrical behaviors mediated by messaging, morphological patterning,
and adaptive apoptosis. We report eight case studies from replicates where
transitions occurred and explore the diverse range of adaptive evolved
multicellular strategies.


  In this research, we study biological networks at different scales: a gene
autoregulatory network at the single-cell level and the gut microbiota at the
population level.
  Proteins are the main actors in cells, they are the building blocks, act as
enzymes and antibodies. The production of proteins is mediated by transcription
factors. In some cases, a protein acts as its own transcription factor, this is
called autoregulation. It is known that autorepression speeds up the response
and that autoactivation can lead to multiple stable equilibria. In this thesis,
we study the effects of the combination of activation and repression in
autoregulation, as a case study we investigate the possible dynamics of the
leucine responsive protein B of the archaeon Sulfolobus solfataricus (Ss-LrpB),
a protein that regulates itself in a unique and non-monotonic way via three
binding boxes. We examine for which conditions this type of network leads to
oscillations or bistability.
  In the second part, much larger biological systems are considered. Ecological
systems, among which the human gut microbiome, are characterized by
heavy-tailed abundance profiles. We study how these distributions can arise
from population-based models by adding saturation effects and linear noise.
Moreover, we examine different characteristics of experimental time series of
microbial communities, such as the noise color and neutrality of the
biodiversity, and look at the influence of the parameters on these
characteristics. With the first research topic we want to lay a foundation for
the understanding of non-monotonic gene regulation and take the first steps
toward synthetic biology in archaea. In the second part of the thesis, we
investigate experimental time series from complex ecosystems and seek
theoretical models reproducing all observed characteristics in view of building
predictive models.


  In a complex community, species continuously adapt to each other. On rare
occasions, the adaptation of a species can lead to the extinction of others,
and even its own. "Adaptive dynamics" is the standard mathematical framework to
describe evolutionary changes in community interactions, and in particular,
predict adaptation driven extinction. Unfortunately, most authors implement the
equations of adaptive dynamics through computer simulations, that require
assuming a large number of questionable parameters and fitness functions. In
this study we present analytical solutions to adaptive dynamics equations,
thereby clarifying how outcomes depend on any computational input. We develop
general formulas that predict equilibrium abundances over evolutionary time
scales. Additionally, we predict which species will go extinct next, and when
this will happen.


  Languages emerge and change over time at the population level though
interactions between individual speakers. It is, however, hard to directly
observe how a single speaker's linguistic innovation precipitates a
population-wide change in the language, and many theoretical proposals exist.
We introduce a very general mathematical model that encompasses a wide variety
of individual-level linguistic behaviours and provides statistical predictions
for the population-level changes that result from them. This model allows us to
compare the likelihood of empirically-attested changes in definite and
indefinite articles in multiple languages under different assumptions on the
way in which individuals learn and use language. We find that accounts of
language change that appeal primarily to errors in childhood language
acquisition are very weakly supported by the historical data, whereas those
that allow speakers to change incrementally across the lifespan are more
plausible, particularly when combined with social network effects.


  The results of Kermack-McKendrick SIR model are planned to be reproduced by
cellular automata (CA) lattice model. The CA algorithms are proposed to study
the model of epidemic, systematically. The basic goal is to capture the effects
of spreading of infection over a scale of length. This CA model can provide the
rate of growth of the infection over the space which was lacking in the
mean-field like SIR model. The motion of the circular front of an infected
cluster shows a linear behaviour in time. The correlation of a particular site
to be infected, with respect to the central site is also studied. The outcomes
of CA model are in good agreement with that obtained from SIR model. The
results of vaccination have been also incorporated in the CA algorithm with a
satisfactory degree of success. The advantage of the present model is that it
can shed considerable amount of light to the physical properties of the spread
of a typical epidemic in a simple, yet robust way.


  Understanding the biological function of knots in proteins and their folding
process is an open and challenging question in biology. Recent studies classify
the topology and geometry of knotted proteins by analysing the distribution of
a protein's planar projections using topological objects called knotoids. We
approach the analysis of proteins with the same topology by introducing a
topologically inspired statistical metric between their knotoid distributions.
We detect geometric differences between trefoil proteins by characterising
their entanglement and we recover a clustering by sequence similarity. By
looking directly at the geometry and topology of their native states, we are
able to probe different folding pathways for proteins forming open-ended
trefoil knots. Interestingly, our pipeline reveals that the folding pathway of
shallow knotted Carbonic Anhydrases involves the creation of a double-looped
structure, differently from what was previously observed for deeply knotted
trefoil proteins. We validate this with Molecular Dynamics simulations.


  We consider an ergodic harvesting problem with model ambiguity that arises
from biology. To account for the ambiguity, the problem is constructed as a
stochastic game with two players: the decision-maker (DM) chooses the `best'
harvesting policy and an adverse player chooses the `worst' probability
measure. The main result is establishing an optimal strategy (also referred to
as a control) of the DM and showing that it is a threshold policy. The optimal
threshold and the optimal payoff are obtained by solving a free-boundary
problem emerging from the Hamilton--Jacobi--Bellman (HJB) equation. As part of
the proof, we fix a gap that appeared in the HJB analysis of [Alvarez and
Hening, {\em Stochastic Process. Appl.}, 2019], a paper that analyzed the
risk-neutral version of the ergodic harvesting problem. Finally, we study the
dependence of the optimal threshold and the optimal payoff on the ambiguity
parameter and show that if the ambiguity goes to 0, the problem converges, to
the risk-neutral problem.


  In epidemiology, the effective reproduction number $R_e$ is used to
characterize the growth rate of an epidemic outbreak. In this paper, we
investigate properties of $R_e$ for a modified SEIR model of COVID-19 in the
city of Houston, TX USA, in which the population is divided into low-risk and
high-risk subpopulations. The response of $R_e$ to two types of control
measures (testing and distancing) applied to the two different subpopulations
is characterized. A nonlinear cost model is used for control measures, to
include the effects of diminishing returns. We propose three types of heuristic
strategies for mitigating COVID-19 that are targeted at reducing $R_e$, and we
exhibit the tradeoffs between strategy implementation costs and number of
deaths. We also consider two variants of each type of strategy: basic
strategies, which consider only the effects of controls on $R_e$, without
regard to subpopulation; and high-risk prioritizing strategies, which maximize
control of the high-risk subpopulation. Results showed that of the three
heuristic strategy types, the most cost-effective involved setting a target
value for $R_e$ and applying sufficient controls to attain that target value.
This heuristic led to strategies that begin with strict distancing of the
entire population, later followed by increased testing. Strategies that
maximize control on high-risk individuals were less cost-effective than basic
strategies that emphasize reduction of the rate of spreading of the disease.
The model shows that delaying the start of control measures past a certain
point greatly worsens strategy outcomes. We conclude that the effective
reproduction can be a valuable real-time indicator in determining
cost-effective control strategies.


  We present artificial neural networks as a feasible replacement for a
mechanistic model of mosquito abundance. We develop a feed-forward neural
network, a long short-term memory recurrent neural network, and a gated
recurrent unit network. We evaluate the networks in their ability to replicate
the spatiotemporal features of mosquito populations predicted by the
mechanistic model, and discuss how augmenting the training data with time
series that emphasize specific dynamical behaviors affects model performance.
We conclude with an outlook on how such equation-free models may facilitate
vector control or the estimation of disease risk at arbitrary spatial scales.


  The dynamics of the spread of contagions such as viruses, infectious diseases
or even rumors/opinions over contact networks (graphs) have effectively been
captured by the well known \textit{Susceptible-Infected-Susceptible} ($SIS$)
epidemic model in recent years. When it comes to competition between two such
contagions spreading on overlaid graphs, their propagation is captured by
so-called \textit{bi-virus} epidemic models. Analysis of such dynamical systems
involve the identification of equilibrium points and its convergence
properties, which determine whether either of the viruses dies out, or both
survive together. We demonstrate how the existing works are unsuccessful in
characterizing a large subset of the model parameter space, including all
parameters for which the competitiveness of the bi-virus system is significant
enough to attain coexistence of the epidemics. In this paper, we fill in this
void and obtain convergence results for the entirety of the model parameter
space; giving precise conditions (necessary and sufficient) under which the
system \textit{globally converges} to a \textit{trichotomy} of possible
outcomes: a virus-free state, a single-virus state, and to a coexistence state
-- the first such result.


  We study the effects of physical distancing measures for the spread of
COVID-19 in regional areas within British Columbia, using the reported cases of
the five provincial Health Authorities. Building on the Bayesian
epidemiological model of Anderson et al. (2020), we propose a hierarchical
regional Bayesian model with time-varying regional parameters between March to
December of 2020. In the absence of COVID-19 variants and vaccinations during
this period, we examine the regionalized basic reproduction number, modelled
prevalence, relative reduction in contact due to physical distancing, and
proportion of anticipated cases that have been tested and reported. We observe
significant differences between the regional and provincial-wide models and
demonstrate the hierarchical regional model can better estimate regional
prevalence, especially in rural regions. These results indicate that it can be
useful to apply similar regional models to other parts of Canada or other
countries.


  Given overlapping subsets of a set of taxa (e.g. species), and posterior
distributions on phylogenetic tree topologies for each of these taxon sets, how
can we infer a posterior distribution on phylogenetic tree topologies for the
entire taxon set? Although the equivalent problem for in the non-Bayesian case
has attracted substantial research, the Bayesian case has not attracted the
attention it deserves. In this paper we develop a variational Bayes approach to
this problem and demonstrate its effectiveness.


  Ion channels are important proteins for physiological information transfer
and functional control. To predict the microscopic origins of their
voltage-conductance characteristics, we here applied dissipation-corrected
targeted Molecular Dynamics in combination with Langevin equation simulations
to potassium diffusion through the Gramicidin A channel as a test system.
Performing a non-equilibrium principal component analysis on backbone dihedral
angles, we find coupled protein-ion dynamics to occur during ion transfer. The
dissipation-corrected free energy profiles correspond well to predictions from
other biased simulation methods. The incorporation of an external electric
field in Langevin simulations enables the prediction of macroscopic observables
in the form of I-V characteristics.


  We introduce an evolutionary game on hypergraphs in which decisions between a
risky alternative and a safe one are taken in social groups of different sizes.
The model naturally reproduces choice shifts, namely the differences between
the preference of individual decision makers and the consensual choice of a
group, that have been empirically observed in choice dilemmas. In particular, a
deviation from the Nash equilibrium towards the risky strategy occurs when the
dynamics takes place on heterogeneous hypergraphs. These results can explain
the emergence of irrational herding and radical behaviours in social groups.


  Group testing allows saving chemical reagents, analysis time, and costs, by
testing pools of samples instead of individual samples. We introduce a class of
group testing protocols with small dilution, suited to operate even at high
prevalence ($5\%-10\%$), and maximizing the fraction of samples classified
positive/negative within the first round of tests. Precisely, if the tested
group has exactly one positive sample then the protocols identify it without
further individual tests. The protocols also detect the presence of two or more
positives in the group, in which case a second round could be applied to
identify the positive individuals. With a prevalence of $5\%$ and maximum
dilution 6, with 100 tests we classify 242 individuals, $92\%$ of them in one
round and $8\%$ requiring a second individual test. In comparison, the
Dorfman's scheme can test 229 individuals with 100 tests, with a second round
for $18.5\%$ of the individuals.


  In the case of SARS-CoV-2 pandemic management, wastewater-based epidemiology
aims to derive information on the infection dynamics by monitoring virus
concentrations in the wastewater. However, due to the intrinsic random
fluctuations of the viral signal in the wastewater (due to e.g., dilution;
transport and fate processes in sewer system; variation in the number of
persons discharging; variations in virus excretion and water consumption per
day) the subsequent prevalence analysis may result in misleading conclusions.
It is thus helpful to apply data filtering techniques to reduce the noise in
the signal. In this paper we investigate 13 smoothing algorithms applied to the
virus signals monitored in four wastewater treatment plants in Austria. The
parameters of the algorithms have been defined by an optimization procedure
aiming for performance metrics. The results are further investigated by means
of a cluster analysis. While all algorithms are in principle applicable,
SPLINE, Generalized Additive Model and Friedman Super Smoother are recognized
as superior methods in this context (with the latter two having a tendency to
over-smoothing). A first analysis of the resulting datasets indicates the
influence of catchment size for wastewater-based epidemiology as smaller
communities both reveal a signal threshold before any relation with infection
dynamics is visible and also a higher sensitivity towards infection clusters.


  The oxDNA model of DNA has been applied widely to systems in biology,
biophysics and nanotechnology. It is currently available via two independent
open source packages. Here we present a set of clearly-documented exemplar
simulations that simultaneously provide both an introduction to simulating the
model, and a review of the model's fundamental properties. We outline how
simulation results can be interpreted in terms of -- and feed into our
understanding of -- less detailed models that operate at larger length scales,
and provide guidance on whether simulating a system with oxDNA is worthwhile.


  Metabolic modeling strives to develop modeling approaches that are robust and
highly predictive. To achieve this, various modeling designs, including hybrid
models, and parameter estimation methods that define the type and number of
parameters used in the model, are adapted. Accurate input data play an
important role so that the selection of experimental methods that provide input
data of the required precision with low measurement errors is crucial. For the
biopharmaceutically relevant protein glycosylation, the most prominent
available models are kinetic models which are able to capture the dynamic
nature of protein N-glycosylation. In this review we focus on how to choose the
most suitable model for a specific research question, as well as on parameters
and considerations to take into account before planning relevant experiments.


  With the invention of the COVID-19 vaccine, shipping and distributing are
crucial in controlling the pandemic. In this paper, we build a mean-field
variational problem in a spatial domain, which controls the propagation of
pandemic by the optimal transportation strategy of vaccine distribution. Here
we integrate the vaccine distribution into the mean-field SIR model designed in
our previous paper arXiv:2006.01249. Numerical examples demonstrate that the
proposed model provides practical strategies in vaccine distribution on a
spatial domain.


  We consider the recent surge of information on the potential benefits of
acid-suppression drugs in the context of COVID-19, with an eye on the
variability (and confusion) across the reported findings--at least as regards
the popular antacid famotidine. The inconsistencies reflect contradictory
conclusions from independent clinical-based studies that took roughly similar
approaches, in terms of experimental design (retrospective, cohort-based, etc.)
and statistical analyses (propensity-score matching and stratification, etc.).
The confusion has significant ramifications in choosing therapeutic
interventions: e.g., do potential benefits of famotidine indicate its use in a
particular COVID-19 case? Beyond this pressing therapeutic issue, conflicting
information on famotidine must be resolved before its integration in
ontological and knowledge graph-based frameworks, which in turn are useful in
drug repurposing efforts. To begin systematically structuring the rapidly
accumulating information, in the hopes of clarifying and reconciling the
discrepancies, we consider the contradictory information along three proposed
'axes': (1) a context-of-disease axis, (2) a degree-of-[therapeutic]-benefit
axis, and (3) a mechanism-of-action axis. We suspect that incongruencies in how
these axes have been (implicitly) treated in past studies has led to the
contradictory indications for famotidine and COVID-19. We also trace the
evolution of information on acid-suppression agents as regards the
transmission, severity, and mortality of COVID-19, given the many literature
reports that have accumulated. By grouping the studies conceptually and
thematically, we identify three eras in the progression of our understanding of
famotidine and COVID-19. Harmonizing these findings is a key goal for both
clinical standards-of-care (COVID and beyond) as well as ontological and
knowledge graph-based approaches.


  We investigated the data-driven relationship between features in the tumor
microenvironment (TME) and the overall and 5-year survival in triple-negative
breast cancer (TNBC) and non-TNBC (NTNBC) patients by using Explainable
Artificial Intelligence (XAI) models. We used clinical information from
patients with invasive breast carcinoma from The Cancer Genome Atlas and from
two studies from the cbioPortal, the PanCanAtlas project and the GDAC Firehose
study. In this study, we used a normalized RNA sequencing data-driven cohort
from 1,015 breast cancer patients, alive or deceased, from the UCSC Xena data
set and performed integrated deconvolution with the EPIC method to estimate the
percentage of seven different immune and stromal cells from RNA sequencing
data. Novel insights derived from our XAI model showed that CD4+ T cells and B
cells are more critical than other TME features for enhanced prognosis for both
TNBC and NTNBC patients. Our XAI model revealed the critical inflection points
(i.e., threshold fractions) of CD4+ T cells and B cells above or below which
5-year survival rates improve. Subsequently, we ascertained the conditional
probabilities of $\geq$ 5-year survival in both TNBC and NTNBC patients under
specific conditions inferred from the inflection points. In particular, the XAI
models revealed that a B-cell fraction exceeding 0.018 in the TME could ensure
100% 5-year survival for NTNBC patients. The findings from this research could
lead to more accurate clinical predictions and enhanced immunotherapies and to
the design of innovative strategies to reprogram the TME of breast cancer
patients.


  Fine-grain patterns produced by juxtacrine signalling, have been studied
using static monolayers as cellular domains. Unfortunately, analytical results
are restricted to a few cells due to the algebraic complexity of nonlinear
dynamical systems. Motivated by concentric patterning of Notch expression
observed in the mammary gland, we combine concepts from graph and control
theory to represent cellular connectivity. The resulting theoretical framework
allows us to exploit the symmetry of multicellular bilayer structures in 2D and
3D, thereby deriving analytical conditions that drive the dynamical system to
form laminar patterns consistent with the formulation of cell polarity.
Critically, the conditions are independent of the precise dynamical details,
thus the framework allows for the utmost generality in understanding the
influence of cellular geometry on patterning in lateral-inhibition systems.
Applying the analytic conditions to mammary organoids suggests that intense
cell signalling polarity is required for the maintenance of stratified
cell-types within a static bilayer using a lateral-inhibition mechanism.
Furthermore, by employing 2D and 3D cell-based models, we highlight that the
cellular polarity conditions derived from static domains have the capacity to
generate laminar patterning in dynamic environments. However, they are
insufficient for the maintenance of patterning when subjected to substantial
morphological perturbations. In agreement with the mathematical implications of
strict signalling polarity induced on the cells, we propose an adhesion
dependent Notch-Delta biological process which has the potential to initiate
bilayer stratification in a developing mammary organoid.


  Mutation and drift play opposite roles in genetics. While mutation creates
diversity, drift can cause gene variants to disappear, especially when they are
rare. In the absence of natural selection and migration, the balance between
the drift and mutation in a well-mixed population defines its diversity. The
Moran model captures the effects of these two evolutionary forces and has a
counterpart in social dynamics, known as the Voter model with external opinion
influencers. Two extreme outcomes of the Voter model dynamics are consensus and
coexistence of opinions, which correspond to low and high diversity in the
Moran model. Here we use a Shannon's information-theoretic approach to
characterize the smooth transition between the states of consensus and
coexistence of opinions in the Voter model. Mapping the Moran into the Voter
model we extend the results to the mutation-drift balance and characterize the
transition between low and high diversity in finite populations. Describing the
population as a network of connected individuals we show that the transition
between the two regimes depends on the network topology of the population and
on the possible asymmetries in the mutation rates.


  This work presents a mathematical model which describes both the genesis and
growth of oxygenic photogranules (OPGs) and the related treatment process. The
photogranule has been modelled as a free boundary domain with radial symmetry,
which evolves over time as a result of microbial growth, attachment and
detachment processes. A system of hyperbolic and parabolic PDEs has been
considered to model the advective transport and growth of sessile biomass and
the diffusive transport and conversion of soluble substrates. The reactor has
been modelled as a sequencing batch reactor (SBR) through a system of first
order IDEs. Phototrophic biomass has been considered for the first time in
granular biofilms, and cyanobacteria and microalgae are taken into account
separately to model their differences in growth rate and light harvesting and
utilization. To describe the key role of cyanobacteria in the photogranules
formation process, the attachment velocity of all suspended microbial species
has been modelled as a function of the cyanobacteria concentration in suspended
form. The model takes into account the main biological aspects and processes
involved in OPGs based systems: heterotrophic and photoautotrophic activities
of cyanobacteria and microalgae, metabolic activity of heterotrophic and
nitrifying bacteria, microbial decay, EPS secrection, diffusion and conversion
of soluble substrates (inorganic and organic carbon, ammonia, nitrate and
oxygen), symbiotic and competitive interactions between the different microbial
species, day-night cycle, light diffusion and attenuation across the granular
biofilm and photoinhibion phenomena. The model has been integrated numerically,
investigating the evolution and microbial composition of photogranules and the
treatment efficiency of the OPGs-based system. The results show the consistency
of the model and confirm the effectiveness of the OPGs technology.


  Deep learning has empowered analysis for single-cell sequencing data in many
ways and has generated deep understanding about a range of complex cellular
systems. As the booming single-cell sequencing technologies brings the surge of
high dimensional data that come from different sources and represent cellular
systems with different features, there is an equivalent rise and challenge of
integrating single-cell sequence across modalities. Here, we present a novel
adversarial approach to integrate single-cell chromatin accessibility and gene
expression data in a semi-supervised manner. We demonstrate that our method
substantially improves data integration from a simple adversarial domain
adaption approach, and it also outperforms two state-of-the-art (SOTA) methods.


  In phylogenetics, it is important for the phylogenetic network model
parameters to be identifiable so that the evolutionary histories of a group of
species can be consistently inferred. However, as the complexity of the
phylogenetic network models grows, the identifiability of network models
becomes increasingly difficult to analyze. As an attempt to analyze the
identifiability of network models, we check whether two networks are
distinguishable. In this paper, we specifically study the distinguishability of
phylogenetic network models associated with level-2 networks. Using an
algebraic approach, namely using discrete Fourier transformation, we present
some results on the distinguishability of some level-2 networks, which
generalize earlier work on the distinguishability of level-1 networks. In
particular, we study simple and semisimple level-2 networks. Simple and
semisimple level-2 networks can be thought as generalizations of level-1 sunlet
and cycle networks, respectively. Moreover, we also compare the varieties
associated with semisimple level-2 and cycle networks.


  The rapid emergence and the disastrous impact of the severe acute respiratory
syndrome coronavirus 2 (SARS-CoV-2) pandemic on public health, societies, and
economies around the world has created an urgent need for understanding the
pathways critical for virus transmission. Airborne virus transmission by
asymptomatic SARS-CoV-2-infected individuals is considered to be a major
contributor to the spread of SARS-CoV-2 and social distancing and wearing of
face masks in public have been implemented as countermeasures in many
countries. However, a comprehensive risk assessment framework for the airborne
transmission of SARS-CoV-2 incorporating realistic assumptions on the
filtration of infectious aerosols (IAs) by face masks is not available yet. In
particular, in most end-to-end models for airborne virus transmission, it is
neglected that the stochastic spread of IAs through imperfectly fitted face
masks depends on the dynamics of the breathing of the wearer. In this paper, we
consider airborne virus transmission from an infected but asymptomatic person
to a healthy person, both wearing imperfectly fitted face masks, in an indoor
environment. By framing the end-to-end virus transmission as a Molecular
Communications (MC) system, we obtain a statistical description of the number
of IAs inhaled by the healthy person subject to the respective configurations
of the face masks of both persons. We demonstrate that the exhalation and
inhalation air flow dynamics have a significant impact on the stochastic
filtering of IAs by the face masks. Furthermore, we show that the fit of the
face mask of the infected person can highly impact the infection probability.
We conclude that the proposed MC model may contribute a valuable assessment
tool to fight the spread of SARS-CoV-2 as it encompasses the randomness of the
transmission process and enables comprehensive risk analysis beyond statistical
averages.


  In Kauffman-s NK model, myopic local search involves flipping one
randomly-chosen bit of an N-bit decision string in every time step and
accepting the new configuration if that has higher fitness. One issue is that,
this algorithm consumes the full extent of computational resources allocated -
given by the number of alternative configurations inspected - even though
search is expected to terminate the moment there are no neighbors having higher
fitness. Otherwise, the algorithm must compute the fitness of all N neighbors
in every time step, consuming a high amount of resources. In order to get
around this problem, I describe an algorithm that allows search to logically
terminate relatively early, without having to evaluate fitness of all N
neighbors at every time step. I further suggest that when the efficacy of two
algorithms need to be compared head to head, imposing a common limit on the
number of alternatives evaluated - metering - provides the necessary level
field.


  Understanding influencing factors is essential for the surveillance and
prevention of infectious diseases, and the factors are likely to vary spatially
and temporally as the disease progresses. Taking daily cases and deaths data
during the coronavirus disease 2019 (COVID-19) outbreak in the U.S. as a case
study, we develop a mobility-augmented geographically and temporally weighted
regression (M-GTWR) model to quantify the spatiotemporal impacts of
social-demographic factors and human activities on the COVID-19 dynamics.
Different from the base GTWR model, we incorporate a mobility-adjusted distance
weight matrix where travel mobility is used in addition to the spatial
adjacency to capture the correlations among local observations. The model
residuals suggest that the proposed model achieves a substantial improvement
over other benchmark methods in addressing the spatiotemporal nonstationarity.
Our results reveal that the impacts of social-demographic and human activity
variables present significant spatiotemporal heterogeneity. In particular, a 1%
increase in population density may lead to 0.63% and 0.71% more daily cases and
deaths, and a 1% increase in the mean commuting time may result in 0.22% and
0.95% increases in daily cases and deaths. Although increased human activities
will, in general, intensify the disease outbreak, we report that the effects of
grocery and pharmacy-related activities are insignificant in areas with high
population density. And activities at the workplace and public transit are
found to either increase or decrease the number of cases and deaths, depending
on particular locations. The results of our study establish a quantitative
framework for identifying influencing factors during a disease outbreak, and
the obtained insights may have significant implications in guiding the
policy-making against infectious diseases.


  Kinase-targeted drug design is challenging. It requires designing inhibitors
that can bind to specific kinases when all kinase catalytic domains share a
common folding scaffold that binds ATP. Thus, obtaining the desired
selectivity, given the whole human kinome, is a fundamental task during
early-stage drug discovery. This begins with deciphering the kinase-ligand
characteristics, analyzing the structure-activity relationships, and
prioritizing the desired drug molecules across the whole kinome. Currently,
there are more than 300 kinases with released PDB structures, which provides a
substantial structural basis to gain these necessary insights. Here, we review
in silico structure-based methods - notably, a function-site interaction
fingerprint approach used in exploring the complete human kinome. In silico
methods can be explored synergistically with multiple cell-based or
protein-based assay platforms such as KINOMEscan. We conclude with new drug
discovery opportunities associated with kinase signaling networks and using
machine/deep learning techniques broadly referred to as structural biomedical
data science.


  Predicting the final folded structure of protein molecules and simulating
their folding pathways is of crucial importance for designing viral drugs and
studying diseases such as Alzheimer's at the molecular level. To this end, this
paper investigates the problem of protein conformation prediction under the
constraint of avoiding high-entropy-loss routes during folding. Using the
well-established kinetostatic compliance (KCM)-based nonlinear dynamics of a
protein molecule, this paper formulates the protein conformation prediction as
a pointwise optimal control synthesis problem cast as a quadratic program (QP).
It is shown that the KCM torques in the protein folding literature can be
utilized for defining a reference vector field for the QP-based control
generation problem. The resulting kinetostatic control torque inputs will be
close to the KCM-based reference vector field and guaranteed to be constrained
by a predetermined bound; hence, high-entropy-loss routes during folding are
avoided while the energy of the molecule is decreased.


  This paper dwells on certain novel game-theoretic investigations in
bio-mimicry, discussed from the perspectives of information asymmetry,
individual utility and its optimization via strategic interactions involving
co-evolving preys (e.g., insects) and predators (e.g., reptiles) who learn.
Formally, we consider a panmictic ecosystem, occupied by species of prey with
relatively short lifespan, which evolve mimicry signals over generations as
they interact with predators with relatively longer lifespans, thus endowing
predators with the ability to learn prey signals. Every prey sends a signal and
provides utility to the predator. The prey can be either nutritious or toxic to
the predator, but the prey may signal (possibly) deceptively without revealing
its true "type." The model is used to study the situation where multi-armed
bandit predators with zero prior information are introduced into the ecosystem.
As a result of exploration and exploitation the predators naturally select the
prey that result in the evolution of those signals. This co-evolution of
strategies produces a variety of interesting phenomena which are subjects of
this paper.


  A recent technology breakthrough in spatial molecular profiling has enabled
the comprehensive molecular characterizations of single cells while preserving
spatial information. It provides new opportunities to delineate how cells from
different origins form tissues with distinctive structures and functions. One
immediate question in spatial molecular profiling data analysis is to identify
genes whose expressions exhibit spatially correlated patterns, called spatially
variable genes. Most current methods to identify spatially variable genes are
built upon the geostatistical model with Gaussian process to capture the
spatial patterns, which rely on ad hoc kernels that could limit the models'
ability to identify complex spatial patterns. In order to overcome this
challenge and capture more types of spatial patterns, we introduce a Bayesian
approach to identify spatially variable genes via a modified Ising model. The
key idea is to use the energy interaction parameter of the Ising model to
characterize spatial expression patterns. We use auxiliary variable Markov
chain Monte Carlo algorithms to sample from the posterior distribution with an
intractable normalizing constant in the model. Simulation studies using both
simulated and synthetic data showed that the energy-based modeling approach led
to higher accuracy in detecting spatially variable genes than those
kernel-based methods. When applied to two real spatial transcriptomics
datasets, the proposed method discovered novel spatial patterns that shed light
on the biological mechanisms. In summary, the proposed method presents a new
perspective for analyzing spatial transcriptomics data.


  More than any other infectious disease epidemic, the COVID-19 pandemic has
been characterized by the generation of large volumes of viral genomic data at
an incredible pace due to recent advances in high-throughput sequencing
technologies, the rapid global spread of SARS-CoV-2, and its persistent threat
to public health. However, distinguishing the most epidemiologically relevant
information encoded in these vast amounts of data requires substantial effort
across the research and public health communities. Studies of SARS-CoV-2
genomes have been critical in tracking the spread of variants and understanding
its epidemic dynamics, and may prove crucial for controlling future epidemics
and alleviating significant public health burdens. Together, genomic data and
bioinformatics methods enable broad-scale investigations of the spread of
SARS-CoV-2 at the local, national, and global scales and allow researchers the
ability to efficiently track the emergence of novel variants, reconstruct
epidemic dynamics, and provide important insights into drug and vaccine
development and disease control. Here, we discuss the tremendous opportunities
that genomics offers to unlock the effective use of SARS-CoV-2 genomic data for
efficient public health surveillance and guiding timely responses to COVID-19.


  The cellular milieu is teeming with biochemical nano-machines whose activity
is a strong source of correlated non-thermal fluctuations termed active noise.
Essential elements of this circuitry are enzymes, catalysts that speed up the
rate of metabolic reactions by orders of magnitude, thereby making life
possible. Here, we examine the possibility that active noise in the cell, or in
vitro, affects enzymatic catalytic rate by accelerating or decelerating the
crossing rate of energy barriers during the reaction. Considering hydrodynamic
perturbations induced by biochemical activity as a source of active noise, we
evaluate their impact on the enzymatic cycle using a combination of analytic
and numerical methods. Our estimates show that the fast component of the active
noise spectrum enhances the rate of enzymes, while reactions remain practically
unaffected by the slow noise spectrum. Revisiting the physics of barrier
crossing under the influence of active hydrodynamic fluctuations suggests that
the biochemical activity of macromolecules such as enzymes is coupled to active
noise. Thus, we propose that enzymatic catalysis is a collective, many-body
process in which enzymes may affect each other's activity via long-range
hydrodynamic interaction, with potential impact on biochemical networks in
living and artificial systems alike.


  Kinetic models of biochemical systems used in the modern literature often
contain hundreds or even thousands of variables. While these models are
convenient for detailed simulations, their size is often an obstacle to
deriving mechanistic insights. One way to address this issue is to perform an
exact model reduction by finding a self-consistent lower-dimensional projection
of the corresponding dynamical system. Recently, a new algorithm CLUE has been
designed and implemented, which allows one to construct an exact linear
reduction of the smallest possible dimension such that the fixed variables of
interest are preserved. It turned out that allowing arbitrary linear
combinations (as opposed to zero-one combinations used in the prior approaches)
may yield a much smaller reduction. However, there was a drawback: some of the
new variables did not have clear physical meaning, thus making the reduced
model harder to interpret. We design and implement an algorithm that, given an
exact linear reduction, re-parametrizes it by performing an invertible
transformation of the new coordinates to improve the interpretability of the
new variables. We apply our algorithm to three case studies and show that
"uninterpretable" variables disappear entirely in all the case studies. The
implementation of the algorithm and the files for the case studies are
available at https://github.com/xjzhaang/LumpingPostiviser.


  We introduce and parameterize a chemomechanical model of microtubule dynamics
on the dimer level, which is based on the allosteric tubulin model and includes
attachment, detachment and hydrolysis of tubulin dimers as well as stretching
of lateral bonds, bending at longitudinal junctions, and the possibility of
lateral bond rupture and formation. The model is computationally efficient such
that we reach sufficiently long simulation times to observe repeated
catastrophe and rescue events at realistic tubulin concentrations and
hydrolysis rates, which allows us to deduce catastrophe and rescue rates. The
chemomechanical model also allows us to gain insight into microscopic features
of the GTP-tubulin cap structure and microscopic structural features triggering
microtubule catastrophes and rescues. Dilution simulations show qualitative
agreement with experiments. We also explore the consequences of a possible
feedback of mechanical forces onto the hydrolysis process and the GTP-tubulin
cap structure.


  In this letter, we investigate the population dynamics in a May-Leonard
formulation of the rock-paper-scissors game in which one or two species, which
we shall refer to as "weak", have a reduced predation or reproduction
probability. We show that in a nonspatial model the stationary solution where
all three species coexist is always unstable, while in a spatial stochastic
model coexistence is possible for a wide parameter space. We find, that a
reduced predation probability results in a significantly higher abundance of
"weak" species, in models with either one or two "weak" species, as long as the
simulation lattices are sufficiently large for coexistence to prevail. On the
other hand, we show that a reduced reproduction probability has a smaller
impact on the abundance of "weak" species, generally leading to a slight
decrease of its population size -- the increase of the population size of one
of the "weak" species being more than compensated by the reduction of the
other, in the two species case. We further show that the species abundances in
models where both predation and reproduction probabilities are simultaneously
reduced may be accurately estimated from the results obtained considering only
a reduction of either the predation or the reproduction probability.