  This manuscript contributes a general and practical framework for casting a
Markov process model of a system at equilibrium as a structural causal model,
and carrying out counterfactual inference. Markov processes mathematically
describe the mechanisms in the system, and predict the system's equilibrium
behavior upon intervention, but do not support counterfactual inference. In
contrast, structural causal models support counterfactual inference, but do not
identify the mechanisms. This manuscript leverages the benefits of both
approaches. We define the structural causal models in terms of the parameters
and the equilibrium dynamics of the Markov process models, and counterfactual
inference flows from these settings. The proposed approach alleviates the
identifiability drawback of the structural causal models, in that the
counterfactual inference is consistent with the counterfactual trajectories
simulated from the Markov process model. We showcase the benefits of this
framework in case studies of complex biomolecular systems with nonlinear
dynamics. We illustrate that, in presence of Markov process model
misspecification, counterfactual inference leverages prior data, and therefore
estimates the outcome of an intervention more accurately than a direct
simulation.


  Approaches to understanding resilience from psychology and sociology
emphasise individuals' agency but obscure systemic factors. Approaches to
understanding resilience stemming from ecology emphasise system dynamics such
as feedbacks but obscure individuals. Approaches from both psychology and
ecology examine the actions or attractors available in the present, but neglect
how actions taken now can affect the configuration of the social-ecological
system in the future. Here, we propose an extension to resilience theory, which
we label 'pathway diversity', that links existing individual, systems and
temporal theories of resilience. In our theory of pathway diversity, resilience
is greater if more actions are currently available and can be maintained or
enhanced into the future. Using a toy model of an agricultural
social-ecological system, we show how pathway diversity could deliver a
context-sensitive method of assessing resilience and guiding planning. Using a
toy state-and-transition model of a poverty trap, we show how pathway diversity
is generally consistent with existing definitions of resilience and can
illuminate long-standing questions about normative and descriptive resilience.
Our results show that pathway diversity advances both theoretical understanding
and practical tools for building resilience.


  The method of Probability Density Profile Analysis has been introduced
previously as a tool to find the best match between a set of experimentally
generated Residual Dipolar Couplings and a set of known protein structures.
While it proved effective on small databases in identifying protein fold
families, and for picking the best result from computational protein folding
tool ROBETTA, for larger data sets, more data is required. Here, the method of
2-D Probability Density Profile Analysis is presented which incorporates paired
RDC data from 2 alignment media for N-H vectors. The method was tested using
synthetic RDC data generated with +/-1 Hz error. The results show that the
addition of information from a second alignment medium makes 2-D PDPA a much
more effective tool that is able to identify a structure from a database of 600
protein fold family representatives.


  Aligning multiple protein structures can yield valuable information about
structural similarities among related proteins, as well as provide insight into
evolutionary relationships between proteins in a family. We have developed an
algorithm (msTALI) for aligning multiple protein structures using biochemical
and biophysical properties, including torsion angles, secondary structure,
hydrophobicity, and surface accessibility. The algorithm is a progressive
alignment algorithm motivated by popular techniques from multiple sequence
alignment. It has demonstrated success in aligning the major structural regions
of a set of proteins from the s/r kinase family. The algorithm was also
successful at aligning functional residues of these proteins. In addition, the
algorithm was also successful in aligning seven members of the acyl carrier
protein family, including both experimentally derived as well as
computationally modeled structures.


  Living organisms exhibit consistent homochirality. It is argued that the
specific, binary choice made is not an accident but is a consequence of parity
violation in the weak interaction expressed by cosmic irradiation. The
secondary muons and pairs are spin- and magnetic moment- polarized and may
introduce a small, net chiral preference when they interact electromagnetically
or quantum mechanically with molecules that have made the transition to
self-replication. Although this preference is likely to be very small, it may
suffice to give a chirally-dominant outcome after billions of replications,
especially if combined with chirally-unbiased conflict between the two choices.
Examples of mechanisms that can manifest the three essential steps of
polarization, preference and domination are presented and some variations and
possible implications are discussed.


  The competitive exclusion principle asserts that coexisting species must
occupy distinct ecological niches (i.e. the number of surviving species can not
exceed the number of resources). An open question is to understand if and how
different resource dynamics affect this bound. Here, we analyze a generalized
consumer resource model with externally supplied resources and show that -- in
contrast to self-renewing resources -- species can occupy only half of all
available environmental niches. This motivates us to construct a new schema for
classifying ecosystems based on species packing properties.


  We have a general knowledge of the principles by which catalysts accelerate
the rate of chemical reactions but no precise understanding of the geometrical
and physical constraints to which their design is subject. To analyze these
constraints, we introduce a minimal model of catalysis based on elastic
networks where the implications of the geometry and flexibility of a catalyst
can be studied systematically. The model demonstrates the relevance and
limitations of the principle of transition-state stabilization: optimal
catalysts are found to have a geometry complementary to the transition state
but a degree of flexibility that non-trivially depends on the parameters of the
reaction as well as on external parameters such as the concentrations of
reactants and products. The results illustrate how simple physical models can
provide valuable insights on the design of catalysts.


  In this paper, we discuss the fitness landscape evolution of permanent
replicator systems using a hypothesis that the specific time of evolutionary
adaptation of the system parameters is much slower than the time of internal
evolutionary dynamics. In other words, we suppose that the extreme principle of
Darwinian evolution based the Fisher's fundamental theorem of natural selection
is valid for the steady-states. Various cases of the evolutionary adaptation
for permanent replicator system are considered.


  In order to respond to environmental signals, cells often use small molecular
circuits to transmit information about their surroundings. Recently, motivated
by concrete examples in signaling and gene regulation, a body of work has
focused on the properties of circuits that function out of equilibrium and
dissipate energy. We briefly review the probabilistic measures of information
and dissipation and use simple models to discuss and illustrate trade-offs
between information and dissipation in biological circuits. We find that
circuits with non-steady state initial conditions can transmit more information
at small readout delays than steady state circuits. The dissipative cost of
this additional information proves marginal compared to the steady state
dissipation. Feedback does not significantly increase the transmitted
information for out of steady state circuits but does decrease dissipative
costs. Lastly, we discuss the case of bursty gene regulatory circuits that even
in the fast switching limit function out of equilibrium.


  Chaperone protein - the most disordered among all protein groups - help RNAs
fold into their functional structure by destabilizing misfolded configurations
or stabilizing the functional ones. But disentangling the mechanism underlying
RNA chaperoning is challenging, mostly due to inherent disorder of the
chaperones and the transient nature of their interactions with RNA. In
particular, it is unclear how specific the interactions are and what role is
played by amino acid charge and polarity patterns. Here, we address these
questions in the RNA chaperone StpA. We adapted direct coupling analysis (DCA)
into the $\alpha\beta$DCA method that can treat in tandem sequences written in
two alphabets, nucleotides and amino acids. With $\alpha\beta$DCA, we could
analyze StpA-RNA interactions and show consistency with a previously proposed
two-pronged mechanism: StpA disrupts specific positions in the group I intron
while globally and loosely binding to the entire structure. Moreover, the
interactions are strongly associated with the charge pattern: negatively
charged regions in the destabilizing StpA N-terminal affect a few specific
positions in the RNA, located in stems and in the pseudoknot. In contrast,
positive regions in the C-terminal contain strongly coupled amino acids that
promote non-specific or weakly-specific binding to the RNA. The present study
opens new avenues to examine the functions of disordered proteins and to design
disruptive proteins based on their charge patterns.


  Many fundamental biological processes are regulated by protein-DNA complexes
called {\it synaptosomes}, which possess multiple interaction sites. Despite
the critical importance of synaptosomes, the mechanisms of their formation
remain not well understood. Because of the multi-site nature of participating
proteins, it is widely believed that their search for specific sites on DNA
involves the formation and breaking of DNA loops and sliding in the looped
configurations. In reality, DNA in live cells is densely covered by other
biological molecules that might interfere with the formation of synaptosomes.
In this work, we developed a theoretical approach to evaluate the role of
obstacles in the target search of multi-site proteins when the formation of DNA
loops and the sliding in looped configurations are possible. Our theoretical
method is based on analysis of a discrete-state stochastic model that uses a
master equations approach and extensive computer simulations. It is found that
the obstacle slows down the search dynamics in the system when DNA loops are
long-lived, but the effect is minimal for short-lived DNA loops. In addition,
the relative positions of the target and the obstacle strongly influence the
target search kinetics. Furthermore, the presence of the obstacle might
increase the noise in the system. These observations are discussed using
physical-chemical arguments. Our theoretical approach clarifies the molecular
mechanisms of formation of protein-DNA complexes with multiple interactions
sites.


  Cortisol, corticosterone and aldosterone activate full-length glucocorticoid
receptor (GR) from elephant shark, a cartilaginous fish belonging to the oldest
group of jawed vertebrates. Activation by aldosterone a mineralocorticoid,
indicates partial divergence of elephant shark GR from the MR. Progesterone
activates elephant shark MR, but not elephant shark GR. Progesterone inhibits
steroid binding to elephant shark GR, but not to human GR. Deletion of the
N-terminal domain (NTD) from elephant shark GR (Truncated GR) reduced the
response to corticosteroids, while truncated and full-length elephant shark MR
had similar responses to corticosteroids. Chimeras of elephant shark GR NTD
fused to MR DBD+LBD had increased activation by corticosteroids and
progesterone compared to full-length elephant shark MR. Elephant shark MR NTD
fused to GR DBD+LBD had similar activation as full-length elephant shark MR,
indicating that activation of human GR by the NTD evolved early in GR
divergence from the MR.


  Inspired by recent experiments on the effects of cytosolic crowders on the
