  Cryo-electron microscopy (cryo-EM) extracts single-particle density
projections of individual biomolecules. Although cryo-EM is widely used for 3D
reconstruction, due to its single-particle nature, it has the potential to
provide information about the biomolecule's conformational variability and
underlying free energy landscape. However, treating cryo-EM as a
single-molecule technique is challenging because of the low signal-to-noise
ratio (SNR) in the individual particles. In this work, we developed the
cryo-BIFE method, cryo-EM Bayesian Inference of Free Energy profiles, that uses
a path collective variable to extract free energy profiles and their
uncertainties from cryo-EM images. We tested the framework over several
synthetic systems, where we controlled the imaging parameters and conditions.
We found that for realistic cryo-EM environments and relevant biomolecular
systems, it is possible to recover the underlying free energy, with the pose
accuracy and SNR as crucial determinants. Then, we used the method to study the
conformational transitions of a calcium-activated channel with real cryo-EM
particles. Interestingly, we recover the most probable conformation (used to
generate a high resolution reconstruction of the calcium-bound state), and we
find two additional meta-stable states, one which corresponds to the
calcium-unbound conformation. As expected for turnover transitions within the
same sample, the activation barriers are of the order of a couple $k_BT$.
Extracting free energy profiles from cryo-EM will enable a more complete
characterization of the thermodynamic ensemble of biomolecules.


  Mutualistic networks have attracted increasing attention in the ecological
literature in the last decades as they play a key role in the maintenance of
biodiversity. Here, we develop an analytical framework to study the structural
stability of these networks including both mutualistic and competitive
interactions. Analytical and numerical analyses show that the structure of the
competitive network fundamentally alters the necessary conditions for species
coexistence in communities. Using 50 real mutualistic networks, we show that
when the relative importance of shared partners is incorporated via weighted
competition, the feasibility area in the parameter space is highly correlated
with May's stability criteria and can be predicted by a functional relationship
between the number of species, the network connectance and the average
interaction strength in the community. Our work reopens a decade-long debate
about the complexity-stability relationship in ecological communities, and
highlights the role of the relative structures of different interaction types.


  Database-search algorithms, that deduce peptides from Mass Spectrometry (MS)
data, have tried to improve the computational efficiency to accomplish larger,
and more complex systems biology studies. Existing serial, and high-performance
computing (HPC) search engines, otherwise highly successful, are known to
exhibit poor-scalability with increasing size of theoretical search-space
needed for increased complexity of modern non-model, multi-species MS-based
omics analysis. Consequently, the bottleneck for computational techniques is
the communication costs of moving the data between hierarchy of memory, or
processing units, and not the arithmetic operations. This post-Moore change in
architecture, and demands of modern systems biology experiments have dampened
the overall effectiveness of the existing HPC workflows. We present a novel
efficient parallel computational method, and its implementation on
memory-distributed architectures for peptide identification tool called HiCOPS,
that enables more than 100-fold improvement in speed over most existing HPC
proteome database search tools. HiCOPS empowers the supercomputing database
search concept for comprehensive identification of peptides, and all their
modified forms within a reasonable time-frame. We demonstrate this by searching
Gigabytes of experimental MS data against Terabytes of databases where HiCOPS
completes peptide identification in few minutes using 72 parallel nodes (1728
cores) compared to several weeks required by existing state-of-the-art tools
using 1 node (24 cores); 100 minutes vs 5 weeks; 500x speedup. Finally, we
formulate a theoretical framework for our overhead-avoiding strategy, and
report superior performance evaluation results for key metrics including
execution time, CPU utilization, speedups, and I/O efficiency. The software
will be made available at: hicops.github.io


  We propose an epidemiological model that includes the mobility patterns of
the individuals, in the spirit to those considered in (Barmak, 2011, 2016) and
(Medus, 2011). We assume that people move around in a city of 120x120 blocks
with 300 inhabitants in each block. The mobility pattern is associated to a
complex network in which nodes represent blocks while the links represent the
traveling path of the individuals. We implemented three confinement strategies
in order to mitigate the disease spreading: 1) global confinement, 2) partial
restriction to mobility, and 3) localized confinement. In the first case, it
was observed that a global isolation policy prevents the massive outbreak of
the disease. In the second case, a partial restriction to mobility could lead
to a massive contagion if this was not complemented with sanitary measures such
as the use of masks and social distancing. Finally, a local isolation policy
was proposed, conditioned to the health status of each block. It was observed
that this mitigation strategy was able to contain and even reduce the outbreak
of the disease by intervening in specific regions of the city according to
their level of contagion. It was also observed that this strategy is capable of
controlling the epidemic in the case that a certain proportion of those
infected are asymptomatic.


  An origin of Earth life on Mars would resolve significant inconsistencies
between the inferred history of life and Earth's geologic history. Life as we
know it utilizes amino acids, nucleic acids, and lipids for the metabolic,
informational, and compartment-forming subsystems of a cell. Such building
blocks may have formed simultaneously from cyanosulfidic chemical precursors in
a planetary surface scenario involving ultraviolet light, wet-dry cycling, and
volcanism. However, early Earth was a water world, and the timing of the rise
of oxygen on Earth is inconsistent with final fixation of the genetic code in
response to oxidative stress. A cyanosulfidic origin of life could have taken
place on Mars via photoredox chemistry, facilitated by orders of magnitude more
sub-aerial crust than early Earth, and an earlier transition to oxidative
conditions. Meteoritic bombardment may have generated transient habitable
environments and ejected and transferred life to Earth. The Mars 2020
Perseverance Rover offers an unprecedented opportunity to confirm or refute
evidence consistent with a cyanosulfidic origin of life on Mars, search for
evidence of ancient life, and constrain the evolution of Mars' oxidation state
over time. We should seek to prove or refute a Martian origin for life on Earth
alongside other possibilities.


  SARS-CoV-2 started propagating worldwide in January 2020 and has now reached
virtually all communities on the planet. This short review provides evidence of
this spread and documents modelling efforts undertaken to understand and
forecast it, including a short section about the new variants that emerged in
late 2020.


  Morphological regeneration is an important feature that highlights the
environmental adaptive capacity of biological systems. Lack of this
regenerative capacity significantly limits the resilience of machines and the
environments they can operate in. To aid in addressing this gap, we develop an
approach for simulated soft robots to regrow parts of their morphology when
being damaged. Although numerical simulations using soft robots have played an
important role in their design, evolving soft robots with regenerative
capabilities have so far received comparable little attention. Here we propose
a model for soft robots that regenerate through a neural cellular automata.
Importantly, this approach only relies on local cell information to regrow
damaged components, opening interesting possibilities for physical regenerable
soft robots in the future. Our approach allows simulated soft robots that are
damaged to partially regenerate their original morphology through local cell
interactions alone and regain some of their ability to locomote. These results
take a step towards equipping artificial systems with regenerative capacities
and could potentially allow for more robust operations in a variety of
situations and environments. The code for the experiments in this paper is
available at: \url{github.com/KazuyaHoribe/RegeneratingSoftRobots}.


  The use of reinforcement learning has proven to be very promising for solving
complex activities without human supervision during their learning process.
However, their successful applications are predominantly focused on fictional
and entertainment problems - such as games. Based on the above, this work aims
to shed light on the application of reinforcement learning to solve this
relevant real-world problem, the genome assembly. By expanding the only
approach found in the literature that addresses this problem, we carefully
explored the aspects of intelligent agent learning, performed by the Q-learning
algorithm, to understand its suitability to be applied in scenarios whose
characteristics are more similar to those faced by real genome projects. The
improvements proposed here include changing the previously proposed reward
system and including state space exploration optimization strategies based on
dynamic pruning and mutual collaboration with evolutionary computing. These
investigations were tried on 23 new environments with larger inputs than those
used previously. All these environments are freely available on the internet
for the evolution of this research by the scientific community. The results
suggest consistent performance progress using the proposed improvements,
however, they also demonstrate the limitations of them, especially related to
the high dimensionality of state and action spaces. We also present, later, the
paths that can be traced to tackle genome assembly efficiently in real
scenarios considering recent, successfully reinforcement learning applications
- including deep reinforcement learning - from other domains dealing with
high-dimensional inputs.


  High-dimensional omics data contains intrinsic biomedical information that is
crucial for personalised medicine. Nevertheless, it is challenging to capture
them from the genome-wide data due to the large number of molecular features
and small number of available samples, which is also called 'the curse of
dimensionality' in machine learning. To tackle this problem and pave the way
for machine learning aided precision medicine, we proposed a unified multi-task
deep learning framework named OmiEmbed to capture biomedical information from
high-dimensional omics data with the deep embedding and downstream task
modules. The deep embedding module learnt an omics embedding that mapped
multiple omics data types into a latent space with lower dimensionality. Based
on the new representation of multi-omics data, different downstream task
modules were trained simultaneously and efficiently with the multi-task
strategy to predict the comprehensive phenotype profile of each sample.
OmiEmbed support multiple tasks for omics data including dimensionality
reduction, tumour type classification, multi-omics integration, demographic and
clinical feature reconstruction, and survival prediction. The framework
outperformed other methods on all three types of downstream tasks and achieved
better performance with the multi-task strategy comparing to training them
individually. OmiEmbed is a powerful and unified framework that can be widely
adapted to various application of high-dimensional omics data and has a great
potential to facilitate more accurate and personalised clinical decision
making.


  We have further extended our compartmental model describing the spread of the
