  We present a new method that enables the identification and analysis of both
transition and metastable conformational states from atomistic or
coarse-grained molecular dynamics (MD) trajectories. Our algorithm is presented
and studied by using both analytical and actual examples from MD simulations of
the helix-forming peptide Ala5, and of a larger system, the epidermal growth
factor receptor (EGFR) protein. In all cases, our method identifies
automatically the corresponding transition states and metastable conformations
in an optimal way, with the input of a set of relevant coordinates, by
capturing accurately the intrinsic slowest relaxation rate. Our approach
provides a general and easy to implement analysis method that provides unique
insight into the molecular mechanism and the rare but crucial rate limiting
conformational pathways occurring in complex dynamical systems such as
molecular trajectories.


  Motivation: RNA thermometers (RNATs) are cis-regulatory ele- ments that
change secondary structure upon temperature shift. Often involved in the
regulation of heat shock, cold shock and virulence genes, RNATs constitute an
interesting potential resource in synthetic biology, where engineered RNATs
could prove to be useful tools in biosensors and conditional gene regulation.
Results: Solving the 2-temperature inverse folding problem is critical for RNAT
engineering. Here we introduce RNAiFold2T, the first Constraint Programming
(CP) and Large Neighborhood Search (LNS) algorithms to solve this problem.
Benchmarking tests of RNAiFold2T against existent programs (adaptive walk and
genetic algorithm) inverse folding show that our software generates two orders
of magnitude more solutions, thus allow- ing ample exploration of the space of
solutions. Subsequently, solutions can be prioritized by computing various
measures, including probability of target structure in the ensemble, melting
temperature, etc. Using this strategy, we rationally designed two thermosensor
internal ribosome entry site (thermo-IRES) elements, whose normalized
cap-independent transla- tion efficiency is approximately 50% greater at 42?C
than 30?C, when tested in reticulocyte lysates. Translation efficiency is lower
than that of the wild-type IRES element, which on the other hand is fully
resistant to temperature shift-up. This appears to be the first purely
computational design of functional RNA thermoswitches, and certainly the first
purely computational design of functional thermo-IRES elements. Availability:
RNAiFold2T is publicly available as as part of the new re- lease RNAiFold3.0 at
https://github.com/clotelab/RNAiFold and http:
//bioinformatics.bc.edu/clotelab/RNAiFold, which latter has a web server as
well. The software is written in C++ and uses OR-Tools CP search engine.


  Single-cell RNA-Sequencing (scRNA-Seq) is a revolutionary technique for
discovering and describing cell types in heterogeneous tissues, yet its
measurement of expression often suffers from large systematic bias. A major
source of this bias is the cell cycle, which introduces large within-cell-type
heterogeneity that can obscure the differences in expression between cell
types. The current method for removing the cell-cycle effect is unable to
effectively identify this effect and has a high risk of removing other
biological components of interest, compromising downstream analysis. We present
ccRemover, a new method that reliably identifies the cell-cycle effect and
removes it. ccRemover preserves other biological signals of interest in the
data and thus can serve as an important pre-processing step for many scRNA-Seq
data analyses. The effectiveness of ccRemover is demonstrated using simulation
data and three real scRNA-Seq datasets, where it boosts the performance of
existing clustering algorithms in distinguishing between cell types.


  Phytochemicals are plant-derived secondary metabolites, which may exert many
biological activities in humans, including anticancer properties. Although
recent findings appear to support their role in cancer prevention and
treatment, this issue is still controversial. Anti-cancer activity of
phytochemicals mainly depends on their multi-target mechanism of action,
including antimutagenic, antioxidant and antiproliferative activities.
Furthermore, they may modulate the host immune response to cancer, reducing
inflammatory microenvironment and enhancing lymphocyte onco-surveillance. Since
carcinogenesis is multi-factorial and involves several signaling pathways,
multi-targeted agents as phytochemicals may represent promising anticancer
compounds. This narrative review aims to analyze the current literature on
phytochemicals highlighting their specific targets on carcinogenic molecular
pathways and their chemopreventive role. A full comprehension of their activity
at molecular and cellular levels will contribute for a better understanding of
phytochemical clinical efficacy, thus promoting the identification of new
effective plant-derived therapeutics.


  We present an end-to-end genome assembly of a female Aedes aegypti mosquito,
which spreads viral diseases such as yellow fever, dengue, chikungunya, and
Zika to humans. The assembly is based on an earlier genome published in 2007
and improved in 2013. The new assembly has a scaffold N50 of 419Mb, with 96.9%
of the ungapped sequence anchored to chromosomes. We used the new assembly to
examine the conservation of A. aegypti chromosomes. Our results suggest that
synteny is strongly conserved between Ae. aegypti and An. gambiae. Comparison
to D. melanogaster highlights the extent to which the identity of entire
chromosome arms is preserved across dipterans.


  A goal in network science is the geometrical characterization of complex
networks. In this direction, we (arXiv:1603.00386; J. Stat. Mech. (2016)
P063206) have recently introduced the Forman's discretization of Ricci
curvature to the realm of undirected networks. Investigation of Forman
curvature in diverse model and real-world undirected networks revealed that
this measure captures several aspects of the organization of complex undirected
networks. However, many important real-world networks are inherently directed
in nature, and the Forman curvature for undirected networks is unsuitable for
analysis of such directed networks. Hence, we here extend the Forman curvature
for undirected networks to the case of directed networks. The simple
mathematical formula for the Forman curvature in directed networks elegantly
incorporates node weights, edge weights and edge direction. By applying the
Forman curvature for directed networks to a variety of model and real-world
directed networks, we show that the measure can be used to characterize the
structure of complex directed networks. Furthermore, our results also hold in
real directed networks which are weighted or spatial in nature. These results
in combination with our previous results suggest that the Forman curvature can
be readily employed to study the organization of both directed and undirected
complex networks.


  Gene co-expression network differential analysis is designed to help
biologists understand gene expression patterns under different condition. By
comparing different gene co-expression networks we may find conserved part as
well as condition specific set of genes. Taking the network as a collection as
modules, we use a sample-saving method to construct condition-specific gene
co-expression network, and identify differentially expressed subnetworks as
conserved or condition specific modules which may be associated with biological
processes. We have implemented the method as an R package which establishes a
pipeline from expression profile to biological explanations. The usefulness of
the method is also demonstrated by synthetic data as well as Daphnia magna gene
expression data under different environmental stresses.


  We propose a new topological characterization of RNA secondary structures
with pseudoknots based on two topological invariants. Starting from the classic
arc-representation of RNA secondary structures, we consider a model that
couples both I) the topological genus of the graph and II) the number of
crossing arcs of the corresponding primitive graph. We add a term proportional
to these topological invariants to the standard free energy of the RNA
molecule, thus obtaining a novel free energy parametrization which takes into
account the abundance of topologies of RNA pseudoknots observed in RNA
databases.


  Recent results have strongly confirmed that multiple supernovae happened at
distances ~100 pc consisting of two main events: one at 1.7 to 3.2 million
years ago, and the other at 6.5 to 8.7 million years ago. These events are said
to be responsible for excavating the Local Bubble in the interstellar medium
and depositing 60Fe on Earth and the Moon. Other events are indicated by
effects in the local cosmic ray (CR) spectrum. Given this updated and refined
picture, we ask whether such supernovae are expected to have had substantial
effects on the terrestrial atmosphere and biota. In a first cut at the most
probable cases, combining photon and cosmic ray effects, we find that a
supernova at 100 pc can have only a small effect on terrestrial organisms from
visible light and that chemical changes such as ozone depletion are weak.
However, tropospheric ionization right down to the ground due to the
penetration of $\geq$TeV cosmic rays will increase by nearly an order of
magnitude for thousands of years, and irradiation by muons on the ground and in
the upper ocean will increase 20-fold, which will approximately triple the
overall radiation load on terrestrial organisms. Such irradiation has been
linked to possible changes in climate and increased cancer and mutation rates.
This may be related to a minor mass extinction around the Pliocene-Pleistocene
boundary, and further research on the effects is needed.


  The availability of databases electronically encoding curated regulatory
networks and of high-throughput technologies and methods to discover regulatory
interactions provides an invaluable source of data to understand the principles
underpinning the organization and evolution of these networks responsible for
cellular regulation. Nevertheless, data on these sources never goes beyond the
regulon level despite the fact that regulatory networks are complex
hierarchical-modular structures still challenging our understanding. This
brings the necessity for an inventory of systems across a large range of
organisms, a key step to rendering feasible comparative systems biology
approaches. In this work, we take the first step towards a global understanding
of the regulatory networks organization by making a cartography of the
functional architectures of diverse bacteria. Abasy (Across-bacteria systems)
Atlas provides a comprehensive inventory of annotated functional systems,
global network properties, and systems-level elements (global regulators,
modular genes shaping functional systems, basal machinery genes, and
intermodular genes) predicted by the natural decomposition approach for
reconstructed and meta-curated regulatory networks across a large range of
bacteria, including pathogenically and biotechnologically relevant organisms.
The meta-curation of regulatory datasets provides the most complete and
reliable set of regulatory interactions currently available. Abasy Atlas
contains systems and system-level elements for 50 regulatory networks
comprising 78,649 regulatory interactions covering 42 bacteria in nine taxa,
containing 3,708 regulons and 1,776 systems. All this brings together a large
corpus of data that will surely inspire studies to generate hypothesis
regarding the principles governing the evolution and organization of systems
and the functional architectures controlling them.


  Many single-stranded (ss) RNA viruses self assemble from capsid protein
