  The development of coalescent theory paved the way to statistical inference
from population genetic data. In the genomic era, however, coalescent models
are limited due to the complexity of the underlying ancestral recombination
graph. The sequentially Markov coalescent (SMC) is a heuristic that enables the
modelling of complete genomes under the coalescent framework. While it empowers
the inference of detailed demographic history of a population from as few as
one diploid genome, current implementations of the SMC make unrealistic
assumptions about the homogeneity of the coalescent process along the genome,
ignoring the intrinsic spatial variability of parameters such as the
recombination rate. Here, I review the historical developments of SMC models
and discuss the evidence for parameter heterogeneity. I then survey approaches
to handle this heterogeneity, focusing on a recently developed extension of the
SMC.


  During the COVID-19 pandemic, many countries opted for strict public health
measures, including closing schools. They have now started relaxing some of
those restrictions. To avoid overwhelming health systems, predictions for the
number of new COVID-19 cases need to be considered when choosing a school
reopening strategy. Using computer simulation, we analyze different strategies
to reopen schools in the Sao Paulo Metropolitan Area, including the official
reopening plan. Our results indicate that reopening schools with all students
at once has a big impact on the number of new COVID-19 cases, which could cause
a collapse of the health system. On the other hand, a controlled school
reopening avoids the collapse of the health system, with a maximum ICU
occupancy between 48.8% and 97.8%, depending on how people attending school
follow sanitary measures. Postponing the schools' return date for after a
vaccine becomes available can save up to 37,753 lives just in the Sao Paulo
Metropolitan Area.


  Mass concentrations of birds, or lack of such, is a phenomenon of great
ecological and domestic significance. Apart from being and indicator for e.g.
food availability, ecological change and population size, it is also a source
of conflict between humans and birds. Moreover, massive gatherings or colonies
of birds also get the attention of the public -- either as a spectacular
phenomenon or as an unwelcome pest -- thereby forming the public perception of
birds and their abundance. In the context of the mass concentration of
bramblings (Fringilla montifringilla) in Sweden the winter 2019-2020, this work
reviews the literature on this striking phenomenon. Winter roosts are found to
amount to on the order of one million birds per hectare of roost area, but the
spread between reports is significant. Support for roosts of up to around 15
million birds was found, but much larger numbers are frequently recited in the
literature. It is argued that these larger numbers are the result of
overestimation or, in some cases, even completely unfounded (potentially
typos). While the difficulties related to the count of large numbers of birds
can explain this state, it is unfortunate that "high numbers" sometimes
displace proper numbers. Since incorrect data, and its persistence, may result
in that incorrect conclusions are drawn from new observations, this matter
deserves attention. As the Brambling is a well-studied species, the matter also
raises concerns regarding numbers for mass concentrations of other species. It
is recommended that very large numbers of birds should be recited and used with
care: underlying data and methods of the original sources should be
scrutinized. Analogously, reporters of large numbers of birds are advised to
describe and document counting methods. In particular, number estimates based
on flock volume and bird density should be avoided.


  Autocatalytic Sets are reaction networks theorised as networks at the basis
of life. Their main feature is the ability of spontaneously emerging and
self-reproducing. The Reflexively and Food-generated theory provides a formal
definition of Autocatalytic Sets in terms of graphs with peculiar topological
properties. This formalisation has been proved to be a powerful tool for the
study of the chemical networks underlying life, and it was able to identify
autocatalytic structures in real metabolic networks. However, the dynamical
behaviour of such networks has not been yet complitely clarified. In this work,
we present a first attempt to connect the topology of an Autocatalytic Set with
its dynamics. For this purpose, we represent Autocatalytic Sets in terms of
Chemical Reaction Networks, and we use the Chemical Reaction Network theory to
detect motifs in the networks'structure, that allow us to determine the
long-term behaviour of the system.


  Deep generative models are increasingly becoming integral parts of the in
silico molecule design pipeline and have dual goals of learning the chemical
and structural features that render candidate molecules viable while also being
flexible enough to generate novel designs. Specifically, Variational Auto
Encoders (VAEs) are generative models in which encoder-decoder network pairs
are trained to reconstruct training data distributions in such a way that the
latent space of the encoder network is smooth. Therefore, novel candidates can
be found by sampling from this latent space. However, the scope of
architectures and hyperparameters is vast and choosing the best combination for
in silico discovery has important implications for downstream success.
Therefore, it is important to develop a principled methodology for
distinguishing how well a given generative model is able to learn salient
molecular features. In this work, we propose a method for measuring how well
the latent space of deep generative models is able to encode structural and
chemical features of molecular datasets by correlating latent space metrics
with metrics from the field of topological data analysis (TDA). We apply our
evaluation methodology to a VAE trained on SMILES strings and show that 3D
topology information is consistently encoded throughout the latent space of the
model.


  Machine learning in drug discovery has been focused on virtual screening of
molecular libraries using discriminative models. Generative models are an
entirely different approach that learn to represent and optimize molecules in a
continuous latent space. These methods have been increasingly successful at
generating two dimensional molecules as SMILES strings and molecular graphs. In
this work, we describe deep generative models of three dimensional molecular
structures using atomic density grids and a novel fitting algorithm for
converting continuous grids to discrete molecular structures. Our models
jointly represent drug-like molecules and their conformations in a latent space
that can be explored through interpolation. We are also able to sample diverse
sets of molecules based on a given input compound and increase the probability
of creating valid, drug-like molecules.


  The protein-ligand residence time, tau, influences molecular function in
biological networks and has been recognized as an important determinant of drug
efficacy. To predict tau, computational methods must overcome the problem that
tau often exceeds the timescales accessible to conventional molecular dynamics
(MD) simulation. Here, we apply the tau-Random Acceleration Molecular Dynamics
(tauRAMD) method to a set of kinetically characterized complexes of T4 lysozyme
mutants with engineered binding cavities. tauRAMD yields relative ligand
dissociation rates in good accordance with experiments and thereby allows a
comprehensive characterization of the ligand egress routes and determinants of
tau. Although ligand dissociation by multiple egress routes is observed, we
find that egress via the predominant route determines the value of tau. We also
find that the presence of metastable states along egress pathways slows down
protein-ligand dissociation. These physical insights could be exploited in the
rational optimization of the kinetic properties of drug candidates.


  Within-host models of COVID-19 infection dynamics enable the merits of
different forms of antiviral therapy to be assessed in individual patients. A
stochastic agent-based model of COVID-19 intracellular dynamics is introduced
here, that incorporates essential steps of the viral life cycle targeted by
treatment options. Integration of model predictions with an intercellular ODE
model of within-host infection dynamics, fitted to patient data, generates a
generic profile of disease progression in patients that have recovered in the
absence of treatment. This is contrasted with the profiles obtained after
variation of model parameters pertinent to the immune response, such as
effector cell and antibody proliferation rates, mimicking disease progression
in immunocompromised patients. These profiles are then compared with disease
progression in the presence of antiviral and convalescent plasma therapy
against COVID-19 infections. The model reveals that using both therapies in
combination can be very effective in reducing the length of infection, but
these synergistic effects decline with a delayed treatment start. Conversely,
early treatment with either therapy alone can actually increase the duration of
infection, with infectious virions still present after the decline of other
markers of infection. This suggests that usage of these treatments should
remain carefully controlled in a clinical environment.


  In the present paper, we model the cumulative number of persons reported to
be infected by the SARS-CoV-2 virus, in a country or a region, by a sum of
logistic functions. For a given logistic function, using Eulerian numbers, we
find the zeros of its successive derivatives and their relationship with the
saturation level of this function. In a given time series, having potentially
the logistic trend, we use its second differences to determine points
corresponding to these zeros. To estimate the parameters of the approximating
logistic function, we define and use logistic wavelets. Then we apply the
theory to the cases of SARS-CoV-2 infections in the United States and the
United Kingdom.


  We propose a mathematical model for the spread of Japanese encephalitis, with
emphasis on environmental effects on the aquatic phase of mosquitoes. The model
is shown to be biologically well-posed and to have a biologically and
ecologically meaningful disease free equilibrium point. Local stability is
analyzed in terms of the basic reproduction number and numerical simulations
presented and discussed.


  Motivation The genotype assignment problem consists of predicting, from the
genotype of an individual, which of a known set of populations it originated
from. The problem arises in a variety of contexts, including wildlife
forensics, invasive species detection and biodiversity monitoring. Existing
approaches perform well under ideal conditions but are sensitive to a variety
of common violations of the assumptions they rely on. Results In this article,
we introduce Mycorrhiza, a machine learning approach for the genotype
assignment problem. Our algorithm makes use of phylogenetic networks to
engineer features that encode the evolutionary relationships among samples.
Those features are then used as input to a Random Forests classifier. The
classification accuracy was assessed on multiple published empirical SNP,
microsatellite or consensus sequence datasets with wide ranges of size,
geographical distribution and population structure and on simulated datasets.
It compared favorably against widely used assessment tests or mixture analysis
methods such as STRUCTURE and Admixture, and against another machine-learning
based approach using principal component analysis for dimensionality reduction.
Mycorrhiza yields particularly significant gains on datasets with a large
average fixation index (FST) or deviation from the Hardy-Weinberg equilibrium.
Moreover, the phylogenetic network approach estimates mixture proportions with
good accuracy.


  Timely prediction of the COVID-19 progression is not possible without a
