  Aggregated health data such as claims data from health insurances become more
and more available for research purposes. Estimates of excess mortality from
prevalence and incidence of a chronic condition have only been possible for
ages 50 years and older and have shown to be unstable in younger ages. The aim
of this article is to explore the reasons why estimates of excess mortality for
younger ages are prone to bias and what can be done to extend the age range to
ages below 50 years.


  Cells adapt to changing environments by sensing ligand concentrations using
specific receptors. The accuracy of sensing is ultimately limited by the finite
number of ligand molecules bound by receptors. Previously derived physical
limits to sensing accuracy have assumed that the concentration was constant and
ignored its temporal fluctuations. We formulate the problem of concentration
sensing in a strongly fluctuating environment as a non-linear field-theoretic
problem, for which we find an excellent approximate Gaussian solution. We
derive a new physical bound on the relative error in concentration $c$ which
scales as $\delta c/c \sim (Dac\tau)^{-1/4}$ with ligand diffusivity $D$,
receptor cross-section $a$, and characteristic fluctuation time scale $\tau$,
in stark contrast with the usual Berg and Purcell bound $\delta c/c \sim
(DacT)^{-1/2}$ for a perfect receptor sensing concentration during time $T$. We
show how the bound can be achieved by a simple biochemical network downstream
the receptor that adapts the kinetics of signaling as a function of the square
root of the sensed concentration.


  Vegetation patterns, a natural phenomenon observed worldwide, are typically
driven by spatially distributed feedback. However, the spatial colonization
mechanisms of clonal plants, driven by the growth of a rhizome, are usually not
considered in prototypical models. Here we propose a general equation for the
vegetation density that includes all main clonal-growth features as well as the
essential ingredients leading to spatial self-organization. This generic model
reproduces the phase diagram of a fully detailed model of clonal growth. The
relation of each term of the model with the mechanisms of clonal growth is
discussed.


  The catalytic reaction system (CRS) formalism by Hordijk and Steel is a
versatile method to model autocatalytic biochemical reaction networks. It is
particularly suited, and has been widely used, to study self-sustainment and
self-generation properties. Its distinguishing feature is the explicit
assignment of a catalytic function to chemicals that are part of the system. In
this work, it is shown that the subsequent and simultaneous catalytic functions
gives rise to an algebraic structure of a semigroup with the additional
compatible operation of idempotent addition and a partial order. The aim of
this article is to demonstrate that such semigroup models are a natural setup
to describe and analyze self-sustaining CRS. The basic algebraic properties of
the models are established and the notion of the function of any set of
chemicals on the whole CRS is made precise. This leads to a natural discrete
dynamics on the network, which results from iteratively considering the
self-action on a set of chemicals by its own function. The fixed points of this
dynamics are proven to correspond to self-sustaining sets of chemicals, which
are functionally closed. Finally, as the main application, a theorem on the
maximal self-sustaining set and a structure theorem on the lattice of
functionally closed self-sustaining sets of chemicals are proven.


  This Thesis presents research at the boundary between Statistical Physics and
Biology. First, we have devised a class of Boolean constraint satisfaction
problems (CSP) whose solutions describe the feasible operational states of a
chemical reaction network. After developing statistical mechanics techniques to
generate solutions and studying the properties of the solution space for both
ensembles and individual instances of random reaction networks, we have applied
this framework to the metabolic network of the bacterium E.Coli. Results
highlight, on one hand, a complex organization of operational states into
"modules" involving different biochemically-defined pathways, and, on the
other, a high degree of cross-talk between modules. In summary, we propose that
this class of CSPs may provide novel and useful quantitative information
linking structure to function in cellular reaction networks.


  Flexibility is an intrinsic essential feature of protein structures, directly
linked to their functions. To this day, most of the prediction methods use the
crystallographic data (namely B-factors) as the only indicator of protein's
inner flexibility and predicts them as rigid or flexible.PredyFlexy stands
differently from other approaches as it relies on the definition of protein
flexibility (i) not only taken from crystallographic data, but also (ii) from
Root Mean Square Fluctuation (RMSFs) observed in Molecular Dynamics
simulations. It also uses a specific representation of protein structures,
named Long Structural Prototypes (LSPs). From Position-Specific Scoring Matrix,
the 120 LSPs are predicted with a good accuracy and directly used to predict
(i) the protein flexibility in three categories (flexible, intermediate and
rigid), (ii) the normalized B-factors, (iii) the normalized RMSFs, and (iv) a
confidence index. Prediction accuracy among these three classes is equivalent
to the best two class prediction methods, while the normalized B-factors and
normalized RMSFs have a good correlation with experimental and in silico
values. Thus, PredyFlexy is a unique approach, which is of major utility for
the scientific community. It support parallelization features and can be run on
a local cluster using multiple cores.The entire project is available under an
open-source license at
http://www.dsimb.inserm.fr/~debrevern/TOOLS/predyflexy_1.3/index.php.


  Post-Translational Modifications (PTMs) are known to play a critical role in
the regulation of the protein functions. Their impact on protein structures,
and their link to disorder regions have already been spotted on the past
decade. Nonetheless, the high diversity of PTMs types, and the multiple schemes
of protein modifications (multiple PTMs, of different types, at different time,
etc) make difficult the direct confrontation of PTM annotations and protein
structures data.We so analyzed the impact of the residue modifications on the
protein structures at local level. Thanks to a dedicated structure database,
namely PTM-SD, a large screen of PTMs have been done and analyze at a local
protein conformation levels using the structural alphabet Protein Blocks (PBs).
We investigated the relation between PTMs and the backbone conformation of
modified residues, of their local environment, and at the level of the complete
protein structure. The two main PTM types (N-glycosylation and phosphorylation)
have been studied in non-redundant datasets, and then, 4 different proteins
were focused, covering 3 types of PTMs: N-glycosylation in renin endopeptidase
and liver carboxylesterase, phosphorylation in cyclin-dependent kinase 2
(CDK2), and methylation in actin. We observed that PTMs could either stabilize
or destabilize the backbone structure, at a local and global scale, and that
these effects depend on the PTM types.


  The aim of this work was to find a minimal set of structurally stable
pentapeptides, which allows forming a polypeptide chain of a required 3D
structure. To search for factors that ensure structural stability of the
pentapeptide, we generated peptide sequences with no more than three functional
groups, based on the alanine pentapeptide AAAAA. We analyzed 44,860 structures
of peptides by the molecular dynamics method and found that 1,225 pentapeptides
over 80% of the simulation time were in a stable conformation. Clustering of
these conformations revealed 54 topological types of conformationally stable
pentapeptides. These conformations relate to different combined elements of the
protein secondary structure. So, we obtained a minimal set of amino acid
structures of conformationally stable pentapeptides, creating a complete set of
different topologies that ensure the formation of pre-folded conformation of
protein structures.


  Protein aggregation on the plasma membrane (PM) is of critical importance to
many cellular processes such as cell adhesion, endocytosis, fibrillar
conformation, and vesicle transport. Lateral diffusion of protein aggregates or
clusters on the surface of the PM plays an important role in governing their
heterogeneous surface distribution. However, the stability behavior of the
surface distribution of protein aggregates remains poorly understood.
Therefore, understanding the spatial patterns that can emerge on the PM solely
through protein-protein interaction, lateral diffusion, and feedback is an
important step towards a complete description of the mechanisms behind protein
clustering on the cell surface. In this work, we investigate the pattern
formation of a reaction-diffusion model that describes the dynamics of a system
of ligand-receptor complexes. The purely diffusive ligand in the cytosol can
bind receptors in the PM, and the resultant ligand-receptor complexes not only
diffuse laterally but can also form clusters resulting in different oligomers.
Finally, the largest oligomers recruit ligands from the cytosol in a positive
feedback. From a methodological viewpoint, we provide theoretical estimates for
diffusion-driven instabilities of the protein aggregates based on the Turing
mechanism. Our main result is a threshold phenomenon, in which a sufficiently
high recruitment of ligands promotes the input of new monomeric components and
consequently drives the formation of a single-patch spatially heterogeneous
steady-state.


  Mobility restrictions - travel advisories, trade and travel bans, border
closures and, in extreme cases, area quarantines or cordons sanitaires - are
among the most widely used measures to control infectious diseases.
Restrictions of this kind were important in the response to epidemics of SARS
(2003), H1N1 influenza (2009), and Ebola (2014). However, they do not always
work as expected. The imposition of a cordon sanitaire to control the 2014 West
African Ebola outbreak, for example, is argued to have led to a
higher-than-expected number of cases in the quarantined area. To determine when
mobility restrictions reduce the size of an epidemic, we use a model of disease
transmission within and between economically heterogeneous locally connected
communities. One community comprises a low-risk, resource-rich, low-density
population with access to effective medical resources. The other comprises a
high-risk, resource-poor, high-density population without access to effective
medical resources. We find that the overall size of an epidemic centered in the
high-risk community is sensitive to the stringency of mobility restrictions
between the two communities. Unrestricted mobility between the two risk
communities increases the number of secondary cases in the low-risk community
but reduces the overall epidemic size. By contrast, the imposition of a cordon
sanitaire around the high-risk community reduces the number of secondary
infections in the low-risk community but increases the overall epidemic size.
The degree to which mobility restrictions increase or decrease the overall
epidemic size depends on the level of risk in each community and the
characteristics of the disease.


  Natural selection has shaped the evolution of cells and multi-cellular