subunits and the nucleic acid to form an infectious virion. It is believed that
the electrostatic interactions between the negatively charged RNA and the
positively charged viral capsid proteins drive the encapsidation, although
there is growing evidence that the sequence of the viral RNA also plays a role
in packaging. In particular the sequence will determine the possible secondary
structures that the ssRNA will take in solution. In this work, we use a mean
field theory to investigate how the secondary structure of the RNA combined
with electrostatic interactions affects the efficiency of assembly and
stability of the assembled virions. We show that the secondary structure of RNA
may result in negative osmotic pressures while a linear polymer causes positive
osmotic pressures for the same conditions. This may suggest that the branched
structure makes the RNA more effectively packaged and the virion more stable.


  This paper gives a systematic method for constructing an N-body potential,
approximating the true potential, that accurately captures meso-scale behavior
of the chemical or biological system using pairwise potentials coming from
experimental data or ab initio methods. The meso-scale behavior is translated
into logic rules for the dynamics. Each pairwise potential has an associated
logic function that is constructed using the logic rules, a class of elementary
logic functions, and AND, OR, and NOT gates. The effect of each logic function
is to turn its associated potential on and off. The N-body potential is
constructed as linear combination of the pairwise potentials, where the
"coefficients" of the potentials are smoothed versions of the associated logic
functions. These potentials allow a potentially low-dimensional description of
complex processes while still accurately capturing the relevant physics at the
meso-scale. We present the proposed formalism to construct coarse-grained
potential models for three examples: an inhibitor molecular system, bond
breaking in chemical reactions, and DNA transcription from biology. The method
can potentially be used in reverse for design of molecular processes by
specifying properties of molecules that can carry them out.


  HIV-1 infection currently cannot be cured because the virus persists as
integrated proviral DNA in long-lived cells despite years of suppressive
antiretroviral therapy (ART). To characterize establishment, turnover, and
evolution of viral DNA reservoirs we deep-sequenced the p17gag region of the
HIV-1 genome from samples obtained after 3-18 years of suppressive ART from 10
patients. For each of these patients, whole genome deep-sequencing data of
HIV-1 RNA populations before onset of ART were available from 6-12 longitudinal
plasma samples spanning 5-8 years of untreated infection. This enabled a
detailed analysis of the dynamics and origin of proviral DNA during ART. A
median of 14% (range 0-42%) of the p17gag DNA sequences were overtly defective
due to G-to-A hypermutation. The remaining sequences were remarkably similar to
previously observed RNA sequences and showed no evidence of evolution over many
years of suppressive ART. Most sequences from the DNA reservoirs were very
similar to viruses actively replicating in plasma (RNA sequences) shortly
before start of ART. The results do not support persistent HIV-1 replication as
a mechanism to maintain the HIV-1 reservoir during suppressive therapy. Rather,
the data indicate that viral DNA variants are turning over as long as patients
are untreated and that suppressive ART halts this turnover.


  Conserved moieties are groups of atoms that remain intact in all reactions of
a metabolic network. Identification of conserved moieties gives insight into
the structure and function of metabolic networks and facilitates metabolic
modelling. All moiety conservation relations can be represented as nonnegative
integer vectors in the left null space of the stoichiometric matrix
corresponding to a biochemical network. Algorithms exist to compute such
vectors based only on reaction stoichiometry but their computational complexity
has limited their application to relatively small metabolic networks. Moreover,
the vectors returned by existing algorithms do not, in general, represent
conservation of a specific moiety with a defined atomic structure. Here, we
show that identification of conserved moieties requires data on reaction atom
mappings in addition to stoichiometry. We present a novel method to identify
conserved moieties in metabolic networks by graph theoretical analysis of their
underlying atom transition networks. Our method returns the exact group of
atoms belonging to each conserved moiety as well as the corresponding vector in
the left null space of the stoichiometric matrix. It can be implemented as a
pipeline of polynomial time algorithms. Our implementation completes in under
five minutes on a metabolic network with more than 4,000 mass balanced
reactions. The scalability of the method enables extension of existing
applications for moiety conservation relations to genome-scale metabolic
networks. We also give examples of new applications made possible by
elucidating the atomic structure of conserved moieties.


  During tissue development, patterns of gene expression determine the spatial
arrangement of cell types. In many cases, gradients of secreted signaling
molecules - morphogens - guide this process. The continuous positional
information provided by the gradient is converted into discrete cell types by
the downstream transcriptional network that responds to the morphogen. A
mechanism commonly used to implement a sharp transition between two adjacent
cell fates is the genetic toggle switch, composed of cross-repressing
transcriptional determinants. Previous analyses emphasize the steady state
output of these mechanisms. Here, we explore the dynamics of the toggle switch
and use exact numerical simulations of the kinetic reactions, the Chemical
Langevin Equation, and Minimum Action Path theory to establish a framework for
studying the effect of gene expression noise on patterning time and boundary
position. This provides insight into the time scale, gene expression
trajectories and directionality of stochastic switching events between cell
states. Taking gene expression noise into account predicts that the final
boundary position of a morphogen-induced toggle switch, although robust to
changes in the details of the noise, is distinct from that of the deterministic
system. Moreover, stochastic switching introduces differences in patterning
time along the morphogen gradient that result in a patterning wave propagating
away from the morphogen source. The velocity of this wave is influenced by
noise; the wave sharpens and slows as it advances and may never reach steady
state in a biologically relevant time. This could explain experimentally
observed dynamics of pattern formation. Together the analysis reveals the
importance of dynamical transients for understanding morphogen-driven
transcriptional networks and indicates that gene expression noise can
qualitatively alter developmental patterning.


  We present a biochemical model of the wall shear stress (WSS)-induced
activation of endothelial nitric oxide synthase (eNOS) in an endothelial cell
(EC). The model includes three key mechanotransducers: mechanosensing ion
channels, integrins and G-protein-coupled receptors. The reaction cascade
consists of two interconnected parts. The first is rapid activation of calcium,
which results in formation of calcium-calmodulin complexes, followed by
recruitment of eNOS from caveolae. The second is phosphoryaltion of eNOS by
protein kinases PKC and AKT. The model also includes a negative feedback loop
due to inhibition of calcium influx into the cell by cyclic guanosine
monophosphate (cGMP). In this feedback, increased nitric oxide (NO) levels
cause an increase in cGMP levels, so that cGMP inhibition of calcium influx can
limit NO production. The model was used to predict the dynamics of NO
production by an EC subjected to a step increase of WSS from zero to a finite
physiologically relevant value. Among several experimentally observed features,
the model predicts a highly nonlinear, biphasic transient behavior of eNOS
activation and NO production: a rapid initial activation due to the very rapid
influx of calcium into the cytosol (occurring within 1 to 5 minutes) is
followed by a sustained period of activation due to protein kinases.


  We previously proposed the existence of DNA to DNA transcription in
eukaryotic cells, but the mechanism by which single-stranded DNA (ssDNA)
transcript is produced and released from the genome remains unknown. We once
speculated that the mechanism of DNA to DNA transcription might be similar to
that of DNA to RNA transcription, but now we propose that endonuclease
dependent transcript cutout may be a possible mechanism of DNA to DNA
transcription, in which a copy of ssDNA fragment (transcript) between two nicks
produced by nicking endonuclease is released from double-stranded DNA (dsDNA)
region in the genome by an unknown ssDNA fragment releasing enzyme. The gap in
the dsDNA will be filled through DNA repair mechanism. Occasionally, multiple
copies of ssDNA transcripts could be produced through multiple rounds of
cutout-repair-cutout cycle.


  The problem of antibiotic resistance poses challenges across many
disciplines. One such challenge is to understand the fundamental science of how
antibiotics work, and how resistance to them can emerge. This is an area where
physicists can make important contributions. Here, we highlight cases where
this is already happening, and suggest directions for further physics
involvement in antimicrobial research.


  During long time, historians and genealogists have interpreted that the elite
that emerged during the Spanish Conquest was almost exclusively European. We
reconstructed a deep matrilineal genealogy which includes recent Costa Rican
ex-presidents and religious authorities back to their ancestors at the early
17th century, and compared their historic ethnic affinities with genetic
mitochondrial evidence of some living descendents. The observed DNA lineage has
an Amerindian ancestry. Such results point out that an Amerindian gene flow had
occurred into the Spanish group during the first generations of colonial
society. This conclusion do not support the current idea that the Spanish elite
avoided interethnic marriages.


  The technological advances in Genetics have given rise to the science of
Molecular Genealogy, by giving us the opportunity to approach the study of
pedigrees from a new perspective, examining alive people at the same time than
collating records from their ancestors. A four generations matrilineal
genealogy is presented from Da. Angela Antunez (b. a. 1881 in Barcelona, Spain)
to their descendants in Cuba and Costa Rica. The observed mitochondrial lineage
has a substitution in the region HVR1 respect to the reference sequence (CRS)
and belongs to H haplogroup. This is consistent -as expected, with the recent
documented European lineage. This is an example of differential female
migration, in which women move to the husband's place of origin. Because of
patrilocal marriages are practiced in most societies, such migration mechanism
has had a major impact since the arrival of immigrant women bring new and
different lineages of mtDNA, which together with the admixture process
contributed to enrich the diversity and regionalization of the gene pool of
human populations in the Americas.


  Traditionally the Costa Rican historians and genealogists have interpreted
