  Just as in eukaryotes, high-throughput chromosome conformation capture (Hi-C)
data have revealed nested organizations of bacterial chromosomes into
overlapping interaction domains. In this chapter, we present a multiscale
analysis framework aiming at capturing and quantifying these properties. These
include both standard tools (e.g. contact laws) and novel ones such as an index
that allows identifying loci involved in domain formation independently of the
structuring scale at play. Our objective is two-fold. On the one hand, we aim
at providing a full, understandable Python/Jupyter-based code which can be used
by both computer scientists as well as biologists with no advanced
computational background. On the other hand, we discuss statistical issues
inherent to Hi-C data analysis, focusing more particularly on how to properly
assess the statistical significance of results. As a pedagogical example, we
analyze data produced in {\it Pseudomonas aeruginosa}, a model pathogenetic
bacterium. All files (codes and input data) can be found on a github
repository. We have also embedded the files into a Binder package so that the
full analysis can be run on any machine through internet.


  The Tree of Life is the graphical structure that represents the evolutionary
process from single-cell organisms at the origin of life to the vast
biodiversity we see today. Reconstructing this tree from genomic sequences is
challenging due to the variety of biological forces that shape the signal in
the data, and many of those processes like incomplete lineage sorting and
hybridization can produce confounding information. Here, we present the
mathematical version of the identifiability proofs of phylogenetic networks
under the pseudolikelihood model in SNaQ. We establish that the ability to
detect different hybridization events depends on the number of nodes on the
hybridization blob, with small blobs (corresponding to closely related species)
being the hardest to be detected. Our work focuses on level-1 networks, but
raises attention to the importance of identifiability studies on phylogenetic
inference methods for broader classes of networks.


  After the blockade that many nations suffered to stop the growth of the
incidence curve of COVID-19 during the first half of 2020, they face the
challenge of resuming their social and economic activity. The rapid airborne
transmissibility of SARS-CoV-2, and the absence of a vaccine, calls for active
containment measures to avoid the propagation of transmission chains. The best
strategy up to date, popularly known as Test-Track-Treat (TTT), consist in
testing the population for diagnosis, track the contacts of those infected, and
treat by quarantine all these cases. The dynamical process that better
describes the combined action of the former mechanisms is that of a contagion
process that competes with the spread of the pathogen, cutting off potential
contagion pathways. Here we propose a compartmental model that couples the
dynamics of the infection with the contact tracing and isolation of cases. We
develop an analytical expression for the effective case reproduction number
$R_c(t)$ that reveals the role of contact tracing in the mitigation and
suppression of the epidemics. We show that there is a trade off between the
infection propagation and the isolation of cases. If the isolation is limited
to symptomatic individuals only, the incidence curve can be flattened but not
bended. However, if contact tracing is applied to asymptomatic individuals too,
the strategy can bend the curve and suppress the epidemics. Quantitative
results are dependent on the network topology. We quantify, the most important
indicator of the effectiveness of contact tracing, namely its capacity to
reverse the increasing tendency of the epidemic curve, causing its bending.


  Efficient design of biological sequences will have a great impact across many
industrial and healthcare domains. However, discovering improved sequences
requires solving a difficult optimization problem. Traditionally, this
challenge was approached by biologists through a model-free method known as
"directed evolution", the iterative process of random mutation and selection.
As the ability to build models that capture the sequence-to-function map
improves, such models can be used as oracles to screen sequences before running
experiments. In recent years, interest in better algorithms that effectively
use such oracles to outperform model-free approaches has intensified. These
span from approaches based on Bayesian Optimization, to regularized generative
models and adaptations of reinforcement learning. In this work, we implement an
open-source Fitness Landscape EXploration Sandbox (FLEXS:
github.com/samsinai/FLEXS) environment to test and evaluate these algorithms
based on their optimality, consistency, and robustness. Using FLEXS, we develop
an easy-to-implement, scalable, and robust evolutionary greedy algorithm
(AdaLead). Despite its simplicity, we show that AdaLead is a remarkably strong
benchmark that out-competes more complex state of the art approaches in a
variety of biologically motivated sequence design challenges.


  While often believed to be a passive agent that merely exploits its host's
metabolism, influenza virus has recently been shown to actively move across
glycan-coated surfaces. This form of enzymatically driven surface motility is
currently not well understood and has been loosely linked to burnt-bridge
Brownian ratchet mechanisms. Starting from known properties of influenza's
spike proteins, we develop a physical model that quantitatively describes the
observed motility. It predicts a collectively emerging dynamics of spike
proteins and surface bound ligands that combined with the virus' geometry give
rise to a self-organized rolling propulsion. We show that in contrast to a
Brownian ratchet, the rotary spike drive is not fluctuation driven but operates
optimally as a macroscopic engine in the deterministic regime. The mechanism
also applies to relatives of influenza and to man-made analogues like
DNA-monowheels and should give guidelines for their optimization.


  Being able to store and transmit human genome sequences is an important part
in genomic research and industrial applications. The complete human genome has
3.1 billion base pairs (haploid), and storing the entire genome naively takes
about 3 GB, which is infeasible for large scale usage.
  However, human genomes are highly redundant. Any given individual's genome
would differ from another individual's genome by less than 1%. There are tools
like DNAZip, which express a given genome sequence by only noting down the
differences between the given sequence and a reference genome sequence. This
allows losslessly compressing the given genome to ~ 4 MB in size.
  In this work, we demonstrate additional improvements on top of the DNAZip
library, where we show an additional ~ 11% compression on top of DNAZip's
already impressive results. This would allow further savings in disk space and
network costs for transmitting human genome sequences.


  Dynamical systems running on the top of complex networks has been extensively
investigated for decades. But this topic still remains among the most relevant
issues in complex network theory due to its range of applicability. The contact
process (CP) and the susceptible-infected-susceptible (SIS) model are used
quite often to describe epidemic dynamics. Despite their simplicity, these
models are robust to predict the kernel of real situations. In this work, we
review concisely both processes that are well-known and very applied examples
of models that exhibit absorbing-state phase transitions. In the epidemic
scenario, individuals can be infected or susceptible. A phase transition
between a disease-free (absorbing) state and an active stationary phase (where
a fraction of the population is infected) are separated by an epidemic
threshold. For the SIS model, the central issue is to determine this epidemic
threshold on heterogeneous networks. For the CP model, the main interest is to
relate critical exponents with statistical properties of the network.


  We propose a mathematical model for the transmission dynamics of SARS-CoV-2
in a homogeneously mixing non constant population, and generalize it to a model
where the parameters are given by piecewise constant functions. This allows us
to model the human behavior and the impact of public health policies on the
dynamics of the curve of active infected individuals during a COVID-19 epidemic
outbreak. After proving the existence and global asymptotic stability of the
disease-free and endemic equilibrium points of the model with constant
parameters, we consider a family of Cauchy problems, with piecewise constant
parameters, and prove the existence of pseudo-oscillations between a
neighborhood of the disease-free equilibrium and a neighborhood of the endemic
equilibrium, in a biologically feasible region. In the context of the COVID-19
pandemic, this pseudo-periodic solutions are related to the emergence of
epidemic waves. Then, to capture the impact of mobility in the dynamics of
COVID-19 epidemics, we propose a complex network with six distinct regions
based on COVID-19 real data from Portugal. We perform numerical simulations for
the complex network model, where the objective is to determine a topology that
minimizes the level of active infected individuals and the existence of
topologies that are likely to worsen the level of infection. We claim that this
methodology is a tool with enormous potential in the current pandemic context,
and can be applied in the management of outbreaks (in regional terms) but also
to manage the opening/closing of borders.


  RNA-sequencing (RNA-seq) has become an exemplar technology in modern biology
and clinical applications over the past decade. It has gained immense
popularity in the recent years driven by continuous efforts of the
bioinformatics community to develop accurate and scalable computational tools.
RNA-seq is a method of analyzing the RNA content of a sample using the modern
sequencing platforms. It generates enormous amounts of transcriptomic data in
the form of nucleotide sequences, known as reads. RNA-seq analysis enables the
probing of genes and corresponding transcripts which is essential for answering
important biological questions, such as detecting novel exons, transcripts,
gene expressions, and studying alternative splicing structure. However,
obtaining meaningful biological signals from raw data using computational
methods is challenging due to the limitations of modern sequencing
technologies. The need to leverage these technological challenges have pushed
the rapid development of many novel computational tools which have evolved and
diversified in accordance with technological advancements, leading to the
current myriad population of RNA-seq tools. Our review provides a systemic
overview of RNA-seq technology and 235 available RNA-seq tools across various
domains published from 2008 to 2020, discussing the interdisciplinary nature of
bioinformatics involved in RNA sequencing, analysis, and software development.


  A defining feature of human culture is that knowledge and technology
continually improve over time. Such cumulative cultural evolution (CCE)
probably depends far more heavily on how reliably information is preserved than
on how efficiently it is refined. Therefore, one possible reason that CCE
appears diminished or absent in other species is that it requires accurate but
specialized forms of social learning at which humans are uniquely adept. Here,
we develop a Bayesian model to contrast the evolution of high-fidelity social
learning, which supports CCE, against low-fidelity social learning, which does
not. We find that high-fidelity transmission evolves when (1) social and (2)
individual learning are inexpensive, (3) traits are complex, (4) individual
learning is abundant, (5) adaptive problems are difficult and (6) behaviour is
flexible. Low-fidelity transmission differs in many respects. It not only
evolves when (2) individual learning is costly and (4) infrequent but also
proves more robust when (3) traits are simple and (5) adaptive problems are
easy. If conditions favouring the evolution of high-fidelity transmission are
stricter (3 and 5) or harder to meet (2 and 4), this could explain why social
learning is common, but CCE is rare.


  Recent studies demonstrate that trends in indicators extracted from measured