infection in Italy. The model is based on the assumption that the time
evolution of all of the observable quantities (number of people still positive
to the infection, hospitalized and fatalities cases, healed people, and total
number of people that has contracted the infection) depend on average
parameters, namely people diffusion coefficient, infection cross-section, and
population density. The model provides precious information on the tight
relationship between the variation of the reported infection cases and a well
defined observable physical quantity: the average number of people that lie
within the daily displacement area of any single person. The extension of the
model now includes self-consistent evaluation of the reproduction index, effect
of immunization due to vaccination, and potential impact of virus variants on
the dynamical evolution of the outbreak. The model fits very well the epidemic
data, and allows us to strictly relate the time evolution of the number of
hospitalized case and fatalities to the change of people mobility, vaccination
rate, and appearance of an initial concentration of people positives for new
variants of the virus.


  Infectious diseases are a significant threat to human society which was over
sighted before the incidence of COVID-19, although according to the report of
the World Health Organisation (WHO) about 4.2 million people die annually due
to infectious disease. Due to recent COVID-19 pandemic, more than 2 million
people died during 2020 and 96.2 million people got affected by this
devastating disease. Recent research shows that applying individual
interactions and movements data could help managing the pandemic though
modelling the spread of infectious diseases on social contact networks.
Infectious disease spreading can be explained with the theories and methods of
diffusion processes where a dynamic phenomena evolves on networked systems. In
the modelling of diffusion process, it is assumed that contagious items spread
out in the networked system through the inter-node interactions. This resembles
spreading of infectious virus, e.g. spread of COVID-19, within a population
through individual social interactions. The evolution behaviours of the
diffusion process are strongly influenced by the characteristics of the
underlying system and the mechanism of the diffusion process itself. Thus,
spreading of infectious disease can be explained how people interact with each
other and by the characteristics of the disease itself. This paper presenters
the relevant theories and methodologies of diffusion process that can be used
to model the spread of infectious diseases.


  In late 2019, a novel coronavirus, the SARS-CoV-2 outbreak was identified in
Wuhan, China and later spread to every corner of the globe. Whilst the number
of infection-induced deaths in Ghana, West Africa are minimal when compared
with the rest of the world, the impact on the local health service is still
significant. Compartmental models are a useful framework for investigating
transmission of diseases in societies. To understand how the infection will
spread and how to limit the outbreak. We have developed a modified SEIR
compartmental model with nine compartments (CoVCom9) to describe the dynamics
of SARS-CoV-2 transmission in Ghana. We have carried out a detailed
mathematical analysis of the CoVCom9, including the derivation of the basic
reproduction number, $\mathcal{R}_{0}$. In particular, we have shown that the
disease-free equilibrium is globally asymptotically stable when
$\mathcal{R}_{0}<1$ via a candidate Lyapunov function. Using the SARS-CoV-2
reported data for confirmed-positive cases and deaths from March 13 to August
10, 2020, we have parametrised the CoVCom9 model. The results of this fit show
good agreement with data. We used Latin hypercube sampling-rank correlation
coefficient (LHS-PRCC) to investigate the uncertainty and sensitivity of
$\mathcal{R}_{0}$ since the results derived are significant in controlling the
spread of SARS-CoV-2. We estimate that over this five month period, the basic
reproduction number is given by $\mathcal{R}_{0} = 3.110$, with the 95\%
confidence interval being $2.042 \leq \mathcal{R}_0 \leq 3.240$, and the mean
value being $\mathcal{R}_{0}=2.623$. Of the 32 parameters in the model, we find
that just six have a significant influence on $\mathcal{R}_{0}$, these include
the rate of testing, where an increasing testing rate contributes to the
reduction of $\mathcal{R}_{0}$.


  A macroscopic theory for describing cellular states during steady-growth is
presented, which is based on the consistency between cellular growth and
molecular replication, as well as the robustness of phenotypes against
perturbations. Adaptive changes in high-dimensional phenotypes were shown to be
restricted within a low-dimensional slow manifold, from which a macroscopic law
for cellular states was derived, which was confirmed by adaptation experiments
on bacteria under stress. Next, the theory was extended to phenotypic
evolution, leading to proportionality between phenotypic responses against
genetic evolution and environmental adaptation. The link between robustness to
noise and mutation, as a result of robustness in developmental dynamics to
perturbations, showed proportionality between phenotypic plasticity by genetic
changes and by environmental noise. Accordingly, directionality and constraint
in phenotypic evolution was quantitatively formulated in terms of phenotypic
fluctuation and the response against environmental change. The evolutionary
relevance of slow modes in controlling high-dimensional phenotypes is
discussed.


  The development of therapeutic targets for COVID-19 treatment is based on the
understanding of the molecular mechanism of pathogenesis. The identification of
genes and proteins involved in the infection mechanism is the key to shed out
light into the complex molecular mechanisms. The combined effort of many
laboratories distributed throughout the world has produced the accumulation of
both protein and genetic interactions. In this work we integrate these
available results and we obtain an host protein-protein interaction network
composed by 1432 human proteins. We calculate network centrality measures to
identify key proteins. Then we perform functional enrichment of central
proteins. We observed that the identified proteins are mostly associated with
several crucial pathways, including cellular process, signalling transduction,
neurodegenerative disease. Finally, we focused on proteins involved in causing
disease in the human respiratory tract. We conclude that COVID19 is a complex
disease, and we highlighted many potential therapeutic targets including RBX1,
HSPA5, ITCH, RAB7A, RAB5A, RAB8A, PSMC5, CAPZB, CANX, IGF2R, HSPA1A, which are
central and also associated with multiple diseases


  Non-genetic forms of antimicrobial drug resistance can result from
cell-to-cell variability that is not encoded in the genetic material. Data from
recent studies also suggest that non-genetic mechanisms can facilitate the
development of genetic drug resistance. In this Perspective article, we
speculate on how the interplay between non-genetic and genetic mechanisms may
affect microbial adaptation and evolution during drug treatment. We argue that
cellular heterogeneity arising from fluctuations in gene expression, epigenetic
modifications, as well as genetic changes contributes to drug resistance at
different timescales, and that the interplay between these mechanisms may
influence the evolutionary dynamics of pathogen resistance. Accordingly,
developing a better understanding of non-genetic mechanisms in drug resistance
and how they interact with genetic mechanisms will enhance our ability to
combat antimicrobial resistance.


  Motivation: The design of enzymes is as challenging as it is consequential
for making chemical synthesis in medical and industrial applications more
efficient, cost-effective and environmentally friendly. While several aspects
of this complex problem are computationally assisted, the drafting of catalytic
mechanisms, i.e. the specification of the chemical steps-and hence intermediate
states-that the enzyme is meant to implement, is largely left to human
expertise. The ability to capture specific chemistries of multi-step catalysis
in a fashion that enables its computational construction and design is
therefore highly desirable and would equally impact the elucidation of existing
enzymatic reactions whose mechanisms are unknown. Results: We use the
mathematical framework of graph transformation to express the distinction
between rules and reactions in chemistry. We derive about 1000 rules for amino
acid side chain chemistry from the M-CSA database, a curated repository of
enzymatic mechanisms. Using graph transformation we are able to propose
hundreds of hypothetical catalytic mechanisms for a large number of unrelated
reactions in the Rhea database. We analyze these mechanisms to find that they
combine in chemically sound fashion individual steps from a variety of known
multi-step mechanisms, showing that plausible novel mechanisms for catalysis
can be constructed computationally.


  The spread of the coronavirus disease 2019 (COVID-19) caused by severe acute
respiratory syndrome coronavirus-2 (SARS-CoV-2) has become a global health
crisis. The binding affinity of SARS-CoV-2 (in particular the receptor binding
domain, RBD) to its receptor angiotensin converting enzyme 2 (ACE2) and the
antibodies is of great importance in understanding the infectivity of COVID-19
and evaluating the candidate therapeutic for COVID-19. In this work, we propose
a new method based on molecular mechanics/Poisson-Boltzmann surface area
(MM/PBSA) to accurately calculate the free energy of SARS-CoV-2 RBD binding to
ACE2 and antibodies. The calculated binding free energy of SARS-CoV-2 RBD to
ACE2 is -13.3 kcal/mol, and that of SARS-CoV RBD to ACE2 is -11.4 kcal/mol,
which agrees well with experimental result (-11.3 kcal/mol and -10.1 kcal/mol,
respectively). Moreover, we take two recently reported antibodies as the
example, and calculate the free energy of antibodies binding to SARS-CoV-2 RBD,
which is also consistent with the experimental findings. Further, within the
framework of the modified MM/PBSA, we determine the key residues and the main
driving forces for the SARS-CoV-2 RBD/CB6 interaction by the computational
alanine scanning method. The present study offers a computationally efficient
and numerically reliable method to evaluate the free energy of SARS-CoV-2
binding to other proteins, which may stimulate the development of the
therapeutics against the COVID-19 disease in real applications.


  The availability of genomic data is essential to progress in biomedical
research, personalized medicine, etc. However, its extreme sensitivity makes it
problematic, if not outright impossible, to publish or share it. As a result,
several initiatives have been launched to experiment with synthetic genomic
data, e.g., using generative models to learn the underlying distribution of the
real data and generate artificial datasets that preserve its salient
characteristics without exposing it. This paper provides the first evaluation
of both utility and privacy protection of six state-of-the-art models for
generating synthetic genomic data. We assess the performance of the synthetic
data on several common tasks, such as allele population statistics and linkage
disequilibrium. We then measure privacy through the lens of membership
inference attacks, i.e., inferring whether a record was part of the training
data. Our experiments show that no single approach to generate synthetic
genomic data yields both high utility and strong privacy across the board.
Also, the size and nature of the training dataset matter. Moreover, while some
combinations of datasets and models produce synthetic data with distributions
close to the real data, there often are target data points that are vulnerable
to membership inference. Looking forward, our techniques can be used by
practitioners to assess the risks of deploying synthetic genomic data in the
wild and serve as a benchmark for future work.


  As Europe is facing the second wave of the CoViD-19 pandemic, each country