comprehensive understanding of environmental factors that may affect the
infection transmissibility. Studies addressing parameters that may influence
COVID-19 progression relied on either the total numbers of detected cases and
similar proxies and/or a small number of analyzed factors, including analysis
of regions that display a narrow range of these parameters. We here apply a
novel approach, exploiting widespread growth regimes in COVID-19 detected case
counts. By applying nonlinear dynamics methods to the exponential regime, we
extract basic reproductive number R0 (i.e., the measure of COVID-19 inherent
biological transmissibility), applying to the completely naive population in
the absence of social distancing, for 118 different countries. We then use
bioinformatics methods to systematically collect data on a large number of
demographics and weather parameters from these countries, and seek their
correlations with the rate of COVID-19 spread. In addition to some of the
already reported tendencies, we show a number of both novel results and those
that help settle existing disputes: the absence of dependence on wind speed and
air pressure, negative correlation with precipitation; significant positive
correlation with society development level (human development index)
irrespective of testing policies, and percent of the urban population, but an
absence of correlation with population density per se. We find a strong
positive correlation of transmissibility on alcohol consumption, and the
absence of correlation on refugee numbers, contrary to some widespread beliefs.
Significant tendencies with health-related factors are reported, including a
detailed analysis of the blood type group showing consistent tendencies on Rh
factor, and a strong positive correlation of transmissibility with cholesterol
levels.


  Lamina-associated domains (LADs) cover a large part of the human genome and
are thought to play a major role in shaping the nuclear architectural
landscape. Here, we perform polymer simulations, microscopy and mass
spectrometry to dissect the roles played by heterochromatin- and
lamina-mediated interactions in nuclear organisation. Our model explains the
conventional organisation of heterochromatin and euchromatin in growing cells
and the pathological organisation found in oncogene-induced senescence and
progeria. We show that the experimentally observed changes in the locality of
contacts in senescent and progeroid cells can be explained as arising due to
phase transitions in the system. Within our simulations LADs are highly
stochastic, as in experiments. Our model suggests that, once established, the
senescent phenotype should be metastable even if lamina-mediated interactions
were reinstated. Overall, our simulations uncover a generic physical mechanism
that can regulate heterochromatin segregation and LAD formation in a wide range
of mammalian nuclei.


  GNNs and chemical fingerprints are the predominant approaches to representing
molecules for property prediction. However, in NLP, transformers have become
the de-facto standard for representation learning thanks to their strong
downstream task transfer. In parallel, the software ecosystem around
transformers is maturing rapidly, with libraries like HuggingFace and BertViz
enabling streamlined training and introspection. In this work, we make one of
the first attempts to systematically evaluate transformers on molecular
property prediction tasks via our ChemBERTa model. ChemBERTa scales well with
pretraining dataset size, offering competitive downstream performance on
MoleculeNet and useful attention-based visualization modalities. Our results
suggest that transformers offer a promising avenue of future work for molecular
representation learning and property prediction. To facilitate these efforts,
we release a curated dataset of 77M SMILES from PubChem suitable for
large-scale self-supervised pretraining.


  One of the most computationally intensive tasks in computational biology is
de novo genome assembly, the decoding of the sequence of an unknown genome from
redundant and erroneous short sequences. A common assembly paradigm identifies
overlapping sequences, simplifies their layout, and creates consensus. Despite
many algorithms developed in the literature, the efficient assembly of large
genomes is still an open problem. In this work, we introduce new
distributed-memory parallel algorithms for overlap detection and layout
simplification steps of de novo genome assembly, and implement them in the
diBELLA 2D pipeline. Our distributed memory algorithms for both overlap
detection and layout simplification are based on linear-algebra operations over
semirings using 2D distributed sparse matrices. Our layout step consists of
performing a transitive reduction from the overlap graph to a string graph. We
provide a detailed communication analysis of the main stages of our new
algorithms. diBELLA 2D achieves near linear scaling with over 80% parallel
efficiency for the human genome, reducing the runtime for overlap detection by
1.2-1.3x for the human genome and 1.5-1.9x for C. elegans compared to the
state-of-the-art. Our transitive reduction algorithm outperforms an existing
distributed-memory implementation by 10.5-13.3x for the human genome and 18-29x
for the C. elegans. Our work paves the way for efficient de novo assembly of
large genomes using long reads in distributed memory.


  We develop a theoretical description of the topological disentanglement
occurring when torus knots reach the ends of a semi-flexible polymer under
tension. These include decays into simpler knots and total unknotting. The
minimal number of crossings and the minimal knot contour length are the
topological invariants playing a key role in the model. The crossings behave as
particles diffusing along the chain and the application of appropriate boundary
conditions at the ends of the chain accounts for the knot disentanglement.
Starting from the number of particles and their positions, suitable rules allow
reconstructing the type and location of the knot moving on the chain. Our
theory is extensively benchmarked with corresponding Molecular Dynamics
simulations and the results show a remarkable agreement between the simulations
and the theoretical predictions of the model.


  We present a symbolic algorithmic approach that allows to compute invariant
manifolds and corresponding reduced systems for differential equations modeling
biological networks which comprise chemical reaction networks for cellular
biochemistry, and compartmental models for pharmacology, epidemiology and
ecology. Multiple time scales of a given network are obtained by scaling, based
on tropical geometry. Our reduction is mathematically justified within a
singular perturbation setting. The existence of invariant manifolds is subject
to hyperbolicity conditions, for which we propose an algorithmic test based on
Hurwitz criteria. We finally obtain a sequence of nested invariant manifolds
and respective reduced systems on those manifolds. Our theoretical results are
generally accompanied by rigorous algorithmic descriptions suitable for direct
implementation based on existing off-the-shelf software systems, specifically
symbolic computation libraries and Satisfiability Modulo Theories solvers. We
present computational examples taken from the well-known BioModels database
using our own prototypical implementations.


  When the hydration shell of a protein is filled with at least 0.6 gram of
water per gram of protein, a significant anti-correlation between the
vibrational free energy and the potential energy of energy-minimized conformers
is observed. This means that low potential energy, well-hydrated, protein
conformers tend to be more rigid than high-energy ones. On the other hand, in
the case of CASP target 624, when its hydration shell is filled, a significant
average energy gap is observed between the crystal structure and the best
conformers proposed during the prediction experiment, strongly suggesting that
including explicit water molecules may help identifying unlikely conformers
among good-looking ones.


  Purpose: This paper proposes a methodology and a computational tool to study
the COVID-19 pandemic throughout the world and to perform a trend analysis to
assess its local dynamics.
  Methods: Mathematical functions are employed to describe the number of cases
and demises in each region and to predict their final numbers, as well as the
dates of maximum daily occurrences and the local stabilization date. The model
parameters are calibrated using a computational methodology for numerical
optimization. Trend analyses are run, allowing to assess the effects of public
policies. Easy to interpret metrics over the quality of the fitted curves are
provided. Country-wise data from the European Centre for Disease Prevention and
Control (ECDC) concerning the daily number of cases and demises around the
world are used, as well as detailed data from Johns Hopkins University and from
the Brasil.io project describing individually the occurrences in United States
counties and in Brazilian states and cities, respectively. U. S. and Brazil
were chosen for a more detailed analysis because they are the current foci of
the pandemic.
  Results: Illustrative results for different countries, U. S. counties and
Brazilian states and cities are presented and discussed.
  Conclusion: The main contributions of this work lie in (i) a straightforward
model of the curves to represent the data, which allows automation of the
process without requiring interventions from experts; (ii) an innovative
approach for trend analysis, whose results provide important information to
support authorities in their decision-making process; and (iii) the developed
computational tool, which is freely available and allows the user to quickly
update the COVID-19 analyses and forecasts for any country, United States
county or Brazilian state or city present in the periodic reports from the
authorities.


  With the advent of high-throughput sequencing technologies, the fields of
immunogenomics and adaptive immune receptor repertoire research are facing both
opportunities and challenges. Adaptive immune receptor repertoire sequencing
(AIRR-seq) has become an increasingly important tool to characterize T and B
cell responses in settings of interest. However, the majority of AIRR-seq
studies conducted so far were performed in individuals of European ancestry,
restricting the ability to identify variation in human adaptive immune
responses across populations and limiting their applications. As AIRR-seq
studies depend on the ability to assign VDJ sequence reads to the correct
germline gene segments, efforts to characterize the genomic loci that encode
adaptive immune receptor genes in different populations are urgently needed.
The availability of comprehensive germline gene databases and further
applications of AIRR-seq studies to individuals of non-European ancestry will
substantially enhance our understanding of human adaptive immune responses,
promote the development of effective diagnostics and treatments, and eventually
advance precision medicine.


  We present an SI epidemic model whereby a continuous variable captures
variability in proliferative potential and resistance to infection among
susceptibles. The occurrence of heritable, spontaneous changes in these
phenotype and the presence of a fitness trade-off between resistance to
infection and proliferative potential are incorporated into the model. The
model comprises an ODE for the number of infected individuals that is coupled
with a partial integrodifferential equation for the population density of
susceptibles through an integral term. The expression for the basic
reproduction number $\mathcal{R}_0$ is derived, the disease-free and endemic
equilibrium of the model are characterised and a threshold theorem is proved.
Analytical results are integrated with numerical simulations of a calibrated
version of the model based on the results of artificial selection experiments
in a host-parasite system. The results of our mathematical study disentangle
the impact of different evolutionary parameters on the spread of infectious
diseases and the consequent phenotypic adaption of susceptible individuals. In
particular, these results provide a theoretical basis for the observation that
infectious diseases exerting stronger selective pressures on susceptibles and
being characterised by higher infection rates are more likely to spread.
Moreover, our results indicate that spontaneous phenotypic changes in
proliferative potential and resistance to infection can either promote or
prevent the spread of diseases depending on the strength of selection acting on
susceptible individuals prior to infection. Finally, we demonstrate that, when
an endemic equilibrium is established, higher levels of resistance to infection
and lower degrees of phenotypic heterogeneity are to be expected in the
presence of infections which are characterised by lower rates of death and
exert stronger selective pressures.


  We study the relative importance of two key control measures for epidemic
