  To contain the propagation of emerging diseases that are transmissible from
human to human, non-pharmaceutical interventions (NPIs) aimed at reducing the
interactions between humans are usually implemented. One example of the latter
kind of measures is social distancing, which can be either policy-driven or can
arise endogenously in the population as a consequence of the fear of infection.
However, if NPIs are lifted before the population reaches herd immunity,
further re-introductions of the pathogen would lead to secondary infections.
Here we study the effects of different social distancing schemes on the large
scale spreading of diseases. Specifically, we generalize metapopulation models
to include social distancing mechanisms at the subpopulation level and model
short- and long-term strategies that are fed with local or global information
about the epidemics. We show that different model ingredients might lead to
very diverse outcomes in different subpopulations. Our results suggest that
there is not a unique answer to the question of whether contention measures are
more efficient if implemented and managed locally or globally and that model
outcomes depends on how the full complexity of human interactions is taken into
account.


  Environmental heterogeneity is a determining factor of the structure of
biological communities. Thus, understanding the distribution of species along
environmental gradients provides assistance to conservation. The goal of this
study was to determine the distribution pattern of the herbaceous community in
three areas of the Southern Amazon. Sampling was conducted in three modules
totaling 39 permanent plots according to the protocol of collection of the
Program for Research in Biodiversity. All herbaceous and ground hemiepiphyte
subjects above 5 cm were recorded. Multivariate analyses were used to summarize
the species composition, multiple regression models were used to determine if
environmental variables and disturbance caused by logging influenced the
composition of the herbaceous community. We recorded 7.965 individuals
representing 70 species. The distance of the watercourse was the main factor
associated with the distribution of the species, interactions between variables
showed that canopy openness and sand content also influence the species
composition, and there was no effect on the number of trees cut. Species
richness increased in areas where canopy cover was higher and it decreases as
it becomes more distant from the watercourse. The occurrences of preferred
habitats for some species have, in addition to an ecological interest, a
practical significance for the conservation and management of these species.
Currently, the area of preservation of streams provided by the Forest Code in
effect is 30 m for rivers up to 10 m wide. However, this study shows that the
range of protection should be extended to at least 100 m wide.


  Replication of SARS-CoV2 virions is an error-prone process which may
eventually generate a percentage of impaired protein copies with complete lack
of functionality. For instance, after RNA mis-replication, a very premature
stop codon (UAG, UAA, UGA) coding for a prematurely truncated
(nonsense-mutated) spike protein may occur.
  In the natural virus replication process via cell infection, the nonsense
genomes are corrected by the proofreading enzymes of the virus, strongly
penalized by natural selection and condemned to a very short life by the host
cell's mRNA watching mechanisms. However, for the very long spike genome of
1273 codons, a truncated non-functional spike protein may potentially still
occur with a high frequency, even in presence of a low mutation rate per single
nucleotide.
  With this paper, a hi-fidelity post-processing of SARS-CoV2 spike sequences
is provided: in ex-vivo samples from patients, an impressively high rate of
26\% of prematurely-stopped (nonsense-mutated) spike genomes sequences due to
insertions/deletions is found, compared with a 9.7\% obtained from in-vitro
cell culture.
  A general warning on the possible high rate of prematurely-stopped spike
protein sequences is also raised for "artificial" de novo DNA synthesis
processes of SARS-CoV2 spike genomes with no associated natural
proofreading/selection, possibly including vaccine preparations.
  Finally, a metric based on the ratio between prematurely stopped and "normal"
genomes is proposed as a potential host-independent variant-watching tool, able
to classify the infectivity of new spike mutations.


  We propose and study a new mathematical model of the human immunodeficiency
virus (HIV). The main novelty is to consider that the antibody growth depends
not only on the virus and on the antibodies concentration but also on the
uninfected cells concentration. The model consists of five nonlinear
differential equations describing the evolution of the uninfected cells, the
infected ones, the free viruses, and the adaptive immunity. The adaptive immune
response is represented by the cytotoxic T-lymphocytes (CTL) cells and the
antibodies with the growth function supposed to be trilinear. The model
includes two kinds of treatments. The objective of the first one is to reduce
the number of infected cells, while the aim of the second is to block free
viruses. Firstly, the positivity and the boundedness of solutions are
established. After that, the local stability of the disease free steady state
and the infection steady states are characterized. Next, an optimal control
problem is posed and investigated. Finally, numerical simulations are performed
in order to show the behavior of solutions and the effectiveness of the two
incorporated treatments via an efficient optimal control strategy.


  Biology offers many examples of large-scale, complex, concurrent systems:
many processes take place in parallel, compete on resources and influence each
other's behavior. The scalable modeling of biological systems continues to be a
very active field of research. In this paper we introduce a new approach based
on Event-B, a state-based formal method with refinement as its central
ingredient, allowing us to check for model consistency step-by-step in an
automated way. Our approach based on functions leads to an elegant and concise
modeling method. We demonstrate this approach by constructing what is, to our
knowledge, the largest ever built Event-B model, describing the ErbB signaling
pathway, a key evolutionary pathway with a significant role in development and
in many types of cancer. The Event-B model for the ErbB pathway describes 1320
molecular reactions through 242 events.


  Airborne diseases can be transmitted by infectious aerosols in the near
field, i.e., in close proximity, or in the far field, i.e., by infectious
aerosols that are well mixed within the indoor air. Is it possible to say which
mode of disease transmission is predominant in large indoor spaces? We
addressed this question by measuring the transport of aerosols equivalent to
the size of human respiratory particles in two large hardware stores (V>10000
m$^3$). We found that aerosol concentrations in both stores decreased rapidly
and almost independently of aerosol size, despite the different ventilation
systems. A persistent and directional airflow on the order of a few cm/s was
observed in both stores. Consequently, aerosol dynamics in such open settings
can be expected to be dominated by turbulent dispersion and sweeping, and the
accumulation of infectious aerosols in the indoor air is unlikely to contribute
significantly to the risk of infection as long as the occupancy of the store is
not too high. Under these conditions, well-fitting face masks are an excellent
means of preventing disease transmission by human aerosols.


  Antibiotic resistance is a growing public health problem. To gain a
fundamental understanding of resistance evolution, a combination of systematic
experimental and theoretical approaches is required. Evolution experiments
combined with next-generation sequencing techniques, laboratory automation, and
mathematical modeling are enabling the investigation of resistance development
at an unprecedented level of detail. Recent work has directly tracked the
intricate stochastic dynamics of bacterial populations in which resistant
mutants emerge and compete. In addition, new approaches have enabled measuring
how prone a large number of genetically perturbed strains are to evolve
resistance. Based on advances in quantitative cell physiology, predictive
theoretical models of resistance are increasingly being developed. Taken
together, a new strategy for observing, predicting, and ultimately controlling
resistance evolution is emerging.


  Knowledge Graphs (KG) and associated Knowledge Graph Embedding (KGE) models
have recently begun to be explored in the context of drug discovery and have
the potential to assist in key challenges such as target identification. In the
drug discovery domain, KGs can be employed as part of a process which can
result in lab-based experiments being performed, or impact on other decisions,
incurring significant time and financial costs and most importantly, ultimately
influencing patient healthcare. For KGE models to have impact in this domain, a
better understanding of not only of performance, but also the various factors
which determine it, is required. In this study we investigate, over the course
of many thousands of experiments, the predictive performance of five KGE models
on two public drug discovery-oriented KGs. Our goal is not to focus on the best
overall model or configuration, instead we take a deeper look at how
performance can be affected by changes in the training setup, choice of
hyperparameters, model parameter initialisation seed and different splits of
the datasets. Our results highlight that these factors have significant impact
on performance and can even affect the ranking of models. Indeed these factors
should be reported along with model architectures to ensure complete
reproducibility and fair comparisons of future work, and we argue this is
critical for the acceptance of use, and impact of KGEs in a biomedical setting.


  We examine a pair of graph generative models for the therapeutic design of