organization of bacterial chromosomes, we consider a "feather-boa" type model
chromosome in the presence of non-additive crowders, encapsulated within a
cylindrical cell. We observe spontaneous emergence of complementary helicity of
the confined polymer and crowders. This feature is reproduced within a
simplified effective model of the chromosome. This latter model further
establishes the occurrence of longitudinal and transverse spatial segregation
transitions between the chromosome and crowders upon increasing crowder size.


  We study the emergence of cooperation in structured populations with any
arrangement of cooperators and defectors on the evolutionary graph. Using
structure coefficients defined for configurations describing such arrangements
of any number of mutants, we provide results for weak selection to favor
cooperation over defection on any regular graph with $N \leq 14$ vertices.
Furthermore, the properties of graphs that particularly promote cooperation are
analyzed. It is shown that the number of graph cycles of certain length is a
good predictor for the values of the structure coefficient, and thus a tendency
to favor cooperation. Another property of particularly cooperation-promoting
regular graphs with a low degree is that they are structured to have blocks
with clusters of mutants that are connected by cut vertices and/or hinge
vertices.


  Learning to optimize has emerged as a powerful framework for various
optimization and machine learning tasks. Current such "meta-optimizers" often
learn in the space of continuous optimization algorithms that are point-based
and uncertainty-unaware. To overcome the limitations, we propose a
meta-optimizer that learns in the algorithmic space of both point-based and
population-based optimization algorithms. The meta-optimizer targets at a
meta-loss function consisting of both cumulative regret and entropy.
Specifically, we learn and interpret the update formula through a population of
LSTMs embedded with sample- and feature-level attentions. Meanwhile, we
estimate the posterior directly over the global optimum and use an uncertainty
measure to help guide the learning process. Empirical results over non-convex
test functions and the protein-docking application demonstrate that this new
meta-optimizer outperforms existing competitors.


  Just like the degrees of human and animal interaction networks, the
distribution of the times between interactions is known to often be
right-skewed and fat-tailed. Both these distributions affect epidemic dynamics
strongly, but, as we show in this Letter, for very different reasons. Whereas
the high degrees of the tail are critical for facilitating epidemics, it is the
small interevent times that control the dynamics of epidemics. We investigate
this effect both analytically and numerically for different versions of the
Susceptible-Infected-Recovered model on different types of networks.


  One of the most interesting, difficult, and potentially useful topics in
computational biology is the inference of gene regulatory networks (GRNs) from
expression data. Although researchers have been working on this topic for more
than a decade and much progress has been made, it remains an unsolved problem
and even the most sophisticated inference algorithms are far from perfect. In
this paper, we review the latest developments in network inference, including
state-of-the-art algorithms like PIDC, Phixer, and more. We also discuss
unsolved computational challenges, including the optimal combination of
algorithms, integration of multiple data sources, and pseudo-temporal ordering
of static expression data. Lastly, we discuss some exciting applications of
network inference in cancer research, and provide a list of useful software
tools for researchers hoping to conduct their own network inference analyses.


  The Jaccard similarity index is an important measure of the overlap of two
sets, widely used in machine learning, computational genomics, information
retrieval, and many other areas. We design and implement SimilarityAtScale, the
first communication-efficient distributed algorithm for computing the Jaccard
similarity among pairs of large datasets. Our algorithm provides an efficient
encoding of this problem into a multiplication of sparse matrices. Both the
encoding and sparse matrix product are performed in a way that minimizes data
movement in terms of communication and synchronization costs. We apply our
algorithm to obtain similarity among all pairs of a set of large samples of
genomes. This task is a key part of modern metagenomics analysis and an
evergrowing need due to the increasing availability of high-throughput DNA
sequencing data. The resulting scheme is the first to enable accurate Jaccard
distance derivations for massive datasets, using largescale distributed-memory
systems. We package our routines in a tool, called GenomeAtScale, that combines
the proposed algorithm with tools for processing input sequences. Our
evaluation on real data illustrates that one can use GenomeAtScale to
effectively employ tens of thousands of processors to reach new frontiers in
large-scale genomic and metagenomic analysis. While GenomeAtScale can be used
to foster DNA research, the more general underlying SimilarityAtScale algorithm
may be used for high-performance distributed similarity computations in other
data analytics application domains.


  Biomedical data, particularly in the field of genomics, has characteristics
which make it challenging for machine learning applications - it can be sparse,
high dimensional and noisy. Biomedical applications also present challenges to
model selection - whilst powerful, accurate predictions are necessary, they
alone are not sufficient for a model to be deemed useful. Due to the nature of
the predictions, a model must also be trustworthy and transparent, empowering a
practitioner with confidence that its use is appropriate and reliable. In this
paper, we propose that this can be achieved through the use of judiciously
built feature sets coupled with Bayesian models, specifically Gaussian
processes. We apply Gaussian processes to drug discovery, using inexpensive
transcriptomic profiles from human cell lines to predict animal kidney and
liver toxicity after treatment with specific chemical compounds. This approach
has the potential to reduce invasive and expensive animal testing during
clinical trials if in vitro human cell line analysis can accurately predict
model animal phenotypes. We compare results across a range of feature sets and
models, to highlight model importance for medical applications.


  The purpose of this roadmap article is to draw attention to a paradigm shift
in our understanding of evolution towards a perspective of
ecological-evolutionary feedback, highlighted through two recent highly
simplified examples of rapid evolution. The first example focuses primarily on
population dynamics: anomalies in population cycles can reflect the influence
of strong selection and the interplay with mutations. The second focuses
primarily on the way in which ecological structure can potentially be
influenced by what is arguably the most powerful source of genetic novelty:
horizontal gene transfer. We review the status of rapid evolution and also
enumerate the current and future challenges of achieving a full understanding
of rapid evolution in all its manifestations.


  Animal movement networks are essential in understanding and containing the
spread of infectious diseases in farming industries. Due to its confidential
nature, movement data for the US swine farming population is not readily
available. Hence, we propose a method to generate such networks from limited
data available in the public domain. As a potentially devastating candidate, we
simulate the spread of African swine fever virus (ASFV) in our generated
network and analyze how the network structure affects the disease spread. We
find that high in-degree farm operations (i.e., markets) play critical roles in
the disease spread. We also find that high in-degree based targeted isolation
and hypothetical vaccinations are more effective for disease control compared
to other centrality-based mitigation strategies. The generated networks can be
made more robust by validation with more data whenever more movement data will
be available.


  To mitigate errors induced by the cell's heterogeneous noisy environment, its
main information channels and production networks utilize the kinetic
proofreading (KPR) mechanism. Here, we examine two extensively-studied KPR
circuits, DNA replication by the T7 DNA polymerase and translation by the E.
coli ribosome. Using experimental data, we analyze the performance of these two
vital systems in light of the fundamental bounds set by the recently-discovered
thermodynamic uncertainty relation (TUR), which places an inherent trade-off
between the precision of a desirable output and the amount of energy
dissipation required. We show that the DNA polymerase operates close to the TUR
lower bound, while the ribosome operates $\sim5$ times farther from this bound.
This difference originates from the enhanced binding discrimination of the
polymerase which allows it to operate effectively as a reduced reaction cycle
prioritizing correct product formation. We show that approaching this limit
also decouples the thermodynamic uncertainty factor from speed and error,
thereby relaxing the accuracy-speed trade-off of the system. Altogether, our
results show that operating near this reduced cycle limit not only minimizes
thermodynamic uncertainty, but also results in global performance enhancement
of KPR circuits.


  Conformational change of a DNA molecule is frequently observed in multiple
biological processes and has been modelled using a chain of strongly coupled
oscillators with a nonlinear bistable potential. While the mechanism and
properties of conformational change in the model have been investigated and
several reduced order models developed, the conformational dynamics as a
function of the length of the oscillator chain is relatively less clear. To
address this, we used a modified Lindstedt-Poincare method and numerical
computations. We calculate a perturbation expansion of the frequency of the
model's nonzero modes, finding that approximating these modes with their
unperturbed dynamics, as in a previous reduced order model, may not hold when
the length of the DNA model increases. We investigate the conformational change
to local perturbation in models of varying lengths, finding that for chosen
input and parameters, there are two regions of DNA length in the model, first
where the minimum energy required to undergo the conformational change
increases with DNA length; and second, where it is almost independent of the
length of the DNA model. We analyze the conformational change in these models
by adding randomness to the local perturbation, finding that the tendency of
the system to remain in a stable conformation against random perturbation
decreases with an increase in the DNA length. These results should help to
understand the role of the length of a DNA molecule in influencing its
conformational dynamics.


  We revisit the notion of gene regulatory code in embryonic development in the
light of recent findings about genome spatial organisation. By analogy with the
genetic code, we posit that the concept of code can only be used if the
corresponding adaptor can clearly be identified. An adaptor is here defined as
an intermediary physical entity mediating the correspondence between codewords
and objects in a gratuitous and evolvable way. In the context of the gene
regulatory code, the encoded objects are the gene expression levels, while the
concentrations of specific transcription factors in the cell nucleus provide
the codewords. The notion of code is meaningful in the absence of direct
physicochemical relationships between the objects and the codewords, when the
mediation by an adaptor is required. We propose that a plausible adaptor for
this code is the gene domain, that is, the genome segment delimited by
topological insulators and comprising the gene and its enhancer regulatory
sequences. We review recent evidences, based on genome-wide chromosome
conformation capture experiments, showing that preferential contact domains
found in metazoan genomes are the physical traces of gene domains. Accordingly,
genome 3D folding plays a direct role in shaping the developmental gene
regulatory code.


  The black and gold howler monkey (Alouatta caraya) is a neotropical primate
