  A network is scale-free if its connectivity density function is proportional
to a power-law distribution. Scale-free networks may provide an explanation for
the robustness observed in certain physical and biological phenomena, since the
presence of a few highly connected hub nodes and a large number of small-
degree nodes may provide alternate paths between any two nodes on average --
such robustness has been suggested in studies of metabolic networks, gene
interaction networks and protein folding. A theoretical justification for why
biological networks are often found to be scale-free may lie in the well-known
fact that expanding networks in which new nodes are preferentially attached to
highly connected nodes tend to be scale-free. In this paper, we provide the
first efficient algorithm to compute the connectivity density function for the
ensemble of all secondary structures of a user-specified length, and show both
by computational and theoretical arguments that preferential attachment holds
when expanding the network from length n to length n + 1 structures. Since
existent power-law fitting software, such as powerlaw, cannot be used to
determine a power-law fit for our exponentially large RNA connectivity data, we
also implement efficient code to compute the maximum likelihood estimate for
the power-law scaling factor and associated Kolmogorov-Smirnov p-value.
Statistical goodness-of-fit tests indicate that one must reject the hypothesis
that RNA connectivity data follows a power-law distribution. Nevertheless, the
power-law fit is visually a good approximation for the tail of connectivity
data, and provides a rationale for investigation of preferential attachment in
the context of macromolecular folding.


  Ecological systems are emergent features of ecological and adaptive dynamics
of a community of interacting species. By natural selection through the abiotic
environment and by co-adaptation within the community, species evolve, thereby
giving rise to the ecological networks we regard as ecosystems. This
reductionist perspective can be contrasted with the view that as species have
to fit in the surrounding system, the system itself exerts selection pressure
on the evolutionary pathways of the species. This interplay of bottom-up and
top-down control in the development and growth of ecological systems has long
been discussed, however empirical ecosystem data is scarce and a comprehensive
mathematical framework is lacking. We present a way of quantifying the relative
weight of natural selection and coadaptation grounded in information theory, to
assess the relative role of bottom-up and top-down control in the evolution of
ecological systems, and analyse the information transfer in an individual based
stochastic complex systems model, the Tangled Nature Model of evolutionary
ecology. We show that ecological communities evolve from mainly bottom-up
controlled early-successional systems to more strongly top-down controlled
late-successional systems, as coadaptation progresses. Species which have a
high influence on selection are also generally more abundant. Hence our
findings imply that ecological communities are shaped by a dialogue of
bottom-up and top-down control, where the role of the systemic selection and
integrity becomes more pronounced the further the ecosystem is developed.


  The capture and translocation of biomolecules through nanometer-scale pores
are processes with a potential large number of applications, and hence they
have been intensively studied in the recent years. The aim of this paper is to
review existing models of the capture process by a nanopore, together with some
recent experimental data of short single- and double-stranded DNA captured by
Cytolysin A (ClyA) nanopore. ClyA is a transmembrane protein of bacterial
origin which has been recently engineered through site-specific mutations, to
allow the translocation of double- and single-stranded DNA. A comparison
between theoretical estimations and experiments suggests that for both cases
the capture is a reaction-limited process. This is corroborated by the observed
salt dependence of the capture rate, which we find to be in quantitative
agreement with the theoretical predictions.


  The nature of dispersal of many invasive pests and pathogens in agricultural
and forestry makes it necessary to consider how the actions of one manager
affect neighbouring properties. In addition to the direct effects of a
potential spread of a pest and the resulting economic loss, there are also
indirect consequences that affect whole regions and that require coordinated
actions to manage and/or to eradicate it (like movement restrictions). In this
paper we address the emergence and stability of cooperation among agents who
respond to a threat of an invasive pest or disease. The model, based on the
weakest-link paradigm, uses repeated multi-participant coordination games where
players' pay-offs depend on management decisions to prevent the invasion on
their own land as well as of their neighbours on a network. We show that for
the basic cooperation game agents select the risk-dominant strategy of a Stag
hunt game over the pay-off dominant strategy of implementing control measures.
However, cooperation can be achieved by the social planner offering a
biosecurity payment. The critical level of this payment depends on the details
of the decision-making process, with higher trust (based on a reputation of
other agents reflecting their past performance) allowing a significant
reduction in necessary payments and slowing down decay in cooperation when the
payment is low. We also find that allowing for uncertainty in decision-making
process can enhance cooperation for low levels of payments. Finally, we show
the importance of industry structure to the emergence of cooperation, with
increase in the average coordination number of network nodes leading to
increase in the critical biosecurity payment.


  The exon junction complex (EJC) deposited upstream of mRNA exon junctions
shapes structure, composition and fate of spliced mRNA ribonucleoprotein
particles (mRNPs). To achieve this, the EJC core nucleates assembly of a
dynamic shell of peripheral proteins that function in diverse
post-transcriptional processes. To illuminate consequences of EJC composition
change, we purified EJCs from human cells via peripheral proteins RNPS1 and
CASC3. We show that EJC originates as an SR-rich mega-dalton sized RNP that
contains RNPS1 but lacks CASC3. After mRNP export to the cytoplasm and before
translation, the EJC undergoes a remarkable compositional and structural
remodeling into an SR-devoid monomeric complex that contains CASC3.
Surprisingly, RNPS1 is important for nonsense-mediated mRNA decay (NMD) in
general whereas CASC3 is needed for NMD of only select mRNAs. The promotion of
switch to CASC3-EJC slows down NMD. Overall, the EJC compositional switch
dramatically alters mRNP structure and specifies two distinct phases of
EJC-dependent NMD.


  Extraordinarily large but short electric field pulses are reported by many
experiments to cause bipolar cancellation (BPC). This unusual cell response
occurs if a first pulse is followed by a second pulse with opposite polarity.
Possibly universal, BPC presently lacks a mechanistic explanation. Multiple
versions of the "standard model" of cell electroporation (EP) fail to account
for BPC. Here we show, for the first time, how an extension of the standard
model can account for a key experimental observation that essentially defines
BPC: the amount of a tracer that enters a cell, and how tracer influx can be
decreased by the second part of a bipolar pulse. The extended model can also
account for the recovery of BPC wherein the extent of BPC is diminished if the
spacing between the first and second pulses is increased. Our approach is
reverse engineering, meaning that we identify and introduce an additional
biophysical mechanism that allows pore transport to change. We hypothesize that
occluding molecules from outside the membrane enter or relocate within a pore.
Significantly, the additional mechanism is fundamental and general, involving a
combination of partitioning and hindrance. Molecules near the membrane can
enter pores to block transport of tracer molecules while still passing small
ions (+/- 1) that govern electrical behavior. Accounting for such behavior
requires an extension of the standard model.


  Various theoretical models have been proposed to understand the basic nature
of epidemics. Recent studies focus on the effects of mobility to epidemic
process. However, uncorrelated random walk is typically assumed as the type of
movement. In our daily life, the movement of people sometimes tends to be
limited to a certain direction, which can be described by biased random walk.
Here, we developed an agent-based model of susceptible-infected-recovered (SIR)
epidemic process in a 2D continuous space where agents tend to move in a
certain direction in addition to random movement. Moreover, we mainly focus on
the effect of the reduced mobility of infected agents. Our model assumes that,
when people are infected, their movement activity is greatly reduced because
they are physically weakened by the disease. By conducting extensive
simulations, we found that when the movement of infected people is limited, the
final epidemic size becomes small. However, that crucially depended on the
movement type of agents. Furthermore, the reduced mobility of infected agents
lengthened the duration of the epidemic because the infection progressed
slowly.


  Nonlinear systems driven by recurrent signals are known to exhibit complex
dynamical responses which, in the physiological context, can have important
functional consequences. One of the simplest biological systems that is exposed
to such repeated stimuli is the intra-cellular signaling network. In this paper
we investigate the periodic activation of an evolutionarily conserved motif of
this network, viz., the mitogen-activated protein kinase (MAPK) signaling
cascade, with a train of pulses. The resulting response of the cascade, which
shows integrative capability over several successive pulses, is characterized
by complex adaptive behavior. These include aspects of non-associative
learning, in particular, habituation and sensitization, which are observed in
response to high- and low-frequency stimulation, respectively. In addition, the
existence of a response threshold of the cascade, an apparent refractory
behavior following stimulation with short inter-pulse interval, and an
alternans-like response under certain conditions suggest an analogy with
excitable media.


  This chapter provides a pedagogical introduction and overview of spatial and
temporal correlation and fluctuation effects resulting from the fundamentally
stochastic kinetics underlying chemical reactions and the dynamics of
populations or epidemics. After reviewing the assumptions and mean-field type
approximations involved in the construction of chemical rate equations for
uniform reactant densities, we first discuss spatial clustering in birth-death
systems, where non-linearities are introduced through either density-limiting
pair reactions, or equivalently via local imposition of finite carrying
capacities. The competition of offspring production, death, and non-linear
inhibition induces a population extinction threshold, which represents a
non-equilibrium phase transition that separates active from absorbing states.
This continuous transition is characterized by the universal scaling exponents
of critical directed percolation clusters. Next we focus on the emergence of
depletion zones in single-species annihilation processes and spatial population
segregation with the associated reaction fronts in two-species pair
annihilation. These strong (anti-)correlation effects are dynamically generated
by the underlying stochastic kinetics. Finally, we address noise-induced and
fluctuation-stabilized spatio-temporal patterns in basic predator-prey systems,
exemplified by spreading activity fronts in the two-species Lotka-Volterra
model as well as spiral structures in the May-Leonard variant of cyclically
competing three-species systems akin to rock-paper-scissors games.


  Despite an explosion in the number of experimentally determined, atomically