should carefully review how it dealt with the first wave of outbreak. Lessons
from the first experience should be useful to avoid indiscriminate closures
and, above all, to determine universal (understandable) parameters to guide the
introduction of containment measures to reduce the spreading of the virus. The
use of few (effective) parameters is indeed of extreme importance to create a
link between authorities and population, allowing the latter to understand the
reason for some restrictions and, consequently, to allow an active
participation in the fight against the pandemic.
  Testing strategies, fitting skew parameters (as mean, mode, standard
deviation, and skewness), mortality rates, and weekly CoViD-19 spreading data,
as more people are getting infected, were used to compare the first wave of the
outbreak in the Italian regions and to determine which parameters have to be
checked before introducing restrictive containment measures. We propose few
\textit{universal} parameters that, once appropriately weighed, could be useful
to correctly differentiate the pandemic situation in the national territory and
to rapidly assign the properly pandemic risk to each region.


  Phylogenetic networks can model more complicated evolutionary phenomena that
trees fail to capture such as horizontal gene transfer and hybridization. The
same Markov models that are used to model evolution on trees can also be
extended to networks and similar questions, such as the identifiability of the
network parameter or the invariants of the model, can be asked. In this paper
we focus on finding the invariants of the Cavendar-Farris-Neyman (CFN) model on
level-1 phylogenetic networks. We do this by reducing the problem to finding
invariants of sunlet networks, which are level-1 networks consisting of a
single cycle with leaves at each vertex. We then determine all quadratic
invariants in the sunlet network ideal which we conjecture generate the full
ideal.


  Networks provide a meaningful way to represent and analyze complex biological
information, but the methodological details of network-based tools are often
described for a technical audience. Graphery is a hands-on tutorial webserver
designed to help biological researchers understand the fundamental concepts
behind commonly-used graph algorithms. Each tutorial describes a graph concept
along with executable Python code that visualizes the concept in a code view
and a graph view. Graphery tutorials help researchers understand graph
statistics (such as degree distribution and network modularity) and classic
graph algorithms (such as shortest paths and random walks). Users navigate each
tutorial using their choice of real-world biological networks, ranging in scale
from molecular interaction graphs to ecological networks. Graphery also allows
users to modify the code within each tutorial or write new programs, which all
can be executed without requiring an account. Discipline-focused tutorials will
be essential to help researchers interpret their biological data. Graphery
accepts ideas for new tutorials and datasets that will be shaped by both
computational and biological researchers, growing into a community-contributed
learning platform. Availability: Graphery is available at
https://graphery.reedcompbio.org/.


  High-throughput sequencing of T- and B-cell receptors makes it possible to
track immune repertoires across time, in different tissues, in acute and
chronic diseases and in healthy individuals. However quantitative comparison
between repertoires is confounded by variability in the read count of each
receptor clonotype due to sampling, library preparation, and expression noise.
We review methods for accounting for both biological and experimental noise and
present an easy-to-use python package NoisET that implements and generalizes a
previously developed Bayesian method. It can be used to learn experimental
noise models for repertoire sequencing from replicates, and to detect
responding clones following a stimulus. We test the package on different
repertoire sequencing technologies and datasets. We review how such approaches
have been used to identify responding clonotypes in vaccination and disease
data. Availability: NoisET is freely available to use with source code at
github.com/statbiophys/NoisET.


  Background: Early detection and isolation of COVID-19 patients are essential
for successful implementation of mitigation strategies and eventually curbing
the disease spread. With a limited number of daily COVID-19 tests performed in
every country, simulating the COVID-19 spread along with the potential effect
of each mitigation strategy currently remains one of the most effective ways in
managing the healthcare system and guiding policy-makers. Methods: We introduce
COVIDHunter, a flexible and accurate COVID-19 outbreak simulation model that
evaluates the current mitigation measures that are applied to a region and
provides suggestions on what strength the upcoming mitigation measure should
be. The key idea of COVIDHunter is to quantify the spread of COVID-19 in a
geographical region by simulating the average number of new infections caused
by an infected person considering the effect of external factors, such as
environmental conditions (e.g., climate, temperature, humidity) and mitigation
measures. Results: Using Switzerland as a case study, COVIDHunter estimates
that if the policy-makers relax the mitigation measures by 50% for 30 days then
both the daily capacity need for hospital beds and daily number of deaths
increase exponentially by an average of 5.1x, who may occupy ICU beds and
ventilators for a period of time. Unlike existing models, the COVIDHunter model
accurately monitors and predicts the daily number of cases, hospitalizations,
and deaths due to COVID-19. Our model is flexible to configure and simple to
modify for modeling different scenarios under different environmental
conditions and mitigation measures. Availability: We release the source code of
the COVIDHunter implementation at https://github.com/CMU- SAFARI/COVIDHunter
and show how to flexibly configure our model for any scenario and easily extend
it for different measures and conditions than we account for.


  Realistic models of biological processes typically involve interacting
components on multiple scales, driven by changing environment and inherent
stochasticity. Such models are often analytically and numerically intractable.
We revisit a dynamic maximum entropy method that combines a static maximum
entropy and a quasi-stationary approximation. This allows us to reduce
stochastic non-equilibrium dynamics expressed by the Fokker-Planck equation to
a simpler low-dimensional deterministic dynamics, without the need to track
microscopic details. Although the method has been previously applied to a few
(rather complicated) applications in population genetics, our main goal here is
to explain and to better understand how the method works. We demonstrate the
usefulness of the method for two widely studied stochastic problems,
highlighting its accuracy in capturing important macroscopic quantities even in
rapidly changing non-stationary conditions. For the Ornstein-Uhlenbeck process,
the method recovers the exact dynamics whilst for a stochastic island model
with migration from other habitats, the approximation retains high macroscopic
accuracy under a wide range of scenarios for a dynamic environment.


  We study local conformational biases in the dynamics of {\alpha}-synuclein by
using all-atom simulations with explicit and implicit solvents. The biases are
related to the frequency of the specific contact formation. In both approaches,
the protein is intrinsically disordered, and its strongest bias is to make bend
and turn local structures. The explicit-solvent conformations can be
substantially more extended which allows for formation of transient trefoil
knots, both deep and shallow, that may last for up to 5 {\mu}s. The two-chain
self-association events, both short- and long-lived, are dominated by formation
of contacts in the central part of the sequence. This part tends to form
helices when bound to a micelle.


  Biomolecular condensates such as membraneless organelles, underpinned by
liquid-liquid phase separation (LLPS), are important for physiological
function, with electrostatics -- among other interaction types -- being a
prominent force in their assembly. Charge interactions of intrinsically
disordered proteins (IDPs) and other biomolecules are sensitive to the aqueous
dielectric environment. Because the relative permittivity of protein is
significantly lower than that of water, the interior of an IDP condensate is a
relatively low-dielectric regime, which, aside from its possible functional
effects on client molecules, should facilitate stronger electrostatic
interactions among the scaffold IDPs. To gain insight into this LLPS-induced
dielectric heterogeneity, addressing in particular whether a low-dielectric
condensed phase entails more favorable LLPS than that posited by assuming IDP
electrostatic interactions are uniformly modulated by the higher dielectric
constant of the pure solvent, we consider a simplified multiple-chain model of
polyampholytes immersed in explicit solvents that are either polarizable or
possess a permanent dipole. Notably, simulated phase behaviors of these systems
exhibit only minor to moderate differences from those obtained using
implicit-solvent models with a uniform relative permittivity equals to that of
pure solvent. Buttressed by theoretical treatments developed here using random
phase approximation and polymer field-theoretic simulations, these observations
indicate a partial compensation of effects between favorable solvent-mediated
interactions among the polyampholytes in the condensed phase and favorable
polyampholyte-solvent interactions in the dilute phase, often netting only a
minor enhancement of overall LLPS propensity from the very dielectric
heterogeneity that arises from the LLPS itself. Further ramifications of this
principle are discussed.


  Motivation: Bacteriophages (aka phages), which mainly infect bacteria, play
key roles in the biology of microbes. As the most abundant biological entities
on the planet, the number of discovered phages is only the tip of the iceberg.
Recently, many new phages have been revealed using high throughput sequencing,
particularly metagenomic sequencing. Compared to the fast accumulation of
phage-like sequences, there is a serious lag in taxonomic classification of
phages. High diversity, abundance, and limited known phages pose great
challenges for taxonomic analysis. In particular, alignment-based tools have
difficulty in classifying fast accumulating contigs assembled from metagenomic
data. Results: In this work, we present a novel semi-supervised learning model,
named PhaGCN, to conduct taxonomic classification for phage contigs. In this
learning model, we construct a knowledge graph by combining the DNA sequence
features learned by convolutional neural network (CNN) and protein sequence
similarity gained from gene-sharing network. Then we apply graph convolutional
network (GCN) to utilize both the labeled and unlabeled samples in training to
enhance the learning ability. We tested PhaGCN on both simulated and real
sequencing data. The results clearly show that our method competes favorably
against available phage classification tools.


  We study effects of strategy-dependent time delays on equilibria of evolving
populations. It is well known that time delays may cause oscillations in
dynamical systems. Here we report a novel behavior. We show that microscopic
models of evolutionary games with strategy-dependent time delays lead to a new
type of replicator dynamics. It describes the time evolution of fractions of
the population playing given strategies and the size of the population. Unlike
in all previous models, stationary states of such dynamics depend continuously
on time delays. We show that in games with an interior stationary state (a
globally asymptotically stable equilibrium in the standard replicator
dynamics), at certain time delays, it may disappear or there may appear another
interior stationary state. In the Prisoner's Dilemma game, for time delays of
cooperation smaller than time delays of defection, there appears an ustable
interior equilibrium and therefore for some initial conditions, the population
converges to the homogeneous state with just cooperators.


  The control of Covid 19 epidemics by public health policy in Italy during the
