  Pituitary gonadotropins play an essential and pivotal role in the control of
human and animal reproduction within the hypothalamic-pituitary-gonadal (HPG)
axis. The computational modeling of pituitary gonadotropin signaling
encompasses phenomena of different natures such as the dynamic encoding of
gonadotropin secretion, and the intracellular cascades triggered by
gonadotropin binding to their cognate receptors, resulting in a variety of
biological outcomes. We overview historical and ongoing issues in modeling and
data analysis related to gonadotropin secretion in the field of both physiology
and neuro-endocrinology. We mention the different mathematical formalisms
involved, their interest and limits. We discuss open statistical questions in
signal analysis associated with key endocrine issues. We also review recent
advances in the modeling of the intracellular pathways activated by
gonadotropins, which yields promising development for innovative approaches in
drug discovery. The greatest challenge to be tackled in computational modeling
of pituitary gonadotropin signaling is the embedding of gonadotropin signaling
within its natural multi-scale environment, from the single cell level, to the
organic and whole HPG level. The development of modeling approaches of G
protein-coupled receptor signaling, together with multicellular systems biology
may lead to unexampled mechanistic understanding with critical expected
fallouts in the therapeutic management of reproduction.


  Biased signaling has recently emerged as an interesting mean to modulate the
function of many G protein-coupled receptors (GPCRs). Previous studies reported
two negative allosteric modulators (NAMs) of follicle-stimulating hormone
receptor (FSHR), ADX68692 and ADX68693, with differential effects on
FSHR-mediated steroidogenesis and ovulation. In this study, we attempted to
pharmacologically profile these NAMs on the closely related
luteinizing/chorionic gonadotropin hormone receptor (LH/CGR) with regards to
its canonical Gs/cAMP pathway as well as \beta-arrestin recruitment in HEK293
cells. The NAMs effects on progesterone and testosterone production were also
assessed in murine Leydig tumor cell line (mLTC-1). We found that both NAMs
strongly antagonized LH/CGR signaling in both HEK293 and mLTC-1 cells. ADX68693
appeared more potent than ADX68692 to inhibit hCG-induced cAMP and
\beta-arrestin 2 in HEK293 and mLTC-1 cells whereas no significant difference
in their efficacy on hCG-promoted \beta-arrestin 2 recruitment. Interestingly,
differential antagonism of the two NAMs on hCG-promoted steroidogenesis in
mLTC-1 cells was observed with significant inhibition of testosterone but not
progesterone production. These observations suggest biased effects of the two
NAMs on LH/CGR-dependent pathways controlling steroidogenesis, which appeared
to be different to that previously shown on FSHR. This also illustrates the
complexity of signaling pathways controlling FSHR- and LH/CGR-mediated
steroidogenesis, suggesting differential implication of cAMP and
\beta-arrestins. Together, our data demonstrate that ADX68692 and ADX68693 are
NAMs at the LH/CGR in addition to FSHR. These pharmacological characteristics
are important to consider for potential contraceptive and therapeutic
applications based on such compounds.


  This work focuses on the question of how identifiability of a mathematical
model, that is, whether parameters can be recovered from data, is related to
identifiability of its submodels. We look specifically at linear compartmental
models and investigate when identifiability is preserved after adding or
removing model components. In particular, we examine whether identifiability is
preserved when an input, output, edge, or leak is added or deleted. Our
approach, via differential algebra, is to analyze specific input-output
equations of a model and the Jacobian of the associated coefficient map. We
clarify a prior determinantal formula for these equations, and then use it to
prove that, under some hypotheses, a model's input-output equations can be
understood in terms of certain submodels we call "output-reachable". Our proofs
use algebraic and combinatorial techniques.


  In 2000, Elowitz and Leibler introduced the repressilator--a synthetic gene
circuit with three genes that cyclically repress transcription of the next
gene--as well as a corresponding mathematical model. Experimental data and
model simulations exhibited oscillations in the protein concentrations across
generations. In 2006, M\"{u}ller \textit{et al.}\ generalized the model to an
arbitrary number of genes and analyzed the resulting dynamics. Their new model
arose from five key assumptions, two of which are restrictive given current
biological knowledge. Accordingly, we propose a new repressilator system that
allows for general functions to model transcription, degradation, and
translation. We prove that, with an odd number of genes, the new model has a
unique steady state and the system converges to this steady state or to a
periodic orbit. We also give a necessary and sufficient condition for stability
of steady states when the number of genes is even and conjecture a condition
for stability for an odd number. Finally, we derive a new rate function
describing transcription that arises under more reasonable biological
assumptions than the widely used single-step binding assumption. With this new
transcription-rate function, we compare the model's amplitude and period with
that of a model with the conventional transcription-rate function. Taken
together, our results enhance our understanding of genetic regulation by
repression.


  Cellular heterogeneity is important to biological processes, including cancer
and development. However, proteome heterogeneity is largely unexplored because
of the limitations of existing methods for quantifying protein levels in single
cells. To alleviate these limitations, we developed Single Cell ProtEomics by
Mass Spectrometry (SCoPE-MS), and validated its ability to identify distinct
human cancer cell types based on their proteomes. We used SCoPE-MS to quantify
over a thousand proteins in differentiating mouse embryonic stem (ES) cells.
The single-cell proteomes enabled us to deconstruct cell populations and infer
protein abundance relationships. Comparison between single-cell proteomes and
transcriptomes indicated coordinated mRNA and protein covariation. Yet many
genes exhibited functionally concerted and distinct regulatory patterns at the
mRNA and the protein levels, suggesting that post-transcriptional regulatory
mechanisms contribute to proteome remodeling during lineage specification,
especially for developmental genes. SCoPE-MS is broadly applicable to measuring
proteome configurations of single cells and linking them to functional
phenotypes, such as cell type and differentiation potentials.


  The consumption rate is a process critically important for the stability of
consumer-resource systems and the persistence, sustainability and biodiversity
of complex food webs. Its mathematical description in the form of functional
response equations is a key problem for describing all trophic interactions.
Because some of the functional response models used in this study presented
redundancy between its parameters two methods were used to check for parameter
redundancy: the Hessian matrix calculation using automatic differentiation (AD)
which calculates derivatives numerically, but does not use finite differences,
and the symbolic method that calculates a derivative matrix and its rank. In
this work, we found that the models that better describe the functional
response of a rotifer is consumer dependant even at low consumer densities, but
their parameters can not be estimated simultaneously because parameter
redundancy. This means that fewer parameters or parameter combination can be
estimated than the original number of parameters in the models. Here, the model
parameters that incorporate intraspecific competition by interference in the
consumer-resource interaction are not identifiable, suggesting that this
problem may be more widespread than is generally appreciated in the literature
of food webs. Including knowledge on competitive interactions in current model
predictions will be a necessity in the next years for ecology as ecological
models are getting more complex and more real. Identifiability of biological
parameters in nonlinear ecological models will be an issue to consider.


  Families of ODE models characterized by a common sign structure of their
Jacobi matrix are investigated within the formalism of qualitative differential
equations. In the context of regulatory networks the sign structure of the
Jacobi matrix carries the information about which components of the network
inhibit or activate each other. Information about constraints on the behavior
of models in this family is stored in a so called qualitative state transition
graph. We showed previously that a similar approach can be used to analyze a
model pool of Boolean functions characterized by a common interaction graph.
Here we show that the opposite approach is fruitful as well. We show that the
qualitative state transition graph can be reduced to a "skeleton" represented
by a Boolean function conserving the reachability properties. This reduction
has the advantage that approaches such as model checking and network inference
methods can be applied to the "skeleton" within the framework of Boolean
networks. Furthermore, our work constitutes an alternative to approach to link
Boolean networks and differential equations.


  Invasion of new territories by alien organisms is of primary concern for
environmental and health agencies and has been a core topic in mathematical
modeling, in particular in the intents of reconstructing the past dynamics of
the alien organisms and predicting their future spatial extents. Partial
differential equations offer a rich and flexible modeling framework that has
been applied to a large number of invasions. In this article, we are
specifically interested in dating and localizing the introduction that led to
an invasion using mathematical modeling, post-introduction data and an adequate
statistical inference procedure. We adopt a mechanistic-statistical approach
grounded on a coupled reaction-diffusion-absorption model representing the
dynamics of an organism in an heterogeneous domain with respect to growth.
Initial conditions (including the date and site of the introduction) and model
parameters related to diffusion, reproduction and mortality are jointly
estimated in the Bayesian framework by using an adaptive importance sampling
algorithm. This framework is applied to the invasion of \textit{Xylella
fastidiosa}, a phytopathogenic bacterium detected in South Corsica in 2015,
France.


  Loss of pod dehiscence is a key step during soybean [Glycine max (L.) Merr.]
domestication. Genome-wide association analysis for soybean shattering
identified loci harboring Pdh1, NST1A and SHAT1-5. Pairwise epistatic
interactions were observed, and the dehiscent Pdh1 overcomes the resistance
conferred by NST1A or SHAT1-5 locus, indicating that Pdh1 predominates pod
dehiscence expression. Further candidate gene association analysis identified a
nonsense mutation in NST1A associated with pod dehiscence. Allele composition
and population differential analyses unraveled that Pdh1 and NST1A, but not
SHAT1-5, underwent domestication and modern breeding selections. Geographic
analysis showed that in Northeast China (NEC), indehiscence at both Pdh1 and
NST1A were required by cultivated soybean; while indehiscent Pdh1 alone is
capable of coping shattering in Huang-Huai-Hai (HHH) valleys where it
originated; and no specific indehiscence was required in Southern China (SC).
Geo-climatic investigation revealed strong correlation between relative
humidity and frequency of indehiscent Pdh1 across China. This study
demonstrates that the epistatic interaction between Pdh1 and NST1A fulfills a
pivotal role in determining the level of resistance against pod dehiscence.
Humidity shapes the distribution of indehiscent alleles. Our results also
suggest that HHH valleys, not NEC, was at least one of the origin centers of
cultivated soybean.


  The paper surveys the evolution of main algorithmic techniques to compare and