that faces the highest capture pressure for illegal trade in Argentina. We
evaluate the applicability of genetic assignment tests based on microsatellite
genotypic data to accurately assign individuals to their site of origin. The
search was conducted on a genetic database to determine the nearest sampled
population or to associate them to three clusters described here for the
Argentinean populations of A. caraya. We correctly assign 73% of the
individuals in the database to nearest population of origin, and 93.3% to their
cluster of origin. With this database, we were able to determine the probable
origin of 17 confiscated individuals, 12 of which were reintroduced in the
province of Misiones and 5 confiscated individuals reintroduced in the province
of Santa Fe. Moreover, we also determined the probable origin of 3 individuals
found dead in cities in northern Argentina. This approach highlights the
relevance of generating genotype indexing databases of species to assist with
in-situ and ex-situ conservation and management programs. Our results
underscore the importance of knowing the origin of individuals for
reintroduction and/or species recovery programs and to pinpoint the hotspots of
illegal capture of various species.


  Both action potentials and mechanosensitive signalling are an important
communication mechanisms in plants. Considering an information theoretic
framework, this paper explores the effective range of multiple action
potentials for a long chain of cells (i.e., up to 100) in different
configurations, and introduces the study of multiple mechanosensitive
activation signals (generated due to a mechanical stimulus) in plants. For both
these signals, we find that the mutual information per cell and information
propagation speed tends to increase up to a certain number of receiver cells.
However, as the number of cells increase beyond 10 to 12, the mutual
information per cell starts to decrease. To validate our model and results, we
include an experimental verification of the theoretical model, using a
PhytlSigns biosignal amplifier, allowing us to measure the magnitude of the
voltage associated with the multiple AP and mechanosensitive activation signals
induced by different stimulus in plants. Experimental data is used to calculate
the mutual information and information propagation speed, which is compared
with corresponding numerical results. Since these signals are used for a
variety of important tasks within the plant, understanding them may lead to new
bioengineering methods for plants.


  Data integration is the problem of combining multiple data groups (studies,
cohorts) and/or multiple data views (variables, features). This task is
becoming increasingly important in many disciplines due to the prevalence of
large and heterogeneous data sets. Data integration commonly aims to identify
structure that is consistent across multiple cohorts and feature sets. While
such joint analyses can boost information from single data sets, it is also
possible that a globally restrictive integration of heterogeneous data may
obscure signal of interest.
  Here, we therefore propose a data adaptive integration method, allowing for
structure in data to be shared across an a priori unknown \emph{subset of
cohorts and views}. The method, Multi-group Multi-view Principal Component
Analysis (MM-PCA), identifies partially shared, sparse low-rank components.
This also results in an integrative bi-clustering across cohorts and views. The
strengths of MM-PCA are illustrated on simulated data and on 'omics data from
The Cancer Genome Atlas. MM-PCA is available as an R-package.
  Key words: Data integration, Multi-view, Multi-group, Bi-clustering


  Next Generation Sequencing can sample the whole genome (WGS) or the 1-2% of
the genome that codes for proteins called the whole exome (WES). Machine
learning approaches to variant calling achieve high accuracy in WGS data, but
the reduced number of training examples causes training with WES data alone to
achieve lower accuracy. We propose and compare three different data
augmentation strategies for improving performance on WES data: 1) joint
training with WES and WGS data, 2) warmstarting the WES model from a WGS model,
and 3) joint training with the sequencing type specified. All three approaches
show improved accuracy over a model trained using just WES data, suggesting the
ability of models to generalize insights from the greater WGS data while
retaining performance on the specialized WES problem. These data augmentation
approaches may apply to other problem areas in genomics, where several
specialized models would each see only a subset of the genome.


  In cheminformatics, network representations of the space of compounds have
been suggested extensively. Among these, the threshold-network consists of
nodes representing molecules. In this network representation, two molecules are
connected by a link if the Tanimoto coefficient, a similarity measure, between
them exceeds a preset threshold. However, the topology of the threshold-network
is affected significantly by the preset threshold. In this study, we collected
the data of biologically relevant compounds and bioactivities. We defined the
weighted network where the weight of each link between the nodes equals the
Tanimoto coefficient between the bioactive compounds (nodes) without using the
threshold. We investigated the relationship between the strength of the link
connection and the bioactivity closeness in the weighted networks. We found
that compounds with significantly high or low bioactivity have a stronger
connection than those in the overall network.


  Reconstructing components of a genomic mixture from data obtained by means of
DNA sequencing is a challenging problem encountered in a variety of
applications including single individual haplotyping and studies of viral
communities. High-throughput DNA sequencing platforms oversample mixture
components to provide massive amounts of reads whose relative positions can be
determined by mapping the reads to a known reference genome; assembly of the
components, however, requires discovery of the reads' origin -- an NP-hard
problem that the existing methods struggle to solve with the required level of
accuracy. In this paper, we present a learning framework based on a graph
auto-encoder designed to exploit structural properties of sequencing data. The
algorithm is a neural network which essentially trains to ignore sequencing
errors and infers the posteriori probabilities of the origin of sequencing
reads. Mixture components are then reconstructed by finding consensus of the
reads determined to originate from the same genomic component. Results on
realistic synthetic as well as experimental data demonstrate that the proposed
framework reliably assembles haplotypes and reconstructs viral communities,
often significantly outperforming state-of-the-art techniques.


  Drug development is a very costly and lengthy process, while repositioned or
repurposed drugs could be brought into clinical practice within a shorter
time-frame and at a much reduced cost. The past decade has observed a massive
growth in the amount of data from genome-wide association studies (GWAS). The
rich information contained in GWAS data has great potential to guide drug
discovery or repositioning. Here we provide an overview of different
computational approaches which employ GWAS data to guide drug repositioning.
These methods include selection of top candidate genes from GWAS as drug
targets, deducing drug candidates based on drug-drug and disease-disease
similarity, searching for reversed expression profiles between drugs and
diseases, pathway-based methods as well as repositioning based on analysis of
biological networks. Each method is illustrated with examples, and their
respective strengths and limitations are discussed. Finally we discussed
several areas for future research.


  In this paper, we propose an original complex network model for an epidemic
problem in an heterogeneous geographical area. The complex network is
constructed by coupling nonidentical instances of a HIV/AIDS epidemiological
model for which a disease-free equilibrium and an endemic equilibrium can
coexist. After proving the existence of a positively invariant region for the
solutions of the complex network problem, we investigate the effect of the
coupling on the dynamics of the network, and establish the existence of a
unique disease-free equilibrium for the whole network, which is globally
asymptotically stable. We prove the existence of an optimal topology that
minimizes the level of infected individuals, and apply the theoretical results
to the case of the Cape Verde archipelago.


  Proteins are the major building blocks of life, and actuators of almost all
chemical and biophysical events in living organisms. Their native structures in
turn enable their biological functions which have a fundamental role in drug
design. This motivates predicting the structure of a protein from its sequence
of amino acids, a fundamental problem in computational biology. In this work,
we demonstrate state-of-the-art protein structure prediction (PSP) results
using embeddings and deep learning models for prediction of backbone atom
distance matrices and torsion angles. We recover 3D coordinates of backbone
atoms and reconstruct full atom protein by optimization. We create a new gold
standard dataset of proteins which is comprehensive and easy to use. Our
dataset consists of amino acid sequences, Q8 secondary structures, position
specific scoring matrices, multiple sequence alignment co-evolutionary
features, backbone atom distance matrices, torsion angles, and 3D coordinates.
We evaluate the quality of our structure prediction by RMSD on the latest
Critical Assessment of Techniques for Protein Structure Prediction (CASP) test
data and demonstrate competitive results with the winning teams and AlphaFold
in CASP13 and supersede the results of the winning teams in CASP12. We make our
data, models, and code publicly available.


  Objectives Lung squamous cell carcinoma (LUSC) often diagnosed as advanced
with poor prognosis. The mechanisms of its pathogenesis and prognosis require
urgent elucidation. This study was performed to screen potential biomarkers
related to the occurrence, development and prognosis of LUSC to reveal unknown
physiological and pathological processes. Materials and Methods Using
bioinformatics analysis, the lung squamous cell carcinoma microarray datasets
from the GEO and TCGA databases were analyzed to identify differentially
expressed genes(DEGs). Furthermore, PPI and WGCNA network analysis were
integrated to identify the key genes closely related to the process of LUSC
development. In addition, survival analysis was performed to achieve a
prognostic model that accomplished a high level of prediction accuracy. Results
and Conclusion Eighty-five up-regulated and 39 down-regulated genes were
identified, on which functional and pathway enrichment analysis was conducted.
GO analysis demonstrated that up-regulated genes were principally enriched in
epidermal development and DNA unwinding in DNA replication. Down-regulated
genes were mainly involved in cell adhesion, signal transduction and positive
regulation of inflammatory response. After PPI and WGCNA network analysis,
eight genes, including AURKA, RAD51, TTK, AURKB, CCNA2, TPX2, KPNA2 and KIF23,
have been found to play a vital role in LUSC development. The prognostic model
contained 20 genes, 18 of which were detrimental to prognosis. The AUC of the
established prognostic model for predicting the survival of patients at 1, 3,
and 5 years was 0.828, 0.826 and 0.824, respectively. To conclude, this study
identified a number of biomarkers of significant interest for additional
investigation of the therapies and methods of prognosis of lung squamous cell
carcinoma.


  Biomolecular feedback systems are now a central application area of interest