detailed structures of biomolecules, many critical tasks in structural biology
remain data-limited. Whether performance in such tasks can be improved by using
large repositories of tangentially related structural data remains an open
question. To address this question, we focused on a central problem in biology:
predicting how proteins interact with one another---that is, which surfaces of
one protein bind to those of another protein. We built a training dataset, the
Database of Interacting Protein Structures (DIPS), that contains biases but is
two orders of magnitude larger than those used previously. We found that these
biases significantly degrade the performance of existing methods on
gold-standard data. Hypothesizing that assumptions baked into the hand-crafted
features on which these methods depend were the source of the problem, we
developed the first end-to-end learning model for protein interface prediction,
the Siamese Atomic Surfacelet Network (SASNet). Using only spatial coordinates
and identities of atoms, SASNet outperforms state-of-the-art methods trained on
gold-standard structural data, even when trained on only 3% of our new dataset.
Code and data available at https://github.com/drorlab/DIPS.


  The central goal of a dynamical theory of evolution is to abstract the mean
evolutionary trajectory in the trait space by considering ecological processes
at the level of the individual. In this work, we develop such a theory for a
new class of deterministic individual based models describing individual births
and deaths, which captures the essential features of standard stochastic
individual-based models and become identical with the latter under maximal
competition. The key motivation is to derive the canonical equation of adaptive
dynamics from this microscopic ecological model, which can be regarded as a
"toy model" for evolution, in a simple way and give it an intuitive geometric
interpretation. Another goal is to study evolution and sympatric speciation
under "maximal" competition. We show that these models, in the deterministic
limit of adaptive dynamics, lead to the same equations that describe the
unraveling of the mean evolutionary trajectory as those obtained from the
standard stochastic models. We further study conditions under which these
models lead to evolutionary branching and find them to be similar with those
obtained from the standard stochastic models. We find that though deterministic
models result in a strong competition that leads to a speed up in the temporal
dynamics of a population cloud in the phenotypic space as well as an increase
in the rate of generation of biodiversity, it does not seem to result in an
absolute increase in biodiversity as far as total number of species are
concerned. Hence, the "toy model" essentially captures all the features of the
standard stochastic model. Interestingly, the notion of a fitness function does
not explicitly enter in our derivation of the canonical equation, thereby
advocating a mechanistic view of evolution.


  Models in evolutionary game theory traditionally assume symmetric
interactions in homogeneous environments. Here, we consider populations
evolving in a heterogeneous environment, which consists of patches of different
qualities that are occupied by one individual each. The fitness of individuals
is not only determined by interactions with others but also by environmental
quality. This heterogeneity results in asymmetric interactions where the
characteristics of the interaction may depend on an individual's location.
Interestingly, in non-varying heterogeneous environments, the long-term
dynamics are the same as for symmetric interactions in an average, homogeneous
environment. However, introducing environmental feedback between an
individual's strategy and the quality of its patch results in rich
eco-evolutionary dynamics. Thus, individuals act as ecosystem engineers. The
nature of the feedback and the rate of ecological changes can relax or
aggravate social dilemmas and promote persistent periodic oscillations of
strategy abundance and environmental quality.


  Heart muscle contraction is normally activated by a synchronized Ca release
from sarcoplasmic reticulum (SR), a major intracellular Ca store. However,
under abnormal conditions Ca leaks from the SR, decreasing heart contraction
amplitude and increasing risk of life-threatening arrhythmia. The mechanisms
and regimes of SR operation generating the abnormal Ca leak remain unclear.
Here we employed both numerical and analytical modeling to get mechanistic
insights into the emergent Ca leak phenomenon. Our numerical simulations using
a detailed realistic model of Ca release unit (CRU) reveal sharp transitions
resulting in Ca leak. The emergence of leak is closely mapped mathematically to
the Ising model from statistical mechanics. The system steady-state behavior is
determined by two aggregate parameters: the analogues of magnetic field ($h$)
and the inverse temperature ($\beta$) in the Ising model, for which we have
explicit formulas in terms of SR Ca and release channel opening/closing rates.
The classification of leak regimes takes the shape of a phase $\beta$-$h$
diagram, with the regime boundaries occurring at $h$=0 and a critical value of
$\beta$ ($\beta*$) which we estimate using a classical Ising model and mean
field theory. Our theory predicts that a synchronized Ca leak will occur when
$h$>0 and $\beta>\beta*$ and a disordered leak occurs when $\beta<\beta*$ and
$h$ is not too negative. The disorder leak is distinguished from synchronized
leak (in long-lasting sparks) by larger Peierls contour lengths, an output
parameter reflecting degree of disorder. Thus, in addition to our detailed
numerical model approach we also offer an instantaneous computational tool
using analytical formulas of the Ising model for respective RyR parameters and
SR Ca load that describe and classify phase transitions and leak emergence.


  Based on statistical analysis of the complete genome sequences, a remote
relationship has been observed between the evolution of the genetic code and
the three domain tree of life. The existence of such a remote relationship need
to be explained. The unity of the living system throughout the history of life
relies on the common features of life: the homochirality, the genetic code and
the universal genome format. The universal genome format has been observed in
the genomic codon distributions as a common feature of life at the sequence
level. A main aim of this article is to reconstruct and to explain the
Phanerozoic biodiversity curve. It has been observed that the exponential
growth rate of the Phanerozoic biodiversity curve is about equal to the
exponential growth rate of genome size evolution. Hence it is strongly
indicated that the expansion of genomes causes the exponential trend of the
Phanerozoic biodiversity curve, where the conservative property during the
evolution of life is guaranteed by the universal genome format at the sequence
level. In addition, a consensus curve based on the climatic and eustatic data
is obtained to explain the fluctuations of the Phanerozoic biodiversity curve.
Thus, the reconstructed biodiversity curve based on genomic, climatic and
eustatic data agrees with Sepkoski's curve based on fossil data. The five mass
extinctions can be discerned in this reconstructed biodiversity curve, which
indicates a tectonic cause of the mass extinctions. And the declining
origination rate and extinction rate throughout the Phanerozoic eon might be
due to the growth trend in genome size evolution.


  With the progressive scale-down of semiconductor's feature size, people are
looking forward to More Moore and More than Moore. In order to offer a possible
alternative implementation process, people are trying to figure out a feasible
transfer from silicon to molecular computing. Such transfer lies on bio-based
modules programming with computer-like logic, aiming at realizing the Turing
machine. To accomplish this, the DNA-based combinational logic is inevitably
the first step we have taken care of. This timely overview paper introduces
combinational logic synthesized in DNA computing from both analog and digital
perspectives separately. State-of-the-art research progress is summarized for
interested readers to quick understand DNA computing, initiate discussion on
existing techniques and inspire innovation solutions. We hope this paper can
pave the way for the future DNA computing synthesis.


  Pareatic snakes possess outstanding asymmetry in the mandibular tooth number,
which has probably been caused by its evolution to improve the feeding on the
predominant dextral snails. Gene mutation can generate chiral inversion on the
snail body. A sinistral snail population can thrive in this ecological context.
The interactions between dextral/sinistral snails and Pareas snakes are modeled
in this paper by using a new generalized functional response of Holling type
II. Distinct Pareas species show different bilateral asymmetry degrees. This
parameter plays an essential role in our model and determines the evolution of
the populations. Stability of the solutions is also analyzed for different
regimes in the space of parameters.


  We study the optimal sustainable harvesting of a population that lives in a
random environment. The novelty of our setting is that we maximize the
asymptotic harvesting yield, both in an expected value and almost sure sense,
for a large class of harvesting strategies and unstructured population models.
We prove under relatively weak assumptions that there exists a unique optimal
harvesting strategy characterized by an optimal threshold below which the
population is maintained at all times by utilizing a local time push-type
policy. We also discuss, through Abelian limits, how our results are related to
the optimal harvesting strategies when one maximizes the expected cumulative
present value of the harvesting yield and establish a simple connection and
ordering between the values and optimal boundaries. Finally, we explicitly
characterize the optimal harvesting strategies in two different cases, one of
which is the celebrated stochastic Verhulst Pearl logistic model of population
growth.


  Studying how diverse human populations are related is of historical and
anthropological interest, in addition to providing a realistic null model for
testing for signatures of natural selection or disease associations.
Furthermore, understanding the demographic histories of other species is
playing an increasingly important role in conservation genetics. A number of
statistical methods have been developed to infer population demographic
histories using whole-genome sequence data, with recent advances focusing on
allowing for more flexible modeling choices, scaling to larger data sets, and
increasing statistical power. Here we review coalescent hidden Markov models, a
powerful class of population genetic inference methods that can effectively
utilize linkage disequilibrium information. We highlight recent advances, give
advice for practitioners, point out potential pitfalls, and present possible
future research directions.


  In permissive environments, E. coli can double its dry mass every 21 minutes.
During this time, ribosomes, RNA polymerases, and the proteome are all doubled.
Yet, the question of how to relate bacterial doubling time to other
biologically relevant time scales in the growth process remains illusive, due
to the complex temporal nesting pattern of these processes. In particular, the
relation between the cell's doubling time and the ribosome assembly time is not
known. Here we develop a model that connects growth rate to ribosome assembly
time and show that the existence of a self-assembly step increases the overall
growth rate, because during ribosome self-assembly existing ribosomes can start
a new round of reproduction, by making a new batch of ribosomal proteins prior
to the completion of the previous round. This overlapping of ribosome
reproduction cycles increases growth rate beyond the serial-limit that is
typically assumed to hold. Using recent data from ribosome profiling and well
known measurements of the average translation rate, rigid bounds on the in-vivo
ribosome self-assembly time are set, which are robust to the assumptions
regarding the biological noises involved. At 21 minutes doubling time, the
ribosome assembly time is found to be approximately 6 minutes --- three fold
larger than the common estimate. We further use our model to explain the
detrimental effect of a recently discovered ribosome assembly inhibitor drug,
and predict the effect of limiting the expression of ribosome assembly
chaperons on the overall growth rate.


  Persistence length of dsDNA is known to decrease with increase in ionic