search biological sequences. We highlight key algorithmic ideas emerged in
response to several interconnected factors: shifts of biological analytical
paradigm, advent of new sequencing technologies, and a substantial increase in
size of the available data. We discuss the expansion of alignment-free
techniques coming to replace alignment-based algorithms in large-scale
analyses. We further emphasize recently emerged and growing applications of
sketching methods which support comparison of massive datasets, such as
metagenomics samples. Finally, we focus on the transition to population
genomics and outline associated algorithmic challenges.


  In combining DNA nanotechnology and high-bandwidth single-molecule detection
in nanopipettes, we demonstrate an all-electric, label-free hybridisation
sensor for short DNA sequences (< 100 nt). Such short fragments are known to
occur as circulating cell-free DNA in various bodily fluids, such as blood
plasma and saliva, and have been identified as disease markers for cancer and
infectious diseases. To this end, we use as a model system a 88-mer target from
the RV1910c gene in Mycobacterium tuberculosis that is associated with
antibiotic (isoniazid) resistance in TB. Upon binding to short probes attached
to long carrier DNA, we show that resistive pulse sensing in nanopipettes is
capable of identifying rather subtle structural differences, such as the
hybridisation state of the probes, in a statistically robust manner. With
significant potential towards multiplexing and high-throughput analysis, our
study points towards a new, single-molecule DNA assay technology that is fast,
easy to use and compatible with point of care environments.


  Many researches demonstrated that the DNA methylation, which occurs in the
context of a CpG, has strong correlation with diseases, including cancer. There
is a strong interest in analyzing the DNA methylation data to find how to
distinguish different subtypes of the tumor. However, the conventional
statistical methods are not suitable for analyzing the highly dimensional DNA
methylation data with bounded support. In order to explicitly capture the
properties of the data, we design a deep neural network, which composes of
several stacked binary restricted Boltzmann machines, to learn the low
dimensional deep features of the DNA methylation data. Experiments show these
features perform best in breast cancer DNA methylation data cluster analysis,
comparing with some state-of-the-art methods.


  The generalized Lotka-Volterra (gLV) equations, a classic model from
theoretical ecology, describe the population dynamics of a set of interacting
species. As the number of species in these systems grow in number, their
dynamics become increasingly complex and intractable. We introduce Steady State
Reduction (SSR), a method that reduces a gLV system of many ecological species
into two-dimensional (2D) subsystems that each obey gLV dynamics and whose
basis vectors are steady states of the high-dimensional model. We apply this
method to an experimentally-derived model of the gut microbiome in order to
observe the transition between "healthy" and "diseased" microbial states.
Specifically, we use SSR to investigate how fecal microbiota transplantation, a
promising clinical treatment for dysbiosis, can revert a diseased microbial
state to health.


  Follicle-stimulating hormone (FSH) and luteinizing hormone (LH) play
essential roles in animal reproduction. They exert their function through
binding to their cognate receptors, which belong to the large family of G
protein-coupled receptors (GPCRs). This recognition at the plasma membrane
triggers a plethora of cellular events, whose processing and integration
ultimately lead to an adapted biological response. Understanding the nature and
the kinetics of these events is essential for innovative approaches in drug
discovery. The study and manipulation of such complex systems requires the use
of computational modeling approaches combined with robust in vitro functional
assays for calibration and validation. Modeling brings a detailed understanding
of the system and can also be used to understand why existing drugs do not work
as well as expected, and how to design more efficient ones.


  Human luteinizing hormone (LH) and chorionic gonadotropin (hCG) have been
considered biologically equivalent because of their structural similarities and
their binding to the same receptor; the LH/CGR. However, accumulating evidence
suggest that LH/CGR differentially responds to the two hormones triggering
differential intracellular signaling and steroidogenesis. The mechanistic basis
of such differential responses remains mostly unknown. Here, we compared the
abilities of recombinant rhLH and rhCG to elicit cAMP, \beta-arrestin 2
activation, and steroidogenesis in HEK293 cells and mouse Leydig tumor cells
(mLTC-1). For this, BRET and FRET technologies were used allowing quantitative
analyses of hormone activities in real-time and in living cells. Our data
indicate that rhLH and rhCG differentially promote cell responses mediated by
LH/CGR revealing interesting divergences in their potencies, efficacies and
kinetics: rhCG was more potent than rhLH in both HEK293 and mLTC-1 cells.
Interestingly, partial effects of rhLH were found on \beta-arrestin recruitment
and on progesterone production compared to rhCG. Such a link was further
supported by knockdown experiments. These pharmacological differences
demonstrate that rhLH and rhCG act as natural biased agonists. The discovery of
novel mechanisms associated with gonadotropin-specific action may ultimately
help improve and personalize assisted reproduction technologies


  Follicle-stimulating hormone (FSH) is produced in the pituitary and is
essential for reproduction. It specifically binds to a membrane receptor (FSHR)
expressed in somatic cells of the gonads. The FSH/FSHR system presents many
peculiarities compared to classical G protein-coupled receptors (GPCRs). FSH is
a large naturally heterogeneous heterodimeric glycoprotein. The FSHR is
characterized by a very large NH2-terminal extracellular domain, which binds
the FSH and participates to the activation/inactivation switch of the receptor.
Once activated, the FSHR couples to G{\alpha}s and, in some instances, to other
G{\alpha} subunits. G protein-coupled receptor kinases and \beta-arrestins are
also recruited to the FSHR and account for its desensitization, the control of
its trafficking and its intracellular signalling. Of note, the FSHR
internalization and recycling are very fast and involve very early endosomes
instead of early endosomes. All the transduction mechanisms triggered upon FSH
stimulation lead to the activation of a complex signalling network that
controls gene expression by acting at multiple levels. The integration of these
mechanisms leads to context-adapted responses from the target gonadal cells,
but also indirectly affects the fate of germ cells. Depending of the
physiological/developmental stage, FSH elicits proliferation, differentiation
or apoptosis in order to maintain the homeostasis of the reproductive system.
Pharmacological tools targeting FSHR recently came to the fore and open
promising prospects both for basic research and therapeutic applications. This
paper provides an updated review of the most salient aspects and peculiarities
of FSHR biology and pharmacology.


  Protein-DNA interactions are critical for the successful functioning of all
natural systems. The key role in these interactions is played by processes of
protein search for specific sites on DNA. Although it has been studied for many
years, only recently microscopic aspects of these processes became more clear.
In this work, we present a review on current theoretical understanding of the
molecular mechanisms of the protein target search. A comprehensive
discrete-state stochastic method to explain the dynamics of the protein search
phenomena is introduced and explained. Our theoretical approach utilizes a
first-passage analysis and it takes into account the most relevant
physical-chemical processes. It is able to describe many fascinating features
of the protein search, including unusually high effective association rates,
high selectivity and specificity, and the robustness in the presence of
crowders and sequence heterogeneity.


  In this article we study the stabilizing of a primitive pattern of behaviour
for the two-species community with chemotaxis due to the short-wavelength
external signal. We use a system of Patlak-Keller-Segel type as a model of the
community. It is well-known that such systems can produce complex unsteady
patterns of behaviour which are usually explained mathematically by
bifurcations of some basic solutions that describe simpler patterns. As far as
we aware, all such bifurcations in the models of the Patlak-Keller-Segel type
had been found for homogeneous (i.e. translationally invariant) systems where
the basic solutions are equilibria with homogeneous distributions of all
species. The model considered in the present paper does not possess the
translational invariance: one of species (the predators) is assumed to be
capable of moving in response to a signal produced externally in addition to
the signal emitted by another species (the prey). For instance, the external
signal may arise from the inhomogeneity of the distribution of an environmental
characteristic such as temperature, salinity, terrain relief, etc. Our goal is
to examine the effect of short-wavelength inhomogeneity. To do this, we employ
a certain homogenization procedure. We separate the short-wavelength and smooth
components of the system response and derive a slow system governing the latter
one. Analysing the slow system and comparing it with the case of homogeneous
environment shows that, generically, a short-wavelength inhomogeneity results
in an exponential decrease in the motility of the predators. The loss of
motility prevents, to a great extent, the occurrence of complex unsteady
patterns and dramatically stabilizes the primitive basic solution. In some
sense, the necessity of dealing with intensive small-scale changes of the
environment makes the system unable to respond to other challenges.


  Oxidative stress, a reaction caused by the imbalance between the reactive
oxygen species of human organism and its ability to detoxify reactive
intermediates and to repair the resulting damage plays an important role in
HIV-infections. On one hand, HIV infection is responsible for the chronic
oxidative stress of the patients. On the other hand, the oxidative stress
contributions to the HIV disease pathogenesis. In this paper, we integrate
oxidative stress into an HIV infection model to investigate its effects on the
virus dynamics. Through mathematical analysis, we obtain the basic reproduction
number R0 of the model which describes the persistence of viruses. In
particular, we show that for R0 > 1, the model has a bistable interval with
virus rebound threshold and elite control threshold. Numerical simulations and
bifurcation analysis are presented to illustrate the viral dynamics under
oxidative stress. Our investigation reveals the interplay between viruses and
the reaction of human organism including immune response and oxidative stress,
and their effects on the health of human being.


  RNA forms elaborate secondary structures through intramolecular base pairing.
These structures perform critical biological functions within each cell. Due to
the availability of a polynomial algorithm to calculate the partition function
over these structures, they are also a suitable system for the statistical
physics of disordered systems. In this model, below the denaturation
temperature, random RNA secondary structures exist in one of two phases: a
strongly disordered, low-temperature glass phase, and a weakly disordered,
high-temperature molten phase. The probability of two bases to pair decays with
their distance with an exponent 3/2 in the molten phase, and about 4/3 in the
glass phase. Inspired by previous results from a renormalized field theory of
the glass transition separating the two phases, we numerically study this
transition. We introduce distinct order parameters for each phase, that both
vanish at the critical point. We finally explore the driving mechanism behind
this transition.


  Groupthink occurs when everyone in a group starts thinking alike, as when
