  The role of dispersal on the stability and synchrony of a metacommunity is a
topic of considerable interest in theoretical ecology. Dispersal is known to
promote both synchrony, which enhances the likelihood of extinction, and
spatial heterogeneity, which favors the persistence of the population. Several
efforts have been made to understand the effect of diverse variants of
dispersal in the spatially distributed ecological community. Despite the
environmental change strongly affect the dispersal, the effects of controlled
dispersal on the metacommunity stability and their persistence remain unknown.
We study the influence of limiting the immigration using two patch
prey-predator metacommunity at both local and spatial scales. We find that the
spread of the inhomogeneous stable steady states (asynchronous states)
decreases monotonically upon limiting the predator dispersal. Nevertheless, at
the local scale, the spread of the inhomogeneous steady states increases up to
a critical value of the limiting factor, favoring the metacommunity
persistence, and then starts decreasing for further decrease in the limiting
factor with varying local interaction. Interestingly, limiting the prey
dispersal promotes inhomogeneous steady states in a large region of the
parameter space, thereby increasing the metacommunity persistence both at
spatial and local scales. Further, we show similar qualitative dynamics in an
entire class of complex networks consisting of a large number of patches. We
also deduce various bifurcation curves and stability condition for the
inhomogeneous steady states, which we find to agree well with the simulation
results. Thus, our findings on the effect of the limiting dispersal can help to
develop conservation measures for ecological community.


  Reputation is one of key mechanisms to maintain human cooperation, but its
analysis gets complicated if we consider the possibility that reputation does
not reach consensus because of erroneous assessment. The difficulty is
alleviated if we assume that reputation and cooperation do not take binary
values but have continuous spectra so that disagreement over reputation can be
analysed in a perturbative way. In this work, we carry out the analysis by
expanding the dynamics of reputation to the second order of perturbation under
the assumption that everyone initially cooperates with good reputation. The
second-order theory clarifies the difference between Image Scoring and Simple
Standing in that punishment for defection against a well-reputed player should
be regarded as good for maintaining cooperation. Moreover, comparison among the
leading eight shows that the stabilizing effect of justified punishment weakens
if cooperation between two ill-reputed players is regarded as bad. Our analysis
thus explains how Simple Standing achieves a high level of stability by
permitting justified punishment and also by disregarding irrelevant information
in assessing cooperation. This observation suggests which factors affect the
stability of a social norm when reputation can be perturbed by noise.


  Design of de novo biological sequences with desired properties, like protein
and DNA sequences, often involves an active loop with several rounds of
molecule ideation and expensive wet-lab evaluations. These experiments can
consist of multiple stages, with increasing levels of precision and cost of
evaluation, where candidates are filtered. This makes the diversity of proposed
candidates a key consideration in the ideation phase. In this work, we propose
an active learning algorithm leveraging epistemic uncertainty estimation and
the recently proposed GFlowNets as a generator of diverse candidate solutions,
with the objective to obtain a diverse batch of useful (as defined by some
utility function, for example, the predicted anti-microbial activity of a
peptide) and informative candidates after each round. We also propose a scheme
to incorporate existing labeled datasets of candidates, in addition to a reward
function, to speed up learning in GFlowNets. We present empirical results on
several biological sequence design tasks, and we find that our method generates
more diverse and novel batches with high scoring candidates compared to
existing approaches.


  The area of population dynamics has a rich history of the development and
analysis of models of biological and social phenomena using ordinary
differential equations. This paper describes a method for understanding the
influence one variable exerts on another in such models as a force, with the
relative effects of these forces providing a narrative explanation of the
curvature in variable behaviour. Using the stock/flow form of a model, a
symbolic notation is developed that identifies the forces with the causal
pathways of the model's feedback loops. A force is measured by its impact,
defined as the ratio of acceleration to rate of change, computed by
differentiation along its associated pathway between variables. Different
phases of force dominance are determined to enhance the standard stability
analysis of the models, providing an explanation of model behaviour in
Newtonian mechanical terms. The concepts developed are applied to well-known
models from mathematical biology: the Spruce Budworm model, where force
dominance identifies scenarios that give clarity to intervention points; and
the Lotka-Volterra predator-prey model where the analysis highlights the
importance of dissipative forces in achieving stability. Conclusions are drawn
on the explanatory power of this approach, with suggestions made for future
work.


  In many organisms, cell division is driven by the constriction of a
cytokinetic ring, which consists of actin filaments and crosslinking proteins.
While it has long been believed that the constriction is driven by motor
proteins, it has recently been discovered that passive crosslinkers that do not
turn over fuel are able to generate enough force to constrict actin filament
rings. To study the ring constriction dynamics, we develop a model that
includes the driving force of crosslinker condensation and the opposing forces
of friction and filament bending. We analyze the constriction force as a
function of ring topology and crosslinker concentration, and predict forces
that are sufficient to constrict an unadorned plasma membrane. Our model also
predicts that actin-filament sliding arises from an interplay between filament
rotation and crosslinker hopping, producing frictional forces that are low
compared to those of crosslinker-mediated microtubule sliding.


  Electrospray ion-beam deposition (ES-IBD) is a versatile tool to study
structure and reactivity of molecules from small metal clusters to large
protein assemblies. It brings molecules gently into the gas phase where they
can be accurately manipulated and purified, followed by controlled deposition
onto various substrates. In combination with imaging techniques, direct
structural information of well-defined molecules can be obtained, which is
essential to test and interpret results from indirect mass spectrometry
techniques. To date, ion-beam deposition experiments are limited to a small
number of custom instruments worldwide, and there are no commercial
alternatives. Here we present a module that adds ion-beam deposition
capabilities to a popular commercial MS platform (Thermo
Scientific$^{\mathrm{TM}}$ Q Exactive$^{\mathrm{TM}}$ UHMR). This combination
significantly reduces the overhead associated with custom instruments, while
benefiting from established high performance and reliability. We present
current performance characteristics including beam intensity, landing-energy
control, and deposition spot size for a broad range of molecules. In
combination with atomic force microscopy (AFM) and transmission electron
microscopy (TEM), we distinguish near-native from unfolded proteins and show
retention of native shape of protein assemblies after dehydration and
deposition. Further, we use an enzymatic assay to quantify activity of an
non-covalent protein complex after deposition an a dry surface. Together, these
results indicate a great potential of ES-IBD for applications in structural
biology, but also outline the challenges that need to be solved for it to reach
its full potential.


  We were interested to read the recent update on recommendations for reporting
of secondary findings in clinical sequencing1, and the accompanying updated
list of genes in which secondary findings should be sought (ACMG SF v3.0)2.
Though the authors discuss challenges around incomplete penetrance in
considerable detail, we are concerned that the recommendations do not fully
convey the degree of uncertainty regarding the penetrance of variants in genes
associated with inherited cardiomyopathies, which make up almost a quarter of
the list. Since penetrance is incomplete and age-related, individuals found to
carry variants will often require surveillance, rather than a one-off
definitive diagnostic assessment. There is a lack of evidence regarding
benefits, harms, and healthcare costs associated with opportunistic screening.


  Multi-task learning for molecular property prediction is becoming
increasingly important in drug discovery. However, in contrast to other
domains, the performance of multi-task learning in drug discovery is still not
satisfying as the number of labeled data for each task is too limited, which
calls for additional data to complement the data scarcity. In this paper, we
study multi-task learning for molecular property prediction in a novel setting,
where a relation graph between tasks is available. We first construct a dataset
(ChEMBL-STRING) including around 400 tasks as well as a task relation graph.
Then to better utilize such relation graph, we propose a method called SGNN-EBM
to systematically investigate the structured task modeling from two
perspectives. (1) In the \emph{latent} space, we model the task representations
by applying a state graph neural network (SGNN) on the relation graph. (2) In
the \emph{output} space, we employ structured prediction with the energy-based
model (EBM), which can be efficiently trained through noise-contrastive
estimation (NCE) approach. Empirical results justify the effectiveness of
SGNN-EBM. Code is available on https://github.com/chao1224/SGNN-EBM.


  Graph neural networks have become a powerful framework for learning complex
structure-property relationships and fast screening of chemical compounds.
Recently proposed methods have demonstrated that using 3D geometry information
of the molecule along with the bonding structure can lead to more accurate
prediction on a wide range of properties. A common practice is to use 3D
geometries computed through density functional theory (DFT) for both training
and testing of models. However, the computational time needed for DFT
calculations can be prohibitively large. Moreover, many of the properties that
we aim to predict can often be obtained with little or no overhead on top of
the DFT calculations used to produce the 3D geometry information, voiding the
need for a predictive model. To be practically useful for high-throughput
chemical screening and drug discovery, it is desirable to work with 3D
geometries obtained using less-accurate but much more efficient non-DFT
methods. In this work we investigate the impact of using non-DFT conformations
in the training and the testing of existing models and propose a data
augmentation method for improving the prediction accuracy of classical
forcefield-derived geometries.


  AlphaFold is a neural-network-based tool for the prediction of 3D structures
of protein. In CASP14, a blind structure prediction challenge, it performed
significantly better than other competitors, which makes it the best available
structure prediction tool. One of the outputs of AlphaFold is the probability
profile of residue-residue distances. This makes it possible to score any
conformation of the studied protein to express its compliance with the
AlphaFold model. Here we show how this score can be used to drive protein
folding simulation by metadynamics and parallel tempering metadynamics. By
parallel tempering metadynamics, we simulated folding of a mini-protein
Trp-cage beta hairpin and predicted their folding equilibria. We see the
potential of AlphaFold-based collective variable in applications beyond
structure prediction, such as in structure refinement or prediction of the
outcome of a mutation.


  Hybrid neural network models combine the advantages of a neural network's