that the Spanish ruling elite emerged after the conquest was exclusively of
European origin. On the other side, recent technological advances in Genetics
give us the opportunity to approach the study of pedigrees from a new
perspective, examining alive people and simultaneously collating the historical
information of their ancestors. In this paper, a complete matrilineal genealogy
was reconstructed from nowadays "white" Costa Ricans to their ancestors in the
early Colonial society (XVI century). It was compared the correlation between
ethnic affiliations deduced from historical records with the genetic
inheritance from maternal lineages. The MtDNA lineage observed corresponds to a
Native American ancestry. These results show that some Amerindian gene flow
into the Spanish group must have occurred since the first generation of
Colonial society, a finding that contrasts with the prevailing ideas that the
Spanish elite avoided the intermarriage with other ethnic groups. Examples like
this one confirm that miscegenation began early in Costa Rica. So, those who
considered themselves "Spaniards" in the late colonial era, were actually
biologically mestizos. It is widely accepted, that the general Costa Rican
population is the result of an admixture process between Europeans, Amerindians
and Africans.


  Seed banks are a common characteristics to many plant species, which allow
storage of genetic diversity in the soil as dormant seeds for various periods
of time. We investigate an above-ground population following a Fisher-Wright
model with selection coupled with a deterministic seed bank assuming the length
of the seed bank is kept constant and the number of seeds is large. To assess
the combined impact of seed banks and selection on genetic diversity, we derive
a general diffusion model. The applied techniques outline a path of
approximating a stochastic delay differential equation by an appropriately
rescaled stochastic differential equation, which is a common issue in
statistical physics. We compute the equilibrium solution of the site-frequency
spectrum and derive the times to fixation of an allele with and without
selection. Finally, it is demonstrated that seed banks enhance the effect of
selection onto the site-frequency spectrum while slowing down the time until
the mutation-selection equilibrium is reached.


  A decrease in system size can induce qualitatively different behavior
compared to the macroscopic behavior of the corresponding large-size system.
The mechanisms of this transition, which is known as the small-size transition,
can be attributed to either a relative increase in the noise intensity or to
the discreteness of the state space due to the small system size. The former
mechanism has been intensively investigated using several toy and realistic
models. However, the latter has rarely been analyzed and is sometimes confused
with the former, because a toy model that extracts the essence of the
discreteness-induced transition mechanism is lacking. In this work, we propose
a 1- and 3-body reaction system as a minimal model of the discreteness-induced
transition and derive the conditions under which this transition occurs in more
complex systems. This work enriches our understanding of the influence of small
system size on system behavior.


  In this work, some phenomenological models, those that are based only on the
population information (macroscopic level), are deduced in an intuitive way.
These models, as for instance Verhulst, Gompertz and Bertalanffy models, are
posted in such a manner that all the parameters involved have physical
interpretation. A model based on the interaction (distance dependent) between
the individuals (microscopic level) is also presented. This microscopic model
reachs the phenomenological models presented as particular cases. In this
approach, the Verhulst model represents the situation in which all the
individuals interact in the same way, regardless of the distance between them.
That means Verhulst model is a kind of mean field model. The other
phenomenological models are reaching from the microscopic model according to
two quantities: i) the relation between the way that the interaction decays
with the distance; and ii) the dimension of the spatial structure formed by the
population. This microscopic model allows understanding population growth by
first principles, because it predicts that some phenomenological models can be
seen as a consequence of the same individual level kind of interaction. In this
sense, the microscopic model that will be discussed here paves the way to
finding universal patterns that are common to all types of growth, even in
systems of very different nature.


  In structure-based models of proteins, one often assumes that folding is
accomplished when all contacts are established. This assumption may frequently
lead to a conceptual problem that folding takes place in a temperature region
of very low thermodynamic stability, especially when the contact map used is
too sparse. We consider six different structure-based models and show that
allowing for a small, but model-dependent, percentage of the native contacts
not being established boosts the folding temperature substantially while
affecting the time scales of folding only in a minor way. We also compare other
properties of the six models. We show that the choice of the description of the
backbone stiffness has a substantial effect on the values of characteristic
temperatures that relate both to equilibrium and kinetic properties. Models
without any backbone stiffness (like the self-organized polymer) are found to
perform similar to those with the stiffness, including in the studies of
stretching.


  Dengue is a vector-borne disease transmitted by the mosquito Aedes aegypti.
It has been observed that its incidence is strongly influenced by temperature
and other abiotic factors like rainfall and humidity. In this work we compare
the effects of seasonality in temperature that affect the entomological
parameters of the mosquito with precipitation that affects the hatchery
capacity. We also analyze its joint action using a dynamical model for the life
cycle of Aedes aegypti and historical weather data from 8 regions of Mexico. We
found that the joint action of different mechanisms can enhance the prevalence
of the disease, but also inhibit it when they act in an asynchronous way. We
found that for the studied regions, the seasonality of the low temperature
rather than mean temperature is the main driving force of Dengue outbreaks. We
also analyzed the role of the diapause in these kinds of outbreaks. The
methodology developed here can be used to discover the underlying mechanism of
Dengue outbreaks in different regions and thus help to apply targeted control
measures.


  Molecular motor proteins serve as an essential component of intracellular
transport by generating forces to haul cargoes along cytoskeletal filaments.
Two species of motors that are directed oppositely (e.g. kinesin, dynein) can
be attached to the same cargo, which is known to produce bidirectional net
motion. Although previous work focuses on the motor number as the driving noise
source for switching, we propose an alternative mechanism: cargo diffusion. A
mean-field mathematical model of mechanical interactions of two populations of
molecular motors with cargo thermal fluctuations (diffusion) is presented to
study this phenomenon. The delayed response of a motor to fluctuations in the
cargo velocity is quantified, allowing for the reduction of the full model a
single "characteristic distance", a proxy for the net force on the cargo. The
system is then found to be metastable, with switching exclusively due to cargo
diffusion between distinct directional transport states. The time to switch
between these states is then investigated using a mean first passage time
analysis. The switching time is found to be non-monotonic in the drag of the
cargo, providing an experimental test of the theory.


  Evolutionary games on networks traditionally involve the same game at each
interaction. Here we depart from this assumption by considering mixed games,
where the game played at each interaction is drawn uniformly at random from a
set of two different games. While in well-mixed populations the random mixture
of the two games is always equivalent to the average single game, in structured
populations this is not always the case. We show that the outcome is in fact
strongly dependent on the distance of separation of the two games in the
parameter space. Effectively, this distance introduces payoff heterogeneity,
and the average game is returned only if the heterogeneity is small. For higher
levels of heterogeneity the distance to the average game grows, which often
involves the promotion of cooperation. The presented results support preceding
research that highlights the favorable role of heterogeneity regardless of its
origin, and they also emphasize the importance of the population structure in
amplifying facilitators of cooperation.


  Evolution occurs in populations of reproducing individuals. The structure of
a biological population affects which traits evolve. Understanding evolutionary
game dynamics in structured populations is difficult. Precise results have been
absent for a long time, but have recently emerged for special structures where
all individuals have the same number of neighbors. But the problem of
determining which trait is favored by selection in the natural case where the
number of neighbors can vary, has remained open. For arbitrary selection
intensity, the problem is in a computational complexity class which suggests
there is no efficient algorithm. Whether there exists a simple solution for
weak selection was unanswered. Here we provide, surprisingly, a general formula
for weak selection that applies to any graph or social network. Our method uses
coalescent theory and relies on calculating the meeting times of random walks.
We can now evaluate large numbers of diverse and heterogeneous population
structures for their propensity to favor cooperation. We can also study how
small changes in population structure---graph surgery---affect evolutionary
outcomes. We find that cooperation flourishes most in societies that are based
on strong pairwise ties.


  Motivated by a non-random but clustered distribution of SNPs, we introduce a
phenomenological model to account for the clustering properties of SNPs in the
human genome. The phenomenological model is based on a preferential mutation to
the closer proximity of existing SNPs. With the Hapmap SNP data, we empirically
demonstrate that the preferential model is better for illustrating the
clustered distribution of SNPs than the random model. Moreover, the model is
applicable not only to autosomes but also to the X chromosome, although the X
chromosome has different characteristics from autosomes. The analysis of the
estimated parameters in the model can explain the pronounced population
structure and the low genetic diversity of the X chromosome. In addition,
correlation between the parameters reveals the population-wise difference of
the mutation probability. These results support the mutational non-independence
hypothesis against random mutation.


  Mathematical models of cholera and waterborne disease vary widely in their
structures, in terms of transmission pathways, loss of immunity, and other
features. These differences may yield different predictions and parameter
estimates from the same data. Given the increasing use of models to inform
public health decision-making, it is important to assess distinguishability
(whether models can be distinguished based on fit to data) and inference
robustness (whether model inferences are robust to realistic variations in
model structure). We examined the effects of uncertainty in model structure in
epidemic cholera, testing a range of models based on known features of cholera
epidemiology. We fit to simulated epidemic and long-term data, as well as data
from the 2006 Angola epidemic. We evaluated model distinguishability based on
data fit, and whether parameter values and forecasts can accurately be inferred
from incidence data. In general, all models were able to successfully fit to
all data sets, even if misspecified. However, in the long-term data, the best
model fits were achieved when the loss of immunity form matched those of the
model that simulated the data. Two transmission and reporting parameters were
accurately estimated across all models, while the remaining showed broad
variation across the different models and data sets. Forecasting efforts were
not successful early, but once the epidemic peak had been achieved, most models
showed similar accuracy. Our results suggest that we are unlikely to be able to
infer mechanistic details from epidemic case data alone, underscoring the need
for broader data collection. Nonetheless, with sufficient data, conclusions
from forecasting and some parameter estimates were robust to variations in the
model structure, and comparative modeling can help determine how variations in
model structure affect conclusions drawn from models and data.


  Maintenance of sexual reproduction and genetic recombination imposes