within control theory. While classical control techniques provide invaluable
insight into the function and design of both natural and synthetic biomolecular
systems, there are certain aspects of biological control that have proven
difficult to analyze with traditional methods. To this end, we describe here
how the recently developed tools of dominance analysis can be used to gain
insight into the nonlinear behavior of the antithetic integral feedback
circuit, a recently discovered control architecture which implements integral
control of arbitrary biomolecular processes using a simple feedback mechanism.
We show that dominance theory can predict both monostability and periodic
oscillations in the circuit, depending on the corresponding parameters and
architecture. We then use the theory to characterize the robustness of the
asymptotic behavior of the circuit in a nonlinear setting.


  The implementation of natural and artificial proteins with designer
properties and functionalities offers unparalleled opportunity for functional
nanoarchitectures formed through self-assembly. However, to exploit the
opportunities offered we require the ability to control protein assembly into
the desired architecture while avoiding denaturation and therefore retaining
protein functionality. Here we address this challenge with a model system of
fluorescent proteins. Using techniques of self-assembly manipulation inspired
by soft matter where interactions between components are controlled to yield
the desired structure, we show that it is possible to assemble networks of
proteins of one species which we can decorate with another, whose coverage we
can tune. Consequently, the interfaces between domains of each component can
also be tuned, with applications for example in energy transfer. Our model
system of fluorescent proteins eGFP and mCherry retain their fluorescence
throughout the assembly process, thus demonstrating that functionality is
preserved.


  Living cells sense their environment through the binding of extra-cellular
molecular ligands to cell surface receptors. Puzzlingly, vast numbers of
signaling pathways exhibit a high degree of cross talk between different
signals whereby different ligands act through the same receptor or shared
components downstream. It remains unclear how a cell can accurately process
information from the environment in such cross-wired pathways. We show that a
feature which commonly accompanies cross talk - signaling pleiotropy (the
ability of a receptor to produce multiple outputs) - offers a solution to the
cross talk problem. In a minimal model we show that a single pleiotropic
receptor can simultaneously identify and accurately sense the concentrations of
arbitrary unknown ligands present individually or in a mixture. We calculate
the fundamental limits of the signaling specificity and accuracy of such
signaling schemes. The model serves as an elementary "building block" towards
understanding more complex cross-wired receptor-ligand signaling networks.


  In this work we have investigated the evolutionary dynamics of a generalist
pathogen, e.g. a virus population, that evolves towards specialisation in an
environment with multiple host types. We have particularly explored under which
conditions generalist viral strains may rise in frequency and coexist with
specialist strains or even dominate the population. By means of a nonlinear
mathematical model and bifurcation analysis, we have determined the theoretical
conditions for stability of nine identified equilibria and provided biological
interpretation in terms of the infection rates for the viral specialist and
generalist strains. By means of a stability diagram we identified stable fixed
points and stable periodic orbits, as well as regions of bistability. For
arbitrary biologically feasible initial population sizes, the probability of
evolving towards stable solutions is obtained for each point of the analyzed
parameter space. This probability map shows combinations of infection rates of
the generalist and specialist strains that might lead to equal chances for each
type becoming the dominant strategy. Furthermore, we have identified infection
rates for which the model predicts the onset of chaotic dynamics. Several
degenerate Bogdanov-Takens and zero-Hopf bifurcations are detected along with
generalized Hopf and zero-Hopf bifurcations. This manuscript provides
additional insights into the dynamical complexity of host-pathogen evolution
towards different infection strategies.


  Is it feasible to create an analysis paradigm that can analyze and then
accurately and quickly predict known drugs from experimental data? PharML.Bind
is a machine learning toolkit which is able to accomplish this feat. Utilizing
deep neural networks and big data, PharML.Bind correlates
experimentally-derived drug affinities and protein-ligand X-ray structures to
create novel predictions. The utility of PharML.Bind is in its application as a
rapid, accurate, and robust prediction platform for discovery and personalized
medicine. This paper demonstrates that graph neural networks (GNNs) can be
trained to screen hundreds of thousands of compounds against thousands of
targets in minutes, a vastly shorter time than previous approaches. This
manuscript presents results from training and testing using the entirety of
BindingDB after cleaning; this includes a test set with 19,708 X-ray structures
and 247,633 drugs, leading to 2,708,151 unique protein-ligand pairings.
PharML.Bind achieves a prodigious 98.3% accuracy on this test set in under 25
minutes. PharML.Bind is premised on the following key principles: 1) speed and
a high enrichment factor per unit compute time, provided by high-quality
training data combined with a novel GNN architecture and use of
high-performance computing resources, 2) the ability to generalize to proteins
and drugs outside of the training set, including those with unknown active
sites, through the use of an active-site-agnostic GNN mapping, and 3) the
ability to be easily integrated as a component of increasingly-complex
prediction and analysis pipelines. PharML.Bind represents a timely and
practical approach to leverage the power of machine learning to efficiently
analyze and predict drug action on any practical scale and will provide utility
in a variety of discovery and medical applications.


  The evolution of drug-resistant microbial species is one of the major
challenges to global health. The development of new antimicrobial treatments
such as antimicrobial peptides needs to be accelerated to combat this threat.
However, the discovery of novel antimicrobial peptides is hampered by
low-throughput biochemical assays. Computational techniques can be used for
rapid screening of promising antimicrobial peptide candidates prior to testing
in the wet lab. The vast majority of existing antimicrobial peptide predictors
are non-targeted in nature, i.e., they can predict whether a given peptide
sequence is antimicrobial, but they are unable to predict whether the sequence
can target a particular microbial species. In this work, we have developed a
targeted antimicrobial peptide activity predictor that can predict whether a
peptide is effective against a given microbial species or not. This has been
made possible through zero-shot and few-shot machine learning. The proposed
predictor called AMP0 takes in the peptide amino acid sequence and any
N/C-termini modifications together with the genomic sequence of a target
microbial species to generate targeted predictions. It is important to note
that the proposed method can generate predictions for species that are not part
of its training set. The accuracy of predictions for novel test species can be
further improved by providing a few example peptides for that species. Our
computational cross-validation results show that the pro-posed scheme is
particularly effective for targeted antimicrobial prediction in comparison to
existing approaches and can be used for screening potential antimicrobial
peptides in a targeted manner especially for cases in which the number of
training examples is small. The webserver of the method is available at
http://ampzero.pythonanywhere.com.


  In this paper, we propose a novel approach for manifold learning that
combines the Earthmover's distance (EMD) with the diffusion maps method for
dimensionality reduction. We demonstrate the potential benefits of this
approach for learning shape spaces of proteins and other flexible
macromolecules using a simulated dataset of 3-D density maps that mimic the
non-uniform rotary motion of ATP synthase. Our results show that EMD-based
diffusion maps require far fewer samples to recover the intrinsic geometry than
the standard diffusion maps algorithm that is based on the Euclidean distance.
To reduce the computational burden of calculating the EMD for all volume pairs,
we employ a wavelet-based approximation to the EMD which reduces the
computation of the pairwise EMDs to a computation of pairwise weighted-$\ell_1$
distances between wavelet coefficient vectors.


  We present an edge-based framework for the study of geometric elastic network
models to model mechanical interactions in physical systems. We use a
formulation in the edge space, instead of the usual node-centric approach, to
characterise edge fluctuations of geometric networks defined in d- dimensional
space and define the edge mechanical embeddedness, an edge mechanical
susceptibility measuring the force felt on each edge given a force applied on
the whole system. We further show that this formulation can be directly related
to the infinitesimal rigidity of the network, which additionally permits three-
and four-centre forces to be included in the network description. We exemplify
the approach in protein systems, at both the residue and atomistic levels of
description.


  The evolution of antimicrobial resistance can be strongly affected by
variations of antimicrobial concentration. Here, we study the impact of
periodic alternations of absence and presence of antimicrobial on resistance
evolution in a microbial population, using a stochastic model that includes
variations of both population composition and size, and fully incorporates
stochastic population extinctions. We show that fast alternations of presence
and absence of antimicrobial are inefficient to eradicate the microbial
population and strongly favor the establishment of resistance, unless the
antimicrobial increases enough the death rate. We further demonstrate that if
the period of alternations is longer than a threshold value, the microbial
population goes extinct upon the first addition of antimicrobial, if it is not
rescued by resistance. We express the probability that the population is
eradicated upon the first addition of antimicrobial, assuming rare mutations.
Rescue by resistance can happen either if resistant mutants preexist, or if
they appear after antimicrobial is added to the environment. Importantly, the
latter case is fully prevented by perfect biostatic antimicrobials that
completely stop division of sensitive microorganisms. By contrast, we show that
the parameter regime where treatment is efficient is larger for biocidal drugs
than for biostatic drugs. This sheds light on the respective merits of
different antimicrobial modes of action.


  Development combines three basic processes asymmetric --- cell division,