novel drug candidates targeting SARS-CoV-2 viral proteins. Due to a sense of
urgency, we chose well-validated models with unique strengths: an autoencoder
that generates molecules with similar structures to a dataset of drugs with
anti-SARS activity and a reinforcement learning algorithm that generates highly
novel molecules. During generation, we explore optimization toward several
design targets to balance druglikeness, synthetic accessability, and anti-SARS
activity based on \icfifty. This generative
framework\footnote{https://github.com/exalearn/covid-drug-design} will
accelerate drug discovery in future pandemics through the high-throughput
generation of targeted therapeutic candidates.


  Short-read DNA sequencing instruments can yield over 1e+12 bases per run,
typically composed of reads 150 bases long. Despite this high throughput, de
novo assembly algorithms have difficulty reconstructing contiguous genome
sequences using short reads due to both repetitive and difficult-to-sequence
regions in these genomes. Some of the short read assembly challenges are
mitigated by scaffolding assembled sequences using paired-end reads. However,
unresolved sequences in these scaffolds appear as "gaps". Here, we introduce
GapPredict, a tool that uses a character-level language model to predict
unresolved nucleotides in scaffold gaps. We benchmarked GapPredict against the
state-of-the-art gap-filling tool Sealer, and observed that the former can fill
65.6% of the sampled gaps that were left unfilled by the latter, demonstrating
the practical utility of deep learning approaches to the gap-filling problem in
genome sequence assembly.


  Changes in human behavior are increasingly recognized as a major determinant
of epidemic dynamics. Although collective activity can be modified through
imposed measures to control epidemic progression, spontaneous changes can also
arise as a result of uncoordinated individual responses to the perceived risk
of contagion. Here we introduce a stochastic epidemic model that implements
population responses driven by individual- and time-dependent risk-taking
propensity. The model reveals an emergent mechanism for the generation of
multiple infection waves of decreasing amplitude without the need to consider
external modulation of parameters. Successive waves tune the effective
reproduction number to its critical value $R=1$. This process is a consequence
of the interplay of the fractions of susceptible and infected population and
the average risk-taking propensity, as shown by a mean-field approach. The
proposed mechanism also shows how, under the threat of contagion, the
distribution of individual risk propensities evolves towards a well-defined
profile. Successive waves trigger selective sweeps of risk-taking propensity at
a pace determined by individual risk reaction rates. This kind of collective,
self-generated pressure, may therefore shape risk-aversion profiles associated
to epidemics in human groups. The final state is self-organized and generic,
independent of the parameter values. We conclude that uncoordinated changes in
human behavior can, by themselves, explain major qualitative and quantitative
features of the epidemic process, as the emergence of multiple waves and the
tendency to remain around $R=1$ observed worldwide after the first few waves of
COVID-19.


  Standard approaches to decision-making under uncertainty focus on sequential
exploration of the space of decisions. However, \textit{simultaneously}
proposing a batch of decisions, which leverages available resources for
parallel experimentation, has the potential to rapidly accelerate exploration.
We present a family of (parallel) contextual bandit algorithms applicable to
problems with bounded eluder dimension whose regret is nearly identical to
their perfectly sequential counterparts -- given access to the same total
number of oracle queries -- up to a lower-order ``burn-in" term. We further
show these algorithms can be specialized to the class of linear reward
functions where we introduce and analyze several new linear bandit algorithms
which explicitly introduce diversity into their action selection. Finally, we
also present an empirical evaluation of these parallel algorithms in several
domains, including materials discovery and biological sequence design problems,
to demonstrate the utility of parallelized bandits in practical settings.


  In this work, we derive a nonstandard finite difference scheme for the SICA
(Susceptible-Infected-Chronic-AIDS) model and analyze the dynamical properties
of the discretized system. We prove that the discretized model is dynamically
consistent with the continuous, maintaining the essential properties of the
standard SICA model, namely, the positivity and boundedness of the solutions,
equilibrium points, and their local and global stability.


  We study real steady state varieties of the dynamics of chemical reaction
networks. The dynamics are derived using mass action kinetics with parametric
reaction rates. The models studied are not inherently parametric in nature.
Rather, our interest in parameters is motivated by parameter uncertainty, as
reaction rates are typically either measured with limited precision or
estimated. We aim at detecting toricity and shifted toricity, using a framework
that has been recently introduced and studied for the non-parametric case over
both the real and the complex numbers. While toricity requires that the variety
specifies a subgroup of the direct power of the multiplicative group of the
underlying field, shifted toricity requires only a coset. In the non-parametric
case these requirements establish real decision problems. In the presence of
parameters we must go further and derive necessary and sufficient conditions in
the parameters for toricity or shifted toricity to hold. Technically, we use
real quantifier elimination methods. Our computations on biological networks
here once more confirm shifted toricity as a relevant concept, while toricity
holds only for degenerate parameter choices.


  Human mobility is an important factor in the propagation of infectious
diseases. In particular, the spatial spread of a disease is a consequence of
human mobility. On the other hand, the control strategies based on mobility
restrictions are generally unpopular and costly. These high social and economic
costs make it very important to design global protocols where the cost is
minimized and effects maximized. In this work, we calculate the percolation
threshold of the spread in a network of a disease. In particular, we found the
number of roads to close and regions to isolate in the Puebla State, Mexico, to
avoid the global spread of COVID-19. Computational simulations taking into
account the proposed strategy show a potential reduction of 94% of infections.
This methodology can be used in broader and different areas to help in the
design of health policies.
  --
  La movilidad de las personas es uno de los principales factores que propician
la propagaci\'on de epidemias. En particular, es el factor que genera el
esparcimiento de la enfermedad en diferentes regiones. Las medidas de control
epidemiol\'ogico basadas en la restricci\'on de movilidad son generalmente poco
populares y las consecuencias econ\'omicas pueden llegar a ser muy grandes.
Debido a los altos costos de estas medidas, es de gran relevancia tener
estrategias globales que optimicen las medidas minimizando los costos. En este
trabajo, se calcula el umbral de percolaci\'on de la propagaci\'on de
enfermedades en redes. De manera particular, se encuentra el n\'umero de
caminos a restringir y localidades que tienen que ser aisladas para limitar la
propagaci\'on global de COVID-19 en el Estado de Puebla, M\'exico. Simulaciones
computacionales donde se implementan las medidas de restricci\'on de movilidad
entre los diferentes municipios, junto con las medidas de confinamiento,
muestran que es posible reducir un 94% de la poblaci\'on afectada comparado con
el caso en el que no se implementa ninguna medida. Esta metodolog\'ia puede ser
aplicada a distintas zonas para ayudar a las autoridades de salud en la toma de
decisiones.


  The isolation of infectious individuals is a key measure of public health for
the control of communicable diseases. However, involving a strong perturbation
of daily life, it often causes psychosocial distress, and severe financial and
social costs. These may act as mechanisms limiting the adoption of the measure
in the first place or the adherence throughout its full duration. In addition,
difficulty of recognizing mild symptoms or lack of symptoms may impact
awareness of the infection and further limit adoption. Here, we study an
epidemic model on a network of contacts accounting for limited adherence and
delayed awareness to self-isolation, along with fatigue causing overhasty
termination. The model allows us to estimate the role of each ingredient and
analyze the tradeoff between adherence and duration of self-isolation. We find
that the epidemic threshold is very sensitive to an effective compliance that
combines the effects of imperfect adherence, delayed awareness and fatigue. If
adherence improves for shorter quarantine periods, there exists an optimal
duration of isolation, shorter than the infectious period. However,
heterogeneities in the connectivity pattern, coupled to a reduced compliance
for highly active individuals, may almost completely offset the effectiveness
of self-isolation measures on the control of the epidemic.


  While recent work has established divergence as a key framework for
understanding evenness, there is currently no research exploring how the
families of measures within the divergence-based framework relate to each
other. This paper uses geometry to show that, holding order and richness
constant, the families of divergence-based evenness measures nest. This
property allows them to be ranked based on their reactivity to changes in
relatively even assemblages or changes in relatively uneven ones. We establish
this ranking and explore how the distance-based measures relate to it for both
order q=2 and q=1. We also derive a new family of distance-based measures that
captures the angular distance between the vector of relative abundances and a
perfectly even vector and is highly reactive to changes in even assemblages.
Finally, we show that if we only require evenness to be a divergence, then any
smooth, monotonically increasing function of diversity can be made into an
evenness measure. A deeper understanding of how to measure evenness will
require empirical or theoretical research that uncovers which kind of
divergence best reflects the underlying concept.


  Although traditional models of epidemic spreading focus on the number of
infected, susceptible and recovered individuals, a lot of attention has been
devoted to integrate epidemic models with population genetics. Here we develop
an individual-based model for epidemic spreading on networks in which viruses
are explicitly represented by finite chains of nucleotides that can mutate
inside the host. Under the hypothesis of neutral evolution we compute
analytically the average pairwise genetic distance between all infecting
viruses over time. We also derive a mean-field version of this equation that
can be added directly to compartmental models such as SIR or SEIR to estimate
the genetic evolution. We compare our results with the inferred genetic
evolution of SARS-CoV-2 at the beginning of the epidemic in China and found
good agreement with the analytical solution of our model. Finally, using
genetic distance as a proxy for different strains, we use numerical simulations
to show that the lower the connectivity between communities, e.g., cities, the
higher the probability of reinfection.


  We evaluate the efficiency of various heuristic strategies for allocating
vaccines against COVID-19 and compare them to strategies found using optimal
control theory. Our approach is based on a mathematical model which tracks the
spread of disease among different age groups and across different geographical
regions, and we introduce a method to combine age-specific contact data to
geographical movement data. As a case study, we model the epidemic in the
population of mainland Finland utilizing mobility data from a major telecom
operator. Our approach allows to determine which geographical regions and age
groups should be targeted first in order to minimize the number of deaths. In
the scenarios that we test, we find that distributing vaccines demographically
and in an age-descending order is not optimal for minimizing deaths and the
burden of disease. Instead, more lives could potentially be saved by using
strategies which emphasize high-incidence regions and distribute vaccines in
parallel to multiple age groups. The level of emphasis that high-incidence
regions should be given depends on the overall transmission rate in the
population. This observation highlights the importance of updating the
vaccination strategy when the effective reproduction number changes due to the
general contact patterns changing and new virus variants entering.


  Disease transmission is studied through disciplines like epidemiology,
applied mathematics, and statistics. Mathematical simulation models for
transmission have implications in solving public and personal health
challenges. The SIR model uses a compartmental approach including dynamic and
nonlinear behavior of transmission through three factors: susceptible,
infected, and removed (recovered and deceased) individuals. Using the Lambert W
Function, we propose a framework to study solutions of the SIR model. This
demonstrates the applications of COVID-19 transmission data to model the spread
of a real-world disease. Different models of disease including the SIR, SIRmp
and SEIRpqr model are compared with respect to their ability to predict disease
spread. Physical distancing impacts and personal protection equipment use are
discussed with relevance to the COVID-19 spread.


  At odds with a traditional view of molecular evolution that seeks a
descent-with-modification relationship between functional sequences, new
functions can emerge {\it de novo} with relative ease. At early times of
molecular evolution, random polymers could have sufficed for the appearance of
incipient chemical activity, while the cellular environment harbors a myriad of
proto-functional molecules. The emergence of function is facilitated by several
mechanisms intrinsic to molecular organization, such as redundant mapping of
sequences into structures, phenotypic plasticity, modularity, or cooperative
associations between genomic sequences. It is the availability of niches in the
molecular ecology that filters new potentially functional proposals. New
phenotypes and subsequent levels of molecular complexity could be attained
through combinatorial explorations of currently available molecular variants.
Natural selection does the rest.


  Prediction of node and graph labels are prominent network science tasks. Data
analyzed in these tasks are sometimes related: entities represented by nodes in
a higher-level (higher-scale) network can themselves be modeled as networks at
a lower level. We argue that systems involving such entities should be
integrated with a "network of networks" (NoN) representation. Then, we ask
whether entity label prediction using multi-level NoN data via our proposed
approaches is more accurate than using each of single-level node and graph data
alone, i.e., than traditional node label prediction on the higher-level network
and graph label prediction on the lower-level networks. To obtain data, we
develop the first synthetic NoN generator and construct a real biological NoN.
We evaluate accuracy of considered approaches when predicting artificial labels
from the synthetic NoNs and proteins' functions from the biological NoN. For
the synthetic NoNs, our NoN approaches outperform or are as good as node- and
network-level ones depending on the NoN properties. For the biological NoN, our
NoN approaches outperform the single-level approaches for just under half of
the protein functions, and for 30% of the functions, only our NoN approaches
make meaningful predictions, while node- and network-level ones achieve random
accuracy. So, NoN-based data integration is important.


  Understanding the causes and consequences of, and devising countermeasures
to, global warming is a profoundly complex problem. Network representations are
sometimes the only way forward, and sometimes able to reduce the complexity of
the original problem. Networks are both necessary and natural elements of
climate science. Furthermore, networks form a mathematical foundation for a
multitude of computational and analytical techniques. We are only beginning to
see the benefits of this connection between the sciences of climate change and
network science. In this review, we cover the wide spectrum of network
applications in the climate-change literature -- what they represent, how they
are analyzed, and what insights they bring. We also discuss network data,
tools, and problems yet to be explored.


  Cholera, a severe gastrointestinal infection caused by the bacterium Vibrio
cholerae, remains a major threat to public health with a yearly estimated
global burden of 2.9 million cases. Although the majority of existing models
for the disease focus on its population dynamics, it's important to link the
multiple scales of the disease to gain better perspectives on its spread and
control. In this study, we formulate an immuno-epidemiological model for
cholera linking the between-host and within-host dynamics of the disease. The
within-host model utilizes time-scale methods to differentiate the pathogen
dynamics from the dynamics of the immune response. Bifurcation analysis of the
within-host system reveals the necessary conditions for the existence of both
the Hopf and saddle-node bifurcations. Contrary to other within-host models,
the current approach allows for the elimination of the pathogen after a finite
time. The epidemic model takes into account the direct human-to-human
transmission route of the infection as well as the transmission via the
environment. It is represented by a dynamical system structured on the immune
status, which is a function derived from the within-host immune response. The
basic reproduction number is derived and the stability of equilibrium points
analysed. Analysis of the endemic equilibrium reveals additional constraints
that lead to its stability. Without loss of immunity, the endemic equilibrium,
if it exists, is globally asymptotically stable.


  We review the development of thermodynamic protein hydropathic scaling
theory, starting from backgrounds in mathematics and statistical mechanics, and
leading to biomedical applications. Darwinian evolution has organized each
protein family in different ways, but dynamical hydropathic scaling theory is
both simple and effective in providing readily transferable dynamical insights
for many proteins represented in the uncounted amino acid sequences, as well as
the 90 thousand static structures contained in the online Protein Data Base.
Critical point theory is general, and recently it has proved to be the most
effective way of describing protein networks that have evolved towards nearly
perfect functionality in given environments, self-organized criticality.
Darwinian evolutionary patterns are governed by common dynamical hydropathic
scaling principles, which can be quantified using scales that have been
developed bioinformatically by studying thousands of static PDB structures. The
most effective dynamical scales involve hydropathic globular sculpting
interactions averaged over length scales centered on domain dimensions. A
central feature of dynamical hydropathic scaling theory is the characteristic
domain length associated with a given protein functionality. Evolution has
functioned in such a way that the minimal critical length scale established so
far is about nine amino acids, but in some cases it is much larger. Some
ingenuity is needed to find this primary length scale, as shown by the examples
discussed here. Often a survey of the Darwinian evolution of a protein sequence
suggests a means of determining the critical length scale. The evolution of
Coronavirus is an interesting application; it identifies critical mutations.


  We construct a family of genealogy-valued Markov processes that are induced
by a continuous-time Markov population process. We derive exact expressions for
the likelihood of a given genealogy conditional on the history of the
underlying population process. These lead to a nonlinear filtering equation
which can be used to design efficient Monte Carlo inference algorithms. We
demonstrate these calculations with several examples. Existing full-information
approaches for phylodynamic inference are special cases of the theory.


  The lack of explainability is one of the most prominent disadvantages of deep
learning applications in omics. This "black box" problem can undermine the
credibility and limit the practical implementation of biomedical deep learning
models. Here we present XOmiVAE, a variational autoencoder (VAE) based
interpretable deep learning model for cancer classification using
high-dimensional omics data. XOmiVAE is capable of revealing the contribution
of each gene and latent dimension for each classification prediction, and the
correlation between each gene and each latent dimension. It is also
demonstrated that XOmiVAE can explain not only the supervised classification
but the unsupervised clustering results from the deep learning network. To the
best of our knowledge, XOmiVAE is one of the first activation level-based
interpretable deep learning models explaining novel clusters generated by VAE.
The explainable results generated by XOmiVAE were validated by both the
performance of downstream tasks and the biomedical knowledge. In our
experiments, XOmiVAE explanations of deep learning based cancer classification
and clustering aligned with current domain knowledge including biological
annotation and academic literature, which shows great potential for novel
biomedical knowledge discovery from deep learning models.


  A protein is a linear chain containing a set of amino acids, which folds on
itself to create a specific native structure, also called the minimum energy
conformation. It is the native structure that determines the functionality of
each protein. The protein folding problem (PFP) remains one of the more
difficult problems in computational and chemical biology. The principal
challenge of PFP is to predict the optimal conformation of a given protein by
considering only its amino acid sequence. As the conformational space contains
a very large number of conformations, even when addressing short sequences,
different simplified models have been developed and applied to make the PFP
less complex. In the last few years, many computational approaches have been
proposed to solve the PFP. They are based on simplified lattice models such as
the hydrophobic-polar model. In this paper, we present a new Hybrid Cuckoo
Search Algorithm (HCSA) to solve the 3D-HP protein folding optimization
problem. Our proposed algorithm consists of combining the Cuckoo Search
Algorithm (CSA) with the Hill Climbing (HC) algorithm. Simulation results on
different benchmark sequences are presented and compared to the
state-of-the-art algorithms.


  The estimate of the remaining time of an ongoing wave of epidemic spreading
is a critical issue. Due to the variations of a wide range of parameters in an
epidemic, for simple models such as Susceptible-Infected-Removed (SIR) model,
it is difficult to estimate such a time scale. On the other hand,
multidimensional data with a large set attributes are precisely what one can
use in statistical learning algorithms to make predictions. Here we show, how
the predictability of the SIR model changes with various parameters using a
supervised learning algorithm. We then estimate the condition in which the
model gives the least error in predicting the duration of the first wave of the
COVID-19 pandemic in different states in India. Finally, we use the SIR model
with the above mentioned optimal conditions to generate a training data set and
use it in the supervised learning algorithm to estimate the end-time of the
ongoing second wave of the pandemic in different states in India.


  Physically, disordered ensembles of non-homopolymeric polypeptides are
expected to be heterogeneous; i.e., they should differ from those homogeneous
ensembles of homopolymers that harbor an essentially unique relationship
between average values of end-to-end distance $R_{\rm EE}$ and radius of
gyration $R_{\rm g}$. It was posited recently, however, that small-angle X-ray
scattering (SAXS) data on conformational dimensions of disordered proteins can
be rationalized almost exclusively by homopolymer ensembles. Assessing this
perspective, chain-model simulations are used to evaluate the discriminatory
power of SAXS-determined molecular form factors (MFFs) with regard to
homogeneous versus heterogeneous ensembles. The general approach adopted here
is not bound by any assumption about ensemble encodability, in that the
postulated heterogeneous ensembles we evaluated are not restricted to those
entailed by simple interaction schemes. Our analysis of MFFs for certain
heterogeneous ensembles with more narrowly distributed $R_{\rm EE}$ and $R_{\rm
g}$ indicates that while they deviates from MFFs of homogeneous ensembles, the
differences can be rather small. Remarkably, some heterogeneous ensembles with
asphericity and $R_{\rm EE}$ drastically different from those of homogeneous
ensembles can nonetheless exhibit practically identical MFFs, demonstrating
that SAXS MFFs do not afford unique characterizations of basic properties of
conformational ensembles in general. In other words, the ensemble to MFF
mapping is practically many-to-one and likely non-smooth. Heteropolymeric
variations of the $R_{\rm EE}$--$R_{\rm g}$ relationship were further showcased
using an analytical perturbation theory developed here for flexible
heteropolymers. Ramifications of our findings for interpretation of
experimental data are discussed.


  Background: Prokaryotic viruses, which infect bacteria and archaea, are the
most abundant and diverse biological entities in the biosphere. To understand
their regulatory roles in various ecosystems and to harness the potential of
bacteriophages for use in therapy, more knowledge of viral-host relationships
is required. High-throughput sequencing and its application to the microbiome
have offered new opportunities for computational approaches for predicting
which hosts particular viruses can infect. However, there are two main
challenges for computational host prediction. First, the empirically known
virus-host relationships are very limited. Second, although sequence similarity
between viruses and their prokaryote hosts have been used as a major feature
for host prediction, the alignment is either missing or ambiguous in many
cases. Thus, there is still a need to improve the accuracy of host prediction.
Results: In this work, we present a semi-supervised learning model, named
HostG, to conduct host prediction for novel viruses. We construct a knowledge
graph by utilizing both virus-virus protein similarity and virus-host DNA
sequence similarity. Then graph convolutional network (GCN) is adopted to
exploit viruses with or without known hosts in training to enhance the learning
ability. During the GCN training, we minimize the expected calibrated error
(ECE) to ensure the confidence of the predictions. We tested HostG on both
simulated and real sequencing data and compared its performance with other
state-of-the-art methods specifcally designed for virus host classification
(VHM-net, WIsH, PHP, HoPhage, RaFAH, vHULK, and VPF-Class). Conclusion: HostG
outperforms other popular methods, demonstrating the efficacy of using a
GCN-based semi-supervised learning approach. A particular advantage of HostG is
its ability to predict hosts from new taxa.


  As the only thiol-bearing amino acid, cysteine (Cys) residues in proteins
have the reactive thiol side chain, which is susceptible to a series of
post-translational modifications (PTMs). These PTMs participate in a wide range
of biological activities including the alteration of enzymatic reactions,
protein-protein interactions and protein stability. Here we summarize the
advance of cysteine PTM identification technologies and the features of the
various kinds of the PTMs. We also discuss in silico approaches for the
prediction of the different types of cysteine modified sites, giving directions
for future study.


  Pavelka et al. (Science, Reports, 7 May 2021) claim that a single round of
population-wide antigen testing in Slovakia reduced the observed COVID-19
prevalence by 58%, and that it played a substantial role in curbing the
pandemic. We argue that this estimate, which is based on incorrect assumptions,
is exaggerated, and that the relief was short-lived with little effect on
mitigating the pandemic.


  The identification of protein-ligand interaction plays a key role in
biochemical research and drug discovery. Although deep learning has recently
shown great promise in discovering new drugs, there remains a gap between deep
learning-based and experimental approaches. Here we propose a novel framework,
named AIMEE, integrating AI Model and Enzymology Experiments, to identify
inhibitors against 3CL protease of SARS-CoV-2, which has taken a significant
toll on people across the globe. From a bioactive chemical library, we have
conducted two rounds of experiments and identified six novel inhibitors with a
hit rate of 29.41%, and four of them showed an IC50 value less than 3 {\mu}M.
Moreover, we explored the interpretability of the central model in AIMEE,
mapping the deep learning extracted features to domain knowledge of chemical
properties. Based on this knowledge, a commercially available compound was
selected and proven to be an activity-based probe of 3CLpro. This work
highlights the great potential of combining deep learning models and
biochemical experiments for intelligent iteration and expanding the boundaries
of drug discovery.


  We study the degree distribution of a randomly chosen vertex in a
duplication--divergence graph, under a variety of different generalizations of
the basic model of Bhan, Galas and Dewey (2002) and V\'azquez, Flammini,
Maritan and Vespignani (2003). We pay particular attention to what happens when
a non-trivial proportion of the vertices have large degrees, establishing a
central limit theorem for the logarithm of the degree distribution. Our
approach, as in Jordan (2018) and Hermann and Pfaffelhuber (2021), relies
heavily on the analysis of related birth--catastrophe processes, and couplings
are used to show that a number of different formulations of the process have
asymptotically similar expected degree distributions.


  The importance of spatial networks in the spread of an epidemic is an
essential aspect in modeling the dynamics of an infectious disease.
Additionally, any realistic data-driven model must take into account the large
uncertainty in the values reported by official sources, such as the amount of
infectious individuals. In this paper we address the above aspects through a
hyperbolic compartmental model on networks, in which nodes identify locations
of interest, such as cities or regions, and arcs represent the ensemble of main
mobility paths. The model describes the spatial movement and interactions of a
population partitioned, from an epidemiological point of view, on the basis of
an extended compartmental structure and divided into commuters, moving on a
suburban scale, and non-commuters, acting on an urban scale. Through a
diffusive rescaling, the model allows us to recover classical diffusion
equations related to commuting dynamics. The numerical solution of the
resulting multiscale hyperbolic system with uncertainty is then tackled using a
stochastic collocation approach in combination with a finite-volume IMEX
method. The ability of the model to correctly describe the spatial
heterogeneity underlying the spread of an epidemic in a realistic city network
is confirmed with a study of the outbreak of COVID-19 in Italy and its spread
in the Lombardy Region.


  Management of invasive species and pathogens requires information about the
traffic of potential vectors. Such information is often taken from vector
traffic models fitted to survey data. Here, user-specific data collected via
mobile apps offer new opportunities to obtain more accurate estimates and to
analyze how vectors' individual preferences affect propagule flows. However,
data voluntarily reported via apps may lack some trip records, adding a
significant layer of uncertainty. We show how the benefits of app-based data
can be exploited despite this drawback.
  Based on data collected via an angler app, we built a stochastic model for
angler traffic in the Canadian province Alberta. There, anglers facilitate the
spread of whirling disease, a parasite-induced fish disease. The model is
temporally and spatially explicit and accounts for individual preferences and
repeating behaviour of anglers, helping to address the problem of missing trip
records.
  We obtained estimates of angler traffic between all subbasins in Alberta. The
model's accuracy exceeds that of direct empirical estimates even when fewer
data were used to fit the model. The results indicate that anglers' local
preferences and their tendency to revisit previous destinations reduce the
number of long inter-waterbody trips potentially dispersing whirling disease.
According to our model, anglers revisit their previous destination in 64% of
their trips, making these trips irrelevant for the spread of whirling disease.
Furthermore, 54% of fishing trips end in individual-specific spatially
contained areas with mean radius of 54.7km. Finally, although the fraction of
trips that anglers report was unknown, we were able to estimate the total
yearly number of fishing trips in Alberta, matching an independent empirical
estimate.


  We propose a risk measure and construct an infectious diseases hazard map for
India. Given an outbreak location, a hazard index is assigned to each city
using an effective distance that depends on inter-city mobilities instead of
geographical distance. We demonstrate its utility using an SIR model augmented
with air, rail, and road data between top 446 cities. Simulations show that the
effective distance from outbreak location reliably predicts the time of arrival
of infection in other cities. The hazard index predictions compare well with
the observed spread of SARS-CoV-2. The hazard map can be useful in other
outbreaks also.


  The ongoing COVID-19 pandemic is being responded with various methods,
applying vaccines, experimental treatment options, total lockdowns or partial
curfews. Weekend curfews is one of the methods to reduce the amount of infected
persons and this method is practically applied in some countries such as
Turkey. In this study, the effect of weekend curfews on reducing the spread of
a contagious disease, such as COVID-19, is modeled using a Monte Carlo
algorithm with a hybrid lattice model. In the simulation setup, a fictional
country with three towns and 26,610 citizens were used as a model. Results
indicate that applying a weekend curfew reduces the active cases significantly
and is one of the efficient ways to fight the epidemic. The results also show
that applying personal precautions such as social distancing is important for
reducing the number of cases and deaths.


  Many of the models used nowadays in mathematical epidemiology, in particular
in COVID-19 research, belong to a certain sub-class of compartmental models
whose classes may be divided into three "(x, y, z)" groups, which we will call
respectively "susceptible/entrance, diseased, and output" (in the classic SIR
case, there is only one class of each type). Roughly, the ODE dynamics of these
models contain only linear terms, with the exception of products between x and
y terms. It has long been noticed that the basic reproduction number R has a
very simple formula (3.3) in terms of the matrices which define the model, and
an explicit first integral formula (3.8) is also available. These results can
be traced back at least to [ABvdD+07] and [Fen07], respectively, and may be
viewed as the "basic laws of SIR-type epidemics"; however many papers continue
to reprove them in particular instances (by the next-generation matrix method
or by direct computations, which are unnecessary). This motivated us to redraw
the attention to these basic laws and provide a self-contained reference of
related formulas for (x, y, z) models. We propose to rebaptize the class to
which they apply as matrix SIR epidemic models, abbreviated as SYR, to
emphasize the similarity to the classic SIR case. For the case of one
susceptible class, we propose to use the name SIR-PH, due to a simple
probabilistic interpretation as SIR models where the exponential infection time
has been replaced by a PH-type distribution. We note that to each SIR-PH model,
one may associate a scalar quantity Y(t) which satisfies "classic SIR
relations", see (3.8). In the case of several susceptible classes, this
generalizes to (5.10); in a future paper, we will show that (3.8), (5.10) may
be used to obtain approximate control policies which compare well with the
optimal control of the original model.


  Context. Homochirality is a generic and unique property of life on Earth and
is considered a universal and agnostic biosignature. Homochirality induces
fractional circular polarization in the incident light that it reflects.
Because this circularly polarized light can be sensed remotely, it can be one
of the most compelling candidate biosignatures in life detection missions.
While there are also other sources of circular polarization, these result in
spectrally flat signals with lower magnitude. Additionally, circular
polarization can be a valuable tool in Earth remote sensing because the
circular polarization signal directly relates to vegetation physiology. Aims.
While high-quality circular polarization measurements can be obtained in the
laboratory and under semi-static conditions in the field, there has been a
significant gap to more realistic remote sensing conditions. Methods. In this
study, we present sensitive circular spectropolarimetric measurements of
various landscape elements taken from a fast-moving helicopter. Results. We
demonstrate that during flight, within mere seconds of measurements, we can
differentiate (S/N>5) between grass fields, forests, and abiotic urban areas.
Importantly, we show that with only nonzero circular polarization as a
discriminant, photosynthetic organisms can even be measured in lakes.
Conclusions. Circular spectropolarimetry can be a powerful technique to detect
life beyond Earth, and we emphasize the potential of utilizing circular
spectropolarimetry as a remote sensing tool to characterize and monitor in
detail the vegetation physiology and terrain features of Earth itself.


  Antibodies are proteins in the immune system which bind to antigens to detect
and neutralise them. The binding sites in an antibody-antigen interaction are
known as the paratope and epitope, respectively, and the prediction of these
regions is key to vaccine and synthetic antibody development. Contrary to prior
art, we argue that paratope and epitope predictors require asymmetric
treatment, and propose distinct neural message passing architectures that are
geared towards the specific aspects of paratope and epitope prediction,
respectively. We obtain significant improvements on both tasks, setting the new
state-of-the-art and recovering favourable qualitative predictions on antigens
of relevance to COVID-19.


  Summary: More sophisticated models are needed to address problems in
bioscience, synthetic biology, and precision medicine. To help facilitate the
collaboration needed for such models, the community developed the Simulation
Experiment Description Markup Language (SED-ML), a common format for describing
simulations. However, the utility of SED-ML has been hampered by limited
support for SED-ML among modeling software tools and by different
interpretations of SED-ML among the tools that support the format. To help
modelers debug their simulations and to push the community to use SED-ML
consistently, we developed a tool for validating SED-ML files. We have used the
validator to correct the official SED-ML example files. We plan to use the
validator to correct the files in the BioModels database so that they can be
simulated. We anticipate that the validator will be a valuable tool for
developing more predictive simulations and that the validator will help
increase the adoption and interoperability of SED-ML.
  Availability: The validator is freely available as a webform, HTTP API,
command-line program, and Python package at
https://run.biosimulations.org/utils/validate and
https://pypi.org/project/biosimulators-utils. The validator is also embedded
into interfaces to 11 simulation tools. The source code is openly available as
described in the Supplementary data.
  Contact: karr@mssm.edu


  Intrinsically disordered proteins (IDPs) constitute a broad set of proteins
with few uniting and many diverging properties. IDPs-and intrinsically
disordered regions (IDRs) interspersed between folded domains-are generally
characterized as having no persistent tertiary structure; instead they
interconvert between a large number of different and often expanded structures.
IDPs and IDRs are involved in an enormously wide range of biological functions
and reveal novel mechanisms of interactions, and while they defy the common
structure-function paradigm of folded proteins, their structural preferences
and dynamics are important for their function. We here discuss open questions
in the field of IDPs and IDRs, focusing on areas where machine learning and
other computational methods play a role. We discuss computational methods aimed
to predict transiently formed local and long-range structure, including methods
for integrative structural biology. We discuss the many different ways in which
IDPs and IDRs can bind to other molecules, both via short linear motifs, as
well as in the formation of larger dynamic complexes such as biomolecular
condensates. We discuss how experiments are providing insight into such
complexes and may enable more accurate predictions. Finally, we discuss the
role of IDPs in disease and how new methods are needed to interpret the
mechanistic effects of genomic variants in IDPs.


  A variety of agent-based models has been proposed to account for the
emergence of coordinated collective behavior of animal groups from simple
interaction rules. A common, simplifying assumption of such collective movement
models, is the consideration of individual agents moving with a constant speed.
In this work we critically re-asses this assumption underlying a vast majority
of collective movement models. First, we show the omnipresent speed variability
observed in different species of live fish and artificial agents (RoboFish).
Based on theoretical considerations accounting for inertia and rotational
friction, we derive a functional dependence of the turning response of
individuals on their instantaneous speed (confirmed by experimental data). We
investigate how the interplay of variable speed and speed-dependent turning
affects self-organized collective behavior by implementing an agent-based model
which accounts for both effects. We show, that besides average speed, the
individual speed variability may have a dramatic impact on the emergent
collective dynamics, as two groups differing only in their speed variability,
and being otherwise identical in all other behavioral parameters, can be in two
fundamentally different stationary states. We find that the local coupling
between group polarization and individual speed is strongest at the
order-disorder transition. Furthermore, we demonstrate a decrease in
polarization with group size for groups of individuals with variable speed, and
a sudden decrease in mean individual speed at a critical group size (N=4 for
Voronoi interactions) linked to a topological transition from an all-to-all to
a distributed spatial interaction network. Overall, our work highlights the
importance to account for fundamental kinematic constraints in general, and
variable speed in particular, when modeling self-organized collective dynamics.


  Cell migration in confining physiological environments relies on the
concerted dynamics of several cellular components, including protrusions,
adhesions with the environment, and the cell nucleus. However, it remains
poorly understood how the dynamic interplay of these components and the cell
polarity determine the emergent migration behavior at the cellular scale. Here,
we combine data-driven inference with a mechanistic bottom-up approach to
develop a model for protrusion and polarity dynamics in confined cell
migration, revealing how the cellular dynamics adapt to confining geometries.
Specifically, we use experimental data of joint protrusion-nucleus migration
trajectories of cells on confining micropatterns to systematically determine a
mechanistic model linking the stochastic dynamics of cell polarity,
protrusions, and nucleus. This model indicates that the cellular dynamics adapt
to confining constrictions through a switch in the polarity dynamics from a
negative to a positive, self-reinforcing feedback loop. Our model further
reveals how this feedback loop leads to stereotypical cycles of
protrusion-nucleus dynamics that drive the migration of the cell through
constrictions. These cycles are disrupted upon perturbation of cytoskeletal
components, indicating that the positive feedback is controlled by cellular
migration mechanisms. Our data-driven theoretical approach therefore identifies
polarity feedback adaptation as a key mechanism in confined cell migration.


  Many biological processes discriminate between correct and incorrect
substrates through the kinetic proofreading mechanism which enables lower error
at the cost of higher energy dissipation. Elucidating physicochemical
constraints for global minimization of dissipation and error is important for
understanding enzyme evolution. Here, we identify theoretically a fundamental
error-cost bound which tightly constrains the performance of proofreading
networks under any parameter variations preserving the rate discrimination
between substrates. The bound is kinetically controlled, i.e. completely
determined by the difference between the transition state energies on the
underlying free energy landscape. The importance of the bound is analyzed for
three biological processes. DNA replication by T7 DNA polymerase is shown to be
nearly optimized, i.e. its kinetic parameters place it in the immediate
proximity of the error-cost bound. The isoleucyl-tRNA synthetase (IleRS) of E.
coli also operates close to the bound, but further optimization is prevented by
the need for reaction speed. In contrast, E. coli ribosome operates in a
high-dissipation regime, potentially in order to speed up protein production.
Together, these findings establish a fundamental error-dissipation relation in
biological proofreading networks and provide a theoretical framework for
studying error-dissipation trade-off in other systems with biological
discrimination.


  We report a droplet microfluidic method to target and sort individual cells
directly from complex microbiome samples, and to prepare these cells for bulk
whole genome sequencing without cultivation. We characterize this approach by
recovering bacteria spiked into human stool samples at a ratio as low as 1:250
and by successfully enriching endogenous Bacteroides vulgatus to the level
required for de-novo assembly of high-quality genomes. While microbiome strains
are increasingly demanded for biomedical applications, the vast majority of
species and strains are uncultivated and without reference genomes. We address
this shortcoming by encapsulating complex microbiome samples directly into
microfluidic droplets and amplify a target-specific genomic fragment using a
custom molecular TaqMan probe. We separate those positive droplets by droplet
sorting, selectively enriching single target strain cells. Finally, we present
a protocol to purify the genomic DNA while specifically removing amplicons and
cell debris for high-quality genome sequencing.


  Health benefits of EGCG are well established; however, the mechanisms of EGCG
action are not completely understood. In our previous study we discovered
solvatochromism of EGCG fluorescence and described the dependence of
fluorescence maxima on solvent polarity. We also noted that fluorescence
intensity (FI) depends on solvent. The goal of this study was to gain insights
into how the protonic properties of the environment affect FI. We demonstrate
that 1) FI of EGCG inversely correlates with the autoprotolysis constant (Kap)
of the solvents, 2) HCl decreases, and NaOH increases FI of EGCG in water, and
3) NaOH evokes slow (~10 min time scale) transient changes in FI of EGCG. We
conclude that EGCG fluorescence depends on Kap, i.e. protonicity of the
environment, that is useful for differentiating EGCG in protic aqueous
environment at physiological pH~7.0{\div}7.4, where EGCG fluorescence is
substantially quenched, from EGCG in aprotic environments, where EGCG
fluorescence significantly increases. Thus, this property of EGCG fluorescence
is useful for investigation of specific EGCG interactions within more aprotic
protein binding sites.


  In this work, the SIR epidemiological model is reformulated so to highlight