physiological costs when compared to parthenogenic reproduction, most
prominently: for maintaining the corresponding (molecular) machinery, for
finding a mating partner, and through the decreased fraction of females in a
population, which decreases the reproductive capacity. Based on principles from
information theory, we have previously developed a new population genetic
model, and applying it in simulations, we have recently hypothesized that all
species maintain the maximum genomic complexity that is required by their niche
and allowed by their mutation rate and selection intensity. Applying this idea
to the complexity overhead of recombination maintenance, its costs must be more
than compensated by an additional capacity for complexity in recombining
populations. Here, we show a simple mechanism, where recombination helps to
maintain larger biases of alleles frequencies in a population, so the
advantageous alleles can have increased frequency. This allows recombining
populations to maintain higher fitness and phenotypic efficiency in comparison
with asexual populations with the same parameters. Random mating alone already
significantly increases the ability to maintain genomic and phenotypic
complexity. Sexual selection provides additional capacity for this complexity.
The model can be considered as a unifying synthesis of previous hypotheses
about the roles of recombination in Muller's ratchet, mutation purging and Red
Queen dynamics, because the introduction of recombination both increases
population frequencies of beneficial alleles and decreases detrimental ones. In
addition, we suggest simple explanations for niche-dependent prevalence of
transient asexuality and the exceptional asexual lineage of Bdelloid rotifers.


  Precise positioning of the mitotic spindle is important for specifying the
plane of cell division, which in turn determines how the cytoplasmic contents
are partitioned into the daughter cells, and how the daughters are positioned
within the tissue. During metaphase in the early C. elegans embryo, the spindle
is aligned and centered on the anterior-posterior axis by a
microtubule-dependent machinery that exerts restoring forces when the spindle
is displaced from the center. To investigate the accuracy and stability of
centering, we tracked the position and orientation of the mitotic spindle
during the first cell division with high temporal and spatial resolution. We
found that the precision is remarkably high: the cell-to-cell variation in the
transverse position of the center of the spindle during metaphase, as measured
by the standard deviation, was only 1.5% of the length of the short axis of the
cell. Spindle position is also very stable: the standard deviation of the
fluctuations in transverse spindle position during metaphase was only 0.5% of
the short axis of the cell. Assuming that stability is limited by fluctuations
in the number of independent motor elements such as microtubules or dyneins
underlying the centering machinery, we infer that the number is on the order of
one thousand, consistent with the several thousand of astral microtubules in
these cells. Astral microtubules grow out from the two spindle poles, make
contact with the cell cortex, and then shrink back shortly thereafter. The high
stability of centering can be accounted for quantitatively if, while making
contact with the cortex, the astral microtubules buckle as they exert
compressive, pushing forces. We thus propose that the large number of
microtubules in the asters provides a highly precise mechanism for positioning
the spindle during metaphase while assembly is completed prior to the onset of
anaphase.


  Selecting regularization parameters in penalized high-dimensional graphical
models in a principled, data-driven, and computationally efficient manner
continues to be one of the key challenges in high-dimensional statistics. We
present substantial computational gains and conceptual generalizations of the
Stability Approach to Regularization Selection (StARS), a state-of-the-art
graphical model selection scheme. Using properties of the Poisson-Binomial
distribution and convex non-asymptotic distributional modeling we propose lower
and upper bounds on the StARS graph regularization path which results in
greatly reduced computational cost without compromising regularization
selection. We also generalize the StARS criterion from single edge to induced
subgraph (graphlet) stability. We show that simultaneously requiring edge and
graphlet stability leads to superior graph recovery performance independent of
graph topology. These novel insights render Gaussian graphical model selection
a routine task on standard multi-core computers.


  We use a recently proposed method called Spectral Gap Optimization of Order
Parameters (SGOOP) (Tiwary and Berne, Proc. Natl. Acad. Sci 2016, 113, 2839
(2016)), to determine an optimal 1-dimensional reaction coordinate (RC) for the
unbinding of a bucky-ball from a pocket in explicit water. This RC is estimated
as a linear combination of the multiple available order parameters that
collectively can be used to distinguish the various stable states relevant for
unbinding. We pay special attention to determining and quantifying the degree
to which water molecules should be included in the RC. Using SGOOP with
under-sampled biased simulations, we predict that water plays a distinct role
in the reaction coordinate for unbinding in the case when the ligand is
sterically constrained to move along an axis of symmetry. This prediction is
validated through extensive calculations of the unbinding times through
metadynamics, and by comparison through detailed balance with unbiased
molecular dynamics estimate of the binding time. However when the steric
constraint is removed, we find that the role of water in the reaction
coordinate diminishes. Here instead SGOOP identifies a good one-dimensional RC
involving various motional degrees of freedom.


  Gene transcription is a highly stochastic and dynamic process. As a result,
the mRNA copy number of a given gene is heterogeneous both between cells and
across time. We present a framework to model gene transcription in populations
of cells with time-varying (stochastic or deterministic) transcription and
degradation rates. Such rates can be understood as upstream cellular drives
representing the effect of different aspects of the cellular environment. We
show that the full solution of the master equation contains two components: a
model-specific, upstream effective drive, which encapsulates the effect of
cellular drives (e.g., entrainment, periodicity or promoter randomness), and a
downstream transcriptional Poissonian part, which is common to all models. Our
analytical framework treats cell-to-cell and dynamic variability consistently,
unifying several approaches in the literature. We apply the obtained solution
to characterise different models of experimental relevance, and to explain the
influence on gene transcription of synchrony, stationarity, ergodicity, as well
as the effect of time-scales and other dynamic characteristics of drives. We
also show how the solution can be applied to the analysis of noise sources in
single-cell data, and to reduce the computational cost of stochastic
simulations.


  Human immunodeficiency virus (HIV) capsid proteins spontaneously assemble
around the genome into a protective protein shell called the capsid, which can
take on a variety of shapes broadly classified as conical, cylindrical and
irregular. The majority of capsids seen in in vivo studies are conical in
shape, while in vitro experiments have shown a preference for cylindrical
capsids. The factors involved in the selection of the unique shape of HIV
capsids are not well understood, and in particular the impact of RNA on the
formation of the capsid is not known. In this work, we study the role of the
genome and its interaction with the capsid protein by modeling the genomic RNA
through a mean-field theory. Our results show that the confinement free energy
for a homopolymeric model genome confined in a conical capsid is lower than
that in a cylindrical capsid, at least when the genome does not interact with
the capsid, which seems to be the case in in vivo experiments. Conversely, the
confinement free energy for the cylinder is lower than for a conical capsid if
the genome is attracted to the capsid proteins as the in vitro experiments.
Understanding the factors that contribute to the formation of conical capsids
may shed light on the infectivity of HIV particles.


  Initial protein structural comparisons were sequence-based. Since amino acids
that are distant in the sequence can be close in the 3-dimensional (3D)
structure, 3D contact approaches can complement sequence approaches.
Traditional 3D contact approaches study 3D structures directly. Instead, 3D
structures can be modeled as protein structure networks (PSNs). Then, network
approaches can compare proteins by comparing their PSNs. Network approaches may
improve upon traditional 3D contact approaches. We cannot use existing PSN
approaches to test this, because: 1) They rely on naive measures of network
topology. 2) They are not robust to PSN size. They cannot integrate 3) multiple
PSN measures or 4) PSN data with sequence data, although this could help
because the different data types capture complementary biological knowledge. We
address these limitations by: 1) exploiting well-established graphlet measures
via a new network approach, 2) introducing normalized graphlet measures to
remove the bias of PSN size, 3) allowing for integrating multiple PSN measures,
and 4) using ordered graphlets to combine the complementary PSN data and
sequence data. We compare both synthetic networks and real-world PSNs more
accurately and faster than existing network, 3D contact, or sequence
approaches. Our approach finds PSN patterns that may be biochemically
interesting.


  One of the most intriguing results of single molecule experiments on proteins
and nucleic acids is the discovery of functional heterogeneity: the observation
that complex cellular machines exhibit multiple, biologically active
conformations. The structural differences between these conformations may be
subtle, but each distinct state can be remarkably long-lived, with random
interconversions between states occurring only at macroscopic timescales,
fractions of a second or longer. Though we now have proof of functional
heterogeneity in a handful of systems---enzymes, motors, adhesion
complexes---identifying and measuring it remains a formidable challenge. Here
we show that evidence of this phenomenon is more widespread than previously
known, encoded in data collected from some of the most well-established single
molecule techniques: AFM or optical tweezer pulling experiments. We present a
theoretical procedure for analyzing distributions of rupture/unfolding forces
recorded at different pulling speeds. This results in a single parameter,
quantifying the degree of heterogeneity, and also leads to bounds on the
equilibration and conformational interconversion timescales. Surveying ten
published datasets, we find heterogeneity in five of them, all with
interconversion rates slower than 10 s$^{-1}$. Moreover, we identify two
systems where additional data at realizable pulling velocities is likely to
find a theoretically predicted, but so far unobserved cross-over regime between
heterogeneous and non-heterogeneous behavior. The significance of this regime
is that it will allow far more precise estimates of the slow conformational
switching times, one of the least understood aspects of functional
heterogeneity.


  Delayed processes are ubiquitous in biological systems and are often
characterized by delay differential equations (DDEs) and their extension to
include stochastic effects. DDEs do not explicitly incorporate intermediate
states associated with a delayed process but instead use an estimated average
delay time. In an effort to examine the validity of this approach, we study
systems with significant delays by explicitly incorporating intermediate steps.
We show by that such explicit models often yield significantly different
equilibrium distributions and transition times as compared to DDEs with
deterministic delay values. Additionally, different explicit models with
qualitatively different dynamics can give rise to the same DDEs revealing
important ambiguities. We also show that DDE-based predictions of oscillatory
behavior may fail for the corresponding explicit model.


  High-diversity assemblages are very common in nature, and yet the factors