people put unlimited faith in a leader. Avoiding this phenomenon is a
ubiquitous challenge to problem-solving enterprises and typical countermeasures
involve the mobility of group members. Here we use an agent-based model of
imitative learning to study the influence of the mobility of the agents on the
time they require to find the global maxima of NK-fitness landscapes. The
agents cooperate by exchanging information on their fitness and use this
information to copy the fittest agent in their influence neighborhoods, which
are determined by face-to-face interaction networks. The influence
neighborhoods are variable since the agents perform random walks in a
two-dimensional space. We find that mobility is slightly harmful for solving
easy problems, i.e. problems that do not exhibit suboptimal solutions or local
maxima. For difficult problems, however, mobility can prevent the imitative
search being trapped in suboptimal solutions and guarantees a better
performance than the independent search for any system size.


  Ecosystems are often under threat by invasive species which, through their
invasion dynamics, create ecological networks to spread. We present preliminary
results using a technique of GIS coupled with complex network analysis to model
the movement and spread of Lantana Camara in Rajaji Tiger Reserve, India, where
prey species are being affected because of habitat degradation due to Lantana
invasion. Understanding spatio-temporal aspects of the spread mechanism are
essential for better management in the region. The objective of the present
study is to develop insight into some key characteristics of the regulatory
mechanism for lantana spread inside RTR. Lantana mapping was carried out by
field observations along multiple transects and plots and the data generated
was used as input for MaxEnt modelling to identify land patches in the study
area that are favourable for lantana growth. The patch information so obtained
is integrated with a raster map generated by identifying different
topographical features in the study area which are favourable for lantana
growth. The integrated data is analysed with a complex network perspective,
where relatively dense potential lantana distribution patches are considered as
vertices, connected by relatively sparse lantana continuities, identified as
edges. The network centrality analysis reveal key patches in the study area
that play specialized roles in the spread of lantana in a large region. Hubs in
the lantana network are primarily identified as dry seasonal river beds and
their management is proposed as a vital strategy to contain lantana invasion.
The lantana network is found to exhibit small-world architecture with a well
formed community structure. We infer that the above properties of the lantana
network have major contribution in regulating the rapid infestation and even
spread of the plant through the entire region of study.


  Finding the largest few principal components of a matrix of genetic data is a
common task in genome-wide association studies (GWASs), both for dimensionality
reduction and for identifying unwanted factors of variation. We describe a
simple random matrix model for matrices that arise in GWASs, showing that the
singular values have a bulk behavior that obeys a Marchenko-Pastur distributed
with a handful of large outliers. We also implement Golub-Kahan-Lanczos (GKL)
bidiagonalization in the Julia programming language, providing thick restarting
and a choice between full and partial reorthogonalization strategies to control
numerical roundoff. Our implementation of GKL bidiagonalization is up to 36
times faster than software tools used commonly in genomics data analysis for
computing principal components, such as EIGENSOFT and FlashPCA, which use dense
LAPACK routines and randomized subspace iteration respectively.


  It has recently been suggested that the fundamental haploid-diploid cycle of
eukaryotic sex exploits a rudimentary form of the Baldwin effect. Thereafter
the other associated phenomena can be explained as evolution tuning the amount
and frequency of learning experienced by an organism. Using the well-known NK
model of fitness landscapes it is here shown that the emergence of sex
determination systems can also be explained under this view of eukaryotic
evolution.


  Direct Coupling Analysis (DCA) is a now widely used method to leverage
statistical information from many similar biological systems to draw meaningful
conclusions on each system separately. DCA has been applied with great success
to sequences of homologous proteins, and also more recently to whole-genome
population-wide sequencing data. We here argue that the use of DCA on the
genome scale is contingent on fundamental issues of population genetics. DCA
can be expected to yield meaningful results when a population is in the
Quasi-Linkage Equilibrium (QLE) phase studied by Kimura and others, but not,
for instance, in a phase of Clonal Competition. We discuss how the exponential
(Potts model) distributions emerge in QLE, and compare couplings to
correlations obtained in a study of about 3,000 genomes of the human pathogen
Streptococcus pneumoniae.


  We study models of biological evolution and investigate a key factor to yield
self-organized criticality (SOC). The Bak-Sneppen (BS) model is the most basic
model that shows an SOC state, which is developed based on minimal and
plausible assumptions of Darwinian competition. Another class of models, which
have population dynamics and simple rules for species migrations, has also been
studied. It turns out that they do not show an SOC state although the
assumptions made in these models are similar to those in the BS model. To
clarify the origin of these differences and to identify a key ingredient of
SOC, we study models that bridge the BS model and the Dynamical Graph model,
which is a representative of the population dynamics models. From a comparative
study of the models, we find that SOC is found when the fluctuations of the
number of species $N$ are suppressed, while it shows off-critical states when
$N$ changes according to its evolutionary dynamics. This indicates that the
assumption of the fixed system size in the BS model plays a pivotal role to
drive the system into an SOC state, and casts doubt on its applicability to
actual evolutionary dynamics.


  Protein-peptide interactions play essential roles in many cellular processes
and their structural characterization is the major focus of current
experimental and theoretical research. Two decades ago, it was proposed to
employ the steered molecular dynamics to assess the strength of protein-peptide
interactions. The idea behind using steered molecular dynamics simulations is
that the mechanical stability can be used as a promising and an efficient
alternative to computationally highly demanding estimation of binding affinity.
However, mechanical stability defined as a peak in force-extension profile
depends on the choice of the pulling direction. Here we propose an uncommon
choice of the pulling direction along resultant dipole moment vector, which has
not been explored in simulations so far. Using explicit solvent all-atom MD
simulations, we apply steered molecular dynamics technique to probe mechanical
resistance of protein-peptide system pulled along two different vectors. A
novel pulling direction, along the resultant dipole moment vector, results in
stronger forces compared to commonly used peptide unbinding along center of
masses vector. Our results demonstrate that resultant dipole moment is one of
the factors influencing the mechanical stability of protein-peptide complex.


  We study existence and stability of stationary solutions of a system of
semilinear parabolic partial differential equations that occurs in population
genetics. It describes the evolution of gamete frequencies in a geographically
structured population of migrating individuals in a bounded habitat. Fitness of
individuals is determined additively by two recombining, diallelic genetic loci
that are subject to spatially varying selection. Migration is modeled by
diffusion. Of most interest are spatially non-constant stationary solutions,
so-called clines. In a two-locus cline all four gametes are present in the
population, i.e., it is an internal stationary solution. We provide conditions
for existence and linear stability of a two-locus cline if recombination is
either sufficiently weak or sufficiently strong relative to selection and
diffusion. For strong recombination, we also prove uniqueness and global
asymptotic stability. For arbitrary recombination, we determine the stability
properties of the monomorphic equilibria, which represent fixation of a single
gamete.


  The objective of a genome-wide association study (GWAS) is to associate
subsequences of individuals' genomes to the observable characteristics called
phenotypes (e.g., high blood pressure). Motivated by the GWAS problem, in this
paper we introduce the information-theoretic problem of \emph{associated
subsequence retrieval}, where a dataset of $N$ (possibly high-dimensional)
sequences of length $G$, and their corresponding observable (binary)
characteristics is given. The sequences are chosen independently and uniformly
at random from $\mathcal{X}^G$, where $\mathcal{X}$ is a finite alphabet. The
observable (binary) characteristic is only related to a specific unknown
subsequence of length $L$ of the sequences, called \textit{associated
subsequence}. For each sequence, if the associated subsequence of it belongs to
a universal finite set, then it is more likely to display the observable
characteristic (i.e., it is more likely that the observable characteristic is
one). The goal is to retrieve the associated subsequence using a dataset of $N$
sequences and their observable characteristics. We demonstrate that as the
parameters $N$, $G$, and $L$ grow, a threshold effect appears in the curve of
probability of error versus the rate which is defined as ${Gh(L/G)}/{N}$, where
$h(\cdot)$ is the binary entropy function. This effect allows us to define the
capacity of associated subsequence retrieval. We develop an achievable scheme
and a matching converse for this problem, and thus characterize its capacity in
two scenarios: the zero-error-rate and the $\epsilon$-error-rate.


  The driving force for active physical and biological systems is determined by
both the underlying landscape and the non-equilibrium curl flux. While
landscape can be quantified in the experiments by the histograms of the
collecting trajectories of the observables, the experimental flux
quantification is still challenging. In this work, we studied the single
molecule enzyme dynamics and observed the deviation in kinetics from the
conventional Michaelis-Menton reaction rate. We identified and quantified the
non-equilibrium flux as the origin of such non-Michaelis-Menton enzyme rate
behavior. This is the first time of rigorous quantification of the flux for the
driving force of the non-equilibrium active dynamics. We also quantified the
corresponding non-equilibrium thermodynamics in terms of chemical potential and
entropy production. We identified and quantified the origin of the flux,
chemical potential and entropy production as the heat absorbed (energy input)
in the enzyme reaction.


  Plants can sense both intracellular and extracellular mechanical forces and
