  Haplotype-resolved de novo assembly is the ultimate solution to the study of
sequence variations in a genome. However, existing algorithms either collapse
heterozygous alleles into one consensus copy or fail to cleanly separate the
haplotypes to produce high-quality phased assemblies. Here we describe hifiasm,
a new de novo assembler that takes advantage of long high-fidelity sequence
reads to faithfully represent the haplotype information in a phased assembly
graph. Unlike other graph-based assemblers that only aim to maintain the
contiguity of one haplotype, hifiasm strives to preserve the contiguity of all
haplotypes. This feature enables the development of a graph trio binning
algorithm that greatly advances over standard trio binning. On three human and
five non-human datasets, including California redwood with a $\sim$30-gigabase
hexaploid genome, we show that hifiasm frequently delivers better assemblies
than existing tools and consistently outperforms others on haplotype-resolved
assembly.


  Ordinary differential equations (ODE) models have a wide variety of
applications in the fields of mathematics, statistics, and the sciences. Though
they are widely used, these models are sometimes viewed as inflexible with
respect to the incorporation of time delays. The Generalized Linear Chain Trick
(GLCT) serves as a way for modelers to incorporate much more flexible delay or
dwell time distribution assumptions than the usual exponential and Erlang
distributions. In this paper we demonstrate how the GLCT can be used to
generate new ODE models by generalizing or approximating existing models to
yield much more general ODEs with phase-type distributed delays or dwell times.


  We provide an overview of the methods that can be used for prediction under
uncertainty and data fitting of dynamical systems, and of the fundamental
challenges that arise in this context. The focus is on SIR-like models, that
are being commonly used when attempting to predict the trend of the COVID-19
pandemic. In particular, we raise a warning flag about identifiability of the
parameters of SIR-like models; often, it might be hard to infer the correct
values of the parameters from data, even for very simple models, making it
non-trivial to use these models for meaningful predictions. Most of the points
that we touch upon are actually generally valid for inverse problems in more
general setups.


  A central question of evolutionary dynamics on graphs is whether or not a
mutation introduced in a population of residents survives and eventually even
spreads to the whole population, or gets extinct. The outcome naturally depends
on the fitness of the mutant and the rules by which mutants and residents may
propagate on the network, but arguably the most determining factor is the
network structure. Some structured networks are transient amplifiers. They
increase for a certain fitness range the fixation probability of beneficial
mutations as compared to a well-mixed population. We study a perturbation
methods for identifying transient amplifiers for death-Birth updating. The
method includes calculating the coalescence times of random walks on graphs and
finding the vertex with the largest remeeting time. If the graph is perturbed
by removing an edge from this vertex, there is a certain likelihood that the
resulting perturbed graph is a transient amplifier. We test all pairwise
nonisomorphic cubic and quartic regular graphs up to a certain size and thus
cover the whole structural range expressible by these graphs. We carry out a
spectral analysis and show that the graphs from which the transient amplifiers
can be constructed share certain structural properties. The graphs are
path-like, have low conductance and are rather easy to divide into subgraphs by
removing edges and/or vertices. This is connected with the subgraphs being
identical (or almost identical) building blocks and the frequent occurrence of
cut and/or hinge vertices. Identifying spectral and structural properties may
promote finding and designing such networks.


  Carbohydrates such as the trisaccharide motif LeX are key constituents of
cell surfaces. Despite intense research, the interactions between carbohydrates
of apposing cells or membranes are not well understood. In this article, we
investigate carbohydrate-carbohydrate interactions in membrane adhesion as well
as in solution with extensive atomistic molecular dynamics simulations that
exceed the simulation times of previous studies by orders of magnitude. For
LeX, we obtain association constants of soluble carbohydrates, adhesion
energies of lipid-anchored carbohydrates, and maximally sustained forces of
carbohydrate complexes in membrane adhesion that are in good agreement with
experimental results in the literature. Our simulations thus appear to provide
a realistic, detailed picture of LeX-LeX interactions in solution and during
membrane adhesion. In this picture, the LeX-LeX interactions are fuzzy, i.e.
LeX pairs interact in a large variety of short-lived, bound conformations. For
the synthetic tetrasaccharide Lac 2, which is composed of two lactose units, we
observe similarly fuzzy interactions and obtain association constants of both
soluble and lipid-anchored variants that are comparable to the corresponding
association constants of LeX. The fuzzy, weak carbohydrate-carbohydrate
interactions quantified in our simulations thus appear to be a generic feature
of small, neutral carbohydrates such as LeX and Lac 2.


  While allostery is of paramount importance for protein regulation, the
underlying dynamical process of ligand (un)binding at one site, resulting time
evolution of the protein structure, and change of the binding affinity at a
remote site is not well understood. Here the ligand-induced conformational
transition in a widely studied model system of allostery, the PDZ2 domain, is
investigated by transient infrared spectroscopy accompanied by molecular
dynamics simulations. To this end, an azobenzene derived photoswitch is linked
to a peptide ligand in a way that its binding affinity to the PDZ2 domain
changes upon switching, thus initiating an allosteric transition in the PDZ2
domain protein. The subsequent response of the protein, covering four decades
of time ranging from $\sim$1~ns to $\sim$10~$\mu$s, can be rationalize by a
remodelling of its rugged free energy landscape, with ver subtle shifts in the
populations of a small number of structurally well defined states. It is
proposed that structurally and dynamically driven allostery, often discussed as
limiting scenarios of allosteric communication, actually go hand-in-hand,
allowing the protein to adapt its free energy landscape to incoming signals.


  COVID-19 epidemics in Japan and Tokyo were analyzed by a fundamental equation
of the dynamic phase transition. As a result, the epidemic was found to be in
good agreement with the random nucleation and linear growth model suggesting
that the epidemic between March 13, 2020 and May 22, 2020 was simply
rate-limited by the three constant-parameters: the initial susceptible, domain
growth rate, and nucleation decay constant. This model provides a good
predictor of the epidemic because it consists of one equation and the initial
specific plot is linear.


  This work is devoted to the study of a stochastic logistic growth model with
and without the Allee effect. Such a model describes the evolution of a
population under environmental stochastic fluctuations and is in the form of a
stochastic differential equation driven by multiplicative Gaussian noise. With
the help of the associated Fokker-Planck equation, we analyze the population
extinction probability and the probability of reaching a large population size
before reaching a small one. We further study the impact of the harvest rate,
noise intensity, and the Allee effect on population evolution. The analysis and
numerical experiments show that if the noise intensity and harvest rate are
small, the population grows exponentially, and upon reaching the carrying
capacity, the population size fluctuates around it. In the stochastic
logistic-harvest model without the Allee effect, when noise intensity becomes
small (or goes to zero), the stationary probability density becomes more acute
and its maximum point approaches one. However, for large noise intensity and
harvest rate, the population size fluctuates wildly and does not grow
exponentially to the carrying capacity. So as far as biological meanings are
concerned, we must catch at small values of noise intensity and harvest rate.
Finally, we discuss the biological implications of our results.


  We revisit the topic of human genome guanine-cytosine content under neutral
evolution. For this study, the de novo mutation data within human is used to
estimate mutational rate instead of using base substitution data between
related species. We then define a new measure of mutation bias which separate
the de novo mutation counts from the background guanine-cytosine content
itself, making comparison between different datasets easier. We derive a new
formula for calculating limiting guanine-cytosine content by separating
CpG-involved mutational events as an independent variable. Using the formula
when CpG-involved mutations are considered, the guanine-cytosine content drops
less severely in the limit of neutral dynamics. We provide evidence, under
certain assumptions, that an isochore-like structure might remain as a limiting
configuration of the neutral mutational dynamics.


  In response to the COVID-19 pandemic, governments have implemented policies
to curb the spread of the novel virus. Little is known about how these policies
impact various groups in society. This paper explores the relationship between
social distancing policies, traffic volumes and air quality and how they impact
various socioeconomic groups. This study aims to understand how disparate
communities respond to Stay-at-Home Orders and other social distancing policies
to understand how human behavior in response to policy may play a part in the
prevalence of COVID-19 positive cases. We collected data on traffic density,
air quality, socio-economic status, and positive cases rates of COVID-19 for
each zip code of Salt Lake County, Utah (USA) between February 17 and June 12,
2020. We studied the impact of social distancing policies across three periods
of policy implementation. We found that wealthier and whiter zip codes
experienced a greater reduction in traffic and air pollution during the
Stay-at-Home period. However, air quality did not necessarily follow traffic
volumes in every case due to the complexity of interactions between emissions
and meteorology. We also found a strong relationship between lower
socioeconomic status and positive COVID-19 rates. This study provides initial
evidence for social distancing's effectiveness in limiting the spread of
COVID-19, while providing insight into how socioeconomic status has compounded
vulnerability during this crisis. Behavior restrictions disproportionately
benefit whiter and wealthier communities both through protection from spread of
COVID-19 and reduction in air pollution. Such findings may be further
compounded by the impacts of air pollution, which likely exacerbate COVID-19
transmission and mortality rates. Policy makers need to consider adapting
social distancing policies to maximize equity in health protection.


  Antibodies are widely used reagents to test for expression of proteins and