allowing for the maintenance of biodiversity remain obscure. The competitive
exclusion principle and May's complexity-diversity puzzle both suggest that a
community can support only a small number of species, turning the spotlight at
the dynamics of local patches or islands, where stable and uninvadable (SU)
subsets of species play a crucial role. Here we map the community SUs question
to the geometric problem of finding maximal cliques of the corresponding graph.
We solve for the number of SUs as a function of the species richness in the
regional pool, $N$, showing that this growth is subexponential, contrary to
long-standing wisdom. We show that symmetric systems relax rapidly to an SU,
where the system stays until a regime shift takes place. In asymmetric systems
the relaxation time grows much faster with $N$, suggesting an excitable
dynamics under noise.


  Gene expression is a complex and quantitative trait that is influenced by
both genetic and non-genetic regulators including environmental factors.
Evaluating the contribution of environment to gene expression regulation and
identifying which genes are more likely to be influenced by environmental
factors are important for understanding human complex traits. We hypothesize
that by living together as couples, there can be commonly co-regulated genes
that may reflect the shared living environment (e.g., diet, indoor air
pollutants, behavioral lifestyle). The lymphoblastoid cell lines (LCLs) derived
from unrelated couples of African ancestry (YRI, Yoruba people from Ibadan,
Nigeria) from the International HapMap Project provided a unique model for us
to characterize gene expression pattern in couples by comparing gene expression
levels between husbands and wives. Strikingly, 778 genes were found to show
much smaller variances in couples than random pairs of individuals at a false
discovery rate (FDR) of 5%. Since genetic variation between unrelated family
members in a general population is expected to be the same assuming a
random-mating society, non-genetic factors (e.g., epigenetic systems) are more
likely to be the mediators for the observed transcriptional similarity in
couples. We thus evaluated the contribution of modified cytosines to those
genes showing transcriptional similarity in couples as well as the
relationships these CpG sites with other gene regulatory elements, such as
transcription factor binding sites (TFBS). Our findings suggested that
transcriptional similarity in couples likely reflected shared common
environment partially mediated through cytosine modifications.


  Social dilemmas are an integral part of social interactions. Cooperative
actions, ranging from secreting extra-cellular products in microbial
populations to donating blood in humans, are costly to the actor and hence
create an incentive to shirk and avoid the costs. Nevertheless, cooperation is
ubiquitous in nature. Both costs and benefits often depend non-linearly on the
number and types of individuals involved -- as captured by idioms such as `too
many cooks spoil the broth' where additional contributions are discounted, or
`two heads are better than one' where cooperators synergistically enhance the
group benefit. Interaction group sizes may depend on the size of the population
and hence on ecological processes. This results in feedback mechanisms between
ecological and evolutionary processes, which jointly affect and determine the
evolutionary trajectory. Only recently combined eco-evolutionary processes
became experimentally tractable in microbial social dilemmas. Here we analyse
the evolutionary dynamics of non-linear social dilemmas in settings where the
population fluctuates in size and the environment changes over time. In
particular, cooperation is often supported and maintained at high densities
through ecological fluctuations. Moreover, we find that the combination of the
two processes routinely reveals highly complex dynamics, which suggests common
occurrence in nature.


  The binding of molecules, ions or proteins to specific target sites is a
generic step for cell activation. However, this step relies on rare events
where stochastic particles located in a large bulk are searching for small and
often hidden targets and thus remains difficult to study. We present here a
hybrid discrete-continuum model where the large ensemble of particles is
described by mass-action laws. The rare discrete binding events are modeled by
a Markov chain for the encounter of a finite number of small targets by few
Brownian particles, for which the arrival time is Poissonian. This model is
applied for predicting the time distribution of vesicular release at neuronal
synapses that remains elusive. This release is triggered by the binding of few
calcium ions that can originate either from the synaptic bulk or from the
transient entry through calcium channels. We report that the distribution of
release time is bimodal although triggered by a single fast action potential:
while the first peak follows a stimulation, the second corresponds to the
random arrival over much longer time of ions located in the bulk to small
binding targets. To conclude, the present multiscale stochastic chemical
reaction modeling allows studying cellular events based on integrating discrete
molecular events over various time scales.


  Evolutionary game dynamics are often studied in the context of different
population structures. Here we propose a new population structure that is
inspired by simple multicellular life forms. In our model, cells reproduce but
can stay together after reproduction. They reach complexes of a certain size,
n, before producing single cells again. The cells within a complex derive
payoff from an evolutionary game by interacting with each other. The
reproductive rate of cells is proportional to their payoff. We consider all
two-strategy games. We study deterministic evolutionary dynamics with
mutations, and derive exact conditions for selection to favor one strategy over
another. Our main result has the same symmetry as the well-known sigma
condition, which has been proven for stochastic game dynamics and weak
selection. For a maximum complex size of n=2 our result holds for any intensity
of selection. For n > 2 it holds for weak selection. As specific examples we
study the prisoner's dilemma and hawk-dove games. Our model advances
theoretical work on multicellularity by allowing for frequency-dependent
interactions within groups.


  We present a maximum-caliber method for inferring transition rates of a
Markov State Model (MSM) with perturbed equilibrium populations, given
estimates of state populations and rates for an unperturbed MSM. It is similar
in spirit to previous approaches but given the inclusion of prior information
it is more robust and simple to implement. We examine its performance in simple
biased diffusion models of kinetics, and then apply the method to predicting
changes in folding rates for several highly non-trivial protein folding systems
for which non-native interactions play a significant role, including (1)
tryptophan variants of GB1 hairpin, (2) salt-bridge mutations of Fs peptide
helix, and (3) MSMs built from ultra-long folding trajectories of FiP35 and GTT
variants of WW domain. In all cases, the method correctly predicts changes in
folding rates, suggesting the wide applicability of maximum-caliber approaches
to efficiently predict how mutations perturb protein conformational dynamics.


  Antibiotics such as the quinolones and fluoroquinolones kill bacterial
pathogens ultimately through DNA damage. They target the essential type IIA
topoisomerases in bacteria by stabilising the normally transient double strand
break state which is created to modify the supercoiling state of the DNA. Here
we discuss the development of these antibiotics and their method of action.
Existing methods for DNA damage visualisation, such as the comet assay and
immunofluorescence imaging can often only be analysed qualitatively and this
analysis is subjective. We describe a putative single-molecule fluorescence
technique for quantifying DNA damage via the total fluorescence intensity of a
DNA origami tile fully saturated with an intercalating dye, along with the
optical requirements for how to implement these into a light microscopy imaging
system capable of single-molecule millisecond timescale imaging. This system
promises significant improvements in reproducibility of the quantification of
DNA damage over traditional techniques.


  Biophysical force spectroscopy tools - for example optical tweezers, magnetic
tweezers, atomic force microscopy, - have been used to study elastic,
mechanical, conformational and dynamic properties of single biological
specimens from single proteins to whole cells to reveal information not
accessible by ensemble average methods such as X-ray crystallography, mass
spectroscopy, gel electrophoresis and so on. Here we review the application of
these tools on a range of infection-related questions from antibody-inhibited
protein processivity to virus-cell adhesion. In each case we focus on how the
instrumental design tailored to the biological system in question translates
into the functionality suitable for that particular study. The unique insights
that force spectroscopy has gained to complement knowledge learned through
population averaging techniques in interrogating biomolecular details prove to
be instrumental in therapeutic innovations such as those in structure-based
drug design.


  The method of action of many antibiotics is to interfere with DNA replication
- quinolones trap DNA gyrase and topoisomerase proteins onto DNA while
metronidazole causes single and double stranded breaks in DNA. To understand
how bacteria respond to these drugs, it is important to understand the repair
processes utilised when DNA replication is blocked. We have used tandem lac
operators inserted into the chromosome bound by fluorescently labelled lac
repressors as a model protein block to replication in E. coli. We have used
dual-colour, alternating-laser, single-molecule narrowfield microscopy to
quantify the amount of operator at the block and simultaneously image
fluorescently labelled DNA polymerase. We anticipate use of this system as a
quantitative platform to study replication stalling and repair proteins.


  The SMC complex, MukBEF, is important for chromosome organization and
segregation in Escherichia coli. Fluorescently tagged MukBEF forms distinct
spots (or 'foci') in the cell, where it is thought to carry out most of its
chromosome associated activities. This chapter outlines the technique of
Fluorescence Recovery After Photobleaching (FRAP) as a method to study the
properties of YFP-tagged MukB in fluorescent foci. This method can provide
important insight into the dynamics of MukB on DNA and be used to study its
biochemical properties in vivo.


  Single-molecule narrowfield microscopy is a versatile tool to investigate a
diverse range of protein dynamics in live cells and has been extensively used
in bacteria. Here, we describe how these methods can be extended to larger
eukaryotic, yeast cells, which contain sub-cellular compartments. We describe
how to obtain single-molecule microscopy data but also how to analyse these
data to track and obtain the stoichiometry of molecular complexes diffusing in
the cell. We chose glucose mediated signal transduction of live yeast cells as
the system to demonstrate these single-molecule techniques as transcriptional
regulation is fundamentally a single molecule problem - a single repressor
protein binding a single binding site in the genome can dramatically alter
behaviour at the whole cell and population level.


  We consider the following frustrated optimization problem: given a prior
probability distribution $q$, find the distribution $p$ minimizing the relative
entropy with respect to $q$ such that $\textrm{mean}(p)$ is fixed and large. We
show that solutions to this problem are asymptotically Gaussian. As an
application we derive an $H$-type theorem for evolutionary dynamics: the
entropy of the (standardized) distribution of fitness of a population evolving
under natural selection is eventually increasing.


  Genomic and metagenomic fields, generating huge sets of short genomic