organisms such that social cooperation can often be preferred over an
individualistic approach to metabolic regulation. This paper extends a
framework for dynamic metabolic resource allocation based on the maximum
entropy principle to spatiotemporal models of metabolism with cooperation. Much
like the maximum entropy principle encapsulates `bet-hedging' behaviour
displayed by organisms dealing with future uncertainty in a fluctuating
environment, its cooperative extension describes how individuals adapt their
metabolic resource allocation strategy to further accommodate limited knowledge
about the welfare of others within a community. The resulting theory explains
why local regulation of metabolic cross-feeding can fulfil a community-wide
metabolic objective if individuals take into consideration an ensemble measure
of total population performance as the only form of global information. The
latter is likely supplied by quorum sensing in microbial systems or signalling
molecules such as hormones in multi-cellular eukaryotic organisms.


  Inspired by observations of beads packed on a thin string in such systems as
sea-grapes and dental plaque, we study the random sequential adsorption of
spheres on a cylinder. We determine the asymptotic fractional coverage of the
cylinder as a function of the sole parameter in the problem, the ratio of the
sphere radius to the cylinder radius (for a very long cylinder) using a
combination of analysis and numerical simulations. Examining the asymptotic
structures, we find weak chiral ordering on sufficiently small spatial scales.
Experiments involving colloidal microspheres that can attach irreversibly to a
silica wire via electrostatic forces or DNA hybridization allow us to verify
our predictions for the asymptotic coverage.


  How specific environmental contexts contribute to the robustness and
variation of developmental trajectories and evolutionary transitions is a
central point in Eco-Evo-Devo. However, the articulation of ecological,
evolutionary and developmental processes into integrative frameworks has been
elusive, partly because standard experimental designs neglect or oversimplify
ecologically meaningful contexts. Microbial models are useful to expose and
discuss two possible sources of bias associated with gene-centered experimental
designs: the use of laboratory strains and laboratory environmental conditions.
We illustrate our point by showing how contrasting developmental phenotypes in
Myxococcus xanthus depend on the joint variation of temperature and substrate
stiffness. Microorganismal development can provide key information for better
understanding the role of environmental conditions in the evolution of
developmental variation, and to overcome some of the limitations associated
with current experimental approaches.


  Biological systems exhibit processes on a wide range of time and length
scales. This work demonstrates that models, wherein the interaction between
system constituents is captured by algebraic operations, inherently allow for
successive coarse-graining operations through quotients of the algebra.
Thereby, the class of model is retained and all possible coarse-graining
operations are encoded in the lattice of congruences of the model. We analyze a
class of algebraic models generated by the subsequent and simultaneous
catalytic functions of chemicals within a reaction network. Our ansatz yields
coarse-graining operations that cover the network with local functional patches
and delete the information about the environment, and complementary operations
that resolve only the large-scale functional structure of the network. Finally,
we present a geometric interpretation of the algebraic models through an
analogy with classical models on vector fields. We then use the geometric
framework to show how a coarse-graining of the algebraic model naturally leads
to a coarse-graining of the state-space. The framework developed here is aimed
at the study of the functional structure of cellular reaction networks spanning
a wide range of scales.


  Critical transitions or regime shifts are sudden and unexpected changes in
the state of an ecosystem, that are usually associated with dangerous levels of
environmental change. However, recent studies show that critical transitions
can also be triggered by dangerous rates of environmental change. In contrast
to classical regime shifts, such rate-induced critical transitions do not
involve any obvious loss of stability, or a bifurcation, and thus cannot be
explained by the linear stability analysis. In this work, we demonstrate that
the well-known Rosenzweig-MacArthur predator-prey model can undergo a
rate-induced critical transition in response to a continuous decline in the
habitat quality, resulting in a collapse of the predator and prey populations.
Rather surprisingly, the collapse occurs even if the environmental change is
slower than the slowest process in the model. To explain this counterintuitive
phenomenon, we combine methods from geometric singular perturbation theory with
the concept of a moving equilibrium, and study critical rates of environmental
change with dependence on the initial state and the system parameters.
Moreover, for a fixed rate of environmental change, we determine the set of
initial states that undergo a rate-induced population collapse. Our results
suggest that ecosystems may be more sensitive to how fast environmental
conditions change than previously assumed. In particular, unexpected critical
transitions with dramatic ecological consequences can be triggered by
environmental changes that (i) do not exceed any dangerous levels, and (ii) are
slower than the natural timescales of the ecosystem. This poses an interesting
research question whether regime shifts observed in the natural world are
predominantly rate-induced or bifurcation-induced.


  Consider stochastic models for the spread of an infection in a structured
community, where this structured community is itself described by a random
network model. Some common network models and transmission models are defined
and large population proporties of them are presented. Focus is then shifted to
statistical methodology: what can be estimated and how, depending on the
underlying network, transmission model and the available data? This survey
paper discusses several different scenarios, also giving references to
publications where more details can be found.


  Traditional anatomical analyses captured only a fraction of real phenomic
information. Here, we apply deep learning to quantify total phenotypic
similarity across 2468 butterfly photographs, covering 38 subspecies from the
polymorphic mimicry complex of $\textit{Heliconius erato}$ and
$\textit{Heliconius melpomene}$. Euclidean phenotypic distances, calculated
using a deep convolutional triplet network, demonstrate significant convergence
between interspecies co-mimics. This quantitatively validates a key prediction
of M\"ullerian mimicry theory, evolutionary biology's oldest mathematical
model. Phenotypic neighbor-joining trees are significantly correlated with wing
pattern gene phylogenies, demonstrating objective, phylogenetically informative
phenome capture. Comparative analyses indicate frequency-dependent, mutual
convergence with coevolutionary exchange of wing pattern features. Therefore,
phenotypic analysis supports reciprocal coevolution, predicted by classical
mimicry theory but since disputed, and reveals mutual convergence as an
intrinsic generator for the surprising diversity of M\"ullerian mimicry. This
demonstrates that deep learning can generate phenomic spatial embeddings which
enable quantitative tests of evolutionary hypotheses previously only testable
subjectively.


  Inverted solubility--a crystal melting upon cooling--is observed in a handful
of proteins, such as carbomonoxy hemoglobin and $\gamma$D-crystallin. In human
$\gamma$D-crystallin, the phenomenon is associated with the mutation of the
23$^\mathrm{rd}$ residue, a proline, to a threonine, serine or valine. One
proposed microscopic mechanism for this effect entails an increase in
hydrophobicity upon mutagenesis. Recent crystal structures of a double mutant
that includes the P23T mutation allows for a more careful investigation of this
proposal. Here, we first measure the surface hydrophobicity of various mutant
structures of this protein and determine that it does not discernibly increase
upon the mutating the 23$^\mathrm{rd}$ residue. We then investigate the
solubility inversion regime with a schematic patchy particle model that
includes one of three models for temperature-dependent patch energies: two of
the hydrophobic effect, and a more generic description. We conclude that while
solubility inversion due to the hydrophobic effect may be possible, microscopic
evidence to support it in $\gamma$D-crystallin is weak. More generally, we find
that solubility inversion requires a fine balance between patch strengths and
the temperature-dependent contribution, which may explain why inverted
solubility is not commonly observed in proteins. In any event, we also find
that the temperature-dependent interaction has only a negligible impact on the
critical properties of the $\gamma$D-crystallin, in line with previous
experimental observations.


  Understanding chromatin organization and dynamics is important since they
crucially affect DNA functions. In this study, we investigate chromatin
dynamics by statistically analyzing single-nucleosome movement in living human
cells. Bi-modal nature of the mean squared displacement distribution of
nucleosomes allows for a natural categorization of the nucleosomes as fast and
slow. Analyses of the nucleosome-nucleosome correlation functions within these
categories along with the density of vibrational modes show that the
nucleosomes form dynamically correlated fluid regions, i.e., dynamic domains of
fast and slow nucleosomes. Perturbed nucleosome dynamics by global histone
acetylation or cohesin inactivation indicate that nucleosome-nucleosome
interactions along with tethering of chromatin chains organize nucleosomes into
fast and slow dynamic domains. A simple polymer model is introduced, which
shows the consistency of this dynamic domain picture. Statistical analyses of
single-nucleosome movement provide rich information on how chromatin is
dynamically organized in a fluid manner in living cells.


  Apart from the knots formed by the main-chain, the proteins can form numerous
topological structures, when included the covalent and ion-mediated
interactions. In this work, we define the protein non-trivial $\theta$-curves
and identify 7 different topologies in all structures known up to date. We
study the correlation of the motif with the function and organism of origin,
and pointing the similarity with main-chain knots, we show that some motifs may
indeed be functional. We also analyze the folding and bridge-induced stability
of an exemplary protein with $\theta$-curve motif and provide a catalogue of
possible $\theta$-curves in proteins.


  Cooperative behaviour constitutes a key aspect of both human society and
non-human animal systems, but explaining how cooperation evolves represents a
major scientific challenge. It is now well established that social network
structure plays a central role for the viability of cooperation. However, not
much is known about the importance of the positions of cooperators in the
networks for the evolution of cooperation. Here, we investigate how cooperation
is affected by correlations between cooperativeness and individual social
connectedness. Using simulation models, we find that the effect of correlation
between cooperativeness and connectedness (degree) depends on the social
network structure, with positive effect in standard scale-free networks and no
effect in standard Poisson networks. Furthermore, when degree assortativity is
increased such that individuals cluster with others of similar social
connectedness, we find that bridge areas between social clusters can act as
barriers to the spread of defection, leading to strong enhancement of
cooperation in particular in Poisson networks. But this effect is sensitive to
the presence of Trojan horses (defectors placed within cooperator clusters).
The study provides new knowledge about the conditions under which cooperation
may evolve and persist, and the results are also relevant to consider in regard
to human cooperation experiments.


  A great number of biological organisms live in aqueous environments. Major
evolutionary transitions, including the emergence of life itself, likely
occurred in such environments. While the chemical aspects of the role of water
in biology are well-studied, the effects of water's physical characteristics on
evolutionary events, such as the control of population structure via its rich
transport properties, are less clear. Evolutionary transitions such as the
emergence of the first cells and of multicellularity, require cooperation among
groups of individuals. However, evolution of cooperation faces challenges in
unstructured "well-mixed" populations, as parasites quickly overwhelm
cooperators. Models that assume population structure to promote cooperation
envision such structure to arise from spatial "lattice" models (e.g. surface
bound individuals) or compartmentalization models, often realized as
protocells. Here we study the effect of turbulent motions in spatial models,
and propose that coherent structures, i.e. flow patterns which trap fluid and
arise naturally in turbulent flows, may serve many of the properties associated
with compartments--collocalization, division, and merging--and thought to play
a key role in the origins of life and other evolutionary transitions. These
results suggest that group selection models may be applicable with fewer
physical and chemical constraints than previously thought, and apply much more
widely in aqueous environments.


  We develop a generalized group-based epidemic model (GgroupEM) framework for
any compartmental epidemic model (for example;
susceptible-infected-susceptible, susceptible-infected-recovered,
susceptible-exposed-infected-recovered). Here, a group consists of a collection
of individual nodes. This model can be used to understand the important dynamic
characteristics of a stochastic epidemic spreading over very large complex
networks, being informative about the state of groups. Aggregating nodes by
groups, the state space becomes smaller than the individual-based approach at
the cost of aggregation error, which is strongly bounded by the isoperimetric
inequality. We also develop a mean-field approximation of this framework to
further reduce the state-space size. Finally, we extend the GgroupEM to
multilayer networks. Since the group-based framework is computationally less
expensive and faster than an individual-based framework, then this framework is
useful when the simulation time is important.


  In real social networks, person-to-person interactions are known to be
heterogeneous, which can affect the way a disease spreads through a population,
reaches a tipping point in the fraction of infected individuals, and becomes an
epidemic. This property, called disorder, is usually associated with contact
times between individuals and can be modeled by a weighted network, where the
weights are related to normalized contact times $\omega$. In this paper, we
study the SIR model for disease spreading when both close and distant types of
interactions are present. We develop a mitigation strategy that reduces only
the time duration of distant contacts, which are easier to alter in practice.
Using branching theory, supported by simulations, we found that the
effectiveness of the strategy increases when the density $f_1$ of close
contacts decreases. Moreover, we found a threshold $\tilde{f}_1 = T_c / \beta$
below which the strategy can bring the system from an epidemic to a
non-epidemic phase, even when close contacts have the longest time durations.


  Different aspects of a clinical sample can be revealed by multiple types of
omics data. Integrated analysis of multi-omics data provides a comprehensive
view of patients, which has the potential to facilitate more accurate clinical
decision making. However, omics data are normally high dimensional with large
number of molecular features and relatively small number of available samples
with clinical labels. The "dimensionality curse" makes it challenging to train
a machine learning model using high dimensional omics data like DNA methylation
and gene expression profiles. Here we propose an end-to-end deep learning model
called OmiVAE to extract low dimensional features and classify samples from
multi-omics data. OmiVAE combines the basic structure of variational
autoencoders with a classification network to achieve task-oriented feature
extraction and multi-class classification. The training procedure of OmiVAE is
comprised of an unsupervised phase without the classifier and a supervised
phase with the classifier. During the unsupervised phase, a hierarchical
cluster structure of samples can be automatically formed without the need for
labels. And in the supervised phase, OmiVAE achieved an average classification
accuracy of 97.49% after 10-fold cross-validation among 33 tumour types and
normal samples, which shows better performance than other existing methods. The
OmiVAE model learned from multi-omics data outperformed that using only one
type of omics data, which indicates that the complementary information from
different omics datatypes provides useful insights for biomedical tasks like
cancer classification.


  Vital biological processes such as genome repair require fast and efficient
binding of selected proteins to specific target sites on DNA. Here we propose
an active target search mechanism based on "chromophoresis", the dynamics of
DNA-binding proteins up or down gradients in the density of epigenetic marks,
or colours (biochemical tags on the genome). We focus on a set of proteins that
deposit marks from which they are repelled---a case which is only encountered
away from thermodynamic equilibrium. For suitable ranges of kinetic parameter
values, chromophoretic proteins can perform unidirectional motion and are
optimally redistributed along the genome. Importantly, they can also locally
unravel a region of the genome which is collapsed due to self-interactions and
"dive" deep into its core, for a striking enhancement of the efficiency of
target search on such an inaccessible substrate. We discuss the potential
relevance of chromophoresis for the location of DNA lesions.


  The prevalence of bacteria in the atmosphere has been well established in
relevant literature, suggesting that airborne bacteria can influence
atmospheric characteristics including the development of clouds. Studies have
also demonstrated that the atmospheric biological profile is influenced by the
underlying terrestrial biomes. An understanding of the complex interplay of
factors that can influence the atmospheric biological profile, not to mention
developing a biological census of the atmosphere, requires a cost-effective
experimental system capable of generating reproducible results with reliable
data. However, as has been demonstrated by payloads launched by space agencies
such as NASA and JAXA, these payloads are both complex and cost prohibitive.
This paper discusses the design and implementation of a biologically oriented
experimental payload for high-altitude ballooning that is within the means of
most student-run experimental programs. The payload highlighted in this
presentation, PHANTOM (Probe for High Altitude Numeration and Tracking of
Microorganisms, which has the goal of capturing aerial microorganisms at
multiple altitudes in order to characterize the biological composition of the
upper atmosphere), has undergone a number of successful flight trials, and
serves to highlight the feasibility and utility of interdisciplinary projects
between aerospace and the biological sciences.


  Late diagnosis and high costs are key factors that negatively impact the care
of cancer patients worldwide. Although the availability of biological markers
for the diagnosis of cancer type is increasing, costs and reliability of tests
currently present a barrier to the adoption of their routine use. There is a
pressing need for accurate methods that enable early diagnosis and cover a
broad range of cancers. The use of machine learning and RNA-seq expression
analysis has shown promise in the classification of cancer type. However,
research is inconclusive about which type of machine learning models are
optimal. The suitability of five algorithms were assessed for the
classification of 17 different cancer types. Each algorithm was fine-tuned and
trained on the full array of 18,015 genes per sample, for 4,221 samples (75 %
of the dataset). They were then tested with 1,408 samples (25 % of the dataset)
for which cancer types were withheld to determine the accuracy of prediction.
The results show that ensemble algorithms achieve 100% accuracy in the
classification of 14 out of 17 types of cancer. The clustering and
classification models, while faster than the ensembles, performed poorly due to
the high level of noise in the dataset. When the features were reduced to a
list of 20 genes, the ensemble algorithms maintained an accuracy above 95% as
opposed to the clustering and classification models.


  Infectious disease models can be of great use for understanding the
underlying mechanisms that influence the spread of diseases and predicting
future disease progression. Modeling has been increasingly used to evaluate the
potential impact of different control measures and to guide public health
policy decisions. In recent years, there has been rapid progress in developing
spatio-temporal modeling of infectious diseases and an example of such recent
developments is the discrete time individual-level models (ILMs). These models
are well developed and provide a common framework for modeling many disease
systems, however, they assume the probability of disease transmission between
two individuals depends only on their spatial separation and not on their
spatial locations. In cases where spatial location itself is important for
understanding the spread of emerging infectious diseases and identifying their
causes, it would be beneficial to incorporate the effect of spatial location in
the model. In this study, we thus generalize the ILMs to a new class of
geographically-dependent ILMs (GD-ILMs), to allow for the evaluation of the
effect of spatially varying risk factors (e.g., education, environmental), as
well as unobserved spatial structure, upon the transmission of infectious
disease. Specifically, we consider a conditional autoregressive model to
capture the effects of unobserved spatially structured latent covariates or
measurement error. This results in flexible infectious disease models that can
be used for formulating etiological hypotheses and identifying geographical
regions of unusually high risk to formulate preventive action. The reliability
of these models are investigated on a combination of simulated epidemic data
and Alberta seasonal influenza outbreak data (2009). This new class of models
is fitted to data within a Bayesian statistical framework using MCMC methods.


  Parametric sensitivity analysis is a critical component in the study of
mathematical models of physical systems. Due to its simplicity, finite
difference methods are used extensively for this analysis in the study of
stochastically modeled reaction networks. Different coupling methods have been
proposed to build finite difference estimators, with the "split coupling," also
termed the "stacked coupling," yielding the lowest variance in the vast
majority of cases. Analytical results related to this coupling are sparse, and
include an analysis of the variance of the coupled processes under the
assumption of globally Lipschitz intensity functions [Anderson, SIAM Numerical
Analysis, Vol. 50, 2012]. Because of the global Lipschitz assumption utilized
in [Anderson, SIAM Numerical Analysis, Vol. 50, 2012], the main result there is
only applicable to a small percentage of the models found in the literature,
and it was conjectured that similar results should hold for a much wider class
of models. In this paper we demonstrate this conjecture to be true by proving
the variance of the coupled processes scales in the desired manner for a large
class of non-Lipschitz models. We further extend the analysis to allow for time
dependence in the parameters. In particular, binary systems with or without
time-dependent rate parameters, a class of models that accounts for the vast
majority of systems considered in the literature, satisfy the assumptions of
our theory.


  We present a computational simulation of evolving entities that includes
symbiosis with shifting levels of selection. Evolution by natural selection
shifts from the level of the original entities to the level of the new
symbiotic entity. In the simulation, the fitness of an entity is measured by a
series of one-on-one competitions in the Immigration Game, a two-player
variation of Conway's Game of Life. Mutation, reproduction, and symbiosis are
implemented as operations that are external to the Immigration Game. Because
these operations are external to the game, we are able to freely manipulate the
operations and observe the effects of the manipulations. The simulation is
composed of four layers, each layer building on the previous layer. The first
layer implements a simple form of asexual reproduction, the second layer
introduces a more sophisticated form of asexual reproduction, the third layer
adds sexual reproduction, and the fourth layer adds symbiosis. The experiments
show that a small amount of symbiosis, added to the other layers, significantly
increases the fitness of the population. We suggest that the model may provide
new insights into symbiosis in biological and cultural evolution.


  This article aims at exploring the short-wavelength stabilization and
destabilization of the advection-diffusion systems formulated using the
Patlak-Keller-Segel cross-diffusion. We study a model of the taxis partly
driven by an external signal. We address the general short-wavelength signal
using the homogenization technique, and then we give a detailed analysis of the
signals emitted as the travelling waves. It turns out that homogenizing
produces the drift of species, which is the main translator of the external
signal effects, in particular, on the stability issues. We examine the
stability of the quasi-equilibria - that is, the simplest short-wavelength
patterns fully imposed by the external signal. Comparing the results to the
case of switching the signal off allows us to estimate the effect of it. For
instance, the effect of the travelling wave turns out to be not single-valued
but depending on the wave speed. Namely, there is an independent threshold
value such that increasing the amplitude of the wave destabilizes the
quasi-equilibria provided that the wave speed is above this value. Otherwise,
the same action exerts the opposite effect. It is worth to note that the effect
is exponential in the amplitude of the wave in both cases.


  On June 25th, 2018, Huang et al. published a computational method SAVER on
Nature Methods for imputing dropout gene expression levels in single cell RNA
sequencing (scRNA-seq) data. Huang et al. performed a set of comprehensive
benchmarking analyses, including comparison with the data from RNA fluorescence
in situ hybridization, to demonstrate that SAVER outperformed two existing
scRNA-seq imputation methods, scImpute and MAGIC. However, their computational
analyses were based on semi-synthetic data that the authors had generated
following the Poisson-Gamma model used in the SAVER method. We have therefore
re-examined Huang et al.'s study. We find that the semi-synthetic data have
very different properties from those of real scRNA-seq data and that the cell
clusters used for benchmarking are inconsistent with the cell types labeled by
biologists. We show that a reanalysis based on real scRNA-seq data and grounded
on biological knowledge of cell types leads to different results and
conclusions from those of Huang et al.


  We present a linear stability analysis of stationary states (or fixed points)
in large dynamical systems defined on random directed graphs with a prescribed
distribution of indegrees and outdegrees. We obtain two remarkable results for
such dynamical systems: First, infinitely large systems on directed graphs can
be stable even when the degree distribution has unbounded support; this result
is surprising since their counterparts on nondirected graphs are unstable when
system size is large enough. Second, we show that the phase transition between
the stable and unstable phase is universal in the sense that it depends only on
a few parameters, such as, the mean degree and a degree correlation
coefficient. In addition, in the unstable regime we characterize the nature of
the destabilizing mode, which also exhibits universal features. These results
follow from an exact theory for the leading eigenvalue of infinitely large
graphs that are locally tree-like and oriented, as well as, for the right and
left eigenvectors associated with the leading eigenvalue. We corroborate
analytical results for infinitely large graphs with numerical experiments on
random graphs of finite size. We discuss how the presented theory can be
extended to graphs with diagonal disorder and to graphs that contain
nondirected links. Finally, we discuss the influence of small cycles and how
they can destabilize large dynamical systems when they induce strong enough
feedback loops.


  Simulating drug binding and unbinding is a challenge, as the rugged energy
landscapes that separate bound and unbound states require extensive sampling
that consumes significant computational resources. Here, we describe the use of
interactive molecular dynamics in virtual reality (iMD-VR) as an accurate
low-cost strategy for flexible protein-ligand docking. We outline an
experimental protocol which enables expert iMD-VR users to guide ligands into
and out of the binding pockets of trypsin, neuraminidase, and HIV-1 protease,
and recreate their respective crystallographic protein-ligand binding poses
within 5 - 10 minutes. Following a brief training phase, our studies shown that
iMD-VR novices were able to generate unbinding and rebinding pathways on
similar timescales as iMD-VR experts, with the majority able to recover binding
poses within 2.15 Angstrom RMSD of the crystallographic binding pose. These
results indicate that iMD-VR affords sufficient control for users to carry out
the detailed atomic manipulations required to dock flexible ligands into
dynamic enzyme active sites and recover crystallographic poses, offering an
interesting new approach for simulating drug docking and generating binding
hypotheses.


  Common and complex traits are the consequence of the interaction and
regulation of multiple genes simultaneously, which work in a coordinated way.
However, the vast majority of studies focus on the differential expression of
one individual gene at a time. Here, we aim to provide insight into the
underlying relationships of the genes expressed in the human brain in cases
with schizophrenia (SCZ) and controls. We introduced a novel approach to
identify differential gene regulatory patterns and identify a set of essential
genes in the brain tissue. Our method integrates genetic, transcriptomic, and
Hi-C data and generates a transcriptomic-causal network. Employing this
approach for analysis of RNA-seq data from CommonMind Consortium, we identified
differential regulatory patterns for SCZ cases and control groups to unveil the
mechanisms that control the transcription of the genes in the human brain. Our
analysis identified modules with a high number of SCZ-associated genes as well
as assessing the relationship of the hubs with their down-stream genes in both,
cases and controls. In addition, the results identified essential genes for
brain function and suggested new genes putatively related to SCZ.


  Dimeric molecular motors walk on polar tracks by binding and hydrolyzing one
ATP per step. Despite tremendous progress, the waiting state for ATP binding in
the well-studied kinesin that walks on microtubule (MT), remains controversial.
One experiment suggests that in the waiting state both heads are bound to the
MT, while the other shows that ATP binds to the leading head after the partner
head detaches. To discriminate between these two scenarios, we developed a
theory to calculate accurately several experimentally measurable quantities as
a function of ATP concentration and resistive force.
  In particular, we predict that measurement of the randomness parameter could
discriminate between the two scenarios for the waiting state of kinesin,
thereby resolving this standing controversy.


  The rise in the availability of bacterial genomes defines a need for
synthesis: abstracting from individual taxa, to see larger patterns of
bacterial lifestyles across systems. A key concept for such synthesis in
ecology is the niche, the set of capabilities that enables a population's
persistence and defines its impact on the environment. The set of possible
niches forms the niche space, a conceptual space delineating ways in which
persistence in a system is possible. Here we use manifold learning to map the
space of metabolic networks representing thousands of bacterial genera. The
results reveal a metabolic niche space with a complex branching geometry, whose
branches constitute major strategies spanning life in different habitats and
hosts. We further demonstrate that communities from similar ecosystem types map
to characteristic regions of this new functional coordinate system, permitting
ecological descriptions of microbiomes in terms of large scale metabolic roles
that may be filled.


  We present the assembly category assessment in the 13th edition of the CASP
community-wide experiment. For the second time, protein assemblies constitute
an independent assessment category. Compared to the last edition we see a clear
uptake in participation, more oligomeric targets released, and consistent,
albeit modest, improvement of the predictions quality. Looking at the tertiary
structure predictions we observe that ignoring the oligomeric state of the
targets hinders modelling success. We also note that some contact prediction
groups successfully predicted homomeric interfacial contacts, though it appears
that these predictions were not used for assembly modelling. Homology modelling
with sizeable human intervention appears to form the basis of the assembly
prediction techniques in this round of CASP. Future developments should see
more integrated approaches to modelling where multiple subunits are a natural
part of the modelling process, which would benefit the structure prediction
field as a whole.


  We investigate a two-strain disease model with amplification to simulate the
prevalence of drug-susceptible (s) and drug-resistant (m) disease strains. We
model the emergence of drug resistance as a consequence of inadequate
treatment, i.e. amplification. We perform a dynamical analysis of the resulting
system and find that the model contains three equilibrium points: a
disease-free equilibrium; a mono-existent disease-endemic equilibrium with
respect to the drug-resistant strain; and a co-existent disease-endemic
equilibrium where both the drug-susceptible and drug-resistant strains persist.
We found two basic reproduction numbers: one associated with the
drug-susceptible strain $R_{0s}$; the other with the drug-resistant strain
$R_{0m}$,and showed that at least one of the strains can spread in a population
if ($R_{0s}$,$R_{0m}$) > 1 (epidemic).Furthermore, we also showed that if
$R_{0m}$ > max($R_{0s}$,1), the drug-susceptible strain dies out but the
drug-resistant strain persists in the population; however if $R_{0s}$ >
max($R_{0m}$,1), then both the drug-susceptible and drug-resistant strains
persist in the population. We conducted a local stability analysis of the
system equilibrium points using the Routh-Hurwitz conditions and a global
stability analysis using appropriate Lyapunov functions. Sensitivity analysis
was used to identify the most important model parameters through the partial
rank correlation coefficient (PRCC) method. We found that the contact rate of
both strains had the largest influence on prevalence. We also investigated the
impact of amplification and treatment rates of both strains on the equilibrium
prevalence of infection; results suggest that poor quality treatment make
coexistence more likely but increase the relative abundance of resistant
infections.


  Here, we study the evolution of specialization using realistic computer
simulations of bacteria that secrete two public goods in a dynamic fluid.
Through this first principles approach, we find physical factors such as
diffusion, flow patterns, and decay rates are as influential as fitness
economics in governing the evolution of community structure, to the extent that
when mechanical factors are taken into account, (1) Generalist communities can
resist becoming specialists, despite the invasion fitness of specialization,
(2) Generalist and specialists can both resist cheaters despite the invasion
fitness of free-riding, (3) Multiple community structures can coexist despite
the opposing force of competitive exclusion. Our results emphasize the role of
spatial assortment and physical forces on niche partitioning and the evolution
of diverse community structures.


  Ecological selection is a major driver of community assembly. Selection is
classified as stabilizing when species with intermediate trait values gain the
highest reproductive success, whereas selection is considered directional when
fitness is highest for species with extreme trait values. Previous studies have
investigated the effects of different selection types on trait distribution,
but the effects of selection on species diversity have remained unclear. Here,
we propose a framework for inferring the type and strength of selection by
studying species diversity and trait distribution together against null
expectations. We use a simulation model to confirm our prediction that
directional selection should lead to lower species diversity than stabilizing
selection despite a similar effect on trait community-weighted variance. We
apply the framework to a mesocosm system of annual plants to test whether
differences in species diversity between two habitats that vary in productivity
are related to differences in selection on seed mass. We show that, in both
habitats, species diversity was lower than the null expectation, but that
species diversity was lower in the more productive habitat. We attribute this
difference to strong directional selection for large-seeded species in the
productive habitat as indicated by trait community-weighted-mean being higher
and community-weighted variance being lower than the null expectations. In the
less productive habitat, we found that community-weighted variance was higher
than expected by chance, suggesting that seed mass could be a driver of niche
partitioning under such conditions. Altogether, our results suggest that
viewing species diversity and trait distribution as interrelated patterns
driven by the same process, ecological selection, is helpful in understanding
community assembly.


  In this paper, we have studied epidemiological models for Ebola infection
using nonlinear ordinary differential equations and optimal control theory. We
considered optimal control analysis of SIR and SEIR models for the deadly Ebola
infection using vaccination, treatment and educational campaign as
time-dependent controls functions. We have applied indirect methods to study
existing deterministic optimal control epidemic models for Ebola virus disease.
These methods in optimal control are based on Hamiltonian function and the
Pontryagin's maximum principle to construct adjoint equations and optimality
systems. The forward-backward sweep numerical scheme with fourth-order
Runge-Kutta method is used to solve the optimality system for the various
control strategies. From our numerical illustrations, we can conclude that,
effective educational campaigns and vaccination of susceptible individuals as
were as effective treatments of infected individuals can help reduce the
disease transmission.


  Background. Human aging is linked to many prevalent diseases. The aging
process is highly influenced by genetic factors. Hence, it is important to
identify human aging-related genes. We focus on supervised prediction of such
genes. Gene expression-based methods for this purpose study genes in isolation
from each other. While protein-protein interaction (PPI) network-based methods
for this purpose account for interactions between genes' protein products,
current PPI network data are context-unspecific, spanning different biological
conditions. Instead, here, we focus on an aging-specific subnetwork of the
entire PPI network, obtained by integrating aging-specific gene expression data
and PPI network data. The potential of such data integration has been
recognized but mostly in the context of cancer. So, we are the first to propose
a supervised learning framework for predicting aging-related genes from an
aging-specific PPI subnetwork.
  Results. In a systematic and comprehensive evaluation, we find that in many
of the evaluation tests: (i) using an aging-specific subnetwork indeed yields
more accurate aging-related gene predictions than using the entire network, and
(ii) predictive methods from our framework that have not previously been used
for supervised prediction of aging-related genes outperform existing prominent
methods for the same purpose.
  Conclusion. These results justify the need for our framework.


  A rearrangement operation makes a small graph-theoretical change to a
phylogenetic network to transform it into another one. For unrooted
phylogenetic trees and networks, popular rearrangement operations are tree
bisection and reconnection (TBR) and prune and regraft (PR) (called subtree
prune and regraft (SPR) on trees). Each of these operations induces a metric on
the sets of phylogenetic trees and networks. The TBR-distance between two
unrooted phylogenetic trees $T$ and $T'$ can be characterised by a maximum
agreement forest, that is, a forest with a minimum number of components that
covers both $T$ and $T'$ in a certain way. This characterisation has
facilitated the development of fixed-parameter tractable algorithms and
approximation algorithms. Here, we introduce maximum agreement graphs as a
generalisations of maximum agreement forests for phylogenetic networks. While
the agreement distance -- the metric induced by maximum agreement graphs --
does not characterise the TBR-distance of two networks, we show that it still
provides constant-factor bounds on the TBR-distance. We find similar results
for PR in terms of maximum endpoint agreement graphs.


  Diversity indices are useful single-number metrics for characterizing a
complex distribution of a set of attributes across a population of interest.
The utility of these different metrics or sets of metrics depend on the context
and application, and whether a predictive mechanistic model exists. In this
topical review, we first summarize the relevant mathematical principles
underlying heterogeneity in a large population before outlining the various
definitions of `diversity' and providing examples of scientific topics in which
its quantification plays an important role. We then review how diversity has
been a ubiquitous concept across multiple fields including ecology, immunology,
cellular barcoding experiments, and socioeconomic studies. Since many of these
applications involve sampling of populations, we also review how diversity in
small samples is related to the diversity in the entire population. Features
that arise in each of these applications are highlighted.


  Ingroup favoritism, the tendency to favor ingroup over outgroup, is often