can respond through morphological changes. The signaling components responsible
for mechanotransduction of the touch response are largely unknown. Here, we
performed a high-throughput SILIA (stable isotope labeling in
Arabidopsis)-based quantitative phosphoproteomics analysis to profile changes
in protein phosphorylation resulting from 40 seconds of force stimulation in
Arabidopsis thaliana. Of the 24 touch-responsive phosphopeptides identified,
many were derived from kinases, phosphatases, cytoskeleton proteins, membrane
proteins and ion transporters. TOUCH-REGULATED PHOSPHOPROTEIN1 (TREPH1) and MAP
KINASE KINASE 2 (MKK2) and/or MKK1 became rapidly phosphorylated in
touch-stimulated plants. Both TREPH1 and MKK2 are required for touch-induced
delayed flowering, a major component of thigmomorphogenesis. The treph1-1 and
mkk2 mutants also exhibited defects in touch-inducible gene expression. A
non-phosphorylatable site-specific isoform of TREPH1 (S625A) failed to restore
touch-induced flowering delay of treph1-1, indicating the necessity of S625 for
TREPH1 function and providing evidence consistent with the possible functional
relevance of the touch-regulated TREPH1 phosphorylation. Bioinformatic analysis
and biochemical subcellular fractionation of TREPH1 protein indicate that it is
a soluble protein. Altogether, these findings identify new protein players in
Arabidopsis thigmomorphogenesis regulation, suggesting that protein
phosphorylation may play a critical role in plant force responses.


  Beta-turn prediction is useful in protein function studies and experimental
design. Although recent approaches using machine-learning techniques such as
SVM, neural networks, and K-NN have achieved good results for beta-turn
pre-diction, there is still significant room for improvement. As previous
predictors utilized features in a sliding window of 4-20 residues to capture
interactions among sequentially neighboring residues, such feature engineering
may result in incomplete or biased features, and neglect interactions among
long-range residues. Deep neural networks provide a new opportunity to address
these issues. Here, we proposed a deep dense inception network (DeepDIN) for
beta-turn prediction, which takes advantages of the state-of-the-art deep
neural network design of the DenseNet and the inception network. A test on a
recent BT6376 benchmark shows that the DeepDIN outperformed the previous best
BetaTPred3 significantly in both the overall prediction accuracy and the
nine-type beta-turn classification. A tool, called MUFold-BetaTurn, was
developed, which is the first beta-turn prediction tool utilizing deep neural
networks. The tool can be downloaded at
http://dslsrv8.cs.missouri.edu/~cf797/MUFoldBetaTurn/download.html.


  Phylodynamics is an area of population genetics that uses genetic sequence
data to estimate past population dynamics. Modern state-of-the-art Bayesian
nonparametric methods for recovering population size trajectories of unknown
form use either change-point models or Gaussian process priors. Change-point
models suffer from computational issues when the number of change-points is
unknown and needs to be estimated. Gaussian process-based methods lack local
adaptivity and cannot accurately recover trajectories that exhibit features
such as abrupt changes in trend or varying levels of smoothness. We propose a
novel, locally-adaptive approach to Bayesian nonparametric phylodynamic
inference that has the flexibility to accommodate a large class of functional
behaviors. Local adaptivity results from modeling the log-transformed effective
population size a priori as a horseshoe Markov random field, a recently
proposed statistical model that blends together the best properties of the
change-point and Gaussian process modeling paradigms. We use simulated data to
assess model performance, and find that our proposed method results in reduced
bias and increased precision when compared to contemporary methods. We also use
our models to reconstruct past changes in genetic diversity of human hepatitis
C virus in Egypt and to estimate population size changes of ancient and modern
steppe bison. These analyses show that our new method captures features of the
population size trajectories that were missed by the state-of-the-art methods.


  A reparametrized Dirichlet-multinomial distribution is introduced, and the
covariance matrix, as well as, the algorithm for calculating the PDF for n
species are provided. The distribution is suited for modelling the joint
distribution of pin-point cover data of spatially aggregated plant species, and
the parametrization ensures that the degree of spatially aggregation is
modelled by a parameter and the mean cover of the species is modelled by the
other parameters. This last property is convenient for using the distribution
in Bayesian hierarchical models where the mean cover of the different species
typically is modelled as latent variables.


  We study periodic, quasi-periodic (Thue-Morse, Fibonacci, Period Doubling,
Rudin-Shapiro), fractal (Cantor, generalized Cantor), Kolakoski and random
binary sequences using a tight-binding wire model, where a site is a monomer
(e.g., in DNA, a base pair). We use B-DNA as our prototype system. All
sequences have purines, guanine (G) or adenine (A) on the same strand, i.e.,
our prototype binary alphabet is (G,A). Our aim is to examine the influence of
sequence intricacy and magnitude of parameters on energy structure,
localization and charge transport. We study quantities such as autocorrelation
function, eigenspectra, density of states, Lyapunov exponents, transmission
coefficients and current-voltage curves. We show that the degree of sequence
intricacy and the presence of correlations decisively affect the aforementioned
physical properties. Periodic segments have enhanced transport properties.
Specifically, in homogeneous sequences transport efficiency is maximum. There
are several deterministic aperiodic sequences that can support significant
currents, depending on the Fermi level of the leads. Random sequences is the
less efficient category.


  We study a biological autoregulation process, involving a protein that
enhances its own transcription, in a parameter region where bistability would
be present in the absence of fluctuations. We calculate the rate of
fluctuation-induced rare transitions between locally-stable states using a path
integral formulation and Master and Chapman-Kolmogorov equations. As in simpler
models for rare transitions, the rate has the form of the exponential of a
quantity $S_0$ (a "barrier") multiplied by a prefactor $\eta$. We calculate
$S_0$ and $\eta$ first in the bursting limit (where the ratio $\gamma$ of the
protein and mRNA lifetimes is very large). In this limit, the calculation can
be done almost entirely analytically, and the results are in good agreement
with simulations. For finite $\gamma$ numerical calculations are generally
required. However, $S_0$ can be calculated analytically to first order in
$1/\gamma$, and the result agrees well with the full numerical calculation for
all $\gamma > 1$. Employing a method used previously on other problems, we find
we can account qualitatively for the way the prefactor $\eta$ varies with
$\gamma$, but its value is 15-20% higher than that inferred from simulations.


  A hierarchy of timescales is ubiquitous in biological systems, where