signaling and gene regulation --- in a multitude of ways to create an
overwhelming diversity of multicellular life-forms. Here, we attempt to chart
this diversity using a generative model. We sample millions of biologically
feasible developmental schemes, allowing us to comment on the statistical
properties of cell-differentiation trajectories they produce. Our results
indicate that, in contrast to common views, cell-type lineage graphs are
unlikely to be tree-like. Instead, they are more likely to be directed acyclic
graphs, with multiple lineages converging on the same terminal cell-type.
Additionally, in line with the hypothesis that whole body regeneration is an
epiphenomenon of development, a majority of the `organism' generated by our
model can regenerate using pluripotent cells. The generative framework is
modular and flexible, and can be adapted to test additional hypotheses about
general features of development.


  One of the major goals of nanobionics and bottom-up synthetic biology is the
development of artificial cell organelles for the creation of cell-like
structures operating similar to biological systems with a minimalistic set of
building blocks. In the present contribution, versatile strategies to develop
artificial reaction centers for novel photoautotrophic processes and to provide
fully biocompatible synthetic enzyme counterparts for stepwise sustrate
photo-conversion in abiotic organelles are summarized. This unique approach has
the potential to create artificial cell-like systems that can be readily
manipulated by light providing both energy input and fundamental regulatory
functions.


  Hydrostatic pressure is a common perturbation to probe the conformations of
proteins. There are two common forms of pressure dependent potentials of mean
force (PMFs) derived from hydrophobic molecules available for the coarse
grained molecular simulations of protein folding and unfolding under
hydrostatic pressure. Although both PMF includes a desolvation barrier
separating the well of a direct contact and the well of a solvent mediated
contact, how these features vary with hydrostatic pressure is still debated.
There is a need of a systematic comparison of these two PMFs on a protein. We
investigated the two different pressure dependencies on the desolvation
potential in a structure based protein model using coarse grained molecular
simulations. We compared them to the known behavior a real protein based on
experimental evidence. We showed that the protein s folding transition curve on
the pressure temperature phase diagram depends on the relationship between the
potential well minima and pressure. For protein that reduces the total volume
under pressure, it is essential for the PMF to carry the feature that the
direct contact well is essential less stable than the water mediated contact
well at high pressure. We also comment on the practicality and importance of
structure based minimalist models for understanding the phenomenological
behavior of a protein under a wide range of phase space.


  In this paper, we construct a model to describe the transmission of HIV in a
homogeneous host population. By considering the specific mechanism of HIV, we
derive a model structured in three successive stages: (i) primary infection,
(ii) long phase of latency without symptoms and (iii) AIDS. Each HIV stage is
stratified by the duration for which individuals have been in the stage,
leading to a continuous age-structure model. In the first part of the paper, we
provide a global analysis of the model depending upon the basic reproduction
number R 0. When R 0 $\le$ 1, then the disease-free equilibrium is globally
asymptotically stable and the infection is cleared in the host population. On
the contrary, if R 0 > 1, we prove the epidemic's persistence with the
asymptotic stability of the endemic equilibrium. By performing the sensitivity
analysis, we then determine the impact 1 of control-related parameters of the
outbreak severity. For the second part, the initial model is extended with
intervention methods. By taking into account ART interventions and the
probability of treatment drop out, we discuss optimal interventions methods
which minimize the number of AIDS cases.


  We consider a collection of Markov chains that model the evolution of
multitype biological populations. The state space of the chains is the positive
orthant, and the boundary of the orthant is absorbing representing the
extinction states of different population types. We are interested in the
long-term behavior of the Markov chain away from extinction, under a small
noise scaling. Under this scaling, the trajectory of the Markov process over
any compact interval converges in distribution to the solution of an ordinary
differential equation (ODE) evolving in the positive orthant. We study the
asymptotic behavior of the quasi-stationary distributions (QSD) in this scaling
regime. Our main result shows that the limit points of the QSD are supported on
the union of interior attractors of the flow determined by the ODE. We also
give lower bounds on expected extinction times which scale exponentially with
the system size. Results of this type when the deterministic dynamical system
obtained under the scaling limit is given by a discrete time evolution equation
and the dynamics are essentially in a compact space (namely, the one step map
is a bounded function) have been studied by Faure and Schreiber (2014). Our
results extend these to a setting of an unbounded state space and continuous
time dynamics. The proofs rely on uniform large deviation results for small
noise stochastic dynamical systems and methods from the theory of dynamical
systems. In general QSD for Markov chains with absorbing states and unbounded
state spaces may not exist. We study one basic family of Binomial-Poisson
models in the positive orthant where one can use Lyapunov function methods to
establish existence of QSD and also to argue the tightness of the QSD of the
scaled sequence of Markov chains. The results from the first part are then used
to characterize the support of limit points of this sequence of QSD.


  The ability to switch on the activity of an enzyme through its spontaneous
reconstitution has proven to be a valuable tool in fundamental studies of
enzyme structure/reactivity relationships or in the design of artificial signal
transduction systems in bioelectronics, synthetic biology, or bioanalytical
applications. In particular, those based on the spontaneous
reconstitution/activation of the apo-PQQ-dependent soluble glucose
dehydrogenase (sGDH) from Acinetobacter calcoaceticus were widely developed.
However, the reconstitution mechanism of sGDH with its two cofactors, i.e. the
pyrroloquinoline quinone (PQQ) and Ca2+, remains unknown. The objective here is
to elucidate this mechanism by stopped-flow kinetics under single-turnover
conditions. The reconstitution of sGDH exhibited biphasic kinetics,
characteristic of a square reaction scheme associated to two activation
pathways. From a complete kinetic analysis, we were able to fully predict the
reconstitution dynamic, but also to demonstrate that when PQQ first binds to
the apo-sGDH, it strongly impedes the access of Ca2+ to its enclosed position
at the bottom of the enzyme binding site, thereby greatly slowing down the
reconstitution rate of sGDH. This slow calcium insertion may purposely be
accelerated by providing more flexibility to the Ca2+ binding loop through the
specific mutation of the calcium coordinating P248 proline residue, reducing
thus the kinetic barrier to calcium ion insertion. The dynamic nature of the
reconstitution process is also supported by the observation of a clear loop
shift and a reorganization of the hydrogen bonding network and van der Waals
interactions observed in both active sites of the apo and holo forms, a
structural change modulation that was revealed from the refined X-ray structure
of apo-sGDH (PDB:5MIN).


  Indoor microbial diversity and composition are suggested to affect the
prevalence and severity of asthma. In this study, we collected floor dust and
environmental characteristics from 21 classrooms, and health data related to
asthma symptoms from 309 students, in junior high schools in Johor Bahru,
Malaysia. Bacterial and fungal composition was characterized by sequencing 16s
rRNA gene and internal transcribed spacer (ITS) region, and the absolute
microbial concentration was quantified by qPCR. In total, 326 bacterial and 255
fungal genera were characterized. Five bacterial (Sphingobium, Rhodomicrobium,
Shimwellia, Solirubrobacter, Pleurocapsa) and two fungal (Torulaspora and
Leptosphaeriaceae) taxa were protective for asthma severity. Two bacterial
taxa, Izhakiella and Robinsoniella, were positively associated with asthma
severity. Several protective bacterial taxa including Rhodomicrobium,
Shimwellia and Sphingobium has been reported as protective microbes in previous
studies, whereas other taxa were first time reported. Environmental
characteristics, such as age of building, size of textile curtain per room
volume, occurrence of cockroaches, concentration of house dust mite allergens
transferred from homes by the occupants, were involved in shaping the overall
microbial community but not asthma-associated taxa; whereas visible dampness
and mold, which did not change the overall microbial community for floor dust,
decreased the concentration of protective bacteria Rhodomicrobium
(\b{eta}=-2.86, p=0.021) of asthma, indicating complex interactions between
microbes, environmental characteristics and asthma symptoms. Overall, this is
the first indoor microbiome study to characterize the asthma-associated
microbes and their environmental determinant in tropical area, promoting the
understanding of microbial exposure and respiratory health in this region.


  Bacterial infectious diseases are a major threat to human health. Timely and
sensitive pathogenic bacteria detection is crucial in identifying the bacterial
contaminations and preventing the spread of infectious diseases. Due to
limitations of conventional bacteria detection techniques there have been
concerted research efforts towards development of new biosensors. Biosensors
offering label free, whole bacteria detection are highly desirable over those
relying on label based or pathogenic molecular components detection. The major
advantage is eliminating the additional time and cost required for labeling or
extracting the desired bacterial components. Here, we demonstrate rapid,
sensitive and label free E. coli detection utilizing interferometric
reflectance imaging enhancement allowing for visualizing individual pathogens
captured on the surface. Enabled by our ability to count individual bacteria on
a large sensor surface, we demonstrate a limit of detection of 2.2 CFU/ml from
a buffer solution with no sample preparation. To the best of our knowledge,
this high level of sensitivity for whole E. coli detection is unprecedented in
label free biosensing. The specificity of our biosensor is validated by
comparing the response to target bacteria E. coli and non target bacteria S.
aureus, K. pneumonia and P. aeruginosa. The biosensor performance in tap water
also proves that its detection capability is unaffected by the sample
complexity. Furthermore, our sensor platform provides high optical
magnification imaging and thus validation of recorded detection events as the
target bacteria based on morphological characterization. Therefore, our
sensitive and label free detection method offers new perspectives for direct
bacterial detection in real matrices and clinical samples.


  Machine learning (ML), deep learning (DL), and artificial intelligence (AI)
are of increasing importance in biomedicine. The goal of this work is to show
progress in ML in digital health, to exemplify future needs and trends, and to
identify any essential prerequisites of AI and ML for precision health.
High-throughput technologies are delivering growing volumes of biomedical data,
such as large-scale genome-wide sequencing assays, libraries of medical images,
or drug perturbation screens of healthy, developing, and diseased tissue.
Multi-omics data in biomedicine is deep and complex, offering an opportunity
for data-driven insights and automated disease classification. Learning from
these data will open our understanding and definition of healthy baselines and
disease signatures. State-of-the-art applications of deep neural networks
include digital image recognition, single cell clustering, and virtual drug
screens, demonstrating breadths and power of ML in biomedicine. Significantly,
AI and systems biology have embraced big data challenges and may enable novel
biotechnology-derived therapies to facilitate the implementation of precision
medicine approaches.


  The search for meaningful structure in biological data has relied on
cutting-edge advances in computational technology and data science methods.
However, challenges arise as we push the limits of scale and complexity in
biological problems. Innovation in massively parallel, classical computing
hardware and algorithms continues to address many of these challenges, but
there is a need to simultaneously consider new paradigms to circumvent current
barriers to processing speed. Accordingly, we articulate a view towards quantum
computation and quantum information science, where algorithms have demonstrated
potential polynomial and exponential computational speedups in certain
applications, such as machine learning. The maturation of the field of quantum
computing, in hardware and algorithm development, also coincides with the
growth of several collaborative efforts to address questions across length and
time scales, and scientific disciplines. We use this coincidence to explore the
potential for quantum computing to aid in one such endeavor: the merging of
insights from genetics, genomics, neuroimaging and behavioral phenotyping. By
examining joint opportunities for computational innovation across fields, we
highlight the need for a common language between biological data analysis and
quantum computing. Ultimately, we consider current and future prospects for the
employment of quantum computing algorithms in the biological sciences.


  The canonical bacteriophage is obligately lytic: the virus infects a
bacterium and hijacks cell functions to produce large numbers of new viruses
which burst from the cell. These viruses are well-studied, but there exist a
wide range of coexisting virus lifestyles that are less understood. Temperate
viruses exhibit both a lytic cycle and a latent (lysogenic) cycle, in which
viral genomes are integrated into the bacterial host. Meanwhile, chronic
(persistent) viruses use cell functions to produce more viruses without killing
the cell; chronic viruses may also exhibit a latent stage in addition to the
productive stage. Here, we study the ecology of these competing viral
strategies. We demonstrate the conditions under which each strategy is
dominant, which aids in control of human bacterial infections using viruses. We
find that low lysogen frequencies provide competitive advantages for both virus
types; however, chronic viruses maximize steady state density by eliminating
lysogeny entirely, while temperate viruses exhibit a non-zero `sweet spot'
lysogen frequency. Viral steady state density maximization leads to coexistence
of temperate and chronic viruses, explaining the presence of multiple viral
strategies in natural environments.


  In this paper we provide a model to describe the dynamics of the species of