explained as a product of intergroup conflict, or correlations between group
tags and behavior. Such accounts assume that group membership is meaningful,
whereas human data show that ingroup favoritism occurs even when it confers no
advantage and groups are transparently arbitrary. Another possibility is that
ingroup favoritism arises due to perceptual biases like outgroup homogeneity,
the tendency for humans to have greater difficulty distinguishing outgroup
members than ingroup ones. We present a prisoner's dilemma model, where
individuals use Bayesian inference to learn how likely others are to cooperate,
and then act rationally to maximize expected utility. We show that, when such
individuals exhibit outgroup homogeneity bias, ingroup favoritism between
arbitrary groups arises through direct reciprocity. However, this outcome may
be mitigated by: (1) raising the benefits of cooperation, (2) increasing
population diversity, and (3) imposing a more restrictive social structure.


  The last decade has witnessed a swiftly increasing interest in the design and
production of novel multivalent molecules as powerful alternatives for
conventional antibodies in the fight against cancer and infectious diseases.
However, while it is widely accepted that large-scale flexibility ($10-100$ nm)
and free/constrained dynamics (100 ns $- \mu$s) control the activity of such
novel molecules, computational strategies at the mesoscale still lag behind
experiments in optimizing the design of crucial features, such as the binding
cooperativity (a.k.a. avidity).
  In this study, we introduced different coarse-grained models of a