first and the second epidemic waves has been driven by using reproductive
number Rt(t) to identify the supercritical (percolative), the subcritical
(arrested), separated by the critical regime. Here we show that to quantify the
Covid-19 spreading rate with containment measures (CSRwCM) there is a need of a
3D expanded parameter space phase diagram built by the combination of Rt(t) and
doubling time Td(t). In this space we identify the dynamics of the Covid-19
dynamics Italy and its administrative Regions. The supercritical regime is
mathematically characterized by i) the power law of Td vs. [Rt(t)-1] and ii)
the exponential behaviour of Td vs. time, either in the first and in the second
wave. The novel 3D phase diagram shows clearly metastable states appearing
before and after the second wave critical regime. for loosening quarantine and
tracing of actives cases. The metastable states are precursors of the abrupt
onset of a next nascent wave supercritical regime. This dynamic description
allows epidemics predictions needed by policymakers to activate
non-pharmaceutical interventions (NPIs), a key issue for avoiding economical
losses, reduce fatalities and avoid new virus variant during vaccination
campaign


  Introduction. Can the infection due to the human immunodeficiency virus type
1 induce a change in the differentiation status or process in T cells?.
Methods. We will consider two stochastic Markov chain models, one which will
describe the T-helper cell differentiation process, and another one describing
that process of infection of the T-helper cell by the virus; in these Markov
chains, we will consider a set of states $\{X_t \}$ comprised of those proteins
involved in each of the processes and their interactions (either
differentiation or infection of the cell), such that we will obtain two
stochastic transition matrices ($A,B$), one for each process; afterwards, the
computation of their eigenvalues shall be performed, in which, should the
eigenvalue $\lambda_i=1$ exist, the computation for the equilibrium
distribution $\pi^n$ will be obtained for each of the matrices, which will
inform us on the trends of interactions amongst the proteins in the long-term.
Results. The stochastic processes considered possess an equilibrium
distribution, when reaching their equilibrium distribution, there exists an
increase in their informational entropy, and their log-rank distributions can
be modeled as discrete beta generalized distributions (DGBD). Discussion. The
equilibrium distributions of both process can be regarded as states in which
the cell is well-differentiated, ergo there exists an induction of a novel
HIV-dependent differentiated state in the T-cell; these processes due to their
DGBD distribution can be considered complex processes; due to the increasing
entropy, the equilibrium states are stable ones. Conclusion. The HIV virus can
promote a novel differentiated state in the T-cell, which can give account for
clinical features seen in patients; this model, notwithstanding does not give
account of YES/NO logical switches involved in the regulatory networks.


  Stochastic hybrid systems involve a coupling between a discrete Markov chain
and a continuous stochastic process. If the latter evolves deterministically
between jumps in the discrete state, then the system reduces to a piecewise
deterministic Markov process (PDMP). Well known examples include stochastic
gene expression, voltage fluctuations in neurons, and motor-driven
intracellular transport. In this paper we use coherent spin states to construct
a new path integral representation of the probability density functional for
stochastic hybrid systems, which holds outside the weak noise regime. We use
the path integral to derive a system of Langevin equations in the
semi-classical limit, which extends previous diffusion approximations based on
a quasi-steady-state reduction. We then show how in the weak noise limit the
path integral is equivalent to an alternative representation that was
previously derived using Doi-Peliti operators. The action functional of the
latter is related to a large deviation principle for stochastic hybrid systems.


  Recent advancements in machine learning techniques for protein folding
motivate better results in its inverse problem -- protein design. In this work
we introduce a new graph mimetic neural network, MimNet, and show that it is
possible to build a reversible architecture that solves the structure and
design problems in tandem, allowing to improve protein design when the
structure is better estimated. We use the ProteinNet data set and show that the
state of the art results in protein design can be improved, given recent
architectures for protein folding.


  Exploring the possible consequences of spatial reciprocity on the evolution
of cooperation is an intensively studied research avenue. Related works assumed
a certain interaction graph of competing players and studied how particular
topologies may influence the dynamical behavior. In this paper we apply a
numerically more demanding off-lattice population approach which could be
potentially relevant especially in microbiological environments. As expected,
results are conceptually similar to those which were obtained for lattice-type
interaction graphs, but some spectacular differences can also be revealed. On
one hand, in off-lattice populations spatial reciprocity may work more
efficiently than for a lattice-based system. On the other hand, competing
strategies may separate from each other in the continuous space concept, which
gives a chance for cooperators to survive even at relatively high temptation
values. Furthermore, the lack of strict neighborhood results in soft borders
between competing patches which jeopardizes the long term stability of
homogeneous domains. We survey the major social dilemma games based on pair
interactions of players and reveal all analogies and differences compared to
on-lattice simulations.


  There is a near consensus view that SARS-CoV-2 has a natural zoonotic origin;
however, several characteristics of SARS-CoV-2 taken together are not easily
explained by a natural zoonotic origin hypothesis. These include: a low rate of
evolution in the early phase of transmission; the lack of evidence of
recombination events; a high pre-existing binding to human ACE2; a novel furin
cleavage site insert; a flat glycan binding domain of the spike protein which
conflicts with host evasion survival patterns exhibited by other coronaviruses,
and high human and mouse peptide mimicry. Initial assumptions against a
laboratory origin, by contrast, have remained unsubstantiated. Furthermore,
over a year after the initial outbreak in Wuhan, there is still no clear
evidence of zoonotic transfer from a bat or intermediate species. Given the
immense social and economic impact of this pandemic, identifying the true
origin of SARS-CoV-2 is fundamental to preventing future outbreaks. The search
for SARS-CoV-2's origin should include an open and unbiased inquiry into a
possible laboratory origin.


  We present an algorithm for building the extended BWT (eBWT) of a string
collection from its grammar-compressed representation. Our technique exploits
the string repetitions captured by the grammar to boost the computation of the
eBWT. Thus, the more repetitive the collection is, the lower are the resources
we use per input symbol. We rely on a new grammar recently proposed at DCC'21
whose nonterminals serve as building blocks for inducing the eBWT. A relevant
application for this idea is the construction of self-indexes for analyzing
sequencing reads -- massive and repetitive string collections of raw genomic
data. Self-indexes have become increasingly popular in Bioinformatics as they
can encode more information in less space. Our efficient eBWT construction
opens the door to perform accurate bioinformatic analyses on more massive
sequence datasets, which are not tractable with current eBWT construction
techniques.


  One pivotal feature of transcriptomics data is the unwanted variations caused
by disparate experimental handling, known as handling effects. Various data
normalization methods were developed to alleviate the adverse impact of
handling effects in the setting of differential expression analysis. However,
little research has been done to evaluate their performance in the setting of
survival outcome prediction, an important analysis goal for transcriptomics
data in biomedical research. Leveraging a unique pair of datasets for the same
set of tumor samples-one with handling effects and the other without, we
developed a benchmarking tool for conducting such an evaluation in microRNA
microarrays. We applied this tool to evaluate the performance of three popular
normalization methods-quantile normalization, median normalization, and
variance stabilizing normalization-in survival prediction using various
approaches for model building and designs for sample assignment. We showed that
handling effects can have a strong impact on survival prediction, and that
quantile normalization, a most popular method in current practice, tends to
underperform median normalization and variance stabilizing normalization. We
demonstrated with a small example the reason for quantile normalization's poor
performance in this setting. Our finding highlights the importance of putting
normalization evaluation in the context of the downstream analysis setting and
the potential of improving the development of survival predictors by applying
median normalization. We make available our benchmarking tool for performing
such evaluation on additional normalization methods in connection with
prediction modeling approaches.


  Nowadays, due to the increasing amount of experimental data obtained by
sequencing, the most interest is focused on determining the functions and
characteristics of its individual parts of the genome instead of determining
the nucleotide sequence of the genome. The genome annotation includes the
identification of coding and non-coding sequences, determining the structure of
the gene and determining the functions of these sequences. Despite the
significant achievements in computational technologies working with sequencing
data, there is no general approach to the functional annotation of the genome
in the reason of the large number of unresolved molecular determination of the
function of some genomes parts. Nevertheless, the scientific community is
trying to solve this problem. This review analyzed existing approaches to
eukaryotic genome annotation. This work includes 3 main parts: introduction,
main body and discussion. The introduction reflects the development of
independent tools and automatic pipelines for annotation of eukaryotic genomes,
which are associated with existing achievements in annotating prokaryotic ones.
The main body consists of two distinguished parts, the first one is devoted to
instructions for annotating genomes of non-model eukaryotes, and the second
block is about recent versions of automatic pipelines that require minimal
user's curation. The question of assessing the quality and completeness of the
annotated genome is noted briefly, and the tools to conduct this analysis are
discussed. Currently, there is no universal automatic software for eukaryotic
genome annotation, covering the whole list of tasks, without manual curation or
using additional external tools and resources. Thus it leads to the task of
developing a wider functional and universal protocol for automatic annotation
of small eukaryotic genomes.


  Reconstructing interactions from observational data is a critical need for
investigating natural biological networks, wherein network dimensionality (i.e.
number of interacting components) is usually high and interactions are
time-varying. These pose a challenge to existing methods that can quantify only
small interaction networks or assume static interactions under steady state.
Here, we proposed a novel approach to reconstruct high-dimensional,
time-varying interaction networks using empirical time series. This method,
named "multiview distance regularized S-map", generalized the state space
reconstruction to accommodate high dimensionality and overcome difficulties in
quantifying massive interactions with limited data. When we evaluated this
method using the time series generated from a large theoretical model involving
hundreds of interacting species, estimated interaction strengths were in good
agreement with theoretical expectations. As a result, reconstructed networks
preserved important topological properties, such as centrality, strength
distribution and derived stability measures. Moreover, our method effectively
forecasted the dynamic behavior of network nodes. Applying this method to a
natural bacterial community helped identify keystone species from the
interaction network and revealed the mechanisms governing the dynamical
stability of bacterial community. Our method overcame the challenge of high
dimensionality and disentangled complex time-varying interactions in large
natural dynamical systems.


  The novel coronavirus disease (COVID-19) pneumonia has posed a great threat