fitting functionality with differential equation models to reflect actual
physical processes and are widely used in analyzing time-series data. Most
related studies have focused on linear hybrid models, but only a few have
examined nonlinear problems. In this work, we use a hybrid nonlinear epidemic
neural network as the entry point to study its power in predicting the correct
infection function of an epidemic model. To achieve this goal, we combine the
bifurcation theory of the nonlinear differential model with the mean-squared
error loss and design a novel loss function to ensure model trainability.
Furthermore, we find the unique existence conditions supporting ordinary
differential equations to estimate the correct infection function. Using the
Runge Kutta method, we perform numerical experiments on our proposed model and
verify its soundness. We also apply it to real COVID-19 data to accurately
discover the change law of its infectivity.


  Crowdsourcing is a mechanism by means of which groups of people are able to
execute a task by sharing ideas, efforts and resources. Thanks to the online
technologies, crowdsourcing has become in the last decade an even more utilized
process in different and diverse fields. An instance of such process is the
so-called "label aggregation problem": in practice, it is the evaluation of an
item by groups of agents, where each agent gives its own judgment of it.
Starting from the individual evaluations of their members, how can the groups
give their global assessment? In this work, by means of a game-theoretical,
evolutionary approach, we show that in most cases the majority rule (the group
evaluation is the evaluation of the majority of its members) is still the best
way to get a reliable group evaluation, even when the agents are not the best
experts of the topic at stake; on the other hand, we also show that noise
(i.e., fortuitous errors, misunderstanding, or every possible source of
non-deterministic outcomes) can undermine the efficiency of the procedure in
non-trivial situations. Therefore, in order to make the process as reliable as
possible, the presence of noise and its effects should be carefully taken into
account.


  Background: The increasing availability of databases containing both magnetic
resonance imaging (MRI) and genetic data allows researchers to utilize
multimodal data to better understand the characteristics of dementia of
Alzheimer's type (DAT). Objective: The goal of this study was to develop and
analyze novel biomarkers that can help predict the development and progression
of DAT. Methods: We used feature selection and ensemble learning classifier to
develop an image/genotype-based DAT score that represents a subject's
likelihood of developing DAT in the future. Three feature types were used: MRI
only, genetic only, and combined multimodal data. We used a novel data
stratification method to better represent different stages of DAT. Using a
pre-defined 0.5 threshold on DAT scores, we predicted whether or not a subject
would develop DAT in the future. Results: Our results on Alzheimer's Disease
Neuroimaging Initiative (ADNI) database showed that dementia scores using
genetic data could better predict future DAT progression for currently normal
control subjects (Accuracy=0.857) compared to MRI (Accuracy=0.143), while MRI
can better characterize subjects with stable mild cognitive impairment
(Accuracy=0.614) compared to genetics (Accuracy=0.356). Combining MRI and
genetic data showed improved classification performance in the remaining
stratified groups. Conclusion: MRI and genetic data can contribute to DAT
prediction in different ways. MRI data reflects anatomical changes in the
brain, while genetic data can detect the risk of DAT progression prior to the
symptomatic onset. Combining information from multimodal data in the right way
can improve prediction performance.


  Network structure or connectivity pattern is critical in determining
collective dynamics among interacting species in ecosystems. Conventional
research on species persistence in spatial populations has focused on static
network structure, though most real network structures change in time, forming
time-varying networks. This raises the question, in metacommunities, how does
the pattern of synchrony vary with temporal evolution in the network structure.
The synchronous dynamics among species are known to reduce metacommunity
persistence. Here, we consider a time-varying metacommunity small-world network
consisting of a chaotic three-species food chain oscillator in each patch/node.
The rate of change in the network connectivity is determined by the natural
frequency or its subharmonics of the constituent oscillator to allow sufficient
time for the evolution of species in between successive rewirings. We find that
over a range of coupling strengths and rewiring periods, even higher rewiring
probabilities drive a network from asynchrony towards synchrony. Moreover, in
networks with a small rewiring period, an increase in average degree (more
connected networks) pushes the asynchronous dynamics to synchrony. On the
contrary, in networks with a low average degree, a higher rewiring period
drives the synchronous dynamics to asynchrony resulting in increased species
persistence. Our results also follow the calculation of synchronization time
and robust across other ecosystem models. Overall, our study opens the
possibility of developing temporal connectivity strategies to increase species
persistence in ecological networks.


  Acylphosphatase (AcP) is a small protein with 98 amino acid residues that
catalyzes the hydrolysis of carboxyl-phosphate bonds. AcP is a typical
two-state protein with slow folding rate due to its relatively large contact
order in the native structure. The mechanical properties and unfolding behavior
of AcP has been studied by atomic force microscope. But the folding and
unfolding dynamics at low forces has not been reported. Here using stable
magnetic tweezers, we measured the force-dependent folding rates within a force
range from 1 pN to 3 pN, and unfolding rates from 15 pN to 40 pN. The obtained
unfolding rates show different force sensitivities at forces below and above
~27 pN, which determines a free energy landscape with two energy barriers. Our
results indicate that the free energy landscape of small globule proteins have
general Bactrian camel shape, and large contact order of the native state
produces a high barrier dominate at low forces.


  Cities become mission-critical zones during pandemics and it is vital to
develop a better understanding of the factors that are associated with
infection levels. The COVID-19 pandemic has impacted many cities severely;
however, there is significant variance in its impact across cities. Pandemic
infection levels are associated with inherent features of cities (e.g.,
population size, density, mobility patterns, socioeconomic condition, and
health environment), which need to be better understood. Intuitively, the
infection levels are expected to be higher in big urban agglomerations, but the
measurable influence of a specific urban feature is unclear. The present study
examines 41 variables and their potential influence on COVID-19 cases and
fatalities. The study uses a multi-method approach to study the influence of
variables, classified as demographic, socioeconomic, mobility and connectivity,
urban form and density, and health and environment dimensions. This study
develops an index dubbed the PVI-CI for classifying the pandemic vulnerability
levels of cities, grouping them into five vulnerability classes, from very high
to very low. Furthermore, clustering and outlier analysis provides insights on
the spatial clustering of cities with high and low vulnerability scores. This
study provides strategic insights into levels of influence of key variables
upon the spread of infections as well as fatalities, along with an objective
ranking for the vulnerability of cities. Thus it provides critical wisdom
needed for urban healthcare policy and resource management. The pandemic
vulnerability index calculation method and the process present a blueprint for
the development of similar indices for cities in other countries, leading to a
better understanding and improved pandemic management for urban areas and
post-pandemic urban planning across the world.


  In this paper, we present a method to forecast the spread of SARS-CoV-2
across regions with a focus on the role of mobility. Mobility has previously
been shown to play a significant role in the spread of the virus, particularly
between regions. Here, we investigate under which epidemiological circumstances
incorporating mobility into transmission models yields improvements in the
accuracy of forecasting, where we take the situation in the Netherlands during
and after the first wave of transmission in 2020 as a case study. We assess the
quality of forecasting on the detailed level of municipalities, instead of on a
nation-wide level. To model transmissions, we use a simple mobility-enhanced
SEIR compartmental model with subpopulations corresponding to the Dutch
municipalities. We use commuter information to quantify mobility, and develop a
method based on maximum likelihood estimation to determine the other relevant
parameters. We show that taking inter-regional mobility into account generally
leads to an improvement in forecast quality. However, at times when policies
are in place that aim to reduce contacts or travel, this improvement is very
small.


  The dynamics of microbial communities is incredibly complex, determined by
competition for metabolic substrates and cross-feeding of byproducts. Species
in the community grow by harvesting energy from chemical reactions that
transform substrates to products. In many anoxic environments, these reactions
are close to thermodynamic equilibrium and growth is slow. To understand the
community structure in these energy-limited environments, we developed a
microbial community consumer-resource model incorporating energetic and
thermodynamic constraints on an interconnected metabolic network. The central
ingredient of the model is product inhibition, meaning that microbial growth
may be limited not only by depletion of metabolic substrates but also by
accumulation of products. We demonstrate that these additional constraints on
microbial growth cause a convergence in the structure and function of the
community metabolic network -- independent of species composition and
biochemical details -- providing a possible explanation for convergence of
community function despite taxonomic variation observed in many natural and
industrial environments. Furthermore, we discovered that the structure of
community metabolic network is governed by the thermodynamic principle of
maximum heat dissipation. Overall, the work demonstrates how universal
thermodynamic principles may constrain community metabolism and explain
observed functional convergence in microbial communities.


  String edit distances have been used for decades in applications ranging from
spelling correction and web search suggestions to DNA analysis. Most string
edit distances are variations of the Levenshtein distance and consider only
single-character edits. In forensic applications polymorphic genetic markers
such as short tandem repeats (STRs) are used. At these repetitive motifs the
DNA copying errors consist of more than just single base differences. More
often the phenomenon of ``stutter'' is observed, where the number of repeated
units differs (by whole units) from the template. To adapt the Levenshtein
distance to be suitable for forensic applications where DNA sequence similarity
is of interest, a generalized string edit distance is defined that accommodates
the addition or deletion of whole motifs in addition to single-nucleotide
edits. A dynamic programming implementation is developed for computing this
distance between sequences. The novelty of this algorithm is in handling the
complex interactions that arise between multiple- and single-character edits.
Forensic examples illustrate the purpose and use of the Restricted Forensic
Levenshtein (RFL) distance measure, but applications extend to sequence
alignment and string similarity in other biological areas, as well as dynamic
programming algorithms more broadly.


  Numerous studies have investigated the differences and similarities between