polymer-linked, two-nanobody composite molecule, with the aim of laying down
the physical bases of a thorough computational drug design protocol at the
mesoscale. We show that the calculation of suitable potentials of mean force
allows one to apprehend the nature, range and strength of the thermodynamic
forces that govern the motion of free and wall-tethered molecules. Furthermore,
we develop a simple computational strategy to quantify the
encounter/dissociation dynamics between the free end of a wall-tethered
molecule and the surface, at the roots of binding cooperativity. This procedure
allows one to pinpoint the role of internal flexibility and weak non-specific
interactions on the kinetic constants of the NB-wall encounter and
dissociation. Finally, we quantify the role and weight of rare events, which
are expected to play a major role in real-life situations, such as in the
immune synapse, where the binding kinetics is likely dominated by fluctuations.


  Resource competition theory predicts coexistence and exclusion patterns based
on species R*s, the minimum resource values required for a species to persist.
A central assumption of the theory is that all species have equal access to
resources. However, many systems are characterized by preemption exploitation,
where some species deplete resources before their competitors can access them
(e.g., asymmetric light competition, contest competition among animals). We
hypothesize that coexistence under preemption requires an R*-preemption
tradeoff, i.e., the species with the priority access should have a higher R*
(lower efficiency). Thus, we developed an extension of resource competition
theory to investigate partial and total preemption (in the latter, the
preemptor is unaffected by species with lower preemption rank). We found that
an R*-preemption tradeoff is a necessary condition for coexistence in all
models. Moreover, under total preemption, the tradeoff alone is sufficient for
coexistence. In contrast, under partial preemption, more conditions are needed,
which restricts the parameter space of coexistence. Finally, we discussed the
implications of our finding for seemingly distinct tradeoffs, which we view as
special cases of R*-preemption tradeoff. These tradeoffs include the
digger-grazer, the competition-colonization, and tradeoffs related to light
competition between trees and understories.


  Interlocus gene conversion (IGC) homogenizes paralogs. Little is known