the important {\em effective reproduction number}, as well as to account for
the {\em generation time}, inverse of the {\em incidence rate}, and the {\em
infectious period} (or {\em removal period}), inverse of the {\em removal
rate}. The aim is to check whether the relationships the model poses among the
various observables are actually found in the data. The study case of the
second through the third wave of the Covid-19 pandemic in Italy is taken. Given
its scale invariance, initially the model is tested with reference to the curve
of swab-confirmed infectious individuals only. It is found to match the data if
the given curve of the {\em removed} (that is healed or deceased) individuals
is assumed underestimated by a factor of about 3 together with other related
curves. Contextually, the {\em generation time} and the {\em removal period},
as well as the {\em effective reproduction number}, are obtained fitting the
SIR equations to the data; the outcomes prove to be in good agreement with
those of other works. Then, using knowledge of the proportion of Covid-19
transmissions likely occurring from individuals who didn't develop symptoms,
thus mainly undetected, an estimate of the {\em "true numbers"} of the epidemic
is obtained, looking also in good agreement with results from other, completely
different works. The line of this work is new and the procedures are
computationally really inexpensive.


  Predicting DNA-protein binding is an important and classic problem in
bioinformatics. Convolutional neural networks have outperformed conventional
methods in modeling the sequence specificity of DNA-protein binding. However,
none of the studies has utilized graph convolutional networks for motif
inference. In this work, we propose to use graph convolutional networks for
motif inference. We build a sequence k-mer graph for the whole dataset based on
k-mer co-occurrence and k-mer sequence relationship and then learn DNA Graph
Convolutional Network (DNA-GCN) for the whole dataset. Our DNA-GCN is
initialized with a one-hot representation for all nodes, and it then jointly
learns the embeddings for both k-mers and sequences, as supervised by the known
labels of sequences. We evaluate our model on 50 datasets from ENCODE. DNA-GCN
shows its competitive performance compared with the baseline model. Besides, we
analyze our model and design several different architectures to help fit
different datasets.


  Cyclic dominance of competing species is an intensively used working