spreading: endogenous social self-distancing and exogenous imposed quarantine.
We use the framework of adaptive networks, moment-closure, and ordinary
differential equations (ODEs) to introduce several novel models based upon
susceptible-infected-recovered (SIR) dynamics. First, we compare
computationally expensive, adaptive network simulations with their
corresponding computationally highly efficient ODE equivalents and find
excellent agreement. Second, we discover that there exists a relatively simple
critical curve in parameter space for the epidemic threshold, which strongly
suggests that there is a mutual compensation effect between the two mitigation
strategies: as long as social distancing and quarantine measures are both
sufficiently strong, large outbreaks are prevented. Third, we study the total
number of infected and the maximum peak during large outbreaks using a
combination of analytical estimates and numerical simulations. Also for large
outbreaks we find a similar compensation effect as for the epidemic threshold.
This suggests that if there is little incentive for social distancing within a
population, drastic quarantining is required, and vice versa. Both pure
scenarios are unrealistic in practice. Our models show that only a combination
of measures is likely to succeed to control epidemic spreading. Fourth, we
analytically compute an upper bound for the total number of infected on
adaptive networks, using integral estimates in combination with the
moment-closure approximation on the level of an observable. This is a
methodological innovation. Our method allows us to elegantly and quickly check
and cross-validate various conjectures about the relevance of different network
control measures.


  The Poisson-Boltzmann (PB) implicit solvent model is a popular framework for
studying the electrostatics of biomolecules immersed in water with dissolved
salt. In this model the dielectric interface between the biomolecule and
solvent is often taken to be the molecular surface or solvent-excluded surface
(SES), and the quality of the SES triangulation is critical in boundary element
simulations of the PB model. In this work we compare the MSMS and NanoShaper
surface triangulation codes for a set of 38 biomolecules. While MSMS produces
triangles of exceedingly small area and large aspect ratio, the two codes yield
comparable values for the SES surface area and electrostatic solvation energy,
where the latter calculations were performed using the treecode-accelerated
boundary integral (TABI) PB solver. However we found that NanoShaper is more
efficient and reliable than MSMS, especially when parameters are set to produce
highly resolved triangulations.


  Machine learning methods are increasingly employed to address challenges
faced by biologists. One area that will greatly benefit from this
cross-pollination is the problem of biological sequence design, which has
massive potential for therapeutic applications. However, significant
inefficiencies remain in communication between these fields which result in
biologists finding the progress in machine learning inaccessible, and hinder
machine learning scientists from contributing to impactful problems in
bioengineering. Sequence design can be seen as a search process on a discrete,
high-dimensional space, where each sequence is associated with a function. This
sequence-to-function map is known as a "Fitness Landscape". Designing a
sequence with a particular function is hence a matter of "discovering" such a
(often rare) sequence within this space. Today we can build predictive models
with good interpolation ability due to impressive progress in the synthesis and
testing of biological sequences in large numbers, which enables model training
and validation. However, it often remains a challenge to find useful sequences
with the properties that we like using these models. In particular, in this
primer we highlight that algorithms for experimental design, what we call
"exploration strategies", are a related, yet distinct problem from building
good models of sequence-to-function maps. We review advances and insights from
current literature -- by no means a complete treatment -- while highlighting
desirable features of optimal model-guided exploration, and cover potential
pitfalls drawn from our own experience. This primer can serve as a starting
point for researchers from different domains that are interested in the problem
of searching a sequence space with a model, but are perhaps unaware of
approaches that originate outside their field.


  We propose a model for demixing of two species by assuming a
density-dependent effective diffusion coefficient of the particles. Both sorts
of microswimmers diffuse as active overdamped Brownian particles with a noise
intensity that is determined by the surrounding density of the respective other
species within a sensing radius $r_s$. A higher concentration of the first
(second) sort will enlarge the diffusion and, in consequence, the intensity of
the noise experienced by the second (first) sort. Numerical and analytical
investigations of steady states of the macroscopic equations prove the demixing
of particles due to this reciprocally concentration-dependent diffusivity. An
ambiguity of the numerical integration scheme for the purely local model
($r_s\to 0$) is resolved by considering nonvanishing sensing radii in a
nonlocal model with $r_s > 0$.


  Epidemic models are useful tools in the fight against infectious diseases, as
they allow policy makers to test and compare various strategies to limit
disease transmission while mitigating collateral damage on the economy.
Epidemic models that are more faithful to the microscopic details of disease
transmission can offer more reliable projections, which in turn can lead to
more reliable control strategies. For example, many epidemic models describe
disease progression via a series of artificial 'stages' or 'compartments' (e.g.
exposed, activated, infectious, etc.) but an epidemic model that explicitly
tracks time since infection (TSI) can provide a more precise description. At
present, epidemic models with 'compartments' are more common than TSI models ,
largely due to higher computational cost and complexity typically associated
with TSI models. Here, however, we show that with the right discretization
scheme a TSI model is not much more difficult to solve than a comparment model
with three or four 'stages' for the infected class. We also provide a new
perspective for adding 'stages' to a TSI model in a way that decouples the
disease transmission dynamics from the residence time distributions at each
stage. These results are also generalized for age-structured TSI models in an
appendix. Finally, as proof-of-principle for the efficiency of the proposed
numerical methods, we provide calculations for optimal epidemic control by
non-pharmaceutical intervention. Many of the tools described in this report are
available through the software package 'pyross'


  In this paper, I study epidemic diffusion in a generalized spatial SEIRD
model, where individuals are initially connected in a social or geographical
network. As the virus spreads in the network, the structure of interactions
between people may endogenously change over time, due to quarantining measures
and/or spatial-distancing policies. I explore via simulations the dynamic
properties of the co-evolutionary process dynamically linking disease diffusion
and network properties. Results suggest that, in order to predict how epidemic
phenomena evolve in networked populations, it is not enough to focus on the
properties of initial interaction structures. Indeed, the co-evolution of
network structures and compartment shares strongly shape the process of
epidemic diffusion, especially in terms of its speed. Furthermore, I show that
the timing and features of spatial-distancing policies may dramatically
influence their effectiveness.


  Liquid-liquid phase separation is now recognized as a common mechanism for
regulating enzyme activity in cells. Insights from studies in cells are
complemented by in vitro studies aimed at developing better understanding of
mechanisms underlying such control. These mechanisms are often based on the
influence of LLPS on the physicochemical properties of the enzyme's
environment. Biochemical mechanisms underlying such regulation include the
potential for concentrating reactants together, tuning reaction rates, and
controlling competing metabolic pathways. LLPS is thus a powerful tool with
extensive utilities for the cell, e.g. for consolidating cell survival under
stress or rerouting metabolic pathways in response to the energy state of the
cell. By examining the evidence of how LLPS affects enzyme catalysis, we can
begin to understand emerging concepts and expand our understanding of enzyme
catalysis in living cells.


  Affinity Maturation (AM) is the process through which the immune system is
able to develop potent antibodies against new pathogens it encounters, and is
at the base of the efficacy of vaccines. At its core AM is analogous to a
Darwinian evolutionary process, where B-cells mutate and are selected on the
base of their affinity for an Antigen (Ag), and Ag availability tunes the
selective pressure. In cases when this selective pressure is high the number of
B-cells might quickly decrease and the population might risk extinction in what
is known as a population bottleneck. Here we study the probability for a B-cell
lineage to survive this bottleneck scenario as a function of the progenitor
affinity for the Ag. Using recursive relations and probability generating
functions we derive expressions for the average extinction time and progeny
size for lineages that go extinct. We then extend our results to the full
population, both in the absence and presence of competition for T-cell help,
and quantify the population survival probability as a function of Ag
concentration and initial population size. Our study suggests the population
bottleneck phenomenology might represent a limit case in the space of
biologically plausible maturation scenarios, whose characterization could help
guide the process of vaccine development.


  The extraction of k-mers from sequencing reads is an important task in many
bioinformatics applications, such as all DNA sequence analysis methods based on
de Bruijn graphs. These methods tend to be more accurate when the used k-mers
are unique in the analyzed DNA, and thus the use of longer k-mers is preferred.
When the read lengths of short read sequencing technologies increase, the error
rate will become the determining factor for the largest possible value of k.
Here we propose LoMeX which uses spaced seeds to extract long k-mers accurately
even in the presence of sequencing errors. Our experiments show that LoMeX can
extract long k-mers from current Illumina reads with a higher recall than a
standard k-mer counting tool. Furthermore, our experiments on simulated data
show that when the read length further increases, the performance of standard
k-mer counters declines, whereas LoMeX still extracts long k-mers successfully.


  Missing genotypes can affect the efficacy of machine learning approaches to
