  Allostery is an intrinsic spatiotemporal property of all proteins, resulting
from long range correlations in the order of several nanometers and time scales
of nanoseconds. Information is carried asymmetrically from one part to another
by entropy transfer. Here, we present a master equation model of allosteric
communication in proteins based on the transfer entropy concept of Schreiber
(PRL, 85, 465, 2000). We show how the model relates the path and velocity of
asymmetric entropy transfer to conformational transitions over the rugged
energy surface of proteins and how this relates to function.


  We show a proof of principle for warping, a method to interpret the inner
working of neural networks in the context of gene expression analysis. Warping
is an efficient way to gain insight to the inner workings of neural nets and
make them more interpretable. We demonstrate the ability of warping to recover
meaningful information for a given class on a samplespecific individual basis.
We found warping works well in both linearly and nonlinearly separable
datasets. These encouraging results show that warping has a potential to be the
answer to neural networks interpretability in computational biology.


  CRISPR is a newly discovered prokaryotic immune system. Bacteria and archaea
with this system incorporate genetic material from invading viruses into their
genomes, providing protection against future infection by similar viruses. The
conditions for coexistence of prokaryots and viruses is an interesting problem
in evolutionary biology. In this work, we show an intriguing phase diagram of
the virus extinction probability, which is more complex than that of the
classical predator-prey model. As the CRISPR incorporates genetic material,
viruses are under pressure to evolve to escape the recognition by CRISPR. When
bacteria have a small rate of deleting spacers, a new parameter region in which
bacteria and viruses can coexist arises, and it leads to a more complex
coexistence patten for bacteria and viruses. For example, when the virus
mutation rate is low, the virus extinction probability changes non-montonically
with the bacterial exposure rate. The virus and bacteria co-evolution not only
alters the virus extinction probability, but also changes the bacterial
population structure. Additionally, we show that recombination is a successful
strategy for viruses to escape from CRISPR recognition when viruses have
multiple proto-spacers, providing support for a recombination-mediated escape
mechanism suggested experimentally. Finally, we suggest that the reentrant
phase diagram, in which phages can progress through three phases of extinction
and two phases of abundance at low spacer deletion rates as a function of
exposure rate to bacteria, is an experimentally testable phenomenon.


  The nonlocal Fisher equation is a diffusion-reaction equation with a nonlocal
quadratic competition, which describes the reaction between distant
individuals. This equation arises in evolutionary biological systems, where the
arena for the dynamics is trait space, diffusion accounts for mutations and
individuals with similar traits compete, resulting in partial niche overlap. It
has been found that the (non-cutoff) deterministic system gives rise to a
spatially inhomogeneous state for a certain class of interaction kernels, while
the stochastic system produces an inhomogeneous state for small enough
population densities. Here we study the problem of front propagation in this
system, comparing the stochastic dynamics to the heuristic approximation of
this system by a deterministic system where the linear growth term is cut off
below some critical density. Of particular interest is the nontrivial pattern
left behind the front. For large population density, or small cutoff, there is
a constant velocity wave propagating from the populated region to the
unpopulated region. As in the local Fisher equation, the spreading velocity is
much lower than the Fisher velocity which is the spreading velocity for
infinite population size. The stochastic simulations give approximately the
same spreading velocity as the deterministic simulation with appropriate birth
cutoff. When the population density is small enough, there is a different
mechanism of population spreading. The population is concentrated on clusters
which divide and separate. This mode of spreading has small spreading velocity,
decaying exponentially with the range of the interaction kernel.


  The mathematical modelling of genetic-metabolic networks is of out most
importance in the field of systems biology. Different formalisms and a huge
variety of classical mathematical tools have been used to describe and analyse
such networks. Michael A. Savageau defined a formalism to model
genetic-metabolic networks called S-Systems, [12], [13], [15]. There is a limit
in the number of nodes that can be analysed when these systems are solved using
classical numerical methods such as non-linear dynamic analysis and linear
optimization algorithms. We propose to use toric algebraic geometry to solve
S-systems. In this work we prove that S-systems are toric varieties and that as
a consequence Hilbert basis can be used to solve them. This is achieved by
applying two theorems, proved here, the theorem about Embedding of S-Systems in
toric varieties and the theorem about Toric Resolution on S-Systems. In
addition, we define the realization of minimal phenotypic polytopes, phenotypic
toric ideals and phenotypic toric varieties as a generalization of the
phenotypic polytopes described by M. Savageau, [8], [9], [10]. Also we studied
a synthetic oscillatory network of two genes under diverse environmental
control modes, we reproduce S. Savageau results and additionally we obtain an
invariant sustained oscillatory region. The implications of the results here
presented is that, in principle, they will facilitate solving large scale
genetic-metabolic networks by means of toric algebraic geometry tools as well
as facilitate elucidating their dynamics.


  Much of our understanding of ecological and evolutionary mechanisms derives
from analysis of low-dimensional models: with few interacting species, or few
axes defining "fitness". It is not always clear to what extent the intuition
derived from low-dimensional models applies to the complex, high-dimensional
reality. For instance, most naturally occurring microbial communities are
strikingly diverse, harboring a large number of coexisting species, each of
which contributes to shaping the environment of others. Understanding the
eco-evolutionary interplay in these systems is an important challenge, and an
exciting new domain for statistical physics. Recent work identified a promising
new platform for investigating highly diverse ecosystems, based on the classic
resource competition model of MacArthur. Here, we describe how the same
analytical framework can be used to study evolutionary questions. Our analysis
illustrates how, at high dimension, the intuition promoted by a one-dimensional
(scalar) notion of fitness can become misleading. Specifically, while the
low-dimensional picture emphasizes organism cost or efficiency, we exhibit a
regime where cost becomes irrelevant for survival, and link this observation to
generic properties of high-dimensional geometry.


  The prevalent method for RNA secondary structure prediction for a single
sequence is free energy minimization based on the nearest neighbor
thermodynamic model (NNTM). One of the least well-developed parts of the model
is the energy function assigned to the multibranch loops. Parametric analysis
can be performed to elucidate the dependance of the prediction on the branching
parameters used in the NNTM. Since the objective function is linear, this boils
down to analyzing the normal fans of the branching polytopes. Here we show that
because of the way the multibranch loops are scored under the NNTM, certain
branching patterns are possible for all sequences. We do this by characterizing
the dominant parts of the parameter space obtained by looking at the relevant
section of the normal fan; therefore, we conclude that the structures that are
normally found in nature are obtained for a relatively small set of parameters.


  MicroRNAs (miRNAs) are small endogenous regulatory molecules that modulate
gene expression post-transcriptionally. Although differential expression of
miRNAs have been implicated in many diseases (including cancers), the
underlying mechanisms of action remain unclear. Because each miRNA can target
multiple genes, miRNAs may potentially have functional implications for the
overall behavior of entire pathways. Here we investigate the functional
consequences of miRNA dysregulation through an integrative analysis of miRNA
and mRNA expression data using a novel approach that incorporates pathway
information a priori. By searching for miRNA-pathway associations that differ
between healthy and tumor tissue, we identify specific relationships at the
systems-level which are disrupted in cancer. Our approach is motivated by the
hypothesis that if a miRNA and pathway are associated, then the expression of
the miRNA and the collective behavior of the genes in a pathway will be
correlated. As such, we first obtain an expression-based summary of pathway
activity using Isomap, a dimension reduction method which can articulate
nonlinear structure in high-dimensional data. We then search for miRNAs that
exhibit differential correlations with the pathway summary between phenotypes
as a means of finding aberrant miRNA-pathway coregulation in tumors. We apply
our method to cancer data using gene and miRNA expression datasets from The
Cancer Genome Atlas (TCGA) and compare ${\sim}10^5$ miRNA-pathway relationships
between healthy and tumor samples from four tissues (breast, prostate, lung,
and liver). Many of the flagged pairs we identify have a biological basis for
disruption in cancer.


  Single-cell gene expression is inherently stochastic; its emergent behavior
can be defined in terms of the chemical master equation describing the
evolution of the mRNA and protein copy numbers as the latter tends to infinity.
We establish two types of "macroscopic limits": the Kurtz limit is consistent
with the classical chemical kinetics, while the L\'{e}vy limit provides a
theoretical foundation for an empirical equation proposed in [Phys. Rev. Lett.
97:168302, 2006]. Furthermore, we clarify the biochemical implications and
ranges of applicability for various macroscopic limits and calculate a
comprehensive analytic expression for the protein concentration distribution in
autoregulatory gene networks. The relationship between our work and modern
population genetics is discussed.


  Bacteria populations rely on mechanisms such as quorum sensing to coordinate
complex tasks that cannot be achieved by a single bacterium. Quorum sensing is
used to measure the local bacteria population density, and it controls
cooperation by ensuring that a bacterium only commits the resources for
cooperation when it expects its neighbors to reciprocate. This paper proposes a
simple model for sharing a resource in a bacterial environment, where knowledge
of the population influences each bacterium's behavior. Game theory is used to
model the behavioral dynamics, where the net payoff (i.e., utility) for each
bacterium is a function of its current behavior and that of the other bacteria.
The game is first evaluated with perfect knowledge of the population. Then, the
unreliability of diffusion introduces uncertainty in the local population
estimate and changes the perceived payoffs. The results demonstrate the
sensitivity to the system parameters and how population uncertainty can
overcome a lack of explicit coordination.


  We propose a model for the human immunodeficiency virus type 1 (HIV-1)
infection with intracellular delay and prove the local asymptotical stability
of the equilibrium points. Then we introduce a control function representing
the efficiency of reverse transcriptase inhibitors and consider the
pharmacological delay associated to the control. Finally, we propose and
analyze an optimal control problem with state and control delays. Through
numerical simulations, extremal solutions are proposed for minimization of the
virus concentration and treatment costs.


  We investigate three aspects of the importance of nodes with respect to