regarding the mutation events that cause IGC and even less is known about the
IGC mutations that experience fixation. To disentangle the rates of fixed IGC
mutations from the tract lengths of these fixed mutations, we employ a
composite likelihood procedure. We characterize the procedure with simulations.
We apply the procedure to duplicated primate introns and to protein-coding
paralogs from both yeast and primates. Our estimates from protein-coding data
concerning the mean length of fixed IGC tracts were unexpectedly low and are
associated with high degrees of uncertainty. In contrast, our estimates from
the primate intron data had lengths in the general range expected from IGC
mutation studies. While it is challenging to separate the rate at which fixed
IGC mutations initiate from the average number of nucleotide positions that
these IGC events affect, all of our analyses indicate that IGC is responsible
for a substantial proportion of evolutionary change in duplicated regions. Our
results suggest that IGC should be considered whenever the evolution of
multigene families is examined.


  Phylogenetic trees are often constructed by using a metric on the set of taxa
that label the leaves of the tree. While there are a number of methods for
constructing a tree using a given metric, such trees will only display the
metric if it satisfies the so-called "four point condition", established by
Buneman in 1971. While this condition guarantees that a unique tree will
display the metric, meaning that the distance between any two leaves can be
found by adding the distances on arcs in the path between the leaves, it
doesn't exclude the possibility that a phylogenetic network might also display
the metric. This possibility was recently pointed out and "tree-metrized"
networks --- that display a tree metric --- with a single reticulation were
characterized. In this paper, we show that in the case of HGT (horizontal gene
transfer) networks, in fact there are tree-metrized networks containing many
reticulations.


  Focal adhesions are the loci of cellular adhesion to the extracellular