the ecosystem after it has been raided by a bad competing specie. The competing
specie invades the native plants for nutrition, carbon dioxide and space. This
affects the population of the native species of the ecosystem. The effect of
the bad biomass on the ecosystem is examined by considering it's equilibrium
points as well as its stability. An optimal control system is developed by
using Pontryagin's maximum principle to construct a Hamiltonian function which
minimizes the spread of the bad biomass. Numerical simulations were conducted
to analyze the results. It was established that the intrinsic growth rate r of
the good biomass was responsible for the sustenance and continuous survival of
the ecosystem. If growth rate of the good biomass is increased, the other
species showed positive growth. In addition to the growth rate of the good
biomass, death rates of both fish and birds also affected the state of the
ecosystem. Simulation results also revealed that the invasive bad species
biomass introduction affected the growth of the good biomass. This was as a
result of the competition for nutrition, carbon dioxide, space between the good
and the bad biomases. After implementing the control units, simulation results
showed an improvement in the growth of the good biomass, fish and birds
populations whereas it showed a decline in the bad biomass growth.


  Determining the primary site of origin for metastatic tumors is one of the
open problems in cancer care because the efficacy of treatment often depends on
the cancer tissue of origin. Classification methods that can leverage tumor
genomic data and predict the site of origin are therefore of great value.
Because tumor DNA point mutation data is very sparse, only limited accuracy
(64.5% for 12 tumor classes) was previously demonstrated by methods that rely
on point mutations as features (1). Tumor classification accuracy can be
greatly improved (to over 90% for 33 classes) by relying on gene expression
data (2). However, this additional data is often not readily available in
clinical setting, because point mutations are better profiled and targeted by
clinical mutational profiling.
  Here we sought to develop an accurate deep transfer learning and fine-tuning
method for tumor sub-type classification, where predicted class is indicative
of the primary site of origin. Our method significantly outperforms the
state-of-the-art for tumor classification using DNA point mutations, reducing
the error by more than 30% at the same time discriminating over many more
classes on The Cancer Genome Atlas (TCGA) dataset. Using our method, we achieve
state-of-the-art tumor type classification accuracy of 78.3% for 29 tumor
classes relying on DNA point mutations in the tumor only.


  Methicillin-resistant Staphylococcus aureus (MRSA) is a significant
contributor to the growing concern of antibiotic resistant bacteria, especially
given its stubborn persistence in hospitals and other health care facility
settings. In combination with this characteristic of S. aureus (colloquially
referred to as staph), MRSA presents an additional barrier to treatment and is
now believed to have colonized two of every 100 people worldwide. According to
the CDC, MRSA prevalence sits as high as 25-50% in countries such as the United
Kingdom and the United States. Given the resistant nature of staph as well as
its capability of evolving to compensate antibiotic treatment, controlling MRSA
levels is more a matter of precautionary and defensive measures. This study
examines the method of "search and isolation," which seeks to isolate MRSA
positive patients in a hospital so as to decrease infection potential. Although
this strategy is straightforward, the question of just whom to screen is of
practical importance. We compare screening at admission to screening at
discharge. To do this, we develop a mathematical model and use simulations to
determine MRSA endemic levels in a hospital with either control measure
implemented. We found that screening at discharge was the more effective method
in controlling MRSA endemicity, but at the cost of a greater number of isolated
patients.


  In many biological systems, the movement of individual agents is commonly
characterized as having multiple qualitatively distinct behaviors that arise
from various biophysical states. This is true for vesicles in intracellular
transport, micro-organisms like bacteria, or animals moving within and
responding to their environment. For example, in cells the movement of
vesicles, organelles and other cargo are affected by their binding to and
unbinding from cytoskeletal filaments such as microtubules through molecular
motor proteins. A typical goal of theoretical or numerical analysis of models
of such systems is to investigate the effective transport properties and their
dependence on model parameters. While the effective velocity of particles
undergoing switching diffusion is often easily characterized in terms of the
long-time fraction of time that particles spend in each state, the calculation
of the effective diffusivity is more complicated because it cannot be expressed
simply in terms of a statistical average of the particle transport state at one
moment of time. However, it is common that these systems are regenerative, in
the sense that they can be decomposed into independent cycles marked by returns
to a base state. Using decompositions of this kind, we calculate effective
transport properties by computing the moments of the dynamics within each cycle
and then applying renewal-reward theory. This method provides a useful
alternative large-time analysis to direct homogenization for linear
advection-reaction-diffusion partial differential equation models. Moreover, it
applies to a general class of semi-Markov processes and certain stochastic
differential equations that arise in models of intracellular transport.
Applications of the proposed framework are illustrated for case studies such as
mRNA transport in developing oocytes and processive cargo movement by teams of
motor proteins.


  We present a novel approach to represent ecological systems using reaction
networks, and show how a particular framework called Chemical Organization
Theory (COT) sheds new light on the longstanding complexity-stability debate.
Namely, COT provides a novel conceptual landscape plenty of analytic tools to
explore the interplay between structure and stability of ecological systems.
Given a large set of species and their interactions, COT identifies, in a
computationally feasible way, each and every sub-collection of species that is
closed and self-maintaining. These sub-collections, called organizations,
correspond to the groups of species that can survive together (co-exist) in the
long-term. Thus, the set of organizations contains all the stable regimes that
can possibly happen in the dynamics of the ecological system. From here, we
propose to conceive the notion of stability from the properties of the
organizations, and thus apply the vast knowledge on the stability of reaction
networks to the Complexity-Stability debate. As an example of the potential of
COT to introduce new mathematical tools, we show that the set of organizations
can be equipped with suitable joint and meet operators, and that for certain
ecological systems the organizational structure is a non-boolean lattice,
providing in this way an unexpected connection between logico-algebraic
structures, popular in the foundations of quantum theory, and ecology.


  Abiotic emergence of ordered information stored in the form of RNA is an
important unresolved problem concerning the origin of life. A polymer longer
than 40--100 nucleotides is necessary to expect a self-replicating activity,
but the formation of such a long polymer having a correct nucleotide sequence
by random reactions seems statistically unlikely. However, our universe,
created by a single inflation event, likely includes more than $10^{100}$
Sun-like stars. If life can emerge at least once in such a large volume, it is
not in contradiction with our observations of life on Earth, even if the
expected number of abiogenesis events is negligibly small within the observable
universe that contains only $10^{22}$ stars. Here, a quantitative relation is
derived between the minimum RNA length $l_{\min}$ required to be the first
biological polymer, and the universe size necessary to expect the formation of
such a long and active RNA by randomly adding monomers. It is then shown that
an active RNA can indeed be produced somewhere in an inflationary universe,
giving a solution to the abiotic polymerization problem. On the other hand,
$l_{\min}$ must be shorter than $\sim$20 nucleotides for the abiogenesis
probability close to unity on a terrestrial planet, but a self-replicating
activity is not expected for such a short RNA. Therefore, if extraterrestrial
organisms of a different origin from those on Earth are discovered in the
future, it would imply an unknown mechanism at work to polymerize nucleotides
much faster than random statistical processes.


  The second domain of gelsolin (G2) hosts mutations responsible for a
hereditary form of amyloidosis. The active form of gelsolin is Ca2+-bound; it
is also a dynamic protein, hence structural biologists often rely on the study
of the isolated G2. However, the wild type G2 structure that have been used so
far in comparative studies is bound to a crystallographic Cd2+, in lieu of the
physiological calcium. Here, we report the wild type structure of G2 in complex
with Ca2+ highlighting subtle ion-dependent differences. Previous findings on
different G2 mutations are also briefly revised in light of these results.


  Mutations in the gelsolin protein are responsible for a rare conformational