Susceptible-Infectious-Removed (SIR) disease dynamics: influence maximization
(the expected outbreak size given a set of seed nodes), the effect of
vaccination (how much deleting nodes would reduce the expected outbreak size)
and sentinel surveillance (how early an outbreak could be detected with sensors
at a set of nodes). We calculate the exact expressions of these quantities, as
functions of the SIR parameters, for all connected graphs of three to seven
nodes. We obtain the smallest graphs where the optimal node sets are not
overlapping. We find that: node separation is more important than centrality
for more than one active node, that vaccination and influence maximization are
the most different aspects of importance, and that the three aspects are more
similar when the infection rate is low.


  Kauffman's NK-model is a paradigmatic example of a class of stochastic models
of genotypic fitness landscapes that aim to capture generic features of
epistatic interactions in multilocus systems. Genotypes are represented as
sequences of $L$ binary loci. The fitness assigned to a genotype is a sum of
contributions, each of which is a random function defined on a subset of $k \le
L$ loci. These subsets or neighborhoods determine the genetic interactions of
the model. Whereas earlier work on the NK model suggested that most of its
properties are robust with regard to the choice of neighborhoods, recent work
has revealed an important and sometimes counter-intuitive influence of the
interaction structure on the properties of NK fitness landscapes. Here we
review these developments and present new results concerning the number of
local fitness maxima and the statistics of selectively accessible (that is,
fitness-monotonic) mutational pathways. In particular, we develop a unified
framework for computing the exponential growth rate of the expected number of
local fitness maxima as a function of $L$, and identify two different
universality classes of interaction structures that display different
asymptotics of this quantity for large $k$. Moreover, we show that the
probability that the fitness landscape can be traversed along an accessible
path decreases exponentially in $L$ for a large class of interaction structures
that we characterize as locally bounded. Finally, we discuss the impact of the
NK interaction structures on the dynamics of evolution using adaptive walk
models.


  The dynamics of gene regulatory networks are often modeled with the
assumption of cellular homogeneity. However, this assumption contradicts the
plethora of experimental results in a variety of systems, which designates that
cell populations are heterogeneous systems in the sense that properties such as
size, shape, and DNA/RNA content are unevenly distributed amongst their
individuals. In order to address the implications of heterogeneity, we utilize
the so-called cell population balance (CPB) models. Here, we solve numerically
multivariable CPB models to study the effect of heterogeneity on populations
carrying the toggle switch network, which features nonlinear behavior at the
single-cell level. In order to answer whether this nonlinear behavior is
inherited to the heterogeneous population level, we perform bifurcation
analysis on the steady-state solutions of the CPB model. We show that
bistability is present at the population level with the pertinent bistability
region shrinking when the impact of heterogeneity is enhanced.


  When the interactions between cooperators (C) and defectors (D) can be
partially avoided within a population, there may be an overall enhancement of
cooperation. One example of such screening mechanism occurs in the presence of
risk-averse agents (loners, L) that are neutral towards others, i.e., both L
and its opponent, whatever its strategy, receive the same payoff. Their
presence in the Prisoner's Dilemma (PD) game sustains the coexistence of
cooperators and defectors far beyond the level attained in their absence.
Another screening mechanism is a heterogeneous landscape obtained, for example,
by site diluting the lattice. In this case, cooperation is enhanced with some
fraction of such inactive, interaction-averse sites. By considering the
interplay of both mechanisms, we show that there is an explosive increase in
the range of densities, just above the percolation threshold, where neutrality
is prevented and loners become extinct, the behavior reverting to the pure PD
game. Interestingly, this occurs despite defectors being usually abundant in
that region. This has to be compared with the corresponding loner-free region
in the undiluted case that, besides being very small, is dominatedby
cooperators.


  The formalism of partial information decomposition provides independent or
non-overlapping components constituting total information content provided by a
set of source variables about the target variable. These components are
recognised as unique information, synergistic information and, redundant
information. The metric of net synergy, conceived as the difference between
synergistic and redundant information, is capable of detecting synergy,
redundancy and, information independence among stochastic variables. And it can
be quantified, as it is done here, using appropriate combinations of different
Shannon mutual information terms. Utilisation of such a metric in network
motifs with the nodes representing different biochemical species, involved in
information sharing, uncovers rich store for interesting results. In the
current study, we make use of this formalism to obtain a comprehensive
understanding of the relative information processing mechanism in a diamond
motif and two of its sub-motifs namely bifurcation and integration motif
embedded within the diamond motif. The emerging patterns of synergy and
redundancy and their effective contribution towards ensuring high fidelity
information transmission are duly compared in the sub-motifs and independent
motifs (bifurcation and integration). In this context, the crucial roles played
by various time scales and activation coefficients in the network topologies
are especially emphasised. We show that the origin of synergy and redundancy in
information transmission can be physically justified by decomposing diamond
motif into bifurcation and integration motif.


  In this paper, a multiscale virtual particle based elastic network model
(MVP-ENM) is proposed for biomolecular normal mode analysis. The multiscale
virtual particle model is proposed for the discretization of biomolecular
density data in different scales. Essentially, the model works as the
coarse-graining of the biomolecular structure, so that a delicate balance
between biomolecular geometric representation and computational cost can be
achieved. To form "connections" between these multiscale virtual particles, a
new harmonic potential function, which considers the influence from both mass
distributions and distance relations, is adopted between any two virtual
particles. Unlike the previous ENMs that use a constant spring constant, a
particle-dependent spring parameter is used in MVP-ENM. Two independent models,
i.e., multiscale virtual particle based Gaussian network model (MVP-GNM) and
multiscale virtual particle based anisotropic network model (MVP-ANM), are
proposed. Even with a rather coarse grid and a low resolution, the MVP-GNM is
able to predict the Debye-Waller factors (B-factors) with considerable good
accuracy. Similar properties have also been observed in MVP-ANM. More
importantly, in B-factor predictions, the mismatch between the predicted
results and experimental ones is predominantly from higher fluctuation regions.
Further, it is found that MVP-ANM can deliver a very consistent low-frequency
eigenmodes in various scales. This demonstrates the great potential of MVP-ANM
in the deformation analysis of low resolution data. With the multiscale
rigidity function, the MVP-ENM can be applied to biomolecular data represented
in density distribution and atomic coordinates. Further, the great advantage of
my MVP-ENM model in computational cost has been demonstrated by using two
poliovirus virus structures. Finally, the paper ends with a conclusion.


  According to implicit ligand theory, the standard binding free energy is an
exponential average of the binding potential of mean force (BPMF), an
exponential average of the interaction energy between the ligand apo ensemble
and a rigid receptor. Here, we use the Fast Fourier Transform (FFT) to
efficiently estimate BPMFs by calculating interaction energies as rigid ligand
configurations from the apo ensemble are discretely translated across rigid
receptor conformations. Results for standard binding free energies between T4
lysozyme and 141 small organic molecules are in good agreement with previous
alchemical calculations based on (1) a flexible complex (R ~ 0.9 for 24
systems) and (2) flexible ligand with multiple rigid receptor configurations (R
~ 0.8 for 141 systems). While the FFT is routinely used for molecular docking,
to our knowledge this is the first time that the algorithm has been used for
rigorous binding free energy calculations.


  Spatially extended population dynamics models that incorporate intrinsic
noise serve as case studies for the role of fluctuations and correlations in
biological systems. Including spatial structure and stochastic noise in
predator-prey competition invalidates the deterministic Lotka-Volterra picture
of neutral population cycles. Stochastic models yield long-lived erratic
population oscillations stemming from a resonant amplification mechanism. In
spatially extended predator-prey systems, one observes noise-stabilized
activity and persistent correlations. Fluctuation-induced renormalizations of
the oscillation parameters can be analyzed perturbatively. The critical
dynamics and the non-equilibrium relaxation kinetics at the predator extinction
threshold are characterized by the directed percolation universality class.
Spatial or environmental variability results in more localized patches which
enhances both species densities. Affixing variable rates to individual
particles and allowing for trait inheritance subject to mutations induces fast
evolutionary dynamics for the rate distributions. Stochastic spatial variants
of cyclic competition with rock-paper-scissors interactions illustrate
connections between population dynamics and evolutionary game theory, and
demonstrate how space can help maintain diversity. In two dimensions,
three-species cyclic competition models of the May-Leonard type are
characterized by the emergence of spiral patterns whose properties are
elucidated by a mapping onto a complex Ginzburg-Landau equation. Extensions to
general food networks can be classified on the mean-field level, which provides
both a fundamental understanding of ensuing cooperativity and emergence of
alliances. Novel space-time patterns emerge as a result of the formation of
competing alliances, such as coarsening domains that each incorporate
rock-paper-scissors competition games.


  Clustering cancer patients into subgroups and identifying cancer subtypes is
an important task in cancer genomics. Clustering based on comprehensive
multi-omic molecular profiling can often achieve better results than those
using a single data type, since each omic data type (representing one view of
patients) may contain complementary information. However, it is challenging to
integrate heterogeneous omic data types directly. Based on one popular method
-- Similarity Network Fusion (SNF), we presented Affinity Network Fusion (ANF)
in this paper, an "upgrade" of SNF with several advantages. Similar to SNF, ANF
treats each omic data type as one view of patients and learns a fused affinity
(transition) matrix for clustering. We applied ANF to a carefully processed
harmonized cancer dataset downloaded from GDC data portals consisting of 2193
patients, and generated promising results on clustering patients into correct
disease types. Our experimental results also demonstrated the power of feature
selection and transformation combined with using ANF in patient clustering.
Moreover, eigengap analysis suggests that the learned affinity matrices of four
cancer types using our proposed framework may have successfully captured
patient group structure and can be used for discovering unknown cancer
subtypes.


  Simple method to improve traditional approaches to name taxa of higher ranks