protein structures determined by solution NMR spectroscopy and those determined
by x-ray crystallography. A fundamental question is whether any observed
differences are due to differing methodologies, or to differences in the
behavior of proteins in solution versus in the crystalline state. Here, we
compare the properties of the hydrophobic cores of high-resolution protein
crystal structures and those in NMR structures, determined using increasing
numbers and types of restraints. Prior studies have reported that many NMR
structures have denser cores compared to those of high-resolution x-ray crystal
structures. Our current work investigates this result in more detail, and finds
that these NMR structures tend to violate basic features of protein
stereochemistry, such as small non-bonded atomic overlaps and few Ramachandran
and side chain dihedral angle outliers. We find that NMR structures solved with
more restraints, and which do not significantly violate stereochemistry, have
hydrophobic cores that have a similar size and packing fraction as their
counterparts determined by x-ray crystallography at high-resolution. These
results lead us to conclude that, at least regarding the core packing
properties, high-quality structures determined by NMR and x-ray crystallography
are the same, and the differences reported earlier are most likely a
consequence of methodology, rather than fundamental differences between the
protein in the two different environments.


  Organisms may respond to local stimuli that benefit or threaten their
fitness. The adaptive movement behaviour may allow individuals to adjust their
speed to maximise the chances of being in comfort zones, where death risk is
minimal. We investigate spatial cyclic models where the rock-paper-scissors
game rules describe the nonhierarchical dominance. We assume that organisms of
one out of the species can control the mobility rate in response to the
information obtained from scanning the environment. Running a series of
stochastic simulations, we quantify the effects of the movement strategy on the
spatial patterns and population dynamics. Our findings show that the ability to
change mobility to adapt to environmental clues is not reflected in an
advantage in cyclic spatial games. The adaptive movement provokes a delay in
the spatial domains occupied by the species in the spiral waves, making the
group more vulnerable to the advance of the dominant species and less efficient
in taking territory from the dominated species. Our outcomes also show that the
effects of adaptive movement behaviour accentuate whether most individuals have
a long-range neighbourhood perception. Our results may be helpful for
biologists and data scientists to comprehend the dynamics of ecosystems where
adaptive processes are fundamental.


  Network science has evolved into an indispensable platform for studying
complex systems. But recent research has identified limits of classical
networks, where links connect pairs of nodes, to comprehensively describe group
interactions. Higher-order networks, where a link can connect more than two
nodes, have therefore emerged as a new frontier in network science. Since group
interactions are common in social, biological, and technological systems,
higher-order networks have recently led to important new discoveries across
many fields of research. We here review these works, focusing in particular on
the novel aspects of the dynamics that emerges on higher-order networks. We
cover a variety of dynamical processes that have thus far been studied,
including different synchronization phenomena, contagion processes, the
evolution of cooperation, and consensus formation. We also outline open
challenges and promising directions for future research.


  Graphs are ubiquitous in encoding relational information of real-world
objects in many domains. Graph generation, whose purpose is to generate new
graphs from a distribution similar to the observed graphs, has received
increasing attention thanks to the recent advances of deep learning models. In
this paper, we conduct a comprehensive review on the existing literature of
deep graph generation from a variety of emerging methods to its wide
application areas. Specifically, we first formulate the problem of deep graph
generation and discuss its difference with several related graph learning
tasks. Secondly, we divide the state-of-the-art methods into three categories
based on model architectures and summarize their generation strategies.
Thirdly, we introduce three key application areas of deep graph generation.
Lastly, we highlight challenges and opportunities in the future study of deep
graph generation. We hope that our survey will be useful for researchers and
practitioners who are interested in this exciting and rapidly-developing field.


  Computational models have great potential to accelerate bioscience,