to the world recent months by causing many deaths and enormous economic damage
worldwide. The first case of COVID-19 in Morocco was reported on 2 March 2020,
and the number of reported cases has increased day by day. In this work, we
extend the well-known SIR compartmental model to deterministic and stochastic
time-delayed models in order to predict the epidemiological trend of COVID-19
in Morocco and to assess the potential role of multiple preventive measures and
strategies imposed by Moroccan authorities. The main features of the work
include the well-posedness of the models and conditions under which the
COVID-19 may become extinct or persist in the population. Parameter values have
been estimated from real data and numerical simulations are presented for
forecasting the COVID-19 spreading as well as verification of theoretical
results.


  Graph neural networks (GNNs) have been proposed for a wide range of
graph-related learning tasks. In particular, in recent years, an increasing
number of GNN systems were applied to predict molecular properties. However, a
direct impediment is to select appropriate hyperparameters to achieve
satisfactory performance with lower computational cost. Meanwhile, many
molecular datasets are far smaller than many other datasets in typical deep
learning applications. Most hyperparameter optimization (HPO) methods have not
been explored in terms of their efficiencies on such small datasets in the
molecular domain. In this paper, we conducted a theoretical analysis of common
and specific features for two state-of-the-art and popular algorithms for HPO:
TPE and CMA-ES, and we compared them with random search (RS), which is used as
a baseline. Experimental studies are carried out on several benchmarks in
MoleculeNet, from different perspectives to investigate the impact of RS, TPE,
and CMA-ES on HPO of GNNs for molecular property prediction. In our
experiments, we concluded that RS, TPE, and CMA-ES have their individual
advantages in tackling different specific molecular problems. Finally, we
believe our work will motivate further research on GNN as applied to molecular
machine learning problems in chemistry and materials sciences.


  We introduce a novel data-driven framework for the design of targeted gene
panels for estimating exome-wide biomarkers in cancer immunotherapy. Our first
goal is to develop a generative model for the profile of mutation across the
exome, which allows for gene- and variant type-dependent mutation rates. Based
on this model, we then propose a new procedure for estimating biomarkers such
as tumour mutation burden and tumour indel nurden. Our approach allows the
practitioner to select a targeted gene panel of a prespecified size, and then
construct an estimator that only depends on the selected genes. Alternatively,
the practitioner may apply our method to make predictions based on an existing
gene panel, or to augment a gene panel to a given size. We demonstrate the
excellent performance of our proposal using data from three non-small cell lung
cancer studies, as well as data from six other cancer types.


  In a recent study, Bryant, Francis and Steel investigated the concept of
\enquote{future-proofing} consensus methods in phylogenetics. That is, they
investigated if such methods can be robust against the introduction of
additional data like extra trees or new species. In the present manuscript, we
analyze consensus methods under a different aspect of introducing new data,
namely concerning the discovery of new clades. In evolutionary biology, often
formerly unresolved clades get resolved by refined reconstruction methods or
new genetic data analyses. In our manuscript we investigate which properties of
consensus methods can guarantee that such new insights do not disagree with
previously found consensus trees but merely refine them. We call consensus
methods with this property \emph{refinement-stable}. Along these lines, we also
study two famous super tree methods, namely Matrix Representation with
Parsimony (MRP) and Matrix Representation with Compatibility (MRC), which have
also been suggested as consensus methods in the literature. While we (just like
Bryant, Francis and Steel in their recent study) unfortunately have to conclude
some negative answers concerning general consensus methods, we also state some
relevant and positive results concerning the majority rule (MR) and strict
consensus methods, which are amongst the most frequently used consensus
methods. Moreover, we show that there exist infinitely many consensus methods
which are refinement-stable and have some other desirable properties.


  While the analytical solution for the marginal distribution of a stochastic
chemical reaction network has been extensively studied, its joint distribution,
i.e. the solution of a high-dimensional chemical master equation, has received
much less attention. Here we develop a novel method of computing the exact
joint distributions of a wide class of first-order stochastic reaction systems
in steady-state conditions. The effectiveness of our method is validated by
applying it to four gene expression models of biological significance,
including models with 2A peptides, nascent mRNA, gene regulation, translational
bursting, and alternative splicing.


  This paper concerns the modeling and numerical simulation of the process of
speciation. In particular, given conditions for which one or more speciation
events within an ecosystem occur, our aim is to develop the necessary modeling
and simulation tools. Care is also taken to establish a solid mathematical
foundation on which our modeling framework is built. This is the subject of the
first half of the paper. The second half is devoted to developing a multi-scale
framework for eco-evolutionary modeling, where the relevant scales are that of
species and individual/population, respectively. Hence, a system of interacting
species can be described at the species level, while for branching species a
population level description is necessary. Our multi-scale framework thus
consists of coupling the species and population level models where speciation
events are detected in advance and then resolved at the population scale until
the branching is complete. Moreover, since the population level model is
formulated as a PDE, we first establish the well-posedness in the time-discrete
setting, and then derive the a posteriori error estimates which provides a
fully computable upper bound on an energy-type error, including also for the
case of general smooth distributions (which will be useful for the detection of
speciation events). Several numerical tests validate our framework in practice.


  Collective migration of cells and animals often relies on a specialised set
of "leaders", whose role is to steer a population of naive followers towards
some target. We formulate a continuous model to understand the dynamics and
structure of such groups, splitting a population into separate follower and
leader types with distinct orientation responses. We incorporate "leader
influence" via three principal mechanisms: a bias in the orientation of leaders
according to the destination, distinct speeds of movement and distinct levels
of conspicuousness. Using a combination of analysis and numerical computation
on a sequence of models of increasing complexity, we assess the extent to which
leaders successfully shepherd the swarm. While all three mechanisms can lead to
a successfully steered swarm, parameter regime is crucial with non successful
choices generating a variety of unsuccessful attempts, including movement away
from the target, swarm splitting or swarm dispersal.


  The outbreak of the coronavirus disease 2019 (COVID-19) caused by severe
acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has become a worldwide
public health crisis. When the SARS-CoV-2 enters the biological fluids in the
human body, different types of biomolecules (in particular proteins) may adsorb
on its surface and alter its infection ability. Although great efforts have
recently been devoted to the interaction of the specific antibodies with the
SARS-CoV-2, it still remains largely unknown how the other serum proteins
affect the infection of the SARS-CoV-2. In this work, we systematically
investigate the interaction of serum proteins with the SARS-CoV-2 RBD by the
molecular docking and the all-atom molecular dynamics simulations. It is found
that the non-specific immunoglobulin (Ig) indeed cannot effectively bind to the
SARS-CoV-2 RBD while the human serum albumin (HSA) may have some potential of
blocking its infection (to ACE2). More importantly, we find that the RBD can
cause the significant structural change of the Apolipoprotein E (ApoE), by
which SARS-CoV-2 may hijack the metabolic pathway of the ApoE to facilitate its
cell entry. The present study enhances the understanding of the role of protein
corona in the bio-behaviors of SARS-CoV-2, which may aid the more precise and
personalized treatment for COVID-19 infection in the clinic.


  Digital contact tracing is a public health intervention. It should be
integrated with local health policy, provide rapid and accurate notifications
to exposed individuals, and encourage high app uptake and adherence to
quarantine. Real-time monitoring and evaluation of effectiveness of app-based
contact tracing is key for improvement and public trust.


  The long-term growth rate of populations in varying environments quantifies
the evolutionary value of processing the information that biological
individuals inherit from their ancestors and acquire from their environment.
Previous models were limited to asexual reproduction with inherited information
coming from a single parent with no recombination. We present a general
extension to sexual reproduction and an analytical solution for a particular
but important case, the infinitesimal model of quantitative genetics which
assumes traits to be normally distributed. We study with this model the
conditions under which sexual reproduction is advantageous and can evolve in
the context of autocorrelated or directionally varying environments, mutational
biases, spatial heterogeneities and phenotypic plasticity. Our results
generalize and unify previous analyses. We also examine the proposal made by
Geodakyan that the presence of two phenotypically distinct sexes permits an
optimal adaptation to varying environments. We verify that conditions exists
where sexual dimorphism is adaptive but find that its evolutionary value does
not generally compensate for the two-fold cost of males.


  The present paper is devoted to the investigation of the long term behavior
of a class of singular multi-dimensional diffusion processes that get absorbed
in finite time with probability one. Our focus is on the analysis of
quasi-stationary distributions (QSDs), which describe the long term behavior of
the system conditioned on not being absorbed. Under natural Lyapunov
conditions, we construct a QSD and prove the sharp exponential convergence to
this QSD for compactly supported initial distributions. Under stronger Lyapunov
conditions ensuring that the diffusion process comes down from infinity, we
show the uniqueness of a QSD and the exponential convergence to the QSD for all
initial distributions. Our results can be seen as the multi-dimensional
generalization of Cattiaux et al (Ann. Prob. 2009) as well as the complement to
Hening and Nguyen (Ann. Appl. Prob. 2018) which looks at the long term behavior
of multi-dimensional diffusions that can only become extinct asymptotically.
  The centerpiece of our approach concerns a uniformly elliptic operator that