concentration of the solution. In contrast to this, here we show that
persistence length of dsDNA increases dramatically as a function of ionic
liquid (IL) concentration. Using all atomic explicit solvent molecular dynamics
simulations and theoretical models we present, for the first time, a systematic
study to determine the mechanical properties of dsDNA in various hydrated ionic
liquids at different concentrations. We find that dsDNA in 50 wt% ILs have
lower persistence length and stretch modulus in comparison to 80 wt% ILs. We
further observe that both persistence length and stretch modulus of dsDNA
increase as we increase the ILs concentration. Present trend of stretch modulus
and persistence length of dsDNA with ILs concentration supports the predictions
of the macroscopic elastic theory, in contrast to the behavior exhibited by
dsDNA in monovalent salt. Our study further suggests the preferable ILs that
can be used for maintaining DNA stability during long-term storage.


  We formulate a multi-group and multi-vector epidemic model in which hosts'
dynamics is captured by staged-progression $SEIR$ framework and the dynamics of
vectors is captured by an $SI$ framework. The proposed model describes the
evolution of a class of zoonotic infections where the pathogen is shared by $m$
host species and transmitted by $p$ arthropod vector species. In each host, the
infectious period is structured into $n$ stages with a corresponding
infectiousness parameter to each vector species. We determine the basic
reproduction number $\mathcal{R}_0^2(m,n,p)$ and investigate the dynamics of
the systems when this threshold is less or greater than one. We show that the
dynamics of the multi-host, multi-stage, and multi-vector system is completely
determined by the basic reproduction number and the structure of the
host-vector network configuration. Particularly, we prove that the disease-free
\mbox{equilibrium} is globally asymptotically stable (GAS) whenever
$\mathcal{R}_0^2(m,n,p)<1$, and a unique strongly endemic equilibrium exists
and is GAS if $\mathcal{R}_0^2(m,n,p)>1$ and the host-vector configuration is
irreducible. That is, either the disease dies out or persists in all hosts and
all vector species.


  Microbial colonies cultured on agar Petri dishes have become a model system
to study biological evolution in populations expanding in space. Processes such
as clonal segregation and gene surfing have been shown to be affected by
interactions between microbial cells and their environment. In this work we
investigate the role of mechanical interactions such as cell-surface adhesion.
We compare two strains of the bacterium E. coli: a wild-type strain and a
"shaved" strain that adheres less to agar. We show that the shaved strain has a
selective advantage over the wild type: although both strains grow with the
same rate in liquid media, the shaved strain produces colonies that expand
faster on agar. This allows the shaved strain outgrow the wild type when both
strains compete for space. We hypothesise that, in contrast to a more common
scenario in which selective advantage results from increased growth rate, the
higher fitness of the shaved strain is caused by reduced adhesion and friction
with the agar surface.


  Bistability and multistationarity are properties of reaction networks linked
to switch-like responses and connected to cell memory and cell decision making.
Determining whether and when a network exhibits bistability is a hard and open
mathematical problem. One successful strategy consists of analyzing small
networks and deducing that some of the properties are preserved upon passage to
the full network. Motivated by this we study chemical reaction networks with
few chemical complexes. Under mass-action kinetics the steady states of these
networks are described by fewnomial systems, that is polynomial systems having
few distinct monomials. Such systems of polynomials are often studied in real
algebraic geometry by the use of Gale dual systems. Using this Gale duality we
give precise conditions in terms of the reaction rate constants for the number
and stability of the steady states of families of reaction networks with one
non-flow reaction.


  We derive an analytical expression for the transition path time (TPT)
distribution for a one-dimensional particle crossing a parabolic barrier. The
solution is expressed in terms of the eigenfunctions and eigenvalues of the
associated Fokker-Planck equation. The particle performs an anomalous dynamics
generated by a power-law memory kernel, which includes memoryless Markovian
dynamics as a limiting case. Our result takes into account absorbing boundary
conditions, extending existing results obtained for free boundaries. We show
that TPT distributions obtained from numerical simulations are in excellent
agreement with analytical results, while the typically employed free boundary
conditions lead to a systematic overestimation of the barrier height. These
findings may be useful in the analysis of experimental results on transition
path times. A web tool to perform this analysis is freely available.


  In the last decades, dispersal studies have benefited from the use of
molecular markers for detecting patterns differing between categories of
individuals and have highlighted sex-biased dispersal in several species. To
explain this phenomenon, several hypotheses implying mating systems,
intrasexual competition or sex-related handicaps have been proposed. In this
context, we investigated sex-biased dispersal in Armadillidium vulgare, a
terrestrial isopod with a promiscuous mating system. As a proxy for effective
dispersal, we performed a fine-scale investigation of the spatial genetic
structure in males and females, using individuals originating from five
sampling points located within 70 meters of each other. Based on microsatellite
markers and spatial autocorrelation analyses, our results revealed that while
males did not present a significant genetic structure at this geographic scale,
females were significantly and genetically more similar to each other when they
were collected in the same sampling point. As females invest more parental care
than males in A. vulgare, but also because this species is promiscuous and
males experience a high intrasexual competition, our results meet the
predictions of most classical hypotheses for sex-biased dispersal. We suggest
that widening dispersal studies to other isopods or crustaceans, differing in
their ecology or mating system and displaying varying levels of parental care,
might shed light on the processes underlying the evolution of sex-biased
dispersal.


  Protein aggregation including the formation of dimers and multimers in
solution, underlies an array of human diseases such as systemic amyloidosis
which is a fatal disease caused by misfolding of native globular proteins
damaging the structure and function of affected organs. Different kind of
interactors can interfere with the formation of protein dimers and multimers in
solution. A very special class of interactors are nanoparticles thanks to the
extremely efficient extension of their interaction surface. In particular
citrate-coated gold nanoparticles (cit-AuNPs) were recently investigated with
amyloidogenic protein $\beta$2-microglobulin ($\beta$2m). Here we present the
computational studies on two challenging models known for their enhanced
amyloidogenic propensity, namely $\Delta$N6 and D76N $\beta$2m naturally
occurring variants, and disclose the role of cit-AuNPs on their
fibrillogenesis. The proposed interaction mechanism lies in the interference of
the cit-AuNPs with the protein dimers at the early stages of aggregation, that
induces dimer disassembling. As a consequence, natural fibril formation can be
inhibited. Relying on the comparison between atomistic simulations at multiple
levels (enhanced sampling molecular dynamics and Brownian dynamics) and protein
structural characterisation by NMR, we demonstrate that the cit-AuNPs
interactors are able to inhibit protein dimer assembling. As a consequence, the
natural fibril formation is also inhibited, as found in experiment.


  Let Sn denote the network of all RNA secondary structures of length n, in
which undirected edges exist between structures s, t such that t is obtained
from s by the addition, removal or shift of a single base pair. Using
context-free grammars, generating functions and complex analysis, we show that
the asymptotic average degree is O(n) and that the asymptotic clustering
coeffcient is O(1/n), from which it follows that the family Sn, n = 1,2,3,...
of secondary structure networks is not small-world.


  Continuous cultures of mammalian cells are complex systems displaying
hallmark phenomena of nonlinear dynamics, such as multi-stability, hysteresis,
as well as sharp transitions between different metabolic states. In this
context mathematical models may suggest control strategies to steer the system
towards desired states. Although even clonal populations are known to exhibit
cell-to-cell variability, most of the currently studied models assume that the
population is homogeneous. To overcome this limitation, we use the maximum
entropy principle to model the phenotypic distribution of cells in a chemostat
as a function of the dilution rate. We consider the coupling between cell
metabolism and extracellular variables describing the state of the bioreactor
and take into account the impact of toxic byproduct accumulation on cell
viability. We present a formal solution for the stationary state of the
chemostat and show how to apply it in two examples. First, a simplified model
of cell metabolism where the exact solution is tractable, and then a
genome-scale metabolic network of the Chinese hamster ovary (CHO) cell line.
Along the way we discuss several consequences of heterogeneity, such as:
qualitative changes in the dynamical landscape of the system, increasing
concentrations of byproducts that vanish in the homogeneous case, and larger
population sizes.


  Understanding systems level behaviour of many interacting agents is
challenging in various ways, here we'll focus on the how the interaction
between components can lead to hierarchical structures with different types of
dynamics, or causations, at different levels. We use the Tangled Nature model
to discuss the co-evolutionary aspects connecting the microscopic level of the
individual to the macroscopic systems level. At the microscopic level the
individual agent may undergo evolutionary changes due to mutations of
strategies. The micro-dynamics always run at a constant rate. Nevertheless, the
system's level dynamics exhibit a completely different type of intermittent
abrupt dynamics where major upheavals keep throwing the system between
meta-stable configurations. These dramatic transitions are described by a
log-Poisson time statistics. The long time effect is a collectively adapted of
the ecological network. We discuss the ecological and macroevolutionary
consequences of the adaptive dynamics and briefly describe work using the
Tangled Nature framework to analyse problems in economics, sociology,
innovation and sustainability


  The metric space of phylogenetic trees defined by Billera, Holmes, and
Vogtmann, which we refer to as BHV space, provides a natural geometric setting
for describing collections of trees on the same set of taxa. However, it is
sometimes necessary to analyze collections of trees on non-identical taxa sets
(i.e., with different numbers of leaves), and in this context it is not evident
how to apply BHV space. Davidson et al. recently approached this problem by
describing a combinatorial algorithm extending tree topologies to regions in
higher dimensional tree spaces, so that one can quickly compute which
topologies contain a given tree as partial data. In this paper, we refine and
adapt their algorithm to work for metric trees to give a full characterization
of the subspace of extensions of a subtree. We describe how to apply our
algorithm to define and search a space of possible supertrees and, for a
collection of tree fragments with different leaf sets, to measure their
compatibility.


  Complementary DNA strands in solution reliably hybridize to form stable