bioengineering, and medicine. However, it remains challenging to reproduce and
reuse simulations, in part, because the numerous formats and methods for
simulating various subsystems and scales remain siloed by different software
tools. For example, each tool must be executed through a distinct interface. To
help investigators find and use simulation tools, we developed BioSimulators
(https://biosimulators.org), a central registry of the capabilities of
simulation tools and consistent Python, command-line, and containerized
interfaces to each version of each tool. The foundation of BioSimulators is
standards, such as CellML, SBML, SED-ML, and the COMBINE archive format, and
validation tools for simulation projects and simulation tools that ensure these
standards are used consistently. To help modelers find tools for particular
projects, we have also used the registry to develop recommendation services. We
anticipate that BioSimulators will help modelers exchange, reproduce, and
combine simulations.


  Background. During 2021, the COVID-19 pandemic was characterized by the
emergence of lineages with increased fitness. For most of these variants,
quantitative information is scarce on epidemiological quantities such as the
incubation period and generation time, which are critical for both public
health decisions and scientific research. Method. We analyzed a dataset
collected during contact tracing activities in the province of Reggio Emilia,
Italy, throughout 2021. We determined the distributions of the incubation
period using information on negative PCR tests and the date of last exposure
from 282 symptomatic cases. We estimated the distributions of the intrinsic
generation time (the time between the infection dates of an infector and its
secondary cases under a fully susceptible population) using a Bayesian
inference approach applied to 4,435 SARS-CoV-2 cases clustered in 1,430
households where at least one secondary case was recorded. Results. We
estimated a mean incubation period of 4.9 days (95% credible intervals, CrI,
4.4-5.4; 95 percentile of the mean distribution: 1-12) for Alpha and 4.5 days
(95%CrI 4.0-5.0; 95 percentile: 1-10) for Delta. The intrinsic generation time
was estimated to have a mean of 6.0 days (95% CrI 5.6-6.4; 95 percentile: 1-15)
for Alpha and of 6.6 days (95%CrI 6.0-7.3; 95 percentile: 1-18) for Delta. The
household serial interval was 2.6 days (95%CrI 2.4-2.7) for Alpha and 2.4 days
(95%CrI 2.2-2.6) for Delta, and the estimated proportion of pre-symptomatic
transmission was 54-55% for both variants. Conclusions. These results indicate
limited differences in the incubation period and intrinsic generation time of
SARS-CoV-2 variants Alpha and Delta compared to ancestral lineages.


  In the brain, both neurons and glial cells work in conjunction with each
other during information processing. Stimulation of neurons can cause calcium
oscillations in astrocytes which in turn can affect neuronal calcium dynamics.
The "glissandi" effect is one such phenomenon, associated with a decrease in
infraslow fluctuations, in which synchronized calcium oscillations propagate as
a wave in hundreds of astrocytes. Nitric oxide molecules released from the
astrocytes contribute to synaptic functions on the basis of the underlying
astrocyte-neuron interaction network. In this study, by defining an
astrocyte-neuronal (A-N) unit as an integrated circuit of one neuron and one
astrocyte, we developed a minimal model of neuronal stimulus-dependent and
nitric oxide-mediated emergence of calcium waves in astrocytes. Incorporating
inter-unit communication via nitric oxide molecules, a coupled network of 1,000
such A-N units is developed in which multiple stable regimes were found to
emerge in astrocytes. We examined the ranges of neuronal stimulus strength and
the coupling strength between A-N units that give rise to such dynamical
behaviors. We also report that there exists a range of coupling strength,
wherein units not receiving stimulus also start showing oscillations and become
synchronized. Our results support the hypothesis that glissandi-like phenomena
exhibiting synchronized calcium oscillations in astrocytes help in efficient
synaptic transmission by reducing the energy demand of the process.


  Amsterdam and other UNAIDS Fast-Track cities aim for zero new HIV infections.
Utilising molecular and clinical data of the ATHENA observational HIV cohort,
our primary aims are to estimate the proportion of undiagnosed HIV infections
and the proportion of locally acquired infections in Amsterdam in 2014-2018,
both in MSM and heterosexuals and Dutch-born and foreign-born individuals.
  We located diagnosed HIV infections in Amsterdam using postcode data at time
of registration to the cohort, and estimated their date of infection using
clinical HIV data. We then inferred the proportion undiagnosed from the
estimated times to diagnosis. To determine sources of Amsterdam infections, we
used HIV sequences of people living with HIV (PLHIV) within a background of
other Dutch and international sequences to phylogenetically reconstruct
transmission chains. Frequent late diagnoses indicate that more recent
phylogenetically observed chains are increasingly incomplete, and we use a
Bayesian model to estimate the actual growth of Amsterdam transmission chains,
and the proportion of locally acquired infections.
  We estimate that 20% [95% CrI 18-22%] of infections acquired among MSM
between 2014-2018 were undiagnosed by the start of 2019, and 44% [37-50%] among
heterosexuals, with variation by place of birth. The estimated proportion of
MSM infections in 2014-2018 that were locally acquired was 68% [61-74%], with
no substantial differences by region of birth. In heterosexuals, this was 57%
[41-71%] overall, with heterogeneity by place of birth.
  The data indicate substantial potential to further curb local transmission,
in both MSM and heterosexual Amsterdam residents. In 2014-2018 the largest
proportion of local transmissions in Amsterdam are estimated to have occurred
in foreign-born MSM, who would likely benefit most from intensified
interventions.


  The problem of molecular generation has received significant attention
recently. Existing methods are typically based on deep neural networks and
require training on large datasets with tens of thousands of samples. In
practice, however, the size of class-specific chemical datasets is usually
limited (e.g., dozens of samples) due to labor-intensive experimentation and
data collection. This presents a considerable challenge for the deep learning
generative models to comprehensively describe the molecular design space.
Another major challenge is to generate only physically synthesizable molecules.
This is a non-trivial task for neural network-based generative models since the
relevant chemical knowledge can only be extracted and generalized from the
limited training data. In this work, we propose a data-efficient generative
model that can be learned from datasets with orders of magnitude smaller sizes
than common benchmarks. At the heart of this method is a learnable graph
grammar that generates molecules from a sequence of production rules. Without
any human assistance, these production rules are automatically constructed from
training data. Furthermore, additional chemical knowledge can be incorporated
in the model by further grammar optimization. Our learned graph grammar yields
state-of-the-art results on generating high-quality molecules for three monomer
datasets that contain only ${\sim}20$ samples each. Our approach also achieves
remarkable performance in a challenging polymer generation task with only $117$
training samples and is competitive against existing methods using $81$k data
points. Code is available at https://github.com/gmh14/data_efficient_grammar.


  Cryo-electron microscopy (cryo-EM) has become a tool of fundamental
importance in structural biology, helping us understand the basic building
blocks of life. The algorithmic challenge of cryo-EM is to jointly estimate the
unknown 3D poses and the 3D electron scattering potential of a biomolecule from
millions of extremely noisy 2D images. Existing reconstruction algorithms,
however, cannot easily keep pace with the rapidly growing size of cryo-EM
datasets due to their high computational and memory cost. We introduce cryoAI,
an ab initio reconstruction algorithm for homogeneous conformations that uses
direct gradient-based optimization of particle poses and the electron
scattering potential from single-particle cryo-EM data. CryoAI combines a
learned encoder that predicts the poses of each particle image with a
physics-based decoder to aggregate each particle image into an implicit
representation of the scattering potential volume. This volume is stored in the
Fourier domain for computational efficiency and leverages a modern coordinate
network architecture for memory efficiency. Combined with a symmetrized loss
function, this framework achieves results of a quality on par with
state-of-the-art cryo-EM solvers for both simulated and experimental data, one
order of magnitude faster for large datasets and with significantly lower
memory requirements than existing methods.


  The principles that govern the organization of genomes, which are needed for
a deeper understanding of how chromosomes are packaged and function in
eukaryotic cells, could be deciphered if the three-dimensional (3D) structures
are known. Recently, single-cell imaging experiments have determined the 3D
coordinates of a number of loci in a chromosome. Here, we introduce a
computational method (Distance Matrix to Ensemble of Structures, DIMES), based
on the maximum entropy principle, with experimental pair-wise distances between
loci as constraints, to generate a unique ensemble of 3D chromatin structures.
Using the ensemble of structures, we quantitatively account for the
distribution of pair-wise distances, three-body co-localization and
higher-order interactions. We demonstrate that the DIMES method can be applied
to both small length-scale and chromosome-scale imaging data to quantify the
extent of heterogeneity and fluctuations in the shapes on various length
scales. We develop a perturbation method that is used in conjunction with DIMES
to predict the changes in 3D structures from structural variations. Our method
also reveals quantitative differences between the 3D structures inferred from
Hi-C and the ones measured in imaging experiments. Finally, the physical
interpretation of the parameters extracted from DIMES provides insights into
the origin of phase separation between euchromatin and heterochromatin domains.


  As the outbreak of COVID-19 enters its third year, we have now enough data to
analyse the behavior of the pandemic with mathematical models over a long
period of time. The pandemic alternates periods of high and low infections, in
a way that sheds a light on the nature of mathematical model that can be used
for reliable predictions. The main hypothesis of the model presented here is
that the oscillatory behavior is a structural feature of the outbreak, even
without postulating a time-dependence of the coefficients. As such, it should
be reflected by the presence of limit cycles as asymptotic solutions. This
stems from the introduction of (i) a non-linear waning immunity based on the
concept of immunity booster (already used for other pathologies); (ii) a fine
description of the compartments with a discrimination between individuals
infected/vaccinated for the first time, and individuals already
infected/vaccinated, undergoing to new infections/doses. We provide a
proof-of-concept that our novel model is capable of reproducing long-term
oscillatory behavior of many infectious diseases, and, in particular, the
periodic nature of the waves of infection. Periodic solutions are inherent to
the model, and achieved without changing parameter values in time. This may
represent an important step in the long-term modeling of COVID-19 and similar
diseases, as the natural, unforced behavior of the solution shows the
qualitative characteristics observed during the COVID-19 pandemic.


  Increased ability to predict protein structures is moving research focus
towards understanding protein dynamics. A promising approach is to represent
protein dynamics through networks and take advantage of well-developed methods
from network science. Most studies build protein dynamics networks from
correlation measures, an approach that only works under very specific
conditions, instead of the more robust inverse approach. Thus, we apply the
inverse approach to the dynamics of protein dihedral angles, a system of
internal coordinates, to avoid structural alignment. Using the
well-characterized adhesion protein, FimH, we show that our method identifies
networks that are physically interpretable, robust, and relevant to the
allosteric pathway sites. We further use our approach to detect dynamical
differences, despite structural similarity, for Siglec-8 in the immune system,
and the SARS-CoV-2 spike protein. Our study demonstrates that using the inverse
approach to extract a network from protein dynamics yields important
biophysical insights.


  During the Covid-19 pandemic a key role is played by vaccination to combat
the virus. There are many possible policies for prioritizing vaccines, and
different criteria for optimization: minimize death, time to herd immunity,
functioning of the health system. Using an age-structured population
compartmental finite-dimensional optimal control model, our results suggest
that the eldest to youngest vaccination policy is optimal to minimize deaths.
Our model includes the possible infection of vaccinated populations. We apply
our model to real-life data from the US Census for New Jersey and Florida,
which have a significantly different population structure. We also provide
various estimates of the number of lives saved by optimizing the vaccine
schedule and compared to no vaccination.


  Spatial transcriptomics (ST) has advanced significantly in the last few
years. Such advancement comes with the urgent need for novel computational
methods to handle the unique challenges of ST data analysis. Many artificial
intelligence (AI) methods have been developed to utilize various machine
learning and deep learning techniques for computational ST analysis. This
review provides a comprehensive and up-to-date survey of current AI methods for
ST analysis.


  Neural codes appear efficient. Naturally, neuroscientists contend that an
efficient process is responsible for generating efficient codes. They argue
that natural selection is the efficient process that generates those codes.
Although natural selection is an adaptive process, evolution itself, is not.
Evolution consists of not only natural selection, but also neutral stochastic
forces that can generate biological inefficiencies. The explanatory power of
natural selection cannot be appealed to, without regards for the remaining
evolutionary forces. In this paper, we aim to reformulate the explanatory role
of evolutionary forces on neural coding, with special attention to neutral
forces. We propose a framework that argues for differing contributions of
adaptive and stochastic evolutionary forces, for different phenotypic `levels',
including those of neural codes. We assert that this framework is of special
interest to neuroscience, because the field has derived much progress from an
efficiency-based worldview. We advocate for a pluralistic neuroscience capable
of appealing to both adaptive and non-adaptive explanations.


  Population decline is projected to become widespread in Europe, with the
continental population set to reverse its longstanding trajectory of growth
within the next five years. This represents unfamiliar demographic territory.
Despite this, literature on decline remains sparse and our understanding
porous. Particular epistemological deficiencies stem from a lack of both
cross-national and temporal analyses of population decline. This study seeks to
address these gapsthrough the novel application of sequence and cluster
analysis techniques to examine variations in population decline trajectories
since 2000 in 696 sub-national areas across 33 European territories. The
methodology allows for a holistic understanding of decline trajectories
capturing differences in the ordering, timing, magnitude and spatial structure
of population decline. We identify a typology of population decline
distinguishing seven distinct pathways to depopulation and chart their
geographies. Results revealed differentiated pathways of depopulation in
continental sub-regions, with consistent and rapid declines in the east,
persistent but moderate declines in central Europe, accelerating declines in
the south and decelerating population declines in the west. Results also
revealed differentiated patterns of depopulation across the rural-urban
continuum, with urban and populous areas experiencing deceleration in
population decline, while population decline accelerates or stabilises in rural
areas. Small and mid-sized areas displayed heterogeneous depopulation
trajectories, highlighting the importance of local contextual factors in
influencing trajectories of population decline.


  In this paper we propose a new mathematical model capable of merging
Darwinian Evolution, Human History and SETI into a single mathematical scheme:
1) Darwinian Evolution over the last 3.5 billion years is defined as one
particular realization of a certain stochastic process called Geometric
Brownian Motion (GBM). This GBM yields the fluctuations in time of the number
of species living on Earth. Its mean value curve is an increasing exponential
curve, i.e. the exponential growth of Evolution. 2) In 2008 this author
provided the statistical generalization of the Drake equation yielding the
number N of communicating ET civilizations in the Galaxy. N was shown to follow
the lognormal probability distribution. 3) We call "b-lognormals" those
lognormals starting at any positive time b ("birth") larger than zero. Then the
exponential growth curve becomes the geometric locus of the peaks of a
one-parameter family of b-lognormals: this is our way to re-define Cladistics.
4) b-lognormals may be also be interpreted as the lifespan of any living being
(a cell, or an animal, a plant, a human, or even the historic lifetime of any
civilization). Applying this new mathematical apparatus to Human History, leads
to the discovery of the exponential progress between Ancient Greece and the
current USA as the envelope of all b-lognormals of Western Civilizations over a
period of 2500 years. 5) We then invoke Shannon's Information Theory. The
b-lognormals' entropy turns out to be the index of "development level" reached
by each historic civilization. We thus get a numerical estimate of the entropy
difference between any two civilizations, like the Aztec-Spaniard difference in
1519. 6) In conclusion, we have derived a mathematical scheme capable of
estimating how much more advanced than Humans an Alien Civilization will be
when the SETI scientists will detect the first hints about ETs.


  A rooted phylogenetic network is a directed acyclic graph with a single root,
whose sinks correspond to a set of species. As such networks are useful for
representing the evolution of species that have undergone reticulate evolution,
there has been great interest in developing the theory behind and algorithms
for constructing them. However, unlike evolutionary trees, these networks can
be highly non-planar, which can make them difficult to visualise and interpret.
Here we investigate properties of planar rooted phylogenetic networks and
algorithms for deciding whether or not rooted networks have certain special
planarity properties. In particular, we introduce three natural subclasses of
planar rooted phylogenetic networks and show that they form a hierarchy. In
addition, for the well-known level-k networks, we show that level-1, -2, -3
networks are always outer, terminal, and upward planar, respectively, and that
level-4 networks are not necessarily planar. Finally, we show that a regular
network is terminal planar if and only if it is pyramidal. Our results make use
of the highly developed field of planar digraphs, and we believe that the link
between phylogenetic networks and planar graphs should prove useful in future
for developing new approaches to both construct and visualise phylogenetic
networks.


  The discovery of new exoplanets makes us wonder where each new exoplanet
stands along its way to develop life as we know it on Earth. Our Evo-SETI
Theory is a mathematical way to face this problem. We describe cladistics and
evolution by virtue of a few statistical equations based on lognormal
probability density functions (pdf) in the time. We call b-lognormal a
lognormal pdf starting at instant b (birth). Then, the lifetime of any living
being becomes a suitable b-lognormal in the time. Next, our "Peak-Locus
Theorem" translates cladistics: each species created by evolution is a
b-lognormal whose peak lies on the exponentially growing number of living
species. This exponential is the mean value of a stochastic process called
"Geometric Brownian Motion" (GBM). Past mass extinctions were all-lows of this
GBM. In addition, the Shannon Entropy (with a reversed sign) of each
b-lognormal is the measure of how evolved that species is, and we call it
EvoEntropy. The "molecular clock" is re-interpreted as the EvoEntropy straight
line in the time whenever the mean value is exactly the GBM exponential. We
were also able to extend the Peak-Locus Theorem to any mean value other than
the exponential. For example, we derive in this paper for the first time the
EvoEntropy corresponding to the Markov-Korotayev (2007) "cubic" evolution: a
curve of logarithmic increase.


  Thermodynamic selection is an indirect competition between agents feeding on
the same energy resource and obeying the laws of thermodynamics. We examine
scenarios of this selection, where the agent is modeled as a heat-engine
coupled to two thermal baths and extracting work from the high-temperature
bath. The agents can apply different work-extracting, game-theoretical
strategies, e.g. the maximum power or the maximum efficiency. They can also
have a fixed structure or be adaptive. Depending on whether the resource (i.e.
the high-temperature bath) is infinite or finite, the fitness of the agent
relates to the work-power or the total extracted work. These two selection
scenarios lead to increasing or decreasing efficiencies of the work-extraction,
respectively. The scenarios are illustrated via plant competition for sunlight,
and the competition between different ATP production pathways. We also show
that certain general concepts of game-theory and ecology--the prisoner's
dilemma and the maximal power principle--emerge from the thermodynamics of
competing agents. We emphasize the role of adaptation in developing efficient
work-extraction mechanisms.


  We study a fundamental problem in structure-based drug design -- generating
molecules that bind to specific protein binding sites. While we have witnessed
the great success of deep generative models in drug design, the existing
methods are mostly string-based or graph-based. They are limited by the lack of
spatial information and thus unable to be applied to structure-based design
tasks. Particularly, such models have no or little knowledge of how molecules
interact with their target proteins exactly in 3D space. In this paper, we
propose a 3D generative model that generates molecules given a designated 3D
protein binding site. Specifically, given a binding site as the 3D context, our
model estimates the probability density of atom's occurrences in 3D space --
positions that are more likely to have atoms will be assigned higher
probability. To generate 3D molecules, we propose an auto-regressive sampling
scheme -- atoms are sampled sequentially from the learned distribution until
there is no room for new atoms. Combined with this sampling scheme, our model
can generate valid and diverse molecules, which could be applicable to various
structure-based molecular design tasks such as molecule sampling and linker
design. Experimental results demonstrate that molecules sampled from our model
exhibit high binding affinity to specific targets and good drug properties such
as drug-likeness even if the model is not explicitly optimized for them.


  We recently provided a recursive construction of the solution of the
selection-recombination equation under the assumption of single-crossover and
with selection acting on a single site. This recursion is intimately connected
to various dual processes that describe the ancestry of a random sample,
backward in time. We now aim to unify and generalise this approach by
introducing a new labelled partitioning process with general Markovian labels.
As a concrete example, we consider the selection-mutation-recombination
equation.


  Gene expression is a biochemical process, where stochastic binding and
un-binding events naturally generate fluctuations and cell-to-cell variability
in gene dynamics. These fluctuations typically have destructive consequences
for proper biological dynamics and function (e.g., loss of timing and synchrony
in biological oscillators). Here, we show that gene expression noise
counter-intuitively accelerates the evolution of a biological oscillator and,
thus, can impart a benefit to living organisms. We used computer simulations to
evolve two mechanistic models of a biological oscillator at different levels of
gene expression noise. We first show that gene expression noise induces
oscillatory-like dynamics in regions of parameter space that cannot oscillate
in the absence of noise. We then demonstrate that these noise-induced
oscillations generate a fitness landscape whose gradient robustly and quickly
guides evolution by mutation towards robust and self-sustaining oscillation.
These results suggest that noise can help dynamical systems evolve or learn new
behavior by revealing cryptic dynamic phenotypes outside the bifurcation point.


  Antigen-specific T cells play an essential role in immunoregulation and
diseases such as cancer. Characterizing the T cell receptor (TCR) sequences
that encode T cell specificity is critical for elucidating the antigenic
determinants of immunological diseases and designing therapeutic remedies.
However, methods of obtaining single-cell TCR sequencing data are labor and
cost intensive, requiring cell sorting and full length single-cell
RNA-sequencing (scRNA-seq). New high-throughput 3' cell-barcoding scRNA-seq
methods can simplify and scale this process; however, they do not routinely
capture TCR sequences during library preparation and sequencing. While 5'
cell-barcoding scRNA-seq methods can be used to examine TCR repertoire at
single-cell resolution, it requires specialized reagents which cannot be
applied to samples previously processed using 3' cell-barcoding methods. Here,
we outline a method for sequencing TCR$\alpha$ and TCR$\beta$ transcripts from
samples already processed using 3' cell-barcoding scRNA-seq platforms, ensuring
TCR recovery at a single-cell resolution. In short, a fraction of the 3'
barcoded whole transcriptome amplification (WTA) product typically used to
generate a massively parallel 3' scRNA-seq library is enriched for TCR
transcripts using biotinylated probes, and further amplified using the same
universal primer sequence from WTA. Primer extension using TCR V-region primers
and targeted PCR amplification results in a 3' barcoded single-cell
CDR3-enriched library that can be sequenced with custom sequencing primers.
Coupled with 3' scRNA-seq of the same WTA, this method enables simultaneous
analysis of single-cell transcriptomes and TCR sequences which can help
interpret inherent heterogeneity among antigen-specific T cells and salient
disease biology. This method can be adapted to enrich other transcripts of
interest from 3' and 5' barcoded WTA libraries.


  Next generation sequencing of B cell receptor (BCR) repertoires has become a
ubiquitous tool for understanding the antibody-mediated immune response: it is
now common to have large volumes of sequence data coding for both the heavy and
light chain subunits of the BCR. However, until the recent development of high
throughput methods of preserving heavy/light chain pairing information, these
samples contained no explicit information on which heavy chain sequence pairs
with which light chain sequence. One of the first steps in analyzing such BCR
repertoire samples is grouping sequences into clonally related families, where
each stems from a single rearrangement event. Many methods of accomplishing
this have been developed, however, none so far has taken full advantage of the
newly-available pairing information. This information can dramatically improve
clustering performance, especially for the light chain. The light chain has
traditionally been challenging for clonal family inference because of its low
diversity and consequent abundance of non-clonal families with
indistinguishable naive rearrangements. Here we present a method of
incorporating this pairing information into the clustering process in order to
arrive at a more accurate partition of the data into clonally related families.
We also demonstrate two methods of fixing imperfect pairing information, which
may allow for simplified sample preparation and increased sequencing depth.
Finally, we describe several other improvements to the partis software package
(https://github.com/psathyrella/partis).


  Chemical reaction prediction, involving forward synthesis and retrosynthesis
prediction, is a fundamental problem in organic synthesis. A popular
computational paradigm formulates synthesis prediction as a
sequence-to-sequence translation problem, where the typical SMILES is adopted
for molecule representations. However, the general-purpose SMILES neglects the
characteristics of chemical reactions, where the molecular graph topology is
largely unaltered from reactants to products, resulting in the suboptimal
performance of SMILES if straightforwardly applied. In this article, we propose
the root-aligned SMILES (R-SMILES), which specifies a tightly aligned
one-to-one mapping between the product and the reactant SMILES for more
efficient synthesis prediction. Due to the strict one-to-one mapping and
reduced edit distance, the computational model is largely relieved from
learning the complex syntax and dedicated to learning the chemical knowledge
for reactions. We compare the proposed R-SMILES with various state-of-the-art
baselines and show that it significantly outperforms them all, demonstrating
the superiority of the proposed method.


  The identification of drug-target binding affinity (DTA) has attracted
increasing attention in the drug discovery process due to the more specific
interpretation than binary interaction prediction. Recently, numerous deep
learning-based computational methods have been proposed to predict the binding
affinities between drugs and targets benefiting from their satisfactory
performance. However, the previous works mainly focus on encoding biological
features and chemical structures of drugs and targets, with a lack of
exploiting the essential topological information from the drug-target affinity
network. In this paper, we propose a novel hierarchical graph representation
learning model for the drug-target binding affinity prediction, namely
HGRL-DTA. The main contribution of our model is to establish a hierarchical
graph learning architecture to incorporate the intrinsic properties of
drug/target molecules and the topological affinities of drug-target pairs. In
this architecture, we adopt a message broadcasting mechanism to integrate the
hierarchical representations learned from the global-level affinity graph and
the local-level molecular graph. Besides, we design a similarity-based
embedding map to solve the cold start problem of inferring representations for
unseen drugs and targets. Comprehensive experimental results under different
scenarios indicate that HGRL-DTA significantly outperforms the state-of-the-art
models and shows better model generalization among all the scenarios.


  A living cell uses long tubular appendages for locomotion and sensory
purposes. Hence, assembling and maintaining a protrusion of correct length is
crucial for its survival and overall performance. Usually the protrusions lack
the machinery for the synthesis of building blocks and imports them from the
cell body. What are the unique features of the transport logistics which
facilitate the exchange of these building blocks between the cell and the
protrusion? What kind of `rulers' and `timers' does the cell use for
constructing its appendages of correct length on time? How do the multiple
appendages coordinate and communicate among themselves during different stages
of their existence? How frequently do the fluctuations drive the length of
these dynamic protrusions out of the acceptable bounds? These questions are
addressed from a broad perspective in this review which is organized in three
parts. In part-I the list of all known cell protrusions is followed by a
comprehensive list of the mechanisms of length control of cell protrusions
reported in the literature. We review not only the dynamics of the genesis of
the protrusions, but also their resorption and regrowth as well as regeneration
after amputation. As a case study in part-II, the specific cell protrusion that
has been discussed in detail is eukaryotic flagellum (also known as cilium);
this choice was dictated by the fact that flagellar length control mechanisms
have been studied most extensively over more than half a century in cells with
two or more flagella. Although limited in scope, brief discussions on a few
non-flagellar cell protrusions in part-III of this review is intended to
provide a glimpse of the uncharted territories and challenging frontiers of
research on subcellular length control phenomena that awaits vigorous
investigations.


  Control of multistable dynamical system has important applications, from
physics to biology. Here, we attack this problem from the perspective of local
sensitivity analysis. We develop sensitivity rules to control properties of
continuous-time multistable dynamics with simple attractors, namely, the
relative size and depth of their basins of attraction. Our parameter control
signal is computationally cheap and provides counter-intuitive information
about the sensitive parameters to be manipulated in an experimental setting.


  We study the role of the adaptive movement strategy in promoting biodiversity
in cyclic models described by the rock-paper-scissors game rules. We assume
that individuals of one out of the species may adjust their movement to escape
hostile regions and stay longer in their comfort zones. Running a series of
stochastic simulations, we calculate the alterations in the spatial patterns
and population densities in scenarios where not all organisms are physically or
cognitively conditioned to perform the behavioural strategy. Although the
adaptive movement strategy is not profitable in terms of territorial dominance
for the species, it may promote biodiversity. Our findings show that if all
individuals are apt to move adaptively, coexistence probability increases for
intermediary mobility. The outcomes also show that even if not all individuals
can react to the signals received from the neighbourhood, biodiversity is still
benefited, but for a shorter mobility range. We find that the improvement in
the coexistence conditions is more accentuated if organisms adjust their
movement intensely and can receive sensory information from longer distances.
We also discover that biodiversity is slightly promoted for high mobility if
the proportion of individuals participating in the strategy is low. Our results
may be helpful for biologists and data scientists to understand adaptive
process learning in system biology.


  De novo peptide sequencing aims to recover amino acid sequences of a peptide
from tandem mass spectrometry (MS) data. Existing approaches for de novo
analysis enumerate MS evidence for all amino acid classes during inference. It
leads to over-trimming on receptive fields of MS data and restricts MS evidence
associated with following undecoded amino acids. Our approach, DPST,
circumvents these limitations with two key components: (1) A confidence value
aggregation encoder to sketch spectrum representations according to
amino-acid-based connectivity among MS; (2) A global-local fusion decoder to
progressively assimilate contextualized spectrum representations with a
predefined preconception of localized MS evidence and amino acid priors. Our
components originate from a closed-form solution and selectively attend to
informative amino-acid-aware MS representations. Through extensive empirical
studies, we demonstrate the superiority of DPST, showing that it outperforms
state-of-the-art approaches by a margin of 12% - 19% peptide accuracy.


  The Lagrangian formalism is developed for the population dynamics of
interacting species that are described by several well-known models. The
formalism is based on standard Lagrangians, which represent differences between
the physical kinetic and potential energy-like terms. A method to derive these
Lagrangians is presented and applied to selected theoretical models of the
population dynamics. The role of the derived Lagrangians and the energy-like
terms in the population dynamics is investigated and discussed. It is suggested
that the obtained standard Lagrangians can be used to identify physical
similarities between different population models.


  We analyse mobile-immobile transport of particles that switch between the
mobile and immobile phases with finite rates. Despite this seemingly simple
assumption of Poissonian switching we unveil a rich transport dynamics
including significant transient anomalous diffusion and non-Gaussian
displacement distributions. Our discussion is based on experimental parameters
for tau proteins in neuronal cells, but the results obtained here are expected
to be of relevance for a broad class of processes in complex systems.
Concretely, we obtain that when the mean binding time is significantly longer
than the mean mobile time, transient anomalous diffusion is observed at short
and intermediate time scales, with a strong dependence on the fraction of
initially mobile and immobile particles. We unveil a Laplace distribution of
particle displacements at relevant intermediate time scales. For any initial
fraction of mobile particles, the respective mean squared displacement displays
a plateau. Moreover, we demonstrate a short-time cubic time dependence of the
mean squared displacement for immobile tracers when initially all particles are
immobile.


  Recently, molecular communication (MC) has been studied as a micro-scale
communication between cells or molecular robots. In previous works, the MC
channels in unbounded environment was analyzed. However, many of the
experimentally implemented MC channels are surrounded by walls, thus the
boundary condition should be explicitly considered to analyze the dynamics of
MC channels. In this paper, we propose a framework to analyze the frequency
response of one-dimensional MC channels based on a diffusion equation with a
boundary. In particular, we decompose the MC channel into the diffusion system
and the boundary system, and show the relation between the cut-off frequency of
the MC channel and the communication distance based on the transfer function.
We then analyze the frequency response of a specific MC channel and reveal that
the boundary can restrict the communication bandwidth of the MC channel.


  The study of control mechanisms of biological systems allows for interesting
applications in bioengineering and medicine, for instance in cell reprogramming
or drug target identification. A control strategy often consists of a set of
interventions that, by fixing the values of some components, ensure that the
long term dynamics of the controlled system is in a desired state. A common
approach to control in the Boolean framework consists in checking how the fixed
values propagate through the network, to establish whether the effect of
percolating the interventions is sufficient to induce the target state.
Although methods based uniquely on value percolation allow for efficient
computation, they can miss many control strategies. Exhaustive methods for
control strategy identification, on the other hand, often entail high
computational costs. In order to increase the number of control strategies
identified while still benefiting from an efficient implementation, we
introduce a method based on value percolation that uses trap spaces, subspaces
of the state space that are closed with respect to the dynamics, and that can
usually be easily computed in biological networks. The approach allows for node
interventions, which fix the value of certain components, and edge
interventions, which fix the effect that one component has on another. The
method is implemented using Answer Set Programming, extending an existing
efficient implementation of value percolation to allow for the use of trap
spaces and edge control. The applicability of the approach is studied for
different control targets in a biological case study, identifying in all cases
new control strategies that would escape usual percolation-based methods.


  A key challenge in metabolomics is annotating measured spectra from a
biological sample with chemical identities. Currently, only a small fraction of
measurements can be assigned identities. Two complementary computational
approaches have emerged to address the annotation problem: mapping candidate
molecules to spectra, and mapping query spectra to molecular candidates. In
essence, the candidate molecule with the spectrum that best explains the query
spectrum is recommended as the target molecule. Despite candidate ranking being
fundamental in both approaches, no prior works utilized rank learning tasks in
determining the target molecule. We propose a novel machine learning model,
Ensemble Spectral Prediction (ESP), for metabolite annotation. ESP takes
advantage of prior neural network-based annotation models that utilize
multilayer perceptron (MLP) networks and Graph Neural Networks (GNNs). Based on
the ranking results of the MLP and GNN-based models, ESP learns a weighting for
the outputs of MLP and GNN spectral predictors to generate a spectral
prediction for a query molecule. Importantly, training data is stratified by
molecular formula to provide candidate sets during model training. Further,
baseline MLP and GNN models are enhanced by considering peak dependencies
through multi-head attention mechanism and multi-tasking on spectral topic
distributions. ESP improves average rank by 41% and 30% over the MLP and GNN
baselines, respectively, demonstrating remarkable performance gain over
state-of-the-art neural network approaches. We show that annotation
performance, for ESP and other models, is a strong function of the number of
molecules in the candidate set and their similarity to the target molecule.


  Molecule design is a fundamental problem in molecular science and has
critical applications in a variety of areas, such as drug discovery, material
science, etc. However, due to the large searching space, it is impossible for
human experts to enumerate and test all molecules in wet-lab experiments.
Recently, with the rapid development of machine learning methods, especially
generative methods, molecule design has achieved great progress by leveraging
machine learning models to generate candidate molecules. In this paper, we
systematically review the most relevant work in machine learning models for
molecule design. We start with a brief review of the mainstream molecule
featurization and representation methods (including 1D string, 2D graph, and 3D
geometry) and general generative methods (deep generative and combinatorial
optimization methods). Then we summarize all the existing molecule design
problems into several venues according to the problem setup, including input,
output types and goals. Finally, we conclude with the open challenges and point
out future opportunities of machine learning models for molecule design in
real-world applications.


  Drug-drug interaction (DDI) prediction provides a drug combination strategy
for systemically effective treatment. Previous studies usually model drug
information constrained on a single view such as the drug itself, leading to
incomplete and noisy information, which limits the accuracy of DDI prediction.
In this work, we propose a novel multi- view drug substructure network for DDI
prediction (MSN-DDI), which learns chemical substructures from both the
representations of the single drug (intra-view) and the drug pair (inter-view)
simultaneously and utilizes the substructures to update the drug representation
iteratively. Comprehensive evaluations demonstrate that MSN-DDI has almost
solved DDI prediction for existing drugs by achieving a relatively improved
accuracy of 19.32% and an over 99% accuracy under the transductive setting.
More importantly, MSN-DDI exhibits better generalization ability to unseen
drugs with a relatively improved accuracy of 7.07% under more challenging
inductive scenarios. Finally, MSN-DDI improves prediction performance for
real-world DDI applications to new drugs.


  Background School closures and distance learning have been extensively
applied to control SARS-CoV-2 transmission. Despite evidence of viral
circulation in schools, the contribution of students and of in-person schooling
to the transmission remains poorly quantified. Methods We analyze 976 exposure
events, involving 460 positive individuals, as identified in early 2021 by
routine surveillance and through an extensive screening conducted on students,
school personnel, and their household members during an outbreak in a small
municipality of Italy. Results From the analysis of potential transmission
chains, we estimated that, on average, 55.1%, 17.3% and 27.6% infection
episodes were linked to household, school, and community contacts,
respectively. Clusters originated from students or school personnel showed a
larger average cluster size (3.32 vs 1.15), a larger average number of
generations in the transmission chain (1.56 vs 1.17) and a larger set of
associated close contacts (11.3 vs 3.15, on average). We found substantial
transmission heterogeneities, with 20% positive individuals seeding 75-80 of
all transmission. A higher proportion of infected individuals causing onward
transmission was found among students (48.8% vs 29.9%, on average), who also
caused a markedly higher number of secondary cases (mean: 1.3 vs 0.5).
Conclusions Uncontrolled transmission at school could disrupt the regular
conduct of teaching activities, likely seeding the transmission into other
settings, and increasing the burden on contact-tracing operations.


  Post-transcriptional modifications are crucial for RNA function and can
affect its structure and dynamics. Force-field based classical molecular
dynamics simulations are a fundamental tool to characterize biomolecular
dynamics and their application to RNA is flourishing. Here we show that the set
of force-field parameters for N$^6$-methyladenosine (m$^6$A) developed for the
commonly used AMBER force field does not reproduce duplex denaturation
experiments and, specifically, cannot be used to describe both paired and
unpaired states. Then we use reweighting techniques to derive new parameters
matching available experimental data. The resulting force field can be used to
properly describe paired and unpaired m$^6$A in both syn and anti conformation,
and thus opens the way to the use of molecular simulations to investigate the
effects of N6 methylations on RNA structural dynamics.


  In pop culture, there are many strategies to curtail zombie apocalypses.
However, it remains unclear what routes we should take to eliminate zombies
effectively and affordably. So, we created a mathematical model to examine
interventions that armed adults with steelhead axes or provided enough
ammunition and a 9mm handgun to kill zombies in either a "single-tap" or
"double-tap". We investigate each case over two years under slow, moderate
speed, and fast zombies scenarios. We quantify health burden by zombies
averted, disability-adjusted life-years averted, and determine
cost-effectiveness using the incremental cost-effectiveness ratio. Our
predictions show the single=tap intervention is the best for stopping the
zombie apocalypse, as it would avert hundreds of millions of zombies and deaths
while also being the most cost-effective intervention. Altogether, this
suggests conserving ammunition and supplying ranged weapons would be an
effective use of limited resources in the event of a zombie uprising.


  This work reports on two results. At first we revisit the Berg and Purcell
calculation that provides a lower bound to the error in concentration
measurement by cells, by considering the realistic case when the cell starts
measuring the moment it comes in contact with the chemoattractants, instead of
measuring after equilibrating with the chemotactic concentration as done in the
classic Berg and Purcell paper. We find that the error in concentration
measurement is still the same as evaluated by Berg and Purcell. We next derive
a lower bound on measurement time below which it is not possible for the cell
to discern extra-cellular chemotactic gradients through spatial sensing
mechanisms. This bound is independent of diffusion rate and concentration of
the chemoattracts and is instead set by detachment rate of ligands from the
cell receptors. The result could help explain experimental observations.


  In recent years, there has been increasing interest in developing models and
tools to address the complex patterns of connectivity found in brain tissue.
Specifically, this is due to a need to understand how emergent properties
emerge from these network structures at multiple spatiotemporal scales. We
argue that computational models are key tools for elucidating the possible
functionalities that can emerge from interactions of heterogeneous neurons
connected by complex networks on multi-scale temporal and spatial domains. Here
we review several classes of models including spiking neurons, integrate and
fire neurons with short term plasticity (STP), conductance based
integrate-and-fire models with STP, and population density neural field (PDNF)
models using simple examples with emphasis on neuroscience applications while
also providing some potential future research directions for AI. These
computational approaches allow us to explore the impact of changing underlying
mechanisms on resulting network function both experimentally as well as
theoretically. Thus we hope these studies will inform future developments in
artificial intelligence algorithms as well as help validate our understanding
of brain processes based on experiments in animals or humans.


  Self-supervised neural language models with attention have recently been
applied to biological sequence data, advancing structure, function and
mutational effect prediction. Some protein language models, including MSA
Transformer and AlphaFold's EvoFormer, take multiple sequence alignments (MSAs)
of evolutionarily related proteins as inputs. Simple combinations of MSA
Transformer's row attentions have led to state-of-the-art unsupervised
structural contact prediction. We demonstrate that similarly simple, and
universal, combinations of MSA Transformer's column attentions strongly
correlate with Hamming distances between sequences in MSAs. Therefore,
MSA-based language models encode detailed phylogenetic relationships. We
further show that these models can separate coevolutionary signals encoding
functional and structural constraints from phylogenetic correlations reflecting
historical contingency. To assess this, we generate synthetic MSAs, either
without or with phylogeny, from Potts models trained on natural MSAs. We find
that unsupervised contact prediction is substantially more resilient to
phylogenetic noise when using MSA Transformer versus inferred Potts models.


  The discovery of nucleic acids and the structure of DNA have brought
considerable advances in the understanding of life. The development of
next-generation sequencing technologies has led to a large-scale generation of
data, for which computational methods have become essential for analysis and
knowledge discovery. In particular, RNAs have received much attention because
of the diversity of their functionalities in the organism and the discoveries
of different classes with different functions in many biological processes.
Therefore, the correct identification of RNA sequences is increasingly
important to provide relevant information to understand the functioning of
organisms. This work addresses this context by presenting a new method for the
classification of biological sequences through complex networks and entropy
maximization. The maximum entropy principle is proposed to identify the most
informative edges about the RNA class, generating a filtered complex network.
The proposed method was evaluated in the classification of different RNA
classes from 13 species. The proposed method was compared to PLEK, CPC2 and
BASiNET methods, outperforming all compared methods. BASiNETEntropy classified
all RNA sequences with high accuracy and low standard deviation in results,
showing assertiveness and robustness. The proposed method is implemented in an
open source in R language and is freely available at
https://cran.r-project.org/web/packages/BASiNETEntropy.


  Successful predictive modeling of epidemics requires an understanding of the
implicit feedback control strategies which are implemented by populations to
modulate the spread of contagion. While this task of capturing endogenous
behavior can be achieved through intricate modeling assumptions, we find that a
population's reaction to case counts can be described through a second order
affine dynamical system with linear control which fits well to the data across
different regions and times throughout the COVID-19 pandemic. The model fits
the data well both in and out of sample across the 50 states of the United
States, with comparable $R^2$ scores to state of the art ensemble predictions.
In contrast to recent models of epidemics, rather than assuming that
individuals directly control the contact rate which governs the spread of
disease, we assume that individuals control the rate at which they vary their
number of interactions, i.e. they control the derivative of the contact rate.
We propose an implicit feedback law for this control input and verify that it
correlates with policies taken throughout the pandemic. A key takeaway of the
dynamical model is that the "stable" point of case counts is non-zero, i.e.
COVID-19 will not be eradicated under the current collection of policies and
strategies, and additional policies are needed to fully eradicate it quickly.
Hence, we suggest alternative implicit policies which focus on making
interventions (such as vaccinations and mobility restrictions) a function of
cumulative case counts, for which our results suggest a better possibility of
eradicating COVID-19.


  Spatial population synchrony has been the focus of theoretical and empirical
studies for decades, in the hopes of understanding mechanisms and interactions
driving ecological dynamics. In many systems, it is well-known that seasonality
plays a critical role in the density dependence structure of the populations,
yet this has hardly received any attention in synchrony studies. Here, we
propose a protocol that allows to elucidate deterministic and stochastic
sources of spatial synchrony, while accounting for geographic- and
season-specific density dependence. We apply our protocol to seasonally-sampled
time series of sub-arctic gray-sided voles, known for marked spatial synchrony.
Dissociating seasonal density-dependence contributions to the total observed
synchrony reveals differential strength and shape of synchrony patterns by
season. Mild winter weather reveals to be an important driver of vole spatial
synchrony, with lagged effects in the fall. This has direct implications to the
future population dynamics of such species when facing climate change.


  More than half of the human population lives in cities meaning that most
people predominantly experience nature in urban greenspace. Nature exposure is
an important contributor to social, mental and physical health. As the world
faces a pandemic which threatens the physical and mental health of billions of
people, it is crucial to understand that all have the possibility to access
nature exposure to alleviate some of these challenges. Here, for the first
time, we integrate data from Facebook, Twitter, and Google Search users to show
that people looked for greenspace during COVID-19 mobility restrictions but may
not have always managed to reach it. People spent more time in areas with
greenspace when they could and that depended on the level of multiple
deprivation in the neighbourhood where the greenspace was embedded.
Importantly, while people sought greenspace throughout the first 20 months of
the pandemic, this preference intensified through the waves of lockdown. Living
in an affluent area conferred a greenspace advantage in London and Paris but we
find that in Berlin more deprived neighbourhoods sought greenspace more,
including outside their neighbourhood. This highlights the need to understand
how greenspace access and deprivation interact to create more sustainable
communities.


  Omics software tools have reshaped the landscape of modern biology and become
an essential component of biomedical research. The increasing dependence of
biomedical scientists on these powerful tools creates a need for easier
installation and greater usability. Packaging, virtualization, and
containerization are different approaches to satisfy this need by wrapping
omics tools in additional software that makes the omics tools easier to install
and use. Here, we systematically review practices across prominent packaging,
virtualization, and containerization platforms. We outline the challenges,
advantages, and limitations of each approach and some of the most widely used
platforms from the perspectives of users, software developers, and system
administrators. We also propose principles to make packaging, virtualization,
and containerization of omics software more sustainable and robust to increase
the reproducibility of biomedical and life science research.


  Correlation between weather and the transmission of SARS-CoV-2 may suggest
its seasonality. Cointegration analysis can avoid spurious correlation among
time series data. We examined the cointegration of virus transmission with
daily temperature, dewpoint, and confounding factors of mobility measurements
during the first year of the pandemic in the United States. We examined the
cointegration of the effective reproductive rate, Rt, of the virus with the
dewpoint at two meters, the temperature at two meters, Apple driving mobility,
and Google workplace mobility measurements. We found that dewpoint and Apple
driving mobility are the best factors to cointegrate with Rt, although
temperature and Google workplace mobility also cointegrate with Rt at
substantial levels. We found that the optimal lag is two days for cointegration
between Rt and weather variables, and three days for Rt and mobility. We
observed clusters of states that share similar cointegration results of Rt,
weather, and mobility, suggesting regional patterns. Our results support the
correlation of weather with the spread of SARS-CoV-2 and its potential
seasonality.


  Pandemic management requires that scientists rapidly formulate and analyze
epidemiological models in order to forecast the spread of disease and the
effects of mitigation strategies. Scientists must modify existing models and
create novel ones in light of new biological data and policy changes such as
social distancing and vaccination. Traditional scientific modeling workflows
detach the structure of a model -- its submodels and their interactions -- from
its implementation in software. Consequently, incorporating local changes to
model components may require global edits to the code-base through a manual,
time-intensive, and error-prone process. We propose a compositional modeling
framework that uses high-level algebraic structures to capture domain-specific
scientific knowledge and bridge the gap between how scientists think about
models and the code that implements them. These algebraic structures, grounded
in applied category theory, simplify and expedite modeling tasks such as model
specification, stratification, analysis, and calibration. With their structure
made explicit, models also become easier to communicate, criticize, and refine
in light of stakeholder feedback.


  In this paper, we propose a molecular communication system to localize an
abnormality in a diffusion based medium. We consider a general setup to perform
joint sensing, communication and localization. This setup consists of three
types of devices, each for a different task: mobile sensors for navigation and
molecule releasing (for communication), fusion centers (FC)s for sampling,
amplifying and forwarding the signal, and a gateway for making decision or
exchanging the information with an external device. The sensors move randomly
in the environment to reach the abnormality. We consider both collaborative and
non-collaborative sensors that simultaneously release their molecules to the
FCs when the number of activated sensors or the moving time reach a certain
threshold, respectively. The FCs amplify the received signal and forward it to
the gateway for decision making using either an ideal or a noisy communication
channel. A practical application of the proposed model is drug delivery in a
tissue of human body, in order to guide the nanomachine bound drug to the exact
location. The decision rules and probabilities of error are obtained for two
considered sensors types in both ideal and noisy communication channels.


  Optimization and expansion are two modes of staged evolution of complex
systems where macroscopic observables change at a decreasing, respectively
increasing, rate.
  A prime example of evolutionary expansion, Gross Domestic Product (GDP) time
series gauge economic activities in changing societal structures,} and the
accelerating trend of their growth probably reflects a manyfold increase of the
human interactions that drive change. We show how optimization and expansion
can coexist by replacing `wall clock time' $t$ as independent variable with a
measure of human interactions intensity $\tau$.
  Our analysis of eight centuries of yearly GDP data from three regions of
Western Europe is carried out in two steps. First, a Monte Carlo algorithm is
used to fit the GDP data to a piecewise continuous function comprising a
sequence of exponentials with different exponents.
  In a second step, GDP data are plotted vs. $\tau$ and shown to display two
logarithmic regimes, both decelerating, that are joined by a power-law
cross-over period.
  We connect the end of the first regime and the beginning of the second with
the dawn of the Industrial Revolution and the societal impact of new transport,
communication and production technologies that became widely available after
World War I.
  We conclude that wealth evolution in terms of $\tau$ is a decelerating
process with the hallmarks of record dynamics optimization.


  Imaging genetic studies aim to find associations between genetic variants and
imaging quantitative traits. Traditional genome-wide association studies (GWAS)
are based on univariate statistical tests, but when multiple traits are
analyzed together they suffer from a multiple-testing problem and from not
taking into account correlations among the traits. An alternative approach to
multi-trait GWAS is to reverse the functional relation between genotypes and
traits, by fitting a multivariate regression model to predict genotypes from
multiple traits simultaneously. However, current reverse genotype prediction
approaches are mostly based on linear models. Here, we evaluated random forest
regression (RFR) as a method to predict SNPs from imaging QTs and identify
biologically relevant associations. We learned machine learning models to
predict 518,484 SNPs using 56 brain imaging QTs. We observed that genotype
regression error is a better indicator of permutation p-value significance than
genotype classification accuracy. SNPs within the known Alzheimer disease (AD)
risk gene APOE had lowest RMSE for lasso and random forest, but not ridge
regression. Moreover, random forests identified additional SNPs that were not
prioritized by the linear models but are known to be associated with
brain-related disorders. Feature selection identified well-known brain regions
associated with AD,like the hippocampus and amygdala, as important predictors
of the most significant SNPs. In summary, our results indicate that non-linear
methods like random forests may offer additional insights into
phenotype-genotype associations compared to traditional linear multi-variate
GWAS methods.


  The prediction of the behavior of the disease, the level of affectation in a
population and the ways to control it are the most important aspects studied by
epidemiology using tools such as historical data and mathematical models. So,
our objective is (1) to provide a tool capable of analyzing epidemiological
phenomena starting from the most common social interactions within a group of
individuals. (2) To provide a methodology to build epidemiological models from
patterns and logical rules. (3) Determine the impact of social interactions on
the spread of a disease. This paper describes the logical construction of two
epidemiological models in cellular automata together with two of their
variations based on topological and dynamical principles.


  We are concerned with a special form of the (stochastic) Allen-Cahn equation,
which can be seen as a model of hybrid zones in population genetics.
Individuals in the population can be of one of three types; $aa$ are fitter
than $AA$, and both are fitter than the $aA$ heterozygotes. The hybrid zone is
the region separating a subpopulation consisting entirely of $aa$ individuals
from one consisting of $AA$ individuals. We investigate the interplay between
the motion of the hybrid zone and the shape of the habitat, both with and
without genetic drift (corresponding to stochastic and deterministic models
respectively). In the deterministic model, we investigate the effect of a wide
opening and provide some explicit sufficient conditions under which the spread
of the advantageous type is halted, and complementary conditions under which it
sweeps through the whole population. As a standing example, we are interested
in the outcome of the advantageous population passing through an isthmus. We
also identify rather precise conditions under which genetic drift breaks down
the structure of the hybrid zone, complementing previous work that identified
conditions on the strength of genetic drift under which the structure of the
hybrid zone is preserved.
  Our results demonstrate that, even in cylindrical domains, it can be
misleading to caricature allele frequencies by one-dimensional travelling
waves, and that the strength of genetic drift plays an important role in
determining the fate of a favoured allele.


  This paper analyses the stability of cycles within a heteroclinic network
lying in a three-dimensional manifold formed by six cycles, for a one-parameter
model developed in the context of game theory. We show the asymptotic stability
of the network for a range of parameter values compatible with the existence of
an interior equilibrium and we describe an asymptotic technique to decide which
cycle (within the network) is visible in numerics. The technique consists of
reducing the relevant dynamics to a suitable one-dimensional map, the so called
\emph{projective map}. Stability of the fixed points of the projective map
determines the stability of the associated cycles. The description of this new
asymptotic approach is applicable to more general types of networks and is
potentially useful in computational dynamics.


  Two alternative views of an economy are combined and studied. The first view
is that of technological evolution as a process of combinatorial innovation.
Recently a simple mathematical model (TAP) was introduced to study such a
combinatorial process. The second view is that of a network of production
functions forming an autocatalytic set. Autocatalytic (RAF) sets have been
studied extensively in the context of chemical reaction networks.
  Here, we combine the two models (TAP and RAF) and show that they are
compatible. In particular, it is shown that production function networks
resulting from the combinatorial TAP model have a high probability of
containing autocatalytic (RAF) sets. We also study the size distribution and
robustness of such "economic autocatalytic sets", and compare our results with
those from the chemical context. These initial results strongly support earlier
claims that the economy can indeed be seen as an autocatalytic set.