we relate to the generator, or the Fokker-Planck operator, associated to the
diffusion process. This operator only has singular coefficients in its
zeroth-order terms and can be handled more easily than the generator. For this
operator, we establish the discreteness of its spectrum, its principal spectral
theory, the stochastic representation of the semigroup generated by it, and the
global regularity for the associated parabolic equation. We show how our
results can be applied to most ecological models, among which cooperative,
competitive, and predator-prey Lotka-Volterra systems.


  Ecology and evolution under changing environments are important in many
subfields of biology with implications for medicine. Here, we explore an
example: the consequences of fluctuating environments on the emergence of
antibiotic resistance, which is an immense and growing problem. Typically, high
doses of antibiotics are employed to eliminate the infection quickly and
minimize the time under which resistance may emerge. However, this strategy may
not be optimal. Since competition can reduce fitness and resistance typically
has a reproductive cost, resistant mutants' fitness can depend on their
environment. Here we show conditions under which environmental varying fitness
can be exploited to prevent the emergence of resistance. We develop a
stochastic Lotka-Volterra model of a microbial system with competing
phenotypes: a wild strain susceptible to the antibiotic, and a mutant strain
that is resistant. We investigate the impact of various pulsed applications of
antibiotics on population suppression. Leveraging competition, we show how a
strategy of environmental switching can suppress the infection while avoiding
resistant mutants. We discuss limitations of the procedure depending on the
microbe and pharmacodynamics and methods to ameliorate them.


  We present a two-dimensional continuum model of tumor growth, which treats
the tissue as a composition of six distinct fluid phases; their dynamics are
governed by the equations of mass and momentum conservation. Our model divides
the cancer cells phase into two sub-phases depending on their maturity state.
The same approach is also applied for the vasculature phase, which is divided
into young sprouts (products of angiogenesis), and fully formed-mature vessels.
The remaining two phases correspond to healthy cells and extracellular material
(ECM). Furthermore, the model foresees the existence of nutrient chemical
species, which are transferred within the tissue through diffusion or supplied
by the vasculature (blood vessels). The model is numerically solved with the
Finite Elements Method and computations are performed with the commercial
software Comsol Multiphysics. The numerical simulations predict that mature
cancer cells are well separated from young cancer cells, which form a
protective shield for the growing tumor. We study the effect of different
mitosis and death rates for mature and young cancer cells on the tumor growth
rate, and predict accelerated rates when the mitosis rate of young cancer cells
is higher compared to mature cancer cells.


  Structure-based virtual screening (SBVS) is a promising in silico technique
that integrates computational methods into drug design. An extensively used
method in SBVS is molecular docking. However, the docking process can hardly be
computationally efficient and accurate simultaneously because classic mechanics
scoring function is used to approximate, but hardly reach, the quantum
mechanics precision in this method. In order to reduce the computational cost
of the protein-ligand scoring process and use data driven approach to boost the
scoring function accuracy, we introduce a docking-based SBVS method and,
furthermore, a deep learning non-docking-based method that is able to avoid the
computational cost of the docking process. Then, we try to integrate these two
methods into an easy-to-use framework, ParaVS, that provides both choices for
researchers. Graph neural network (GNN) is employed in ParaVS, and we explained
how our in-house GNN works and how to model ligands and molecular targets. To
verify our approaches, cross validation experiments are done on two datasets,
an open dataset Directory of Useful Decoys: Enhanced (DUD.E) and an in-house
proprietary dataset without computational generated artificial decoys
(NoDecoy). On DUD.E we achieved a state-of-the-art AUC of 0.981 and a
state-of-the-art enrichment factor at 2% of 36.2; on NoDecoy we achieved an AUC
of 0.974. We further finish inference of an open database, Enamine REAL
Database (RDB), that comprises over 1.36 billion molecules in 4050 core-hours
using our ParaVS non-docking method (ParaVS-ND). The inference speed of
ParaVS-ND is about 3.6e5 molecule / core-hour, while this number of a
conventional docking-based method is around 20, which is about 16000 times
faster. The experiments indicate that ParaVS is accurate, computationally
efficient and can be generalized to different molecular.


  Cryo-electron microscopy (cryo-EM) has become a major experimental technique
to determine the structures of large protein complexes and molecular
assemblies, as evidenced by the 2017 Nobel Prize. Although cryo-EM has been
drastically improved to generate high-resolution three-dimensional (3D) maps
that contain detailed structural information about macromolecules, the
computational methods for using the data to automatically build structure
models are lagging far behind. The traditional cryo-EM model building approach
is template-based homology modeling. Manual de novo modeling is very
time-consuming when no template model is found in the database. In recent
years, de novo cryo-EM modeling using machine learning (ML) and deep learning
(DL) has ranked among the top-performing methods in macromolecular structure
modeling. Deep-learning-based de novo cryo-EM modeling is an important
application of artificial intelligence, with impressive results and great
potential for the next generation of molecular biomedicine. Accordingly, we
systematically review the representative ML/DL-based de novo cryo-EM modeling
methods. And their significances are discussed from both practical and
methodological viewpoints. We also briefly describe the background of cryo-EM
data processing workflow. Overall, this review provides an introductory guide
to modern research on artificial intelligence (AI) for de novo molecular
structure modeling and future directions in this emerging field.


  We investigate, through a fractional mathematical model, the effects of
physical distance on the SARS-CoV-2 virus transmission. Two controls are
considered in our model for eradication of the spread of COVID-19: media
education, through campaigns explaining the importance of social distancing,
use of face masks, etc., towards all population, while the second one is
quarantine social isolation of the exposed individuals. A general fractional
order optimal control problem, and associated optimality conditions of
Pontryagin type, are discussed, with the goal to minimize the number of
susceptible and infected while maximizing the number of recovered. The
extremals are then numerically obtained.


  Radiomic models have been shown to outperform clinical data for outcome
prediction in glioblastoma (GBM). However, clinical implementation is limited
by lack of parameters standardization. We aimed to compare nine machine
learning classifiers, with different optimization parameters, to predict
overall survival (OS), isocitrate dehydrogenase (IDH) mutation,
O-6-methylguanine-DNA-methyltransferase (MGMT) promoter methylation, epidermal
growth factor receptor (EGFR) VII amplification and Ki-67 expression in GBM
patients, based on radiomic features from conventional and advanced MR. 156
adult patients with pathologic diagnosis of GBM were included. Three tumoral
regions were analyzed: contrast-enhancing tumor, necrosis and non-enhancing
tumor, selected by manual segmentation. Radiomic features were extracted with a
custom version of Pyradiomics, and selected through Boruta algorithm. A Grid
Search algorithm was applied when computing 4 times K-fold cross validation
(K=10) to get the highest mean and lowest spread of accuracy. Once optimal
parameters were identified, model performances were assessed in terms of Area
Under The Curve-Receiver Operating Characteristics (AUC-ROC). Metaheuristic and
ensemble classifiers showed the best performance across tasks. xGB obtained
maximum accuracy for OS (74.5%), AB for IDH mutation (88%), MGMT methylation
(71,7%), Ki-67 expression (86,6%), and EGFR amplification (81,6%). Best
performing features shed light on possible correlations between MR and tumor
histology.


  A model-based signal processing framework is proposed for pandemic trend
forecasting and control, by using non-pharmaceutical interventions (NPI) at
regional and country levels worldwide. The control objective is to prescribe
quantifiable NPI strategies at different levels of stringency, which balance
between human factors (such as new cases and death rates) and cost of
intervention per region/country. Due to infrastructural disparities and
differences in priorities of regions and countries, strategists are given the
flexibility to weight between different NPIs and to select the desired balance
between the human factor and overall NPI cost.
  The proposed framework is based on a \textit{finite-horizon optimal control}
(FHOC) formulation of the bi-objective problem and the FHOC is numerically
solved by using an ad hoc \textit{extended Kalman filtering/smoothing}
framework for optimal NPI estimation and pandemic trend forecasting. The
algorithm enables strategists to select the desired balance between the human
factor and NPI cost with a set of weights and parameters. The parameters of the
model are partially selected by epidemiological facts from COVID-19 studies,
and partially trained by using machine learning techniques. The developed
algorithm is applied on ground truth data from the Oxford COVID-19 Government
Response Tracker project, which has categorized and quantified the regional
responses to the pandemic for more than 300 countries and regions worldwide,
since January 2020. The dataset was used for NPI-based prediction and
prescription during the XPRIZE Pandemic Response Challenge.


  An astonishingly diverse biomolecular circuitry orchestrates the functioning
machinery underlying every living cell. These biomolecules and their circuits
have been engineered not only for various industrial applications but also to
perform other atypical functions that they were not evolved for - including
computation. Various kinds of computational challenges, such as solving
NP-complete problems with many variables, logical computation, neural network
operations, and cryptography, have all been attempted through this
unconventional computing paradigm. In this review, we highlight key experiments
across three different eras of molecular computation, beginning with molecular
solutions, transitioning to logic circuits and ultimately, more complex
molecular networks. We also discuss a variety of applications of molecular
computation, from solving NP-hard problems to self-assembled nanostructures for
delivering molecules, and provide a glimpse into the exciting potential that
molecular computing holds for the future.


  In recent years, several global events have severely disrupted economies and