duplexes. We study the kinetics of the hybridization process and the mechanisms
by which two initially isolated strands come together to form a stable double
helix. We adopt a multi-step computational approach. First, we perform a large
number of Brownian dynamics simulations of the hybridization process using the
coarse-grained oxDNA2 model. Second, we use these simulations to construct a
Markov State Model of DNA dynamics that uses a state decomposition based on the
inter-strand hydrogen bonding pattern. Third, we take advantage of Transition
Path Theory to obtain quantitative information about the thermodynamic and
dynamic properties of the hybridization process. We find that while there is a
large ensemble of possible hybridization pathways there is a single dominant
mechanism in which an initial base pair forms close to either end of the
nascent double helix, and the remaining bases pair sequentially in a
zipper-like fashion. We also show that the number of formed base pairs by
itself is insufficient to describe the transition state of the hybridization
process.


  Aim: The Akaike information Criterion (AIC) is widely used science to make
predictions about complex phenomena based on an entire set of models weighted
by Akaike weights. This approach (AIC model averaging; hereafter AvgAICc) is
often preferable than alternatives based on the selection of a single model.
Surprisingly, AvgAICc has not yet been introduced in ecological niche modeling
(ENM). We aimed to introduce AvgAICc in the context of ENM to serve both as an
optimality criterion in analyses that tune-up model parameters and as a
multi-model prediction strategy.
  Innovation: Results from the AvgAICc approach differed from those of
alternative approaches with respect to model complexity, contribution of
environmental variables, and predicted amount and geographic location of
suitable conditions for the focal species. Two theoretical properties of the
AvgAICc approach might justify that future studies will prefer its use over
alternative approaches: (1) it is not limited to make predictions based on a
single model, but it also uses secondary models that might have important
predictive power absent in a given single model favored by alternative
optimality criteria; (2) it balances goodness of fit and model accuracy, this
being of critical importance in applications of ENM that require model
transference.
  Main conclusions: Our introduction of the AvgAICc approach in ENM; its
theoretical properties, which are expected to confer advantages over
alternatives approaches; and the differences we found when comparing the
AvgAICc approach with alternative ones; should eventually lead to a wider use
of the AvgAICc approach. Our work should also promote further methodological
research comparing properties of the AvgAICc approach with respect to those of
alternative procedures.


  Motivation: High-throughput sequencing of large immune repertoires has
enabled the development of methods to predict the probability of generation by
V(D)J recombination of T- and B-cell receptors of any specific nucleotide
sequence. These generation probabilities are very non-homogeneous, ranging over
20 orders of magnitude in real repertoires. Since the function of a receptor
really depends on its protein sequence, it is important to be able to predict
this probability of generation at the amino acid level. However, brute-force
summation over all the nucleotide sequences with the correct amino acid
translation is computationally intractable. The purpose of this paper is to
present a solution to this problem.
  Results: We use dynamic programming to construct an efficient and flexible
algorithm, called OLGA (Optimized Likelihood estimate of immunoGlobulin
Amino-acid sequences), for calculating the probability of generating a given
CDR3 amino acid sequence or motif, with or without V/J restriction, as a result
of V(D)J recombination in B or T cells. We apply it to databases of
epitope-specific T-cell receptors to evaluate the probability that a typical
human subject will possess T cells responsive to specific disease-associated
epitopes. The model prediction shows an excellent agreement with published
data. We suggest that OLGA may be a useful tool to guide vaccine design.
  Availability: Source code is available at https://github.com/zsethna/OLGA


  The need to harmonise apparently irreconcilable arrangements in an ecosystem
--nestedness and segregation-- has triggered so far different strategies.
Methodological refinements, or the inclusion of behavioural preferences to the
network dynamics offer a limited approach to the problem, since one remains
constrained within a 2-dimensional view of an ecosystem, i.e. a
presence-absence matrix. Here we question this partial-view paradigm, and show
that nestedness and segregation may coexist across a varied range of scenarios.
To do so, we rely on an upscaled representation of an ecological community as
an $n$-partite hypergraph, inspired by Hutchinson's high-dimensional niche
concept and the latest trends on ecological multilayer networks. This yields an
inclusive description of an ecological system, for which we develop a natural
extension of the definition of nestedness to larger dimensional spaces,
revealing how competitive exclusion may operate regardless of a highly nested
bipartite system.


  An approach to estimate the influence of the treatment-type controls on the
basic reproduction number, R 0 , is proposed and elaborated. The presented
approach allows one to estimate the effect of a given treatment strategy or to
compare a number of different treatment strategies on the basic reproduction
number. All our results are valid for sufficiently small values of the control.
However, in many cases it is possible to extend this analysis to larger values
of the control as was illustrated by examples.


  Axons are linear processes of nerve cells that can range from a few tens of
micrometers up to meters in length. In addition to external cues, the length of
an axon is also regulated by unknown internal mechanisms. Molecular motors have
been suggested to generate oscillations with an axon length-dependent frequency
that could be used to measure an axon's extension. Here, we present a mechanism
that depends on the spectral decomposition of the oscillatory signal to
determine the axon length.


  The fitness landscape metaphor plays a central role on the modeling of
optimizing principles in many research fields, ranging from evolutionary
biology, where it was first introduced, to management research. Here we
consider the ensemble of trajectories of the imitative learning search, in
which agents exchange information on their fitness and imitate the fittest
agent in the population aiming at reaching the global maximum of the fitness
landscape. We assess the degree to which the starting and ending points
determine the learning trajectories using two measures, namely, the
predictability that yields the probability that two randomly chosen
trajectories are the same, and the mean path divergence that gauges the
dissimilarity between two learning trajectories. We find that the
predictability is greater in rugged landscapes than in smooth ones. The mean
path divergence, however, is strongly affected by the search parameters --
population size and imitation propensity -- that obliterate the influence of
the underlying landscape. The learning trajectories become more deterministic,
in the sense that there are fewer distinct trajectories and those trajectories
are more similar to each other, with increasing population size and imitation
propensity. In addition, we find that the roughness of the learning
trajectories, which measures the deviation from additivity of the fitness
function, is always greater than the roughness estimated over the entire
fitness landscape.


  It is desirable for detection and classification algorithms to generalize to
unfamiliar environments, but suitable benchmarks for quantitatively studying
this phenomenon are not yet available. We present a dataset designed to measure
recognition generalization to novel environments. The images in our dataset are
harvested from twenty camera traps deployed to monitor animal populations.
Camera traps are fixed at one location, hence the background changes little
across images; capture is triggered automatically, hence there is no human
bias. The challenge is learning recognition in a handful of locations, and
generalizing animal detection and classification to new locations where no
training data is available. In our experiments state-of-the-art algorithms show
excellent performance when tested at the same location where they were trained.
However, we find that generalization to new locations is poor, especially for
classification systems.


  Given a real sparse polynomial system, we present a general framework to find
explicit coefficients for which the system has more than one positive solution,
based on the recent article by Bihan, Santos and Spaenlehauer. We apply this
approach to find explicit reaction rate constants and total conservation
constants in biochemical reaction networks for which the associated dynamical
system is multistationary.


  Compartmental equations are primary tools in disease spreading studies. Their
predictions are accurate for large populations but disagree with empirical and
simulated data for finite populations, where uncertainties become a relevant
factor. Starting from the agent-based approach, we investigate the role of
uncertainties and autocorrelation functions in SIS epidemic model, including
their relationship with epidemiological variables. We find new differential
equations that take uncertainties into account. The findings provide improved
predictions to the SIS model and it can offer new insights for emerging
diseases.


  In this work we present new scalable, information theory-based variational
methods for the efficient model reduction of high-dimensional deterministic and
stochastic reaction networks. The proposed methodology combines, (a)
information theoretic tools for sensitivity analysis that allow us to identify
the proper coarse variables of the reaction network, with (b) variational
approximate inference methods for training a best-fit reduced model. This
approach takes advantage of both physicochemical modeling and data-based
approaches and allows to construct optimal parameterized reduced dynamics in
the number of variables, reactions and parameters, while controlling the
information loss due to the reduction. We demonstrate the effectiveness of our
model reduction method on several complex, high-dimensional chemical reaction
networks arising in biochemistry.


  Utilizing common resources is always a dilemma for community members. While
cooperator players restrain themselves and consider the proper state of
resources, defectors demand more than their supposed share for a higher payoff.
To avoid the tragedy of the common state, punishing the latter group seems to
be an adequate reaction. This conclusion, however, is less straightforward when
we acknowledge the fact that resources are finite and even a renewable resource
has limited growing capacity. To clarify the possible consequences, we consider
a coevolutionary model where beside the payoff-driven competition of cooperator
and defector players the level of a renewable resource depends sensitively on
the fraction of cooperators and the total consumption of all players. The
applied feedback-evolving game reveals that beside a delicately adjusted
punishment it is also fundamental that cooperators should pay special attention
to the growing capacity of renewable resources. Otherwise, even the usage of
tough punishment cannot save the community from an undesired end.


  Gene expression (GE) data capture valuable condition-specific information
("condition" can mean a biological process, disease stage, age, patient, etc.)
However, GE analyses ignore physical interactions between gene products, i.e.,
proteins. Since proteins function by interacting with each other, and since
biological networks (BNs) capture these interactions, BN analyses are
promising. However, current BN data fail to capture condition-specific
information. Recently, GE and BN data have been integrated using network
propagation (NP) to infer condition-specific BNs. However, existing NP-based
studies result in a static condition-specific subnetwork, even though cellular
processes are dynamic. A dynamic process of our interest is human aging. We use
prominent existing NP methods in a new task of inferring a dynamic rather than
static condition-specific (aging-related) subnetwork. Then, we study evolution
of network structure with age - we identify proteins whose network positions
significantly change with age and predict them as new aging-related candidates.
We validate the predictions via e.g., functional enrichment analyses and
literature search. Dynamic network inference via NP yields higher prediction
quality than the only existing method for inferring a dynamic aging-related BN,
which does not use NP.


  We present a computational method to evaluate the end-to-end and the contour