disease known as AGel amyloidosis. Four of these mutations are hosted by the
second domain of the protein (G2): D187N/Y, G167R and N184K. The impact of the
latter has been so far evaluated only by studies on the isolated G2. Here we
report the characterization of full-length gelsolin carrying the N184K mutation
and compare the findings with those obtained on the wild type and the other
variants. The crystallographic structure of the N184K variant in the Ca2+-free
conformation shows remarkable similarities with the wild type protein. Only
minimal local rearrangements can be observed and the mutant is as efficient as
the wild type in severing filamentous actin. However, thermal stability of the
pathological variant is compromised in the Ca2+-free conditions. These data
suggest that the N to K substitution causes a local disruption of the H-bond
network in the core of the G2 domain. Such a subtle rearrangement of the
connections does not lead to significant conformational changes but severely
affects the stability of the protein.


  We analyze the long term behavior of interacting populations which can be
controlled through harvesting. The dynamics is assumed to be discrete in time
and stochastic due to the effect of environmental fluctuations. We present
extinction and coexistence criteria when there are one or two interacting
species. We then use these tools in order to see when harvesting leads to
extinction or persistence of species, as well as what the optimal harvesting
strategies, which maximize the expected long term yield, look like. For single
species systems, we show under certain conditions that the optimal harvesting
strategy is of bang-bang type: there is a threshold under which there is no
harvesting, while everything above this threshold gets harvested.
  The second part of the paper is concerned with the analysis of ecosystems
that have two interacting species which can be harvested. In particular, we
carefully study predator-prey and competitive Ricker models when there are two
species. In this setting we show how to find the optimal proportional
harvesting strategy. If the system is of predator-prey type the optimal
proportional harvesting strategy is, depending on the interaction parameters
and the price of predators relative to prey, either to harvest the predator to
extinction and maximize the asymptotic yield of the prey or to not harvest the
prey and to maximize the asymptotic harvesting yield of the predators. If the
system is competitive, in certain instances it is optimal to drive one species
extinct and to harvest the other one. In other cases it is best to let the two
species coexist and harvest both species while maintaining coexistence.
  In the setting of the competitive Ricker model we show that if one competitor
is dominant and pushes the other species to extinction, the harvesting of the
dominant species can lead to coexistence.


  PURPOSE: The popularity of germline genetic panel testing has led to a vast
accumulation of variant-level data. Variant names are not always consistent
across laboratories and not easily mappable to public variant databases such as
ClinVar. A tool that can automate the process of variants harmonization and
mapping is needed to help clinicians ensure their variant interpretations are
accurate. METHODS: We present a Python-based tool, Ask2Me VarHarmonizer, that
incorporates data cleaning, name harmonization, and a four-attempt mapping to
ClinVar procedure. We applied this tool to map variants from a pilot dataset
collected from 11 clinical practices. Mapping results were evaluated with and
without the transcript information. RESULTS: Using Ask2Me VarHarmonizer, 4728
out of 6027 variant entries (78%) were successfully mapped to ClinVar,
corresponding to 3699 mappable unique variants. With the addition of 1099
unique unmappable variants, a total of 4798 unique variants were eventually
identified. 427 (9%) of these had multiple names, of which 343 (7%) had
multiple names within-practice. 99% mapping consistency was observed with and
without transcript information. CONCLUSION: Ask2Me VarHarmonizer aggregates and
structures variant data, harmonizes names, and maps variants to ClinVar.
Performing harmonization removes the ambiguity and redundancy of variants from
different sources.


  Many native structures of proteins accomodate complex topological motifs such
as knots, lassos, and other geometrical entanglements. How proteins can fold
quickly even in the presence of such topological obstacles is a debated
question in structural biology. Recently, the hypothesis that energetic
frustration might be a mechanism to avoid topological frustration has been put
forward based on the empirical observation that loops involved in entanglements
are stabilized by weak interactions between amino-acids at their extrema. To
verify this idea, we use a toy lattice model for the folding of proteins into
two almost identical structures, one entangled and one not. As expected, the
folding time is longer when random sequences folds into the entangled
structure. This holds also under an evolutionary pressure simulated by
optimizing the folding time. It turns out that optmized protein sequences in
the entangled structure are in fact characterized by frustrated interactions at
the closures of entangled loops. This phenomenon is much less enhanced in the
control case where the entanglement is not present. Our findings, which are in
agreement with experimental observations, corroborate the idea that an
evolutionary pressure shapes the folding funnel to avoid topological and
kinetic traps.


  We investigate the complexity of local search based on steepest ascent. We
show that even when all variables have domains of size two and the underlying
constraint graph of variable interactions has bounded treewidth (in our
construction, treewidth 7), there are fitness landscapes for which an
exponential number of steps may be required to reach a local optimum. This is
an improvement on prior recursive constructions of long steepest ascents, which
we prove to need constraint graphs of unbounded treewidth.


  Traditional approaches to elucidation of protein structures by NMR
spectroscopy rely on distance restraints also known as nuclear Overhauser
effects (NOEs). The use of NOEs as the primary source of structure
determination by NMR spectroscopy is time consuming and expensive. Residual
Dipolar Couplings (RDCs) have become an alternate approach for structure
calculation by NMR spectroscopy. In previous works, the software package
REDCRAFT has been presented as a means of harnessing the information containing
in RDCs for structure calculation of proteins. In this work, we present
significant improvements to the REDCRAFT package including: refinement of the
decimation procedure, the inclusion of graphical user interface, adoption of
NEF standards, and addition of scripts for enhanced protein modeling options.
The improvements to REDCRAFT have resulted in the ability to fold proteins that
the previous versions were unable to fold. For instance, we report the results
of folding of the protein 1A1Z in the presence of highly erroneous data.


  Large scale initiatives such as the Human Genome Project, Structural
Genomics, and individual research teams have provided large deposits of genomic
and proteomic data. The transfer of data to knowledge has become one of the
existing challenges, which is a consequence of capturing data in databases that
are optimally designed for archiving and not mining. In this research, we have
targeted the Protein Databank (PDB) and demonstrated a transformation of its
content, named PDBMine, that reduces storage space by an order of magnitude,
and allows for powerful mining in relation to the topic of protein structure
determination. We have demonstrated the utility of PDBMine in exploring the
prevalence of dimeric and trimeric amino acid sequences and provided a
mechanism of predicting protein structure.


  Gene set enrichment analyses of 8,405 genes linked with 35,074 human-specific
(hs) regulatory single-nucleotide changes (SNCs) revealed the staggering
breadth of significant associations with morphological structures,
physiological processes, and pathological conditions of Modern Humans.
Significant enrichment traits include more than 1,000 anatomically-distinct
regions of the adult human brain, many different types of human cells and
tissues, more than 200 common human disorders and more than 1,000 records of
rare diseases. Thousands of genes connected with regulatory hsSNCs have been
identified in this contribution, which represent essential genetic elements of
the autosomal inheritance and survival of species phenotypes: a total of 1,494
genes linked with either autosomal dominant or recessive inheritance as well as
2,273 genes associated with premature death, embryonic lethality, as well as
pre-, peri-, neo-, and post-natal lethality of both complete and incomplete
penetrance. Therefore, thousands of heritable traits and critical genes
impacting the offspring survival appear under the human-specific regulatory
control in genomes of Modern Humans. These observations highlight the
remarkable translational opportunities afforded by the discovery of genetic
regulatory loci harboring hsSNCs that are fixed in humans, distinct from other
primates, and located in differentially-accessible (DA) chromatin regions
during human brain development.


  Protein domains are highly conserved functional units of proteins. Because
they carry functionally significant information, the majority of the coding
disease variants are located on domains. Additionally, domains are specific
units of the proteins that can be targeted for drug delivery purposes. Here,
using information about variants sites associated with diseases, a disease
network was built, based on their sharing the same domain and domain variation
site. The result was 49,990 disease pairs linked by domain variant site and
533,687 disease pairs that share the same mutated domain. These pairs were
compared to disease pairs made using previous methods such as gene identity and
gene variant site identity, which revealed that over 8,000 of these pairs were
not only missing from the gene pairings but also not found commonly together in
literature. The disease network was analyzed from their disease subject
categories, which when compared to the gene-based disease network revealed that
the domain method results in higher number of connections across disease
categories versus within a disease category. Further, a study into the drug
repurposing possibilities of the disease network created using domain revealed
that 16,902 of the disease pairs had a drug reported for one disease but not
the other, highlighting the drug repurposing potential of this new methodology.


  Over 50% of the world population is at risk of mosquito-borne diseases.
Female Ae. aegypti mosquito species transmit Zika, Dengue, and Chikungunya. The
spread of these diseases correlate positively with the vector population, and
this population depends on biotic and abiotic environmental factors including
temperature, vegetation condition, humidity and precipitation. To combat virus
outbreaks, information about vector population is required. To this aim, Earth
observation (EO) data provide fast, efficient and economically viable means to
estimate environmental features of interest. In this work, we present a
temporal distribution model for adult female Ae. aegypti mosquitoes based on
the joint use of the Normalized Difference Vegetation Index, the Normalized
Difference Water Index, the Land Surface Temperature (both at day and night
time), along with the precipitation information, extracted from EO data. The
model was applied separately to data obtained during three different vector
control and field data collection condition regimes, and used to explain the
differences in environmental variable contributions across these regimes. To
this aim, a random forest (RF) regression technique and its nonlinear features
importance ranking based on mean decrease impurity (MDI) were implemented. To
prove the robustness of the proposed model, we compared it to other machine
learning techniques and statistical models. Our results show that machine
learning techniques perform better than linear statistical models for the task
at hand, and RF performs best. By ranking the importance of all features based
on MDI in RF and selecting the subset comprising the most informative ones, a
more parsimonious but equally effective and explainable model can be
obtained...


  Genome-wide association studies (GWAS) provide a means of examining the