hypothesis to explain biodiversity in certain living systems, where the
evolutionary selection principle would dictate a single victor otherwise.
Technically the May--Leonard models offer a mathematical framework to describe
the mentioned non-transitive interaction of competing species when individual
movement is also considered in a spatial system. Emerging rotating spirals
composed by the competing species are frequently observed character of the
resulting patterns. But how do these spiraling patterns change when we vary the
external environment which affects the general vitality of individuals?
Motivated by this question we suggest an off-lattice version of the tradition
May--Leonard model which allows us to change the actual state of the
environment gradually. This can be done by introducing a local carrying
capacity parameter which value can be varied gently in an off-lattice
environment. Our results support a previous analysis obtained in a more
intricate metapopulation model and we show that the well-known rotating spirals
become evident in a benign environment when the general density of the
population is high. The accompanying time-dependent oscillation of competing
species can also be detected where the amplitude and the frequency show a
scaling law of the parameter that characterizes the state of the environment.
These observations highlight that the assumed non-transitive interaction alone
is insufficient condition to maintain biodiversity safely, but the actual state
of the environment, which characterizes the general living conditions, also
plays a decisive role on the evolution of related systems.


  Models at various levels of resolution are commonly used, both for forest
management and in ecological research. They all have comparative advantages and
disadvantages, making desirable a better understanding of the relationships
between the various approaches. It is found that accounting for crown and root
plasticity creates more realistic links between spatial and non-spatial models
than simply ignoring spatial structure. The article reviews also the connection
between distance-independent models and size distributions, and how
distributions evolve over time and relate to whole-stand descriptions. In
addition, some ways in which stand-level knowledge feeds back into detailed
individual-tree formulations are demonstrated. The presentation intends to be
accessible to non-specialists.
  Study implications: Introducing plasticity improves the representation of