social structures, undermining confidence in the resilience of modern
societies. While empirical evidence on the dynamics and drivers of past
societal collapse is mounting, a process-based understanding of these dynamics
is still in its infancy. Here we aim to identify and illustrate the underlying
drivers of such societal instability or even collapse. The inspiration for this
work is Joseph Tainter's theory of the "collapse of complex societies", which
postulates that the complexity of societies increases as they solve problems,
leading to diminishing returns on complexity investments, and ultimately to
collapse. In this work we have abstracted this theory into a low-dimensional
and stylised model of two classes of networked agents, hereafter referred to as
"laborers" and "administrators". We numerically modeled the dynamics of
societal complexity, measured as the fraction of "administrators", which is
assumed to affect the productivity of connected energy-producing "laborers". We
show that collapse becomes increasingly likely as the complexity of the model
society continuously increases in response to external stresses that emulate
Tainter's abstract notion of problems that societies must solve. We also
provide an analytical approximation of the system's dominant dynamics, which
matches well with the numerical experiments, and use it to study the influence
on network link density, social mobility and productivity. Our work advances
the understanding of social-ecological collapse and illustrates its potentially
direct link to an ever-increasing societal complexity in response to external
shocks or stresses via a self-reinforcing feedback.


  High-throughput RNA-sequencing (RNA-seq) technologies are powerful tools for
understanding cellular state. Often it is of interest to quantify and summarize
changes in cell state that occur between experimental or biological conditions.
Differential expression is typically assessed using univariate tests to measure
gene-wise shifts in expression. However, these methods largely ignore changes
in transcriptional correlation. Furthermore, there is a need to identify the
low-dimensional structure of the gene expression shift to identify collections
of genes that change between conditions. Here, we propose contrastive latent
variable models designed for count data to create a richer portrait of
differential expression in sequencing data. These models disentangle the
sources of transcriptional variation in different conditions, in the context of
an explicit model of variation at baseline. Moreover, we develop a model-based
hypothesis testing framework that can test for global and gene subset-specific
changes in expression. We test our model through extensive simulations and
analyses with count-based gene expression data from perturbation and
observational sequencing experiments. We find that our methods can effectively
summarize and quantify complex transcriptional changes in case-control
experimental sequencing data.


  Protein fold classification is a classic problem in structural biology and
bioinformatics. We approach this problem using persistent homology. In
particular, we use alpha shape filtrations to compare a topological
representation of the data with a different representation that makes use of
knot-theoretic ideas. We use the statistical method of Angle-based Joint and
Individual Variation Explained (AJIVE) to understand similarities and
differences between these representations.


  The COVID-19 pandemic has proved to be one of the most disruptive public
health emergencies in recent memory. Among non-pharmaceutical interventions,
social distancing and lockdown measures are some of the most common tools
employed by governments around the world to combat the disease. While
mathematical models of COVID-19 are ubiquitous, few have leveraged network
theory in a general way to explain the mechanics of social distancing. In this
paper, we build on existing network models for heterogeneous, clustered
networks with random link activation/deletion dynamics to put forth realistic
mechanisms of social distancing using piecewise constant activation/deletion
rates. We find our models are capable of rich qualitative behavior, and offer
meaningful insight with relatively few intervention parameters. In particular,
we find that the severity of social distancing interventions and when they
begin have more impact than how long it takes for the interventions to take
full effect.


  The outbreak of COVID-19 in 2020 has led to a surge in interest in the
mathematical modeling of infectious diseases. Such models are usually defined
as compartmental models, in which the population under study is divided into
compartments based on qualitative characteristics, with different assumptions
about the nature and rate of transfer across compartments. Though most commonly
formulated as ordinary differential equation (ODE) models, in which the
compartments depend only on time, recent works have also focused on partial
differential equation (PDE) models, incorporating the variation of an epidemic
in space. Such research on PDE models within a Susceptible, Infected, Exposed,
Recovered, and Deceased (SEIRD) framework has led to promising results in
reproducing COVID-19 contagion dynamics. In this paper, we assess the
robustness of this modeling framework by considering different geometries over
more extended periods than in other similar studies. We first validate our code
by reproducing previously shown results for Lombardy, Italy. We then focus on
the U.S. state of Georgia and on the Brazilian state of Rio de Janeiro, one of
the most impacted areas in the world. Our results show good agreement with
real-world epidemiological data in both time and space for all regions across
major areas and across three different continents, suggesting that the modeling
approach is both valid and robust.


  The outbreak of COVID-19 led to a record-breaking race to develop a vaccine.
However, the limited vaccine capacity creates another massive challenge: how to
distribute vaccines to mitigate the near-end impact of the pandemic? In the
United States in particular, the new Biden administration is launching mass
vaccination sites across the country, raising the obvious question of where to
locate these clinics to maximize the effectiveness of the vaccination campaign.
This paper tackles this question with a novel data-driven approach to optimize
COVID-19 vaccine distribution. We first augment a state-of-the-art
epidemiological model, called DELPHI, to capture the effects of vaccinations
and the variability in mortality rates across age groups. We then integrate
this predictive model into a prescriptive model to optimize the location of
vaccination sites and subsequent vaccine allocation. The model is formulated as
a bilinear, non-convex optimization model. To solve it, we propose a coordinate
descent algorithm that iterates between optimizing vaccine distribution and
simulating the dynamics of the pandemic. As compared to benchmarks based on
demographic and epidemiological information, the proposed optimization approach
increases the effectiveness of the vaccination campaign by an estimated $20\%$,
saving an extra $4000$ extra lives in the United States over a three-month
period. The proposed solution achieves critical fairness objectives -- by
reducing the death toll of the pandemic in several states without hurting
others -- and is highly robust to uncertainties and forecast errors -- by
achieving similar benefits under a vast range of perturbations.


  We propose a stochastic SIR model, specified as a system of stochastic
differential equations, to analyse the data of the Italian COVID-19 epidemic,
taking also into account the under-detection of infected and recovered
individuals in the population. We find that a correct assessment of the amount
of under-detection is important to obtain reliable estimates of the critical
model parameters. Moreover, a single SIR model over the whole epidemic period
is unable to correctly describe the behaviour of the pandemic. Then, the
adaptation of the model in every time-interval between relevant government
decrees that implement contagion mitigation measures, provides short-term
predictions and a continuously updated assessment of the basic reproduction
number.


  As 2021 dawns, the COVID-19 pandemic is still raging strongly as vaccines
finally appear and hopes for a return to normalcy start to materialize. There
is much to be learned from the pandemic's first year data that will likely
remain applicable to future epidemics and possible pandemics. With only minor
variants in virus strain, countries across the globe have suffered roughly the
same pandemic by first glance, yet few locations exhibit the same patterns of
viral spread, growth, and control as the state of Hawai'i. In this paper, we
examine the data and compare the COVID-19 spread statistics between the
counties of Hawai'i as well as examine several locations with similar
properties to Hawai'i.


  We develop a model-independent reduction method of chemical reaction systems
based on the stoichiometry, which determines their network topology. A
subnetwork can be eliminated systematically to give a reduced system with fewer
degrees of freedom. This subnetwork removal is accompanied by rewiring of the
network, which is prescribed by the Schur complement of the stoichiometric
matrix. Using homology and cohomology groups to characterize the topology of
chemical reaction networks, we can track the changes of the network topology
induced by the reduction through the changes in those groups. We prove that,
when certain topological conditions are met, the steady-state chemical
concentrations and reaction rates of the reduced system are ensured to be the
same as those of the original system. This result holds regardless of the
modeling of the reactions, namely chemical kinetics, since the conditions only
involve topological information. This is advantageous because the details of
reaction kinetics and parameter values are difficult to identify in many
practical situations. The method allows us to reduce a reaction network while
preserving its original steady-state properties, thereby complex reaction
systems can be studied efficiently. We demonstrate the reduction method in
hypothetical networks and the central carbon metabolism of Escherichia coli.


  Precision medicine is a paradigm shift in healthcare relying heavily on