common genetic variation underlying a range of traits and disorders. In
addition, it is hoped that GWAS may provide a means of differentiating affected
from unaffected individuals. This has potential applications in the area of
risk prediction. Current attempts to address this problem focus on using the
polygene risk score (PRS) to predict case-control status on the basis of GWAS
data. However this approach has so far had limited success for complex traits
such as schizophrenia (SZ). This is essentially a classification problem.
Artificial neural networks (ANNs) have been shown in recent years to be highly
effective in such applications. Here we apply an ANN to the problem of
distinguishing SZ patients from unaffected controls. We compare the
effectiveness of the ANN with the PRS in classifying individuals by
case-control status based only on genetic data from a GWAS. We use the
schizophrenia dataset from the Psychiatric Genomics Consortium (PGC) for this
study. Our analysis indicates that the ANN is more sensitive to sample size
than the PRS. As larger and larger sample sizes become available, we suggest
that ANNs are a promising alternative to the PRS for classification and risk
prediction for complex genetic disorders.


  Analysis of somatic mutation profiles from cancer patients is essential in
the development of cancer research. However, the low frequency of most
mutations and the varying rates of mutations across patients makes the data
extremely challenging to statistically analyze as well as difficult to use in
classification problems, for clustering, visualization or for learning useful
information. Thus, the creation of low dimensional representations of somatic
mutation profiles that hold useful information about the DNA of cancer cells
will facilitate the use of such data in applications that will progress
precision medicine. In this paper, we talk about the open problem of learning
from somatic mutations, and present Flatsomatic: a solution that utilizes
variational autoencoders (VAEs) to create latent representations of somatic
profiles. The work done in this paper shows great potential for this method,
with the VAE embeddings performing better than PCA for a clustering task, and
performing equally well to the raw high dimensional data for a classification
task. We believe the methods presented herein can be of great value in future
research and in bringing data-driven models into precision oncology.


  Often the development of novel functional peptides is not amenable to high
throughput or purely computational screening methods. Peptides must be
synthesized one at a time in a process that does not generate large amounts of
data. One way this method can be improved is by ensuring that each experiment
provides the best improvement in both peptide properties and predictive
modeling accuracy. Here, we study the effectiveness of active learning,
optimizing experiment order, and meta-learning, transferring knowledge between
contexts, to reduce the number of experiments necessary to build a predictive
model. We present a multi-task benchmark database of peptides designed to
advance these methods for experimental design. Each task is binary
classification of peptides represented as a sequence string. We find neither
active learning method tested to be better than random choice. The
meta-learning method Reptile was found to improve average accuracy across
datasets. Combining meta-learning with active learning offers inconsistent
benefits.


  In this work, we employ autoregressive models developed in financial
engineering for modeling of forest dynamics. Autoregressive models have some
theoretical advantage over currently employed forest modeling approaches such
as Markov chains and individual-based models, as autoregressive models are both
analytically tractable and operate with continuous state space. We perform time
series statistical analysis of forest biomass and basal area recorded in Quebec
provincial forest inventories in 1970-2007. The geometric random walk model
adequately describes the yearly average dynamics. For individual patches, we
fit an AR(1) process capable to model negative feedback (mean-reversion).
Overall, the best fit also turns out to be geometric random walk, however, the
normality tests for residuals fail. In contrast, yearly means are adequately
described by normal fluctuations, with annual growth, on average, 2.3%, but
with standard deviation of order 40%. We use Bayesian analysis to account for
uneven number of observations per year. This work demonstrates that
autoregressive models represent a valuable tool for modeling of forest
dynamics. In particular, they quantify stochastic effects of environmental
disturbances and develop predictive empirical models on short and intermediate
temporal scales.


  Cyclic dominance between species may yield spiral waves that are known to
provide a mechanism enabling persistent species coexistence. This observation
holds true even in presence of spatial heterogeneity in the form of quenched
disorder. In this work we study the effects on spatio-temporal patterns and
species coexistence of structured spatial heterogeneity in the form of habitats
that locally provide one of the species with an advantage. Performing extensive
numerical simulations of systems with three and six species we show that these
structured habitats destabilize spiral waves. Analyzing extinction events, we
find that species extinction probabilities display a succession of maxima as
function of time, that indicate a periodically enhanced probability for species
extinction. Analysis of the mean extinction time reveals that as a function of
the parameter governing the advantage of one of the species a transition
between stable coexistence and unstable coexistence takes place. We also
investigate how efficiency as a predator or a prey affects species coexistence.


  The adhesion of biomembranes is mediated by the binding of membrane-anchored
receptor and ligand proteins. The proteins can only bind if the separation
between apposing membranes is sufficiently close to the length of the protein
complexes, which leads to an interplay between protein binding and membrane
shape. In this article, we review current models of biomembrane adhesion and
novel insights obtained from the models. Theory and simulations with
elastic-membrane and coarse-grained molecular models of biomembrane adhesion
indicate that the binding of proteins in membrane adhesion strongly depends on
nanoscale shape fluctuations of the apposing membranes, which results in
binding cooperativity. A length mismatch between protein complexes leads to
repulsive interactions that are caused by membrane bending and act as a driving
force for the length-based segregation of proteins during membrane adhesion.


  In animal groups, individual decisions are best characterised by
probabilistic rules. Furthermore, animals of many species live in small groups.
Probabilistic interactions among small numbers of individuals lead to a so
called intrinsic noise at the group level. Theory predicts that the strength of
intrinsic noise is not a constant but often depends on the collective state of
the group; hence, it is also called a state-dependent noise or a multiplicative
noise. Surprisingly, such noise may produce collective order. However, only a
few empirical studies on collective behaviour have paid attention to such
effects due to the lack of methods that enable us to connect data with theory.
Here, we demonstrate a method to characterise the role of stochasticity
directly from high-resolution time-series data of collective dynamics. We do
this by employing two well-studied individual-based toy models of collective
behaviour. We argue that the group-level noise may encode important information
about the underlying processes at the individual scale. In summary, we describe
a method that enables us to establish connections between empirical data of
animal (or cellular) collectives with the phenomenon of noise-induced states, a
field that is otherwise largely limited to the theoretical literature.


  Many different processes take place to facilitate lubrication of the joints
functioning in human locomotion system. The main purpose of this is to avoid
destroying the articular cartilage. Viscoelastic properties of the joints
system are very sensitive on both temperature and concentration changes because
of the change in conformation presented in the system proteins and protein
network formation. We are searching for an answer to the question on how
changes in temperature and concentration influence the conformational entropy
of mucin protein which is a part of one of the key components, lubricin, which
is believed to be responsible for gel formation inside synovial fluid. We are
using molecular dynamic technique to obtain the information about dihedral
(phi, psi) angles of the mucin during protein self - assembly by means of the
computer simulation with a time duration up to 50 ns, parameterized by six
temperatures ranged between 300 - 315 K, and six concentrations 10.68 - 267.1
g/L. The results show that between c3 and c4 (160 g/L and 214 g/L) a transition
exists where crowding begins affecting the dynamics of protein network
formation. In such a concentration ranges mucin has a chance to change the
frictional properties of the system. Simultaneously there were no significant
changes in conformations of the mucins molecules even after they created
networks. The temperature changes also did not affect much of mucins
conformations but it introduced slightly modifications in dihedral angles and
after some critical value T=306 K it changed conformational entropy trend from
decreasing to raising.


  Many aspects of the study of protein folding and dynamics have been affected
by the recent advances in machine learning. Methods for the prediction of
protein structures from their sequences are now heavily based on machine
learning tools. The way simulations are performed to explore the energy
landscape of protein systems is also changing as force-fields are started to be
designed by means of machine learning methods. These methods are also used to
extract the essential information from large simulation datasets and to enhance
the sampling of rare events such as folding/unfolding transitions. While
significant challenges still need to be tackled, we expect these methods to
play an important role on the study of protein folding and dynamics in the near
future. We discuss here the recent advances on all these fronts and the
questions that need to be addressed for machine learning approaches to become
mainstream in protein simulation.


  Fatty acid sugar esters represent an important class of non-ionic bio-based
surfactants. They can be synthesized from vinyl fatty acids and sugars with
enzyme as a catalyst. Herein, the influence of the solvent, the lipase and the
temperature on a model reaction between vinyl palmitate and glucose via
enzymatic catalysis has been investigated and the reaction conditions
optimized. Full conversion into 6-O-glucose palmitate was reached in 40 hours
in acetonitrile starting from a reactant ratio 1:1, at only 5%-wt loading of
lipase from Candida antarctica B (CALB) without the presence of molecular
sieves.


  PPI network alignment aims to find topological and functional similarities
between networks of different species. Several alignment approaches have been
proposed. Each of these approaches relies on a different alignment method and
uses different biological information during the alignment process such as the
topological structure of the networks and the sequence similarities between the
proteins, but less of them integrate the functional similarities between
proteins. In this context, we present our algorithm PPINA (Protein-Protein
Interaction Network Aligner), which is an extension of the NETAL algorithm. The
latter aligns two networks based on the sequence, functional and network
topology similarity of the proteins. PPINA has been tested on real PPI
networks. The results show that PPINA has outperformed other alignment
algorithms where it provides biologically meaningful results.