has been proposed. Instead of base name + postfixes (which vary between codes
of nomenclature and sometimes not standardized at all), use numerical prefixes,
where numbers designate the rank of taxon, from lowest (species = 1) to highest
(kingdom = 7). This "numerical ranks" method is not only more simple and
straightforward then current practice, but also easy to extend for numerous
additional uses. There is a hope that numerical ranks will facilitate creation
of typified names for all higher-level taxomonic groups.


  An efficient structural identifiability analysis algorithm is developed in
this study for a broad range of network structures. The proposed method adopts
the Wright's path coefficient method to generate identifiability equations in
forms of symbolic polynomials, and then converts these symbolic equations to
binary matrices (called identifiability matrix). Several matrix operations are
introduced for identifiability matrix reduction with system equivalency
maintained. Based on the reduced identifiability matrices, the structural
identifiability of each parameter is determined. A number of benchmark models
are used to verify the validity of the proposed approach. Finally, the network
module for influenza A virus replication is employed as a real example to
illustrate the application of the proposed approach in practice. The proposed
approach can deal with cyclic networks with latent variables. The key advantage
is that it intentionally avoids symbolic computation and is thus highly
efficient. Also, this method is capable of determining the identifiability of
each single parameter and is thus of higher resolution in comparison with many
existing approaches. Overall, this study provides a basis for systematic
examination and refinement of graphical models of biological networks from the
identifiability point of view, and it has a significant potential to be
extended to more complex network structures or high-dimensional systems.


  Motivation: In transcriptomic and proteomics-based studies, the abundance of
genes is often compared to functional pathways such as the Kyoto Encyclopaedia
at Genes and Genomes (KEGG) to identify active metabolic processes. Even though
a plethora of tools allow to analyze and to compare omics data in respect to
KEGG pathways, the analysis of multiple conditions is often limited to only a
defined set of conditions. Furthermore, for transcriptomic datasets, it is
crucial to compare the entire set of pathways in order to obtain a global
overview of the species' metabolic functions. Results: Here, we present the
tool KEGGexpressionMapper, a module, that is implemented in the programming
language R. The module allows to highlight the expression of transcriptomic or
proteomic measurements in various conditions on pathways and incorporates
methods to analyze gene enrichment analyses and expression clustering in time
series data. KEGGexpressionMapper supports time series data from transcriptomic
or proteomics measurements from different individuals. As the tool is
implemented in the scripting language R, it can be integrated into existing
analysis pipelines to obtain a global overview of the dataset. The R package
can be downloaded from https://github.com/nthomasCUBE/KEGGexpressionmapper.


  In the protein sequence space, natural proteins form clusters of families
which are characterized by their unique native folds whereas the great majority
of random polypeptides are neither clustered nor foldable to unique structures.
Since a given polypeptide can be either foldable or unfoldable, a kind of
"folding transition" is expected at the boundary of a protein family in the
sequence space. By Monte Carlo simulations of a statistical mechanical model of
protein sequence alignment that coherently incorporates both short-range and
long-range interactions as well as variable-length insertions to reproduce the
statistics of the multiple sequence alignment of a given protein family, we
demonstrate the existence of such transition between natural-like sequences and
random sequences in the sequence subspaces for 15 domain families of various
folds. The transition was found to be highly cooperative and two-state-like.
Furthermore, enforcing or suppressing consensus residues on a few of the
well-conserved sites enhanced or diminished, respectively, the natural-like
pattern formation over the entire sequence. In most families, the key sites
included ligand binding sites. These results suggest some selective pressure on
the key residues, such as ligand binding activity, may cooperatively facilitate
the emergence of a protein family during evolution. From a more practical
aspect, the present results highlight an essential role of long-range effects
in precisely defining protein families, which are absent in conventional
sequence models.


  One of the greatest challenges in biophysical models of translation is to
identify coding sequences features that affect the rate of translation and
therefore the overall protein production in the cell. We propose an analytic
method to solve a translation model based on the inhomogeneous totally
asymmetric simple exclusion process, which allows us to unveil simple design
principles of nucleotide sequences determining protein production rates. Our
solution shows an excellent agreement when compared to numerical genome-wide
simulations of S. cerevisiae transcript sequences and predicts that the first
10 codons, together with the value of the initiation rate, are the main
determinants of protein production rate. Finally, we interpret the obtained
analytic results based on the evolutionary role of codons' choice for
regulating translation rates and ribosome densities.


  We have extended our analytically derived PDB-NMA formulation, ATMAN [1], to
include protein dimers using mixed internal and Cartesian coordinates. A test
case on a 1.3\AA\ resolution model of a small homodimer, ActVA-ORF6, consisting
of two 112-residue subunits identically folded in a compact 50\AA\ sphere,
reproduces the distinct experimental Debye-Waller motility asymmetry for the
two chains, demonstrating that structure sensitively selects vibrational
signatures. The vibrational analysis of this PDB entry, together with
biochemical and crystallographic data, demonstrates the cooperative nature of
the dimeric interaction of the two subunits and suggests a mechanical model for
subunit interconversion during the catalytic cycle.


  Leucine zippers are alpha helical monomers dimerized to a coiled coil
structure. Various scenarios for dimerization span interaction of unstructured
monomers that form alpha helices in the process of dimerization to preformed
alpha helical monomers dimerizing.Another suggested possibility, a trigger
sequence hypothesis (M. O. Steinmetz et al., Proc. Natl. Acad. Sci. U. S. A.
2007, 104, 7062-7067), is that a C terminal (CT) part of each monomer is a
trigger sequence that is dimerized and, subsequently, the remainder N terminal
(NT) part of each monomer (that is partially ordered) zips together to form the
coiled coil. In this work, methods are developed to computationally explore the
trigger sequence hypothesis based on an extension of a previously introduced
(R. I. Cukier, J. Chem. Phys. 2011, 134, 045104) Hamiltonian Temperature
Replica Exchange Method Mean Square (HTREM_MS) procedure, which scales the
Hamiltonian in both potential and kinetic energies, to enhance sampling
generically, along with additional restraints to enhance sampling in desired
regions of conformational space. The method is applied to a 31-residue
truncation of the 33-residue leucine zipper (GCN4-p1) of the yeast
transcriptional activator GCN4. Using a variety of HTREM_MS simulations, we
find that the NT of one monomer of the dimer becomes disordered when the NT
parts of the monomers are separated. In contrast, when the CT parts of the
monomers are separated, both monomers remain alpha helical. These simulations
suggest that the CT does act as a trigger sequence and are consistent with a
disordered NT. We also investigate whether the disordered NT can be induced to
re-form the dimerized leucine zipper. It does, but with some lack of recovery
of alpha helical hydrogen bonding structure in the NT.


  The variational principle for conformational dynamics has enabled the
systematic construction of Markov state models through the optimization of
hyperparameters by approximating the transfer operator. In this note we discuss
why lag time of the operator being approximated must be held constant in the
variational approach.


  This work introduces a number of algebraic topology approaches, such as
multicomponent persistent homology, multi-level persistent homology and
electrostatic persistence for the representation, characterization, and
description of small molecules and biomolecular complexes. Multicomponent
persistent homology retains critical chemical and biological information during
the topological simplification of biomolecular geometric complexity.
Multi-level persistent homology enables a tailored topological description of
inter- and/or intra-molecular interactions of interest. Electrostatic
persistence incorporates partial charge information into topological
invariants. These topological methods are paired with Wasserstein distance to
characterize similarities between molecules and are further integrated with a
variety of machine learning algorithms, including k-nearest neighbors, ensemble
of trees, and deep convolutional neural networks, to manifest their descriptive
and predictive powers for chemical and biological problems. Extensive numerical
experiments involving more than 4,000 protein-ligand complexes from the PDBBind
database and near 100,000 ligands and decoys in the DUD database are performed
to test respectively the scoring power and the virtual screening power of the
proposed topological approaches. It is demonstrated that the present approaches
outperform the modern machine learning based methods in protein-ligand binding
affinity predictions and ligand-decoy discrimination.


  This review discusses research developments and applications of
isotachophoresis (ITP) to the initiation, control, and acceleration of chemical
reactions, emphasizing reactions involving biomolecular reactants such as
nucleic acids, proteins, and live cells. ITP is a versatile technique which
requires no specific geometric design or material, and is compatible with a
wide range of microfluidic and automated platforms. Though ITP has
traditionally been used as a purification and separation technique, recent
years have seen its emergence as a method to automate and speed up chemical
reactions. ITP has been used to demonstrate up to 14,000-fold acceleration of
nucleic acid assays, and has been used to enhance lateral flow and other
immunoassays, and even whole bacterial cell detection assays. We here classify
these studies into two categories: homogeneous (all reactants in solution) and
heterogeneous (at least one reactant immobilized on a solid surface) assay
configurations. For each category, we review and describe physical modeling and
scaling of ITP-aided reaction assays, and elucidate key principles in ITP assay
design. We summarize experimental advances, and identify common threads and
approaches which researchers have used to optimize assay performance. Lastly,
we propose unaddressed challenges and opportunities that could further improve
these applications of ITP.


  Computational elucidation of membrane protein (MP) structures is challenging