physio-ecological processes in spatial modelling. Plasticity explains in part
the practical success of distance-independent models. The nature of size
distributions and their relationship to individual-tree and whole-stand models
are discussed. I point out limitations of various approaches and questions for
future research.


  Enzyme is the major workhorse to carry out the diverse cellular functions. It
catalyzes the biological reactions with a high specificity, with its topology
playing a crucial role. For ecologically safe production of numerous
bioproducts including drugs and chemicals, we have been striving to design the
industrially useful enzyme molecules with highly improved catalytic capability.
As the sequence space is enormous for an enzyme, its quick and effective
exploration is quite improbable for the mutagenesis studies whose accuracy is
greatly reliant on the prior information of the mutated sites and the extent of
rigorous screening of the mutant libraries. Although directed evolution methods
significantly aid the construction of a functionally improved molecule, their
credibility depends on the successful excavation of the functionally similar
sequence space in the available databases, encompassing billions of proteins.
As deep learning methods aid us to extensively uncover the underlying network
of all the key catalytic positions without any experimental data, their
implementation has reliably increased the accuracy of directed evolution. The
chapter comprehensively explains data mining and deep learning methods to
further showcase their importance in enzyme engineering methods. The key
biological and algorithmic limitations of these deep learning methodologies are
lastly highlighted.


  Tyshkovskiy and Panchin have recently published a commentary on our paper in