identify the risk genetic variants of common diseases and traits. The problem
occurs when genotypic data are collected from different experiments with
different DNA microarrays, each being characterised by its pattern of uncalled
(missing) genotypes. This can prevent the machine learning classifier from
assigning the classes correctly. To tackle this issue, we used well-developed
notions of object-attribute biclusters and formal concepts that correspond to
dense subrelations in the binary relation $\textit{patients} \times
\textit{SNPs}$. The paper contains experimental results on applying a
biclustering algorithm to a large real-world dataset collected for studying the
genetic bases of ischemic stroke. The algorithm could identify large dense
biclusters in the genotypic matrix for further processing, which in return
significantly improved the quality of machine learning classifiers. The
proposed algorithm was also able to generate biclusters for the whole dataset
without size constraints in comparison to the In-Close4 algorithm for
generation of formal concepts.


  Epidemiological forecasts are beset by uncertainties about the underlying
epidemiological processes, and the surveillance process through which data are
acquired. We present a Bayesian inference methodology that quantifies these
uncertainties, for epidemics that are modelled by (possibly) non-stationary,
continuous-time, Markov population processes. The efficiency of the method
derives from a functional central limit theorem approximation of the
likelihood, valid for large populations. We demonstrate the methodology by
analysing the early stages of the COVID-19 pandemic in the UK, based on
age-structured data for the number of deaths. This includes maximum a
posteriori estimates, MCMC sampling of the posterior, computation of the model
evidence, and the determination of parameter sensitivities via the Fisher
information matrix. Our methodology is implemented in PyRoss, an open-source
platform for analysis of epidemiological compartment models.


  The outbreak of COVID-19 in 2020 has led to a surge in the interest in the