genomics data. However, the complexity of biological interactions, the large
number of genes as well as the lack of comparisons on the analysis of data,
remain a tremendous bottleneck regarding clinical adoption. In this paper, we
introduce a novel, automatic and unsupervised framework to discover
low-dimensional gene biomarkers. Our method is based on the LP-Stability
algorithm, a high dimensional center-based unsupervised clustering algorithm,
that offers modularity as concerns metric functions and scalability, while
being able to automatically determine the best number of clusters. Our
evaluation includes both mathematical and biological criteria. The recovered
signature is applied to a variety of biological tasks, including screening of
biological pathways and functions, and characterization relevance on tumor
types and subtypes. Quantitative comparisons among different distance metrics,
commonly used clustering methods and a referential gene signature used in the
literature, confirm state of the art performance of our approach. In
particular, our signature, that is based on 27 genes, reports at least $30$
times better mathematical significance (average Dunn's Index) and 25% better
biological significance (average Enrichment in Protein-Protein Interaction)
than those produced by other referential clustering methods. Finally, our
signature reports promising results on distinguishing immune inflammatory and
immune desert tumors, while reporting a high balanced accuracy of 92% on tumor
types classification and averaged balanced accuracy of 68% on tumor subtypes
classification, which represents, respectively 7% and 9% higher performance
compared to the referential signature.


  The study of epidemiological systems has generated deep interest in exploring
the dynamical complexity of common infectious diseases driven by seasonally
varying contact rates. Mathematical modeling and field observations have shown
that, under seasonal variation, the incidence rates of some endemic infectious
diseases fluctuate dramatically and the dynamics is often characterized by
chaotic oscillations in the absence of specific vaccination programs. In fact,
the existence of chaotic behavior has been precisely stated in the literature
as a noticeable feature in the dynamics of the classical
Susceptible-Infected-Recovered (SIR) seasonally forced epidemic model. However,
in the context of epidemiology, chaos is often regarded as an undesirable
phenomenon associated with the unpredictability of infectious diseases. As a
consequence, the problem of converting chaotic motions into regular motions
becomes particularly relevant. In this article, we consider the phase control
technique applied to the seasonally forced SIR epidemic model to suppress
chaos. Interestingly, this method of controlling chaos takes on a clear meaning
as a weak perturbation on a seasonal vaccination strategy. Numerical
simulations show that the phase difference between the two periodic forces -
contact rate and vaccination - plays a very important role in controlling
chaos.


  Polygenic risk score (PRS) analysis is a powerful method been used to
estimate an individual's genetic risk towards targeted traits. PRS analysis
could be used to obtain evidence of a genetic effect beyond Genome-Wide
Association Studies (GWAS) results i.e. when there are no significant markers.
PRS analysis has been widely applied to investigate the genetic basis of
several traits including rare diseases. However, the accuracy of PRS analysis
depends on the genomic data of the underlying population. For instance, several
studies showed that obtaining higher prediction power of PRS analysis is
challenging for non-Europeans. In this manuscript, we reviewed the conventional
PRS methods and their application to sub-saharan Africa communities. We
concluded that the limiting factor of applying PRS analysis to sub-saharan
populations is the lack of sufficient GWAS data. Also, we recommended
developing African-specific PRS tools


  This paper focuses on a core task in computational sustainability and
statistical ecology: species distribution modeling (SDM). In SDM, the
occurrence pattern of a species on a landscape is predicted by environmental
features based on observations at a set of locations. At first, SDM may appear
to be a binary classification problem, and one might be inclined to employ
classic tools (e.g., logistic regression, support vector machines, neural
networks) to tackle it. However, wildlife surveys introduce structured noise
(especially under-counting) in the species observations. If unaccounted for,
these observation errors systematically bias SDMs. To address the unique
challenges of SDM, this paper proposes a framework called StatEcoNet.
Specifically, this work employs a graphical generative model in statistical
ecology to serve as the skeleton of the proposed computational framework and
carefully integrates neural networks under the framework. The advantages of
StatEcoNet over related approaches are demonstrated on simulated datasets as
well as bird species data. Since SDMs are critical tools for ecological science
and natural resource management, StatEcoNet may offer boosted computational and
analytical powers to a wide range of applications that have significant social
impacts, e.g., the study and conservation of threatened species.


  As the worldwide population gets increasingly aged, in-home telemedicine and
mobile-health solutions represent promising services to promote active and
independent aging and to contribute to a paradigm shift towards patient-centric
healthcare. In this work, we present ACTA (Advanced Cognitive Training for
Aging), a prototype mobile-health solution to provide advanced cognitive
training for senior citizens with mild cognitive impairments. We disclose here
the conceptualization of ACTA as the integration of two promising
rehabilitation strategies: the "Nudge theory", from the cognitive domain, and
the neurofeedback, from the neuroscience domain. Moreover, in ACTA we exploit
the most advanced machine learning techniques to deliver customized and fully
adaptive support to the elderly, while training in an ecological environment.
ACTA represents the next-step beyond SENIOR, an earlier mobile-health project
for cognitive training based on Nudge theory, currently ongoing in Lombardy
Region. Beyond SENIOR, ACTA represents a highly-usable, accessible, low-cost,
new-generation mobile-health solution to promote independent aging and
effective motor-cognitive training support, while empowering the elderly in
their own aging.


  Despite the progress in medical care, combined population-wide interventions
(such as physical distancing, testing and contact tracing) are still crucial to
manage the SARS-CoV-2 pandemic, aggravated by the emergence of new highly
transmissible variants. We combine the compartmental SIDARTHE model, predicting
the course of COVID-19 infections, with a new data-based model that projects
new cases onto casualties and healthcare system costs. Based on the Italian
case study, we outline several scenarios: mass vaccination campaigns with
different paces, different transmission rates due to new variants, and
different enforced countermeasures, including the alternation of opening and
closure phases. Our results demonstrate that non-pharmaceutical interventions
(NPIs) have a higher impact on the epidemic evolution than vaccination, which
advocates for the need to keep containment measures in place throughout the
vaccination campaign. We also show that, if intermittent open-close strategies
are adopted, deaths and healthcare system costs can be drastically reduced,
without any aggravation of socioeconomic losses, as long as one has the
foresight to start with a closing phase rather than an opening one.


  Biological processes at the cellular level are stochastic in nature, and the
immune response system is no different. Therefore, models that attempt to
explain this system need to also incorporate noise or fluctuations that can
account for the observed variability. In this work, a stochastic model of the
immune response system is presented in terms of the dynamics of the T cells and
the virus particles. Making use of the Green's function and the Wilemski-Fixman
approximation, this model is then solved to obtain the analytical expression
for the joint probability density function of these variables in the early and
late stages of infection. This is then also used to calculate the average level
of virus particles in the system. Upon comparing the theoretically predicted
average virus levels to those of COVID-19 patients, it is hypothesized that the
long lived dynamics that are characteristic of such viral infections are due to
the long range correlations in the temporal fluctuations of the virions. This
model therefore provides an insight into the effects of noise on viral
dynamics.


  We are interested in describing the infected size of the SIS Epidemic model
using Birth-Death Markov process. The Susceptible-Infected-Susceptible (SIS)
model is defined within a population of constant size $M$; the size is kept
constant by replacing each death with a newborn healthy individual. The life
span of each individual in the population is modelled by an exponential
distribution with parameter $\alpha$; and the disease spreads within the
population is modelled by a Poisson process with a rate $\lambda_{I}$.
$\lambda_{I}=\beta I(1-\frac{I}{M}) $ is similar to the instantaneous rate in
the logistic population growth model. The analysis is focused on the disease
outbreak, where the reproduction number $R=\frac{\beta}{\alpha} $ is greater
than one. As methodology, we use both numerical and analytical approaches. The
analysis relies on the stationary distribution for Birth and Death Markov
process. The numerical approach creates sample path simulations into order show
the infected size dynamics, and the relationship between infected size and $R$.
As $M$ becomes large, some stable statistical characteristics of the infected
size distribution can be deduced. And the infected size is shown analytically
to follow a normal distribution with mean $(1-\frac{1}{R}) M$ and Variance
$\frac{M}{R} $.


  Recent work has highlighted roles for thermodynamic phase behavior in diverse
cellular processes. Proteins and nucleic acids can phase separate into
three-dimensional liquid droplets in the cytoplasm and nucleus and the plasma
membrane of animal cells appears tuned close to a two-dimensional liquid-liquid
critical point. In some examples, cytoplasmic proteins aggregate at plasma
membrane domains, forming structures such as the post-synaptic density and
diverse signaling clusters. Here we examine the physics of these surface
densities, employing minimal simulations of co-acervating polymers coupled to
an Ising membrane surface in conjunction with a complementary Landau theory. We
argue that these surface densities are a novel phase reminiscent of
pre-wetting, in which a molecularly thin three-dimensional liquid forms on a
usually solid surface. However, in surface densities the solid surface is
replaced by a membrane with an independent propensity to phase separate. We
show that proximity to criticality in the membrane dramatically increases the
parameter regime in which a pre-wetting-like transition occurs, leading to a
broad region where coexisting surface phases can form even when a bulk phase is
unstable. Our simulations naturally exhibit three surface phase coexistence
even though both the membrane and the polymer bulk can only display two phase
coexistence on their own. We argue that the physics of these surface densities
enables diverse functions seen in Eukaryotic cells.


  This paper is concerned with exploring the microscopic basis for the discrete
versions of the standard replicator equation and the adjusted replicator
equation. To this end, we introduce frequency-dependent selection -- as a
result of competition fashioned by game-theoretic consideration -- into the
Wright--Fisher process, a stochastic birth-death process. The process is
further considered to be active in a generation-wise nonoverlapping finite
population where individuals play a two-strategy bimatrix population game.
Subsequently, connections among the corresponding master equation, the
Fokker--Planck equation, and the Langevin equation are exploited to arrive at
the deterministic discrete replicator maps in the limit of infinite population
size.


  Therapeutics machine learning is an emerging field with incredible
opportunities for innovatiaon and impact. However, advancement in this field
requires formulation of meaningful learning tasks and careful curation of
datasets. Here, we introduce Therapeutics Data Commons (TDC), the first
unifying platform to systematically access and evaluate machine learning across
the entire range of therapeutics. To date, TDC includes 66 AI-ready datasets
spread across 22 learning tasks and spanning the discovery and development of
safe and effective medicines. TDC also provides an ecosystem of tools and
community resources, including 33 data functions and types of meaningful data
splits, 23 strategies for systematic model evaluation, 17 molecule generation
oracles, and 29 public leaderboards. All resources are integrated and
accessible via an open Python library. We carry out extensive experiments on
selected datasets, demonstrating that even the strongest algorithms fall short
of solving key therapeutics challenges, including real dataset distributional
shifts, multi-scale modeling of heterogeneous data, and robust generalization
to novel data points. We envision that TDC can facilitate algorithmic and
scientific advances and considerably accelerate machine-learning model
development, validation and transition into biomedical and clinical
implementation. TDC is an open-science initiative available at
https://tdcommons.ai.


  Many biological processes have been the source of inspiration for heuristic
methods that generate high-quality solutions to solve optimization and search
problems. This thesis presents an epigenetic technique for Evolutionary
Algorithms, inspired by the epigenetic regulation process, a mechanism to
better understand the ability of individuals to adapt and learn from the
environment. Epigenetic regulation comprises biological mechanisms by which
small molecules, also known as epigenetic tags, are attached to or removed from
a particular gene, affecting the phenotype. Five fundamental elements form the
basis of the designed technique: first, a metaphorical representation of
Epigenetic Tags as binary strings; second, a layer on chromosome top structure
used to bind the tags (the Epigenotype layer); third, a Marking Function to
add, remove, and modify tags; fourth, an Epigenetic Growing Function that acts
like an interpreter, or decoder of the tags located over the alleles, in such a
way that the phenotypic variations can be reflected when evaluating the
individuals; and fifth, a tags inheritance mechanism. A set of experiments are
performed for determining the applicability of the proposed approach.