matrix. At these sites, various integrins forge connections between the
intracellular cytoskeleton and the outside world: large patches of multiple
types of integrins together grip hold of collagen, fibronectin and other
extracellular matrix components. The mixture of integrins composing the FA
will, in general, contain both slip bond integrins and catch bond
integrins---bonds whose lifetime increases with applied load and bonds for whom
it decreases when forced. Prior work suggests that catch bonds are essential
for proper FA stability and mechanosensory functionality. In the present work,
we investigate, numerically, the interplay between the two distinct types of
bonds and ask how the presence, in the same FA cluster, of slip bonds augments
the behavior of the catch bonds. We show, that mixing the two components m
low-force mechanical integrity, lacking in purely catch systems, while
preserving the potential to strengthen the FA bond by force as well as the
mechanosensory qualities of the catch bonds. We investigate the spatial
distribution in mixed-integrin FA's and show that the differential response to
loading leads, via an excluded volume interaction, to a dependence of the
individual integrin diffusivities on the applied load, an effect that has been
reported in experiments.


  We study the extinction risk of a fragmented population residing on a network
of patches coupled by migration, where the local patch dynamics include the
Allee effect. We show that mixing between patches dramatically influences the
population's viability. Slow migration is shown to always increase the
population's global extinction risk compared to the isolated case. At fast
migration, we demonstrate that synchrony between patches minimizes the
population's extinction risk. Moreover, we discover a critical migration rate
that maximizes the extinction risk of the population, and identify an
early-warning signal when approaching this state. Our theoretical results are
confirmed via the highly-efficient weighted ensemble method. Notably, our
analysis can also be applied to studying switching in gene regulatory networks
with multiple transcriptional states.


  Genomes in living organisms consist of the nucleotides adenine (A), guanine
(G), cytosine (C) and thymine (T). All prokaryotes have genomes consisting of
double-stranded DNA, where the A's and G's (purines) of one strand bind
respectively to the T's and C's (pyrimidines) of the other. As such, the number
of A's on one strand nearly equals the number of T's on the other, and the same
is true of one strand's G's and the other's C's. Globally, this relationship is
formalized as Chargaff's first parity rule; its strandwise equivalent is
Chargaff's second parity rule. Therefore, the GC content of any double-stranded
DNA genome can be expressed as %GC=100%-%AT. Variation in prokaryotic GC
content can be substantial between taxa but is generally small within microbial
genomes. This variation has been found to correlate with both phylogeny and
environmental factors. Since novel single-nucleotide polymorphisms (SNPs)
within genomes are at least partially linked to the environment, SNP GC content
can be considered a compound measure of an organism's environmental influences,
lifestyle and phylogeny. We present a mathematical model that describes how SNP
GC content in microbial genomes evolves over time as a function of the AT->GC
and GC->AT mutation rates with Gaussian white noise disturbances. The model
suggests that, in non-recombining bacteria, mutations can first accumulate
unnoticeably and then abruptly fluctuate out of control. Thus, minuscule
variations in mutation rates can suddenly become unsustainable, ultimately
driving a species to extinction if not counteracted early enough. This model,
which is suited specifically to symbiotic prokaryotes, conforms to scenarios
predicted by Muller's ratchet and may suggest that this is not always a
gradual, degrading process.


  Medical statistics reveal a significant dependence of hospitalized dengue
patient on the patient's age. To incorporate an age-dependence into a
mathematical model, we extend the classical ODE system of disease dynamics to a
PDE system. The equilibrium distribution is then determined by the fixed points
of resulting integro-differential equations. In this paper we use an extension
of the concept of the basic reproductive number to characterize parameter
regimes, where either only the disease-free or an endemic equilibrium exists.
Using rather general and minimal assumptions on the population distribution and
on the age-dependent transmission rate, we prove the existence of those
equilibria. Furthermore, we are able to prove the convergence of an iteration
scheme to compute the endemic equilibrium. To validate our model, we use
existing data from the city of Semarang, Indonesia for comparison and to
identify the model parameters.


  Deciphering gene regulatory networks is a central problem in computational
biology. Here, we explore the use of multi-modal neural networks to learn
predictive models of gene expression that include cis and trans regulatory
components. We learn models of stress response in the budding yeast
Saccharomyces cerevisiae. Our models achieve high performance and substantially
outperform other state-of-the-art methods such as boosting algorithms that use
pre-defined cis-regulatory features. Our model learns several cis and trans
regulators including well-known master stress response regulators. We use our
models to perform in-silico TF knock-out experiments and demonstrate that
in-silico predictions of target gene changes correlate with the results of the
corresponding TF knockout microarray experiment.


  We consider certain analytical features of a stochastic model that can
explain among other things competition among species and simultaneous predation
on the competing species from a geometric perspective which allows for a
systematic description of models admitting singular Lagrangians. The model
equations are shown to admit a Jacobi Last Multiplier which in turn allows for
the construction of a Lagrangian. The Lagrangian is of singular nature so that
construction of the Hamiltonian via a Legendre transformation is not possible.
A Hamiltonian description of the model therefore requires the introduction of
Dirac brackets. Explicit results are presented for the "Kill the winner" model
and its reductions.


  In systems biology, it is common to measure biochemical entities at different
levels of the same biological system. One of the central problems for the data
fusion of such data sets is the heterogeneity of the data. This thesis
discusses two types of heterogeneity. The first one is the type of data, such
as metabolomics, proteomics and RNAseq data in genomics. These different omics
data reflect the properties of the studied biological system from different
perspectives. The second one is the type of scale, which indicates the
measurements obtained at different scales, such as binary, ordinal, interval
and ratio-scaled variables. In this thesis, we developed several statistical
methods capable to fuse data sets of these two types of heterogeneity. The
advantages of the proposed methods in comparison with other approaches are
assessed using comprehensive simulations as well as the analysis of real
biological data sets.


  The physical chemistry of liquid-liquid phase separation (LLPS) of polymer
solutions bears directly on the assembly of biologically functional
droplet-like bodies from proteins and nucleic acids. These biomolecular
condensates include certain extracellular materials, and intracellular
compartments that are characterized as "membraneless organelles". Analytical
theories are a valuable, computationally efficient tool for addressing general
principles. LLPS of neutral homopolymers are quite well described by theory;
but it has been a challenge to develop general theories for the LLPS of
heteropolymers involving charge-charge interactions. Here we present a novel
theory that combines a random-phase-approximation treatment of polymer density
fluctuations and an account of intrachain conformational heterogeneity based
upon renormalized Kuhn lengths to provide predictions of LLPS properties as a
function of pH, salt, and charge patterning along the chain sequence. Advancing
beyond more limited analytical approaches, our LLPS theory is applicable to a
wide variety of charged sequences ranging from highly charged polyelectrolytes
to neutral or nearly neutral polyampholytes. The new theory should be useful in
high-throughput screening of protein and other sequences for their LLPS
propensities and can serve as a basis for more comprehensive theories that
incorporate non-electrostatic interactions. Experimental ramifications of our
theory are discussed.


  Ecologists have long suspected that species are more likely to interact if
their traits match in a particular way. For example, a pollination interaction
may be more likely if the proportions of a bee's tongue fit a plant's flower
shape. Empirical estimates of the importance of trait-matching for determining
species interactions, however, vary significantly among different types of
ecological networks. Here, we show that ambiguity among empirical
trait-matching studies may have arisen at least in parts from using overly
simple statistical models. Using simulated and real data, we contrast
conventional generalized linear models (GLM) with more flexible Machine
Learning (ML) models (Random Forest, Boosted Regression Trees, Deep Neural
Networks, Convolutional Neural Networks, Support Vector Machines, naive Bayes,
and k-Nearest-Neighbor), testing their ability to predict species interactions
based on traits, and infer trait combinations causally responsible for species
interactions. We find that the best ML models can successfully predict species
interactions in plant-pollinator networks, outperforming GLMs by a substantial
margin. Our results also demonstrate that ML models can better identify the
causally responsible trait-matching combinations than GLMs. In two case
studies, the best ML models successfully predicted species interactions in a
global plant-pollinator database and inferred ecologically plausible
trait-matching rules for a plant-hummingbird network, without any prior
assumptions. We conclude that flexible ML models offer many advantages over
traditional regression models for understanding interaction networks. We
anticipate that these results extrapolate to other ecological network types.
More generally, our results highlight the potential of machine learning and
artificial intelligence for inference in ecology, beyond standard tasks such as
image or pattern recognition.


  One of the most important challenges in the analysis of high-throughput
genetic data is the development of efficient computational methods to identify
statistically significant Single Nucleotide Polymorphisms (SNPs). Genome-wide
association studies (GWAS) use single-locus analysis where each SNP is
independently tested for association with phenotypes. The limitation with this
approach, however, is its inability to explain genetic variation in complex
diseases. Alternative approaches are required to model the intricate
relationships between SNPs. Our proposed approach extends GWAS by combining
deep learning stacked autoencoders (SAEs) and association rule mining (ARM) to
identify epistatic interactions between SNPs. Following traditional GWAS
quality control and association analysis, the most significant SNPs are
selected and used in the subsequent analysis to investigate epistasis. SAERMA
controls the classification results produced in the final fully connected
multi-layer feedforward artificial neural network (MLP) by manipulating the
interestingness measures, support and confidence, in the rule generation
process. The best classification results were achieved with 204 SNPs compressed
to 100 units (77% AUC, 77% SE, 68% SP, 53% Gini, logloss=0.58, and MSE=0.20),
although it was possible to achieve 73% AUC (77% SE, 63% SP, 45% Gini,
logloss=0.62, and MSE=0.21) with 50 hidden units - both supported by close
model interpretation.


  Synthetic biology is an interdisciplinary field aiming to design biochemical