other antigens. However, they might not always reliably produce results when
they do not specifically bind to the target proteins that their providers
designed them for, leading to unreliable research results. While many proposals
have been developed to deal with the problem of antibody specificity, it is
still challenging to cover the millions of antibodies that are available to
researchers. In this study, we investigate the feasibility of automatically
generating alerts to users of problematic antibodies by extracting statements
about antibody specificity reported in the literature. The extracted alerts can
be used to construct an "Antibody Watch" knowledge base containing supporting
statements of problematic antibodies. We developed a deep neural network system
and tested its performance with a corpus of more than two thousand articles
that reported uses of antibodies. We divided the problem into two tasks. Given
an input article, the first task is to identify snippets about antibody
specificity and classify if the snippets report that any antibody exhibits
non-specificity, and thus is problematic. The second task is to link each of
these snippets to one or more antibodies mentioned in the snippet. The
experimental evaluation shows that our system can accurately perform both
classification and linking tasks with weighted F-scores over 0.925 and 0.923,
respectively, and 0.914 overall when combined to complete the joint task. We
leveraged Research Resource Identifiers (RRID) to precisely identify antibodies
linked to the extracted specificity snippets. The result shows that it is
feasible to construct a reliable knowledge base about problematic antibodies by
text mining.


  The spreading dynamics of an epidemic and the collective behavioral pattern
of the population over which it spreads are deeply intertwined and the latter
can critically shape the outcome of the former. Motivated by this, we design a
parsimonious game-theoretic behavioral--epidemic model, in which an interplay
of realistic factors shapes the co-evolution of individual decision-making and
epidemics on a network. Although such a co-evolution is deeply intertwined in
the real-world, existing models schematize population behavior as
instantaneously reactive, thus being unable to capture human behavior in the
long term. Our model offers a unified framework to model and predict complex
emergent phenomena, including successful collective responses, periodic
oscillations, and resurgent epidemic outbreaks. The framework also allows to
assess the effectiveness of different policy interventions on ensuring a
collective response that successfully eradicates the outbreak. Two case
studies, inspired by real-world diseases, are presented to illustrate the
potentialities of the proposed model.


  This paper introduces a mathematical framework for determining second surge
behavior of COVID-19 cases in the United States. Within this framework, a
flexible algorithmic approach selects a set of turning points for each state,
computes distances between them, and determines whether each state is in (or
over) a first or second surge. Then, appropriate distances between normalized
time series are used to further analyze the relationships between case
trajectories on a month-by-month basis. Our algorithm shows that 31 states are
experiencing second surges, while 4 of the 10 largest states are still in their
first surge, with case counts that have never decreased. This analysis can aid
in highlighting the most and least successful state responses to COVID-19.


  This study reports the development of a pattern recognition search engine for
a World Wide Web-based database of gas chromatography-electron impact mass
spectra (GC-EIMS) of partially methylated Alditol Acetates (PMAAs). Here, we
also report comparative results for two pattern recognition techniques that
were employed for this study. The first technique is a statistical technique
using Bayesian classifiers and Parzen density estimators. The second technique
involves an artificial neural network module trained with reinforcement
learning. We demonstrate here that both systems perform well in identifying
spectra with small amounts of noise. Both system's performance degrades with
degrading signal-to-noise ratio (SNR). When dealing with partial spectra
(missing data), the artificial neural network system performs better. The
developed system is implemented on the world wide web, and is intended to
identify PMAAs using submitted spectra of these molecules recorded on any
GC-EIMS instrument. The system, therefore, is insensitive to instrument and
column dependent variations in GC-EIMS spectra.


  In this paper we find common features between the equations that are used for
pandemic or epidemic modeling and the renormalization group equations that are
used in high energy physics. Some of these features include the relation of
contact matrices in pandemic modeling and operator mixing in the
renormalization group equations. Another common feature are the use of flow
diagrams and the study of fixed points both in pandemic modeling and in
evolution under renormalization group equations. We illustrate these relations
through the study of some cases of interest to the current COVID-19 pandemic.
These include pandemic modeling with mixing between different age groups and
also contact matrices associated with contact between countries. For the final
example we study the effect on mortality of waning from nonspecific vaccines
which are designed to combat different pathogens but nevertheless may lessen
the severity and mortality of COVID-19 infections.


  Currently, Brazil has one of the fastest increasing COVID-19 epidemics in the
world, that has caused at least 94 thousand confirmed deaths until now. The
city of Sao Paulo is particularly vulnerable because it is the most populous in
the country. Analyzing the spatiotemporal dynamics of COVID-19 is important to
help the urgent need to integrate better actions to face the pandemic. Thus,
this study aimed to analyze the COVID-19 mortality, from March to July 2020,
considering the spatio-time architectures, the socio-economic context of the
population, and using a fine granular level, in the most populous city in
Brazil. For this, we conducted an ecological study, using secondary public data
from the mortality information system. We describe mortality rates for each
epidemiological week and the entire period by sex and age. We modelled the
deaths using spatiotemporal and spatial architectures and Poisson probability
distributions in a latent Gaussian Bayesian model approach. We obtained the
relative risks for temporal and spatiotemporal trends and socio-economic
conditions. To reduce possible sub notification, we considered the confirmed
and suspected deaths. Our findings showed an apparent stabilization of the
temporal trend, at the end of the period, but that may change in the future.
Mortality rate increased with increasing age and was higher in men. The risk of
death was greater in areas with the worst social conditions throughout the
study period. However, this was not a uniform pattern over time, since we
identified a shift from the high risk in the areas with best socio-economic
conditions to the worst ones. Our study contributed by emphasizing the
importance of geographic screening in areas with a higher risk of death, and,
currently, worse socio-economic contexts, as a crucial aspect to reducing
disease mortality and health inequities, through integrated public health
actions.


  Automated identification of protein conformational states from simulation of