which they outline several "points of disagreement with the Segreto/Deigin
hypothesis". As our paper is titled "The genetic structure of SARS-CoV-2 does
not rule out a laboratory origin", points of disagreement should provide
evidence that rules out a laboratory origin. However, Tyshkovskiy and Panchin
provide no such evidence and instead attempt to criticize our arguments that
highlight aspects of SARS-CoV-2 that could be consistent with the lab leak
hypothesis. Strikingly, Tyshkovskiy and Panchin's main point of criticism is
based on a false premise that we have claimed RaTG13 to be a direct progenitor
of SARS-CoV-2, and their other points of criticism are either incorrect or
irrelevant to our hypotheses. Thus, the genetic structure of SARS-CoV-2 remains
consistent with both natural or laboratory origin, which means that both the
zoonotic and the lab leak hypothesis need to be investigated equally
thoroughly.


  Infectious diseases are caused by pathogenic microorganisms and can spread
through different ways. Mathematical models and computational simulation have
been used extensively to investigate the transmission and spread of infectious
diseases. In other words, mathematical model simulation can be used to analyse
the dynamics of infectious diseases, aiming to understand the effects and how
to control the spread. In general, these models are based on compartments,
where each compartment contains individuals with the same characteristics, such
as susceptible, exposed, infected, and recovered. In this paper, we cast
further light on some classical epidemic models, reporting possible outcomes
from numerical simulation. Furthermore, we provide routines in a repository for
simulations.


  We developed Distilled Graph Attention Policy Network (DGAPN), a
reinforcement learning model to generate novel graph-structured chemical
representations that optimize user-defined objectives by efficiently navigating
a physically constrained domain. The framework is examined on the task of
generating molecules that are designed to bind, noncovalently, to functional
sites of SARS-CoV-2 proteins. We present a spatial Graph Attention (sGAT)
mechanism that leverages self-attention over both node and edge attributes as
well as encoding the spatial structure -- this capability is of considerable
interest in synthetic biology and drug discovery. An attentional policy network
is introduced to learn the decision rules for a dynamic, fragment-based
chemical environment, and state-of-the-art policy gradient techniques are
employed to train the network with stability. Exploration is driven by the
stochasticity of the action space design and the innovation reward bonuses
learned and proposed by random network distillation. In experiments, our
framework achieved outstanding results compared to state-of-the-art algorithms,
while reducing the complexity of paths to chemical synthesis.


  Spatial context is central to understanding health and disease. Yet reference
protein interaction networks lack such contextualization, thereby limiting the
study of where protein interactions likely occur in the human body and how they
may be altered in disease. Contextualized protein interactions could better
characterize genes with disease-specific interactions and elucidate diseases'
manifestation in specific cell types. Here, we introduce AWARE, a graph neural
message passing approach to inject cellular and tissue context into protein
embeddings. AWARE optimizes for a multi-scale embedding space, whose structure
reflects network topology at a single-cell resolution. We construct a
multi-scale network of the Human Cell Atlas and apply AWARE to learn protein,
cell type, and tissue embeddings that uphold cell type and tissue hierarchies.
We demonstrate AWARE's utility on the novel task of predicting whether a
protein is altered in disease and where that association most likely manifests
in the human body. To this end, AWARE outperforms generic embeddings without
contextual information by at least 12.5%, showing AWARE's potential to reveal
context-dependent roles of proteins in disease.


  Statistical Physics has proved essential to analyze multi-agent environments.
Motivated by the empirical observation of various non-equilibrium features in
Barro Colorado and other ecological systems, we analyze a plant-species
abundance model, presenting analytical evidence of scale-invariant plant
clusters and non-trivial emergent modular correlations. Such first theoretical
confirmation of a scale-invariant region, based on percolation processes,
reproduces the key features in actual ecological ecosystems and can confer the
most stable equilibrium for ecosystems with vast biodiversity.


  During their evolution, proteins explore sequence space via an interplay
between random mutations and phenotypic selection. Here we build upon recent
progress in reconstructing data-driven fitness landscapes for families of
homologous proteins, to propose stochastic models of experimental protein
evolution. These models predict quantitatively important features of
experimentally evolved sequence libraries, like fitness distributions and
position-specific mutational spectra. They also allow us to efficiently
simulate sequence libraries for a vast array of combinations of experimental
parameters like sequence divergence, selection strength and library size. We
showcase the potential of the approach in re-analyzing two recent experiments
to determine protein structure from signals of epistasis emerging in
experimental sequence libraries. To be detectable, these signals require
sufficiently large and sufficiently diverged libraries. Our modeling framework
offers a quantitative explanation for the variable success of recently
published experiments. Furthermore, we can forecast the outcome of time- and
resource-intensive evolution experiments, opening thereby a way to
computationally optimize experimental protocols.


  Since November 6$^{\mathrm{th}}$, 2020, Italian regions have been classified
according to four levels, corresponding to specific risk scenarios, for which
specific restrictive measures have been foreseen. By analyzing the time
evolution of the reproduction number $R_t$, we estimate how much different
restrictive measures affect $R_t$, and we quantify the combined effect of the
diffusion of virus variants and the beginning of the vaccination campaign upon
the $R_t$ trend. We also compute the time delay between implementation of
restrictive measures and the resulting effects. Three different models to
describe the effects of restrictive measures are discussed and the results are
cross-checked with two different algorithms for the computation of $R_t$.


  Motivation: SBML is the most widespread language for the definition of