sequences, brought their own share of high performance problems. To extract
relevant pieces of information from the huge data sets generated by current
sequencing techniques, one must rely on extremely scalable methods and
solutions. Indexing billions of objects is a task considered too expensive
while being a fundamental need in this field. In this paper we propose a
straightforward indexing structure that scales to billions of element and we
propose two direct applications in genomics and metagenomics. We show that our
proposal solves problem instances for which no other known solution scales-up.
We believe that many tools and applications could benefit from either the
fundamental data structure we provide or from the applications developed from
this structure.


  What can we learn about controlling a system solely from its underlying
network structure? Here we adapt a recently developed framework for control of
networks governed by a broad class of nonlinear dynamics that includes the
major dynamic models of biological, technological, and social processes. This
feedback-based framework provides realizable node overrides that steer a system
towards any of its natural long term dynamic behaviors, regardless of the
specific functional forms and system parameters. We use this framework on
several real networks, identify the topological characteristics that underlie
the predicted node overrides, and compare its predictions to those of
structural controllability in control theory. Finally, we demonstrate this
framework's applicability in dynamic models of gene regulatory networks and
identify nodes whose override is necessary for control in the general case, but
not in specific model instances.


  A theory for qualitative models of gene regulatory networks has been
developed over several decades, generally considering transcription factors to
regulate directly the expression of other transcription factors, without any
intermediate variables. Here we explore a class of models that explicitly
includes both transcription and translation, keeping track of both mRNA and
protein concentrations. We mainly deal with transcription regulation functions
that are steep sigmoids or step functions, as is often done in protein-only
models, though translation is governed by a linear term. We extend many aspects
of the protein-only theory to this new context, including properties of fixed
points, description of trajectories by mappings between switching points,
qualitative analysis via a state-transition diagram, and a result on periodic
orbits for negative feedback loops. We find that while singular behaviour in
switching domains is largely avoided, non-uniqueness of solutions can still
occur in the step-function limit.


  Multilevel selection is an important organizing principle that crucially
underlies evolutionary processes from the emergence of cells to eusociality and
the economics of nations. Previous studies on multilevel selection assumed that
the effective higher-level selection emerges from lower-level reproduction.
This leads to selection among groups, although only individuals reproduce. We
introduce selective group extinction, where groups die with a probability
inversely proportional to their group fitness. When accounting for this the
critical benefit-to-cost ratio is substantially lowered. Because in game theory
and evolutionary dynamics the degree of cooperation crucially depends on this
ratio above which cooperation emerges previous studies may have substantially
underestimated the establishment and maintenance of cooperation.


  We explore the emergence of persistent infection in a patch of population,
where the disease progression of the individuals is given by the SIRS model and
an individual becomes infected on contact with another infected individual. We
investigate the persistence of contagion qualitatively and quantitatively,
under varying degrees of heterogeneity in the initial population. We observe
that when the initial population is uniform, consisting of individuals at the
same stage of disease progression, infection arising from a contagious seed
does not persist. However when the initial population consists of randomly
distributed refractory and susceptible individuals, a single source of
infection can lead to sustained infection in the population, as heterogeneity
facilitates the de-synchronization of the phases in the disease cycle of the
individuals. We also show how the average size of the window of persistence of
infection depends on the degree of heterogeneity in the initial composition of
the population. In particular, we show that the infection eventually dies out
when the entire initial population is susceptible, while even a few
susceptibles among an heterogeneous refractory population gives rise to a large
persistent infected set.


  We present a simple and robust technique to extract kinetic rate models and
thermodynamic quantities from single molecule time traces. SMACKS (Single
Molecule Analysis of Complex Kinetic Sequences) is a maximum likelihood
approach that works equally well for long trajectories as for a set of short
ones. It resolves all statistically relevant rates and also their
uncertainties. This is achieved by optimizing one global kinetic model based on
the complete dataset, while allowing for experimental variations between
individual trajectories. In particular, neither a priori models nor equilibrium
have to be assumed. The power of SMACKS is demonstrated on the kinetics of the
multi-domain protein Hsp90 measured by smFRET (single molecule F\"orster
resonance energy transfer). Experiments in and out of equilibrium are analyzed
and compared to simulations, shedding new light on the role of Hsp90's ATPase
function. SMACKS pushes the boundaries of single molecule kinetics far beyond
current methods.


  Population genetics struggles to model extinction; standard models track the
relative rather than absolute fitness of genotypes, while the exceptions
describe only the short-term transition from imminent doom to evolutionary
rescue. But extinction can result from failure to adapt not only to
catastrophes, but also to a backlog of environmental challenges. We model
long-term evolution to long series of small challenges, where fitter
populations reach higher population sizes. The population's long-term fitness
dynamic is well approximated by a simple stochastic Markov chain model.
Long-term persistence occurs when the rate of adaptation exceeds the rate of
environmental deterioration for some genotypes. Long-term persistence times are
consistent with typical fossil species persistence times of several million
years. Immediately preceding extinction, fitness declines rapidly, appearing as
though a catastrophe disrupted a stably established population, even though
gradual evolutionary processes are responsible. New populations go through an
establishment phase where, despite being demographically viable, their
extinction risk is elevated. Should the population survive long enough,
extinction risk later becomes constant over time.


  Background: A deterministic model is developed for the spatial spread of an
epidemic disease in a geographical setting. The disease is borne by vectors to
susceptible hosts through criss-cross dynamics. The model is focused on an
epidemic outbreak that initiates from a small number of cases in a small
sub-region of the geographical setting. Methods: Partial differential equations
are formulated to describe the interaction of the model compartments. Results:
The partial differential equations of the model are analyzed and proven to be
well-posed. The epidemic outcomes of the model are correlated to the spatially
dependent parameters and initial conditions of the model. Conclusions: A
version of the model is applied to the 2015-2016 Zika outbreak in the Rio de
Janeiro Municipality in Brazil.


  Spatial big data have the "velocity," "volume," and "variety" of big data
sources and additional geographic information about the record. Digital data
sources, such as medical claims, mobile phone call data records, and geo-tagged
tweets, have entered infectious disease epidemiology as novel sources of data
to complement traditional infectious disease surveillance. In this work, we
provide examples of how spatial big data have been used thus far in
epidemiological analyses and describe opportunities for these sources to
improve public health coordination and disease mitigation strategies. In
addition, we consider the technical, practical, and ethical challenges with the
use of spatial big data in infectious disease surveillance and inference.
Finally, we discuss the implications of the rising use of spatial big data in
epidemiology to health risk communications, across-scale public health
coordination, and public health policy recommendation.


  Stochastic Process Model has many applications in analysis of longitudinal
biodemographic data. Such data contain various physiological variables
(sometimes known as covariates). It also can potentially contain genetic
information available for all or a part of participants. Taking advantage from
both genetic and non-genetic information can provide future insights into a
broad range of processes describing aging-related changes in the organism. In
this paper, we implemented a multi-dimensional Genetic Stochastic Process Model
(GenSPM) in newly developed software tool, R-package stpm, which allows
researchers performing such kind of analysis.


  In genome-wide association studies (GWASs) of common diseases/traits, we
often analyze multiple GWASs with the same phenotype together to discover
associated genetic variants with higher power. Since it is difficult to access
data with detailed individual measurements, summary-statistics-based
meta-analysis methods have become popular to jointly analyze data sets from
multiple GWASs. In this paper, we propose a novel summary-statistics-based
joint analysis method based on controlling the joint local false discovery rate
(Jlfdr). We prove that our method is the most powerful summary-statistics-based
joint analysis method when controlling the false discovery rate at a certain
level. In particular, the Jlfdr-based method achieves higher power than
commonly used meta-analysis methods when analyzing heterogeneous data sets from
multiple GWASs. Simulation experiments demonstrate the superior power of our
method over meta-analysis methods. Also, our method discovers more associations
than meta-analysis methods from empirical data sets of four phenotypes. The
R-package is available at: http://bioinformatics.ust.hk/Jlfdr.html.


  Measurement of the life time of attachments formed by a single microtubule