enzymatic reactions play an important role because they can hasten the
relaxation to equilibrium. We introduced a statistical physics model of
interacting spins that also incorporates enzymatic reactions to extend the
classic model for allosteric regulation. Through Monte Carlo simulations, we
found that the relaxation dynamics are much slower than the elementary
reactions and are logarithmic in time with several plateaus, as is commonly
observed for glasses. This is because of the kinetic constraints from the
cooperativity via the competition for an enzyme, which has different affinity
for molecules with different structures. Our model showed symmetry breaking in
the relaxation trajectories that led to inherently kinetic transitions without
any correspondence to the equilibrium state. In this paper, we discuss the
relevance of these results for diverse responses in biology.


  This article is dedicated to the late Giorgio Israel. R{\'e}sum{\'e}. The aim
of this article is to propose on the one hand a brief history of modeling
starting from the works of Fibonacci, Robert Malthus, Pierre Francis Verhulst
and then Vito Volterra and, on the other hand, to present the main hypotheses
of the very famous but very little known predator-prey model elaborated in the
1920s by Volterra in order to solve a problem posed by his son-in-law, Umberto
D'Ancona. It is thus shown that, contrary to a widely-held notion, Volterra's
model is realistic and his seminal work laid the groundwork for modern
population dynamics and mathematical ecology, including seasonality, migration,
pollution and more. 1. A short history of modeling 1.1. The Malthusian model.
If the rst scientic view of population growth seems to be that of Leonardo
Fibonacci [2], also called Leonardo of Pisa, whose famous sequence of numbers
was presented in his Liber abaci (1202) as a solution to a population growth
problem, the modern foundations of population dynamics clearly date from Thomas
Robert Malthus [20]. Considering an ideal population consisting of a single
homogeneous animal species, that is, neglecting the variations in age, size and
any periodicity for birth or mortality, and which lives alone in an invariable
environment or coexists with other species without any direct or indirect
inuence, he founded in 1798, with his celebrated claim Population, when
unchecked, increases in a geometrical ratio, the paradigm of exponential
growth. This consists in assuming that the increase of the number N (t) of
individuals of this population, during a short interval of time, is
proportional to N (t). This translates to the following dierential equation :
(1) dN (t) dt = $\epsilon$N (t) where $\epsilon$ is a constant factor of
proportionality that represents the growth coe-cient or growth rate. By
integrating (1) we obtain the law of exponential growth or law of Malthusian
growth (see Fig. 1). This law, which does not take into account the limits
imposed by the environment on growth and which is in disagreement with the
actual facts, had a profound inuence on Charles Darwin's work on natural
selection. Indeed, Darwin [1] founded the idea of survival of the ttest on the
1. According to Frontier and Pichod-Viale [3] the correct terminology should be
population kinetics, since the interaction between species cannot be
represented by forces. 2. A population is dened as the set of individuals of
the same species living on the same territory and able to reproduce among
themselves.


  The free rider problem is one of the most well-studied problems in economics.
The solution proposed mainly is punitive in order to deter people from free
riding. In this work, we introduce quantum strategies and study the problem in
the thermodynamic limit by drawing analogies with the 1D Ising model. We
observe that for maximum entanglement, irrespective of the payoffs, quantum
strategy is the equilibrium solution, solving the free rider problem.


  Reaction networks are mainly used to model the time-evolution of molecules of
interacting chemical species. Stochastic models are typically used when the
counts of the molecules are low, whereas deterministic models are used when the
counts are in high abundance. In 2011, the notion of `tiers' was introduced to
study the long time behavior of deterministically modeled reaction networks
that are weakly reversible and have a single linkage class. This `tier' based
argument was analytical in nature. Later, in 2014, the notion of a strongly
endotactic network was introduced in order to generalize the previous results
from weakly reversible networks with a single linkage class to this wider
family of networks. The point of view of this later work was more geometric and
algebraic in nature. The notion of strongly endotactic networks was later used
in 2018 to prove a large deviation principle for a class of stochastically
modeled reaction networks.
  We provide an analytical characterization of strongly endotactic networks in
terms of tier structures. By doing so, we shed light on the connection between
the two points of view, and also make available a new proof technique for the
study of strongly endotactic networks. We show the power of this new technique
in two distinct ways. First, we demonstrate how the main previous results
related to strongly endotactic networks, both for the deterministic and
stochastic modeling choices, can be quickly obtained from our characterization.
Second, we demonstrate how new results can be obtained by proving that a
sub-class of strongly endotactic networks, when modeled stochastically, is
positive recurrent. Finally, and similarly to recent independent work by Agazzi
and Mattingly, we provide an example which closes a conjecture in the negative
by showing that stochastically modeled strongly endotactic networks can be
transient (and even explosive).


  The homochirality of biomolecules remains one of the outstanding puzzles
concerning the beginning of life. Chiral amplification of a randomly perturbed
racemic mixture of chiral molecules is a well-accepted prerequisite for all
routes to biological homochirality. Some models have suggested that such
amplification occurred due to asymmetric discrimination of chiral biotic or
prebiotic molecules when they adsorbed onto crystalline surfaces. While chiral
amplification has been demonstrated on surfaces of both chiral and achiral
crystals, the mechanism that would produce an enantiomeric imbalance in the
chiral surfaces themselves has not been addressed. Here we report strong chiral
amplification in the colloidal synthesis of intrinsically chiral lanthanide
phosphate nanocrystals, quantitatively measured via the circularly polarized
luminescence of the lanthanide ions within the nanocrystals. The amplification
involves spontaneous symmetry breaking into either left- or right-handed
nanocrystals below a critical temperature. Furthermore, chiral tartaric acid
molecules in the solution act as an external chiral field, sensitively
directing the amplified nanocrystal handedness through a discontinuous
transition between left- and right-handed excess. These characteristics suggest
a conceptual framework for chiral amplification, based on the statistical
thermodynamics of critical phenomena, which we use to quantitatively account
for the observations. Our results demonstrate how chiral minerals with high
enantiomeric excess could have grown locally in a primordial racemic aqueous
environment.


  Randomized trials of infectious disease interventions, such as vaccines,
often focus on groups of connected or potentially interacting individuals. When
the pathogen of interest is transmissible between study subjects, interference
may occur: individual infection outcomes may depend on treatments received by
others. Epidemiologists have defined the primary causal effect of interest --
called the "susceptibility effect" -- as a contrast in infection risk under
treatment versus no treatment, while holding exposure to infectiousness
constant. A related quantity -- the "direct effect" -- is defined as an
unconditional contrast between the infection risk under treatment versus no
treatment. The purpose of this paper is to show that under a widely recommended
randomization design, the direct effect may fail to recover the sign of the
true susceptibility effect of the intervention in a randomized trial when
outcomes are contagious. The analytical approach uses structural features of
infectious disease transmission to define the susceptibility effect. A new
probabilistic coupling argument reveals stochastic dominance relations between
potential infection outcomes under different treatment allocations. The results
suggest that estimating the direct effect under randomization may provide
misleading inferences about the effect of an intervention -- such as a vaccine
-- when outcomes are contagious.


  Resolving the relationship between biodiversity and ecosystem functioning has
been one of the central goals of modern ecology. Early debates about the
relationship were finally resolved with the advent of a statistical
partitioning scheme that decomposed the biodiversity effect into a "selection"
effect and a "complementarity" effect. We prove that both the biodiversity
effect and its statistical decomposition into selection and complementarity are
fundamentally flawed because these methods use a na\"ive null expectation based
on neutrality, likely leading to an overestimate of the net biodiversity
effect, and they fail to account for the nonlinear abundance-ecosystem
functioning relationships observed in nature. Furthermore, under such
nonlinearity no statistical scheme can be devised to partition the biodiversity
effects. We also present an alternative metric providing a more reasonable
estimate of biodiversity effect. Our results suggest that all studies conducted
since the early 1990s likely overestimated the positive effects of biodiversity
on ecosystem functioning.


  We study the dynamics of knotted deformable closed chains sedimenting in a
viscous fluid. We show experimentally that trefoil and other torus knots often
attain a remarkably regular horizontal toroidal structure while sedimenting,
with a number of intertwined loops, oscillating periodically around each other.
We then recover this motion numerically and find out that it is accompanied by
a very slow rotation around the vertical symmetry axis. We analyze the
dependence of the characteristic time scales on the chain flexibility and
aspect ratio. It is observed in the experiments that this oscillating mode of
the dynamics can spontaneously form even when starting from a qualitatively
different initial configuration. In numerical simulations, the oscillating
modes are usually present as transients or final stages of the evolution,
depending on chain aspect ratio and flexibility, and the number of loops.


  When three species compete cyclically in a well-mixed, stochastic system of
$N$ individuals, extinction is known to typically occur at times scaling as the
system size $N$. This happens, for example, in rock-paper-scissors games or
conserved Lotka-Volterra models in which every pair of individuals can interact
on a complete graph. Here we show that if the competing individuals also have a
"social temperament" to be either introverted or extroverted, leading them to
cut or add links respectively, then long-living state in which all species
coexist can occur when both introverts and extroverts are present. These states
are non-equilibrium quasi-steady states, maintained by a subtle balance between
species competition and network dynamcis. Remarkably, much of the phenomena is
embodied in a mean-field description. However, an intuitive understanding of
why diversity stabilizes the co-evolving node and link dynamics remains an open
issue.


  As infectious disease outbreaks emerge, public health agencies often enact
vaccination and social distancing measures to slow transmission. Their success
depends on not only strategies and resources, but also public adherence.
Individual willingness to take precautions may be influenced by global factors,
such as news media, or local factors, such as infected family members or
friends. Here, we compare three modes of epidemiological decision-making in the
midst of a growing outbreak. Individuals decide whether to adopt a recommended
intervention based on overall disease prevalence, the proportion of social
contacts infected, or the number of social contacts infected. While all
strategies can substantially mitigate transmission, vaccinating (or self
isolating) based on the number of infected acquaintances is expected to achieve
the greatest herd immunity and number of infections averted, while requiring
the fewest intervention resources.


  Chemical clock reactions are characterised by a relatively long induction
period followed by a rapid `switchover' during which the concentration of a
\emph{clock chemical} rises rapidly. In addition to their interest in chemistry
education, these reactions are relevant to industrial and biochemical
applications. A substrate-depletive, non-autocatalytic clock reaction involving
household chemicals (vitamin C, iodine, hydrogen peroxide and starch) is
modelled mathematically via a system of nonlinear ordinary differential
equations. Following dimensional analysis the model is analysed in the phase
plane and via matched asymptotic expansions. Asymptotic approximations are
found to agree closely with numerical solutions in the appropriate time
regions. Asymptotic analysis also yields an approximate formula for the
dependence of switchover time on initial concentrations and the rate of the
slow reaction. This formula is tested via `kitchen sink chemistry' experiments,
and is found to enable a good fit to experimental series varying in initial
concentrations of both iodine and vitamin C. The vitamin C clock reaction
provides an accessible model system for mathematical chemistry.


  Even in the steady-state, the number of biomolecules in living cells
fluctuates dynamically; and the frequency spectrum of this chemical fluctuation
carries valuable information about the mechanism and the dynamics of the
intracellular reactions creating these biomolecules. Although recent advances
in single-cell experimental techniques enable the direct monitoring of the
time-traces of the biological noise in each cell, the development of the
theoretical tools needed to extract the information encoded in the stochastic
dynamics of intracellular chemical fluctuation is still in its adolescence.
Here, we present a simple and general equation that relates the power-spectrum
of the product number fluctuation to the product lifetime and the reaction
dynamics of the product creation process. By analyzing the time traces of the
protein copy number using this theory, we can extract the power spectrum of the
mRNA number, which cannot be directly measured by currently available
experimental techniques. From the power spectrum of the mRNA number, we can
further extract quantitative information about the transcriptional regulation
dynamics. Our power spectrum analysis of gene expression noise is demonstrated
for the gene network model of luciferase expression under the control of the
Bmal 1a promoter in mouse fibroblast cells. Additionally, we investigate how
the non-Poisson reaction dynamics and the cell-to-cell heterogeneity in
transcription and translation affect the power-spectra of the mRNA and protein
number.


  Biological systems are typically highly open, non-equilibrium systems that
are very challenging to understand from a statistical mechanics perspective.
While statistical treatments of evolutionary biological systems have a long and
rich history, examination of the time-dependent non-equilibrium dynamics has
been less studied. In this paper we first derive a generalized master equation
in the genotype space for diploid organisms incorporating the processes of
selection, mutation, recombination, and reproduction. The master equation is
defined in terms of continuous time and can handle an arbitrary number of gene
loci and alleles, and can be defined in terms of an absolute population or
probabilities. We examine and analytically solve several prototypical cases
which illustrate the interplay of the various processes and discuss the
timescales of their evolution. The entropy production during the evolution
towards steady state is calculated and we find that it agrees with predictions
from non-equilibrium statistical mechanics where it is large when the
population distribution evolves towards a more viable genotype. The stability
of the non-equilibrium steady state is confirmed using the Glansdorff-Prigogine
criterion.


  The discovery of Paranthropus deyiremeda in 3.3-3.5 million year old fossil
sites in Afar (Haile-Selassie, 2015), together with 30% of the gorilla genome
showing lineage sorting between humans and chimpanzees (Scally, 2012), and a
NUMT ("nuclear mitochondrial DNA segment") that is shared by both gorillas,
humans and chimpanzees, and that dates back to 6 million years ago (Popadin,
2017), is conclusive evidence that introgression from the gorilla lineage
caused the speciation of both the Australopithecus lineage and the Paranthropus
lineage, providing a lens into the gorilla-like features within Paranthropus,
as well as traits within Homo that originate from the gorilla branch, such as a
high opposable thumb index (Alm\'ecija, 2015), an adducted great toe (Tocheri,
2011, McHenry, 2006), and large deposits of subcutaneous fat.


  The importance of microscopic details on cooperation level is an intensively
studied aspect of evolutionary game theory. Interestingly, these details become
crucial on heterogeneous populations where individuals may possess diverse
traits. By introducing a coevolutionary model in which not only strategies but
also individual dynamical features may evolve we revealed that the formerly
established conclusion is not necessarily true when different updating rules
are on stage. In particular, we apply two strategy updating rules, imitation
and Death-Birth rule, which allow local selection in a spatial system. Our
observation highlights that the microscopic feature of dynamics, like the level
of learning activity, could be a fundamental factor even if all players share
the same trait uniformly.


  The transition from early long-term potentiation (E-LTP) to late LTP (L-LTP)
involves protein synthesis and degradation. L-LTP is blocked by inhibiting
either protein synthesis or proteasome-dependent degradation prior to and
during a tetanic stimulus, but paradoxically, L-LTP is not blocked when
synthesis and degradation are inhibited simultaneously, suggesting
counter-acting positive and negative proteins regulate L-LTP. To investigate
this paradox, we modeled LTP at the Schaffer collateral synapse. Nine
differential equations describe the levels of positive and negative regulator
proteins (PP and NP) and transitions among five discrete synaptic states, a
basal state (BAS), three E-LTP states (EP1, EP2, ED), and a L-LTP state (LP). A
stimulus initiates the transition from BAS to EP1 and from EP1 to EP2,
initiates the synthesis of PP and NP, and activates the ubiquitin-proteasome
system (UPS). UPS mediates transitions of EP1 and EP2 to ED and the degradation
of NP. The conversion of E-LTP to L-LTP is mediated by a PP-dependent
transition from ED to LP. NP mediates reversal of EP2 to BAS. This model
simulates empirical observations: 1) normal L-LTP, 2) block by either
proteasome inhibitor or protein synthesis inhibitor alone, and 3) preservation
of L-LTP when both inhibitors are applied together. Elements of this abstract
model can be correlated with specific molecules and processes. Moreover, the
model makes testable predictions, such as a unique synaptic state ED that
precedes the transition to L-LTP, and a time window for the action of the UPS
(during the transitions from EP1 and EP2 to ED). Tests of these predictions
will provide insights into the processes of long-term synaptic plasticity.


  We use dynamical generating functionals to study the stability and size of