length distribution functions of short DNA molecules described by a mesoscopic
Hamiltonian. The method generates a large statistical ensemble of possible
configurations for each dimer in the sequence, selects the global equilibrium
twist conformation for the molecule and determines the average base pair
distances along the molecule backbone. Integrating over the base pair radial
and angular fluctuations, we derive the room temperature distribution functions
as a function of the sequence length. The obtained values for the most probable
end-to-end distance and contour length distance, providing a measure of the
global molecule size, are used to examine the DNA flexibility at short length
scales. It is found that, also in molecules with less than $\sim 60$ base
pairs, coiled configurations maintain a large statistical weight and,
consistently, the persistence lengths may be much smaller than in kilo-base
DNA.


  Similarly as in (Blancas et al. 2018) where nested coalescent processes are
studied, we generalize the definition of partition-valued homogeneous Markov
fragmentation processes to the setting of nested partitions, i.e. pairs of
partitions $(\zeta,\xi)$ where $\zeta$ is finer than $\xi$. As in the classical
univariate setting, under exchangeability and branching assumptions, we
characterize the jump measure of nested fragmentation processes, in terms of
erosion coefficients and dislocation measures. Among the possible jumps of a
nested fragmentation, three forms of erosion and two forms of dislocation are
identified - one of which being specific to the nested setting and relating to
a bivariate paintbox process.


  Artificial Intelligence techniques such as agent-based modeling and
probabilistic reasoning have shown promise in modeling complex biological
systems and testing ecological hypotheses through simulation. We develop an
agent-based model of Arctic foxes from Medniy Island while utilizing
Probabilistic Graphical Models to capture the conditional dependencies between
the random variables. Such models provide valuable insights in analyzing
factors behind catastrophic degradation of this population and in revealing
evolutionary mechanisms of its persistence in high-density environment. Using
empirical data from studies in Medniy Island, we create a realistic model of
Arctic foxes as agents, and study their survival and population dynamics under
a variety of conditions.


  We investigate the effects of solvent specificities on the stability of the
native structure (NS) of a protein on the basis of our free-energy function
(FEF). We use CPB-bromodomain (CBP-BD) and apoplastocyanin (apoPC) as
representatives of the protein universe and water, methanol, ethanol, and
cyclohexane as solvents. The NSs of CBP-BD and apoPC consist of 66$\%$
$\alpha$-helices and of 35$\%$ $\beta$-sheets and 4$\%$ $\alpha$-helices,
respectively. In order to assess the structural stability of a given protein
immersed in each solvent, we contrast the FEF of its NS against that of a
number of artificially created, misfolded decoys possessing the same amino-acid
sequence but significantly different topology and $\alpha$-helix and
$\beta$-sheet contents. We find that for both CBP-BD and apoPC, the energetic
component dominates in methanol, ethanol, and cyclohexane, with the most stable
structures in these solvents sharing the same characteristics described as an
association of a-helices. In water, the entropic component is as strong as or
even stronger than the energetic one, with a large gain of translational,
configurational entropy of water becoming crucially important, so that the
relative contents of $\alpha$-helix and $\beta$-sheet and the content of total
secondary structures are carefully selected to achieve sufficiently close
packing of side chains. Our analysis, which requires minimal computational
effort, can be applied to any protein immersed in any solvent and provides
robust predictions that are quite consistent with the experimental observations
for proteins in different solvent environments, thus paving the way toward a
more detailed understanding of the folding process.


  We analyze the movement of a starving forager on a one-dimensional periodic
lattice, where each location contains one unit of food. As the forager lands on
sites with food, it consumes the food, leaving the sites empty. If the forager
lands consecutively on $s$ empty sites, then it will starve. The forager has
two modes of movement: it can either diffuse, by moving with equal probability
to adjacent sites on the lattice, or it can jump to a uniformly randomly chosen
site on the lattice. We show that the lifetime $T$ of the forager in either
paradigm can be approximated by the sum of the cover time $\tau_{\rm cover}$
and the starvation time $s$, when $s$ far exceeds the number $n$ of lattice
sites. Our main findings focus on the hybrid model, where the forager has a
probability of either jumping or diffusing. The lifetime of the forager varies
non-monotonically according to $p_j$, the probability of jumping. By examining
a small system, analyzing a heuristic model, and using direct numerical
simulation, we explore the tradeoff between jumps and diffusion, and show that
the strategy that maximizes the forager lifetime is a mixture of both modes of
movement.


  Similarity of competitors has been proposed to facilitate coexistence of
species because it slows down competitive exclusion, thus making it easier for
equalizing mechanisms to maintain diverse communities. On the other hand, chaos
can promote coexistence of species. Here we link these two previously unrelated
findings, by analyzing the dynamics of food web models. We show that
near-neutrality of competition of prey, in the presence of predators, increases
the chance of developing chaotic dynamics. Moreover we confirm that this
results in a higher biodiversity. Our results suggest that near-neutrality may
promote biodiversity in two ways: through reducing the rates of competitive
displacement and through promoting non-equilibrium dynamics.


  In a neglected 2005 paper, Nobel Laureate Paul Lauterbur proposed that
molecular imprinting in amorphous materials -- a phenomenon with an extensive
experimental literature -- played a key role in abiogenesis. The present paper
builds on Lauterbur's idea to propose imprint-mediated templating (IMT), a
mechanism for prebiotic peptide replication that could potentially avoid a
range of difficulties arising in classic gene-first and metabolism-first models
of abiogenesis. Unlike models that propose prebiotic RNA synthesis, activation,
and polymerization based on unknown chemistries, peptide/IMT models are
compatible with demonstrably realistic prebiotic chemistries: synthesis of
dilute mixtures of racemic amino acids from atmospheric gases, and
polymerization of unactivated amino acids on hot, intermittently-wetted
surfaces. Starting from a peptide/IMT-based genetics, plausible processes could
support the elaboration of genetic and metabolic complexity in an early-Earth
environment, both explaining the emergence of homochirality and providing a
potential bridge to nucleic acid metabolism. Peptide/IMT models suggest
directions for both theoretical and experimental inquiry.


  Spectral clustering uses a graph Laplacian spectral embedding to enhance the
cluster structure of some data sets. When the embedding is one dimensional, it
can be used to sort the items (spectral ordering). A number of empirical
results also suggests that a multidimensional Laplacian embedding enhances the
latent ordering of the data, if any. This also extends to circular orderings, a
case where unidimensional embeddings fail. We tackle the task of retrieving
linear and circular orderings in a unifying framework, and show how a latent
ordering on the data translates into a filamentary structure on the Laplacian
embedding. We propose a method to recover it, illustrated with numerical
experiments on synthetic data and real DNA sequencing data. The code and
experiments are available at https://github.com/antrec/mdso.


  Computational procedures to foresee the 3D structure of aptamers are in
continuous progress. They constitute a crucial input to research, mainly when
the crystallographic counterpart of the structures in silico produced is not
present. At now, many codes are able to perform structure and binding
prediction, although their ability in scoring the results remains rather weak.
In this paper, we propose a novel procedure to complement the ranking outcomes
of free docking code, by applying it to a set of anti-angiopoietin aptamers,
whose performances are known. We rank the in silico produced configurations,
adopting a maximum likelihood estimate, based on their topological and
electrical properties. From the analysis, two principal kinds of conformers are
identified, whose ability to mimick the binding features of the natural
receptor is discussed. The procedure is easily generalizable to many biological
biomolecules, useful for increasing chances of success in designing
high-specificity biosensors (aptasensors).


  We provide a general mathematical framework based on the theory of graphical
models to study admixture graphs. Admixture graphs are used to describe the
ancestral relationships between past and present populations, allowing for
population merges and migration events, by means of gene flow. We give various
mathematical properties of admixture graphs with particular focus on properties
of the so-called $F$-statistics. Also the Wright-Fisher model is studied and a
general expression for the loss of heterozygosity is derived.


  We construct a one-bead-per-residue coarse-grained dynamical model to
describe intrinsically disordered proteins at significantly longer timescales
than in the all-atom models. In this model, inter-residue contacts form and
disappear during the course of the time evolution. The contacts may arise
between the sidechains, the backbones or the sidechains and backbones of the
interacting residues. The model yields results that are consistent with many
all-atom and experimental data on these systems. We demonstrate that the
geometrical properties of various homopeptides differ substantially in this
model. In particular, the average radius of gyration scales with the sequence
length in a residue-dependent manner.


  This paper outlines the implications of neural-level accounts of insight, and
models of the conceptual interactions that underlie creativity, for a theory of
cultural evolution. Since elements of human culture exhibit cumulative,
adaptive, open-ended change, it seems reasonable to view culture as an
evolutionary process, one fueled by creativity. Associative memory models of
creativity and mathematical models of how concepts combine and transform
through interaction with a context, support a view of creativity that is
incompatible with a Darwinian (selectionist) framework for cultural evolution,
but compatible with a non-Darwinian (Self-Other Reorganization) framework. A
theory of cultural evolution in which creativity is centre stage could provide
the kind of integrative framework for the behavioral sciences that Darwin
provided for the life sciences.


  Chase-Escape is a simple stochastic model that describes a predator-prey
interaction. In this model, there are two types of particles, red and blue. Red
particles colonize adjacent empty sites at an exponential rate $\lambda_{R}$,
whereas blue particles take over adjacent red sites at exponential rate
$\lambda_{B}$, but can never colonize empty sites directly. Numerical
simulations suggest that there is a critical value $p_{c}$ for the relative
growth rate $p:=\lambda_{R}/\lambda_{B}$. When $p<p_{c}$, mutual survival of
both types of particles has zero probability, and when $p>p_{c}$ mutual
survival occurs with positive probability. In particular, $p_{c} \approx 0.50$
for the square lattice case ($\mathbb Z^{2}$). Our simulations provide a
plausible explanation for the critical value. Near the critical value, the set
of occupied sites exhibits a fractal nature, and the hole sizes approximately
follow a power-law distribution.


  We consider cascades of enzymatic Goldbeter-Koshland loops with any number n