time series can indicate approaching to an impending transition. Kendall's
{\tau} coefficient is often used to study the trend of statistics related to
the critical slowing down phenomenon and other methods to forecast critical
transitions. Because statistics are estimated from time series, the values of
Kendall's {\tau} are affected by parameters such as window size, sample rate
and length of the time series, resulting in challenges and uncertainties in
interpreting results. In this study, we examine the effects of different
parameters on the distribution of the trend obtained from Kendall's {\tau}, and
provide insights into how to choose these parameters. We also suggest the use
of the non-parametric Mann-Kendall test to evaluate the significance of a
Kendall's {\tau} value. The non-parametric test is computationally much faster
compared to the traditional parametric ARMA test.


  Drawing and visualisation of molecular structures are some of the most common
tasks carried out in structural glycobiology, typically using various software.
In this perspective article, we outline developments in the computational tools
for the sketching, visualisation and modelling of glycans. The article also
provides details on the standard representation of glycans, and
glycoconjugates, which helps the communication of structure details within the
scientific community. We highlight the comparative analysis of the available
tools which could help researchers to perform various tasks related to
structure representation and model building of glycans. These tools can be
useful for glycobiologists or any researcher looking for a ready to use, simple
program for the sketching or building of glycans.


  Chromatin loop extrusion is a popular model for the formation of CTCF loops
and topological domains. Recent HiC data have revealed a strong bias in favour
of a particular arrangement of the CTCF binding motifs that stabilize loops,
and extrusion is the only model to date which can explain this. However, the
model requires a motor to generate the loops, and although cohesin is a strong
candidate for the extruding factor, a suitable motor protein (or a motor
activity in cohesin itself) has yet to be found. Here we explore a new
hypothesis: that there is no motor, and thermal motion within the nucleus
drives extrusion. Using theoretical modelling and computer simulations we ask
whether such diffusive extrusion could feasibly generate loops. Our simulations
uncover an interesting ratchet effect (where an osmotic pressure promotes loop
growth), and suggest, by comparison to recent in vitro and in vivo
measurements, that diffusive extrusion can in principle generate loops of the
size observed in the data.
  Extra View on : C. A. Brackley, J. Johnson, D. Michieletto, A. N. Morozov, M.
Nicodemi, P. R. Cook, and D. Marenduzzo "Non-equilibrium chromosome looping via
molecular slip-links", Physical Review Letters 119, 138101 (2017)


  Delays in protein synthesis cause a confounding effect when constructing Gene
Regulatory Networks (GRNs) from RNA-sequencing time-series data. Accurate GRNs
can be very insightful when modelling development, disease pathways, and drug
side-effects. We present a model which incorporates translation delays by
combining mechanistic equations and Bayesian approaches to fit to experimental
data. This enables greater biological interpretability, and the use of Gaussian
processes enables non-linear expressivity through kernels as well as naturally
accounting for biological variation.


  A simple Monte-Carlo method will be put forward herein, to enable the