an ensemble of structures is a hard problem because it requires teaching a
computer to recognize shapes. We adapt the naive Bayes classifier from the
machine learning community for use on atom-to-atom pairwise contacts. The
result is an unsupervised learning algorithm that samples a `distribution' over
potential classification schemes. We apply the classifier to a series of test
structures and one real protein, showing that it identifies the conformational
transition with > 95% accuracy in most cases. A nontrivial feature of our
adaptation is a new connection to information entropy that allows us to vary
the level of structural detail without spoiling the categorization. This is
confirmed by comparing results as the number of atoms and time-samples are
varied over 1.5 orders of magnitude. Further, the method's derivation from
Bayesian analysis on the set of inter-atomic contacts makes it easy to
understand and extend to more complex cases.


  A dynamic model of non-lineal time-dependent ordinary differential equations
(ODE) has been applied to the interactions of a HIV infection with the immune
system cells. This model has been simplified into two compartments: lymph node
and peripheral blood. The model includes CD4 T-lymphocytes in several states
(quiescent Q, naive N and activated T), cytotoxic CD8 T-cells, B-cells and
dendritic cells. Cytokines and immunoglobulins specific for each antigen (i.e.
gp41 or p24) have been also included in the model, modelling the atraction
effect of CD4 T-cells to the infected area and the reduction of virus
concentration by immunoglobulins. HIV virus infection of CD4 T-lymphocytes is
modelled in several stages: before fusion as HIV-attached (H) and after fusion
as non-permissive / abortively infected (M), and permissive / latently infected
(L) and permissive / actively infected (I). These equations have been
implemented in a C++/Python interface application, called Immune System app,
which runs Open Modelica software to solve the ODE system through a 4th order
Runge-Kutta numerical approximation. Results of the simulation show that
although HIV virus concentration in both compartments is lower than $10^{-10}$
virus/$\mu L$ after t=2 years, quiescent lymphocytes reach an equilibrium with
a concentration lower than the initial conditions, due to the latency state,
which serves as a reservoir in time of virus production. As a conclusion, this
model can provide reliable results in other conditions, such as antiviral
therapies.


  In a given country, the cumulative death toll of the first wave of the
COVID-19 epidemic follows a sigmoid curve as a function of time. In most cases,
the curve is well described by the Gompertz function, which is characterized by
two essential parameters, the initial growth rate and the decay rate as the
first epidemic wave subsides. These parameters are determined by socioeconomic
factors and the countermeasures to halt the epidemic. The Gompertz model
implies that the total death toll depends exponentially, and hence very
sensitively, on the ratio between these rates. The remarkably different
epidemic curves for the first epidemic wave in Sweden and Norway and many other
countries are classified and discussed in this framework, and their usefulness
for the planning of mitigation strategies is discussed.


  Some recent works reveal that there are models of differential equations for
the mean and variance of infected individuals that reproduce the SIS epidemic
model at some point. This stochastic SIS epidemic model can be interpreted as a
Hamiltonian system, therefore we wondered if it could be geometrically handled
through the theory of Lie--Hamilton systems, and this happened to be the case.
The primordial result is that we are able to obtain a general solution for the
stochastic/ SIS-epidemic model (with fluctuations) in form of a nonlinear
superposition rule that includes particular stochastic solutions and certain
constants to be related to initial conditions of the contagion process. The
choice of these initial conditions will be crucial to display the expected
behavior of the curve of infections during the epidemic. We shall limit these
constants to nonsingular regimes and display graphics of the behavior of the
solutions. As one could expect, the increase of infected individuals follows a
sigmoid-like curve.
  Lie--Hamiltonian systems admit a quantum deformation, so does the stochastic
SIS-epidemic model. We present this generalization as well. If one wants to
study the evolution of an SIS epidemic under the influence of a constant heat
source (like centrally heated buildings), one can make use of quantum
stochastic differential equations coming from the so-called quantum
deformation.


  Broadly neutralizing antibodies are promising candidates for treatment and
prevention of HIV-1 infections. Such antibodies can temporarily suppress viral
load in infected individuals; however, the virus often rebounds by escape
mutants that have evolved resistance. In this paper, we map an in vivo fitness
landscape of HIV-1 interacting with broadly neutralizing antibodies, using data
from a recent clinical trial. We identify two fitness factors, antibody dosage
and viral load, that determine viral reproduction rates reproducibly across
different hosts. The model successfully predicts the escape dynamics of HIV-1
in the course of an antibody treatment, including a characteristic frequency
turnover between sensitive and resistant strains. This turnover is governed by
a dosage-dependent fitness ranking, resulting from an evolutionary tradeoff
between antibody resistance and its collateral cost in drug-free growth. Our
analysis suggests resistance-cost tradeoff curves as a measure of antibody
performance in the presence of resistance evolution.


  We investigate the classical two species ODE and PDE Lotka-Volterra
competition models, where one of the competitors could potentially go extinct
in finite time. We show that in this setting, classical theories and intuitions
do not hold, and various counter intuitive dynamics are possible. In
particular, the weaker competitor could avoid competitive exclusion, and the
slower diffuser may not win. Numerical simulations are performed to verify our
analytical findings.


  In marine management, fish stocks are often managed on a stock-by-stock basis
using single-species models. Many of these models are based upon statistical
techniques and are good at assessing the current state and making short-term
predictions; however, as they do not model interactions between stocks, they
lack predictive power on longer timescales. Additionally, there are mechanistic
multi-species models that represent key biological processes and consider
interactions between stocks such as predation and competition for resources.
Due to the complexity of these models, they are difficult to fit to data, and
so many mechanistic multi-species models depend upon single-species models
where they exist, or ad hoc assumptions when they don't, for parameters such as
annual fishing mortality. In this paper we demonstrate that by taking a
state-space approach, many of the uncertain parameters can be treated
dynamically, allowing us to fit, with quantifiable uncertainty, mechanistic
multi-species models directly to data. We demonstrate this by fitting uncertain
parameters, including annual fishing mortality, of a size-based multi-species
model of the Celtic Sea, for species with and without single-species
stock-assessments. Consequently, errors in the single-species models no longer
propagate through the multi-species model and underlying assumptions are more
transparent. Building mechanistic multi-species models that are internally
consistent, with quantifiable uncertainty, will improve their credibility and
utility for management. This may lead to their uptake by being either used to
corroborate single-species models; directly in the advice process to make
predictions into the future; or used to provide a new way of managing
data-limited stocks.


  Population survival depends on a large set of factors that includes
environment structure. Due to landscape heterogeneity, species can occupy
particular regions that provide the ideal scenario for development, working as
a refuge from harmful environmental conditions. Survival occurs if population
growth overcomes the losses caused by adventurous individuals that cross the
patch edge. In this work, we consider a single species dynamics in a bounded
domain with a space-dependent diffusion coefficient. We investigate the impact
of heterogeneous diffusion on the minimal patch size that allows population
survival and show that, typically, this critical size is smaller than the one
for a homogeneous medium with the same average diffusivity.


  Liquid-liquid phase separation is the mechanism underlying the formation of
biomolecular condensates. Disordered protein regions often drive phase
separation, but molecular interactions of disordered protein regions are not
well understood, sometimes leading to the conflation that all disordered
protein regions drive phase separation. Given the critical role of phase
separation in many cellular processes, and that dysfunction of phase separation
can lead to debilitating diseases, it is important that we understand the
interactions and sequence properties underlying phase behavior. A conceptual
framework that divides IDRs into interacting and solvating regions has proven
particularly useful, and analytical instantiations and coarse-grained models
can test our understanding of the driving forces against experimental phase
behavior. Validated simulation paradigms enable the exploration of sequence
space to help our understanding of how disordered protein regions can encode
phase behavior, which IDRs may mediate phase separation in cells, and which
IDRs are in contrast highly soluble.


  Alchemical free energy calculations are a useful tool for predicting free
energy differences associated with the transfer of molecules from one
environment to another. The hallmark of these methods is the use of "bridging"
potential energy functions representing \emph{alchemical} intermediate states
that cannot exist as real chemical species. The data collected from these
bridging alchemical thermodynamic states allows the efficient computation of
transfer free energies (or differences in transfer free energies) with orders
of magnitude less simulation time than simulating the transfer process
directly. While these methods are highly flexible, care must be taken in
avoiding common pitfalls to ensure that computed free energy differences can be
robust and reproducible for the chosen force field, and that appropriate
corrections are included to permit direct comparison with experimental data. In
this paper, we review current best practices for several popular application
domains of alchemical free energy calculations, including relative and absolute
small molecule binding free energy calculations to biomolecular targets.


  In this article, we give visualization and different machine learning
technics for two weeks and 40 days ahead forecast based on public data. On July
15, 2020, Senegal reopened its airspace doors, while the number of confirmed
cases is still increasing. The population no longer respects hygiene measures,
social distancing as at the beginning of the contamination. Negligence or
tiredness to always wear the masks? We make forecasting on the inflection point
and possible ending time.


  We present two different approaches for modeling the spread of the COVID-19
pandemic. Both approaches are based on the population classes susceptible,
exposed, infectious, quarantined, and recovered and allow for an arbitrary
number of subgroups with different infection rates and different levels of
testing. The first model is derived from a set of ordinary differential
equations that incorporates the rates at which population transitions take
place among classes. The other is a particle model, which is a specific case of
crowd simulation model, in which the disease is transmitted through particle
collisions and infection rates are varied by adjusting the particle velocities.
The parameters of these two models are tuned using information on COVID-19 from
the literature and country-specific data, including the effect of restrictions
as they were imposed and lifted. We demonstrate the applicability of both
models using data from Cyprus, for which we find that both models yield very
similar results, giving confidence in the predictions.


  This work constructs, analyzes, and simulates a new compartmental SEIR-type
model for the dynamics and potential control of the current COVID-19 pandemic.
The novelty in this work is two-fold. First, the population is divided
according to its compliance with disease control directives (lockdown,
shelter-in-place, masks/face coverings, physical distancing, etc.) into those
who fully comply and those who follow the directives partially, or are
necessarily mobile (such as medical staff). This split, indirectly, reflects on
the quality and consistency of these measures. This allows the assessment of
the overall effectiveness of the control measures and the impact of their
relaxing or tightening on the disease spread. Second, the adequate contact
rate, which directly affects the infection rate, is one of the model unknowns,
as it keeps track of the changes in the population behavior and the
effectiveness of various disease treatment modalities via a differential
inclusion. Existence, uniqueness and positivity results are proved using a
nonstandard convex analysis-based approach. As a case study, the pandemic
outbreak in the Republic of Korea (South Korea) is simulated. The model
parameters were found by minimizing the deviation of the model prediction from
the reported data over the first 100 days of the pandemic in South Korea.The
simulations show that the model captures accurately the pandemic dynamics in
the subsequent 75 days, which provides confidence in the model predictions and
its future use. In particular, the model predicts that about 40% of the
infections were not documented, which implies that asymptomatic infections
contribute silently but substantially to the spread of the disease indicating
that more widespread asymptomatic testing is necessary.


  In this paper, we develop a probabilistic mathematical model for the spread
of coronavirus disease (COVID-19). It takes into account the known special
characteristics of this disease such as the existence of infectious undetected
cases and the different social and infectiousness conditions of infected
people. In particular, it considers the social structure and governmental
measures in a country, the fraction of detected cases over the real total
infected cases, and the influx of undetected infected people from outside the
borders. Although the model is simple and allows a reasonable identification of
its parameters, using the data provided by local authorities on this pandemic,
it is also complex enough to capture the most important effects. We study the
particular case of Lebanon and use its reported data to estimate the model
parameters, which can be of interest for estimating the spread of COVID-19 in
other countries. We show a good agreement between the reported data and the
estimations given by our model. We also simulate several scenarios that help
policy makers in deciding how to loosen different measures without risking a
severe wave of COVID-19. We are also able to identify the main factors that
lead to specific scenarios which helps in a better understanding of the spread
of the virus.


  Using a simple economic model in which social-distancing reduces contagion,
we study the implications of waning immunity for the epidemiological dynamics
and social activity. If immunity wanes, we find that COVID-19 likely becomes
endemic and that social-distancing is here to stay until the discovery of a
vaccine or cure. But waning immunity does not necessarily change optimal
actions on the onset of the pandemic. Decentralized equilibria are virtually
independent of waning immunity until close to peak infections. For centralized
equilibria, the relevance of waning immunity decreases in the probability of
finding a vaccine or cure, the costs of infection (e.g., infection-fatality
rate), and the presence of other NPIs that lower contagion (e.g., quarantining
and mask use). In simulations calibrated to July 2020, our model suggests that
waning immunity is virtually unimportant for centralized equilibria until at
least 2021. This provides vital time for individuals and policymakers to learn
about immunity against SARS-CoV-2 before it becomes critical.


  We adapt the hybrid mechanistic-statistical approach of Ref. [1] to estimate
the total number of undocumented Covid-19 infections in Cuba. This scheme is
based on the maximum likelihood estimation of a SIR-like model parameters for
the infected population, assuming that the detection process matches a
Bernoulli trial. Our estimations show that (a) 60% of the infections were
undocumented, (b) the real epidemics behind the data peaked ten days before the
reports suggested, and (c) the reproduction number swiftly vanishes after 80
epidemic days.


  Bacterial quorum sensing is the communication that takes place between
bacteria as they secrete certain molecules into the intercellular medium that
later get absorbed by the secreting cells themselves and by others. Depending
on cell density, this uptake has the potential to alter gene expression and
thereby affect global properties of the community. We consider the case of
multiple bacterial species coexisting, referring to each one of them as a
genotype and adopting the usual denomination of the molecules they collectively
secrete as public goods. A crucial problem in this setting is characterizing
the coevolution of genotypes as some of them secrete public goods (and pay the
associated metabolic costs) while others do not but may nevertheless benefit
from the available public goods. We introduce a network model to describe
genotype interaction and evolution when genotype fitness depends on the
production and uptake of public goods. The model comprises a random graph to
summarize the possible evolutionary pathways the genotypes may take as they
interact genetically with one another, and a system of coupled differential
equations to characterize the behavior of genotype abundance in time. We study
some simple variations of the model analytically and more complex variations
computationally. Our results point to a simple trade-off affecting the
long-term survival of those genotypes that do produce public goods. This
trade-off involves, on the producer side, the impact of producing and that of
absorbing the public good. On the non-producer side, it involves the impact of
absorbing the public good as well, now compounded by the molecular
compatibility between the producer and the non-producer. Depending on how these
factors turn out, producers may or may not survive.


  One of the key steps in building deep learning systems for drug
classification and generation is the choice of featurization for the molecules.
Previous featurization methods have included molecular images, binary strings,
graphs, and SMILES strings. This paper proposes the creation of molecular
images captioned with binary vectors that encode information not contained in
or easily understood from a molecular image alone. Specifically, we use Morgan
fingerprints, which encode higher level structural information, and MACCS keys,
which encode yes or no questions about a molecules properties and structure. We
tested our method on the HIV dataset published by the Pande lab, which consists
of 41,127 molecules labeled by if they inhibit the HIV virus. Our final model
achieved a state of the art AUC ROC on the HIV dataset, outperforming all other
methods. Moreover, the model converged significantly faster than most other
methods, requiring dramatically less computational power than unaugmented
images.


  One of the central goals in precision health is the understanding and
interpretation of high-dimensional biological data to identify genes and
markers associated with disease initiation, development, and outcomes. Though
significant effort has been committed to harness gene expression data for
multiple analyses while accounting for time-to-event modeling by including
survival times, many traditional analyses have focused separately on
non-negative matrix factorization (NMF) of the gene expression data matrix and
survival regression with Cox proportional hazards model. In this work, Cox
proportional hazards regression is integrated with NMF by imposing survival
constraints. This is accomplished by jointly optimizing the Frobenius norm and
partial log likelihood for events such as death or relapse. Simulation results
on synthetic data demonstrated the superiority of the proposed method, when
compared to other algorithms, in finding survival associated gene clusters. In
addition, using human cancer gene expression data, the proposed technique can
unravel critical clusters of cancer genes. The discovered gene clusters reflect
rich biological implications and can help identify survival-related biomarkers.
Towards the goal of precision health and cancer treatments, the proposed
algorithm can help understand and interpret high-dimensional heterogeneous
genomics data with accurate identification of survival-associated gene
clusters.


  For three typical sets of small reaction networks (networks with two
reactions, one irreversible and one reversible reaction, or two
reversible-reaction pairs), we completely answer the challenging question: what
is the smallest subset of all multistable networks such that any multistable
network outside of the subset contains either more species or more reactants
than any network in this subset?


  Contact Tracing (CT) is one of the measures taken by government and health
officials to mitigate the spread of the novel coronavirus. In this paper, we
investigate its efficacy by developing a compartmental model for assessing its
impact on mitigating the spread of the virus. We describe the impact on the
reproduction number ${\mathcal{R}_0}$ of COVID-19. In particular, we discuss
the importance and relevance of parameters of the model such as the number of
reported cases, effectiveness of tracking and monitoring policy, and the
transmission rates to contact tracing. We describe the terms "perfect
tracking", "perfect monitoring" and "perfect reporting" to indicate that traced
contacts will be tracked while incubating, tracked contacts are efficiently
monitored so that they do not cause secondary infections, and all infected
persons are reported, respectively. We consider three special scenarios: (1)
perfect monitoring and perfect tracking of contacts of a reported case, (2)
perfect reporting of cases and perfect monitoring of tracked reported cases and
(3) perfect reporting and perfect tracking of contacts of reported cases.
Furthermore, we gave a lower bound on the proportion of contacts to be traced
to ensure that the effective reproduction, ${\mathcal{R}_c}$, is below one and
describe ${\mathcal{R}_c}$ in terms of observable quantities such as the
proportion of reported and traced cases. Model simulations using the COVID-19
data obtained from John Hopkins University for some selected states in the US
suggest that even late intervention of CT may reasonably reduce the
transmission of COVID-19 and reduce peak hospitalizations and deaths.


  This study investigates whether coupling crop modeling and machine learning
(ML) improves corn yield predictions in the US Corn Belt. The main objectives
are to explore whether a hybrid approach (crop modeling + ML) would result in
better predictions, investigate which combinations of hybrid models provide the
most accurate predictions, and determine the features from the crop modeling
that are most effective to be integrated with ML for corn yield prediction.
Five ML models (linear regression, LASSO, LightGBM, random forest, and XGBoost)
and six ensemble models have been designed to address the research question.
The results suggest that adding simulation crop model variables (APSIM) as
input features to ML models can decrease yield prediction root mean squared
error (RMSE) from 7 to 20%. Furthermore, we investigated partial inclusion of
APSIM features in the ML prediction models and we found soil moisture related
APSIM variables are most influential on the ML predictions followed by
crop-related and phenology-related variables. Finally, based on feature
importance measure, it has been observed that simulated APSIM average drought
stress and average water table depth during the growing season are the most
important APSIM inputs to ML. This result indicates that weather information
alone is not sufficient and ML models need more hydrological inputs to make
improved yield predictions.


  Pyroptosis is an inflammatory mode of cell death that can contribute to the
cytokine storm associated with severe cases of coronavirus disease 2019
(COVID-19). The formation of the NLRP3 inflammasome is central to pyroptosis,
which may be induced by severe acute respiratory syndrome coronavirus 2
(SARS-CoV-2). Inflammasome formation, and by extension pyroptosis, may be
inhibited by certain anti-inflammatory drugs. In this study, we present a
single-cell mathematical model that captures the formation of the NLRP3
inflammasome, pyroptotic cell death and responses to anti-inflammatory
intervention that hinder the formation of the NLRP3 inflammasome. The model is
formulated in terms of a system of ordinary differential equations (ODEs) that
describe the dynamics of the biological components involved in pyroptosis. Our
results demonstrate that an anti-inflammatory drug can delay the formation of
the NLRP3 inflammasome, and thus may alter the mode of cell death from
inflammatory (pyroptosis) to non-inflammatory e.g., apoptosis). The single-cell
model is being implemented in a SARS-CoV-2 Tissue Simulator, in collaboration
with a multidisciplinary coalition investigating within host-dynamics of
COVID-19. In this paper, we provide an overview of the SARS-CoV-2 Tissue
Simulator and highlight the effects of pyroptosis on a cellular level.


  Objective: The COVID-19 pandemic has created many challenges that need
immediate attention. Various epidemiological and deep learning models have been
developed to predict the COVID-19 outbreak, but all have limitations that
affect the accuracy and robustness of the predictions. Our method aims at
addressing these limitations and making earlier and more accurate pandemic
outbreak predictions by (1) using patients' EHR data from different counties
and states that encode local disease status and medical resource utilization
condition; (2) considering demographic similarity and geographical proximity
between locations; and (3) integrating pandemic transmission dynamics into deep
learning models. Materials and Methods: We proposed a spatio-temporal attention
network (STAN) for pandemic prediction. It uses an attention-based graph
convolutional network to capture geographical and temporal trends and predict
the number of cases for a fixed number of days into the future. We also
designed a physical law-based loss term for enhancing long-term prediction.
STAN was tested using both massive real-world patient data and open source
COVID-19 statistics provided by Johns Hopkins university across all U.S.
counties. Results: STAN outperforms epidemiological modeling methods such as
SIR and SEIR and deep learning models on both long-term and short-term
predictions, achieving up to 87% lower mean squared error compared to the best
baseline prediction model. Conclusions: By using information from real-world
patient data and geographical data, STAN can better capture the disease status
and medical resource utilization information and thus provides more accurate
pandemic modeling. With pandemic transmission law based regularization, STAN
also achieves good long-term prediction performance.


  Viral spread has been intermittently threatening human life over time.
Characterizing the viral concentration and modelling the viral transmission
are, therefore, considered major milestones for enhancing viral detection
capabilities. This paper addresses the problem of viral aerosol detection based
on the exhaled breath in a bounded environment, e.g., a bounded room. The paper
models the exhaled breath as a cloud which is emitted through the room
continuously, and analyzes the temporal-spatial virus concentration by
accounting for partial absorption and reflection at each side of the room. The
paper first derives a closed form expression of the temporal-spatial virus
concentration. It then considers the deployment of a receiver composed of an
air sampler and a bio-sensor to detect the viral existence of a specific virus.
We, therefore, assess the detection capabilities of the proposed system via
evaluating the viral miss-detection probability as a function of the sampling
volume and the detection time-instance at the receiver side. Our numerical
simulations verify the validity of the analytical results, and illustrate the
ability of the proposed system to detect viruses in indoor environments. The
results further characterize the impacts of several system parameters on the
miss-detection probability.


  As the ravages caused by COVID-19 pandemic are becoming inevitable with every
moment, monitoring and understanding of transmission and fatality rate has
become even more paramount for containing its spread. The key purpose of this
analysis is to report the real-time effective reproduction rate ($R_t$ ) and
case fatality rates (CFR) of COVID-19 in South Asia region. Data for this study
are extracted from JHU CSSE COVID-19 Data source up to July 31, 2020. $R_t$ is
estimated using exponential growth and time-dependent methods. R0 package in
R-language is employed to estimate $R_t$ by fitting the existing epidemic
curve. Case fatality rate is estimated by using Naive and Kaplan-Meier methods.
Owing to exponential increase in cases of COVID-19, the pandemic will ensue in
India, Maldives and in Nepal as $R_t$ was estimated greater than 1 for these
countries. Although case fatality rates are found lesser as compared to other
highly affected regions in the world, strict monitoring of deaths for better
health facilities and care of patients is emphasized. More regional level
cooperation and efforts are the need of time to minimize the detrimental
effects of the virus.


  The synthesis of complex oligosaccharides is desired for their potential as
prebiotics, and their role in the pharmaceutical and food industry.
Levansucrase (LS, EC 2.4.1.10), a fructosyl-transferase, can catalyze the
synthesis of these compounds. LS acquires a fructosyl residue from a donor
molecule and performs a non-Lenoir transfer to an acceptor molecule, via
$\beta$-(2$\to$6)-glycosidic linkages. Genome mining was used to uncover new LS
enzymes with increased transfructosylating activity and wider acceptor
promiscuity, with an initial screening revealing five LS enzymes. The product
profiles and activities of these enzymes were examined after their incubation
with sucrose. Alternate acceptor molecules were also incubated with the enzymes
to study their consumption. LSs from Gluconobacter oxydans and Novosphingobium
aromaticivorans synthesized fructooligosaccharides (FOSs) with up to 13 units
in length. Alignment of their amino acid sequences and substrate docking with
homology models identified structural elements causing differences in their
product spectra. Raffinose, over sucrose, was the preferred donor molecule for
the LS from Vibrio natriegens, N. aromaticivorans, and Paraburkolderia
graminis. The LSs examined were found to have wide acceptor promiscuity,
utilizing monosaccharides, disaccharides, and two alcohols to a high degree.


  Structure determination is key to understanding protein function at a
molecular level. Whilst significant advances have been made in predicting
structure and function from amino acid sequence, researchers must still rely on
expensive, time-consuming analytical methods to visualise detailed protein
conformation. In this study, we demonstrate that it is possible to make
accurate ($\geq$80%) predictions of protein class and architecture from
structures determined at low ($>$3A) resolution, using a deep convolutional
neural network trained on high-resolution ($\leq$3A) structures represented as
2D matrices. Thus, we provide proof of concept for high-speed, low-cost protein
structure classification at low resolution, and a basis for extension to
prediction of function. We investigate the impact of the input representation
on classification performance, showing that side-chain information may not be
necessary for fine-grained structure predictions. Finally, we confirm that
high-resolution, low-resolution and NMR-determined structures inhabit a common
feature space, and thus provide a theoretical foundation for boosting with
single-image super-resolution.


  We highlight the usefulness of city-scale agent-based simulators in studying
various non-pharmaceutical interventions to manage an evolving pandemic. We
ground our studies in the context of the COVID-19 pandemic and demonstrate the
power of the simulator via several exploratory case studies in two
metropolises, Bengaluru and Mumbai. Such tools become common-place in any city
administration's tool kit in our march towards digital health.


  This study describes the dynamics of COVID-19 deaths and infections via a
Monte Carlo approach. The analyses include death's data from USA, Brazil,
Mexico, UK, India and Russia, which comprise the four countries with the
highest number of deaths/confirmed cases, as of Aug 07, 2020, according to the
WHO. The Gompertz functions were fitted to the data of weekly averaged
confirmed deaths per day by mapping the $\chi^2$ values. The uncertainties,
variances and covariances of the model parameters were calculated by
propagation. The fitted functions for the average deaths per day for USA and
India have an upward trend, with the former having a higher growth rate and
quite huge uncertainties. For Mexico, UK and Russia, the fits are consistent
with a slope down pattern. For Brazil we found a subtle trend down, but with
significant uncertainties. The USA, UK and India data shown a first peak with a
higher growth rate when compared to the second one, demonstrating the benefits
of non-pharmaceutical interventions of sanitary measures and social distance
flattening the curve. For USA, a third peak seems quite plausible, most likely
related with the recent relaxation policies. Brazil's data are satisfactorily
described by two highly overlapped Gompertz functions with similar growth
rates, suggesting a two-steps process for the pandemic spreading. The 95% CI
for the total number of deaths ($\times 10^3$) predicted by the model for Aug
31, 2020 are 160 to 220, 110 to 130, 59 to 62, 46.6 to 47.3, 54 to 63 and 16.0
to 16.7 for USA, Brazil, Mexico, UK, India and Russia, respectively. Our
estimates for the prevalences of infections are in reasonable agreement with
some preliminary reports from serological studies carried out in USA and
Brazil. The method represents an effective framework to estimate the line-shape
of the infection curves and the uncertainties of the relevant parameters based
on the actual data.


  Standard epidemic models based on compartmental differential equations are
investigated under continuous parameter change as external forcing. We show
that seasonal modulation of the contact parameter superimposed a monotonic
decay needs a different description than that of the standard chaotic dynamics.
The concept of snapshot attractors and their natural probability distribution
has been adopted from the field of the latest climate-change-research to show
the importance of transient effect and ensemble interpretation of disease
spread. After presenting the extended bifurcation diagram of measles, the
temporal change of the phase space structure is investigated. By defining
statistical measures over the ensemble, we can interpret the internal
variability of the epidemic as the onset of complex dynamics even for those
values of contact parameter where regular behavior is expected. We argue that
anomalous outbreaks of infectious class cannot die out until transient chaos is
presented for various parameters. More important, that this fact becomes
visible by using of ensemble approach rather than single trajectory
representation. These findings are applicable generally in nonlinear dynamical
systems such as standard epidemic models regardless of parameter values.


  A general stochastic model for susceptible -> infective -> recovered (SIR)
epidemics in non homogeneous populations is considered. The heterogeneity is a
very important aspect here since it allows more realistic but also more complex
models. The basic reproduction number $R_0$, an indication of the probability
of an outbreak for homogeneous populations does not indicate the probability of
an outbreak for non homogeneous models anymore, because it changes with the
initially infected case. Therefore, we use "individual $R_0$" that is the
expected number of secondary cases for a given initially infected individual.
Thus, the effectiveness of intervention strategies can be assessed by their
capability to reduce individual $R_0$ values. Also an intelligent vaccination
plan for fully heterogeneous populations is proposed. It is based on the
recursive calculation of individual R0 values.


  In this Ansatz we consider theoretical constructions of RNA polymers into
automata, a form of computational structure. The basis for transitions in our
automata are plausible RNA enzymes that may perform ligation or cleavage.
Limited to these operations, we construct RNA automata of increasing
complexity; from the Finite Automaton (RNA-FA) to the Turing Machine equivalent
2-stack PDA (RNA-2PDA) and the universal RNA-UPDA. For each automaton we show
how the enzymatic reactions match the logical operations of the RNA automaton.
A critical theme of the Ansatz is the self-reference in RNA automata
configurations which exploits the program-data duality but results in
computational undecidability. We describe how computational undecidability is
exemplified in the self-referential Liar paradox that places a boundary on a
logical system, and by construction, any RNA automata. We argue that an
expansion of the evolutionary space for RNA-2PDA automata can be interpreted as
a hierarchical resolution of computational undecidability by a meta-system
(akin to Turing's oracle), in a continual process analogous to Turing's ordinal
logics and Post's extensible recursively generated logics. On this basis, we
put forward the hypothesis that the resolution of undecidable configurations in
RNA automata represent a novelty generation mechanism and propose avenues for
future investigation of biological automata.


  Currently, the number of patients with COVID-19 has significantly increased.
Thus, there is an urgent need for developing treatments for COVID-19. Drug
repurposing, which is the process of reusing already-approved drugs for new
medical conditions, can be a good way to solve this problem quickly and
broadly. Many clinical trials for COVID-19 patients using treatments for other
diseases have already been in place or will be performed at clinical sites in
the near future. Additionally, patients with comorbidities such as diabetes
mellitus, obesity, liver cirrhosis, kidney diseases, hypertension, and asthma
are at higher risk for severe illness from COVID-19. Thus, the relationship of
comorbidity disease with COVID-19 may help to find repurposable drugs. To
reduce trial and error in finding treatments for COVID-19, we propose building
a network-based drug repurposing framework to prioritize repurposable drugs.
First, we utilized knowledge of COVID-19 to construct a disease-gene-drug
network (DGDr-Net) representing a COVID-19-centric interactome with components
for diseases, genes, and drugs. DGDr-Net consisted of 592 diseases, 26,681
human genes and 2,173 drugs, and medical information for 18 common
comorbidities. The DGDr-Net recommended candidate repurposable drugs for
COVID-19 through network reinforcement driven scoring algorithms. The scoring
algorithms determined the priority of recommendations by utilizing graph-based
semi-supervised learning. From the predicted scores, we recommended 30 drugs,
including dexamethasone, resveratrol, methotrexate, indomethacin, quercetin,
etc., as repurposable drugs for COVID-19, and the results were verified with
drugs that have been under clinical trials. The list of drugs via a data-driven
computational approach could help reduce trial-and-error in finding treatment
for COVID-19.


  We prove that stochastic replicator dynamics can be interpreted as intrinsic
Brownian motion on the simplex equipped the Aitchison geometry. As an immediate
consequence we derive three approximation results in the spirit of Wong-Zakai
approximation, Donsker's invariance principle and a JKO-scheme. Finally, using
the Fokker-Planck equation and Wasserstein-contraction estimates, we study the
long time behavior of the stochastic replicator equation, as an example of a
non-gradient drift diffusion on the Aitchison simplex.


  Measuring the distance between two bacterial genomes under the inversion
process is usually done by assuming all inversions to occur with equal
probability. Recently, an approach to calculating inversion distance using
group theory was introduced, and is effective for the model in which only very
short inversions occur. In this paper, we show how to use the group-theoretic
framework to establish minimal distance for any weighting on the set of
inversions, generalizing previous approaches. To do this we use the theory of
rewriting systems for groups, and exploit the Knuth--Bendix algorithm, the
first time this theory has been introduced into genome rearrangement problems.
  The central idea of the approach is to use existing group theoretic methods
to find an initial path between two genomes in genome space (for instance using
only short inversions), and then to deform this path to optimality using a
confluent system of rewriting rules generated by the Knuth--Bendix algorithm.


  Chlamydia trachomatis is the most commonly reported infectious disease in the
United States and causes important reproductive morbidity in women. The Centers
for Disease Control and Prevention have recommended routine screening of
sexually active women under age 25 but have not recommended screening among
men. Consequently, untested and untreated men may serve as a reservoir of
infection in women. Despite three decades of screening women, the chlamydia
prevalence has continued to increase. Moreover, chlamydia is five times more
common in African American (AA) youth compared to Whites, constituting an
important health disparity. The Check It program is a bundled Ct intervention
targeting AA men aged 15-24 who have sex with women. We created an
individual-based network model to simulate a realistic chlamydia epidemic on
sexual contact networks for the target population. Based on the practice in
Check It, we quantified the impact of screening young AA men on the chlamydia
prevalence in women. We used sensitivity analysis to quantify the relative
importance of each Check It intervention component, and the significance ranked
from high to low was venue-based screening, expedited index treatment,
expedited partner treatment, rescreening. We estimated that by annually
screening 7.5% of the target male population, the chlamydia prevalence would be
reduced by 8.1% and 8.8% in men and women, respectively. The findings suggested
that male-screening has the potential to significantly reduce the prevalence
among women.


  Newly-translated glycoproteins in the endoplasmic reticulum (ER) often
undergo cycles of chaperone binding and release in order to assist in folding.
Quality control is required to distinguish between proteins that have completed
native folding, those that have yet to fold, and those that have misfolded.
Using quantitative modeling, we explore how the design of the quality-control
pathway modulates its efficiency. Our results show that an energy-consuming
cyclic quality-control process, similar to the observed physiological system,
outperforms alternative designs. The kinetic parameters that optimize the
performance of this system drastically change with protein production levels,
while remaining relatively insensitive to the protein folding rate. Adjusting
only the degradation rate, while fixing other parameters, allows the pathway to
adapt across a range of protein production levels, aligning with in vivo
measurements that implicate the release of degradation-associated enzymes as a
rapid-response system for perturbations in protein homeostasis. The
quantitative models developed here elucidate design principles for effective
glycoprotein quality control in the ER, improving our mechanistic understanding
of a system crucial to maintaining cellular health.


  Evolution can be broadly described in terms of mutations of the genotype and
the subsequent selection of the phenotype. The full enumeration of a given
genotype-phenotype (GP) map is therefore a powerful technique in examining
evolutionary landscapes. However, because the number of genotypes typically
grows exponentially with genome length, such calculations rapidly become
intractable. Here I apply graphics processing unit(GPU) techniques to the
hydrophobic-polar (HP)model for protein folding. This GP map is a simple and
well-studied model for the complex process of protein folding. Prior studies on
relatively small 2D and 3D lattices have been exclusively carried out using
conventional central processing unit (CPU) approaches. By using GPU techniques,
I was able to reproduce the pioneering calculations of Li et al.[1] with a
speed up of 580-700 fold over a CPU. I was also able to perform the largest
enumeration to date of the 6x6 lattice. These novel calculations provide
evidence that a popular "plum-pudding" metaphor that suggests that phenotypes
are disconnected in genotype space does not describe the data. Instead a
"spaghetti" metaphor of connected genotype networks may be more suitable.
Furthermore, the data allows the relationships between designability and
complexity within GP space to be explored. GPU approaches appear extremely well
suited toGP mapping and the success of this work provides a promising
introduction for its wider application in this field.


  In this work, we developed a deep learning model-based approach to forecast
the spreading trend of SARS-CoV-2 in the United States. We implemented the
designed model using the United States to confirm cases and state demographic
data and achieved promising trend prediction results. The model incorporates
demographic information and epidemic time-series data through a Gated Recurrent
Unit structure. The identification of dominating demographic factors is
delivered in the end.


  Since the 1960s, Democrats and Republicans in U.S. Congress have taken
increasingly polarized positions, while the public's policy positions have
remained centrist and moderate. We explain this apparent contradiction by
developing a dynamical model that predicts ideological positions of political
parties. Our approach tackles the challenge of incorporating bounded
rationality into mathematical models and integrates the empirical finding of
satisficing decision making---voters settle for candidates who are "good
enough" when deciding for whom to vote. We test the model using data from the
U.S. Congress over the past 150 years, and find that our predictions are
consistent with the two major political parties' historical trajectory. In
particular, the model explains how polarization between the Democrats and
Republicans since the 1960s could be a consequence of increasing ideological
homogeneity within the parties.


  To improve understanding of the physiology, ecology and regeneration
requirements of black poplar (Populus nigra), a severely endangered tree
species in Europe, we analysed the biochemical composition of seeds from 20
families derived from open pollination of 20 trees in two contrasting
environments in France, along the Dr\^ome and Loire rivers. Significant
between-family differences were detected. Total lipid contents differed
qualitatively between seed of the two populations and were higher in Dr\^ome
population (214 more or less 29.8 versus 172.72 more or less 19.77 mg.g\-1 dw).
Soluble sugars were less abundant in seeds of the Dr\^ome population (78.78
more or less 10.01 versus 104.10 more or less 25.32 mg.g\-1 dw). The protein
content (322 more or less 74.3 mg.g\-1 dw) and quantity of reserve proteins did
not differ between seeds of the two populations. Measurements in three
consecutive years showed that seeds of the Loire population had significantly
higher dry weight (0.69 more or less 0.10 versus 0.45 more or less 0.10 mg.g\-1
dw), but this did not significantly affect the germination rate or root growth
of the seedlings measured 5 days after sowing. However, one group of Loire
seedlings had substantially longer juvenile roots than another group. Overall,
the results suggest that intrinsic qualities of the seeds are not critical
factors in early stages and in selection pressure of the species life cycle and
recruitment.


  Chlamydia trachomatis (Ct) is the most reported sexually transmitted
infection in the United States with a major cause of infertility and pelvic
inflammatory disease among women. Despite decades of screening women for Ct,
rates increase among young African Americans (AA). We create and analyze an
agent-based network model to understand the spread of Ct. We calibrate the
model parameters to agree with survey data showing Ct prevalence of 12% of the
women and 10% of the men in the 15-25 year-old AA in New Orleans, Louisiana.
Our model accounts for long-term and casual partnerships. The network captures
assortative mixing of individuals by preserving the joint-degree distributions
observed in the data. We compare the efficiency of intervention strategies of
randomly screening men, partner notification, which includes partner treatment,
partner screening, and rescreening for infection. We compare the difference
between treating partners of an infected person both with and without testing
them. We observe that although increased Ct screening, rescreening and treating
most of the partners of infected people will reduce the prevalence, these
mitigations alone are not sufficient to control the epidemic. The current
practice is to treat the partners of an infected individual, without first
testing them for infection. The model predicts that if a sufficient number of
the partners of all infected people are tested and treated, then there is a
threshold condition where the epidemic can be mitigated. This threshold results
from the expanded treatment network created by treating the partners of the
infected partners of an individual. Although these conclusions can help design
future Ct mitigation studies, we caution the reader that these conclusions are
for the mathematical model, not the real world, and are contingent on the
validity of the model assumptions.


  In this study, we present a new epidemiological model, with contamination
from confirmed and unreported. We also compute equilibria and study their
stability without intervention strategies. Optimal control theory has proven to
be a successful tool in understanding ways to curtail the spread of infectious
diseases by devising the optimal disease intervention strategies. We
investigate the impact of distancing, case finding, and case holding controls
while at the same time, we minimize the number of infected and dead
individuals. The method consists of minimizing the cost functional related to
infectious, death, and controls through some strategies to reduce the spread of
the COVID19 epidemic.


  The poly-tRNA model was recently presented for the origin and evolution of
genetic coding. This model has led to a rather precise description of what
might have occurred at the beginning of protein synthesis in the first life
form. Here, we further discuss some interesting implications of this model.
First, the system of encoded peptide/protein synthesis appears to have started
and developed on the breeding ground of a rich RNA world, responsible for the
infancy of life existence and complexity. Furthermore, once protein synthesis
was fully established and apparently superseding the RNA world, and at a very
early stage of life beginnings, we already see what has been a recurrent theme
in the likely interpretation of modern comparative molecular studies on
species: the full ability of this nascent life entity to develop itself by
tinkering, using all kinds of available pieces to improve itself. Lastly, and
very instructively, it is deduced from this model that the first peptides to be
produced probably had unique properties, not shared by the arsenal of molecules
present hitherto, which allowed a functional connection between the then
omnipotent RNA world and the lipid membrane vesicles containing it. These
specific functions might have initially been at the origin of the Darwinian
selection of the full-blown protein-synthesis machinery.


  More and more diseases have been found to be strongly correlated with
disturbances in the microbiome constitution, e.g., obesity, diabetes, or some
cancer types. Thanks to modern high-throughput omics technologies, it becomes
possible to directly analyze human microbiome and its influence on the health
status. Microbial communities are monitored over long periods of time and the
associations between their members are explored. These relationships can be
described by a time-evolving graph. In order to understand responses of the
microbial community members to a distinct range of perturbations such as
antibiotics exposure or diseases and general dynamical properties, the
time-evolving graph of the human microbial communities has to be analyzed. This
becomes especially challenging due to dozens of complex interactions among
microbes and metastable dynamics. The key to solving this problem is the
representation of the time-evolving graphs as fixed-length feature vectors
preserving the original dynamics. We propose a method for learning the
embedding of the time-evolving graph that is based on the spectral analysis of
transfer operators and graph kernels. We demonstrate that our method can
capture temporary changes in the time-evolving graph on both created synthetic
data and real-world data. Our experiments demonstrate the efficacy of the
method. Furthermore, we show that our method can be applied to human microbiome
data to study dynamic processes.


  Objective: We study the influence of local reopening policies on the
composition of the infectious population and their impact on future
hospitalization and mortality rates. Materials and Methods: We collected
datasets of daily reported hospitalization and cumulative morality of COVID 19
in Houston, Texas, from May 1, 2020 until June 29, 2020. These datasets are
from multiple sources (USA FACTS, Southeast Texas Regional Advisory Council
COVID 19 report, TMC daily news, and New York Times county level mortality
reporting). Our model, risk stratified SIR HCD uses separate variables to model
the dynamics of local contact (e.g., work from home) and high contact (e.g.,
work on site) subpopulations while sharing parameters to control their
respective $R_0(t)$ over time. Results: We evaluated our models forecasting
performance in Harris County, TX (the most populated county in the Greater
Houston area) during the Phase I and Phase II reopening. Not only did our model
outperform other competing models, it also supports counterfactual analysis to
simulate the impact of future policies in a local setting, which is unique
among existing approaches. Discussion: Local mortality and hospitalization are
significantly impacted by quarantine and reopening policies. No existing model
has directly accounted for the effect of these policies on local trends in
infections, hospitalizations, and deaths in an explicit and explainable manner.
Our work is an attempt to close this important technical gap to support
decision making. Conclusion: Despite several limitations, we think it is a
timely effort to rethink about how to best model the dynamics of pandemics
under the influence of reopening policies.


  The generalized logistic equation is used to interpret the COVID-19 epidemic
data in several countries: Austria, Switzerland, the Netherlands, Italy, Turkey
and South Korea. The model coefficients are calculated: the growth rate and the
expected number of infected people, as well as the exponent indexes in the
generalized logistic equation. It is shown that the dependence of the number of
the infected people on time is well described on average by the logistic curve
(within the framework of a simple or generalized logistic equation) with a
determination coefficient exceeding 0.8. At the same time, the dependence of
the number of the infected people per day on time has a very uneven character
and can be described very roughly by the logistic curve. To describe it, it is
necessary to take into account the dependence of the model coefficients on time
or on the total number of cases. Variations, for example, of the growth rate
can reach 60%. The variability spectra of the coefficients have characteristic
peaks at periods of several days, which corresponds to the observed serial
intervals. The use of the stochastic logistic equation is proposed to estimate
the number of probable peaks in the coronavirus incidence.


  Bacterial genomes have been shown to be partitioned into several kilobases
long chromosomal domains that are topologically independent from each other,
meaning that change of DNA superhelicity in one domain does not propagate to
neighbors. Both in vivo and in vitro experiments have been performed to
question the nature of the topological barriers at play, leading to several
predictions on possible molecular actors. Here, we address the question of
topological barriers using polymer models of supercoiled DNA chains. More
specifically, we determine under which conditions DNA-bridging proteins may act
as topological barriers. To this end, we developed a coarse-grained
bead-and-spring model and investigated its properties through Brownian dynamics
simulations. As a result, we find that DNA-bridging proteins must exert rather
strong constraints on their binding sites: they must block the diffusion of the
excess of twist through the two binding sites on the DNA molecule and,
simultaneously, prevent the rotation of one DNA segment relative to the other
one. Importantly, not all DNA-bridging proteins satisfy this second condition.
For example, single bridges formed by proteins that bind DNA non-specifically,
like H-NS dimers, are expected to fail with this respect. Our findings might
also explain, in the case of specific DNA-bridging proteins like LacI, why
multiple bridges are required to create stable independent topological domains.
Strikingly, when the relative rotation of the DNA segments is not prevented,
relaxation results in complex intrication of the two domains. Moreover, while
the value of the torsional stress in each domain may vary, their differential
is preserved. Our work also predicts that nucleoid associated proteins known to
wrap DNA must form higher protein-DNA complexes to efficiently work as
topological barriers.


  During pandemic events, strategies such as social distancing can be
fundamental to curb viral spreading. Such actions can reduce the number of
simultaneous infections and mitigate the disease spreading, which is relevant
to the risk of a healthcare system collapse. Although these strategies can be
suggested, their actual implementation may depend on the population perception
of the disease risk. The current COVID-19 crisis, for instance, is showing that
some individuals are much more prone than others to remain isolated, avoiding
unnecessary contacts. With this motivation, we propose an epidemiological SIR
model that uses evolutionary game theory to take into account dynamic
individual quarantine strategies, intending to combine in a single process
social strategies, individual risk perception, and viral spreading. The disease
spreads in a population whose agents can choose between self-isolation and a
lifestyle careless of any epidemic risk. The strategy adoption is individual
and depends on the perceived disease risk compared to the quarantine cost. The
game payoff governs the strategy adoption, while the epidemic process governs
the agent's health state. At the same time, the infection rate depends on the
agent's strategy while the perceived disease risk depends on the fraction of
infected agents. Results show recurrent infection waves, which were seen in
previous epidemic scenarios with quarantine. Notably, the risk perception is
found to be fundamental for controlling the magnitude of the infection peak,
while the final infection size is mainly dictated by the infection rates. Low
awareness leads to a single and strong infection peak, while a greater disease
risk leads to shorter, although more frequent, peaks. The proposed model
spontaneously captures relevant aspects of a pandemic event, highlighting the
fundamental role of social strategies.


  We apply topological data analysis, specifically the Mapper algorithm, to the
U.S. COVID-19 data. The resulting Mapper graphs provide visualizations of the
pandemic that are more complete than those supplied by other, more standard
methods. They encode a variety of geometric features of the data cloud created
from geographic information, time progression, and the number of COVID-19
cases. They reflect the development of the pandemic across all of the U.S. and
capture the growth rates as well as the regional prominence of hot-spots. The
Mapper graphs allow for easy comparisons across time and space and have the
potential of becoming a useful predictive tool for the spread of the
coronavirus.


  In this paper, we introduce a novel modeling framework for incorporating fear
of infection and frustration with social distancing into disease dynamics. We
show that the resulting SEIR behavior-perception model has three principal
modes of qualitative behavior---no outbreak, controlled outbreak, and
uncontrolled outbreak. We also demonstrate that the model can produce transient
and sustained waves of infection consistent with secondary outbreaks. We fit
the model to cumulative COVID-19 case and mortality data from several regions.
Our analysis suggests that regions which experience a significant decline after
the first wave of infection, such as Canada and Israel, are more likely to
contain secondary waves of infection, whereas regions which only achieve
moderate success in mitigating the disease's spread initially, such as the
United States, are likely to experience substantial secondary waves or
uncontrolled outbreaks.


  The coronavirus SARS-CoV-2 remains an extant threat against public health on
a global scale. Cell infection begins when the spike protein of SARS-CoV-2
binds with the cell receptor, angiotensin-converting enzyme 2 (ACE2). Here, we
address the role of Tetracycline as an inhibitor for the receptor-binding
domain (RBD) of the spike protein. Targeted molecular investigation show that
Tetracycline binds more favorably to the RBD (-9.40 kcal/mol) compared to
Chloroquine (-6.31 kcal/mol) or Doxycycline (-8.08 kcal/mol) and inhibits
attachment to ACE2 to a greater degree (binding efficiency of 2.98
$\frac{\text{kcal}}{\text{mol}\cdot \text{nm}^2}$ for Tetracycline-RBD, 5.59
$\frac{\text{kcal}}{\text{mol}\cdot \text{nm}^2}$ for Chloroquine-RBD, 5.16
$\frac{\text{kcal}}{\text{mol}\cdot \text{nm}^2}$ for Doxycycline-RBD).
Stronger Tetracycline inhibition is verified with nonequilibrium PMF
calculations, for which the Tetracycline-RBD complex exhibits the lowest free
energy profile along the dissociation pathway from ACE2. Tetracycline appears
to target viral residues that are usually involved in significant hydrogen
bonding with ACE2; this inhibition of cellular infection complements the
anti-inflammatory and cytokine suppressing capability of Tetracycline, and may
further reduce the duration of ICU stays and mechanical ventilation induced by
the coronavirus SARS-CoV-2.


  Governments around the world have implemented preventive measures against the
spread of the coronavirus disease (COVID-19). In this study, we consider a
multivariate discrete-time Markov model to analyze the propagation of COVID-19
across 33 provincial regions in China. This approach enables us to evaluate the
effect of mobility restriction policies on the spread of the disease. We use
data on daily human mobility across regions and apply the Bayesian framework to
estimate the proposed model. The results show that the spread of the disease in
China was predominately driven by community transmission within regions and the
lockdown policy introduced by local governments curbed the spread of the
pandemic. Further, we document that Hubei was only the epicenter of the early
epidemic stage. Secondary epicenters, such as Beijing and Guangdong, had
already become established by late January 2020, and the disease spread out to
connected regions. The transmission from these epicenters substantially
declined following the introduction of human mobility restrictions across
regions.


  COVID-19 is highly transmissible and containing outbreaks requires a rapid
and effective response. Because infection may be spread by people who are
pre-symptomatic or asymptomatic, substantial undetected transmission is likely
to occur before clinical cases are diagnosed. Thus, when outbreaks occur there
is a need to anticipate which populations and locations are at heightened risk
of exposure. In this work, we evaluate the utility of aggregate human mobility
data for estimating the geographic distribution of transmission risk. We
present a simple procedure for producing spatial transmission risk assessments
from near-real-time population mobility data. We validate our estimates against
three well-documented COVID-19 outbreak scenarios in Australia. Two of these
were well-defined transmission clusters and one was a community transmission
scenario. Our results indicate that mobility data can be a good predictor of
geographic patterns of exposure risk from transmission centres, particularly in
scenarios involving workplaces or other environments associated with habitual
travel patterns. For community transmission scenarios, our results demonstrate
that mobility data adds the most value to risk predictions when case counts are
low and spatially clustered. Our method could assist health systems in the
allocation of testing resources, and potentially guide the implementation of
geographically-targeted restrictions on movement and social interaction.


  Epidemiological models are the mathematical models that capture the dynamics
of epidemics. The spread of the virus has two routes - exogenous and
endogenous. The exogenous spread is from outside the population under study,
and endogenous spread is within the population under study. Although some of
the models consider the exogenous source of infection, they have not studied
the interplay between exogenous and endogenous spreads. In this paper, we
introduce a novel model - the Exo-SIR model that captures both the exogenous
and endogenous spread of the virus. We analyze to find out the relationship
between endogenous and exogenous infections during the Covid19 pandemic. First,
we simulate the Exo-SIR model without assuming any contact network for the
population. Second, simulate it by assuming that the contact network is a scale
free network. Third, we implemented the Exo-SIR model on a real dataset
regarding Covid19. We found that endogenous infection is influenced by even a
minimal rate of exogenous infection. Also, we found that in the presence of
exogenous infection, the endogenous infection peak becomes higher, and the peak
occurs earlier. This means that if we consider our response to a pandemic like
Covid19, we should be prepared for an earlier and higher number of cases than
the SIR model suggests if there are the exogenous source(s) of infection.


  We propose a mathematical model to analyze the time evolution of the total
number of infected population with Covid-19 disease at a region in the ongoing
pandemic. Using the available data of Covid-19 infected population on various
countries we formulate a model which can successfully track the time evolution
from early days to the saturation period in a given wave of this infectious
disease. It involves a set of effective parameters which can be extracted from
the available data. Using those parameters the future trajectories of the
disease spread can also be projected. A set of differential equations is also
proposed whose solutions are these time evolution trajectories. Using such a
formalism we project the future time evolution trajectories of infection spread
for a number of countries where the Covid-19 infection is still rapidly rising.


  In this paper, we introduce a new control-theoretic paradigm for mitigating
the spread of a virus. To this end, our discrete-time controller, aims to
reduce the number of new daily deaths, and consequently, the cumulative number
of deaths. In contrast to much of the existing literature, we do not rely on a
potentially complex virus transmission model whose equations must be customized
to the "particulars" of the pandemic at hand. For new viruses such as Covid-19,
the epidemiology driving the modelling process may not be well known and model
estimation with limited data may be unreliable. With this motivation in mind,
the new paradigm described here is data-driven and, to a large extent, we avoid
modelling difficulties by concentrating on just two key quantities which are
common to pandemics: the doubling time, denoted by $d(k)$ and the peak day
denoted by $\theta(k)$. Our numerical studies to date suggest that our
appealingly simple model can provide a reasonable fit to real data. Given that
time is of the essence during the ongoing global health crisis, the intent of
this paper is to introduce this new paradigm to control practitioners and
describe a number of new research directions suggested by our current results.


  As the question of the origins of eukaryotes comes closer to an answer, the
interactions between the archaeal ancestor to eukaryotes and environmental,
molecular, bacterial, other influences and symbionts have become increasingly
relevant. Various studies have pointed towards the Lokiarchaeota as a close
archaeal relative to the eukarya. A study by Imachi et al., published in early
2020, was able to decipher the proteome and genome of the Loki and reported an
entirely new model for eukaryogenesis due to its findings, thus shifting
research of eukaryogenesis away from models focusing on endocytosis and the
rapid acquisition of the alphaproteobacteria that would become the
mitochondria. Imachi and colleagues proposed that ancestors of the Loki lived
symbiotically with other prokaryotes due to its shape, and facilitated
molecular transfer between the symbionts as well as the membrane manipulation
necessary to allow such transfer, before acquiring a proto-mitochondria over
time. This paper will lay out the history of endosymbiosis and the prevailing
theory of endocytosis between archaea and alphaproteobacteria, and consolidates
the knowns and unknowns up to before the publication of Imachi et al. The
findings of their novel study will be discussed, and the knowns and unknowns
with respect to these new findings are re-consolidated. The two models of
eukaryogenesis are compared and, in moving forward with Imachi and colleagues
model, probable areas of research to further validate this approach are
discussed, as well as potential usage of research tools of other fields to
answer the persistent unknowns of how eukaryotes came to be.


  Repeated asymptomatic screening for SARS-CoV-2 promises to control spread of
the virus but would require too many resources to implement at scale. Group
testing is promising for screening more people with fewer test resources:
multiple samples tested together in one pool can be excluded with one negative
test result. Existing approaches to group testing design for SARS-CoV-2
asymptomatic screening, however, do not consider dilution effects: that false
negatives become more common with larger pools. As a consequence, they may
recommend pool sizes that are too large or misestimate the benefits of
screening. Modeling dilution effects, we derive closed-form expressions for the
expected number of tests and false negative/positives per person screened under
two popular group testing methods: the linear and square array methods. We find
that test error correlation induced by a common viral load across an
individual's samples results in many fewer false negatives than would be
expected from less realistic but more widely assumed independent errors. This
insight also suggests that false positives can be controlled through repeated
tests without significantly increasing false negatives. Using these closed-form
expressions to trace a Pareto frontier over error rates and tests, we design
testing protocols for repeated asymptomatic screening of a large population. We
minimize disease prevalence by optimizing a time-varying pool sizes and
screening frequency constrained by daily test capacity and a false positive
limit. This provides a testing protocol practitioners can use for mitigating
COVID-19. In a case study, we demonstrate the effectiveness of this methodology
in controlling spread.


  In this paper we analyze the effects of commuting and social inequalities for
the epidemic development of the novel coronavirus (COVID-19). With this aim we
consider a SEIRD (susceptible, exposed, infected, recovered and dead by
disease) model without vital dynamics in a population divided into patches that
have different economic resources and in which the individuals can commute from
one patch to another (bilaterally). In the modeling we choose the social and
commuting parameters arbitrarily. We calculate the basic reproductive number
$R_0$ with the next generation approach and analyze the sensitivity of $R_0$
with respect to the parameters. Furthermore, we run numerical simulations
considering a population divided into two patches to bring some conclusions on
the number of total infected individuals and cumulative deaths for our model
considering heterogeneous populations.


  In this paper we use asymptotic analysis to determine the steady-state mean
number of resources in each of $N$ small interior targets within a
three-dimensional bounded domain. The accumulation of resources is based on
multiple rounds of search-and-capture events; whenever a searcher finds a
target it delivers a resource packet to the target, after which it escapes and
returns to its initial position (resetting after capture). The searcher is then
resupplied with cargo and a new search process is initiated after a random
delay. Assuming that the accumulation of resources is counterbalanced by
degradation, one can derive general expressions for the moments of the resource
distribution. We use this to show that the mean number of resources in a target
is proportional to its effective "shape capacitance." We then extend the
analysis to the case of diffusive search with stochastic resetting before
capture, where the position of the searcher is reset to its initial position at
a random sequence of times that is statistically independent of the ongoing
search process, in contrast to the sequence of resetting times after capture.


  Reproduction numbers, like the basic reproduction number $\mathcal{R}_0$,
play an important role in the analysis and application of dynamic models,
including contagion models and ecological population models. One difficulty in
deriving these quantities is that they must be computed on a model-by-model
basis, since it is typically impractical to obtain general reproduction number
expressions applicable to a family of related models, especially if these are
of different dimensions. For example, this is typically the case for SIR-type
infectious disease models derived using the linear chain trick (LCT). Here we
show how to find general reproduction number expressions for such models
families (which vary in their number of state variables) using the next
generation operator approach in conjunction with the generalized linear chain
trick (GLCT). We further show how the GLCT enables modelers to draw insights
from these results by leveraging theory and intuition from continuous time
Markov chains (CTMCs) and their absorption time distributions (i.e., phase-type
probability distributions). To do this, we first review the GLCT and other
connections between mean-field ODE model assumptions, CTMCs, and phase-type
distributions. We then apply this technique to find reproduction numbers for
two sets of models: a family of generalized SEIRS models of arbitrary finite
dimension, and a generalized family of finite dimensional predator-prey
(Rosenzweig-MacArthur type) models. These results highlight the utility of the
GLCT for the derivation and analysis of mean field ODE models, especially when
used in conjunction with theory from CTMCs and their associated phase-type
distributions.


  To describe the dynamics of social distancing during pandemics, we follow
previous efforts to combine basic epidemiology models (e.g. SIR - Susceptible,
Infected, and Recovered) with game and economy theory tools. We present an
extension of the SIR model that predicts a series of discontinuous transitions
in social distancing. Each transition resembles a phase transition of the
second-order (Ginzburg-Landau instability) and, therefore, potentially a
general phenomenon. The first wave of COVID-19 led to social distancing around
the globe: severe lockdowns to stop the pandemic were followed by a series of
lockdown lifts. Data analysis of the first wave in Austria, Israel, and Germany
corroborates the soundness of the model. Furthermore, this work presents
analytical tools to analyze pandemic waves, which may be extended to calculate
derivatives of giant components in network percolation transitions and may also
be of interest in the context of crisis formation theories.


  treePL uses a penalised likelihood approach to produce a dated phylogeny in a
maximum likelihood framework. Since its publication in 2012, few resources have
been developed to explain how to use it properly. In this guide, I provide a
step-by-step protocol for producing a dated phylogeny using treePL, based on my
experience building a large dated phylogeny with it and conducting additional
tests on a smaller phylogeny. I also provide the necessary data to reproduce
one of the example phylogenies presented. I compare these treePL phylogenies to
BEAST2-built counterparts. Even though I cannot explain precisely how treePL
works, the evidence discussed in this guide suggest that the empirical protocol
presented is reliable.