of layers, for which there exist two layers involving the same phosphatase.
Even if the number of variables and the number of conservation laws grow
linearly with n, we find explicit regions in reaction rate constant and total
conservation constant space for which the associated mass-action kinetics
dynamical system is multistationary. Our computations are based on the
theoretical results of our companion paper in arXiv:1807.05157, which are
inspired by results in real algebraic geometry by Bihan, Santos and
Spaenlehauer.


  Complex biological functions are carried out by the interaction of genes and
proteins. Uncovering the gene regulation network behind a function is one of
the central themes in biology. Typically, it involves extensive experiments of
genetics, biochemistry and molecular biology. In this paper, we show that much
of the inference task can be accomplished by a deep neural network (DNN), a
form of machine learning or artificial intelligence. Specifically, the DNN
learns from the dynamics of the gene expression. The learnt DNN behaves like an
accurate simulator of the system, on which one can perform in-silico
experiments to reveal the underlying gene network. We demonstrate the method
with two examples: biochemical adaptation and the gap-gene patterning in fruit
fly embryogenesis. In the first example, the DNN can successfully find the two
basic network motifs for adaptation - the negative feedback and the incoherent
feed-forward. In the second and much more complex example, the DNN can
accurately predict behaviors of essentially all the mutants. Furthermore, the
regulation network it uncovers is strikingly similar to the one inferred from
experiments. In doing so, we develop methods for deciphering the gene
regulation network hidden in the DNN "black box". Our interpretable DNN
approach should have broad applications in genotype-phenotype mapping.


  The growing number of large databases of animal tracking provides an
opportunity for analyses of movement patterns at the scales of populations and
even species. We used analytical approaches, developed to cope with big data,
that require no a priori assumptions about the behaviour of the target agents,
to analyse a pooled tracking dataset of 272 elephant seals (Mirounga leonina)
in the Southern Ocean, that was comprised of >500,000 location estimates
collected over more than a decade. Our analyses showed that the displacements
of these seals were described by a truncated power law distribution across
several spatial and temporal scales, with a clear signature of directed
movement. This pattern was evident when analysing the aggregated tracks despite
a wide diversity of individual trajectories. We also identified marine
provinces that described the migratory and foraging habitats of these seals.
Our analysis provides evidence for the presence of intrinsic drivers of
movement, such as memory, that cannot be detected using common models of
movement behaviour. These results highlight the potential for big data
techniques to provide new insights into movement behaviour when applied to
large datasets of animal tracking.


  We use optical trapping to continuously bend an isolated microtubule while
simultaneously measuring the applied force and the resulting filament strain,
thus allowing us to determine its elastic properties over a wide range of
applied strains. We find that, while in the low-strain regime, microtubules may
be quantitatively described in terms of the classical Euler-Bernoulli elastic
filament, above a critical strain they deviate from this simple elastic model,
showing a softening response with increasingdeformations. A three-dimensional
thin-shell model, in which the increased mechanical compliance is caused by
flattening and eventual buckling of the filament cross-section, captures this
softening effect in the high strain regime and yields quantitative values of
the effective mechanical properties of microtubules. Our results demonstrate
that properties of microtubules are highly dependent on the magnitude of the
applied strain and offer a new interpretation for the large variety in
microtubule mechanical data measured by different methods.


  Hypervariable T-cell receptors (TCR) play a key role in adaptive immunity,
recognising a vast diversity of pathogen-derived antigens. High throughput
sequencing of TCR repertoires (RepSeq) produces huge datasets of T-cell
receptor sequences from blood and tissue samples. However, our ability to
extract clinically relevant information from RepSeq data is limited, mainly
because little is known about TCR-disease associations. Here we present a
statistical approach called ALICE (Antigen-specific Lymphocyte Identification
by Clustering of Expanded sequences) that identifies TCR sequences that are
actively involved in the current immune response from a single RepSeq sample,
and apply it to repertoires of patients with a variety of disorders -
autoimmune disease (ankylosing spondylitis), patients under cancer
immunotherapy, or subject to an acute infection (live yellow fever vaccine).
The method's robustness is demonstrated by the agreement of its predictions
with independent assays, and is supported by its ability to selectively detect
responding TCR in the memory but not in the naive subset. ALICE requires no
longitudinal data collection nor large cohorts, and is thus directly applicable
to most RepSeq datasets. Its results facilitate the identification of TCR
variants associated with a wide variety of diseases and conditions, which can
be used for diagnostics, rational vaccine design and evaluation of the adaptive
immune system state.


  Functional protein-protein interactions are crucial in most cellular
processes. They enable multi-protein complexes to assemble and to remain
stable, and they allow signal transduction in various pathways. Functional
interactions between proteins result in coevolution between the interacting
partners, and thus in correlations between their sequences. Pairwise
maximum-entropy based models have enabled successful inference of pairs of
amino-acid residues that are in contact in the three-dimensional structure of
multi-protein complexes, starting from the correlations in the sequence data of
known interaction partners. Recently, algorithms inspired by these methods have
been developed to identify which proteins are functional interaction partners
among the paralogous proteins of two families, starting from sequence data
alone. Here, we demonstrate that a slightly higher performance for partner
identification can be reached by an approximate maximization of the mutual
information between the sequence alignments of the two protein families. Our
mutual information-based method also provides signatures of the existence of
interactions between protein families. These results stand in contrast with
structure prediction of proteins and of multi-protein complexes from sequence
data, where pairwise maximum-entropy based global statistical models
substantially improve performance compared to mutual information. Our findings
entail that the statistical dependences allowing interaction partner prediction
from sequence data are not restricted to the residue pairs that are in direct
contact at the interface between the partner proteins.


  In the present work, we develop a delayed Logistic growth model to study the
effects of decontamination on the bacterial population in the ambient
environment. Using the linear stability analysis, we study different case
scenarios, where bacterial population may establish at the positive equilibrium
or go extinct due to increased decontamination. The results are verified using
numerical simulation of the model.


  We consider a Moran model with recombination in a haploid population of size
$N$. At each birth event, with probability $1-\rho_N R$ the offspring copies
one parent's chromosome, and with probability $\rho_N R$ she inherits a
chromosome that is a mosaic of both parental chromosomes. We assume that at
time $0$ each individual has her chromosome painted in a different color and we
study the color partition of the chromosome that is asymptotically fixed in a
large population, when we look at a portion of the chromosome such that $\rho
:= \lim_{N\to \infty} \frac{\rho_N N}{2}\to \infty$. To do so, we follow
backwards in time the ancestry of the chromosome of a randomly sampled
individual. This yields a Markov process valued in the color partitions of the
half-line, that was introduced by Esser et al. (2016), in which blocks can
merge and split, called the partitioning process. Its stationary distribution
is closely related to the fixed chromosome in our Moran model with
recombination. We are able to provide an approximation of this stationary
distribution when $\rho \gg 1$ and an error bound. This allows us to show that
the distribution of the (renormalised) length of the leftmost block of the
partition (i.e. the region of the chromosome that carries the same color as 0)
converges to an exponential distribution. In addition, the geometry of this
block can be described in terms of a Poisson point process with an explicit
intensity measure.


  In many species, females are hypothesised to obtain 'good genes' for their
offspring by mating with males in good condition. However, female preferences
might deplete genetic variance and make choice redundant. Additionally,
high-condition males sometime produce low-fitness offspring, for example
because of environmental turnover and gene-by-environment interactions (GEIs)
for fitness, or because fit males carry sexually-antagonistic alleles causing
them to produce unfit daughters. Here, we extend previous theory by
investigating the evolution of female mate choice in a spatially explicit
evolutionary simulation implementing both GEIs and intralocus sexual conflict
(IASC), under sex-specific hard or soft selection. We show that IASC can weaken
female preferences for high-condition males or even cause a preference for
males in low condition, depending on the relative benefits of producing
well-adapted sons versus daughters, which in turn depends on the relative
hardness of selection on males and females. We discuss the relevance of our
results to conservation genetics and empirical evolutionary biology.


  The nested Kingman coalescent describes the dynamics of particles (called
genes) contained in larger components (called species), where pairs of species
coalesce at constant rate and pairs of genes coalesce at constant rate provided
they lie within the same species. We prove that starting from $rn$ species, the
empirical distribution of species masses (numbers of genes$/n$) at time $t/n$
converges as $n\to\infty$ to a solution of the deterministic
coagulation-transport equation $$
  \partial_t d \ = \ \partial_x ( \psi d ) \ + \ a(t)\left(d\star d - d
\right), $$ where $\psi(x) = cx^2$, $\star$ denotes convolution and $a(t)=
1/(t+\delta)$ with $\delta=2/r$. The most interesting case when $\delta =0$
corresponds to an infinite initial number of species. This equation describes
the evolution of the distribution of species of mass $x$, where pairs of
species can coalesce and each species' mass evolves like $\dot x = -\psi(x)$.
We provide two natural probabilistic solutions of the latter IPDE and address
in detail the case when $\delta=0$. The first solution is expressed in terms of
a branching particle system where particles carry masses behaving as
independent continuous-state branching processes. The second one is the law of
the solution to the following McKean-Vlasov equation $$ dx_t \ = \ - \psi(x_t)
\,dt \ + \ v_t\,\Delta J_t $$ where $J$ is an inhomogeneous Poisson process
with rate $1/(t+\delta)$ and $(v_t; t\geq0)$ is a sequence of independent rvs
such that $ {\mathcal L}(v_t) = {\mathcal L}(x_t)$. We show that there is a
unique solution to this equation and we construct this solution with the help
of a marked Brownian coalescent point process. When $\psi(x)=x^\gamma$, we show
the existence of a self-similar solution for the PDE which relates when
$\gamma=2$ to the speed of coming down from infinity of the nested Kingman
coalescent.


  We investigated frictional effects on the folding rates of a human telomerase