partially due to lack of sufficient solved structures for homology modeling.
Here we describe a high-throughput deep transfer learning method that first
predicts MP contacts by learning from non-membrane proteins (non-MPs) and then
predicting three-dimensional structure models using the predicted contacts as
distance restraints. Tested on 510 non-redundant MPs, our method has contact
prediction accuracy at least 0.18 better than existing methods, predicts
correct folds for 218 MPs (TMscore at least 0.6), and generates
three-dimensional models with RMSD less than 4 Angstrom and 5 Angstrom for 57
and 108 MPs, respectively. A rigorous blind test in the continuous automated
model evaluation (CAMEO) project shows that our method predicted
high-resolution three-dimensional models for two recent test MPs of 210
residues with RMSD close to 2 Angstrom. We estimated that our method could
predict correct folds for between 1,345 and 1,871 reviewed human multi-pass MPs
including a few hundred new folds, which shall facilitate the discovery of
drugs targeting at membrane proteins.


  The two dominant approaches for the analysis of species-habitat associations
in animals have been shown to reach divergent conclusions. Models fitted from
the viewpoint of an individual (step selection functions), once scaled up, do
not agree with models fitted from a population viewpoint (resource selection
functions). We explain this fundamental incompatibility, and propose a solution
by introducing to the animal movement field a novel use for the well-known
family of Markov chain Monte Carlo (MCMC) algorithms. By design, the step
selection rules of MCMC lead to a steady-state distribution that coincides with
a given underlying function: the target distribution. We therefore propose an
analogy between the movements of an animal and the movements of a MCMC sampler,
to guarantee convergence of the step selection rules to the parameters
underlying the population's utilisation distribution. We introduce a
rejection-free MCMC algorithm, the local Gibbs sampler, that better resembles
real animal movement, and discuss the wide range of biological assumptions that
it can accommodate. We illustrate our method with simulations on a known
utilisation distribution, and show theoretically and empirically that locations
simulated from the local Gibbs sampler give rise to the correct resource
selection function. Using simulated data, we demonstrate how this framework can
be used to estimate resource selection and movement parameters.


  This work deals with a system of three distinct species that changes in time
under the presence of mobility, selection, and reproduction, as in the popular
rock-paper-scissors game. The novelty of the current study is the modification
of the mobility rule to the case of directional mobility, in which the species
move following the direction associated to a larger (averaged) number density
of selection targets in the surrounding neighborhood. Directional mobility can
be used to simulate eyes that see or a nose that smells, and we show how it may
contribute to reduce the probability of coexistence.


  Random matrix theory successfully connects the structure of interactions of
large ecological communities to their ability to respond to perturbations. One
of the most debated aspects of this approach is the missing role of population
abundances. Despite being one of the most studied patterns in ecology, and one
of the most empirically accessible quantities, population abundances are always
neglected in random matrix approaches and their role in determining stability
is still not understood. Here, we tackle this question by explicitly including
population abundances in a random matrix framework. We obtain an analytical
formula that describes the spectrum of a large community matrix for arbitrary
feasible species abundance distributions. The emerging picture is remarkably
simple: while population abundances affect the rate to return to equilibrium
after a perturbation, the stability of large ecosystems is uniquely determined
by the interaction matrix. We confirm this result by showing that the
likelihood of having a feasible and unstable solution in the Lotka-Volterra
system of equations decreases exponentially with the number of species for
stable interaction matrices.


  Environment plays a fundamental role in the competition for resources, and
hence in the evolution of populations. Here, we study a well-mixed, finite
population consisting of two strains competing for the limited resources
provided by an environment that randomly switches between states of abundance
and scarcity. Assuming that one strain grows slightly faster than the other, we
consider two scenarios--one of pure resource competition, and one in which one
strain provides a public good--and investigate how environmental randomness
(external noise) coupled to demographic (internal) noise determines the
population's fixation properties and size distribution. By analytical means and
simulations, we show that these coupled sources of noise can significantly
enhance the fixation probability of the slower-growing species. We also show
that the population size distribution can be unimodal, bimodal or multimodal
and undergoes noise-induced transitions between these regimes when the rate of
switching matches the population's growth rate.


  Phylogenetic species trees typically represent the speciation history as a
bifurcating tree. Speciation events that simultaneously create more than two
descendants, thereby creating polytomies in the phylogeny, are possible.
Moreover, the inability to resolve relationships is often shown as a (soft)
polytomy. Both types of polytomies have been traditionally studied in the
context of gene tree reconstruction from sequence data. However, polytomies in
the species tree cannot be detected or ruled out without considering gene tree
discordance. In this paper, we describe a statistical test based on properties
of the multi-species coalescent model to test the null hypothesis that a branch
in an estimated species tree should be replaced by a polytomy. On both
simulated and biological datasets, we show that the null hypothesis is rejected
for all but the shortest branches, and in most cases, it is retained for true
polytomies. The test, available as part of the ASTRAL package, can help
systematists decide whether their datasets are sufficient to resolve specific
relationships of interest.


  Modern biological techniques enable very dense genetic sampling of unfolding
evolutionary histories, and thus frequently sample some genotypes multiple
times. This motivates strategies to incorporate genotype abundance information
in phylogenetic inference. In this paper, we synthesize a stochastic process
model with standard sequence-based phylogenetic optimality, and show that tree
estimation is substantially improved by doing so. Our method is validated with
extensive simulations and an experimental single-cell lineage tracing study of
germinal center B cell receptor affinity maturation.


  This report describes the development of an aptamer for sensing azole
antifungal drugs for therapeutic drug monitoring. Modified Synthetic Evolution
of Ligands through Exponential Enrichment (SELEX) was used to discover a DNA
aptamer recognizing azole class antifungal drugs. This aptamer undergoes a
secondary structural change upon binding to its target molecule as shown
through fluorescence anisotropy-based binding measurements. Experiments using
circular dichroism spectroscopy, revealed a unique double G-quadruplex
structure that was essential and specific for binding to the azole antifungal
target. Aptamer-functionalized Graphene Field Effect Transistor (GFET) devices
were created and used to measure the binding of strength of azole antifungals
to this surface. In total this aptamer and the supporting sensing platform
could provide a valuable tool for improving the treatment of patients with
invasive fungal infections.


  Heterogeneity in gene expression across isogenic cell populations can give
rise to phenotypic diversity, even when cells are in homogenous environments.
This diversity arises from the discrete, stochastic nature of biochemical
reactions, which naturally arise due to the very small numbers of genes, RNA,
or protein molecules in single cells. Modern measurements of single
biomolecules have created a vast wealth of information about the fluctuations
of these molecules, but a quantitative understanding of these complex,
stochastic systems requires precise computational tools. In this article, we
present modern tools necessary to model variability in biological system and to
compare model results to experimental data. We review the Chemical Master
Equation and approaches to solve for probability distributions of discrete
numbers of biomolecules. We discuss how to fit probability distributions to
single-cell data using likelihood based approaches. Finally, we provide
examples of fitting discrete stochastic models to single-molecule fluorescent
in-situ hybridization data that quantifies RNA levels across populations of
cells in bacteria and yeast.


  Due to their antimicrobial properties, silver nanoparticles (AgNPs) are being
used in non-edible and edible consumer products. It is not clear though if
exposure to these chemicals can exert toxic effects on the host and gut
microbiome. Conflicting studies have been reported on whether AgNPs result in
gut dysbiosis and other changes within the host. We sought to examine whether
exposure of Sprague-Dawley male rats for two weeks to different shapes of
AgNPs, cube (AgNC) and sphere (AgNS) affects gut microbiota, select behaviors,
and induces histopathological changes in the gastrointestinal system and brain.
In the elevated plus maze (EPM), AgNS-exposed rats showed greater number of
entries into closed arms and center compared to controls and those exposed to
AgNC. AgNS and AgNC treated groups had select reductions in gut microbiota
relative to controls. Clostridium spp., Bacteroides uniformis,
Christensenellaceae, and Coprococcus eutactus were decreased in AgNC exposed
group, whereas, Oscillospira spp., Dehalobacterium spp., Peptococcaeceae,
Corynebacterium spp., Aggregatibacter pneumotropica were reduced in AgNS
exposed group. Bacterial reductions correlated with select behavioral changes
measured in the EPM. No significant histopathological changes were evident in
the gastrointestinal system or brain. Findings suggest short-term exposure to
AgNS or AgNC can lead to behavioral and gut microbiome changes.


  Gene regulatory networks (GRNs) control cellular function and decision making
during tissue development and homeostasis. Mathematical tools based on
dynamical systems theory are often used to model these networks, but the size
and complexity of these models mean that their behaviour is not always
intuitive and the underlying mechanisms can be difficult to decipher. For this
reason, methods that simplify and aid exploration of complex networks are
necessary. To this end we develop a broadly applicable form of the Zwanzig-Mori
projection. By first converting a thermodynamic state ensemble model of gene
regulation into mass action reactions we derive a general method that produces
a set of time evolution equations for a subset of components of a network. The
influence of the rest of the network, the bulk, is captured by memory functions
that describe how the subnetwork reacts to its own past state via components in
the bulk. These memory functions provide probes of near-steady state dynamics,
revealing information not easily accessible otherwise. We illustrate the method
on a simple cross-repressive transcriptional motif to show that memory
functions not only simplify the analysis of the subnetwork but also have a
natural interpretation. We then apply the approach to a GRN from the vertebrate
neural tube, a well characterised developmental transcriptional network
composed of four interacting transcription factors. The memory functions reveal
the function of specific links within the neural tube network and identify
features of the regulatory structure that specifically increase the robustness
of the network to initial conditions. Taken together, the study provides
evidence that Zwanzig-Mori projections offer powerful and effective tools for
simplifying and exploring the behaviour of GRNs.


  We consider stochastically modeled reaction networks and prove that if a