mathematical modeling of infectious diseases. Disease transmission may be
modeled as compartmental models, in which the population under study is divided
into compartments and has assumptions about the nature and time rate of
transfer from one compartment to another. Usually, they are composed of a
system of ordinary differential equations (ODE's) in time. A class of such
models considers the Susceptible, Exposed, Infected, Recovered, and Deceased
populations, the SEIRD model. However, these models do not always account for
the movement of individuals from one region to another. In this work, we extend
the formulation of SEIRD compartmental models to diffusion-reaction systems of
partial differential equations to capture the continuous spatio-temporal
dynamics of COVID-19. Since the virus spread is not only through diffusion, we
introduce a source term to the equation system, representing exposed people who
return from travel. We also add the possibility of anisotropic non-homogeneous
diffusion. We implement the whole model in \texttt{libMesh}, an open finite
element library that provides a framework for multiphysics, considering
adaptive mesh refinement and coarsening. Therefore, the model can represent
several spatial scales, adapting the resolution to the disease dynamics. We
verify our model with standard SEIRD models and show several examples
highlighting the present model's new capabilities.


  Trusting in others and reciprocating that trust with trustworthy actions are
crucial to successful and prosperous societies. The Trust Game has been widely
used to quantitatively study trust and trustworthiness, involving a sequential
exchange between an investor and a trustee. The deterministic evolutionary game
theory predicts no trust and no trustworthiness whereas the behavioural
experiments with the one-shot anonymous Trust Game show that people
substantially trust and respond trustworthily. To explain these discrepancies,
previous works often turn to additional mechanisms, which are borrowed from
other games such as Prisoner's Dilemma. Although these mechanisms lead to the
evolution of trust and trustworthiness to an extent, the optimal or the most
common strategy often involves no trustworthiness. In this paper, we study the
impact of asymmetric demographic parameters (e.g. different population sizes)
on game dynamics of the Trust Game. We show that, in weak-mutation limit,
stochastic evolutionary dynamics with the asymmetric parameters can lead to the
evolution of high trust and high trustworthiness without any additional
mechanisms in well-mixed finite populations. Even full trust and near full
trustworthiness can be the most common strategy. These results are
qualitatively different from those of the previous works. Our results thereby
demonstrate rich evolutionary dynamics of the asymmetric Trust Game.


  Genomic Epidemiology (genEpi) is a branch of public health that uses many
different data types including tabular, network, genomic, and geographic, to
identify and contain outbreaks of deadly diseases. Due to the volume and
variety of data, it is challenging for genEpi domain experts to conduct data
reconnaissance; that is, have an overview of the data they have and make
assessments toward its quality, completeness, and suitability. We present an
algorithm for data reconnaissance through automatic visualization
recommendation, GEViTRec. Our approach handles a broad variety of dataset types
and automatically generates coordinated combinations of charts, in contrast to
existing systems that primarily focus on singleton visual encodings of tabular
datasets. We automatically detect linkages across multiple input datasets by
analyzing non-numeric attribute fields, creating an entity graph within which
we analyze and rank paths. For each high-ranking path, we specify chart
combinations with spatial and color alignments between shared fields, using a
gradual binding approach to transform initial partial specifications of
singleton charts to complete specifications that are aligned and oriented
consistently. A novel aspect of our approach is its combination of
domain-agnostic elements with domain-specific information that is captured
through a domain-specific visualization prevalence design space. Our
implementation is applied to both synthetic data and real data from an Ebola
outbreak. We compare GEViTRec's output to what previous visualization
recommendation systems would generate, and to manually crafted visualizations
used by practitioners. We conducted formative evaluations with ten genEpi
experts to assess the relevance and interpretability of our results.


  Effective control of biological systems can often be achieved through the
control of a surprisingly small number of distinct variables. We bring clarity
to such results using the formalism of Boolean dynamical networks, analyzing
the effectiveness of external control in selecting a desired final state when
that state is among the original attractors of the dynamics. Analyzing 49
existing biological network models, we find strong numerical evidence that the
average number of nodes that must be forced scales logarithmically with the
number of original attractors. This suggests that biological networks may be
typically easy to control even when the number of interacting components is
large. We provide a theoretical explanation of the scaling by separating
controlling nodes into three types: those that act as inputs, those that
distinguish among attractors, and any remaining nodes. We further identify
characteristics of dynamics that can invalidate this scaling, and speculate
about how this relates more broadly to non-biological systems.


  Trees corresponding to $\Lambda$- and $\Xi$-$n$-coalescents can be both quite
similar and fundamentally different compared to bifurcating tree models based
on Kingman's $n$-coalescent. This has consequences for inference of a
well-fitting gene genealogy as well as for assessing biological properties of
species having such sample genealogies. Here, mathematical properties
concerning clade sizes in the tree as well as changes of the tree when the
samples are enlargened are highlighted. To be used as realistic genealogy
models for real populations, an extension for changing population sizes is
discussed.


  Non-pharmaceutical measures such as social distancing, can play an important
role to control an epidemic in the absence of vaccinations. In this paper, we
study the impact of social distancing on epidemics for which it is executable.
We use a mathematical model combining human mobility and disease spreading. For
the mobility dynamics, we design an agent based model consisting of pedestrian
dynamics with a novel type of force to resemble social distancing in crowded
sites. For the spreading dynamics, we consider the compartmental SIE dynamics
plus an indirect transmission with the footprints of the infectious pedestrians
being the contagion factor. We show that the increase in the intensity of
social distancing has a significant effect on the exposure risk. By classifying
the population into social distancing abiders and non-abiders, we conclude that
the practice of social distancing, even by a minority of potentially infectious
agents, results in a drastic change on the population exposure risk, but
reduces the effectiveness of the protocols when practiced by the rest of the
population. Furthermore, we observe that for contagions which the indirect
transmission is more significant, the effectiveness of social distancing would
be reduced. This study can provide a quantitative guideline for policy-making
on exposure risk reduction.


  The exclusive presence of $\beta$-D-ribofuranose in nucleic acids is still a
conundrum in prebiotic chemistry, given that pyranose species are substantially
more stable at equilibrium. However, a precise characterisation of the relative
furanose/pyranose fraction at temperatures higher than about 50$^{\,\rm o}$C is
still lacking. Here, we employ a combination of NMR measurements and
statistical mechanics modelling to predict a population inversion between
furanose and pyranose at equilibrium at high temperatures. More importantly, we
show that a steady temperature gradient may steer an open isomerisation network
into a non-equilibrium steady state where furanose is boosted beyond the limits
set by equilibrium thermodynamics. Moreover, we demonstrate that nonequilibrium
selection of furanose is maximum at optimal dissipation, as gauged by the
temperature gradient and energy barriers for isomerisation. The predicted
optimum is compatible with temperature drops found in hydrothermal vents
associated with extremely fresh lava flows on the seafloor.


  This paper is concerned with the utilization of deterministically modeled
chemical reaction networks for the implementation of (feed-forward) neural
networks. We develop a general mathematical framework and prove that the
ordinary differential equations (ODEs) associated with certain reaction network
implementations of neural networks have desirable properties including (i)
existence of unique positive fixed points that are smooth in the parameters of
the model (necessary for gradient descent), and (ii) fast convergence to the
fixed point regardless of initial condition (necessary for efficient
implementation). We do so by first making a connection between neural networks
and fixed points for systems of ODEs, and then by constructing reaction
networks with the correct associated set of ODEs. We demonstrate the theory by
constructing a reaction network that implements a neural network with a
smoothed ReLU activation function, though we also demonstrate how to generalize
the construction to allow for other activation functions (each with the
desirable properties listed previously). As there are multiple types of
"networks" utilized in this paper, we also give a careful introduction to both
reaction networks and neural networks, in order to disambiguate the overlapping
vocabulary in the two settings and to clearly highlight the role of each
network's properties.


  Various methods have been developed to analyze the association between
organisms and their genomic sequences. Among them, sequence alignment is the
most frequently used for comparative analysis of biological genomes. However,
the traditional sequence alignment method is considerably complicated in
proportion to the sequences' length, and it is significantly challenging to
align long sequences such as a human genome. Currently, several multiple
sequence alignment algorithms are available that can reduce the complexity and
improve the alignment performance of various genomes. However, there have been
relatively fewer attempts to improve the alignment performance of the pairwise
alignment algorithm. After grasping these problems, we intend to propose a new
sequence alignment method using deep reinforcement learning. This research
shows the application method of the deep reinforcement learning to the sequence
alignment system and the way how the deep reinforcement learning can improve
the conventional sequence alignment method.


  The coil to globule transition of the polypeptide chain is the physical
phenomenon behind the folding of globular proteins. Globular proteins with a
single domain usually consist of about 30 to 100 amino acid residues, and this
finite size extends the transition interval of the coil-globule phase
transition. Based on the pedantic derivation of the two-state model, we
introduce the number of amino acid residues of a polypeptide chain as a
parameter in the expressions for two cooperativity measures and reveal their
physical significance. We conclude that the $k_2$ measure, related to the
degeneracy of the denatured state, describes the number of cooperative units
involved in the transition; additionally is found that the famous condition
$k_2=1$ is just the necessary condition to classify the protein as the
two-state folder. We also find that $\Omega_c$ is simply proportional to the
square of the transition interval. This fact allows us to perform the classical
size scaling analysis of the coil-globule phase transition. Moreover, these two
measures are shown to describe different characteristics of protein folding.


  In this paper, we obtain finite estimates and asymptotic formulas for the
expected number of real roots of two classes of random polynomials arising from
evolutionary game theory. As a consequence of our analysis, we achieve an
asymptotic formula for the expected number of internal equilibria in
multi-player two-strategy random evolutionary games. Our results contribute
both to evolutionary game theory and random polynomial theory.


  We investigate spreading and recovery of disease in a square lattice, and in
particular, emphasize the role of the initial distribution of infected patches
in the network, on the progression of an endemic and initiation of a recovery
process, if any, due to migration of both the susceptible and infected hosts.
The disease starts in the lattice with three possible initial distribution
patterns of infected and infection-free sites, infected core patches (ICP),
infected peripheral patches (IPP) and randomly distributed infected patches
(RDIP). Our results show that infection spreads monotonically in the lattice
with increasing migration without showing any sign of recovery in the ICP case.
In the IPP case, it follows a similar monotonic progression with increasing
migration, however, a self-organized healing process starts for higher
migration, leading the lattice to full recovery at a critical rate of
migration. Encouragingly, for the initial RDIP arrangement, chances of recovery
are much higher with a lower rate of critical migration. An eigenvalue based
semi-analytical study is made to determine the critical migration rate for
realizing a stable infection-free lattice. The initial fraction of infected
patches and the force of infection play significant roles in the self-organized
recovery. They follow an exponential law, for the RDIP case, that governs the
recovery process. For the frustrating case of ICP arrangement, we propose a
random rewiring of links in the lattice allowing long-distance migratory paths
that effectively initiate a recovery process. Global prevalence of infection
thereby declines and progressively improves with the rewiring probability that
follows a power law with the critical migration and leads to the birth of
emergent infection-free networks.


  In this work we propose a model that represents the relation between fish
ponds, the mosquito population and the transmission of malaria.
  It has been observed that in the Amazonic region of Acre, in the North of
Brazil, fish farming is correlated to the transmission of malaria when carried
out in artificial ponds that become breeding sites. Evidence has been found
indicating that cleaning the vegetation from the edges of the crop tanks helps
to control the size of the mosquito population.
  We use our model to determine the effective contribution of fish farming
practices on malaria transmission dynamics. The model consists of a nonlinear
system of ordinary differential equations with jumps at the cleaning time,
which act as impulsive controls. We study the asymptotic behaviour of the
system in function of the intensity and periodicity of the cleaning, and the
value of the parameters. In particular, we state sufficient conditions under
which the mosquito population is eliminated or persists, and under which the
malaria is eliminated or becomes endemic. We prove our conditions by applying
results for cooperative systems with concave nonlinearities.


  The antioxidant activity of baked foods is of utmost interest when
envisioning enhancing their health benefits. Incorporating functional
ingredients is challenging since their bioactivity naturally declines during
baking. In this study, 3D food printing and design of experiments are employed
to clarify how the antioxidant activity of cookies enriched with encapsulated
polyphenols can be maximized. A synergistic effect between encapsulation, time,
temperature, number of layers, and infill of the printed cookies was observed
on the moisture and antioxidant activity. Four-layer cookies with 30 % infill
provided the highest bioactivity and phenolic content if baked for 10 min and
at 180 {\deg}C. The bioacitivity and total phenolic content improved by 115 %
and 173 %, respectively, comparing to free extract cookies. Moreover, the
proper combination of the design and baking variables allowed to vary the
bioactivity of cooked cookies (moisture 3-5 %) between 300 to 700
{\mu}molTR/gdry. The additive manufacture of foods with interconnected pores
could accelerate baking and browning, or reduce thermal degradation. This
represents a potential approach to enhance the functional and healthy
properties of cookies or other thermal treated bioactive food products.


  Deep generative models have been applied with increasing success to the
generation of two dimensional molecules as SMILES strings and molecular graphs.
In this work we describe for the first time a deep generative model that can
generate 3D molecular structures conditioned on a three-dimensional (3D)
binding pocket. Using convolutional neural networks, we encode atomic density
grids into separate receptor and ligand latent spaces. The ligand latent space
is variational to support sampling of new molecules. A decoder network
generates atomic densities of novel ligands conditioned on the receptor.
Discrete atoms are then fit to these continuous densities to create molecular
structures. We show that valid and unique molecules can be readily sampled from
the variational latent space defined by a reference `seed' structure and
generated structures have reasonable interactions with the binding site. As
structures are sampled farther in latent space from the seed structure, the
novelty of the generated structures increases, but the predicted binding
affinity decreases. Overall, we demonstrate the feasibility of conditional 3D
molecular structure generation and provide a starting point for methods that
also explicitly optimize for desired molecular properties, such as high binding
affinity.


  Initially emerged in the Chinese city Wuhan and subsequently spread almost
worldwide causing a pandemic, the SARS-CoV-2 virus follows reasonably well the
SIR (Susceptible-Infectious-Recovered) epidemic model on contact networks in
the Chinese case. In this paper, we investigate the prediction accuracy of the
SIR model on networks also for Italy. Specifically, the Italian regions are a
metapopulation represented by network nodes and the network links are the
interactions between those regions. Then, we modify the network-based SIR model
in order to take into account the different lockdown measures adopted by the
Italian Government in the various phases of the spreading of the COVID-19. Our
results indicate that the network-based model better predicts the daily
cumulative infected individuals when time-varying lockdown protocols are
incorporated in the classical SIR model.


  We consider phylogeny estimation under a two-state model of sequence
evolution by site substitution on a tree. In the asymptotic regime where the
sequence lengths tend to infinity, we show that for any fixed $k$ no
statistically consistent phylogeny estimation is possible from $k$-mer counts
over the full leaf sequences alone. Formally, we establish that the joint
distribution of $k$-mer counts over the entire leaf sequences on two distinct
trees have total variation distance bounded away from $1$ as the sequence
length tends to infinity. Our impossibility result implies that statistical
consistency requires more sophisticated use of $k$-mer count information, such
as block techniques developed in previous theoretical work.


  A recent article in Nature Physics unified key results from thermodynamics,
statistics, and information theory. The unification arose from a general
equation for the rate of change in the information content of a system. The
general equation describes the change in the moments of an observable quantity
over a probability distribution. One term in the equation describes the change
in the probability distribution. The other term describes the change in the
observable values for a given state. We show the equivalence of this general
equation for moment dynamics with the widely known Price equation from
evolutionary theory, named after George Price. We introduce the Price equation
from its biological roots, review a mathematically abstract form of the
equation, and discuss the potential for this equation to unify diverse
mathematical theories from different disciplines. The new work in Nature
Physics and many applications in biology show that this equation also provides
the basis for deriving many novel theoretical results within each discipline.


  In order to accommodate the empirical fact that population structures are
rarely simple, modern studies of evolutionary dynamics allow for complicated
and highly-heterogeneous spatial structures. As a result, one of the most
difficult obstacles lies in making analytical deductions, either qualitative or
quantitative, about the long-term outcomes of evolution. The
"structure-coefficient" theorem is a well-known approach to this problem for
mutation-selection processes under weak selection, but a general method of
evaluating the terms it comprises is lacking. Here, we provide such a method
for populations of fixed (but arbitrary) size and structure, using easily
interpretable demographic measures. This method encompasses a large family of
evolutionary update mechanisms and extends the theorem to allow for asymmetric
contests to provide a better understanding of the mutation-selection balance
under more realistic circumstances. We apply the method to study social goods
produced and distributed among individuals in spatially-heterogeneous
populations, where asymmetric interactions emerge naturally and the outcome of
selection varies dramatically depending on the nature of the social good, the
spatial topology, and frequency with which mutations arise.


  Infectious diseases that incorporate pre-symptomatic transmission are
challenging to monitor, model, predict and contain. We address this scenario by
studying a variant of a stochastic susceptible-exposed-infected-recovered model
on arbitrary network instances using an analytical framework based on the
method of dynamic message-passing. This framework provides a good estimate of
the probabilistic evolution of the spread on both static and contact networks,
offering a significantly improved accuracy with respect to individual-based
mean-field approaches while requiring a much lower computational cost compared
to numerical simulations. It facilitates the derivation of epidemic thresholds,
which are phase boundaries separating parameter regimes where infections can be
effectively contained from those where they cannot. These have clear
implications on different containment strategies through topological (reducing
contacts) and infection parameter changes (e.g., social distancing and wearing
face masks), with relevance to the recent COVID-19 pandemic.


  We propose a simple model of spreading of some infection in an originally
healthy population which is different from other models existing in the
literature. In particular, we use an operator technique which allows us to
describe in a natural way the possible interactions between healthy and
un-healthy populations, and their transformation into recovered and to dead
people. After a rather general discussion, we apply our method to the analysis
of Chinese data for the SARS-2003 (Severe acute respiratory syndrome;
SARS-CoV-1) and the Coronavirus COVID-19 (Corona Virus Disease; SARS-CoV-2 )
and we show that the model works very well in reproducing the long-time
behaviour of the disease, and in particular in finding the number of affected
and dead people in the limit of large time. Moreover, we show how the model can
be easily modified to consider some lockdown measure, and we deduce that this
procedure drastically reduces the asymptotic value of infected individuals, as
expected, and observed in real life.


  Understanding how to limit biological invasion is critical, especially in the
context of accelerating anthropogenic ecological changes. Although biological
invasion success could be explained by the lack of natural enemies in new
regions, recent studies have revealed that resident herbivores often do have a
substantial effect on both native and invasive plants. Very few studies have
included consideration of native plant resistance while estimating methods of
controlling invasion; hence, it is unclear to what extent the interactive
effects of controlling approaches and native plants' resistance could slow down
or even inhibit biological invasion. We developed a spatial modeling framework,
using a paired logistic equation model, with considerations of the dispersal
processes, to capture the dynamics change of native and invasive plants under
various strategies of control. We found that when biocontrol agents could have
a strong effect on invasive plant, that could almost completely limit the
invasion, together with a high native plant resistance. However, a high
application frequency is needed make an efficient impact, whereas, a low
frequency treatment leads to nearly the same outcome as the no treatment case.
Lastly, we showed that evenly controlling a larger area with a weaker effect
still lead to a better outcome than focusing on small patches with a stronger
effect. Overall, this study has some management implications, such as how to
determine the optimal allocation strategy.


  The purpose of this work is to extend the Geant4-DNA Monte Carlo toolkit to
include electron interactions with the four DNA bases using a set of cross
sections recently implemented in Geant-DNA CPA100 models and available for
liquid water. Electron interaction cross sections for elastic scattering,
ionisation, and electronic excitation were calculated in the four DNA bases
adenine, thymine, guanine and cytosine. The electron energy range is extended
to include relativistic electrons. Elastic scattering cross sections were
calculated using the independent atom model with amplitude derived from ELSEPA
code. Relativistic Binary Encounter Bethe Vriens model was used to calculate
ionisation cross sections. The electronic excitation cross sections
calculations were based on the water cross sections following the same strategy
used in CPA100 code. These were implemented within the Geant4-DNA option6
physics constructor to extend its capability of tracking electrons in DNA
material in addition to liquid water. Since DNA nucleobases have different
molecular structure than water it is important to perform more accurate
simulations especially because DNA is considered the most radiosensitive
structure in cells. Differential and integrated cross sections calculations
were in good agreement with data from the literature for all DNA bases.
Stopping power, range and inelastic mean free path calculations in the four DNA
bases using this new extension of Geant4-DNA option6 are in good agreement with
calculations done by other studies, especially for high energy electrons. Some
deviations are shown at the low electron energy range, which could be
attributed to the different interaction models. Comparison with water
simulations shows obvious difference which emphasizes the need to include DNA
bases cross sections in track structure codes for better estimation of
radiation effects on biological material.


  The Basic Local Alignment Search Tool (BLAST) is currently the most popular
method for searching databases of biological sequences. BLAST compares
sequences via similarity defined by a weighted edit distance, which results in
it being computationally expensive. As opposed to working with edit distance, a
vector similarity approach can be accelerated substantially using modern
hardware or hashing techniques. Such an approach would require fixed-length
embeddings for biological sequences. There has been recent interest in learning
fixed-length protein embeddings using deep learning models under the hypothesis
that the hidden layers of supervised or semi-supervised models could produce
potentially useful vector embeddings. We consider transformer (BERT) protein
language models that are pretrained on the TrEMBL data set and learn
fixed-length embeddings on top of them with contextual lenses. The embeddings
are trained to predict the family a protein belongs to for sequences in the
Pfam database. We show that for nearest-neighbor family classification,
pretraining offers a noticeable boost in performance and that the corresponding
learned embeddings are competitive with BLAST. Furthermore, we show that the
raw transformer embeddings, obtained via static pooling, do not perform well on
nearest-neighbor family classification, which suggests that learning embeddings
in a supervised manner via contextual lenses may be a compute-efficient
alternative to fine-tuning.


  Many biomolecules have flexible structures, requiring distributional
estimates of their conformations. Experiments to acquire distributional data
typically measure pairs of labels separately, losing information on the joint
distribution. These data are assumed independent when estimating the
conformational ensemble. We developed a method to estimate the true joint
distribution from separately acquired measurements, testing it on two
biological systems. This method accurately reproduces the joint distribution
where known and generates testable predictions about complex conformational
ensembles.


  Ecological theory predicts that the soil seed bank stabilises the composition
of annual plant communities in the face of environmental variability. However,
long-term data on the community dynamics in the seed bank and the standing
vegetation are needed to test this prediction. We tested the hypothesis that
the composition of the seed bank undergoes lower temporal variability than the
standing vegetation in a nine-year study in Mediterranean, semi-arid, and arid
ecosystems. The composition of the seed bank was estimated by collecting soil
cores from the studied sites on an annual basis. Seedling emergence under
optimal watering conditions was measured in each soil core for three
consecutive years, to account for seed dormancy. In all sites, the composition
of the seed bank differed from the vegetation throughout the years.
Small-seeded and dormant-seeded species had a higher frequency in the seed bank
than in the standing vegetation. In contrast, functional group membership
(grasses vs. forbs) did not explain differences in species frequency between
the seed bank and the vegetation after controlling for differences between
grasses and forbs in seed mass and seed dormancy. Contrary to predictions, the
magnitude of year-to-year variability (the mean compositional dissimilarity
between consecutive years) was not lower in the seed bank than in the
vegetation in all sites. However, long-term compositional trends in the seed
bank were weaker than in the vegetation in the Mediterranean and semi-arid
sites. In the arid site where year-to-year variability was highest, no
long-term trends were observed. Overall, the effect of the seed bank on the
temporal variability of the vegetation in annual communities depends on site
conditions and time scale. While the year-to-year variability of the seed bank
is similar to the vegetation, the soil seed bank can buffer long-term trends.


  De novo design of molecules has recently enjoyed the power of generative deep
neural networks. Current approaches aim to generate molecules either resembling
the properties of the molecules of the training set or molecules that are
optimized with respect to specific physicochemical properties. None of the
methods generates molecules specific to a target protein. In the approach
presented here, we introduce a method which is conditioned on the protein
target sequence to generate de novo molecules that are relevant to the target.
We use an implementation adapted from Google's "Show and Tell" image caption
generation method, to generate SMILES strings of molecules from protein
sequence embeddings generated by a deep bi-directional language model ELMo.
ELMo is used to generate contextualized embedding vectors of the protein
sequence. Using reinforcement learning, the trained model is further optimized
through augmented episodic likelihood to increase the diversity of the
generated compounds compared to the training set. We used the model to generate
compounds for two major drug target families, i.e. for GPCRs and Tyrosine
Kinase targets. The model generated compounds which are structurally different
form the training set, while also being more similar to compounds known to bind
to the two families of drug targets compared to a random set of molecules. The
compounds further display reasonable synthesizability and drug-likeness scores.


  The dynamics of the structured water molecules in the hydration shell of the
DNA double helix is of paramount importance for the understanding of many
biological mechanisms. In particular, the vibrational dynamics of a water spine
that is formed in the DNA minor groove is the aim of the present study. Within
the framework of the developed phenomenological model, based on the approach of
DNA conformational vibrations, the modes of H-bond stretching, backbone
vibrations, and water translational vibrations have been established. The
calculated frequencies of translation vibrations of water molecules vary from
167 to 205 cm$^{-1}$ depending on the nucleotide sequence. The mode of water
vibrations higher than the modes of internal conformational vibrations of DNA.
The calculated frequencies of water vibrations have shown a sufficient
agreement with the experimental low-frequency vibrational spectra of DNA. The
obtained modes of water vibrations are observed in the same region of the
vibrational spectra of DNA as translation vibrations of water molecules in the
bulk phase. To distinguish the vibrations of water molecules in the DNA minor
groove from those in the bulk water, the dynamics of DNA with heavy water was
also considered. The results have shown that in the case of heavy water the
frequencies of vibrations decrease for about 10 cm$^{-1}$ that may be used in
the experiment to identify the mode of water vibrations in the spine of
hydration in DNA minor groove.


  A growing body of evidence suggests that the structure of mitochondrial
networks is poised near criticality, an intermediate regime lying in between
order and disorder. Such description fits well with the idea that biological
systems, in general, may benefit from the long-range correlations and large
flexibility conferred by a critical regime. Despite the attractiveness of this
proposal, a clear understanding of the possible scenarios leading these
networks to criticality is still lacking. In this work, we compared the
behavior of mitochondrial networks emerging from a dimensionless agent-based
(AB) model and a spatially explicit (SE) model, in which nodes are embedded on
a 2D lattice. In both scenarios, we described the position of the control
parameter at which mitochondrial networks exhibit a dynamical phase transition
as well as the size-dependency of several network features. Furthermore, we
showed that the mitochondrial networks from mouse embryonic fibroblasts
presented similar topologies to the ones generated using the AB model, while
their universal behavior is better described by a SE model. Using finite-size
scaling analysis conducted on models and empirical data we defined the
universality classes they belong and provided the theoretical boundaries for
the mechanisms governing mitochondrial network formation. Our findings predict
the full repertoire of dynamical behavior expected for real mitochondrial
networks under physiological and pathological conditions.


  Understanding protein structure-function relationships is a key challenge in
computational biology, with applications across the biotechnology and
pharmaceutical industries. While it is known that protein structure directly
impacts protein function, many functional prediction tasks use only protein
sequence. In this work, we isolate protein structure to make functional
annotations for proteins in the Protein Data Bank in order to study the
expressiveness of different structure-based prediction schemes. We present
PersGNN - an end-to-end trainable deep learning model that combines graph
representation learning with topological data analysis to capture a complex set
of both local and global structural features. While variations of these
techniques have been successfully applied to proteins before, we demonstrate
that our hybridized approach, PersGNN, outperforms either method on its own as
well as a baseline neural network that learns from the same information.
PersGNN achieves a 9.3% boost in area under the precision recall curve (AUPR)
compared to the best individual model, as well as high F1 scores across
different gene ontology categories, indicating the transferability of this
approach.


  A biologically motivated individual-based framework for evolution in
network-structured populations is developed that can accommodate
eco-evolutionary dynamics. This framework is used to construct a network birth
and death model. The evolutionary graph theory model, which considers
evolutionary dynamics only, is derived as a special case, highlighting
additional assumptions that diverge from real biological processes. This is
achieved by introducing a negative ecological feedback loop that suppresses
ecological dynamics by forcing births and deaths to be coupled. We also
investigate how fitness, a measure of reproductive success used in evolutionary
graph theory, is related to the life-history of individuals in terms of their
birth and death rates. In simple networks, these ecologically motivated
dynamics are used to provide new insight into the spread of adaptive mutations,
both with and without clonal interference. For example, the star network, which
is known to be an amplifier of selection in evolutionary graph theory, can
inhibit the spread of adaptive mutations when individuals can die naturally.


  Understanding the behaviour of hosts of SARS-CoV-2 is crucial to our
understanding of the virus. A comparison of environmental features related to
the incidence of SARS-CoV-2 with those of its potential hosts is critical. We
examine the distribution of coronaviruses among bats. We analyse the
distribution of SARS-CoV-2 in a nine-week period following lockdown in Italy,
Spain, and Australia. We correlate its incidence with environmental variables
particularly ultraviolet radiation, temperature, and humidity. We establish a
clear negative relationship between COVID-19 and ultraviolet radiation,
modulated by temperature and humidity. We relate our results with data showing
that the bat species most vulnerable to coronavirus infection are those which
live in environmental conditions that are similar to those that appear to be
most favourable to the spread of COVID-19. The SARS-CoV-2 ecological niche has
been the product of long-term coevolution of coronaviruses with their host
species. Understanding the key parameters of that niche in host species allows
us to predict circumstances where its spread will be most favourable. Such
conditions can be summarised under the headings of nocturnality and
seasonality. High ultraviolet radiation, in particular, is proposed as a key
limiting variable. We therefore expect the risk of spread of COVID-19 to be
highest in winter conditions, and in low light environments. Human activities
resembling those of highly social cave-dwelling bats (e.g. large nocturnal
gatherings or high density indoor activities) will only serve to compound the
problem of COVID-19.


  In computer science, we can theoretically neatly separate transmission and
processing of information, hardware and software, and programs and their
inputs. This is much more intricate in biology, Nevertheless, I argue that
Shannon's concept of information is useful in biology, although its application
is not as straightforward as many people think. In fact, the recently developed
theory of information decomposition can shed much light on the complementarity
between coding and regulatory, or internal and environmental information. The
key challenge that we formulate in this contribution is to understand how
genetic information and external factors combine to create an organism, and
conversely, how the genome has learned in the course of evolution how to
harness the environment, and analogously, how coding, regulation and spatial
organization interact in cellular processes.


  Coronavirus disease 2019 (COVID-19) poses a great threat to public health and
the economy worldwide. Currently, COVID-19 evolves in many countries to a
second stage, characterized by the need for the liberation of the economy and
relaxation of the human psychological effects. To this end, numerous countries
decided to implement adequate deconfinement strategies. After the first
prolongation of the established confinement, Morocco moves to the deconfinement
stage on May 20, 2020. The relevant question concerns the impact on the
COVID-19 propagation by considering an additional degree of realism related to
stochastic noises due to the effectiveness level of the adapted measures. In
this paper, we propose a delayed stochastic mathematical model to predict the
epidemiological trend of COVID-19 in Morocco after the deconfinement. To ensure
the well-posedness of the model, we prove the existence and uniqueness of a
positive solution. Based on the large number theorem for martingales, we
discuss the extinction of the disease under an appropriate threshold parameter.
Moreover, numerical simulations are performed in order to test the efficiency
of the deconfinement strategies chosen by the Moroccan authorities to help the
policy makers and public health administration to make suitable decisions in
the near future.


  We examine a modification of the Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP)
process in which the diffusing substance requires a parent density field for
reproduction. A biological example would be the density of diffusing spores
(propagules) and the density of a stationary fungus (parent). The parent
produces propagules at a certain rate, and the propagules turn into the parent
substance at another rate. We model this evolution by the FKPP process with
delay, which reflects a finite time typically required for a new parent to
mature before it begins to produce propagules. While the FKPP process with
other types of delays have been considered in the past as a pure mathematical
construct, in our work a delay in the FKPP model arises in a natural science
setting. The speed of the resulting density fronts is shown to decrease with
increasing delay time, and has a non-trivial dependence on the rate of
conversion of propagules into the parent substance. Remarkably, the fronts in
this model are always slower than Fisher waves of the classical FKPP model. The
largest speed is half of the classical value, and it is achieved at zero delay
and when the two rates are matched.


  The relations, rather than the elements, constitute the structure of
networks. We therefore develop a systematic approach to the analysis of
networks, modelled as graphs or hypergraphs, that is based on structural
properties of (hyper)edges, instead of vertices. For that purpose, we utilize
so-called network curvatures. These curvatures quantify the local structural
properties of (hyper)edges, that is, how, and how well, they are connected to
others. In the case of directed networks, they assess the input they receive
and the output they produce, and relations between them. With those tools, we
can investigate biological networks. As examples, we apply our methods here to
protein-protein interaction, transcriptional regulatory and metabolic networks.


  Aggregating transcriptomics data across hospitals can increase sensitivity
and robustness of differential expression analyses, yielding deeper clinical
insights. As data exchange is often restricted by privacy legislation,
meta-analyses are frequently employed to pool local results. However, if class
labels are inhomogeneously distributed between cohorts, their accuracy may
drop. Flimma (https://exbio.wzw.tum.de/flimma/) addresses this issue by
implementing the state-of-the-art workflow limma voom in a privacy-preserving
manner, i.e. patient data never leaves its source site. Flimma results are
identical to those generated by limma voom on combined datasets even in
imbalanced scenarios where meta-analysis approaches fail.


  The assessment of microbial species biodiversity is essential in ecology and
evolutionary biology (Reaka-Kudla et al. 1996), but especially challenging for
communities of microorganisms found in the environment (Das et al. 2006,
Hillebrand et al. 2018). Beyond providing a census of organisms in the ocean,
assessing marine microbial biodiversity can reveal how microbes respond to
environmental change (Salazar et al. 2017), clarify the ecological roles of
community members (Hehemann et al. 2016), and lead to biotechnology discoveries
(Das et al. 2006). Computational approaches to characterize taxonomic diversity
and phylogeny based on the quality of available data for environmental sequence
datasets is fundamental for advancing our understanding of the role of these
organisms in the environment. Even more pressing is the need for comprehensive
and consistent methods to assign taxonomy to environmentally-relevant microbial
eukaryotes. Here, we present EUKulele, an open-source software tool designed to
assign taxonomy to microeukaryotes detected in meta-omic samples, and
complement analysis approaches in other domains by accommodating assembly
output and providing concrete metrics reporting the taxonomic completeness of
each sample.


  In this paper, we study an optimal control problem of a communicable disease
in a prison population. In order to control the spread of the disease inside a
prison, we consider an active case-finding strategy, consisting on screening a
proportion of new inmates at the entry point, followed by a treatment depending
on the results of this procedure. The control variable consists then in the
coverage of the screening applied to new inmates. The disease dynamics is
modeled by a SIS (susceptible-infected-susceptible) model, typically used to
represent diseases that do not confer immunity after infection. We determine
the optimal strategy that minimizes a combination between the cost of the
screening/treatment at the entrance and the cost of maintaining infected
individuals inside the prison, in a given time horizon. Using the Pontryagin
Maximum Principle and Hamilton-Jacobi-Bellman equation, among other tools, we
provide the complete synthesis of an optimal feedback control, consisting in a
bang-bang strategy with at most two switching times and no singular arc
trajectory, characterizing different profiles depending on model parameters.


  Compartmental transmission models have become an invaluable tool to study the
dynamics of infectious diseases. The Susceptible-Infectious-Recovered (SIR)
model is known to have an exact semi-analytical solution. In the current study,
the approach of Harko et al. (2014) is generalised to obtain an approximate
semi-analytical solution of the Susceptible-Exposed-Infectious-Recovered (SEIR)
model. The SEIR model curves have nearly the same shapes as the SIR ones, but
with a stretch factor applied to them across time that is related to the ratio
of the incubation to infectious periods. This finding implies an approximate
characteristic timescale, scaled by this stretch factor, that is universal to
all SEIR models, which only depends on the basic reproduction number and
initial fraction of the population that is infectious.


  Best match graphs (BMG) are a key intermediate in graph-based orthology
detection and contain a large amount of information on the gene tree. We
provide a near-cubic algorithm to determine whether a BMG is
binary-explainable, i.e., whether it can be explained by a fully resolved gene
tree and, if so, to construct such a tree. Moreover, we show that all such
binary trees are refinements of the unique binary-resolvable tree (BRT), which
in general is a substantial refinement of the also unique least resolved tree
of a BMG. Finally, we show that the problem of editing an arbitrary
vertex-colored graph to a binary-explainable BMG is NP-complete and provide an
integer linear program formulation for this task.


  Gliomas are lethal type of central nervous system tumors with a poor
prognosis. Recently, with the advancements in the micro-array technologies
thousands of gene expression related data of glioma patients are acquired,
leading for salient analysis in many aspects. Thus, genomics are been emerged
into the field of prognosis analysis. In this work, we identify survival
related 7 gene signature and explore two approaches for survival prediction and
risk estimation. For survival prediction, we propose a novel probabilistic
programming based approach, which outperforms the existing traditional machine
learning algorithms. An average 4 fold accuracy of 74% is obtained with the
proposed algorithm. Further, we construct a prognostic risk model for risk
estimation of glioma patients. This model reflects the survival of glioma
patients, with high risk for low survival patients.


  At the time of writing, the ongoing COVID-19 pandemic, caused by severe acute
respiratory syndrome coronavirus 2 (SARS-CoV-2), had already resulted in more
than thirty-two million cases infected and more than one million deaths
worldwide.
  Given the fact that the pandemic is still threatening health and safety, it
is in the urgency to understand the COVID-19 contagion process and know how it
might be controlled. With this motivation in mind, in this paper, we consider a
version of a stochastic discrete-time
Susceptible-Infected-Recovered-Death~(SIRD)-based epidemiological model with
two uncertainties: The uncertain rate of infected cases which are undetected or
asymptomatic, and the uncertain effectiveness rate of control. Our aim is to
study the effect of an epidemic control policy on the uncertain model in a
control-theoretic framework. We begin by providing the closed-form solutions of
states in the modified SIRD-based model such as infected cases, susceptible
cases, recovered cases, and deceased cases. Then, the corresponding expected
states and the technical lower and upper bounds for those states are provided
as well. Subsequently, we consider two epidemic control problems to be
addressed: One is almost sure epidemic control problem and the other average
epidemic control problem. Having defined the two problems, our main results are
a set of sufficient conditions on a class of linear control policy which
assures that the epidemic is "well-controlled"; i.e., both of the infected
cases and deceased cases are upper bounded uniformly and the number of infected
cases converges to zero asymptotically. Our numerical studies, using the
historical COVID-19 contagion data in the United States, suggest that our
appealingly simple model and control framework can provide a reasonable
epidemic control performance compared to the ongoing pandemic situation.


  Antibiotic resistance is a worldwide challenge. A potential approach to block
resistance is to simultaneously inhibit WT and known escape variants of the
target bacterial protein. Here we applied an integrated computational and
experimental approach to discover compounds that inhibit both WT and
trimethoprim (TMP) resistant mutants of E. coli dihydrofolate reductase (DHFR).
We identified a novel compound (CD15-3) that inhibits WT DHFR and its TMP
resistant variants L28R, P21L and A26T with IC50 50-75 micromoles against WT
and TMP-resistant strains. Resistance to CD15-3 was dramatically delayed
compared to TMP in in vitro evolution. Whole genome sequencing of CD15-3
resistant strains showed no mutations in the target folA locus. Rather, gene
duplication of several efflux pumps gave rise to weak (about twofold increase
in IC50) resistance against CD15-3. Altogether, our results demonstrate the
promise of strategy to develop evolution drugs - compounds which block
evolutionary escape routes in pathogens.


  The metazoan body plan is established during early embryogenesis via
collective cell rearrangements and evolutionarily conserved gene networks, as
part of a process commonly referred to as gastrulation. While substantial
progress has been achieved in terms of characterizing the embryonic development
of several model organisms, underlying principles of many early patterning
processes nevertheless remain enigmatic. Despite the diversity of
(pre-)gastrulating embryo and adult body shapes across the animal kingdom, the
body axes, which are arguably the most fundamental features, generally remain
identical between phyla. Recently there has been a renewed appreciation of ex
vivo and in vitro embryo-like systems to model early embryonic patterning
events. Here, we briefly review key examples and propose that similarities in
morphogenesis as well as associated gene expression dynamics may reveal an
evolutionarily conserved developmental mode as well as provide further insights
into the role of external or extraembryonic cues in shaping the early embryo.
In summary, we argue that embryo-like systems can be employed to inform
previously uncharted aspects of animal body plan evolution as well as
associated patterning rules.


  In this work of speculative science, scientists from a distant star system
explain the emergence and consequences of triparentalism, when three
individuals are required for sexual reproduction, which is the standard form of
mating on their home world. The report details the evolution of their
reproductive system--that is, the conditions under which triparentalism and
three self-avoiding mating types emerged as advantageous strategies for sexual
reproduction. It also provides an overview of the biological consequences of
triparental reproduction with three mating types, including the genetic
mechanisms of triparental reproduction, asymmetries between the three mating
types, and infection dynamics arising from their different mode of sexual
reproduction. The report finishes by discussing how central aspects of their
society, such as short-lasting unions among individuals and the rise of a
monoculture, might have arisen as a result of their triparental system.


  Humans, domestic animals, orchard crops, and ornamental plants are commonly
treated with antibiotics in response to bacterial infection. By curing
infectious individuals, antibiotic therapy might limit the spread of contagious
disease among hosts. But an antibiotic`s suppression of within-host pathogen
density might also reduce the probability that the host is otherwise removed
from infectious status prior to therapeutic recovery. When rates of both
recovery via treatment and other removal events (e.g., isolation or mortality)
depend directly on within-host pathogen density, antibiotic treatment can relax
the overall removal rate sufficiently to increase between-host disease
transmission. To explore this dependence, a deterministic within-host dynamics
drives the infectious host's time-dependent probability of disease
transmission, as well as the probabilistic duration of the infectious period.
At the within-host scale, the model varies (1) inoculum size, (2) bacterial
self-regulation, (3) the time between infection and initiation of therapy, and
(4) antibiotic efficacy. At the between-host scale the model varies (5) the
size/susceptibility of groups randomly encountered by an infectious host.
Results identify conditions where antibiotic treatment can increase duration of
a host`s infectiousness, and consequently increase the expected number of new
infections. At lower antibiotic efficacy, treatment might convert a rare,
serious bacterial disease into a common, but treatable infection.


  In this paper we study the spatial spread of the COVID-19 infection in
Lebanon. We inspect the spreading of the daily new infections across the 26
administrative districts of the country, and implement Moran's $I$ statistics
in order to analyze the tempo-spatial clustering of the infection in relation
to various variables parameterized by adjacency, proximity, population,
population density, poverty rate and poverty density, and we find out that
except for the poverty rate, the spread of the infection is clustered and
associated to those parameters with varying magnitude for the time span between
July (geographic adjacency and proximity) or August (population, population
density and poverty density) through October. We also determine the temporal
dynamics of geographic location of the mean center of new and cumulative
infections since late March. The results obtained allow for regionally and
locally adjusted health policies and measures that would provide higher levels
of public health safety in the country.


  A large body of empirical evidence suggests that humans are willing to engage
in costly punishment of defectors in public goods games. Based on such pieces
of evidence, it is suggested that punishment serves an important role in
promoting cooperation in humans, and possibly other species. Nevertheless,
theoretical work has been unable to show how this is possible. The problem
originates from the fact that punishment, being costly, is an altruistic act
and its evolution is subject to the same problem that it tries to address. To
suppress this so-called second-order free-rider problem, known theoretical
models on the evolution of punishment resort to one of the few established
mechanisms for the evolution of cooperation. This leaves the question that
whether altruistic punishment can evolve and give rise to the evolution of
cooperation, unaddressed. Here, by considering a population of individuals who
play a public goods game, followed by a public punishing game, introduced here,
we show that altruistic punishment indeed evolves and promotes cooperation, in
a general environment and in the absence of a cooperation favoring mechanism.
Besides, our analysis shows, being close to a physical phase transition
facilitates the evolution of altruistic punishment.


  Machine learning based methods have shown potential for optimizing existing
molecules with more desirable properties, a critical step towards accelerating
new chemical discovery. Here we propose QMO, a generic query-based molecule
optimization framework that exploits latent embeddings from a molecule
autoencoder. QMO improves the desired properties of an input molecule based on
efficient queries, guided by a set of molecular property predictions and
evaluation metrics. We show that QMO outperforms existing methods in the
benchmark tasks of optimizing small organic molecules for drug-likeness and
solubility under similarity constraints. We also demonstrate significant
property improvement using QMO on two new and challenging tasks that are also
important in real-world discovery problems: (i) optimizing existing potential
SARS-CoV-2 Main Protease inhibitors toward higher binding affinity; and (ii)
improving known antimicrobial peptides towards lower toxicity. Results from QMO
show high consistency with external validations, suggesting effective means to
facilitate material optimization problems with design constraints.