communities evolving in Lotka-Volterra systems with random interaction
coefficients. The size of the eco-system is not set from the beginning.
Instead, we start from a set of possible species, which may undergo extinction.
How many species survive depends on the properties of the interaction matrix;
the size of the resulting food web at stationarity is a property of the system
itself in our model, and not a control parameter as in most studies based on
random matrix theory. We find that prey-predator relations enhance stability,
and that variability of species interactions promotes instability. Complexity
of inter-species couplings leads to reduced sizes of ecological communities.
Dynamically evolved community size and stability are hence positively
correlated.


  The central dogma of molecular biology is the principal framework for
understanding how nucleic acid information is propagated and used by living
systems to create complex biomolecules. Here, by integrating the structural and
dynamic paradigms of DNA nanotechnology, we present a rationally designed
synthetic platform which functions in an analogous manner to create complex DNA
nanostructures. Starting from one type of DNA nanostructure, DNA strand
displacement circuits were designed to interact and pass along the information
encoded in the initial structure to mediate the self-assembly of a different
type of structure, the final output structure depending on the type of circuit
triggered. Using this concept of a DNA structure "trans-assembling" a different
DNA structure through non-local strand displacement circuitry, four different
schemes were implemented. Specifically, 1D ladder and 2D double-crossover (DX)
lattices were designed to kinetically trigger DNA circuits to activate
polymerization of either ring structures or another type of DX lattice under
enzyme-free, isothermal conditions. In each scheme, the desired multilayer
reaction pathway was activated, among multiple possible pathways, ultimately
leading to the downstream self-assembly of the correct output structure.


  This paper introduces a unified framework called cooperative extensive form
games, which (i) generalizes standard non-cooperative games, and (ii) allows
for more complex coalition formation dynamics than previous concepts like
coalition-proof Nash equilibrium. Central to this framework is a novel solution
concept called cooperative equilibrium system (CES). CES differs from Nash
equilibrium in two important respects. First, a CES is immune to both
unilateral and multilateral `credible' deviations. Second, unlike Nash
equilibrium, whose stability relies on the assumption that the strategies of
non-deviating players are held fixed, CES allows for the possibility that
players may regroup and adjust their strategies in response to a deviation. The
main result establishes that every cooperative extensive form game, possibly
with imperfect information, possesses a CES. For games with perfect
information, the proof is constructive. This framework is broadly applicable in
contexts such as oligopolistic markets and dynamic political bargaining.


  A large class of mesoscopic or macroscopic flocking theories are
coarse-grained from microscopic models that feature binary interactions as the
chief aligning mechanism. However while such theories seemingly predict the
existence of polar order with just binary interactions, actomyosin motility
assay experiments show that binary interactions are insufficient to obtain
polar order, especially at high densities. To resolve this paradox, here we
introduce a solvable one-dimensional flocking model and derive its stochastic
hydrodynamics. We show that two-body interactions are insufficient to generate
polar order unless the noise is non-Gaussian. We show that noisy three-body
interactions in the microscopic theory allow us to capture all essential
dynamical features of the flocking transition, in systems that achieve
orientational order above a critical density.


  One of the main aims of phylogenetics is the reconstruction of the correct
evolutionary tree when data concerning the underlying species set are given.
These data typically come in the form of DNA, RNA or protein alignments, which
consist of various characters (also often referred to as sites). Often,
however, tree reconstruction methods based on criteria like maximum parsimony
may fail to provide a unique tree for a given dataset, or, even worse,
reconstruct the `wrong' tree (i.e. a tree that differs from the one that
generated the data). On the other hand it has long been known that if the
alignment consists of all the characters that correspond to edges of a
particular tree, i.e. they all require exactly $k=1$ substitution to be
realized on that tree, then this tree will be recovered by maximum parsimony
methods. This is based on Buneman's theorem in mathematical phylogenetics. It
is the goal of the present manuscript to extend this classic result as follows:
We prove that if an alignment consists of all characters that require exactly
$k=2$ substitutions on a particular tree, this tree will always be the unique
maximum parsimony tree (and we also show that this can be generalized to
characters which require at most $k=2$ substitutions). In particular, this also
proves a conjecture based on a recently published observation by Goloboff et
al. affirmatively for the special case of $k=2$.


  Receptor-ligand interactions are essential for biological function and their
binding strength is commonly explained in terms of static lock-and-key models
based on molecular complementarity. However, detailed information of the full
unbinding pathway is often lacking due, in part, to the static nature of atomic
structures and ensemble averaging inherent to bulk biophysics approaches. Here
we combine molecular dynamics and high-speed force spectroscopy on the
streptavidin-biotin complex to determine the binding strength and unbinding
pathways over the widest dynamic range. Experiment and simulation show
excellent agreement at overlapping velocities and provided evidence of the
unbinding mechanisms. During unbinding, biotin crosses multiple energy barriers
and visits various intermediate states far from the binding pocket while
streptavidin undergoes transient induced fits, all varying with loading rate.
This multistate process slows down the transition to the unbound state and
favors rebinding, thus explaining the long lifetime of the complex. We provide
an atomistic, dynamic picture of the unbinding process, replacing a simple
two-state picture with one that involves many routes to the lock and
rate-dependent induced-fit motions for intermediates, which might be relevant
for other receptor-ligand bonds.


  Competition between biological species in marine environments is affected by
the motion of the surrounding fluid. An effective 2D compressibility can arise,
for example, from the convergence and divergence of water masses at the depth
at which passively traveling photosynthetic organisms are restricted to live.
In this report, we seek to quantitatively study genetics under flow. To this
end, we couple an off-lattice agent-based simulation of two populations in 1D
to a weakly compressible velocity field--first a sine wave and then a shell
model of turbulence. We find for both cases that even in a regime where the
overall population structure is approximately unaltered, the flow can
significantly diminish the effect of a selective advantage on fixation
probabilities. We understand this effect in terms of the enhanced survival of
organisms born at sources in the flow and the influence of Fisher genetic
waves.


  Haemostasis is governed by a highly complex system of interacting proteins.
Due to the central role of thrombin, thrombin generation and specifically the
thrombin generation curve (TGC) is commonly used as an indicator of haemostatic
activity. Functional characteristics of the haemostatic system in neonates and
children are significantly different compared with adults; at the same time
plasma levels of haemostatic proteins vary considerably with age. However,
relating one to the other has been difficult, both due to significant
inter-individual differences for individuals of similar age and the complexity
of the biochemical reactions underlying haemostasis. Mathematical modelling has
been very successful at representing the biochemistry of blood clotting. In
this study we address the challenge of large inter-individual variability by
parameterising the Hockin-Mann model with data from individual patients, across
different age groups from neonates to adults. Calculating TGCs for each patient
of a specific age group provides us with insight into the variability of
haemostatic activity across that age group. From our model we observe that two
commonly used metrics for haemostatic activity are significantly lower in
neonates than in older patients. Because both metrics are strongly determined
by prothrombin and prothrombin levels are considerably lower in neonates we
conclude that decreased haemostatic activity in neonates is due to lower
prothrombin availability.


  Statistical analysis of alignments of large numbers of protein sequences has
revealed "sectors" of collectively coevolving amino acids in several protein
families. Here, we show that selection acting on any functional property of a
protein, represented by an additive trait, can give rise to such a sector. As
an illustration of a selected trait, we consider the elastic energy of an
important conformational change within an elastic network model, and we show
that selection acting on this energy leads to correlations among residues. For
this concrete example and more generally, we demonstrate that the main
signature of functional sectors lies in the small-eigenvalue modes of the
covariance matrix of the selected sequences. However, secondary signatures of
these functional sectors also exist in the extensively-studied large-eigenvalue
modes. Our simple, general model leads us to propose a principled method to
identify functional sectors, along with the magnitudes of mutational effects,
from sequence data. We further demonstrate the robustness of these functional
sectors to various forms of selection, and the robustness of our approach to
the identification of multiple selected traits.


  The determination of rates of body growth is the first step in many aquatic