extraction of an estimate for the quarantine duration, applicable to visitors
to high-risk regions. Results will be obtained on the basis of an analysis of
the upper tail of the cumulative distribution function of the time span between
the departure of the travellers from the place where the infection occurs and
the time instant when COVID-19 infections may currently be detected. As
expected, the quarantine duration is a decreasing function of the fraction of
the infected travellers, which one is prepared to identify as `acceptable
risk'. The analysis suggests that a maximal $5~\%$ risk (of new infections
originating from subjects who become infective after their quarantine is
lifted) may be associated with a minimal quarantine duration of about eight
days, $1~\%$ with about twelve, and $0.1~\%$ with about sixteen. Unless the
distribution of the duration of short (typically, up to three weeks) travels
departs significantly from the shape assumed in this study, the results of the
present analysis do not provide support for the plans to shorten the quarantine
duration of about ten days to two weeks, which currently applies to travellers
entering most European countries from regions with a high risk of infection.


  The classic SIR model of epidemic dynamics is solved completely by
quadratures, including a time integral transform expanded in a series of
incomplete gamma functions. The model is also generalized to arbitrary
time-dependent infection rates and solved explicitly when the control parameter
depends on the accumulated infections at time $t$. Numerical results are
presented by way of comparison. Autonomous and non-autonomous generalizations
of SIR for interacting regions are also considered, including non-separability
for two or more interacting regions. A reduction of simple SIR models to one
variable leads us to a generalized logistic model, Richards model, which we use
to fit Mexico's COVID-19 data up to day number 134. Forecasting scenarios
resulting from various fittings are discussed. A critique to the applicability
of these models to current pandemic outbreaks in terms of robustness is
provided. Finally, we obtain the bifurcation diagram for a discretized version
of Richards model, displaying period doubling bifurcation to chaos.


  In the idealized Morgan model of crossover, we study the probability
distributions of shared DNA (identical by descent) between individuals having a
wide range of relationships (not just lineal descendants), especially cases for
which previous work produces inaccurate results. Using Monte Carlo simulation,
we show that a particular, complicated functional form with just one continuous
fitted parameter accurately approximates the distributions in all cases tried.
Analysis of that functional form shows that it is close to a normal
distribution, not in shared fraction f, but in the square-root of f. We
describe a multivariate normal model in this variable for use as a practical
framework for several general tasks in forensic genealogy that are currently
done by less-accurate and less well-founded methods.


  The selection pressures that have shaped the evolution of complex traits in
humans remain largely unknown, and in some contexts highly contentious, perhaps
above all where they concern mean trait differences among groups. To date, the
discussion has focused on whether such group differences have any genetic
basis, and if so, whether they are without fitness consequences and arose via
random genetic drift, or whether they were driven by selection for different
trait optima in different environments. Here, we highlight a plausible
alternative, that many complex traits evolve under stabilizing selection in the
face of shifting environmental effects. Under this scenario, there will be
rapid evolution at the loci that contribute to trait variation, even when the
trait optimum remains the same. These considerations underscore the strong
assumptions about environmental effects that are required in ascribing trait
differences among groups to genetic differences.


  A generalisation of the Susceptible-Infectious model is made to include a
time-dependent transmission rate, which leads to a close analytical expression
in terms of a logistic function. The solution can be applied to any continuous
function chosen to describe the evolution of the transmission rate with time.
Taking inspiration from real data of the Covid-19, for the case of cumulative
confirmed positives and deaths, we propose an exponentially decaying
transmission rate with two free parameters, one for its initial amplitude and
another one for its decaying rate. The resultant time-dependent SI model, which
under extra conditions recovers the standard Gompertz functional form, is then
compared with data from selected countries and its parameters fit using
Bayesian inference. We make predictions about the asymptotic number of
confirmed positives and deaths, and discuss the possible evolution of the
disease in each country in terms of our parametrisation of the transmission
rate.


  In order to fulfill cell proliferation and differentiation through cellular
hierarchy, stem cells can undergo either asymmetric or symmetric divisions.
Recent studies pay special attention to the effect of different modes of stem
cell division on the lifetime risk of cancer, and report that symmetric
division is more beneficial to delay the onset of cancer. The fate uncertainty
of symmetric division is considered to be the reason for the cancer-delaying
effect. In this paper we compare asymmetric and symmetric divisions of stem
cells via studying stochastic stem cell models with mutations. Specially, by
using rigorous mathematical analysis we find that both asymmetric and symmetric
models show the same statistical average, but symmetric model shows higher
fluctuation than asymmetric model. We further show that the difference between
the two models would be more remarkable for lower mutation rates. Our work
quantifies the uncertainty of cell division and highlights the significance of
stochasticity for distinguishing between different modes of stem cell division.


  Spatiotemporal modelling of infectious diseases such as COVID-19 involves
using a variety of epidemiological metrics such as regional proportion of cases
or regional positivity rates. Although observing their changes over time is
critical to estimate the regional disease burden, the dynamical properties of
these measures as well as cross-relationships are not systematically explained.
Here we provide a spatiotemporal framework composed of six commonly used and
newly constructed epidemiological metrics and conduct a case study evaluation.
We introduce a refined risk model that is biased neither by the differences in
population sizes nor by the spatial heterogeneity of testing. In particular,
the proposed methodology is useful for the unbiased identification of time
periods with elevated COVID-19 risk, without sensitivity to spatial
heterogeneity of neither population nor testing. Our results also provide
insights regarding regional prioritization of testing and the consequences of
potential synchronization of epidemics between regions.


  Metabolic heterogeneity is widely recognised as the next challenge in our
understanding of non-genetic variation. A growing body of evidence suggests
that metabolic heterogeneity may result from the inherent stochasticity of
intracellular events. However, metabolism has been traditionally viewed as a
purely deterministic process, on the basis that highly abundant metabolites
tend to filter out stochastic phenomena. Here we bridge this gap with a general
method for prediction of metabolite distributions across single cells. By
exploiting the separation of time scales between enzyme expression and enzyme
kinetics, our method produces estimates for metabolite distributions without
the lengthy stochastic simulations that would be typically required for large
metabolic models. The metabolite distributions take the form of Gaussian
mixture models that are directly computable from single-cell expression data
and standard deterministic models for metabolic pathways. The proposed mixture
models provide a systematic method to predict the impact of biochemical
parameters on metabolite distributions. Our method lays the groundwork for
identifying the molecular processes that shape metabolic heterogeneity and its
functional implications in disease.


  Isolation of symptomatic individuals, tracing and testing of their
nonsymptomatic contacts are fundamental strategies for mitigating the current
COVID-19 pandemic. The breaking of contagion chains relies on two complementary
strategies: manual reconstruction of contacts based on interviews and a digital
(app-based) privacy-preserving contact tracing. We compare their effectiveness
using model parameters tailored to describe SARS-CoV-2 diffusion within the
activity-driven model, a general empirically validated framework for network
dynamics. We show that, even for equal probability of tracing a contact, manual
tracing robustly performs better than the digital protocol, also taking into
account the intrinsic delay and limited scalability of the manual procedure.
This result is explained in terms of the stochastic sampling occurring during
the case-by-case manual reconstruction of contacts, contrasted with the
intrinsically prearranged nature of digital tracing, determined by the decision
to adopt the app or not by each individual. The better performance of manual
tracing is enhanced by heterogeneity in agent behavior: superspreaders not
adopting the app are completely invisible to digital contact tracing, while
they can be easily traced manually, due to their multiple contacts. We show
that this intrinsic difference makes the manual procedure dominant in realistic
hybrid protocols.


  By characterising the time evolution of COVID-19 in term of its "velocity"
(log of the new cases per day) and its rate of variation, or "acceleration", we
show that in many countries there has been a deceleration even before lockdowns
were issued. This feature, possibly due to the increase of social awareness,
can be rationalised by a susceptible-hidden-infected-recovered (SHIR) model
introduced by Barnes, in which a hidden (isolated from the virus) compartment
$H$ is gradually populated by susceptible people, thus reducing the
effectiveness of the virus spreading. By introducing a partial hiding
mechanism, for instance due to the impossibility for a fraction of the
population to enter the hidden state, we obtain a model that, although still
sufficiently simple, faithfully reproduces the different deceleration trends
observed in several major countries.


  Genomic data I used in many fields but, it has become known that most of the
platforms used in the sequencing process produce significant errors. This means
that the analysis and inferences generated from these data may have some errors
that need to be corrected. On the two main types of genome errors -
substitution and indels - our work is focused on correcting indels. A deep
learning approach was used to correct the errors in sequencing the chosen
dataset


  Predicting the effect of mutations in proteins is one of the most critical
challenges in protein engineering; by knowing the effect a substitution of one
(or several) residues in the protein's sequence has on its overall properties,
could design a variant with a desirable function. New strategies and
methodologies to create predictive models are continually being developed.
However, those that claim to be general often do not reach adequate
performance, and those that aim to a particular task improve their predictive
performance at the cost of the method's generality. Moreover, these approaches
typically require a particular decision to encode the amino acidic sequence,
without an explicit methodological agreement in such endeavor. To address these
issues, in this work, we applied clustering, embedding, and dimensionality
reduction techniques to the AAIndex database to select meaningful combinations
of physicochemical properties for the encoding stage. We then used the chosen
set of properties to obtain several encodings of the same sequence, to
subsequently apply the Fast Fourier Transform (FFT) on them. We perform an
exploratory stage of Machine-Learning models in the frequency space, using
different algorithms and hyperparameters. Finally, we select the best
performing predictive models in each set of properties and create an assembled
model. We extensively tested the proposed methodology on different datasets and
demonstrated that the generated assembled model achieved notably better
performance metrics than those models based on a single encoding and, in most
cases, better than those previously reported. The proposed method is available
as a Python library for non-commercial use under the GNU General Public License
(GPLv3) license.


  We introduce a fundamentally different paradigm for contact tracing: for each
positive case, do not only ask direct contacts to quarantine; instead, tell
everyone how many relationships away the disease just struck (so, "2" is a
close physical contact of a close physical contact). This new approach, which
has already been deployed in a publicly downloadable app, brings a new tool to
bear on pandemic control, powered by network theory. Like a weather satellite
providing early warning of incoming hurricanes, it empowers individuals to see
transmission approaching from far away, and incites behavior change to directly
avoid exposure. This flipped perspective engages natural self-interested
instincts of self-preservation, reducing reliance on altruism, and the
resulting caution reduces pandemic spread in the social vicinity of each
infection. Consequently, our new system solves the behavior coordination
problem which has hampered many other app-based interventions to date. We also
provide a heuristic mathematical analysis that shows how our system already
achieves critical mass from the user perspective at very low adoption
thresholds (likely below 10% in some common types of communities as indicated
empirically in the first practical deployment); after that point, the design of
our system naturally accelerates further adoption, while also alerting even
non-users of the app. This article seeks to lay the theoretical foundation for
our approach, and to open the area for further research along many dimensions.


  A deterministic mathematical model is formulated and analyzed to study the
transmission dynamics of Nipah virus both qualitatively and numerically.
Existence and stability of equilibria were investigated and the model was
rigorously analyzed. We then incorporated time dependent controls on the model,
using Pontryagin's Maximum Principle to derive necessary conditions for the
optimal control of the disease. We examined various combination strategies so
as to investigate the impact of the controls on the spread of the disease.
Through the numerical results, we found out that the optimal control strategies
help to reduce the burden of the diseases significantly


  Nowadays, coronavirus disease 2019 (COVID-19) poses a great threat to public
health and economy worldwide. Unfortunately, there is yet no effective drug for
this disease. For this, several countries have adopted multiple preventive
interventions to avoid the spread of COVID-19. Here, we propose a delayed
mathematical model to predict the epidemiological trend of COVID-19 in Morocco.
Parameter estimation and sensitivity analysis of the proposed model are
rigorously studied. Moreover, numerical simulations are presented in order to
test the effectiveness of the preventive measures and strategies that were
imposed by the Moroccan authorities and also help policy makers and public
health administration to develop such strategies.


  The SEIR model is a compartmental model used to simulate the dynamics of an
epidemic. In this chapter, we introduce two control functions in the
compartmental SEIR model representing vaccination and plasma transfusion.
Optimal control problems are proposed to study the effects of these two control
measures, on the reduction of infected individuals and increase of recovered
ones, with minimal costs. Up to our knowledge, the plasma transfusion treatment
has never been considered as a control strategy for epidemics mitigation. The
proposed vaccination and treatment strategies may have a real application in
the challenging and hard problem of controlling the COVID-19 pandemic.


  We consider chemical reaction networks modeled by a discrete state and
continuous in time Markov process for the vector copy number of the species and
provide a novel particle filter method for state and parameter estimation based
on exact observation of some of the species in continuous time. The conditional
probability distribution of the unobserved states is shown to satisfy a system
of differential equations with jumps. We provide a method of simulating a
process that is a proxy for the vector copy number of the unobserved species
along with a weight. The resulting weighted Monte Carlo simulation is then used
to compute the conditional probability distribution of the unobserved species.
We also show how our algorithm can be adapted for a Bayesian estimation of
parameters and for the estimation of a past state value based on observations
up to a future time.


  We propose an alternative delayed population growth difference equation model
based on a modification of the Beverton-Holt recurrence, assuming a delay only
in the growth contribution that takes into account that those individuals that
die during the delay, do not contribute to growth. The model introduced differs
from existing delay difference equations in population dynamics, such as the
delayed logistic difference equation, which was formulated as a discretization
of the Hutchinson model. The analysis of our delayed difference equation model
identifies an important critical delay threshold. If the time delay exceeds
this threshold, the model predicts that the population will go extinct for all
non-negative initial conditions and if it is below this threshold, the
population survives and its size converges to a positive globally
asymptotically stable equilibrium that is decreasing in size as the delay
increases. Firstly, we obtain the local stability results by exploiting the
special structure of powers of the Jacobian matrix. Secondly, we show that
local stability implies global stability using two different techniques. For
one set of parameter values, a contraction mapping result is applied, while for
the remaining set of parameter values, we show that the result follows by first
proving that the recurrence structure is eventually monotonic in each of its
arguments.


  Epidemiologists model the dynamics of epidemics in order to propose control
strategies based on pharmaceutical and non-pharmaceutical interventions
(contact limitation, lock down, vaccination, etc). Hand-designing such
strategies is not trivial because of the number of possible interventions and
the difficulty to predict long-term effects. This task can be cast as an
optimization problem where state-of-the-art machine learning algorithms such as
deep reinforcement learning, might bring significant value. However, the
specificity of each domain -- epidemic modelling or solving optimization
problem -- requires strong collaborations between researchers from different
fields of expertise.
  This is why we introduce EpidemiOptim, a Python toolbox that facilitates
collaborations between researchers in epidemiology and optimization.
EpidemiOptim turns epidemiological models and cost functions into optimization
problems via a standard interface commonly used by optimization practitioners
(OpenAI Gym). Reinforcement learning algorithms based on Q-Learning with deep
neural networks (DQN) and evolutionary algorithms (NSGA-II) are already
implemented. We illustrate the use of EpidemiOptim to find optimal policies for
dynamical on-off lock-down control under the optimization of death toll and
economic recess using a Susceptible-Exposed-Infectious-Removed (SEIR) model for
COVID-19. Using EpidemiOptim and its interactive visualization platform in
Jupyter notebooks, epidemiologists, optimization practitioners and others (e.g.
economists) can easily compare epidemiological models, costs functions and
optimization algorithms to address important choices to be made by health
decision-makers.


  Aligning protein interaction networks (PPI) of two or more organisms consists
of finding a mapping of the nodes (proteins) of the networks that captures
important structural and functional associations (similarity). It is a well
studied but difficult problem. It is provably NP-hard in some instances thus
computationally very demanding. The problem comes in several versions: global
versus local alignment; pairwise versus multiple alignment; one-to-one versus
many-to-many alignment. Heuristics to address the various instances of the
problem abound and they achieve some degree of success when their performance
is measured in terms of node and/or edges conservation. However, as the
evolutionary distance between the organisms being considered increases the
results tend to degrade. Moreover, poor performance is achieved when the
considered networks have remarkably different sizes in the number of nodes
and/or edges. Here we address the challenge of analyzing and comparing
different approaches to global network alignment, when a one-to-one mapping is
sought. We consider and propose various measures to evaluate the agreement
between alignments obtained by existing approaches. We show that some such
measures indicate an agreement that is often about the same than what would be
obtained by chance. That tends to occur even when the mappings exhibit a good
performance based on standard measures.


  The COVID-19 pandemic has hit different parts of the world differently: some
regions are still in the rise of the first wave, other regions are now facing a
decline after a first wave, and yet other regions have started to see a second
wave. The current immunity level $\hat i$ in a region is closely related to the
cumulative fraction infected, which primarily depends on two factors: a) the
initial potential for COVID-19 in the region (often quantified by the basic
reproduction number $R_0$), and b) the timing, amount and effectiveness of
preventive measures put in place. By means of a mathematical model including
heterogeneities owing to age, social activity and susceptibility, and allowing
for time-varying preventive measures, the risk for a new epidemic wave and its
doubling time, and how they depend on $R_0$, $\hat i$ and the overall effect of
the current preventive measures, are investigated. Focus lies on quantifying
the minimal overall effect of preventive measures $p_{Min}$ needed to prevent a
future outbreak. The first result shows that the current immunity level $\hat
i$ plays a more influential roll than when immunity is obtained from
vaccination. Secondly, by comparing regions with different $R_0$ and $\hat i$
it is shown that regions with lower $R_0$ and low $\hat i$ may now need higher
preventive measures ($p_{Min}$) compared with other regions having higher $R_0$
but also higher $\hat i$, even when such immunity levels are far from herd
immunity.


  An unsolved fundamental problem in biology and ecology is to predict
observable traits (phenotypes) from a new genetic constitution (genotype) of an
organism under environmental perturbations (e.g., drug treatment). The
emergence of multiple omics data provides new opportunities but imposes great
challenges in the predictive modeling of genotype-phenotype associations.
Firstly, the high-dimensionality of genomics data and the lack of labeled data
often make the existing supervised learning techniques less successful.
Secondly, it is a challenging task to integrate heterogeneous omics data from
different resources. Finally, the information transmission from DNA to
phenotype involves multiple intermediate levels of RNA, protein, metabolite,
etc. The higher-level features (e.g., gene expression) usually have stronger
discriminative power than the lower level features (e.g., somatic mutation). To
address above issues, we proposed a novel Cross-LEvel Information Transmission
network (CLEIT) framework. CLEIT aims to explicitly model the asymmetrical
multi-level organization of the biological system. Inspired by domain
adaptation, CLEIT first learns the latent representation of high-level domain
then uses it as ground-truth embedding to improve the representation learning
of the low-level domain in the form of contrastive loss. In addition, we adopt
a pre-training-fine-tuning approach to leveraging the unlabeled heterogeneous
omics data to improve the generalizability of CLEIT. We demonstrate the
effectiveness and performance boost of CLEIT in predicting anti-cancer drug
sensitivity from somatic mutations via the assistance of gene expressions when
compared with state-of-the-art methods.


  In this work, we developed an efficient approach to compute ensemble averages
in systems with pairwise-additive energetic interactions between the entities.
Methods involving full enumeration of the configuration space result in
exponential complexity. Sampling methods such as Markov Chain Monte Carlo
(MCMC) algorithms have been proposed to tackle the exponential complexity of
these problems; however, in certain scenarios where significant energetic
coupling exists between the entities, the efficiency of the such algorithms can
be diminished. We used a strategy to improve the efficiency of MCMC by taking
advantage of the cluster structure in the interaction energy matrix to bias the
sampling. We pursued two different schemes for the biased MCMC runs and show
that they are valid MCMC schemes. We used both synthesized and real-world
systems to show the improved performance of our biased MCMC methods when
compared to the regular MCMC method. In particular, we applied these algorithms
to the problem of estimating protonation ensemble averages and titration curves
of residues in a protein.


  In recent years, it became clear that super-spreader events play an important
role, particularly in the spread of airborne infections. We investigate a novel
model for super-spreader events, not based on a heterogeneous contact graph but
on a random contact rate: Many individuals become infected synchronously in
single contact events. We use the branching-process approach for contact
tracing to analyze the impact of super-spreader events on the effect of contact
tracing. Here we neglect a tracing delay. Roughly speaking, we find that
contact tracing is more efficient in the presence of super-spreaders if the
fraction of symptomatics is small, the tracing probability is high, or the
latency period is distinctively larger than the incubation period. In other
cases, the effect of contact tracing can be decreased by super-spreaders.
Numerical analysis with parameters suited for SARS-CoV-2 indicates that
super-spreaders do not decrease the effect of contact tracing crucially in case
of that infection.


  The coronavirus disease 2019 (COVID-19) has become a public health emergency
of international concern affecting more than 200 countries and territories
worldwide. As of September 30, 2020, it has caused a pandemic outbreak with
more than 33 million confirmed infections and more than 1 million reported
deaths worldwide. Several statistical, machine learning, and hybrid models have
previously tried to forecast COVID-19 confirmed cases for profoundly affected
countries. Due to extreme uncertainty and nonstationarity in the time series
data, forecasting of COVID-19 confirmed cases has become a very challenging
job. For univariate time series forecasting, there are various statistical and
machine learning models available in the literature. But, epidemic forecasting
has a dubious track record. Its failures became more prominent due to
insufficient data input, flaws in modeling assumptions, high sensitivity of
estimates, lack of incorporation of epidemiological features, inadequate past
evidence on effects of available interventions, lack of transparency, errors,
lack of determinacy, and lack of expertise in crucial disciplines. This chapter
focuses on assessing different short-term forecasting models that can forecast
the daily COVID-19 cases for various countries. In the form of an empirical
study on forecasting accuracy, this chapter provides evidence to show that
there is no universal method available that can accurately forecast pandemic
data. Still, forecasters' predictions are useful for the effective allocation
of healthcare resources and will act as an early-warning system for government
policymakers.


  The concept of the effective infection opportunity population (EIOP) was
incorporated into the SIQR model, and it was assumed that this EIOP would
change with the spread of infection, and this was named as the effective SIQR
model. When calculated with this model, the uninfected population S decreases
with the passage of time. However, when the EIOP N increases because of any
reason, the infection threshold becomes larger than 1. Even after the first
wave seems to have subsided, the infection begins to spread again. Firstly, we
find the curve of EIOP change so that the calculation result by this model
matches the data of the first and second waves. Then, we use this curve to fit
with only the data of the second wave alone, and the third wave is predicted.
In the case of new coronavirus infection, there are various restrictions on
data collection to identify individual coefficients of mathematical models, and
the true value is almost unknown. Therefore, the discussion in this paper is
only about data fitting for predictive calculation. Therefore, the simulation
on the true value is not aimed. However, since the data of infected persons
reflect the true values, the results of data fitting can be used for the
prediction of infected persons, isolated care recipients, inpatients, and
severely ill persons. They are useful for a qualitative understanding of
infection. The idea of EIOP is important in the sense that it connects the
mean-field and the non-mean field, but the existence of data is essential, and
the theory alone cannot simulate the non-mean field. We have developed two
methods for treating the non-mean field cases where we don't have enough data.
We have briefly introduced them.


  The evolution of cooperation has remained an important problem in
evolutionary theory and social sciences. In this regard, a curious question is
why consistent cooperative and defective personalities exist and if they serve
a role in the evolution of cooperation? To shed light on these questions, here,
I consider a population of individuals who possibly play two consecutive rounds
of public goods game, with different enhancement factors. Importantly,
individuals have independent strategies in the two rounds. However, their
strategy in the first round affects the game they play in the second round. I
consider two different scenarios where either only first-round cooperators play
a second public goods game, or both first-round cooperators and first-round
defectors play a second public goods game, but in different groups. The first
scenario can be considered a reward dilemma, and the second can be considered
an assortative public goods game but with independent strategies of the
individuals in the two rounds. Both models show cooperators can survive either
in a fixed point or through different periodic orbits. Interestingly, due to
the emergence of a correlation between the strategies of the individuals in the
two rounds, individuals develop consistent personalities during the evolution.
This, in turn, helps cooperation to flourish. These findings shed new light on
the evolution of cooperation and show how consistent cooperative and defective
personalities can evolve and play a positive role in solving social dilemmas.


  Chase-escape percolation is a variation of the standard epidemic spread
models. In this model, each site can be in one of three states: unoccupied,
occupied by a single prey, or occupied by a single predator. Prey particles
spread to neighboring empty sites at rate $p$, and predator particles spread
only to neighboring sites occupied by prey particles at rate $1$, killing the
prey particle that existed at that site. It was found that the prey can survive
with non-zero probability, if $p>p_c$ with $p_c<1$. Using Monte Carlo
simulations on the square lattice, we estimate the value of $p_c = 0.49451 \pm
0.00001$, and the critical exponents are consistent with the undirected
percolation universality class. We define a discrete-time parallel-update
version of the model, which brings out the relation between chase-escape and
undirected bond percolation. For all $p < p_c$ in $D$-dimensions, the number of
predators in the absorbing configuration has a stretched-exponential
distribution in contrast to the exponential distribution in the standard
percolation theory. We also study the problem starting from the line initial
condition with predator particles on all lattice points of the line $y=0$ and
prey particles on the line $y=1$. In this case, for $p_c<p < 1$, the center of
mass of the fluctuating prey and predator fronts travel at the same speed. This
speed is strictly smaller than the speed of an Eden front with the same value
of $p$, but with no predators. At $p=1$, the fronts undergo a depinning
transition. The fluctuations of the front follow Kardar-Parisi-Zhang scaling
both above and below this depinning transition.


  We develop a model of infection spread that takes into account the existence
of a vulnerable group as well as the variability of the social relations of
individuals. We develop a compartmentalized power-law model, with power-law
connections between the vulnerable and the general population, considering
these connections as well as the connections among the vulnerable as parameters
that we vary in our tests. We use the model to study a number of vaccination
strategies under two hypotheses: first, we assume a limited availability of
vaccine but an infinite vaccination capacity, so that all the available doses
can be administered in a short time (negligible with respect to the evolution
of the epidemic). Then we assume a limited vaccination capacity, so that the
doses are administered in a time non-negligible with respect to the evolution
of the epidemic. We develop optimal strategies for the various social
parameters, where a strategy consists of (1) the fraction of vaccine that is
administered to the vulnerable population and (2) the criterion that is used to
administer it to the general population. In the case of a limited vaccination
capacity, the fraction (1) is a function of time, and we study how to optimize
it to obtain a maximal reduction in the number of victims.


  With the advancements in Artificial intelligence (AI) and the accumulation of
healthrelated big data, it has become increasingly feasible and commonplace to
leverage machine learning technologies to analyze clinical and omics metadata
to assess the possibility of adverse drug reactions or events (ADRs) in the
course of drug discovery. Here, we have described a novel approach that
combined drug-induced gene expression profile from Open TG-GATEs
(Toxicogenomics Project-Genomics Assisted Toxicity Evaluation Systems) and ADR
occurrence information from FAERS (FDA [Food and Drug Administration] Adverse
Events Reporting System) database to predict the likelihood of ADRs. We
generated a total of 14 models using Deep Neural Networks (DNN) to predict
different ADRs; in the validation tests, our models achieved a mean accuracy of
85.71%, indicating that our approach successfully and consistently predicted
ADRs for a wide range of drugs. As an example, we have described the ADR model
in the context of Duodenal ulcer. We believe that our models will help predict
the likelihood of ADRs while testing novel pharmaceutical compounds, and will
be useful for researchers in drug discovery.


  We review the mechanism and consequences of the "bridging-induced
attraction", a generic biophysical principle which underpins some existing
models for chromosome organisation in 3-D. This attraction, which was revealed
in polymer physics-inspired computer simulations, is a generic clustering
tendency arising in multivalent chromatin-binding proteins, and it provides an
explanation for the biogenesis of nuclear bodies and transcription factories
via microphase separation. Including post-translational modification reactions
involving these multivalent proteins can account for the fast dynamics of the
ensuing clusters, as is observed via microscopy and photobleaching experiments.
The clusters found in simulations also give rise to chromatin domains which
conform well with the observation of A/B compartments in HiC experiments.


  Mitochondrial genetic material is widely used for phylogenetic reconstruction
and as a barcode for species identification. Here we study how mito-nuclear
interactions affect the accuracy of species identification by mtDNA, as well as
the speciation process itself. We simulate the evolution of a population of
individuals who carry a recombining nuclear genome and a mitochondrial genome
inherited maternally. We compare a null model fitness landscape that lacks any
mito-nuclear interaction against a scenario in which interactions influence
fitness. Fitness is assigned to individuals according to their mito-nuclear
compatibility, which drives the coevolution of the nuclear and mitochondrial
genomes. When the population breaks into distinct species we analyze the
accuracy of mtDNA barcode for species identification. Remarkably, we find that
species identification by mtDNA is equally accurate in the presence or absence
of mito-nuclear coupling and that the success of the DNA barcode derives mainly
from population geographical isolation during speciation. Nevertheless,
selection imposed by mito-nuclear compatibility influences the diversification
process and leaves signatures in the genetic content and spatial distribution
of the populations, in three ways: phylogenetic trees are more balanced; clades
correlate strongly with the spatial distribution and; there is a substantial
increase in the intraspecies mtDNA similarity. We compare the evolutionary
patterns observed in our model to empirical data from copepods (\textit{T.
californicus}). We find good qualitative agreement in the geographic patterns
and the topology of the phylogenetic tree, provided the model includes
selection based on mito-nuclear interactions. These results highlight the role
of mito-nuclear compatibility in the speciation process and its reconstruction
from genetic data.


  The observation that phenotypic variability is ubiquitous in isogenic
populations has led to a multitude of experimental and theoretical studies
seeking to probe the causes and consequences of this variability. Whether it be
in the context of antibiotic treatments or exponential growth in constant
environments, non-genetic variability has shown to have significant effects on
population dynamics. Here, we review research that elucidates the relationship
between cell-to-cell variability and population dynamics. After summarizing the
relevant experimental observations, we discuss models of bet-hedging and
phenotypic switching. In the context of these models, we discuss how switching
between phenotypes at the single-cell level can help populations survive in
uncertain environments. Next, we review more fine-grained models of phenotypic
variability where the relationship between single-cell growth rates, generation
times and cell sizes is explicitly considered. Variability in these traits can
have significant effects on the population dynamics, even in a constant
environment. We show how these effects can be highly sensitive to the
underlying model assumptions. We close by discussing a number of open
questions, such as how environmental and intrinsic variability interact and
what the role of non-genetic variability in evolutionary dynamics is.


  How does a cell self-organize so that its appendages attain specific lengths
that are convenient for their respective functions? What kind of 'rulers' does
a cell use to measure the length of these appendages? How does a cell transport
structure building materials between the cell body and distal tips of these
appendages so as to regulate their dynamic lengths during various stages of its
lifetime? Some of these questions are addressed here in the context of a
specific cell appendage called flagellum (also called cilium). A "time of
flight" (ToF) mechanism, adapted from the pioneering idea of Galileo, has been
used successfully very recently to explain the length control of flagella by a
biflagellate green algae. Using the same ToF mechanism, here we develop a
stochastic model for the dynamics of flagella in two different types of
monoflagellate unicellular organisms. A unique feature of these monoflagellates
is that these become transiently multi-flagellated during a short span of their
life time. The mean length of the flagella in our model reproduce the trend of
their temporal variation observed in experiments. Moreover, for probing the
intracellular molecular communication between the dynamic flagella of a given
cell, we have computed the correlation in the fluctuations of their lengths
during the multiflagellated stage of the cell cycle by Monte Carlo simulation.


  A second wave of SARS-CoV-2 is unfolding in dozens of countries. However,
this second wave manifests itself strongly in new reported cases, but less in
death counts compared to the first wave. Over the past three months in Germany,
the reported cases increased by a factor five or more, whereas the death counts
hardly grew. This discrepancy fueled speculations that the rise of reported
cases would not reflect a second wave but only wider testing. We find that this
apparent discrepancy can be explained to a large extent by the age structure of
the infected, and predict a pronounced increase of death counts in the near
future, as the spread once again expands into older age groups. To re-establish
control, and to avoid the tipping point when TTI capacity is exceeded, case
numbers have to be lowered. Otherwise the control of the spread and the
protection of vulnerable people will require more restrictive measures latest
when the hospital capacity is reached.


  High-throughput molecular profiling technologies have produced
high-dimensional multi-omics data, enabling systematic understanding of living
systems at the genome scale. Studying molecular interactions across different
data types helps reveal signal transduction mechanisms across different classes
of molecules. In this paper, we develop a novel Bayesian representation
learning method that infers the relational interactions across multi-omics data
types. Our method, Bayesian Relational Learning (BayReL) for multi-omics data
integration, takes advantage of a priori known relationships among the same
class of molecules, modeled as a graph at each corresponding view, to learn
view-specific latent variables as well as a multi-partite graph that encodes
the interactions across views. Our experiments on several real-world datasets
demonstrate enhanced performance of BayReL in inferring meaningful interactions
compared to existing baselines.


  Genes communicate with each other through different regulatory effects, which
lead to the emergence of complex structures in cells, and such structures are
expected to be different for normal and cancerous cells. To study breast cancer
differences, we have investigated the Gene Regulatory Network (GRN) of cells as
inferred from RNA-sequencing data. The GRN is a signed weighted network
corresponding to the inductive or inhibitory interactions. Here we focus on a
particular of motifs in the GRN, the triangles, which are imbalanced if the
number of negative interactions are odd. By studying the stability of
imbalanced triangles in the GRN, we show that the network of cancerous cells
has fewer imbalanced triangles compared to normal. Moreover, in the normal
cells, imbalanced triangles are isolated from the main part of the network,
while such motifs are part of the network's giant component in cancerous cells.
Our result demonstrates that due to genes' collective behavior the complex
structures are different in cancerous cells from those in normal ones.


  Non-transitive dominance and the resulting cyclic loop of three or more
competing species provide a fundamental mechanism to explain biodiversity in
biological and ecological systems. Both Lotka-Volterra and May-Leonard type
model approaches agree that heterogeneity of invasion rates within this loop
does not hazard the coexistence of competing species. While the resulting
abundances of species become heterogeneous, the species who has the smallest
invasion power benefits the most from unequal invasions. Nevertheless, the
effective invasion rate in a predator and prey interaction can also be modified
by breaking the direction of dominance and allowing reversed invasion with a
smaller probability. While this alteration has no particular consequence on the
behavior within the framework of Lotka-Volterra models, the reactions of
May-Leonard systems are highly different. In the latter case, not just the
mentioned "survival of the weakest" effect vanishes, but also the coexistence
of the loop cannot be maintained if the reversed invasion exceeds a threshold
value. Interestingly, the extinction to a uniform state is characterized by a
non-monotonous probability function. While the presence of reversed invasion
does not fully diminish the evolutionary advantage of the original predator
species, but this weakened effective invasion rate helps the related prey
species to collect larger initial area for the final battle between them. The
competition of these processes determines the likelihood in which uniform state
the system terminates.


  To date, the only effective means to respond to the spreading of COVID-19
pandemic are non-pharmaceutical interventions (NPIs), which entail policies to
reduce social activity and mobility restrictions. Quantifying their effect is
difficult, but it is key to reduce their social and economical consequences.
Here, we introduce a meta-population model based on temporal networks,
calibrated on the COVID-19 outbreak data in Italy and apt to evaluate the
outcomes of these two types of NPIs. Our approach combines the advantages of
granular spatial modelling of meta-population models with the ability to
realistically describe social contacts via activity-driven networks. We provide
a valuable framework to assess the viability of different NPIs, varying with
respect to their timing and severity. Results suggest that the effects of
mobility restrictions largely depend on the possibility to implement timely
NPIs in the early phases of the outbreak, whereas activity reduction policies
should be prioritised afterwards.


  We study the multi-strategy stochastic evolutionary game with death-birth
updating in expanding spatial populations of size $N\to \infty$. The model is a
voter model perturbation. For typical populations, we require perturbation
strengths satisfying $1/N\ll w\ll 1$. Under appropriate conditions on the
space, the limiting density processes of strategy are proven to obey the
replicator equation, and the normalized fluctuations converge to a Gaussian
process with the Wright-Fisher covariance function in the limiting densities.
As an application, we resolve in the positive a conjecture from the biological
literature that the expected density processes approximate the replicator
equation on many non-regular graphs.


  Inference of the evolutionary histories of species, commonly represented by a
species tree, is complicated by the divergent evolutionary history of different
parts of the genome. Different loci on the genome can have different histories
from the underlying species tree (and each other) due to processes such as
incomplete lineage sorting (ILS), gene duplication and loss, and horizontal
gene transfer. The multispecies coalescent is a commonly used model for
performing inference on species and gene trees in the presence of ILS. This
paper introduces Lily-T and Lily-Q, two new methods for species tree inference
under the multispecies coalescent. We then compare them to two frequently used
methods, SVDQuartets and ASTRAL, using simulated and empirical data. Both
methods generally showed improvement over SVDQuartets, and Lily-Q was superior
to Lily-T for most simulation settings. The comparison to ASTRAL was more mixed
- Lily-Q tended to be better than ASTRAL when the length of recombination-free
loci was short, when the coalescent population parameter {\theta} was small, or
when the internal branch lengths were longer.


  For prebiotic chemistry to succeed in producing a starting metastable,
autocatalytic and reproducing system subject to evolutionary selection it must
satisfy at least two apparently contradictory requirements: Because such
systems are rare, a search among vast numbers of molecular combinations must
take place naturally, requiring rapid rearrangement and breaking of covalent
bonds. But once a relevant system is found, such rapid disruption and
rearrangement would be very likely to destroy the system before much evolution
could take place. In this paper we explore the possibility, using a model
developed previously, that the search process could occur under different
environmental conditions than the subsequent fixation and growth of a lifelike
chemical system. We use the example of a rapid change in temperature to
illustrate the effect and refer to the rapid change as a `quench' borrowing
terminology from study of the physics and chemistry of glass formation. The
model study shows that interrupting a high temperature nonequilibrium state
with a rapid quench to lower temperatures can substantially increase the
probability of producing a chemical state with lifelike characteristics of
nonequilibrium metastability, internal dynamics and exponential population
growth in time. Previously published data on the length distributions of
proteomes of prokaryotes may be consistent with such an idea and suggest a
prebiotic high temperature `search' phase near the boiling point of water. A
rapid change in pH could have a similar effect. We discuss possible scenarios
on early earth which might have allowed frequent quenches of the sort
considered here to have occurred. The models show a strong dependence of the
effect on the number of chemical monomers available for bond formation.


  The problem of reconstructing a string from its error-prone copies, the trace
reconstruction problem, was introduced by Vladimir Levenshtein two decades ago.
While there has been considerable theoretical work on trace reconstruction,
practical solutions have only recently started to emerge in the context of two
rapidly developing research areas: immunogenomics and DNA data storage. In
immunogenomics, traces correspond to mutated copies of genes, with mutations
generated naturally by the adaptive immune system. In DNA data storage, traces
correspond to noisy copies of DNA molecules that encode digital data, with
errors being artifacts of the data retrieval process. In this paper, we
introduce several new trace generation models and open questions relevant to
trace reconstruction for immunogenomics and DNA data storage, survey
theoretical results on trace reconstruction, and highlight their connections to
computational biology. Throughout, we discuss the applicability and
shortcomings of known solutions and suggest future research directions.


  Antibody therapeutics and vaccines are among our last resort to end the
raging COVID-19 pandemic. They, however, are prone to over 5,000 mutations on
the spike (S) protein uncovered by a Mutation Tracker based on over 200,000
genome isolates. It is imperative to understand how mutations would impact
vaccines and antibodies in the development. In this work, we study the
mechanism, frequency, and ratio of mutations on the S protein. Additionally, we
use 56 antibody structures and analyze their 2D and 3D characteristics.
Moreover, we predict the mutation-induced binding free energy (BFE) changes for
the complexes of S protein and antibodies or ACE2. By integrating genetics,
biophysics, deep learning, and algebraic topology, we reveal that most of 462
mutations on the receptor-binding domain (RBD) will weaken the binding of S
protein and antibodies and disrupt the efficacy and reliability of antibody
therapies and vaccines. A list of 31 vaccine escape mutants is identified,
while many other disruptive mutations are detailed as well. We also unveil that
about 65\% existing RBD mutations, including those variants recently found in
the United Kingdom (UK) and South Africa, are binding-strengthen mutations,
resulting in more infectious COVID-19 variants. We discover the disparity
between the extreme values of RBD mutation-induced BFE strengthening and
weakening of the bindings with antibodies and ACE2, suggesting that SARS-CoV-2
is at an advanced stage of evolution for human infection, while the human
immune system is able to produce optimized antibodies. This discovery implies
the vulnerability of current vaccines and antibody drugs to new mutations. Our
predictions were validated by comparison with more than 1,400 deep mutations on
the S protein RBD. Our results show the urgent need to develop new
mutation-resistant vaccines and antibodies and to prepare for seasonal
vaccinations.


  Binary trees are fundamental objects in models of evolutionary biology and
population genetics. Here, we discuss some of their combinatorial and
structural properties as they depend on the tree class considered. Furthermore,
the process by which trees are generated determines the probability
distribution in tree space. Yule trees, for instance, are generated by a pure
birth process. When considered as unordered, they have neither a closed-form
enumeration nor a simple probability distribution. But their ordered siblings
have both. They present the object of choice when studying tree structure in
the framework of evolving genealogies.


  Experimental approaches have been applied to address questions in
understanding three-dimensional chromatin organisation and function. As
datasets increase in size and complexity, it becomes a challenge to reach a
mechanistic interpretation of experimental results. Polymer simulations and
mechanistic modelling have been applied to explain experimental observations,
and the links to different aspects of genome function. Here, we provide a guide
for biologists, explaining different simulation approaches and the contexts in
which they have been used.


  Network alignment aims to uncover topologically similar regions in the
protein-protein interaction (PPI) networks of two or more species under the
assumption that topologically similar regions tend to perform similar
functions. Although there exist a plethora of both network alignment algorithms
and measures of topological similarity, currently no gold standard exists for
evaluating how well either is able to uncover functionally similar regions.
Here we propose a formal, mathematically and statistically rigorous method for
evaluating the statistical significance of shared GO terms in a global, 1-to-1
alignment between two PPI networks. We use combinatorics to precisely count the
number of possible network alignments in which $k$ proteins share a particular
GO term. When divided by the number of all possible network alignments, this
provides an explicit, exact $p$-value for a network alignment with respect to a
particular GO term. Just as with BLAST's p-values and bit-scores, this method
is designed not to guide the formation of any particular alignment, but instead
to provide an after-the-fact evaluation of a fixed, given alignment.


  Models in systems biology are mathematical descriptions of biological
processes that are used to answer questions and gain a better understanding of
biological phenomena. Dynamic models represent the network through rates of the
production and consumption for the individual species. The ordinary
differential equations that describe rates of the reactions in the model
include a set of parameters. The parameters are important quantities to
understand and analyze biological systems. Moreover, the perturbation of the
kinetic parameters are correlated with upregulation of the system by
cell-intrinsic and cell-extrinsic factors, including mutations and the
environment changes. Here, we aim at using well-established models of
biological pathways to identify parameter values and point their potential
perturbation/deviation. We present our population-based optimization framework
that is able to identify kinetic parameters in the dynamic model based on only
input and output data (i.e., timecourses of selected metabolites). Our approach
can deal with the identification of the non-measurable parameters as well as
with discovering deviation of the parameters. We present our proposed
optimization framework on the example of the well-studied glycolytic pathway in
Saccharomyces cerevisiae.


  One of the main problems in controlling COVID-19 epidemic spread is the delay
in confirming cases. Having information on changes in the epidemic evolution or
outbreaks rise before lab-confirmation is crucial in decision making for Public
Health policies. We present an algorithm to estimate on-stream the number of
COVID-19 cases using the data from telephone calls to a COVID-line. By modeling
the calls as background (proportional to population) plus signal (proportional
to infected), we fit the calls in Province of Buenos Aires (Argentina) with
coefficient of determination $R^2 > 0.85$. This result allows us to estimate
the number of cases given the number of calls from a specific district, days
before the lab results are available. We validate the algorithm with real data.
We show how to use the algorithm to track on-stream the epidemic, and present
the Early Outbreak Alarm to detect outbreaks in advance to lab results. One key
point in the developed algorithm is a detailed track of the uncertainties in
the estimations, since the alarm uses the significance of the observables as a
main indicator to detect an anomaly. We present the details of the explicit
example in Villa Azul (Quilmes) where this tool resulted crucial to control an
outbreak on time. The presented tools have been designed in urgency with the
available data at the time of the development, and therefore have their
limitations which we describe and discuss. We consider possible improvements on
the tools, many of which are currently under development.


  With the recent advances in machine learning for quantum chemistry, it is now
possible to predict the chemical properties of compounds and to generate novel
molecules. Existing generative models mostly use a string- or graph-based
representation, but the precise three-dimensional coordinates of the atoms are
usually not encoded. First attempts in this direction have been proposed, where
autoregressive or GAN-based models generate atom coordinates. Those either lack
a latent space in the autoregressive setting, such that a smooth exploration of
the compound space is not possible, or cannot generalize to varying chemical
compositions. We propose a new approach to efficiently generate molecular
structures that are not restricted to a fixed size or composition. Our model is
based on the variational autoencoder which learns a translation-, rotation-,
and permutation-invariant low-dimensional representation of molecules. Our
experiments yield a mean reconstruction error below 0.05 Angstrom,
outperforming the current state-of-the-art methods by a factor of four, and
which is even lower than the spatial quantization error of most chemical
descriptors. The compositional and structural validity of newly generated
molecules has been confirmed by quantum chemical methods in a set of
experiments.


  Organismal phenotypes emerge from a complex set of genotypic interactions.
While technological advances in sequencing provide a quantitative description
of an organism's genotype, characterization of an organism's physical phenotype
lags far behind. Here, we relate genotype to the complex and multi-dimensional
phenotype of an anatomical structure using the Drosophila wing as a model
system. We develop a mathematical approach that enables a robust description of
biologically salient phenotypic variation. Analysing natural phenotypic
variation, and variation generated by weak perturbations in genetic and
environmental conditions during development, we observe a highly constrained
set of wing phenotypes. In a striking ex-ample of dimensionality reduction, the
nature of varieties produced by the Drosophila developmental pro-gram is
constrained to a single integrated mode of variation in the wing. Our strategy
demonstrates the emergent simplicity manifest in the genotype-to-phenotype map
in the Drosophila wing and may represent a general approach for interrogating a
variety of genotype-phenotype relationships.


  We study the spread of Covid-19 across neighbourhoods of cities in the
developing world and find that small numbers of neighbourhoods account for a
majority of cases (k-index~0.7). We also find that the countrywide distribution
of cases across states/provinces in these nations also displays similar
inequality, indicating self-similarity across scales. Neighbourhoods with slums
are found to contain the highest density of cases across all cities under
consideration, revealing that slums constitute the most at-risk urban locations
in this epidemic. We present a stochastic network model to study the spread of
a respiratory epidemic through physically proximate and accidental daily human
contacts in a city, and simulate outcomes for a city with two kinds of
neighbourhoods - slum and non-slum. The model reproduces observed empirical
outcomes for a broad set of parameter values - reflecting the potential
validity of these findings for epidemic spread in general, especially across
cities of the developing world. We also find that distribution of cases becomes
less unequal as the epidemic runs its course, and that both peak and cumulative
caseloads are worse for slum neighbourhoods than non-slums at the end of an
epidemic. Large slums in the developing world therefore contain the most
vulnerable populations in an outbreak, and the continuing growth of
metropolises in Asia and Africa presents significant challenges for future
respiratory outbreaks from perspectives of public health and socioeconomic
equity.


  The dramatic outbreak of the coronavirus disease 2019 (COVID-19) pandemics
and its ongoing progression boosted the scientific community's interest in
epidemic modeling and forecasting. The SIR (Susceptible-Infected-Removed) model
is a simple mathematical model of epidemic outbreaks, yet for decades it evaded
the efforts of the community to derive an explicit solution. The present work
demonstrates that this is a non-trivial task. Notably, it is proven that the
explicit solution of the model requires the introduction of a new
transcendental special function, related to the Wright's Omega function. The
present manuscript reports new analytical results and numerical routines
suitable for parametric estimation of the SIR model. The manuscript introduces
iterative algorithms approximating the incidence variable, which allows for
estimation of the model parameters from the numbers of observed cases. The
numerical approach is exemplified with data from the European Centre for
Disease Prevention and Control (ECDC) for several European countries in the
period Jan 2020 -- Jun 2020.


  A method of Transitional Automorphic Mapping of the Genome on Itself (TAMGI)
is aimed at combining detection and reconstruction of correlational and
quasi-periodic motifs in the viral genomic RNA/DNA sequences. The motifs
reconstructed by TAMGI are robust with respect to indels and point mutations
and can be tried as putative therapeutic targets. We developed and tested the
relevant theory and statistical criteria for TAMGI applications. The
applications of TAMGI are illustrated by the study of motifs in the genomes of
the severe acute respiratory syndrome coronaviruses SARS-CoV and SARS-CoV-2
(the latter coronavirus SARS-CoV-2 being responsible for the COVID-19 pandemic)
packaged within filament-like helical capsid. Such ribonucleocapsid is
transported into spherical membrane envelope with incorporated spike
glycoproteins. Two other examples concern the genomes of viruses with
icosahedral capsids, satellite tobacco mosaic virus (STMV) and bacteriophage
PHIX174. A part of the quasi-periodic motifs in these viral genomes was evolved
due to weakly specific cooperative interaction between genomic ssRNA/ssDNA and
nucleocapsid proteins. The symmetry of the capsids leads to the natural
selection of specific quasi-periodic motifs in the related genomic sequences.
Generally, TAMGI provides a convenient tool for the study of numerous molecular
mechanisms with participation of both quasi-periodic motifs and complete
repeats, the genome organization, contextual analysis of cis/trans regulatory
elements, data mining, and correlations in the genomic sequences.


  We discuss a polymer model for the 3D organization of human chromosomes. A
chromosome is represented by a string of beads, with each bead being "colored"
according to 1D bioinformatic data (e.g., chromatin state, histone
modification, GC content). Individual spheres (representing bi- and
multi-valent transcription factors) can bind reversibly and selectively to
beads with the appropriate color. During molecular dynamics simulations, the
factors bind, and the string spontaneously folds into loops, rosettes, and
topologically-associating domains (TADs). This organization occurs in the
absence of any specified interactions between distant DNA segments, or between
transcription factors. A comparison with Hi-C data shows that simulations
predict the location of most boundaries between TADs correctly. The model is
"fitting-free" in the sense that it does not use Hi-C data as an input;
consequently, one of its strengths is that it can -- in principle -- be used to
predict the 3D organization of any region of interest, or whole chromosome, in
a given organism, or cell line, in the absence of existing Hi-C data. We
discuss how this simple model might be refined to include more transcription
factors and binding sites, and to correctly predict contacts between convergent
CTCF binding sites.


  The Covid-19 pandemic is ongoing worldwide, and the damage it has caused is
unprecedented. For prevention, South Korea has adopted a local quarantine
strategy rather than a global lockdown. This approach not only minimizes
economic damage, but it also efficiently prevents the spread of the disease. In
this work, the spread of COVID-19 under local quarantine measures is modeled
using the Susceptible-Exposed-Infected-Recovered model on complex networks. In
this network approach, the links connected to isolated people are disconnected
and then reinstated when they are released. This link dynamics leads to
time-dependent reproduction number. Numerical simulations are performed on
networks with reaction rates estimated from empirical data. The temporal
pattern of the cumulative number of confirmed cases is then reproduced. The
results show that a large number of asymptomatic infected patients are detected
as they are quarantined together with infected patients. Additionally, possible
consequences of the breakdowns of local quarantine measures and social
distancing are considered.


  Deficiency zero is an important network structure and has been the focus of
many celebrated results within reaction network theory. In our previous paper
\textit{Prevalence of deficiency zero reaction networks in an Erd\H os-R\'enyi
framework}, we provided a framework to quantify the prevalence of deficiency
zero among randomly generated reaction networks. Specifically, given a randomly
generated binary reaction network with $n$ species, with an edge between two
arbitrary vertices occurring independently with probability $p_n$, we
established the threshold function $r(n)=\frac{1}{n^3}$ such that the
probability of the random network being deficiency zero converges to 1 if
$\frac{p_n}{r(n)}\to 0$ and converges to 0 if $\frac{p_n}{r(n)}\to\infty$, as
$n \to \infty$.
  With the base Erd\H os-R\'enyi framework as a starting point, the current
paper provides a significantly more flexible framework by weighting the edge
probabilities via control parameters $\alpha_{i,j}$, with $i,j\in \{0,1,2\}$
enumerating the types of possible vertices (zeroth, first, or second order).
The control parameters can be chosen to generate random reaction networks with
a specific underlying structure, such as "closed" networks with very few inflow
and outflow reactions, or "open" networks with abundant inflow and outflow.
Under this new framework, for each choice of control parameters
$\{\alpha_{i,j}\}$, we establish a threshold function $r(n,\{\alpha_{i,j}\})$
such that the probability of the random network being deficiency zero converges
to 1 if $\frac{p_n}{r(n,\{\alpha_{i,j}\})}\to 0$ and converges to 0 if
$\frac{p_n}{r(n,\{\alpha_{i,j}\})}\to \infty$.


  Social insect colonies routinely face large vertebrate predators, against
which they need to mount a collective defense. To do so, honeybees use an alarm
pheromone that recruits nearby bees into mass stinging of the perceived threat.
This alarm pheromone is carried directly on the stinger, hence its
concentration builds up during the course of the attack. Here, we investigate
how individual bees react to different alarm pheromone concentrations, and how
this evolved response-pattern leads to better coordination at the group level.
We first present an individual dose-response curve to the alarm pheromone,
obtained experimentally. Second, we apply Projective Simulation to model each
bee as an artificial learning agent that relies on the pheromone concentration
to decide whether to sting or not. If the emergent collective performance
benefits the colony, the individual reactions that led to it are enhanced via
reinforcement learning, thus emulating natural selection. Predators are modeled
in a realistic way so that the effect of factors such as their resistance,
their killing rate or their frequency of attacks can be studied. We are able to
reproduce the experimentally measured response-pattern of real bees, and to
identify the main selection pressures that shaped it. Finally, we apply the
model to a case study: by tuning the parameters to represent the environmental
conditions of European or African bees, we can predict the difference in
aggressiveness observed between these two subspecies.


  Causal gene networks model the flow of information within a cell, but
reconstructing them from omics data is challenging because correlation does not
imply causation. Combining genomics and transcriptomics data from a segregating
population allows to orient the direction of causality between gene expression
traits using genomic variants. Instrumental-variable methods (IV) use a local
expression quantitative trait locus (eQTL) as a randomized instrument for a
gene's expression level, and assign target genes based on distal eQTL
associations. Mediation-based methods (ME) additionally require that distal
eQTL associations are mediated by the source gene. Here we used Findr, a
software providing uniform implementations of IV, ME, and coexpression-based
methods, a recent dataset of 1,012 segregants from a cross between two budding
yeast strains, and the YEASTRACT database of known transcriptional interactions
to compare causal gene network inference methods. We found that causal
inference methods result in a significant overlap with the ground-truth,
whereas coexpression did not perform better than random. A subsampling analysis
revealed that the performance of ME decreases at large sample sizes, due to a
loss of sensitivity when residual correlations become significant. IV methods
contain false positive predictions, due to genomic linkage between eQTL
instruments. IV and ME methods also have complementary roles for identifying
causal genes underlying transcriptional hotspots. IV methods correctly
predicted STB5 targets for a hotspot centred on the transcription factor STB5,
whereas ME failed due to Stb5p auto-regulating its own expression. ME suggests
a new candidate gene, DNM1, for a hotspot on Chr XII, where IV methods could
not distinguish between multiple genes located within the hotspot.


  The prediction of structure dependent molecular properties, such as collision
cross sections as measured using ion mobility spectrometry, are crucially
dependent on the selection of the correct population of molecular conformers.
Here, we report an in-depth evaluation of multiple conformation selection
techniques, including simple averaging, Boltzmann weighting, lowest energy
selection, low energy threshold reductions, and similarity reduction.
Generating 50,000 conformers each for 18 molecules, we used the In Silico
Chemical Library Engine (ISiCLE) to calculate the collision cross sections for
the entire dataset. First, we employed Monte Carlo simulations to understand
the variability between conformer structures as generated using simulated
annealing. Then we employed Monte Carlo simulations to the aforementioned
conformer selection techniques applied on the simulated molecular property -
the ion mobility collision cross section. Based on our analyses, we found
Boltzmann weighting to be a good tradeoff between precision and theoretical
accuracy. Combining multiple techniques revealed that energy thresholds and
root-mean-squared deviation-based similarity reductions can save considerable
computational expense while maintaining property prediction accuracy. Molecular
dynamic conformer generation tools like AMBER can continue to generate new
lowest energy conformers even after tens of thousands of generations,
decreasing precision between runs. This reduced precision can be ameliorated
and theoretical accuracy increased by running density functional theory
geometry optimization on carefully selected conformers.


  Objective: To introduce a method for tracking results and utilization of
Artificial Intelligence (tru-AI) in radiology. By tracking both large-scale
utilization and AI results data, the tru-AI approach is designed to calculate
surrogates for measuring important disease-related observational quantities
over time, such as the prevalence of intracranial hemorrhage during the
COVID-19 pandemic outbreak. Methods: To quantitatively investigate the clinical
applicability of the tru-AI approach, we analyzed service requests for
automatically identifying intracranial hemorrhage (ICH) on head CT using a
commercial AI solution. This software is typically used for AI-based
prioritization of radiologists' reading lists for reducing turnaround times in
patients with emergent clinical findings, such as ICH or pulmonary embolism.We
analyzed data of N=9,421 emergency-setting non-contrast head CT studies at a
major US healthcare system acquired from November 1, 2019 through June 2, 2020,
and compared two observation periods, namely (i) a pre-pandemic epoch from
November 1, 2019 through February 29, 2020, and (ii) a period during the
COVID-19 pandemic outbreak, April 1-30, 2020. Results: Although daily CT scan
counts were significantly lower during (40.1 +/- 7.9) than before (44.4 +/-
7.6) the COVID-19 outbreak, we found that ICH was more likely to be observed by
AI during than before the COVID-19 outbreak (p<0.05), with approximately one
daily ICH+ case more than statistically expected. Conclusion: Our results
suggest that, by tracking both large-scale utilization and AI results data in
radiology, the tru-AI approach can contribute clinical value as a versatile
exploratory tool, aiming at a better understanding of pandemic-related effects
on healthcare.


  Computational docking methods can provide structural models of
protein-protein complexes, but protein backbone flexibility upon association
often thwarts accurate predictions. In recent blind challenges, medium or high
accuracy models were submitted in less than 20% of the "difficult" targets
(with significant backbone change or uncertainty). Here, we describe recent
developments in protein-protein docking and highlight advances that tackle
backbone flexibility. In molecular dynamics and Monte Carlo approaches,
enhanced sampling techniques have reduced time-scale limitations. Internal
coordinate formulations can now capture realistic motions of monomers and
complexes using harmonic dynamics. And machine learning approaches adaptively
guide docking trajectories or generate novel binding site predictions from deep
neural networks trained on protein interfaces. These tools poise the field to
break through the longstanding challenge of correctly predicting complex
structures with significant conformational change.


  We analyze the dynamics of a population of independent random walkers on a
graph and develop a simple model of epidemic spreading. We assume that each
walker visits independently the nodes of a finite ergodic graph in a
discrete-time markovian walk governed by his specific transition matrix. With
this assumption, we first derive an upper bound for the reproduction numbers.
Then we assume that a walker is in one of the states: susceptible, infectious,
or recovered. An infectious walker remains infectious during a certain
characteristic time. If an infectious walker meets a susceptible one on the
same node there is a certain probability for the susceptible walker to get
infected. By implementing this hypothesis in computer simulations we study the
space-time evolution of the emerging infection patterns. Generally, random walk
approaches seem to have a large potential to study epidemic spreading and to
identify the pertinent parameters in epidemic dynamics.


  The COVID-19 pandemic has led to unprecedented efforts to identify drugs that
can reduce its associated morbidity/mortality rate. Computational chemistry
approaches hold the potential for triaging potential candidates far more
quickly than their experimental counterparts. These methods have been widely
used to search for small molecules that can inhibit critical proteins involved
in the SARS-CoV-2 replication cycle. An important target is the SARS-CoV-2 main
protease Mpro, an enzyme that cleaves the viral polyproteins into individual
proteins required for viral replication and transcription. Unfortunately,
standard computational screening methods face difficulties in ranking diverse
ligands to a receptor due to disparate ligand scaffolds and varying charge
states. Here, we describe full density functional quantum mechanical (DFT/QM)
simulations of Mpro in complex with various ligands to obtain absolute ligand
binding energies. Our calculations are enabled by a new cloud-native parallel
DFT/QM implementation running on computational resources from Amazon Web
Services (AWS). The results we obtain are promising: the approach is quite
capable of scoring a very diverse set of existing drug compounds for their
affinities to Mpro and suggest the DFT/QM approach is potentially more broadly
applicable to repurpose screening against this target. In addition, each DFT/QM
simulation required only ~1 hour (wall clock time) per ligand. The fast
turnaround time raises the practical possibility of a broad application of
large-scale quantum mechanics in the drug discovery pipeline at stages where
ligand diversity is essential.


  Measures of phylogenetic balance, such as the Colless and Sackin indices,
play an important role in phylogenetics. Unfortunately, these indices are
specifically designed for phylogenetic trees, and do not extend naturally to
phylogenetic networks (which are increasingly used to describe reticulate
evolution). This led us to consider a lesser-known balance index, whose
definition is based on a probabilistic interpretation that is equally
applicable to trees and to networks. This index, known as the $B_2$ index, was
first proposed by Shao and Sokal in 1990. Surprisingly, it does not seem to
have been studied mathematically since. Likewise, it is used only sporadically
in the biological literature, where it tends to be viewed as arcane. In this
paper, we study mathematical properties of $B_2$ such as its expectation and
variance under the most common models of random trees and its extremal values
over various classes of phylogenetic networks. We also assess its relevance in
biological applications, and find it to be comparable to that of the Colless
and Sackin indices. Altogether, our results call for a reevaluation of the
status of this somewhat forgotten measure of phylogenetic balance.


  Exploiting the SIQR model for COVID-19, I show that the wavy infection curve
in Japan is the result of fluctuation of policy on isolation measure imposed by
the government and obeyed by citizens. Assuming the infection coefficient be a
two-valued function of the number of daily confirmed new cases, I show that
when the removal rate of infected individuals is between these two values, the
wavy infection curve is self-organized. On the basis of the infection curve, I
classify the outbreak of COVID-19 into five types and show that these
differences can be related to the relative magnitude of the transmission
coefficient and the quarantine rate of infected individuals.


  We introduced a more general predator-prey model to analyze the paradox of
enrichment. We hope the results obtained for the model can guide us on
identifying real field paradox of enrichment.


  Despite recent advancements in deep learning methods for protein structure
prediction and representation, little focus has been directed at the
simultaneous inclusion and prediction of protein backbone and sidechain
structure information. We present SidechainNet, a new dataset that directly
extends the ProteinNet dataset. SidechainNet includes angle and atomic
coordinate information capable of describing all heavy atoms of each protein
structure. In this paper, we provide background information on the availability
of protein structure data and the significance of ProteinNet. Thereafter, we
argue for the potentially beneficial inclusion of sidechain information through
SidechainNet, describe the process by which we organize SidechainNet, and
provide a software package (https://github.com/jonathanking/sidechainnet) for
data manipulation and training with machine learning models.