hairpin (hTR HP) and H-type pseudoknot from the Beet Western Yellow Virus (BWYV
PK) using simulations of the Three Interaction Site (TIS) model for RNA. The
heat capacity from TIS model simulations, calculated using temperature replica
exchange simulations, reproduces nearly quantitatively the available
experimental data for the hTR HP. The corresponding results for BWYV PK serve
as predictions. We calculated the folding rates ($k_\mathrm{F}$) from more than
100 folding trajectories for each value of the solvent viscosity ($\eta$) at a
fixed salt concentration of 200 mM. By using the theoretical estimate
($\propto$$\sqrt{N}$ where $N$ is the number of nucleotides) for folding free
energy barrier, $k_\mathrm{F}$ data for both the RNAs are quantitatively fit
using one-dimensional Kramers' theory with two parameters specifying the
curvatures in the unfolded basin and the barrier top. In the high-friction
regime ($\eta\gtrsim10^{-5}\,\textrm{Pa\ensuremath{\cdot}s}$), for both HP and
PK, $k_\mathrm{F}$s decrease as $1/\eta$ whereas in the low friction regime,
$k_\mathrm{F}$ values increase as $\eta$ increases, leading to a maximum
folding rate at a moderate viscosity
($\sim10^{-6}\,\textrm{Pa\ensuremath{\cdot}s}$), which is the Kramers turnover.
From the fits, we find that the speed limit to RNA folding at water viscosity
is between 1 and 4 $\mathrm{\mu s}$, which is in accord with our previous
theoretical prediction as well as results from several single molecule
experiments. Both the RNA constructs fold by parallel pathways. Surprisingly,
we find that the flux through the pathways could be altered by changing solvent
viscosity, a prediction that is more easily testable in RNA than in proteins.


  Are "turn-on" and "turn-off" functions in protein-protein interaction
networks exact opposites of each other? To answer this question, we implement a
minimal model for the evolution of functional protein-interaction networks
using a sequence-based mutational algorithm, and apply the model to study
neutral drift in networks that yield oscillatory dynamics. We study the roles
of activators and deactivators, two core components of oscillatory protein
interaction networks, and find a striking asymmetry in the roles of activating
and deactivating proteins, where activating proteins tend to be synergistic and
deactivating proteins tend to be competitive.


  An essential feature of the adaptive immune system is the proliferation of
antigen-specific lymphocytes during an immune reaction to form a large pool of
effector cells. This proliferation must be regulated to ensure an effective
response to infection while avoiding immunopathology. Recent experiments in
mice have demonstrated that the expansion of a specific clone of T cells in
response to cognate antigen obeys a striking inverse power law with respect to
the initial number of T cells. Here, we show that such a relationship arises
naturally from a model in which T cell expansion is limited by decaying levels
of presented antigen. The same model also accounts for the observed dependence
of T cell expansion on affinity for antigen and on the kinetics of antigen
administration. Extending the model to address expansion of multiple T cell
clones competing for antigen, we find that higher affinity clones can suppress
the proliferation of lower affinity clones, thereby promoting the specificity
of the response. Employing the model to derive optimal vaccination protocols,
we find that exponentially increasing antigen doses can achieve a nearly
optimized response. We thus conclude that the dynamics of presented antigen is
a key regulator of both the size and specificity of the adaptive immune
response.


  Adaptation mechanism of cells on the basis of stochastic gene expression and
epigenetic modification is proposed. From marginally stable states generated by
epigenetic modification, a gene expression pattern that achieves greater cell
growth is selected, as confirmed by simulations and analysis of several models.
The mechanism does not require any design of gene regulation networks and is
shown to be generic in a stochastic system with marginal stability. General
relevance of the mechanism to cell biology is also discussed.


  DNA Methylation has been the most extensively studied epigenetic mark.
Usually a change in the genotype, DNA sequence, leads to a change in the
phenotype, observable characteristics of the individual. But DNA methylation,
which happens in the context of CpG (cytosine and guanine bases linked by
phosphate backbone) dinucleotides, does not lead to a change in the original
DNA sequence but has the potential to change the phenotype. DNA methylation is
implicated in various biological processes and diseases including cancer. Hence
there is a strong interest in understanding the DNA methylation patterns across
various epigenetic related ailments in order to distinguish and diagnose the
type of disease in its early stages. In this work, the relationship between
methylated versus unmethylated CpG regions and cancer types is explored using
Convolutional Neural Networks (CNNs). A CNN based Deep Learning model that can
classify the cancer of a new DNA methylation profile based on the learning from
publicly available DNA methylation datasets is then proposed.


  Motivation: In this paper we present the latest release of EBIC, a
next-generation biclustering algorithm for mining genetic data. The major
contribution of this paper is adding support for big data, making it possible
to efficiently run large genomic data mining analyses. Additional enhancements
include integration with R and Bioconductor and an option to remove influence
of missing value on the final result.
  Results: EBIC was applied to datasets of different sizes, including a large
DNA methylation dataset with 436,444 rows. For the largest dataset we observed
over 6.6 fold speedup in computation time on a cluster of 8 GPUs compared to
running the method on a single GPU. This proves high scalability of the
algorithm.
  Availability: The latest version of EBIC could be downloaded from
http://github.com/EpistasisLab/ebic . Installation and usage instructions are
also available online.


  How fast a post-translational nascent protein escapes from the ribosomal exit
tunnel is relevant to its folding and protection against aggregation. Here, by
using Langevin molecular dynamics, we show that non-local native interactions
help decreasing the escape time, and foldable proteins generally escape much
faster than same-length self-repulsive homopolymers at low temperatures. The
escape process, however, is slowed down by the local interactions that
stabilize the {\alpha}-helices. The escape time is found to increase with both
the tunnel length and the concentration of macromolecular crowders outside the
tunnel. We show that a simple diffusion model described by the Smoluchowski
equation with an effective linear potential can be used to map out the escape
time distribution for various tunnel lengths and various crowder
concentrations. The consistency between the simulation data and the diffusion
model however is found only for the tunnel length smaller than a cross-over
length of 90 {\AA} to 110 {\AA}, above which the escape time increases much
faster with the tunnel length. It is suggested that the length of ribosomal
exit tunnel has been selected by evolution to facilitate both the efficient
folding and efficient escape of single domain proteins. We show that
macromolecular crowders lead to an increase of the escape time, and attractive
crowders are unfavorable for the folding of nascent polypeptide.


  The attractors of Boolean networks and their basins have been shown to be
highly relevant for model validation and predictive modelling, e.g., in systems
biology. Yet there are currently very few tools available that are able to
compute and visualise not only attractors but also their basins. In the realm
of asynchronous, non-deterministic modeling not only is the repertoire of
software even more limited, but also the formal notions for basins of
attraction are often lacking. In this setting, the difficulty both for theory
and computation arises from the fact that states may be ele- ments of several
distinct basins. In this paper we address this topic by partitioning the state
space into sets that are committed to the same attractors. These commitment
sets can easily be generalised to sets that are equivalent w.r.t. the long-term
behaviours of pre-selected nodes which leads us to the notions of markers and
phenotypes which we illustrate in a case study on bladder tumorigenesis. For
every concept we propose equivalent CTL model checking queries and an extension
of the state of the art model checking software NuSMV is made available that is
capa- ble of computing the respective sets. All notions are fully integrated as
three new modules in our Python package PyBoolNet, including functions for
visualising the basins, commitment sets and phenotypes as quotient graphs and
pie charts.


  Micromanipulations, perfusions and measurements performed using glass
microelectrodes filled with an electrolyte is a conventional technique for
experimental morphological and membrane electrophysiological studies at a
single cell and membrane surface level. The typical (effective) diameter of the
end of the glass microelectrode is from 500 up to less than 100 nm, which
prevents one from observing it using a standard optical microscope in
accordance with the optical resolution criteria, since the diameter less than
500 nm is indistinguishable within the interference zone. Microprocessor
programming of the puller (microforge) that provides pulling and tearing allows
to obtain in certain regimes the adjusted diameter and shape of the
micropipette tip, although this result is not fully controlled due to the above
limitations. In this connection it is necessary to design the control devices
for the micropipette tips both at the preparation and operation stages
(intracellular or extracellular insertion). This method also should provide
visualization of the processes occurring upon interaction of the microelectrode
tip with the cell in real time, depending on the electrode type and state,
which allows to level the artifacts arising with the systematic error frequency
from the uncontrolled operation of the micropipette tip after different ways of
the microelectrode filling with the electrolyte. We propose an installation
scheme that solves the above problems by means of introducing an
interferometric device for microscopic control of the microelectrode and
micromanipulator or microperfusor, for the first time for a given type of
optical instruments combined with the interferometric optical scheme.


  The prospects of application of graphene and related structures as the
membrane mimetic materials, capable of reproducing several biomembrane
functions up to the certain limit, are analyzed in the series of our papers.
This paper considers the possibility of the ion channel function modeling using
graphene and its derivatives. The physical mechanisms providing selective
permeability for different membrane mimetic materials, as well as the limits of
the adequate simulation of the transport, catalytic, sensing and electrogenic
properties of the cell membrane ion channels using bilayered graphene-based
structures are discussed.


  All known terrestrial proteins are coded as continuous strings of ~20 amino
acids. The patterns formed by the repetitions of elements in groups of finite
sequences describes the natural architectures of protein families. We present a
method to search for patterns and groupings of patterns in protein sequences
using a mathematically precise definition for 'repetition', an efficient
algorithmic implementation and a robust scoring system with no adjustable
parameters. We show that the sequence patterns can be well-separated into
disjoint classes according to their recurrence in nested structures. The
statistics of pattern occurrences indicate that short repetitions are enough to
account for the differences between natural families and randomized groups by
more than 10 standard deviations, while patterns shorter than 5 residues are
effectively random. A small subset of patterns is sufficient to account for a
robust ''familiarity'' definition of arbitrary sets of sequences.


  In this article, we use the illness-death model to present a mathematical