systems with desired behaviors. To this end, molecular controllers have been
developed which, when embedded into a pre-existing ambient biochemical network,
control the dynamics of the underlying target molecular species. When
integrated into smaller compartments, such as biological cells in vivo, or
vesicles in vitro, controllers have to be calibrated to factor in the intrinsic
noise. In this context, molecular controllers put forward in the literature
have focused on manipulating the mean (first moment), and reducing the variance
(second moment), of the target species. However, many critical biochemical
processes are realized via higher-order moments, particularly the number and
configuration of the modes (maxima) of the probability distributions. To bridge
the gap, a controller called stochastic morpher is put forward in this paper,
inspired by gene-regulatory networks, which, under suitable time-scale
separations, morphs the probability distribution of the target species into a
desired predefined form. The morphing can be performed at the lower-resolution,
allowing one to achieve desired multi-modality/multi-stability, and at the
higher-resolution, allowing one to achieve arbitrary probability distributions.
Properties of the controller, such as robust perfect adaptation and
convergence, are rigorously established, and demonstrated on various examples.
Also proposed is a blueprint for an experimental implementation of stochastic
morpher.


  The World Health Organization describes zoonotic diseases as a major pandemic
threat, and modeling the behavior of such diseases is a key component of their
control. Many emerging zoonoses, such as SARS, Nipah, and Hendra, mutated from
their wild type while circulating in an intermediate host population, usually a
domestic species, to become more transmissible among humans, and moreover, this
transmission route will only become more likely as agriculture and trade
intensifies around the world. Passage through an intermediate host enables many
otherwise rare diseases to become better adapted to humans, and so
understanding this process with mathematical epidemiological models is
necessary to prevent epidemics of emerging zoonoses, guide policy interventions
in public health, and predict the behavior of an epidemic. In this paper, we
account for spillovers of a zoonotic disease mutating in an intermediate host
by means of modeling transmission dynamics within and between three host
species, namely, wild reservoir, intermediate domestic animals, and humans. We
calculate the basic reproductive number of the pathogen, present critical
conditions for the emergence dynamics of zoonosis, and perform stability
analysis of admissible disease equilibria. Our analytical results agree well
with long-term simulations of the system. We find that in the presence of
biologically realistic interspecies transmission parameters, a zoonotic disease
can establish itself in humans even if it fails to persist in its reservoir and
intermediate host species. Our model and results can be used to understand the
dynamic behavior of any zoonosis with intermediate hosts and assist efforts to
protect public health.


  Cell polarization requires redistribution of specific proteins to the nascent
front and back of a eukarytotic cell. Among these proteins are Rac and Rho,
members of the small GTPase family that regulate the actin cytoskeleton. Rac
promotes actin assembly and protrusion of the front edge, whereas Rho activates
myosin-driven contraction at the back. Mathematical models of cell polarization
at many levels of detail have appeared. One of the simplest based on
"wave-pinning", consists of a pair of reaction-diffusion equations for a single
GTPase. Mathematical analysis of wave-pinning so far is largely restricted to
static domains in one spatial dimension. Here we extend the analysis to cells
that change in size, showing that both shrinking and growing cells can lose
polarity. We further consider the feedback between mechanical tension, GTPase
activation, and cell deformation in both static, growing, shrinking, and moving
cells. Special cases (spatially uniform cell chemistry, absence or presence of
mechanical feedback) are analyzed, and the full model is explored by
simulations in 1D. We find a variety of novel behaviors, including
"dilution-induced" oscillations of Rac activity and cell size, as well as gain
or loss of polarization and motility in the model cell.


  Biochemical signal transduction, a form of molecular communication, can be
modeled using graphical Markov channels with input-modulated transition rates.
Such channel models are strongly non-Gaussian. In this paper we use a linear
noise approximation to construct a novel class of Gaussian additive white noise
channels that capture essential features of fully- and partially-observed
intensity-driven signal transduction. When channel state transitions that are
sensitive to the input signal are directly observable, high-frequency
information is transduced more efficiently than low-frequency information, and
the mutual information rate per bandwidth (spectral efficiency) is
significantly greater than when sensitive transitions and observable
transitions are disjoint. When both observable and hidden transitions are
input-sensitive, we observe a superadditive increase in spectral efficiency.


  The human gut microbiome is a complex ecosystem, in which hundreds of
microbial species and metabolites coexist, in part due to an extensive network
of cross-feeding interactions. However, both the large-scale trophic
organization of this ecosystem, and its effects on the underlying metabolic
flow, remain unexplored. Here, using a simplified model, we provide
quantitative support for a multi-level trophic organization of the human gut
microbiome, where microbes consume and secrete metabolites in multiple
iterative steps. Using a manually-curated set of metabolic interactions between
microbes, our model suggests about four trophic levels, each characterized by a
high level-to-level metabolic transfer of byproducts. It also quantitatively
predicts the typical metabolic environment of the gut (fecal metabolome) in
approximate agreement with the real data. To understand the consequences of
this trophic organization, we quantify the metabolic flow and biomass
distribution, and explore patterns of microbial and metabolic diversity in
different levels. The hierarchical trophic organization suggested by our model
can help mechanistically establish causal links between the abundances of
microbes and metabolites in the human gut.


  The cortical microtubule array of dark-grown hypocotyl cells of plant
seedlings undergoes a striking, and developmentally significant, reorientation
upon exposure to light. This process is driven by the exponential amplification
of a population of longitudinal microtubules, created by severing events
localized at crossovers with the microtubules of the pre-existing transverse
array. We present a dynamic one-dimensional model for microtubule amplification
through this type of templated severing. We focus on the role of the
probability of immediate rescue-after-severing of the newly-created lagging
microtubule, observed to be a characteristic feature of the reorientation
process. Employing stochastic simulations, we show that in the dynamic regime
of unbounded microtubule growth, a finite value of this probability is not
required for amplification to occur, but does strongly influence the degree of
amplification, and hence the speed of the reorientation process. In contrast,
in the regime of bounded microtubule growth, we show that amplification only
occurs above a critical threshold. We construct an approximate analytical
theory, based on a priori limiting the number of crossover events considered,
which allows us to predict the observed critical value of the
rescue-after-severing probability with reasonable accuracy.


  Healthy human prostate epithelial cells have the unique ability to produce
and secrete large amounts of citrate into the lumen of the prostate. Citrate is
a Krebs cycle metabolite produced in the condensation reaction between
acetyl-CoA and oxaloacetate in the mitochondria of the cell. With the
application of $^{13}$C enriched substrates, such as $^{13}$C glucose or
pyruvate, to prostate cells or tissues, it is possible to identify the
contributions of different metabolic pathways to this production and secretion
of citrate. In this work we present a quantitative model describing the
mitochondrial production and the secretion of citrate by prostatic epithelial
cells employing the $^{13}$C labeling pattern of secreted citrate as readout.
We derived equations for the secretion fraction of citrate and the contribution
of pyruvate dehydrogenase complex (PDC) versus the anaplerotic pyruvate
carboxylase (PC) pathways in supplying the Krebs cycle with carbons from
pyruvate for the production of citrate. These measures are independent of
initial $^{13}$C-enrichment of the administered supplements and of $^{13}$C
J-coupling patterns, making this method robust also if SNR is low. We propose
use of these equations to distinguish between citrate metabolism in healthy and
diseased prostate tissue, in particular upon malignant transformation.


  Virtual reality is a powerful tool with the ability to immerse a user within
a completely external environment. This immersion is particularly useful when
visualizing and analyzing interactions between small organic molecules,
molecular inorganic complexes, and biomolecular systems such as redox proteins
and enzymes. A common tool used in the biomedical community to analyze such
interactions is the APBS software, which was developed to solve the equations
of continuum electrostatics for large biomolecular assemblages. Numerous
applications exist for using APBS in the biomedical community including
analysis of protein ligand interactions and APBS has enjoyed widespread
adoption throughout the biomedical community. Currently, typical use of the
full APBS toolset is completed via the command line followed by visualization
using a variety of two-dimensional external molecular visualization software.
This process has inherent limitations: visualization of three-dimensional
objects using a two-dimensional interface masks important information within
the depth component. Herein, we have developed a single application,
UnityMol-APBS, that provides a dual experience where users can utilize the full
range of the APBS toolset, without the use of a command line interface, by use
of a simple \ac{GUI} for either a standard desktop or immersive virtual reality
experience.


  Open-source 3D data processing libraries originally developed for computer