constant solution to the Kolmogorov forward equation decays fast enough
relatively to the transition rates, then the model is non-explosive. In
particular, complex balanced reaction networks are non-explosive.


  Biochemical reactions typically occur at low copy numbers, but at once in
crowded and diverse environments. Space and stochasticity therefore play an
essential role in biochemical networks. Spatial-stochastic simulations have
become a prominent tool for understanding how stochasticity at the microscopic
level influences the macroscopic behavior of such systems. However, while
particle-based models guarantee the level of detail necessary to accurately
describe the microscopic dynamics at very low copy numbers, the algorithms used
to simulate them oftentimes imply trade-offs between computational efficiency
and accuracy. eGFRD (enhanced Green's Function Reaction Dynamics) is an exact
algorithm that evades such trade-offs by partitioning the N-particle system
into M<N analytically tractable one- and two-particle systems; the analytical
solutions (Green's functions) then are used to implement an event-driven
particle-based scheme that allows particles to make large jumps in time and
space while retaining access to their state variables at any moment. Here we
present "eGFRD2", a new eGFRD version that implements the principle of eGFRD in
all dimensions, enabling efficient simulation of biochemical reaction-diffusion
processes in the 3D cytoplasm, on 2D planes representing membranes, and on 1D
elongated cylinders representative of, e.g., cytoskeletal tracks or DNA; in 1D,
it also incorporates convective motion used to model active transport. We find
that, for low particle densities, eGFRD2 is up to 3 orders of magnitude faster
than optimized Brownian Dynamics. We exemplify the capabilities of eGFRD2 by
simulating an idealized model of Pom1 gradient formation, which involves 3D
diffusion, active transport on microtubules, and autophosphorylation on the
membrane, confirming recent results on this system and demonstrating that it
can efficiently operate under genuinely stochastic conditions.


  The honey bees play a role of unquestioned relevance in nature and the
comprehension of the mechanisms affecting their population dynamic is of
fundamental importance. As experimentally documented, the proper development of
a colony is related to the nest temperature, whose value is maintained around
the optimal value if the colony population is sufficiently large. Then, the
environmental temperature, the way in which this influence the nest temperature
and the colony population size, are variables closely linked to each other and
deserve to be taken into account in a model that aims to describe the
population dynamics. In the present study, as first step, the continuous-time
autonomous system proposed by Khoury, Myerscoug and Barron (KMB) in 2011 was
approximated by means a Nonstandard finite difference (NSFD) scheme in order to
obtain a set of autonomous difference equations. Subsequently, with the aim to
introduce the seasonal effects, a nonautonomous version (NAKMB) was proposed
and formulated in discrete-time domain via a NSFD scheme, by introducing a
time-dependent formulation for the queen bee laying rate and the recruitment
rate coefficients. By means the phase-plane analysis was possible to deduce
that, with an appropriate choice of the parameters, the NAKMB model admits both
a limit cycle at nonzero population size and an equilibrium point marking the
colony collapse, depending on the initial population size.


  We will study a mathematical model of the human immunodeficiency virus (HIV)
infection in the presence of combination therapy that includes within-host
infectious dynamics. The deterministic model requires us to analyze asymptotic
stability of two distinct steady states, disease-free and endemic equilibria.
Previous results have focused on investigating the global asymptotic stability
of the trivial steady state using an implicit finite-difference method which
generates a system of difference equations. We, instead, provide analytic
solutions and long term attractive behavior for the endemic steady state using
the theory of difference equations. The dynamics of estimated model is
appropriately determined by a certain quantity threshold maintaining the immune
response to a sufficient level. The result also indicates that a forward
bifurcation in the model happens when the disease-free equilibrium loses its
stability and a stable endemic equilibrium appears as the basic reproduction
number exceeds unity. In this scenario, the classical requirement of the
reproduction number being less than unity becomes a necessary and sufficient
condition for disease mitigation. When the associated reproduction number is in
excess of unity, a stable endemic equilibrium emerges with an unstable
disease-free equilibrium (leading to the persistence and existence of HIV
within the infected individuals). The attractivity of the model reveals that
the disease-free equilibrium is globally asymptotically stable under certain
assumptions. A comparison between the continuous and estimated discrete models
is also provided to have a clear perception in understanding the behavioral
dynamics of disease modelling. Finally, we show that the associated estimation
method is very robust in the sense of numerical stability since the equilibria
and the stability conditions are independent of the time step.


  Demographic (shot) noise in population dynamics scales with the square root
of the population size. This process is very important, as it yields an
absorbing state at zero field, but simulating it, especially on spatial
domains, is a non-trivial task. Here we compare the results of two
operator-splitting techniques suggested for simulating the corresponding
Langevin equation, one by Pechenik and Levine (PL) and the other by Dornic,
Chat\'e and Mu\~noz (DCM). We identify an anomalously strong bias toward the
active phase in the numerical scheme of DCM, a bias which is not present in the
alternative scheme of PL. This bias strongly distorts the phase diagram
determined via the DCM procedure for the range of time-steps used in such
simulations. We pinpoint the underlying cause in the inclusion of the
diffusion, treated as an on-site decay with a constant external source, in the
stochastic part of the algorithm. Treating the diffusion deterministically is
shown to remove this unwanted bias while keeping the simulation algorithm
stable, thus a hybrid numerical technique, in which the DCM approach to
diffusion is applied but the diffusion is simulated deterministically, appears
to be optimal.


  We study the fixation probability of a mutant type when introduced into a
resident population. As opposed to the usual assumption of constant pop-
ulation size, we allow for stochastically varying population sizes. This is
implemented by a stochastic competitive Lotka-Volterra model. The compe- tition
coefficients are interpreted in terms of inverse payoffs emerging from an
evolutionary game. Since our study focuses on the impact of the competition
values, we assume the same birth and death rates for both types. In this gen-
eral framework, we derive an approximate formula for the fixation probability
{\phi} of the mutant type under weak selection. The qualitative behavior of
{\phi} when compared to the neutral scenario is governed by the invasion
dynamics of an initially rare type. Higher payoffs when competing with the
resident type yield higher values of {\phi}. Additionally, we investigate the
influence of the remaining parameters and find an explicit dependence of {\phi}
on the mixed equilibrium value of the corresponding deterministic system (given
that the parameter values allow for its existence).


  Mathematical modelling and numerical simulations of interaction populations
are crucial topics in systems biology. The interactions of ecological models
may occur among individuals of the same species or individuals of different
species. Describing the dynamics of such models occasionally requires some
techniques of model analysis. Choosing appropriate techniques of model analysis
is often a difficult task. We define a prey (mouse) and predator (cat) model.
The system is modelled by a pair of non-linear ordinary differential equations
using mass action law, under constant rates. A proper scaling is suggested to
minimize the number of parameters. More interestingly, we propose a homotopy
technique with n expanding parame- ters for finding some analytical approximate
solutions. Furthermore, using the local sensitivity method is another important
step forward in this study because it helps to identify critical model
parameters. Numerical simulations are provided using Matlab for different
parameters and initial conditions.


  Evolutionary games on graphs have been extensively studied. A variety of
graph structures, graph dynamics, and behaviours of replicators have been
explored. These models have primarily been studied in the framework of
facilitation of cooperation, and much previous research has shed light on this
field of study. However, there has been little attention devoted to truncation
selection as most models employ proportional selection (such as in the
replicator equation) or `imitate the best.' Here we systematically explore
truncation selection on periodic square lattices, where replicators below a
fitness threshold are culled and the reproduction probabilities are equal for
all survivors. We employ two variations of this method: independent truncation,
where the threshold is fixed; and dependent truncation, which is a
generalization of `imitate the best.' Further, we explore the effects of
diffusion in our networks in the following orders of operation:
contest-diffusion-offspring (CDO), and diffusion-contest-offspring (DCO). CDO
and DCO frequently facilitate and diminish cooperation, respectively. For
independent truncation, we find three regimes determined by the fitness
threshold: cooperation decreases as we raise the threshold; polymorphisms and
extinction can occur; and the entire population goes extinct. Further, we show
how an intermediate sucker's payoff maximizes cooperation in the DCO
independent truncation model. We find that dependent truncation affects games
differently; lower levels reduce cooperation for the Hawk Dove game and
increase it for the Stag Hunt, and higher levels produce the opposite effects.
We compare these truncation methods to proportional selection, and show that
they can facilitate cooperation. We conclude that truncation selection can
impact the prevalence of cooperation in complex ways, and therefore merit
further study.


  We present kleuren, a novel assembly-free method to reconstruct phylogenetic
trees using the Colored de Bruijn Graph. kleuren works by constructing the
Colored de Bruijn Graph and then traversing it, finding bubble structures in
the graph that provide phylogenetic signal. The bubbles are then aligned and
concatenated to form a supermatrix, from which a phylogenetic tree is inferred.
We introduce the algorithms that kleuren uses to accomplish this task, and show
its performance on reconstructing the phylogenetic tree of 12 Drosophila
species. kleuren reconstructed the established phylogenetic tree accurately,
and is a viable tool for phylogenetic tree reconstruction using whole genome
sequences. Software package available at: https://github.com/Colelyman/kleuren


  Natural selection explains how life has evolved over millions of years from
more primitive forms. The speed at which this happens, however, has sometimes
defied formal explanations when based on random (uniformly distributed)
mutations. Here we investigate the application of a simplicity bias based on a
natural but algorithmic distribution of mutations (no recombination) in various
examples, particularly binary matrices in order to compare evolutionary
convergence rates. Results both on synthetic and on small biological examples
indicate an accelerated rate when mutations are not statistical uniform but
\textit{algorithmic uniform}. We show that algorithmic distributions can evolve
modularity and genetic memory by preservation of structures when they first
occur sometimes leading to an accelerated production of diversity but also
population extinctions, possibly explaining naturally occurring phenomena such
as diversity explosions (e.g. the Cambrian) and massive extinctions (e.g. the
End Triassic) whose causes are currently a cause for debate. The natural
approach introduced here appears to be a better approximation to biological
evolution than models based exclusively upon random uniform mutations, and it
also approaches a formal version of open-ended evolution based on previous
formal results. These results validate some suggestions in the direction that
computation may be an equally important driver of evolution. We also show that
inducing the method on problems of optimization, such as genetic algorithms,
has the potential to accelerate convergence of artificial evolutionary
algorithms.


  Biological systems reach organizational complexity that far exceeds the
complexity of any known inanimate objects. Biological entities undoubtedly obey
the laws of quantum physics and statistical mechanics. However, is modern
physics sufficient to adequately describe, model and explain the evolution of
biological complexity? Detailed parallels have been drawn between statistical
thermodynamics and the population-genetic theory of biological evolution. Based
on these parallels, we outline new perspectives on biological innovation and
major transitions in evolution, and introduce a biological equivalent of
thermodynamic potential that reflects the innovation propensity of an evolving
population. Deep analogies have been suggested to also exist between the
properties of biological entities and processes, and those of frustrated states
in physics, such as glasses. We extend such analogies by examining
frustration-type phenomena, such as conflicts between different levels of
selection, in biological evolution. We further address evolution in
multidimensional fitness landscapes from the point of view of percolation
theory and suggest that percolation at level above the critical threshold
dictates the tree-like evolution of complex organisms. Taken together, these
multiple connections between fundamental processes in physics and biology imply
that construction of a meaningful physical theory of biological evolution might
not be a futile effort.


  Sand--bubblers are crabs of the genera Dotilla and Scopimera which are known
to produce remarkable patterns and structures at tropical beaches. From these
pattern-making abilities, we may draw inspiration for digital visual art. A
simple mathematical model is proposed and an algorithm is designed that may
create such sand-bubbler patterns artificially. In addition, design parameters
to modify the patterns are identified and analyzed by computational aesthetic
measures. Finally, an extension of the algorithm is discussed that may enable
controlling and guiding generative evolution of the art-making process.


  We present and explore a general method for deriving a Lie-Markov model from
a finite semigroup. If the degree of the semigroup is $k$, the resulting model
is a continuous-time Markov chain on $k$ states and, as a consequence of the
product rule in the semigroup, satisfies the property of multiplicative
closure. This means that the product of any two probability substitution
matrices taken from the model produces another substitution matrix also in the
model. We show that our construction is a natural generalization of the concept
of group-based models.


  In this work we construct individual-based models that give rise to the
generalized logistic model at the mean-field deterministic level and that allow
us to interpret the parameters of these models in terms of individual
interactions. We also study the effect of internal fluctuations on the
long-time dynamics for the different models that have been widely used in the
literature, such as the theta-logistic and Savageau models. In particular, we
determine the conditions for population extinction and calculate the mean time
to extinction. If the population does not become extinct, we obtain analytical
expressions for the population abundance distribution. Our theoretical results
are based on WKB theory and the probability generating function formalism and
are verified by numerical simulations.


  Understanding under which conditions the increase of systems complexity is
evolutionary advantageous, and how this trend is related to the modulation of
the intrinsic noise, are fascinating issues of utmost importance for synthetic
and systems biology. To get insights into these matters, we analyzed chemical
reaction networks with different topologies and degrees of complexity,
interacting or not with the environment. We showed that the global level of
fluctuations at the steady state, as measured by the sum of the Fano factors of
the number of molecules of all species, is directly related to the topology of
the network. For systems with zero deficiency, this sum is constant and equal
to the rank of the network. For higher deficiencies, we observed an increase or
decrease of the fluctuation levels according to the values of the reaction
fluxes that link internal species, multiplied by the associated stoichiometry.
We showed that the noise is reduced when the fluxes all flow towards the
species of higher complexity, whereas it is amplified when the fluxes are
directed towards lower complexity species.


  A key parameter in models for the spread of infectious diseases is the basic
reproduction number $R_0$, which is the expected number of secondary cases a
typical infected primary case infects during its infectious period in a large
mostly susceptible population. In order for this quantity to be meaningful, the
initial expected growth of the number of infectious individuals in the
large-population limit should be exponential.
  We investigate to what extent this assumption is valid by performing repeated
simulations of epidemics on selected empirical networks, viewing each epidemic
as a random process in discrete time. The initial phase of each epidemic is
analyzed by fitting the number of infected people at each time step to a
generalised growth model, allowing for estimating the shape of the growth. For
reference, similar investigations are done on some elementary graphs such as
integer lattices in different dimensions and configuration model graphs, for
which the early epidemic behaviour is known.
  We find that for the empirical networks tested in this paper, exponential
growth characterizes the early stages of the epidemic, except when the network
is restricted by a strong low-dimensional spacial constraint, such as is the
case for the two-dimensional square lattice. However, on finite integer
lattices of sufficiently high dimension, the early development of epidemics
shows exponential growth.


  Short-scale interactions yield large-scale vegetation patterns that, in turn,
shape ecosystem function across landscapes. Fairy circles, which are circular
patches bare of vegetation within otherwise continuous landscapes, are
characteristic features of semiarid grasslands. We report the occurrence of
submarine fairy circle seascapes in seagrass meadows and propose a simple model
that reproduces the diversity of seascapes observed in these ecosystems as
emerging from plant interactions within the meadow. These seascapes include two
extreme cases, a continuous meadow and a bare landscape, along with
intermediate states that range from the occurrence of persistent but isolated
fairy circles, or solitons, to seascapes with multiple fairy circles, banded
vegetation, and "leopard skin" patterns consisting of bare seascapes patterns
consisting of bare seascapes dotted with plant patches. The model predicts that
these intermediate seascapes extending across kilometers emerge as a
consequence of local demographic imbalances along with facilitative and
competitive interactions among the plants with a characteristic spatial scale
of 20 to 30 m, consistent with known drivers of seagrass performance. The
model, which can be extended to clonal growth plants in other landscapes
showing fairy rings, reveals that the different seascapes observed hold
diagnostic power as to the proximity of seagrass meadows to extinction points
that can be used to identify ecosystems at risks.


  Background. A large number of algorithms is being developed to reconstruct
evolutionary models of individual tumours from genome sequencing data. Most
methods can analyze multiple samples collected either through bulk multi-region
sequencing experiments or the sequencing of individual cancer cells. However,
rarely the same method can support both data types.
  Results. We introduce TRaIT, a computational framework to infer mutational
graphs that model the accumulation of multiple types of somatic alterations
driving tumour evolution. Compared to other tools, TRaIT supports multi-region
and single-cell sequencing data within the same statistical framework, and
delivers expressive models that capture many complex evolutionary phenomena.
TRaIT improves accuracy, robustness to data-specific errors and computational
complexity compared to competing methods.
  Conclusions. We show that the application of TRaIT to single-cell and
multi-region cancer datasets can produce accurate and reliable models of
single-tumour evolution, quantify the extent of intra-tumour heterogeneity and
generate new testable experimental hypotheses.


  Sheep pox is a highly transmissible disease which can cause serious loss of
livestock and can therefore have major economic impact. We present data from
sheep pox epidemics which occurred between 1994 and 1998. The data include
weekly records of infected farms as well as a number of covariates. We
implement Bayesian stochastic regression models which, in addition to various
explanatory variables like seasonal and environmental/meteorological factors,
also contain serial correlation structure based on variants of the
Ornstein-Uhlenbeck process. We take a predictive view in model selection by
utilizing deviance-based measures. The results indicate that seasonality and
the number of infected farms are important predictors for sheep pox incidence.


  Reconciling the stability of epigenetic patterns with the rapid turnover of
histone modifications and their adaptability to external stimuli is an
outstanding challenge. Here, we propose a new biophysical mechanism that can
establish and maintain robust yet plastic epigenetic domains via genomic
bookmarking (GBM). We model chromatin as a recolourable polymer whose segments
bear non-permanent histone marks (or colours) which can be modified by "writer"
proteins. The three-dimensional chromatin organisation is mediated by protein
bridges, or "readers", such as Polycomb Repressive Complexes and Transcription
Factors. The coupling between readers and writers drives spreading of
biochemical marks and sustains the memory of local chromatin states across
replication and mitosis. In contrast, GBM-targeted perturbations destabilise
the epigenetic patterns. Strikingly, we demonstrate that GBM alone can explain
the full distribution of Polycomb marks in a whole Drosophila chromosome. We
finally suggest that our model provides a starting point for an understanding
of the biophysics of cellular differentiation and reprogramming.


  Various approaches to alignment-free sequence comparison are based on the
length of exact or inexact word matches between two input sequences. Haubold
{\em et al.} (2009) showed how the average number of substitutions between two
DNA sequences can be estimated based on the average length of exact common
substrings. In this paper, we study the length distribution of $k$-mismatch
common substrings between two sequences. We show that the number of
substitutions per position that have occurred since two sequences have evolved
from their last common ancestor, can be estimated from the position of a local
maximum in the length distribution of their $k$-mismatch common substrings.


  Noncoding RNA (ncRNA) and long noncoding RNA (lncRNA) are scientifically
invalid terms because they define molecular entities according to properties
they do not possess and functions they do not perform. Here, I suggest retiring
these two terms. Instead, I suggest using an evolutionary classification of
genomic function, in which every RNA molecule is classified as either
"functional" or "junk' according to its selected effect function. Dealing with
RNA molecules whose functional status is unknown require us to phrase Popperian
nomenclatures that spell out the conditions for their own refutation. Thus, in
the absence of falsifying evidence, RNA molecules of unknown function must be
considered junk RNA (jRNA).


  Biological network alignment is currently in a state of disarray, with more
than two dozen network alignment tools having been introduced in the past
decade, with no clear winner, and other new tools being published almost
quarterly. Part of the problem is that almost every new tool proposes both a
new objective function and a new search algorithm to optimize said objective.
These two aspects of alignment are orthogonal, and confounding them makes it
difficult to evaluate them separately. A more systematic approach is needed. To
this end, we bring these two orthogonal issues into sharp focus in two
companion papers. In Part 1 (this paper) we show that simulated annealing, as
implemented by SANA, far outperforms all other existing search algorithms
across a wide range of objectives. Part 2 (our companion paper) then uses SANA
to compare over a dozen objectives in terms of the biology they recover,
demonstrating that some objective functions recover useful biological
information while others do not. We propose that further work should focus on
improving objective functions, with SANA the obvious choice as the search
algorithm.


  Traveling waves describe diverse natural phenomena from crystal growth in
physics to range expansions in biology. Two classes of waves exist with very
different properties: pulled and pushed. Pulled waves are driven by high growth
rates at the expansion edge, where the number of organisms is small and
fluctuations are large. In contrast, fluctuations are suppressed in pushed
waves because the region of maximal growth is shifted towards the population
bulk. Although it is commonly believed that expansions are either pulled or
pushed, we found an intermediate class of waves with bulk-driven growth, but
exceedingly large fluctuations. These waves are unusual because their
properties are controlled by both the leading edge and the bulk of the front.


  We begin with the extraordinary observation that the length distribution of
80 million proteins in UniProt, the Universal Protein Resource, measured in
amino acids, is qualitatively identical to the length distribution of large
collections of computer functions measured in programming language tokens, at
all scales. That two such disparate discrete systems share important structural
properties suggests that yet other apparently unrelated discrete systems might
share the same properties, and certainly invites an explanation.
  We demonstrate that this is inevitable for all discrete systems of components
built from tokens or symbols. Departing from existing work by embedding the
Conservation of Hartley-Shannon information (CoHSI) in a classical statistical
mechanics framework, we identify two kinds of discrete system, heterogeneous
and homogeneous. Heterogeneous systems contain components built from a unique
alphabet of tokens and yield an implicit CoHSI distribution with a sharp
unimodal peak asymptoting to a power-law. Homogeneous systems contain
components each built from just one kind of token unique to that component and
yield a CoHSI distribution corresponding to Zipf's law.
  This theory is applied to heterogeneous systems, (proteome, computer
software, music); homogeneous systems (language texts, abundance of the
elements); and to systems in which both heterogeneous and homogeneous behaviour
co-exist (word frequencies and word length frequencies in language texts). In
each case, the predictions of the theory are tested and supported to high
levels of statistical significance. We also show that in the same heterogeneous
system, different but consistent alphabets must be related by a power-law. We
demonstrate this on a large body of music by excluding and including note
duration in the definition of the unique alphabet of notes.


  The folding of a protein towards its native state is a rather complicated
process. However there are empirical evidences that the folding time correlates
with the contact order, a simple measure of the spatial organisation of the
native state of the protein. Contact order is related to the average length of
the main chain loops formed by amino acids which are in contact. Here we argue
that folding kinetics can be influenced also by the entanglement that loops may
undergo within the overall three dimensional protein structure. In order to
explore such possibility, we introduce a novel descriptor, which we call
"maximum intrachain contact entanglement". Specifically, we measure the maximum
Gaussian entanglement between any looped portion of a protein and any other
non-overlapping subchain of the same protein, which is easily computed by
discretized line integrals on the coordinates of the $C_{\alpha}$ atoms. By
analyzing experimental data sets of two-state and multistate folders, we show
that also the new index is a good predictor of the folding rate. Moreover,
being only partially correlated with previous methods, it can be integrated
with them to yield more accurate predictions.


  In this paper we present a model to estimate the density of aedes mosquitoes
in a community affected by dengue. The model is based on the fitting of a
continuous function to the incidence of dengue infections, from which the
density of infected mosquitoes is derived straightforwardly. Further
derivations allows the calculation of the latent and susceptible mosquitoes'
densities, the sum of the three equals the total mosquitoes' density. The model
is illustrated with the case of the risk of urban yellow fever resurgence in
dengue infested areas but the same methods apply for other aedes-transmitted
infections like Zika and chikungunya viruses.


  Biological networks provide insight into the complex organization of
biological processes in a cell at the system level. They are an effective tool
for understanding the comprehensive map of functional interactions, finding the
functional modules and pathways. Reconstruction and comparative analysis of
these networks provide useful information to identify functional modules,
prioritization of disease causing genes and also identification of drug
targets. The talk will consist of two parts. I will discuss several methods for
protein-protein interaction network alignment and investigate their preferences
to other existing methods. Further, I briefly talk about reconstruction of
protein-protein interaction networks by using deep learning.


  Starting from the most general formulation of stochastic
thermodynamics---i.e. a thermodynamically consistent nonautonomous stochastic
dynamics describing systems in contact with several reservoirs---, we define a
procedure to identify the conservative and the minimal set of nonconservative
contributions in the entropy production. The former is expressed as the
difference between changes caused by time-dependent drivings and a generalized
potential difference. The latter is a sum over the minimal set of flux--force
contributions controlling the dissipative flows across the system. When the
system is initially prepared at equilibrium (e.g. by turning off drivings and
forces), a finite-time detailed fluctuation theorem holds for the different
contributions. Our approach relies on identifying the complete set of conserved
quantities and can be viewed as the extension of the theory of generalized
Gibbs ensembles to nonequilibrium situations.


  We propose a probabilistic model for interpreting gene expression levels that
are observed through single-cell RNA sequencing. In the model, each cell has a
low-dimensional latent representation. Additional latent variables account for
technical effects that may erroneously set some observations of gene expression
levels to zero. Conditional distributions are specified by neural networks,
giving the proposed model enough flexibility to fit the data well. We use
variational inference and stochastic optimization to approximate the posterior
distribution. The inference procedure scales to over one million cells, whereas
competing algorithms do not. Even for smaller datasets, for several tasks, the
proposed procedure outperforms state-of-the-art methods like ZIFA and
ZINB-WaVE. We also extend our framework to account for batch effects and other
confounding factors, and propose a Bayesian hypothesis test for differential
expression that outperforms DESeq2.


  Laboratory models are often used to understand the interaction of related
pathogens via host immunity. For example, recent experiments where ferrets were
exposed to two influenza strains within a short period of time have shown how
the effects of cross-immunity vary with the time between exposures and the
specific strains used. On the other hand, studies of the workings of different
arms of the immune response, and their relative importance, typically use
experiments involving a single infection. However, inferring the relative
importance of different immune components from this type of data is
challenging. Using simulations and mathematical modelling, here we investigate
whether the sequential infection experiment design can be used not only to
determine immune components contributing to cross-protection, but also to gain
insight into the immune response during a single infection.
  We show that virological data from sequential infection experiments can be
used to accurately extract the timing and extent of cross-protection. Moreover,
the broad immune components responsible for such cross-protection can be
determined. Such data can also be used to infer the timing and strength of some
immune components in controlling a primary infection, even in the absence of
serological data. By contrast, single infection data cannot be used to reliably
recover this information. Hence, sequential infection data enhances our
understanding of the mechanisms underlying the control and resolution of
infection, and generates new insight into how previous exposure influences the
time course of a subsequent infection.


  Background: Convolutional Neural Networks can be effectively used only when
data are endowed with an intrinsic concept of neighbourhood in the input space,
as is the case of pixels in images. We introduce here Ph-CNN, a novel deep
learning architecture for the classification of metagenomics data based on the
Convolutional Neural Networks, with the patristic distance defined on the
phylogenetic tree being used as the proximity measure. The patristic distance
between variables is used together with a sparsified version of
MultiDimensional Scaling to embed the phylogenetic tree in a Euclidean space.
Results: Ph-CNN is tested with a domain adaptation approach on synthetic data
and on a metagenomics collection of gut microbiota of 38 healthy subjects and
222 Inflammatory Bowel Disease patients, divided in 6 subclasses.
Classification performance is promising when compared to classical algorithms
like Support Vector Machines and Random Forest and a baseline fully connected
neural network, e.g. the Multi-Layer Perceptron. Conclusion: Ph-CNN represents
a novel deep learning approach for the classification of metagenomics data.
Operatively, the algorithm has been implemented as a custom Keras layer taking
care of passing to the following convolutional layer not only the data but also
the ranked list of neighbourhood of each sample, thus mimicking the case of
image data, transparently to the user. Keywords: Metagenomics; Deep learning;
Convolutional Neural Networks; Phylogenetic trees


  RNA function is intimately related to its structural dynamics. Molecular
dynamics simulations are useful for exploring biomolecular flexibility but are
severely limited by the accessible timescale. Enhanced sampling methods allow
this timescale to be effectively extended in order to probe
biologically-relevant conformational changes and chemical reactions. Here, we
review the role of enhanced sampling techniques in the study of RNA systems. We
discuss the challenges and promises associated with the application of these
methods to force-field validation, exploration of conformational landscapes and
ion/ligand-RNA interactions, as well as catalytic pathways. Important technical
aspects of these methods, such as the choice of the biased collective variables
and the analysis of multi-replica simulations, are examined in detail. Finally,
a perspective on the role of these methods in the characterization of RNA
dynamics is provided.


  The increasing demand for renewable energy is projected to result in a
40-fold increase in offshore wind electricity in the European Union by 2030.
Despite a great number of local impact studies for selected marine populations,
the regional ecosystem impacts of offshore wind farm structures are not yet
well assessed nor understood. Our study investigates whether the accumulation
of epifauna, dominated by the filter feeder Mytilus edulis (blue mussel), on
turbine structures affects pelagic primary productivity and ecosystem
functioning in the southern North Sea. We estimate the anthropogenically
increased potential distribution based on the current projections of turbine
locations and reported patterns of M. edulis settlement. This distribution is
integrated through the Modular Coupling System for Shelves and Coasts to
state-of-the-art hydrodynamic and ecosystem models. Our simulations reveal
non-negligible potential changes in regional annual primary productivity of up
to 8% within the offshore wind farm area, and induced maximal increases of the
same magnitude in daily productivity also far from the wind farms. Our setup
and modular coupling are effective tools for system scale studies of other
environmental changes arising from large-scale offshore wind-farming such as
ocean physics and distributions of pelagic top predators.


  We study a membrane tube with unidirectional ion pumps driving an osmotic
pressure difference. A pressure driven peristaltic instability is identified,
qualitatively distinct from similar tension-driven Rayleigh type instabilities
on membrane tubes. We discuss how this instability could be related to the
function and biogenesis of membrane bound organelles, in particular the
contractile vacuole complex. The unusually long natural wavelength of this
instability is in agreement with that observed in cells.


  Many mathematical models of evolution assume that all individuals experience
the same environment. Here, we study the Moran process in heterogeneous
environments. The population is of finite size with two competing types, which
are exposed to a fixed number of environmental conditions. Reproductive rate is
determined by both the type and the environment. We first calculate the
condition for selection to favor the mutant relative to the resident wild type.
In large populations, the mutant is favored if and only if the mutant's spatial
average reproductive rate exceeds that of the resident. But environmental
heterogeneity elucidates an interesting asymmetry between the mutant and the
resident. Specifically, mutant heterogeneity suppresses its fixation
probability; if this heterogeneity is strong enough, it can even completely
offset the effects of selection (including in large populations). In contrast,
resident heterogeneity has no effect on a mutant's fixation probability in
large populations and can amplify it in small populations.


  Single molecule force spectroscopy of DNA strands adsorbed at surfaces is a
powerful technique used in air or liquid environments to quantify their
mechanical properties. Although the force responses are limited to unfolding
events so far, single base detection might be possible in more drastic
cleanliness conditions such as ultra high vacuum. Here, we report on high
resolution imaging and pulling attempts at low temperature (5K) of a single
strand DNA (ssDNA) molecules composed of 20 cytosine bases adsorbed on Au(111)
by scanning probe microscopy and numerical calculations. Using electrospray
deposition technique, the ssDNA were successfully transferred from solution
onto a surface kept in ultra high vacuum. Real space characterizations reveal
that the ssDNA have an amorphous structure on gold in agreement with numerical
calculations. Subsequent substrate annealing promotes the desorption of solvent
molecules, DNA as individual molecules as well as the formation of DNA self
assemblies. Furthermore, pulling experiments by force spectroscopy have been
conducted to measure the mechanical response of the ssDNA while detaching. A
periodic pattern of 0.2 to 0.3nm is observed in the force curve which arises
from the stick slip of single nucleotide bases over the gold. Although an intra
molecular response is obtained in the force curve, a clear distinction of each
nucleotide detachment is not possible due the complex structure of ssDNA
adsorbed on gold.


  Coextinction models are useful to understand community robustness to species
loss and resilience to disturbances. We simulated pollinator extinctions in
pollination networks by using a hybrid model that combined a recently developed
stochastic coextinction model (SCM) for plant extinctions and a topological
model (TCM) for animal extinctions. Our model accounted for variation in
interaction strengths and included empirical estimates of plant dependence on
pollinators to set seeds. The stochastic nature of such model allowed us
determining plant survival to single (and multiple) extinction events, and
identifying which pollinators (keystone species) were more likely to trigger
secondary extinctions. Consistently across three different pollinator removal
sequences, plant robustness was lower than in a pure TCM, and plant survival
was more determined by dependence on the mutualism than by interaction
strength. As expected, highly connected and dependent plants were the most
sensitive to pollinator loss and collapsed faster in extinction cascades. We
predict that the relationship between dependence and plant connectivity is
crucial to determine network robustness to interaction loss. Finally, we showed
that honeybees and several beetles were keystone species in our communities.
This information is of great value to foresee consequences of pollinator losses
facing current global change and to identify target species for effective
conservation.


  Many mathematical frameworks of evolutionary game dynamics assume that the
total population size is constant and that selection affects only the relative
frequency of strategies. Here, we consider evolutionary game dynamics in an
extended Wright-Fisher process with variable population size. In such a
scenario, it is possible that the entire population becomes extinct. Survival
of the population may depend on which strategy prevails in the game dynamics.
Studying cooperative dilemmas, it is a natural feature of such a model that
cooperators enable survival, while defectors drive extinction. Although
defectors are favored for any mixed population, random drift could lead to
their elimination and the resulting pure-cooperator population could survive.
On the other hand, if the defectors remain, then the population will quickly go
extinct because the frequency of cooperators steadily declines and defectors
alone cannot survive. In a mutation-selection model, we find that (i) a steady
supply of cooperators can enable long-term population survival, provided
selection is sufficiently strong, and (ii) selection can increase the abundance
of cooperators but reduce their relative frequency. Thus, evolutionary game
dynamics in populations with variable size generate a multifaceted notion of
what constitutes a trait's long-term success.


  Genome-wide association studies (GWA studies or GWAS) investigate the
relationships between genetic variants such as single-nucleotide polymorphisms
(SNPs) and individual traits. Recently, incorporating biological priors
together with machine learning methods in GWA studies has attracted increasing
attention. However, in real-world, nucleotide-level bio-priors have not been
well-studied to date. Alternatively, studies at gene-level, for example,
protein--protein interactions and pathways, are more rigorous and legitimate,
and it is potentially beneficial to utilize such gene-level priors in GWAS. In
this paper, we proposed a novel two-level structured sparse model, called
Sparse Group Lasso with Group-level Graph structure (SGLGG), for GWAS. It can
be considered as a sparse group Lasso along with a group-level graph Lasso.
Essentially, SGLGG penalizes the nucleotide-level sparsity as well as takes
advantages of gene-level priors (both gene groups and networks), to identifying
phenotype-associated risk SNPs. We employ the alternating direction method of
multipliers algorithm to optimize the proposed model. Our experiments on the
Alzheimer's Disease Neuroimaging Initiative whole genome sequence data and
neuroimage data demonstrate the effectiveness of SGLGG. As a regression model,
it is competitive to the state-of-the-arts sparse models; as a variable
selection method, SGLGG is promising for identifying Alzheimer's
disease-related risk SNPs.


  This paper investigates the evolution of strategic play where players drawn
from a finite well-mixed population are offered the opportunity to play in a
public goods game. All players accept the offer. However, due to the
possibility of unforeseen circumstances, each player has a fixed probability of
being unable to participate in the game, unlike similar models which assume
voluntary participation. We first study how prescribed stochastic opting-out
affects cooperation in finite populations. Moreover, in the model, cooperation
is favored by natural selection over both neutral drift and defection if return
on investment exceeds a threshold value defined solely by the population size,
game size, and a player's probability of opting-out. Ultimately, increasing the
probability that each player is unable to fulfill her promise of participating
in the public goods game facilitates natural selection of cooperators. We also
use adaptive dynamics to study the coevolution of cooperation and opting-out
behavior. However, given rare mutations minutely different from the original
population, an analysis based on adaptive dynamics suggests that the over time
the population will tend towards complete defection and non-participation, and
subsequently, from there, participating cooperators will stand a chance to
emerge by neutral drift. Nevertheless, increasing the probability of
non-participation decreases the rate at which the population tends towards
defection when participating. Our work sheds light on understanding how
stochastic opting-out emerges in the first place and its role in the evolution
of cooperation.


  Drug resistance and strong contacts actually play crucial roles in epidemic
spread in complex systems. Nevertheless, neither theoretical model or
methodology is proposed to address this. We thus consider an edge-based
epidemic spread model considering the two key ingredients, in which the
contacts are grouped into two classes: strong contacts and normal ones. Next,
we present a unified edge-based compartmental approach to the spread dynamics
on Erd\"{o}s-R\'{e}nyi (ER) networks and validate its results by extensive
numerical simulations. In case that epidemic is totally drug-resistant, we both
numerically and theoretically show a slow outbreak (continuous transition) of
epidemics when number of strong contacts is not enough for the emergence of
null threshold. If the epidemic owns partial resistance, we would observe
evident faster-growing outbreaks (discontinuous transitions) and larger final
epidemic sizes for few strong contacts, instead of emergence of null threshold
with increase of strong contacts. Inhibiting effect of infection threshold,
positive roles of strong contacts and strength of strong contacts in promoting
outbreaks are also approved. Throughout this paper, we could drive exact
predictions through the analytical approach, showing good agreements with
numerical simulations.


  In this study, we assess the influence of solvation on the accuracy and
reliability of nuclear shielding calculations for amino acids in comparison to
experimental data. We focus particularly on the performance of solvation
methods for different protonation states, as biological molecules occur almost
exclusively in aqueous solution and are subject to protonation with pH. We
identify significant shortcomings of current implicit solvent models and
present a hybrid solvation approach that improves agreement with experimental
data by taking into account the presence of direct interactions between amino
acid protonation state and water molecules.