framework for studying the compression of morbidity (COM) hypothesis. It turns
out that questions about COM are completely determined by the transition rates
in the illness-death model and a closely related partial differential equation.
By this, the COM hypothesis is analytically tractable. To demonstrate the
usefulness of the mathematical framework, an example is given, which has been
motivated by empirical findings from Germany.


  The construction of a novel surgical instrument is considered, which is also
a probing device providing a signal to the measuring equipment, which after its
interpretation allows to obtain useful information about the section quality
and the biomaterial properties. We propose here some formalized considerations
on the possibility of its implementation for different variables registration.
The idea is also extrapolated into the field of micrurgy which refers to the
microelectrode techniques and the local potential registration in situ.


  Reconstituted nicotinic acetylcholine receptors (nAChRs) exhibit significant
gain-of-function upon addition of cholesterol to reconstitution mixtures, and
cholesterol affects organization of nAChRs within domain-forming membranes, but
whether nAChR partitions to cholesterol-rich liquid-ordered ("raft" or $l_o$)
domains or cholesterol-poor liquid-disordered ($l_{do}$) domains is unknown. We
use coarse-grained molecular dynamics simulations to observe spontaneous
interactions of cholesterol, saturated lipids, and polyunsaturated (PUFA)
lipids with nAChRs. In binary Dipalmitoylphosphatidylcholine:Cholesterol
(DPPC:CHOL) mixtures, both CHOL and DPPC acyl chains were observed
spontaneously entering deep "non-annular" cavities in the nAChR TMD,
particularly at the subunit interface and the $\beta$ subunit center,
facilitated by the low amino acid density in the cryo-EM structure of nAChR in
a native membrane. Cholesterol was highly enriched in the annulus around the
TMD, but this effect extended over (at most) 5-10\AA. In domain-forming ternary
mixtures containing PUFAs, the presence of a single receptor did not
significantly affect the likelihood of domain formation. nAChR partitioned to
any cholesterol-poor $l_{do}$ domain that was present, regardless of whether
the $l_{do}$ or $l_o$ domain lipids had PC or PE headgroups. Enrichment of
PUFAs among boundary lipids was positively correlated with their propensity for
demixing from cholesterol-rich phases. Long $n - 3$ chains (tested here with
Docosahexaenoic Acid, DHA) were highly enriched in annular and non-annular
embedded sites, partially displacing cholesterol and completely displacing
DPPC, and occupying sites even deeper within the bundle. Shorter $n - 6$ chains
were far less effective at displacing cholesterol from non-annular sites.


  The genetic instructions stored in the genome require an additional layer of
information to robustly determine cell fate. This additional regulation is
provided by the interplay between chromosome-patterning biochemical
("epigenetic") marks and three-dimensional genome folding. Yet, the physical
principles underlying the dynamical coupling between three-dimensional genomic
organisation and one-dimensional epigenetic patterns remain elusive. To shed
light on this issue, here we study by mean field theory and Brownian dynamics
simulations a magnetic polymer model for chromosomes, where each monomer
carries a dynamic epigenetic mark. At the single chromosome level, we show that
a first order transition describes the unlimited spreading of epigenetic marks,
a phenomenon that is often observed in vivo. At the level of the whole nucleus,
experiments suggest chromosomes form micro-phase separated compartments with
distinct epigenetic marks. We here discover that for a melt of magnetic
polymers such a morphology is thermodynamically unstable, but can be stabilised
by a non- equilibrium and ATP-mediated epigenetic switch between different
monomer states.


  Social and political polarization is a significant source of conflict and
poor governance in many societies. Thus, understanding its causes has become a
priority of scholars across many disciplines. Here we demonstrate that shifts
in socialization strategies analogous to political polarization and identity
politics can arise as a locally-beneficial response to both rising wealth
inequality and economic decline. Adopting a perspective of cultural evolution,
we develop a framework to study the emergence of polarization under shifting
economic environments. In many contexts, interacting with diverse out-groups
confers benefits from innovation and exploration greater than those that arise
from interacting exclusively with a homogeneous in-group. However, when the
economic environment favors risk-aversion, a strategy of seeking low-risk
interactions can be important to maintaining individual solvency. To capture
this dynamic, we assume that in-group interactions have a lower expected
outcome, but a more certain one. Thus in-group interactions are less risky than
out-group interactions. Our model shows that under conditions of economic
decline or increasing wealth inequality, some members of the population benefit
from adopting a risk-averse, in-group favoring strategy. Moreover, we show that
such in-group polarization can spread rapidly to the whole population and
persist even when the conditions that produced it have reversed. Finally we
offer empirical support for the role of income inequality as a driver of
affective polarization in the United States, mirroring findings on a panel of
developed democracies. Our work provides a framework for studying how disparate
forces interplay, via cultural evolution, to shape patterns of identity, and
unifies what are often seen as conflicting explanations for political
polarization: identity threat versus economic anxiety.


  Gene expression represents a fundamental interface between genes and
environment in the development and ongoing plasticity of the human organism.
Individual differences in gene expression are likely to underpin much of human
diversity, including psychiatric illness. Gene expression shows a distinct
regulatory pattern in different tissues. Therefore, brain tissue analysis
provides insights into brain disorder mechanisms. Furthermore, mechanistic
understanding of gene regulatory pattern can be provided through studying the
underlying relationships as a complex network. Identification of brain specific
gene relationships provides a complementary framework in which to tackle the
complex dysregulations that occur in neuropsychiatric and other neurological
disorders. Using a systems approach established in Mendelian randomization and
Bayesian Network, we integrated genetic and transcriptomic data from the
common-mind consortium and identified transcriptomic causal networks in
observational studies. Focusing on Schizophrenia disorder, we identified high
impact genes and revealed their underlying pathways in brain tissue. In
addition, we generated novel hypotheses including genes as causes of the
schizophrenia-associated genes and new genes associated with Schizophrenia.
This approach may facilitate a better understanding of the disease mechanism
that is complementary to molecular experimental studies especially for complex
systems and large-scale data sets.


  Time-varying network topologies can deeply influence dynamical processes
mediated by them. Memory effects in the pattern of interactions among
individuals are also known to affect how diffusive and spreading phenomena take
place. In this paper we analyze the combined effect of these two ingredients on
epidemic dynamics on networks. We study the susceptible-infected-susceptible
(SIS) and the susceptible-infected-removed (SIR) models on the recently
introduced activity-driven networks with memory. By means of an activity-based
mean-field approach we derive, in the long time limit, analytical predictions
for the epidemic threshold as a function of the parameters describing the
distribution of activities and the strength of the memory effects. Our results
show that memory reduces the threshold, which is the same for SIS and SIR
dynamics, therefore favouring epidemic spreading. The theoretical approach
perfectly agrees with numerical simulations in the long time asymptotic regime.
Strong aging effects are present in the preasymptotic regime and the epidemic
threshold is deeply affected by the starting time of the epidemics. We discuss
in detail the origin of the model-dependent preasymptotic corrections, whose
understanding could potentially allow for epidemic control on correlated
temporal networks.


  Dynamic models of signaling networks allow the formulation of hypotheses on
the topology and kinetic rate laws characterizing a given molecular network,
in-depth exploration and confrontation with kinetic biological data. Despite
its standardization, dynamic modeling of signaling networks still requires
successive technical steps that need to be carefully performed. Here, we detail
these steps by going through the mathematical and statistical framework. We
explain how it can be applied to the understanding of \beta-arrestin-dependent
signaling networks. We illustrate our methodology through the modeling of
\beta-arrestin recruitment kinetics at the Follicle Stimulating Hormone (FSH)
receptor supported by in-house Bioluminescence Resonance Energy Transfer (BRET)
data.


  Since the pioneering works of the founder of membrane mimetic chemistry Janos
H. Fendler it is known that a number of atomic or molecular clusters and films
(including nanoscale ones) are capable of mimicking the membrane functions.
Membrane mimetic materials can be either soft matter or solid state materials.
Conducting films (including those with magnetic properties) and semiconductors
are also known to possess membrane mimetic properties. If we consider the agent
exchange through the membrane in the operator form, the chemical composition of
the membranes and their models, as well as the difference between the atomic
and molecular clusters or layers become not so essential, and hence, membrane
mimetic chemistry of nano- and mesostructures do not differ significantly
within the agent-based approach. This invited review containing several parts
reflects the main aspects of the author's report at the conference "Graphene: a
molecule and 2D-crystal" (September 8-12, 2015, Novosibirsk, Russia) and
considers various aspects of the similarity between the graphene
nanostructures, membranes and bionic membrane-like nanomaterials.


  Recent research on the network modeling of complex systems has led to a
convenient representation of numerous natural, social, and engineered systems
that are now recognized as networks of interacting parts. Such systems can
exhibit a wealth of phenomena that not only cannot be anticipated from merely
examining their parts, as per the textbook definition of complexity, but also
challenge intuition even when considered in the context of what is now known in
network science. Here we review the recent literature on two major classes of
such phenomena that have far-reaching implications: (i) antagonistic responses
to changes of states or parameters and (ii) coexistence of seemingly
incongruous behaviors or properties -- both deriving from the collective and
inherently decentralized nature of the dynamics. They include effects as
diverse as negative compressibility in engineered materials, rescue
interactions in biological networks, negative resistance in fluid networks, and
the Braess paradox occurring across transport and supply networks. They also
include remote synchronization, chimera states and the converse of symmetry
breaking in brain, power-grid and oscillator networks as well as remote control
in biological and bio-inspired systems. By offering a unified view of these
various scenarios, we suggest that they are representative of a yet broader
class of unprecedented network phenomena that ought to be revealed and
explained by future research.