biochemical models. Although dozens of SBML simulators are available, there is
a general lack of support to the integration of SBML models within
open-standard general-purpose simulation ecosystems. This hinders co-simulation
and integration of SBML models within larger model networks, in order to, e.g.
enable in silico clinical trials of drugs, pharmacological protocols, or
engineering artefacts such as biomedical devices against Virtual Physiological
Human models. Modelica is one of the most popular existing open-standard
general-purpose simulation languages, supported by many simulators. Modelica
models are especially suited for the definition of complex networks of
heterogeneous models from virtually all application domains. Models written in
Modelica (and in 100+ other languages) can be readily exported into black-box
Functional Mock-Up Units (FMUs), and seamlessly co-simulated and integrated
into larger model networks within open-standard language-independent simulation
ecosystems.
  Results: In order to enable SBML model integration within heterogeneous model
networks, we present SBML2Modelica, a software system translating SBML models
into well-structured, user-intelligible, easily modifiable Modelica models.
SBML2Modelica is SBML Level 3 Version 2-compliant and succeeds on 96.47% of the
SBML Test Suite Core (with a few rare, intricate and easily avoidable
combinations of constructs unsupported and cleanly signalled to the user). Our
experimental campaign on 613 models from the BioModels database (with up to
5438 variables) shows that the major open-source (general-purpose) Modelica and
FMU simulators achieve performance comparable to state-of-the-art specialized
SBML simulators.
  Availability and implementation: https://bitbucket.org/mclab/sbml2modelica


  Rooted and ranked binary trees are mathematical objects of great importance
used to model hierarchical data and evolutionary relationships with
applications in many fields including evolutionary biology and genetic
epidemiology. Bayesian phylogenetic inference usually explore the posterior
distribution of trees via Markov Chain Monte Carlo methods, however assessing
uncertainty and summarizing distributions or samples of such trees remains
challenging. While labelled phylogenetic trees have been extensively studied,
relatively less literature exists for unlabelled trees which are increasingly
useful, for example when one seeks to summarize samples of trees obtained with
different methods, or from different samples and environments, and wishes to
assess stability and generalizability of these summaries. In our paper, we
exploit recently proposed distance metrics of unlabelled ranked binary trees
and unlabelled ranked genealogies (equipped with branch lengths) to define the
Frechet mean and variance as summaries of these tree distributions. We provide
an efficient combinatorial optimization algorithm for computing the Frechet
mean from a sample of or distribution on unlabelled ranked tree shapes and
unlabelled ranked genealogies. We show the applicability of our summary
statistics for studying popular tree distributions and for comparing the
SARS-CoV-2 evolutionary trees across different locations during the COVID-19
epidemic in 2020.


  The relationship of network structure and dynamics is one of most extensively
investigated problems in the theory of complex systems of the last years.
Understanding this relationship is of relevance to a range of disciplines --
from Neuroscience to Geomorphology. A major strategy of investigating this
relationship is the quantitative comparison of a representation of network
architecture (structural connectivity) with a (network) representation of the
dynamics (functional connectivity). Analysing such SC/FC relationships has over
the past years contributed substantially to our understanding of the functional
role of network properties, such as modularity, hierarchical organization, hubs
and cycles.
  Here, we show that one can distinguish two classes of functional connectivity
-- one based on simultaneous activity (co-activity) of nodes the other based on
sequential activity of nodes. We delineate these two classes in different
categories of dynamical processes -- excitations, regular and chaotic
oscillators -- and provide examples for SC/FC correlations of both classes in
each of these models. We expand the theoretical view of the SC/FC
relationships, with conceptual instances of the SC and the two classes of FC
for various application scenarios in Geomorphology, Freshwater Ecology, Systems
Biology, Neuroscience and Social-Ecological Systems.
  Seeing the organization of a dynamical processes in a network either as
governed by co-activity or by sequential activity allows us to bring some order
in the myriad of observations relating structure and function of complex
networks.


  We report on three pre-registered studies testing whether people in the
position of describing a decision problem to decision-makers exploit this
opportunity for their benefit, by choosing descriptions that may be potentially
beneficial for themselves. In Study 1, recipients of an extreme dictator game
(where dictators can either take the whole pie for themselves or give it
entirely to the receiver) are asked to choose the instructions used to
introduce the game to dictators, among six different instructions that are
known from previous research to affect dictators' decisions. The results
demonstrate that some dictator game recipients tend to choose instructions that
make them more likely to receive a higher payoff. Study 2 shows that people who
choose descriptions that make them more likely to receive a higher payoff
indeed believe that they will receive a higher payoff. Study 3 shows that
receivers are more likely than dictators to choose these descriptions. In sum,
our work suggests that some people choose descriptions that are beneficial to
themselves; we also found some evidence that deliberative thinking and young
age are associated with this tendency.


  Networks of genetic expression can be modelled by hypergraphs with the
additional structure that real coefficients are given to each vertex-edge
incidence. The spectra, i.e. the multiset of the eigenvalues, of such
hypergraphs, are known to encode structural information of the data. We show
how these spectra can be used, in particular, in order to give an estimation of
cellular redundancy, a novel measure of gene expression heterogeneity, of the
network. We analyze some simulated and real data sets of gene expression for
illustrating the new method proposed here.


  Natural selection is general and powerful concept not only to explain
evolutionary processes of biological organisms but also to design engineering
systems such as genetic algorithms and particle filters. There is a surge of
interest, both from biology and engineering, in considering natural selection
of intellectual agents that can learn individually. Learning by individual
agents of better behaviors for survival may accelerate the evolutionary
processes by natural selection. We have accumulating pieces of evidence that
organisms can transmit its information to the next generation via epigenetic
states or memes. Also, such idea is important for engineering applications. To
accelerate the evolutionary process, an agent should change their strategy so
that the population fitness increases the most. Equivalently, an agent should
update the strategy towards a gradient of the population fitness. However, it
has not yet been clarified whether and how an agent can estimate the gradient
and accelerate the evolutionary process. We also lack methodology to quantify
the acceleration to understand and predict the impact of learning. In this
paper, we address these problems. We show that an learning agent can accelerate
the evolutionary process by proposing ancestral learning, which uses the
information transmitted from the ancestor (ancestral information). We next show
that the ancestral information is sufficient to estimate the gradient. In
particular, learning can accelerate the evolutionary process without
communications between agents. Finally, to quantify the acceleration, we extend
the Fisher's fundamental theorem (FF-thm) for natural selection to ancestral
learning. Our extended FF-thm relates the acceleration of the evolutionary
process to the variety of individual fitness of the agent. By the theorem, we
can quantitatively understand when and why learning is beneficial.


  This paper is a short extension of our previous paper [arXiv:2004.06033]
about the use of the Test-Negative design to study risk factors for COVID-19
[See: PubMed and ArXiv reference below] Reason for the extension is that the
conditions under which people refer themselves for testing have greatly
changed: originally, in most countries priority was given to people with
symptoms, but nowadays people without symptoms are also tested for different
reasons, e.g., during contact tracing, or to be allowed on an (international)
flight. Interestingly, this opens new possibilities to separately investigate
risk factors for infection and risk factors for becoming diseased. To use this
new situation to best effect, one has to think carefully about how to elucidate
the different reasons for testing and what analyses one might do with the
different groups.


  The rigidity and flexibility of homologous psychrophilic(P), mesophilic(M)
and thermophilic(T) proteins have been investigated at the global and local
levels in terms of packing factor and atomic fluctuations obtained from
B-factors. For comparison of atomic fluctuations, correction of errors by
considering errors in B-factors from all sources in a consolidated manner and
conversion of the fluctuations to the same temperature have been suggested and
validated. Results indicate no differences in the global values like average
packing factor among the three classes of protein homologs but at local levels
there are differences. Comparison of homologous proteins triplets show that the
average atomic fluctuations at a given temperature obey the order P>M>T.
Packing factors and the atomic fluctuations are anti-correlated suggesting that
altering the rigidity of the active site might be a potential strategy to make
tailor made psychrophilic or thermophilic proteins from their mesophilic
homologs.


  Representing and reasoning about 3D structures of macromolecules is emerging
as a distinct challenge in machine learning. Here, we extend recent work on
geometric vector perceptrons and apply equivariant graph neural networks to a
wide range of tasks from structural biology. Our method outperforms all
reference architectures on three out of eight tasks in the ATOM3D benchmark, is
tied for first on two others, and is competitive with equivariant networks
using higher-order representations and spherical harmonic convolutions. In
addition, we demonstrate that transfer learning can further improve performance
on certain downstream tasks. Code is available at
https://github.com/drorlab/gvp-pytorch.


  The temporal dynamics of membrane voltage changes in neurons is controlled by
ionic currents. These currents are characterized by two main properties:
conductance and kinetics. The hyperpolarization-activated current ($I_{\rm h}$)
strongly modulates subthreshold potential changes by shortening the excitatory
postsynaptic potentials and decreasing their temporal summation. Whereas the
shortening of the synaptic potentials caused by the $I_{\rm h}$ conductance is
well understood, the role of the $I_{\rm h}$ kinetics remains unclear. Here, we
use a model of the $I_{\rm h}$ current model with either fast or slow kinetics
to determine its influence on the membrane time constant ($\tau_m$) of a CA1
pyramidal cell model. Our simulation results show that the $I_{\rm h}$ with
fast kinetics decreases $\tau_m$ and attenuates and shortens the excitatory
postsynaptic potentials more than the slow $I_{\rm h}$. We conclude that the
$I_{\rm h}$ activation kinetics is able to modulate $\tau_m$ and the temporal
properties of excitatory postsynaptic potentials (EPSPs) in CA1 pyramidal
cells. In order to elucidate the mechanisms by which $I_{\rm h}$ kinetics
controls $\tau_m$, we propose a new concept called "time scaling factor". Our
main finding is that the $I_{\rm h}$ kinetics influences $\tau_m$ by modulating
the contribution of the $I_{\rm h}$ derivative conductance to $\tau_m$.


  The lactose uptake-pathway of E. coli is a paradigmatic example of