(MT) with a single kinetochore (kt) {\it in-vitro} under force-clamp conditions
had earlier revealed a catch-bond-like behavior. In the past the physical
origin of this apparently counter-intuitive phenomenon was traced to the nature
of the force-dependence of the (de-)polymerization kinetics of the MTs. Here
first the same model MT-kt attachment is subjected to external tension that
increases linearly with time until rupture occurs. In our {\it force-ramp}
experiments {\it in-silico}, the model displays the well known `mechanical
signatures' of a catch-bond probed by molecular force spectroscopy. Exploiting
this new evidence, we have further strengthened the analogy between MT-kt
attachments and common ligand-receptor bonds in spite of the crucial
differences in their underlying physical mechanisms. We then extend the
formalism to model the stochastic kinetics of an attachment formed by a bundle
of multiple parallel microtubules with a single kt considering the effect of
rebinding under force-clamp and force-ramp conditions. From numerical studies
of the model we predict the trends of variation of the mean life time and mean
rupture force with the increasing number of MTs in the bundle. Both the mean
life time and the mean rupture force display nontrivial nonlinear dependence on
the maximum number of MTs that can attach simultaneously to the same kt.


  Liquid-liquid phase separation of charge/aromatic-enriched intrinsically
disordered proteins (IDPs) is critical in the biological function of
membraneless organelles. Much of the physics of this recent discovery remains
to be elucidated. Here we present a theory in the random phase approximation to
account for electrostatic effects in polyampholyte phase separations, yielding
predictions consistent with recent experiments on the IDP Ddx4. The theory is
applicable to any charge pattern and thus provides a general analytical
framework for studying sequence dependence of IDP phase separation.


  To understand the behaviour of complex systems it is often necessary to use
models that describe the dynamics of subnetworks. It has previously been
established using projection methods that such subnetwork dynamics generically
involves memory of the past, and that the memory functions can be calculated
explicitly for biochemical reaction networks made up of unary and binary
reactions. However, many established network models involve also
Michaelis-Menten kinetics, to describe e.g. enzymatic reactions. We show that
the projection approach to subnetwork dynamics can be extended to such
networks, thus significantly broadening its range of applicability. To derive
the extension we construct a larger network that represents enzymes and enzyme
complexes explicitly, obtain the projected equations, and finally take the
limit of fast enzyme reactions that gives back Michaelis-Menten kinetics. The
crucial point is that this limit can be taken in closed form. The outcome is a
simple procedure that allows one to obtain a description of subnetwork
dynamics, including memory functions, starting directly from any given network
of unary, binary and Michaelis-Menten reactions. Numerical tests show that this
closed form enzyme elimination gives a much more accurate description of the
subnetwork dynamics than the simpler method that represents enzymes explicitly,
and is also more efficient computationally.


  Variant detection from high-throughput sequencing data is an essential step
in identification of alleles involved in complex diseases and cancer. To deal
with these massive data, elaborated sequence analysis pipelines are employed. A
core component of such pipelines is a read mapping module whose accuracy
strongly affects the quality of resulting variant calls.
  We propose a dynamic read mapping approach that significantly improves read
alignment accuracy. The general idea of dynamic mapping is to continuously
update the reference sequence on the basis of previously computed read
alignments. Even though this concept already appeared in the literature, we
believe that our work provides the first comprehensive analysis of this
approach.
  To evaluate the benefit of dynamic mapping, we developed a software pipeline
(http://github.com/karel-brinda/dymas) that mimics different dynamic mapping
scenarios. The pipeline was applied to compare dynamic mapping with the
conventional static mapping and, on the other hand, with the so-called
iterative referencing - a computationally expensive procedure computing an
optimal modification of the reference that maximizes the overall quality of all
alignments. We conclude that in all alternatives, dynamic mapping results in a
much better accuracy than static mapping, approaching the accuracy of iterative
referencing.
  To correct the reference sequence in the course of dynamic mapping, we
developed an online consensus caller named OCOCO
(http://github.com/karel-brinda/ococo). OCOCO is the first consensus caller
capable to process input reads in the online fashion.
  Finally, we provide conclusions about the feasibility of dynamic mapping and
discuss main obstacles that have to be overcome to implement it. We also review
a wide range of possible applications of dynamic mapping with a special
emphasis on variant detection.


  We investigate the replicator dynamics with "sparse" symmetric interactions
which represent specialist-specialist interactions in ecological communities.
By considering a large self interaction $u$, we conduct a perturbative
expansion which manifests that the nature of the interactions has a direct
impact on the species abundance distribution. The central results are all
species coexistence in a realistic range of the model parameters and that a
certain discrete nature of the interactions induces multiple peaks in the
species abundance distribution, providing the possibility of theoretically
explaining multiple peaks observed in various field studies. To get more
quantitative information, we also construct a non-perturbative theory which
becomes exact on tree-like networks if all the species coexist, providing exact
critical values of $u$ below which extinct species emerge. Numerical
simulations in various different situations are conducted and they clarify the
robustness of the presented mechanism of all species coexistence and multiple
peaks in the species abundance distributions.


  Intracellular transport is vital for the proper functioning and survival of a
cell. Cargo (proteins, vesicles, organelles, etc.) is transferred from its
place of creation to its target locations via molecular motor assisted
transport along cytoskeletal filaments. The transport efficiency is strongly
affected by the spatial organization of the cytoskeleton, which constitutes an
inhomogeneous, complex network. In cells with a centrosome microtubules grow
radially from the central microtubule organizing center towards the cell
periphery whereas actin filaments form a dense meshwork, the actin cortex,
underneath the cell membrane with a broad range of orientations. The emerging
ballistic motion along filaments is frequently interrupted due to constricting
intersection nodes or cycles of detachment and reattachment processes in the
crowded cytoplasm. In order to investigate the efficiency of search strategies
established by the cell's specific spatial organization of the cytoskeleton we
formulate a random velocity model with intermittent arrest states. With
extensive computer simulations we analyze the dependence of the mean first
passage times for narrow escape problems on the structural characteristics of
the cytoskeleton, the motor properties and the fraction of time spent in each
state. We find that an inhomogeneous architecture with a small width of the
actin cortex constitutes an efficient intracellular search strategy.


  The protein-protein interactions (PPIs) are crucial for understanding the
majority of cellular processes. PPIs play important role in gene transcription
regulation, cellular signaling, molecular basis of immune response and more.
Moreover, a disruption of hese mechanisms is frequently postulated as a
possible cause of diseases such as Alzheimer's or cancer. For many of
biologically relevant cases the structure of protein-protein complexes remain
unknown. Therefore computational techniques, including molecular docking, have
become a valuable part of drug discovery pipelines. Unfortunately, none of the
widely used protein-protein docking tools is free from serious limitations.
Typically, in docking simulations the protein flexibility is either completely
neglected or very limited. Additionally, some knowledge of the approximate
location and/or the shape of the active site is also required. Such limitations
arise mostly from the enormous number of degrees of freedom of protein-protein
systems. In this paper, an efficient computational method for protein-protein
docking is proposed and initially tested on a single docking case. The proposed
method is based on a two-step procedure. In the first step, CABS-dock web
server for protein-peptide docking is used to dock a peptide, which is the
appropriate protein fragment responsible for the protein-protein interaction,
to the other protein partner. During peptide docking, no knowledge about the
binding site, nor the peptide structure, is used and the peptide is allowed to
be fully flexible. In the second step, the docked peptide is used in the
structural adjustment of protein complex partners. The proposed method allowed
us to obtain a high accuracy model, therefore it provides a promising framework
for further advances.


  Despite considerable efforts, structural prediction of protein-peptide
complexes is still a very challenging task, mainly due to two reasons: high
flexibility of the peptides and transient character of their interactions with
proteins. Recently we have developed an automated web server CABS-dock
(http://biocomp.chem.uw.edu.pl/CABSdock), which conducts flexible
protein-peptide docking without any knowledge about the binding site. Our
method allows for full flexibility of the peptide, whereas the flexibility of
the receptor is restricted to near native conformations considering the main
chain, and full flexibility of the side chains. Performance of the CABS-dock
server was thoroughly tested on a benchmark of 171 test cases, both bound and
unbound. Evaluation of the obtained results showed overall good performance of
the method, especially that no information of the binding site was used. From
unsuccessful experiments we learned that the accuracy of docking might be
significantly improved, if only little information of the binding site was
considered. In fact, in real-life applications user typically has access to
some data indicating the location and/or structure of the binding site. In the
current work, we test and demonstrate the performance of the CABS-dock server
with two new features. The first one allows to utilize the knowledge about
receptor residue(s) constituting the binding site, and the second one allows to
enforce the desired secondary structure on the peptide structure. Based on the
given example, we observe significant improvement of the docking accuracy in
comparison to the default CABS-dock mode.


  Protein-peptide molecular docking is a difficult modeling problem. It is even
more challenging when significant conformational changes that may occur during
the binding process need to be predicted. In this chapter, we demonstrate the
capabilities and features of the CABS-dock server for flexible protein-peptide
docking. CABS-dock allows highly efficient modeling of full peptide flexibility
and significant flexibility of a protein receptor. During CABS-dock docking,
the peptide folding and binding process is explicitly simulated and no
information about the peptide binding site or its structure, is used. This
chapter presents a successful CABS-dock use for docking a potentially
therapeutic peptide to a protein target. Moreover, simulation contact maps, a
new CABS-dock feature, are described and applied to the docking test case.
Finally, a tutorial for running CABS-dock from the command line or command line
scripts is provided. The CABS-dock web server is available from
http://biocomp.chem.uw.edu.pl/CABSdock/


  We assess the impact of cell cycle noise on gene circuit dynamics. For
bistable genetic switches and excitable circuits, we find that transitions
between metastable states most likely occur just after cell division and that
this concentration effect intensifies in the presence of transcriptional delay.
We explain this concentration effect with a 3-states stochastic model. For
genetic oscillators, we quantify the temporal correlations between daughter
cells induced by cell division. Temporal correlations must be captured properly
in order to accurately quantify noise sources within gene networks.


  The speed of site-specific binding of transcription factor (TFs) proteins
with genomic DNA seems to be strongly retarded by the randomly occurring
sequence traps. Traps are those DNA sequences sharing significant similarity
with the original specific binding sites. It is an intriguing question how the
naturally occurring TFs and their specific binding sites are designed to manage
the retarding effects of such randomly occurring traps. We develop a simple
random walk model on the site-specific binding of TFs with genomic DNA in the
presence of sequence traps. Our dynamical model predicts that (a) the retarding
effects of traps will be minimum when the traps are arranged around the
specific binding site such that there is a negative correlation between the
binding strength of TFs with traps and the distance of traps from the specific
binding site and (b) the retarding effects of sequence traps can be appeased by
the condensed conformational state of DNA. Our computational analysis results
on the distribution of sequence traps around the putative binding sites of
various TFs in mouse and human genome clearly agree well the theoretical
predictions. We propose that the distribution of traps can be used as an
additional metric to efficiently identify the specific binding sites of TFs on
genomic DNA.


  Most disorders are caused by a combination of multiple genetic and/or
environmental factors. If two diseases are caused by the same molecular
mechanism, they tend to co-occur in patients. Here we provide a quantitative
method to disentangle how much genetic or environmental risk factors contribute
to the pathogenesis of 358 individual diseases, respectively. We pool data on
genetic, pathway-based, and toxicogenomic disease-causing mechanisms with
disease co-occurrence data obtained from almost two million patients. From this
data we construct a multilayer network where nodes represent disorders that are
connected by links that either represent phenotypic comorbidity of the patients
or the involvement of a certain molecular mechanism. From the similarity of
phenotypic and mechanism-based networks for each disorder we derive measure
that allows us to quantify the relative importance of various molecular
mechanisms for a given disease. We find that most diseases are dominated by
genetic risk factors, while environmental influences prevail for disorders such
as depressions, cancers, or dermatitis. Almost never we find that more than one
type of mechanisms is involved in the pathogenesis of diseases.


  In this paper we study the long term evolution of a continuous time Markov
chain formed by two interacting birth-and-death processes. The interaction
between the processes is modelled by transition rates which are functions with
suitable monotonicity properties. This is in line with the approach proposed by
Gauss G.F. and Kolmogorov A.N. for modelling interaction between species in
ecology. We obtain conditions for transience/recurrence of the Markov chain and
describe in detail its asymptotic behaviour in special transient cases. In
particular, we find that in some of these cases the Markov chain escapes to
infinity in an unusual way, and the corresponding trajectories can be rather
precisely described.


  Understanding the main determinants of species coexistence across space and
time is a central question in ecology. However, ecologists still know little
about the scales and conditions at which biotic interactions matter and how
these interact with the environment to structure species assemblages. Here we
use recent theory developments to analyze plant distribution and trait data
across Europe and find that plant height clustering is related to both
evapotranspiration and gross primary productivity. This clustering is a signal
of interspecies competition between plants, which is most evident in
mid-latitude ecoregions, where conditions for growth (reflected in actual
evapotranspiration rates and gross primary productivities) are optimal. Away
from this optimum, climate severity likely overrides the effect of competition,
or other interactions become increasingly important. Our approach bridges the
gap between species-rich competition theories and large-scale species
distribution data analysis.


  The article analyzes the contribution of stochastic thermal fluctuations in
the attachment times of the immature T-cell receptor TCR:
peptide-major-histocompatibility-complex pMHC immunological synapse bond. The
key question addressed here is the following: how does a synapse bond remain
stabilized in the presence of high frequency thermal noise that potentially
equates to a strong detaching force? Focusing on the average time persistence
of an immature synapse, we show that the high frequency nodes accompanying
large fluctuations are counterbalanced by low frequency nodes that evolve over
longer time periods. Our analysis shows that such a behavior could be easily
explained from the fact that the survival probability distribution is governed
by two distinct phases, for the two different time regimes. The relatively
shorter time scales correspond to the cohesion:adhesion induced immature bond
formation whereas the larger time reciprocates the association:dissociation
regime leading to TCR:pMHC signaling. From an estimation of the bond survival
probability, we show that at shorter time scales, this probability
$P_{\Delta}(\tau)$ scales with time $\tau$ as an universal function of a
rescaled noise amplitude $\frac{D}{\Delta^2}$, such that $P_{\Delta}(\tau)\sim
\tau^{-(\frac{\Delta}{\sqrt{D}}+\frac{1}{2})}$, $\Delta$ being the distance
from the mean inter-membrane (T cell:Antigen Presenting Cell) separation
distance. The crossover from this shorter to a longer time regime leads to an
universality in the dynamics, at which point the survival probability shows a
different power-law scaling compared to the one at shorter time scales. In
biological terms, such a crossover indicates that the TCR:pMHC bond has a
survival probability with a slower decay rate than the longer LFA-1:ICAM-1 bond
justifying its stability.


  Data from a long time evolution experiment with Escherichia Coli and from a
large study on copy number variations in subjects with european ancestry are
analyzed in order to argue that mutations can be described as Levy flights in
the mutation space. These Levy flights have at least two components: random
single-base substitutions and large DNA rearrangements. From the data, we get
estimations for the time rates of both events and the size distribution
function of large rearrangements.


  Allosteric regulation is central to many biochemical processes. Allosteric
sites provide a target to fine-tune protein activity, yet we lack computational
methods to predict them. Here, we present an efficient graph-theoretical
approach for identifying allosteric sites and the mediating interactions that
connect them to the active site. Using an atomistic graph with edges weighted
by covalent and non-covalent bond energies, we obtain a bond-to-bond propensity
that quantifies the effect of instantaneous bond fluctuations propagating
through the protein. We use this propensity to detect the sites and
communication pathways most strongly linked to the active site, assessing their
significance through quantile regression and comparison against a reference set
of 100 generic proteins. We exemplify our method in detail with three
well-studied allosteric proteins: caspase-1, CheY, and h-Ras, correctly
predicting the location of the allosteric site and identifying key allosteric
interactions. Consistent prediction of allosteric sites is then attained in a
further set of 17 proteins known to exhibit allostery. Because our propensity
measure runs in almost linear time, it offers a scalable approach to
high-throughput searches for candidate allosteric sites.


  A genome-scale reconstruction of human metabolism, Recon 2, is available but
no interface exists to interactively visualise its content integrated with
omics data and simulation results. We manually drew a comprehensive map,
ReconMap 2.0, that is consistent with the content of Recon 2. We present it
within a web interface that allows content query, visualization of custom
datasets and submission of feedback to manual curators. ReconMap can be
accessed via http://vmh.uni.lu, with network export in a Systems Biology
Graphical Notation compliant format. A Constraint-Based Reconstruction and
Analysis (COBRA) Toolbox extension to interact with ReconMap is available via
https://github.com/opencobra/cobratoolbox.


  Constraint-Based Reconstruction and Analysis (COBRA) is currently the only
methodology that permits integrated modeling of Metabolism and macromolecular
Expression (ME) at genome-scale. Linear optimization computes steady-state flux
solutions to ME models, but flux values are spread over many orders of
magnitude. Standard double-precision solvers may return inaccurate solutions or
report that no solution exists. Exact simplex solvers are extremely slow and
hence not practical for ME models that currently have 70,000 constraints and
variables and will grow larger. We have developed a quadruple-precision version
of our linear and nonlinear optimizer MINOS, and a solution procedure (DQQ)
involving Double and Quad MINOS that achieves efficiency and reliability for ME
models. DQQ enables extensive use of large, multiscale, linear and nonlinear
models in systems biology and many other applications.


  The capacity of cells and organisms to respond to challenging conditions in a
repeatable manner is limited by a finite repertoire of pre-evolved adaptive
responses. Beyond this capacity, cells can use exploratory dynamics to cope
with a much broader array of conditions. However, the process of adaptation by
exploratory dynamics within the lifetime of a cell is not well understood. Here
we demonstrate the feasibility of exploratory adaptation in a high-dimensional
network model of gene regulation. Exploration is initiated by failure to comply
with a constraint and is implemented by random sampling of network
configurations. It ceases if and when the network reaches a stable state
satisfying the constraint. We find that successful convergence (adaptation) in
high dimensions requires outgoing network hubs and is enhanced by their
auto-regulation. The ability of these empirically-validated features of gene
regulatory networks to support exploratory adaptation without fine-tuning,
makes it plausible for biological implementation.


  MicroRNAs play critical roles in many physiological processes. Their
dysregulations are also closely related to the development and progression of
various human diseases, including cancer. Therefore, identifying new microRNAs
that are associated with diseases contributes to a better understanding of
pathogenicity mechanisms. MicroRNAs also represent a tremendous opportunity in
biotechnology for early diagnosis. To date, several in silico methods have been
developed to address the issue of microRNA-disease association prediction.
However, these methods have various limitations. In this study, we investigate
the hypothesis that information attached to miRNAs and diseases can be revealed
by distributional semantics. Our basic approach is to represent distributional
information on miRNAs and diseases in a high-dimensional vector space and to
define associations between miRNAs and diseases in terms of their vector
similarity. Cross validations performed on a dataset of known miRNA-disease
associations demonstrate the excellent performance of our method. Moreover, the
case study focused on breast cancer confirms the ability of our method to
discover new disease-miRNA associations and to identify putative false
associations reported in databases.


  The cellular adaptive immune response plays a key role in resolving influenza
infection. Experiments where individuals are successively infected with
different strains within a short timeframe provide insight into the underlying
viral dynamics and the role of a cross-reactive immune response in resolving an
acute infection. We construct a mathematical model of within-host influenza
viral dynamics including three possible factors which determine the strength of
the cross-reactive cellular adaptive immune response: the initial naive T cell
number, the avidity of the interaction between T cells and the epitopes
presented by infected cells, and the epitope abundance per infected cell. Our
model explains the experimentally observed shortening of a second infection
when cross-reactivity is present, and shows that memory in the cellular
adaptive immune response is necessary to protect against a second infection.


  It is very common to observe nonlinear features in agricultural and
ecological systems. For example, in tree crop production, alternate bearing is
well known phenomena caused by nonlinear dynamics. Production of single tree of
Citrus Unshiu trees is recognized to be driven by a mechanistic process modeled
with the so-called "resource budget model", which demonstrates phenomenon of
alternate bearing. However, the term of alternative bearing is used not only
for an individual tree's production but also for total production of a large
sized population of trees. In this paper, we developed noise induced uncoupled
dynamics model for population alternate bearing based on Isagi's resource
budget model. Based on numerical experiments with the developed model,
theoretical possibility of a substantial alternate bearing effect even in a
national market was proposed.


  Scaling theory generates transferable (even universal) algebraic and
geometrical relations between the amino acid sequences and the aggregation
functions of the three titled radically different proteins. In addition to the
two hydropathicity scales and beta strand scales used in earlier p53 work, a
second beta strand Hot Spot scale is shown to yield very accurate results for
oligomerization of p53, the tumor suppressor. These algebraic and geometrical
relations could be caused topologically by the dominance of protein-protein
aggregation by interactions in a membrane catalytic surface layer.