population studies and fisheries stock assessments. ELEFAN (Electronic LEngth
Frequency ANalysis) is a widely used method to fit a growth curve to
length-frequency distribution (LFD) data. However, up to now, it was not
possible to assess its accuracy or the uncertainty inherent of this method, or
to obtain confidence intervals for growth parameters within an unconstrained
search space. In this study, experiments were conducted to assess the precision
and accuracy of bootstrapped and single-fit ELEFAN-based curve fitting methods,
using synthetic LFDs with known input parameters and a real data set of Abra
alba shell lengths. The comparison of several types of bootstrap experiments
and their outputs (95% confidence intervals and confidence contour plots)
provided a first glimpse into the accuracy of modern ELEFAN-based fit methods.
The main components of uncertainty (precision and reproducibility of fit
algorithms, seed effects, sample size and matrix information content) could be
assessed from partial bootstraps. Uncertainty was mainly determined by LFD
matrix size, total number of non-zero bins and the sampling of large-sized
individuals. A new pseudo-Rsquared index for the goodness-of-fit of VBGF models
to LFD data is proposed. For a large, perfect synthetic data set,
pseudo-RsquaredPhi was very high (88 to 100%), indicating an excellent fit of
the VBGF model. The small Abra alba data set showed a low pseudo-RsquaredPhi,
from to 54% to 68%, indicating the need for more samples and a larger LFD data
matrix. New, robust, bootstrap-based methods for curve fitting are presented
and discussed. This study demonstrates a promising new path for length-based
analyses of growth and mortality in natural populations, which are the basis
for a new suite of methods that are included in the new fishboot package.


  Background Synonymous codon choice is mainly affected by mutation and
selection. For the majority of genes within a genome, mutational pressure is
the major driving force, but selective strength can be strong and dominant for
specific set of genes or codons. More specifically, the selective strength on
translational efficiency and accuracy increases with the gene's expression
level. Many statistical approaches have been developed to evaluate and quantify
the selective profile in codon usage, including S index and Akashi's test, but
no program or pipeline has been developed that includes these tests and
automates the calculation.
  Results In this study, we release an R package SSCU (selective strength for
codon usage, v2.4.0), which includes tools for codon usage analyses. The
package identifies optimal codons using two approaches (comparative and
correlative methods), implements well-established statistics for detecting
codon selection, such as S index, Akashi's test, and estimates standard genomic
statistics, such as genomic GC3, RSCU and Nc.
  Conclusions The package is useful for researchers working on the codon usage
analysis, and thus has general interest to the biological research community.
The package is deposited and curated at the Bioconductor site, and has
currently been downloaded for more than 2000 times and ranked as top 50%
packages.


  We present a mechanistic formalism for the study of evolutionary dynamics
models based on the diffusion approximation described by the Kimura Equation.
In this formalism, the central component is the fitness potential, from which
we obtain an expression for the amount of work necessary for a given type to
reach fixation. In particular, within this interpretation, we develop a
graphical analysis --- similar to the one used in classical mechanics ---
providing the basic tool for a simple heuristic that describes both the short
and long term dynamics. As a by-product, we provide a new definition of an
evolutionary stable state in finite populations that includes the case of mixed
populations. We finish by showing that our theory -- rigorous for two types
evolution without mutations-- is also consistent with the multi-type case, and
with the inclusion of rare mutations.


  We report the analysis of activation by corticosteroids and progesterone of
full-length mineralocorticoid receptor (MR) from elephant shark, a
cartilaginous fish belonging to the oldest group of jawed vertebrates. Based on
their measured activities, aldosterone, cortisol, 11-deoxycorticosterone,
corticosterone, 11-deoxcortisol, progesterone and 19-norprogesterone are
potential physiological mineralocorticoids. However, aldosterone, the
physiological mineralocorticoid in humans and other terrestrial vertebrates, is
not found in cartilaginous or ray-finned fishes. Because progesterone is a
precursor for corticosteroids that activate elephant shark MR, we propose that
progesterone was an ancestral ligand for elephant shark MR. Although
progesterone activates ray-finned fish MRs, progesterone does not activate
human, amphibian or alligator MRs, suggesting that during the transition to
terrestrial vertebrates, progesterone lost the ability to activate the MR.
Comparison of RNA-sequence analysis of elephant shark MR with that of human MR
suggests that MR expression in the human brain, heart, ovary, testis and other
non-epithelial tissues evolved in cartilaginous fishes. Together, these data
suggest that progesterone-activated MR may have unappreciated functions in
elephant shark ovary and testis.


  This work addresses whether a reaction network, taken with mass-action
kinetics, is multistationary, that is, admits more than one positive steady
state in some stoichiometric compatibility class. We build on previous work on
the effect that removing or adding intermediates has on multistationarity, and
also on methods to detect multistationarity for networks with a binomial steady
state ideal. In particular, we provide a new determinant criterion to decide
whether a network is multistationary, which applies when the network obtained
by removing intermediates has a binomial steady state ideal. We apply this
method to easily characterize which subsets of complexes are responsible for
multistationarity; this is what we call the \emph{multistationarity structure}
of the network. We use our approach to compute the multistationarity structure
of the $n$-site sequential distributive phosphorylation cycle for arbitrary n.


  A very stable binding site for the interaction between an pentameric
oligothiophene and an amyloid-$\beta$(1-42) fibril has been identified by means
of non-biased molecular dynamics simulations. In this site, the probe is locked
in an all-trans conformation with a Coulombic binding energy of 1,200 kJ/mol
due to the interactions between the anionic carboxyl groups of the probe and
the cationic $\epsilon$-amino groups in the lysine side chain. Upon binding,
the conformationally restricted probes show a pronounced increase in molecular
planarity. This is in-line with the observed changes in luminescence properties
that serve as the foundation for their use as biomarkers.


  Mathematical modelers have long known of a "rule of thumb" referred to as the
Linear Chain Trick (LCT; aka the Gamma Chain Trick): a technique used to
construct mean field ODE models from continuous-time stochastic state
transition models where the time an individual spends in a given state (i.e.,
the dwell time) is Erlang distributed (i.e., gamma distributed with integer
shape parameter). Despite the LCT's widespread use, we lack general theory to
facilitate the easy application of this technique, especially for complex
models. This has forced modelers to choose between constructing ODE models
using heuristics with oversimplified dwell time assumptions, using time
consuming derivations from first principles, or to instead use non-ODE models
(like integro-differential equations or delay differential equations) which can
be cumbersome to derive and analyze. Here, we provide analytical results that
enable modelers to more efficiently construct ODE models using the LCT or
related extensions. Specifically, we 1) provide novel extensions of the LCT to
various scenarios found in applications; 2) provide formulations of the LCT and
it's extensions that bypass the need to derive ODEs from integral or stochastic
model equations; and 3) introduce a novel Generalized Linear Chain Trick (GLCT)
framework that extends the LCT to a much broader family of distributions,
including the flexible phase-type distributions which can approximate
distributions on $\mathbb{R}^+$ and be fit to data. These results give modelers
more flexibility to incorporate appropriate dwell time assumptions into mean
field ODEs, including conditional dwell time distributions, and these results
help clarify connections between individual-level stochastic model assumptions
and the structure of corresponding mean field ODEs.


  We analyze an epidemiological model to evaluate the effectiveness of multiple
means of control in malaria-endemic areas. The mathematical model consists of a
system of several ordinary differential equations, and is based on a
multicompartment representation of the system. The model takes into account the
mutliple resting-questing stages undergone by adult female mosquitos during the
period in which they function as disease vectors. We compute the basic
reproduction number $\mathcal R_0$, and show that that if $\mathcal R_0<1$, the
disease free equilibrium (DFE) is globally asymptotically stable (GAS) on the
non-negative orthant. If $\mathcal R_0>1$, the system admits a unique endemic
equilibrium (EE) that is GAS. We perform a sensitivity analysis of the
dependence of $\mathcal R_0$ and the EE on parameters related to control
measures, such as killing effectiveness and bite prevention. Finally, we
discuss the implications for a comprehensive, cost-effective strategy for
malaria control.


  Multiple sequence alignment is a basic procedure in molecular biology, and it
is often treated as being essentially a solved computational problem. However,
this is not so, and here I review the evidence for this claim, and outline the
requirements for a solution. The goal of alignment is often stated to be to
juxtapose nucleotides (or their derivatives, such as amino acids) that have
been inherited from a common ancestral nucleotide (although other goals are
also possible). Unfortunately, this is not an operational definition, because
homology (in this sense) refers to unique and unobservable historical events,
and so there can be no objective mathematical function to optimize.
Consequently, almost all algorithms developed for multiple sequence alignment
are based on optimizing some sort of compositional similarity (similarity =
homology + analogy). As a result, many, if not most, practitioners either
manually modify computer-produced alignments or they perform de novo manual
alignment, especially in the field of phylogenetics. So, if homology is the
goal, then multiple sequence alignment is not yet a solved computational
problem. Several criteria have been developed by biologists to help them
identify potential homologies (compositional, ontogenetic, topographical and
functional similarity, plus conjunction and congruence), and these criteria can
be applied to molecular data, in principle. Current computer programs do
implement one (or occasionally two) of these criteria, but no program
implements them all. What is needed is a program that evaluates all of the
evidence for the sequence homologies, optimizes their combination, and thus
produces the best hypotheses of homology. This is basically an inference
problem not an optimization problem.


  We consider the stochastic gene expression process with inter-state