multistability in gene-regulatory circuits. In the induced state of the
lac-pathway, the genes comprising the lac-operon are transcribed, leading to
the production of proteins which import and metabolize lactose. In the
uninduced state, a stable repressor-DNA loop frequently blocks the
transcription of the lac-genes. Transitions from one phenotypic state to the
other are driven by fluctuations, which arise from the random timing of the
binding of ligands and proteins. This stochasticity affects transcription and
translation, and ultimately molecular copy numbers. Our aim is to understand
the transition from the induced to the uninduced state of the lac-operon. We
use a detailed computational model to show that repressor-operator
binding/unbinding, fluctuations in the total number of repressors, and
inducer-repressor binding/unbinding all play a role in this transition. Based
on the timescales on which these processes operate, we construct a minimal
model of the transition to the uninduced state and compare the results with
simulations and experimental observations. The induced state turns out to be
very stable, with a transition rate to the uninduced state lower than $2 \times
10^{-9}$ per minute. In contrast to the transition to the induced state, the
transition to the uninduced state is well described in terms of a 2D diffusive
system crossing a barrier, with the diffusion rates emerging from a model of
repressor unbinding.


  We explored the dynamical and structural effects of actin-related proteins
2/3 (Arp2/3) on actomyosin networks using mechanochemical simulations of active
matter networks. At a nanoscale, the Arp2/3 complex alters the topology of
actomyosin by nucleating a daughter filament at an angle to a mother filament.
At a subcellular scale, they orchestrate the formation of branched actomyosin
network. Using a coarse-grained approach, we sought to understand how an
actomyosin network temporally and spatially reorganizes itself by varying the
concentration of the Arp2/3 complexes. Driven by the motor dynamics, the
network stalls at a high concentration of Arp2/3 and contracts at a low Arp2/3
concentration. At an intermediate Arp2/3 concentration, however, the actomyosin
network is formed by loosely connected clusters that may collapse suddenly when
driven by motors. This physical phenomenon is called an "avalanche" largely due
to the marginal instability inherent from the morphology of a branched
actomyosin network when the Arp2/3 complex is present. While embracing the data
science approaches, we unveiled the higher-order patterns in the branched
actomyosin networks and discovered a sudden change in the "social" network
topology of the actomyosin. This is a new type of avalanches in addition to the
two types of avalanches associated with a sudden change in the size or the
shape of the whole actomyosin network as shown in the previous investigation.
Our new finding promotes the importance of using network theory and machine
learning models to forecast avalanches in actomyosin networks. The mechanisms
of the Arp2/3 complexes in shaping the architecture of branched actomyosin
networks obtained in this paper will help us better understand the emergent
reorganization of the topology in dense actomyosin networks that are difficult
to detect in experiments.


  This document aims to estimate and describe the effects of the social
distancing measures implemented in several countries with the expectancy of
controlling the spread of COVID-19. The procedure relies on the classic
Susceptible-Infected-Removed (SIR) model, which is modified to incorporate a
permissiveness index, representing the isolation achieved by the social
distancing and the future development of vaccination campaigns and allowing the
math model to reproduce more than one infection wave. The adjusted SIR models
are used to study the compromise between the economy's reactivation and the
resulting infection spreading increase. The document presents a
graphical-abacus that describes the convenience of progressively relax social
distancing measures while a feasible vaccination campaign develops


  We consider an epidemic model of SIR type set on a homogeneous tree and
investigate the spreading properties of the epidemic as a function of the
degree of the tree, the intrinsic basic reproduction number and the strength of
the interactions within the population of infected individuals. When the degree
is one, the homogeneous tree is nothing but the standard lattice on the
integers and our model reduces to a SIR model with discrete diffusion for which
the spreading properties are very similar to the continuous case. On the other
hand, when the degree is larger than two, we observe some new features in the
spreading properties. Most notably, there exists a critical value of the
strength of interactions above which spreading of the epidemic in the tree is
no longer possible.


  How and where proteins interface with one another can ultimately impact the
proteins' functions along with a range of other biological processes. As such,
precise computational methods for protein interface prediction (PIP) come
highly sought after as they could yield significant advances in drug discovery
and design as well as protein function analysis. However, the traditional
benchmark dataset for this task, Docking Benchmark 5 (DB5), contains only a
modest 230 complexes for training, validating, and testing different machine
learning algorithms. In this work, we expand on a dataset recently introduced
for this task, the Database of Interacting Protein Structures (DIPS), to
present DIPS-Plus, an enhanced, feature-rich dataset of 42,112 complexes for
geometric deep learning of protein interfaces. The previous version of DIPS
contains only the Cartesian coordinates and types of the atoms comprising a
given protein complex, whereas DIPS-Plus now includes a plethora of new
residue-level features including protrusion indices, half-sphere amino acid
compositions, and new profile hidden Markov model (HMM)-based sequence features
for each amino acid, giving researchers a large, well-curated feature bank for
training protein interface prediction methods. We demonstrate through rigorous
benchmarks that training an existing state-of-the-art (SOTA) model for PIP on
DIPS-Plus yields SOTA results, surpassing the performance of all other models
trained on residue-level and atom-level encodings of protein complexes to date.


  The conflict between individual and collective interests is in the heart of
every social dilemmas established by evolutionary game theory. We cannot avoid
these conflicts but sometimes we may choose which interaction framework to use
as a battlefield. For instance some people like to be part of a larger group
while other persons prefer to interact in a more personalized, individual way.
Both attitudes can be formulated via appropriately chosen traditional games. In
particular, the prisoner's dilemma game is based on pair interaction while the
public goods game represents multi-point interactions of group members. To
reveal the possible advantage of a certain attitude we extend these models by
allowing players not simply to change their strategies but also let them to
vary their attitudes for a higher individual income. We show that both
attitudes could be the winner at a specific parameter value. Interestingly,
however, the subtle interplay between different states may result in a
counterintuitive evolutionary outcome where the increase of the multiplication
factor of public goods game drives the population to a fully defector state. We
point out that the accompanying pattern formation can only be understood via
the multipoint or multi-player interactions of different microscopic states
where the vicinity of a particular state may influence the relation of two
other competitors.


  The titled subject has attracted much interest. Here we summarize the
substantial results obtained by a physical model of protein evolution based on
hydropathic domain dynamics. In a recent Letter eighteen biologists suggested
that the titled subject should be studied in a way inclusive of broad expertise
(1). There is an even broader view that has been developed over several decades
by physicists (2,3). This view is based on analyzing amino acid sequences of
proteins. These sequences are now available on-line at Uniprot, and represent a
treasure-trove of data (4).


  We describe the population-based SEIR (susceptible, exposed, infected,
removed) model developed by the Irish Epidemiological Modelling Advisory Group
(IEMAG), which advises the Irish government on COVID-19 responses. The model
assumes a time-varying effective contact rate (equivalently, a time-varying
reproduction number) to model the effect of non-pharmaceutical interventions. A
crucial technical challenge in applying such models is their accurate
calibration to observed data, e.g., to the daily number of confirmed new cases,
as the past history of the disease strongly affects predictions of future
scenarios. We demonstrate an approach based on inversion of the SEIR equations
in conjunction with statistical modelling and spline-fitting of the data, to
produce a robust methodology for calibration of a wide class of models of this
type.


  COVID-19 pandemic has been raging all around the world for almost a year now,
as of November 1, 2020. In this paper, we try to analyze the variation of the
COVID-19 pandemic in different countries in the light of some modifications to
the susceptible-infected-recovered (SIR) model. The SIR model was modified by
taking time-dependent rate parameters. From this modified SIR model, the basic
reproduction number, effective reproduction number, herd immunity, and herd
immunity threshold are redefined. The re-outbreak of the COVID-19 is a real
threat to various countries. We have used the above-mentioned quantities to
find the reasons behind the re-outbreak of this disease. Also, the
effectiveness of herd immunity to prevent an epidemic has been analyzed with
respect to this model. We have also tried to show that there are certain
universal aspects in the spread and containment of the disease in various
countries for a short period of time. Finally, we have also analyzed the
current pandemic situation in India and have attempted to discuss the
possibilities in order to predict its future behaviour using our model.


  Organelles of optimum size are crucial for proper functioning of a living
cell. The cell employs various mechanisms for actively sensing and controlling
the size of its organelles. Recently Bauer et al have opened a new research
frontier in the field of subcellular size control by shedding light on the
noise and fluctuations of organelles of controlled size. Taking eukaryotic
flagellum as a model organelle, which is quite popular for such studies because
of its linear geometry and dynamic nature, Bauer et al have analysed the nature
of fluctuations of its length. Here we summarize the key questions and the
fundamental importance of the recent developments. Although our attention is
focussed here mainly on the experimental and theoretical works on eukaryotic
flagellum, the ideas are general and applicable to wide varieties of cell
organelle.


  After the introduction of drastic containment measures aimed at stopping the
epidemic contagion from SARS-CoV2, many governments have adopted a strategy
based on a periodic relaxation of such measures in the face of a severe
economic crisis caused by lockdowns. Assessing the impact of such openings in
relation to the risk of a resumption of the spread of the disease is an
extremely difficult problem due to the many unknowns concerning the actual
number of people infected, the actual reproduction number and infection
fatality rate of the disease. In this work, starting from a compartmental model
with a social structure and stochastic inputs, we derive models with multiple
feedback controls depending on the social activities that allow to assess the
impact of a selective relaxation of the containment measures in the presence of
uncertain data. Specific contact patterns in the home, work, school and other
locations have been considered. Results from different scenarios concerning the
first wave of the epidemic in some major countries, including Germany, France,
Italy, Spain, the United Kingdom and the United States, are presented and
discussed.