vision and pattern recognition are used to align and compare molecular shapes
and sub-shapes. Here, a shape is represented by a set of points distributed on
the van der Waals surface of molecules. Each point is colored by its closest
atom, which itself belongs to a user defined class. The strength of this
representation is that it allows for comparisons of point clouds of different
kind of chemical entities: small molecules, peptides, proteins or cavities (the
negative image of the


  Quasi-steady state (QSS) reduction is a commonly used method to lower the
dimension of a differential equation model of a chemical reaction network. From
a mathematical perspective, QSS reduction is generally interpreted as a special
type of singular perturbation reduction, but in many instances the
correspondence is not worked out rigorously, and the QSS reduction may yield
incorrect results. The present paper contains a thorough discussion of QSS
reduction and its relation to singular perturbation reduction for the special,
but important, case when the right hand side of the differential equation is
linear in the variables to be eliminated. For this class we give necessary and
sufficient conditions for a singular perturbation reduction (in the sense of
Tikhonov and Fenichel) to exist, and to agree with QSS reduction. We then apply
the general results to chemical reaction networks with non-interacting species,
generalizing earlier results and methods for steady states to quasi-steady
state scenarios. We provide easy-to-check graphical conditions to select
parameter values yielding to singular perturbation reductions and additionally,
we identify a choice of parameters for which the corresponding singular
perturbation reduction agrees with the QSS reduction. Finally we consider a
number of examples.


  Many biological functions are executed by molecular machines, which consume
energy and convert it into mechanical work. Biological machines have evolved to
transport cargo, facilitate folding of proteins and RNA, remodel chromatin and
replicate DNA. A common aspect of these machines is that their functions are
driven out of equilibrium. It is a challenge to provide a general framework for
understanding the functions of biological machines, such as molecular motors,
molecular chaperones, and helicases. Using these machines as prototypical
examples, we describe a few general theoretical methods providing insights into
their functions. Although the theories rely on coarse-graining of these complex
systems they have proven useful in not only accounting for many in vitro
experiments but also addressing questions such as how the trade-off between
precision, energetic costs and optimal performances are balanced. However, many
complexities associated with biological machines will require one to go beyond
current theoretical methods. Simple point mutations in the enzyme could
drastically alter functions, making the motors bi-directional or result in
unexpected diseases or dramatically restrict the capacity of molecular
chaperones to help proteins fold. These examples are reminders that while the
search for principles of generality in biology is intellectually stimulating,
one also ought to keep in mind that molecular details must be accounted for to
develop a deeper understanding of processes driven by biological machines.
Going beyond generic descriptions of in vitro behavior to making genuine
understanding of in vivo functions will likely remain a major challenge for
some time to come. The combination of careful experiments and the use of
physical principles will be useful in elucidating the rules governing the
workings of biological machines.


  This paper presents an algebraic framework to study sign-sensitivities for
reaction networks modeled by means of systems of ordinary differential
equations. Specifically, we study the sign of the derivative of the
concentrations of the species in the network at steady state with respect to a
small perturbation on the parameter vector. We provide a closed formula for the
derivatives that accommodates common perturbations, and illustrate its form
with numerous examples. We argue that, mathematically, the study of the
response to the system with respect to changes in total amounts is not well
posed, and that one should rather consider perturbations with respect to the
initial conditions. We find a sign-based criterion to determine, without
computing the sensitivities, whether the sign depends on the steady state and
parameters of the system. This is based on earlier results of so-called
injective networks. Finally, we address systems with multiple steady states and
the restriction to stable steady states.


  The theory of elastic rods can be used to describe certain geometric and
topological properties of the DNA molecules. A similar effective field theory
approach was previously suggested to describe the conformations and dynamics of
proteins. In this letter we report a detailed study of the basic features of a
version of the proposed model, which assumes proteins to be very long
continuous curves. In the most appealing case, the model is based on a
potential with a pair of minima corresponding to helical and strand-like
configurations of the curves. It allows to derive several predictions about the
geometric features of the molecules, and we show that the predictions are
compatible with the phenomenology. While the helices represent the ground state
configurations, the abundance of beta strands is controlled by a parameter,
which can either completely suppress their presence in a molecule, or make them
abundant. The few-parameter model investigated in the letter rather represents
a universality class of protein molecules. Generalizations accounting for the
discrete nature and inhomogeneity of the molecules presumably allow to model
realistic cases.


  The article presents a 3D agent-based model to simulate the dynamics of
redistribution of bicoid mRNA. The model is developed on the basis of the
Skeledyne software [Odell and Foe, 2008]. The model reproduces active transport
by molecular motors over an undirected network of short microtubules. The model
is intended to simulate the data presented in [Spirov et al., 2009; Little et
al., 2011; Fahmy et al., 2014; Ali-Murthy and Kornberg, 2016; Cai, Spirov,
Akber et al., 2017].


  Varying environmental conditions affect relations between interacting
individuals in social dilemmas, thus affecting also the evolution of
cooperation. Oftentimes these environmental variations are seasonal and can
therefore be mathematically described as periodic changes. Accordingly, we here
study how periodic shifts between different manifestations of social dilemmas
affect cooperation. We observe a non-trivial interplay between the inherent
spatiotemporal dynamics that characterizes the spreading of cooperation in a
particular social dilemma type and the frequency of payoff changes. In
particular, we show that periodic changes between two available games with
global ordering best be fast, while periodic changes between global and local
ordering games best be slow for cooperation to thrive. We also show that the
frequency of periodic changes between two local ordering social dilemmas is
irrelevant, because then the process is fast and simply the average cooperation
level of the two is returned. The structure of the interaction network plays an
important role too in that lattices promote local ordering, whilst random
graphs hinder the formation of compact cooperative clusters. Conversely, for
local ordering the regular structure of the interaction network is only
marginally relevant as role-separating checkerboard patterns do not rely on
long-range order.


  We review current thinking about, and draw connections between, brain
energetics and metabolism, mitochondria and traumatic brain injury. In addition
to summarizing current thinking in these disciplines, our goal is to suggest a
framework for mechanisms and pathways based on optimal energetic decisions.


  Stochastic gene regulatory networks with bursting dynamics can be modeled
mesocopically as a generalized density-dependent Markov chain (GDDMC) or
macroscopically as a piecewise-deterministic Markov process (PDMP). Here we
prove a limit theorem showing that each family of GDDMCs will converge to a
PDMP as the system size tends to infinity. Moreover, under a simple dissipative
condition, we prove the existence and uniqueness of the stationary distribution
and the exponential ergodicity for the PDMP limit via the coupling method.
Further extensions and applications to single-cell stochastic gene expression
kinetics and bursty stochastic gene regulatory networks are also discussed and
the convergence of the stationary distribution of the GDDMC model to that of
the PDMP model is also proved.


  Here we investigate single-cell stochastic gene expression kinetics in a
minimal coupled gene circuit with positive-plus-negative feedback. A triphasic
stochastic bifurcation upon the increasing ratio of the positive and negative
feedback strengths is observed, which reveals a strong synergistic interaction
between positive and negative feedback loops. We discover that coupled
positive-plus-negative feedback amplifies gene expression mean but reduces gene
expression noise over a wide range of feedback strengths when promoter
switching is relatively slow, stabilizing gene expression around a relatively
high level. In addition, we study two types of macroscopic limits of the
discrete chemical master equation model: the Kurtz limit applies to proteins
with large burst frequencies and the L\'{e}vy limit applies to proteins with
large burst sizes. We derive the analytic steady-state distributions of the
protein abundance in a coupled gene circuit for both the discrete model and its
two macroscopic limits, generalizing the results obtained in [Chaos 26:043108,
2016]. We also obtain the analytic time-dependent protein distribution for the
classical Friedman-Cai-Xie random bursting model proposed in [Phys. Rev. Lett.
97:168302, 2006]. Our analytic results are further applied to study the
structure of gene expression noise in a coupled gene circuit and a complete
decomposition of noise in terms of five different biophysical origins is
provided.


  Molecule property prediction is a fundamental problem for computer-aided drug
discovery and materials science. Quantum-chemical simulations such as density
functional theory (DFT) have been widely used for calculating the molecule
properties, however, because of the heavy computational cost, it is difficult
to search a huge number of potential chemical compounds. Machine learning
methods for molecular modeling are attractive alternatives, however, the
development of expressive, accurate, and scalable graph neural networks for
learning molecular representations is still challenging. In this work, we
propose a simple and powerful graph neural networks for molecular property
prediction. We model a molecular as a directed complete graph in which each
atom has a spatial position, and introduce a recursive neural network with
simple gating function. We also feed input embeddings for every layers as skip
connections to accelerate the training. Experimental results show that our
model achieves the state-of-the-art performance on the standard benchmark
dataset for molecular property prediction.


  We report a neural architecture search framework, BioNAS, that is tailored
for biomedical researchers to easily build, evaluate, and uncover novel
knowledge from interpretable deep learning models. The introduction of
knowledge dissimilarity functions in BioNAS enables the joint optimization of
predictive power and biological knowledge through searching architectures in a
model space. By optimizing the consistency with existing knowledge, we
demonstrate that BioNAS optimal models reveal novel knowledge in both simulated
data and in real data of functional genomics. BioNAS provides a useful tool for
domain experts to inject their prior belief into automated machine learning and
therefore making deep learning easily accessible to practitioners. BioNAS is
available at https://github.com/zj-zhang/BioNAS-pub.


  In the self-assembly process which drives the formation of cellular
membranes, micelles, and capsids, a collection of separated subunits
spontaneously binds together to form functional and more ordered structures. In
this work, we study the statistical physics of self-assembly in a simpler
scenario: the formation of dimers from a system of monomers. The properties of
the model allow us to frame the microstate counting as a combinatorial problem
whose solution leads to an exact partition function. From the associated
equilibrium conditions, we find that such dimer systems come in two types:
"search-limited" and "combinatorics-limited", only the former of which has
states where partial assembly can be dominated by correct contacts. Using
estimates of biophysical quantities in systems of single-stranded DNA
dimerization, transcription factor and DNA interactions, and protein-protein
interactions, we find that all of these systems appear to be of the
search-limited type, i.e., their fully correct dimerization regimes are more
limited by the ability of monomers to find one another in the constituent
volume than by the combinatorial disadvantage of correct dimers. We derive the
parameter requirements for fully correct dimerization and find that rather than
the ratio of particle number and volume (number density) being the relevant
quantity, it is the product of particle diversity and volume that is
constrained. Ultimately, this work contributes to an understanding of
self-assembly by using the simple case of a system of dimers to analytically
study the combinatorics of assembly.