flip-flops. An exact steady-state solution to the master equation is
calculated. One of the main goals in this paper is to investigate whether the
probability distribution of gene copies contains even-odd number oscillations.
A master equation previously derived in the adiabatic limit of fast switching
by Kepler and Elston \cite{kepler} suggests that the oscillations should be
present. However our analysis demonstrates that the oscillations do not happen
not only in the adiabatic case but they are entirely absent. We discuss the
adiabatic approximation in detail. The other goal is to establish the master
equation that takes into account external fluctuations that is similar to the
master equation in the adiabatic approximation. The equation allows even-odd
oscillations. The reason the behaviour occurs is an underlying interference of
Poisson and Gaussian processes. The master equation contains an extra term that
describes the gene copy number unconventional diffusion and is responsible for
the oscillations. We also point out to a similar phenomenon in quantum physics.


  Is evolution always gradual or can it make leaps? We examine a mathematical
model of an evolutionary process on a fitness landscape and obtain analytic
solutions for the probability of multi-mutation leaps, that is, several
mutations occurring simultaneously, within a single generation in one genome,
and being fixed all together in the evolving population. The results indicate
that, for typical, empirically observed combinations of the parameters of the
evolutionary process, namely, effective population size, mutation rate, and
distribution of selection coefficients of mutations, the probability of a
multi-mutation leap is low, and accordingly, the contribution of such leaps is
minor at best. However, we show that, taking sign epistasis into account, leaps
could become an important factor of evolution in cases of substantially
elevated mutation rates, such as stress-induced mutagenesis in microbes. We
hypothesize that stress-induced mutagenesis is an evolvable adaptive strategy.


  The dynamics of species' densities depend both on internal and external
variables. Internal variables include frequencies of individuals exhibiting
different phenotypes or living in different spatial locations. External
variables include abiotic factors or non-focal species. These internal or
external variables may fluctuate due to stochastic fluctuations in
environmental conditions. We prove theorems for stochastic persistence and
exclusion for stochastic ecological difference equations accounting for
internal and external variables. Specifically, we use a stochastic analog of
average Lyapunov functions to develop sufficient and necessary conditions for
(i) all population densities spending little time at low densities, and (ii)
population trajectories asymptotically approaching the extinction set with
positive probability. For (i) and (ii), respectively, we provide quantitative
estimates on the fraction of time that the system is near the extinction set,
and the probability of asymptotic extinction as a function of the initial state
of the system. Furthermore, we provide lower bounds for the expected time to
escape neighborhoods of the extinction set. To illustrate the applicability of
our results, we analyze stochastic models of evolutionary games, Lotka-Volterra
dynamics, trait evolution, and spatially structured disease dynamics. Our
analysis of these models demonstrates environmental stochasticity facilitates
coexistence of strategies in the hawk-dove game, but inhibits coexistence in
the rock-paper-scissors game and a Lotka-Volterra predator-prey model.
Furthermore, environmental fluctuations with positive auto-correlations can
promote persistence of evolving populations and persistence of diseases in
patchy landscapes. While our results help close the gap between the persistence
theories for deterministic and stochastic systems, we highlight challenges for
future research.


  Photosynthetic eukaryotes show a remarkable variability in photosynthesis,
including large differences in light harvesting proteins and pigment
composition. In vivo circular spectropolarimetry enables us to probe the
molecular architecture of photosynthesis in a non-invasive and non-destructive
way and, as such, can offer a wealth of physiological and structural
information. In the present study we have measured the circular polarizance of
several multicellular green, red and brown algae and higher plants, which show
large variations in circular spectropolarimetric signals with differences in
both spectral shape and magnitude. Many of the algae display spectral
characteristics not previously reported, indicating a larger variation in
molecular organization than previously assumed. As the strengths of these
signals vary by three orders of magnitude, these results also have important
implications in terms of detectability for the use of circular polarization as
a signature of life.


  A Markovian SIR (Susceptible-Infectious-Recovered) model is considered for
the spread of an epidemic on a configuration model network, in which
susceptible individuals may take preventive measures by dropping edges to
infectious neighbours. An effective degree formulation of the model is used in
conjunction with the theory of density dependent population processes to obtain
a law of large numbers and a functional central limit theorem for the epidemic
as the population size $N \to \infty$, assuming that the degrees of individuals
are bounded. A central limit theorem is conjectured for the final size of the
epidemic. The results are obtained for both the Molloy-Reed (in which the
degrees of individuals are deterministic) and Newman-Strogatz-Watts (in which
the degrees of individuals are independent and identically distributed)
versions of the configuration model. The two versions yield the same limiting
deterministic model but the asymptotic variances in the central limit theorema
are greater in the Newman-Strogatz-Watts version. The basic reproduction number
$R_0$ and the process of susceptible individuals in the limiting deterministic
model, for the model with dropping of edges, are the same as for a
corresponding SIR model without dropping of edges but an increased recovery
rate, though, when $R_0>1$, the probability of a major outbreak is greater in
the model with dropping of edges. The results are specialised to the model
without dropping of edges to yield conjectured central limit theorems for the
final size of Markovian SIR epidemics on configuration-model networks, and for
the giant components of those networks. The theory is illustrated by numerical
studies, which demonstrate that the asymptotic approximations are good, even
for moderate $N$.


  Dengue control in India is a challenging task due to complex healthcare
settings. In yesteryears, an amplification of dengue infections in India posed
the need for introspection of existing dengue control policies. Prior
understanding of the impacts of control interventions is necessary for their
future implementation. In this paper, we propose and analyze a compartmental
model of dengue to assess the impact of active case finding (ACF) on dengue
disease transmission. Currently, primary prevention of dengue is possible only
with vector control and personal protection from the bites of infected
mosquitoes. Although a few experimental studies are performed to assess ACF in
dengue disease, but this is the first attempt to represent and study the
dynamics of disease using ACF as a control strategy. Local and global dynamics
of the system are studied. We use sensitivity analysis to see the effects of
controllable parameters of the model on the basic reproduction number and total
number of infective population. We find that decrease in the biting rate of
mosquitoes, and increase in the rate of hospitalization and/or notification,
death rate of mosquitoes and ACF for asymptomatic and symptomatic individuals
play crucial role for the reduction of disease prevalence. We calibrate our
model to the yearly dengue cases in eight dengue endemic states of India. The
results of our study show that ACF of symptomatic individuals will have
significant effect on dengue case reduction but ACF of asymptomatic individuals
cannot be ignored. Our findings indicate that the healthcare organizations must
focus on ACF of symptomatic as well as asymptomatic individuals along with
personal protection and mosquitoes control to achieve rapid reduction of dengue
cases in India.


  One of the focus areas of modern scientific research is to reveal mysteries
related to genes and their interactions. The dynamic interactions between genes
can be encoded into a gene regulatory network (GRN), which can be used to gain
understanding on the genetic mechanisms behind observable phenotypes. GRN
inference from time series data has recently been a focus area of systems
biology. Due to low sampling frequency of the data, this is a notoriously
difficult problem. We tackle the challenge by introducing the so-called
continuous-time Gaussian process dynamical model, based on Gaussian process
framework that has gained popularity in nonlinear regression problems arising
in machine learning. The model dynamics are governed by a stochastic
differential equation, where the dynamics function is modelled as a Gaussian
process. We prove the existence and uniqueness of solutions of the stochastic
differential equation. We derive the probability distribution for the Euler
discretised trajectories and establish the convergence of the discretisation.
We develop a GRN inference method called BINGO, based on the developed
framework. BINGO is based on MCMC sampling of trajectories of the GPDM and
estimating the hyperparameters of the covariance function of the Gaussian
process. Using benchmark data examples, we show that BINGO is superior in
dealing with poor time resolution and it is computationally feasible.


  Symmetry is closely intertwined with the function, genetics, and chemical
properties of multiprotein complexes. Here, we explore the relation between
structural symmetry and the ability of membrane proteins to sense and induce
membrane curvature, which is a key factor for modulating the shape and
organization of cell membranes. Using coarse-grained simulations and elasticity
theory, we show that the potential for direction-dependent membrane curvature
sensing is limited to asymmetric proteins, dimers, and tetramers, and argue
that one should expect this anisotropy to be strongest for dimers. Odd and
higher-order symmetries strongly suppress directional curvature sensing. This
classification gives a new perspective on the structure-function relation for
membrane proteins, and simplifies the task of translating between molecular
sensing mechanisms and their large-scale cellular consequences.


  Peroxisome proliferator-activated receptors gamma (PPAR{\gamma}) are
ligand-activated controllers of various metabolic actions and insulin
sensitivity. PPAR{\gamma} is thus considered as an important target to treat
type 2 diabetes. Available PPAR{\gamma} drugs (full agonists) have robust
insulin-sensitizing properties but are accompanied by severe side effects
leading to complicated health problems. Here, we have used molecular docking
and a molecular dynamics simulation study to find a novel PPAR{\gamma} ligand
from a natural product. Our study suggests that the inhibition of ceramicine B
in the PPAR{\gamma} ligand-binding domain (LBD) could act as a partial agonist
and block cdk5-mediated phosphorylation. This result may provide an opportunity
for the development of new anti-diabetic drugs by targeting PPAR{\gamma} while
avoiding the side effects associated with full agonists.


  We are increasingly using nature for tourism and recreation, an economic
sector now generating more than 10% of the global GDP and 10% of global total
employment. This growth though has come at a cost and we now have 5930 species
for which tourism and recreation are conservation threats. For the first time
we use global social media data to estimate where people go to experience
nature and determine how this tourism and recreation pressure overlap with the
distribution of threatened species. The more people seek interactions with
nature in an area, the larger the number of species threatened by those
interactions is. Clear crisis areas emerge where many species sensitive to
tourism are exposed to high tourism pressures and those are mainly coastal
marine regions. Our current tourism management approaches are not achieving
biodiversity conservation. The global increase in nature tourism and recreation
is set to continue and we need a global consistent response to mitigate its
biodiversity impact. Like with other extractive industries, we must prioritise
our efforts to diverge tourism away from crisis areas.