  Enhanced stabilisation of protein structures via the presence of inert
excipients is a key mechanism adopted both by physiological systems and in
biotechnological applications. While the intrinsic stability of proteins is
ultimately fixed by their amino acid composition and organisation, the
interactions between excipients and proteins together with their concentrations
introduce an additional layer of complexity and in turn, method of modulating
protein stability. Here, we combined experimental measurements with molecular
dynamics simulations and graph-theory based analyses to assess the
stabilising/destabilising effects of different kinds of osmolytes on proteins
during heat-mediated denaturation. We found that (i) proteins in solution with
stability-enhancing osmolytes tend to have more compact interaction networks
than those assumed in presence of destabilising excipients; (ii) a strong
negative correlation (R = -0.85) characterises the relationship between the
melting temperature Tm and the preferential interaction coefficient defined by
the radial distribution functions of osmolytes and water around the protein and
(iii) a positive correlation exists between osmolyte-osmolyte clustering and
the extent of preferential exclusion from the local domain of the protein,
suggesting that exclusion may be driven by enhanced steric hindrance of
aggregated osmolytes.


  Molecular carriers represent an increasingly common strategy in the field of
nanopore sensing to use secondary molecules to selectively report on the
presence of target analytes in solution, allowing for sensitive assays of
otherwise hard-to-detect molecules such as small, weakly-charged proteins.
However, existing carrier designs can often introduce drawbacks to nanopore
experiments including higher levels of cost/complexity and carrier-pore
interactions that lead to ambiguous signals and elevated clogging rates. In
this work, we present a simple method of carrier production based on
sticky-ended DNA molecules that emphasizes ease-of-synthesis and compatibility
with nanopore sensing and analysis. In particular, our method incorporates the
ability to flexibly control the length of the DNA carriers produced, enhancing
the multiplexing potential of this carrier system through the separable
nanopore signals they could generate for distinct targets. A proof-of-concept
nanopore experiment is also presented, involving carriers produced by our
method with multiple lengths and attached to DNA nanostructure targets, in
order to validate the capabilities of the system. As the breadth of
applications for nanopore sensors continues to expand, the availability of
tools such as those presented here to help translate the outcomes of these
applications into robust nanopore signals will be of major importance.


  We constructed a modular, biomimetic red panda paw with which to
experimentally investigate the evolutionary reason for the existence of the
false thumbs of red pandas. These thumbs were once believed to have shared a
common origin with the similar false thumbs of giant pandas; however, the
discovery of a carnivorous fossil ancestor of the red panda that had false
thumbs implies that the red panda did not evolve its thumbs to assist in eating
bamboo, as the giant panda did, but rather evolved its thumbs for some other
purpose. The leading proposal for this purpose is that the thumbs developed to
aid arboreal locomotion. To test this hypothesis, we conducted grasp tests on
rods 5-15 mm in diameter using a biomimetic paw with 0-16 mm interchangeable
thumb lengths. The results of these tests demonstrated an optimal thumb length
of 7 mm, which is just above that of the red panda's true thumb length of 5.5
mm. Given trends in the data that suggest that smaller thumbs are better suited
to grasping larger diameter rods, we conclude that the red panda's thumb being
sized below the optimum length suggests an adaptation toward grasping branches
as opposed to relatively thinner food items, supporting the new proposal that
the red panda's thumbs are an adaptation primary to climbing rather than food
manipulation.


  Type I CRISPR-Cas systems are the most common among six types of CRISPR-Cas
systems, however, non-self-targeting genome editing based on a single Cas3 of
type I CRISPR-Cas systems has not been reported. Here, we present the subtype
I-B-Svi CRISPR-Cas system (with three confirmed CRISPRs and a cas gene cluster)
and genome editing based on this system found in Streptomyces virginiae IBL14.
Importantly, like the animal-derived bacterial protein SpCas9 (1368
amino-acids), the single, compact, non-animal-derived bacterial protein SviCas3
(771 amino-acids) can also direct template-based microbial genome editing
through the target cell's own homology-directed repair system, which breaks the
view that the genome editing based on type I CRISPR-Cas systems requires a full
Cascade. Notably, no off-target changes or indel-formation were detected in the
analysis of potential off-target sites. This discovery broadens our
understanding of the diversity of type I CRISPR-Cas systems and will facilitate
new developments in genome editing tools.


  This study explores the vaccine prioritization strategy to reduce the overall
burden of the pandemic when the supply is limited. Existing methods conduct
macro-level or simplified micro-level vaccine distribution by assuming the
homogeneous behavior within subgroup populations and lacking mobility dynamics
integration. Directly applying these models for micro-level vaccine allocation
leads to sub-optimal solutions due to the lack of behavioral-related details.
To address the issue, we first incorporate the mobility heterogeneity in
disease dynamics modeling and mimic the disease evolution process using a
Trans-vaccine-SEIR model. Then we develop a novel deep reinforcement learning
to seek the optimal vaccine allocation strategy for the high-degree
spatial-temporal disease evolution system. The graph neural network is used to
effectively capture the structural properties of the mobility contact network
and extract the dynamic disease features. In our evaluation, the proposed
framework reduces 7% - 10% of infections and deaths than the baseline
strategies. Extensive evaluation shows that the proposed framework is robust to
seek the optimal vaccine allocation with diverse mobility patterns in the
micro-level disease evolution system. In particular, we find the optimal
vaccine allocation strategy in the transit usage restriction scenario is
significantly more effective than restricting cross-zone mobility for the top
10% age-based and income-based zones. These results provide valuable insights
for areas with limited vaccines and low logistic efficacy.


  Experimental work suggests that biological soil crusts, dominant primary
producers in drylands and tundra, are particularly vulnerable to disturbances
that cause reverse ecological succession. To model successional transitions in
biocrust communities, we propose a resource-firing game that captures
succession dynamics without specifying detailed function forms. The model is
evaluated in idealized terrestrial ecosystems, where disturbances are modeled
as a reduction in available resources that triggers inter-species competition.
The resource-firing game is executed on a finite graph with nodes representing
species in the community and a sink node that becomes active when every species
is depleted of resources. First, we discuss the theoretical basis of the
resource-firing game, evaluate it in the light of existing literature, and
consider the characteristics of a biocrust community that has evolved to
equilibrium. We then examine the dependence of resource-firing and game
stability on species richness, showing that high species richness increases the
probability of very short and long avalanches, but not those of intermediate
length. Indeed, this result suggests that the response of the community to
disturbance is both directional and episodic, proceeding towards reverse
succession in bursts of variable length. Finally, we incorporate the spatial
structure of the biocrust community into a Cayley Tree and derive a formula for
the probability that a disturbance, modeled as a random attack, initiates a
large species-death event.


  The probability that life spontaneously emerges in a suitable environment
(abiogenesis) is one of the major unknowns in astrobiology. Assessing its value
is impeded by the lack of an accepted theory for the origin of life, and is
further complicated by the existence of selection biases. Appealing
uncritically to some version of the ``Principle of Mediocrity'' -- namely, the
supposed typicality of what transpired on Earth -- is problematic on empirical
or logical grounds. In this paper, we adopt a Bayesian statistical approach to
put on rigorous footing the inference of lower bounds for the probability of
abiogenesis, based on current and future evidence. We demonstrate that the
single datum that life has appeared at least once on Earth merely sets weak
constraints on the minimal probability of abiogenesis. In fact, the {\it a
priori} probability assigned to this event (viz., optimistic, pessimistic or
agnostic prior) exerts the strongest influence on the final result. We also
show that the existence of a large number of habitable worlds does not
necessarily imply, by itself, a high probability that life should be common in
the universe. Instead, as delineated before, the choice of prior, which is
subject to uncertainty (i.e., admits multiple scenarios), strongly influences
the likelihood of life being common. If habitable worlds are uncommon, for an
agnostic prior, a deterministic scenario for the origin of life might be
favoured over one where abiogenesis is a fluke event.


  Cell lineage statistics is a powerful tool for inferring cellular parameters,
such as division rate, death rate or the population growth rate. Yet, in
practice such an analysis suffers from a basic problem: how should we treat
incomplete lineages that do not survive until the end of the experiment? Here,
we develop a model-independent theoretical framework to address this issue. We
show how to quantify fitness landscape, survivor bias and selection for
arbitrary cell traits from cell lineage statistics in the presence of death,
and we test this method using an experimental data set in which a cell
population is exposed to a drug that kills a large fraction of the population.
This analysis reveals that failing to properly account for dead lineages can
lead to misleading fitness estimations. For simple trait dynamics, we prove and
illustrate numerically that the fitness landscape and the survivor bias can in
addition be used for the non-parametric estimation of the division and death
rates, using only lineage histories. Our framework provides universal bounds on
the population growth rate, and a fluctuation-response relation which
quantifies the reduction of population growth rate due to the variability in
death rate. Further, in the context of cell size control, we obtain
generalizations of Powell's relation that link the distributions of generation
times with the population growth rate, and show that the survivor bias can
sometimes conceal the adder property, namely the constant increment of volume
between birth and division.


  The statistical field theory of information dynamics on complex networks
concerns the dynamical evolution of large classes of models of complex systems.
Previous work has focused on networks where nodes carry an information field,
which describes the internal state of the node, and its dynamical evolution. In
this work, we propose a more general mathematical framework to model
information dynamics on complex networks, where the internal node states are
vector valued, thus allowing each node to carry multiple types of information.
This setup is relevant for many biophysical and socio-technological models of
complex systems, ranging from viral dynamics on networks to models of opinion
dynamics and social contagion. The full dynamics of these systems is described
in the space of all possible network configurations, as opposed to a node-based
perspective. Here, we illustrate the mathematical framework presented in an
accompanying letter, while focusing on an exemplary application of epidemic
spreading on a low-dimensional network.


  Spatiotemporal coordination of chromatin and subnuclear compartments is
crucial for cells. Numerous enzymes act inside nucleus\textemdash some of those
transiently link and pass two chromatin segments. Here we study how such an
active perturbation affects fluctuating dynamics of an inclusion in the
chromatic medium. Using numerical simulations and a versatile effective model,
we categorize inclusion dynamics into three distinct modes. The
transient-link-and-pass activity speeds up inclusion dynamics by affecting a
slow mode related to chromatin remodeling, viz., size and shape of the
chromatin meshes.


  Mitochondria are critical organelles in eukaryotes that produce the energy
currency ATP. In nerve axons, mitochondria are known to align at almost regular
intervals to maintain a constant ATP concentration, but little is known about
the mechanism. In this letter, we show theoretically that ATP production and
ATP-dependent non-directional movement of mitochondria are sufficient for
alignment, even in the absence of an explicit repulsive force between them.
This is similar to thermodynamic forces driven by thermal fluctuations, even
generated by non-equilibrium processes, and demonstrates the diversity of
mechanisms governing the motion of biological matter.


  Protons in the gap between base pairs of the double helix store the code of
life by breaking the chiral symmetry that swaps the sense strand with its
complementary partner. When these hydrogen bonds break during replication and
transcription, pairs of protons switch sides restoring chiral symmetry and
destroying genetic information. Using time-independent second-order
perturbation theory, we show that the observed rate of such spontaneous
mutations follows in the sudden approximation for bond breaking provided
protons in bonds between bases tunnel across the gap with terahertz
frequencies.


  We investigate the long-time dynamics of a SIR epidemic model with infinitely
many pathogen variants infecting a homogeneous host population. We show that
the basic reproduction number $\mathcal{R}_0$ of the pathogen can be defined in
that case and corresponds to a threshold between the persistence
($\mathcal{R}_0>1$) and the extinction ($\mathcal{R}_0\leq 1$) of the pathogen.
When $\mathcal{R}_0>1$ and the maximal fitness is attained by at least one
variant, we show that the systems reaches an equilibrium state that can be
explicitly determined from the initial data. When $\mathcal{R}_0>1$ but none of
the variants attain the maximal fitness, the situation is more intricate.
  We show that, in general, the pathogen is uniformly persistent and any family
of variants that have a fitness which is uniformly lower than the optimal
fitness, eventually gets extinct. We derive a condition under which the total
pathogen population converges to a limit which can be computed explicitly.
  We also find counterexamples that show that, when our condition is not met,
the total pathogen population may converge to an unexpected value, or the
system can even reach an eternally transient behavior where the total pathogen
population between several values. We illustrate our results with numerical
simulations that emphasize the wide variety of possible dynamics.


  Model selection is critical in the modern statistics and machine learning
community. However, most existing works do not apply to heavy-tailed data,
which are commonly encountered in real applications, such as the single-cell
multiomics data. In this paper, we propose a rank-sum based approach that
outputs a confidence set containing the optimal model with guaranteed
probability. Motivated by conformal inference, we developed a general method
that is applicable without moment or tail assumptions on the data. We
demonstrate the advantage of the proposed method through extensive simulation
and a real application on the COVID-19 genomics dataset (Stephenson et al.,
2021). To perform the inference on rank-sum statistics, we derive a general
Gaussian approximation theory for high dimensional two-sample U-statistics,
which may be of independent interest to the statistics and machine learning
community.


  RNA-guided gene editing based on the CRISPR-Cas system is currently the most
effective genome editing technique. Here, we report that the SviCas3 from the
subtype I-B-Svi Cas system in Streptomyces virginiae IBL14 is an RNA-guided and
DNA-guided DNA endonuclease suitable for the HDR-directed gene and/or base
editing of eukaryotic cell genomes. The genome editing efficiency of SviCas3
guided by DNA is no less than that of SviCas3 guided by RNA. In particular,
t-DNA, as a template and a guide, does not require a proto-spacer-adjacent
motif, demonstrating that CRISPR, as the basis for crRNA design, is not
required for the SviCas3-mediated gene and base editing. This discovery will
broaden our understanding of enzyme diversity in CRISPR-Cas systems, will
provide important tools for the creation and modification of living things and
the treatment of human genetic diseases, and will usher in a new era of
DNA-guided gene editing and base editing.


  Traditional drug discovery programs are being transformed by the advent of
machine learning methods. Among these, Generative AI methods (GM) have gained
attention due to their ability to design new molecules and enhance specific
properties of existing ones. However, current GM methods have limitations, such
as low affinity towards the target, unknown ADME/PK properties, or the lack of
synthetic tractability. To improve the applicability domain of GM methods, we
have developed a workflow based on a variational autoencoder coupled with
active learning steps. The designed GM workflow iteratively learns from
molecular metrics, including drug likeliness, synthesizability, similarity, and
docking scores. In addition, we also included a hierarchical set of criteria
based on advanced molecular modeling simulations during a final selection step.
We tested our GM workflow on two model systems, CDK2 and KRAS. In both cases,
our model generated chemically viable molecules with a high predicted affinity
toward the targets. Particularly, the proportion of high-affinity molecules
inferred by our GM workflow was significantly greater than that in the training
data. Notably, we also uncovered novel scaffolds significantly dissimilar to
those known for each target. These results highlight the potential of our GM
workflow to explore novel chemical space for specific targets, thereby opening
up new possibilities for drug discovery endeavors.


  How predictable is the next move of an animal? Specifically, which factors
govern the short- and long-term motion patterns and the overall dynamics of
landbound, plant-eating animals and ruminants in particular? To answer this
question, we here study the movement dynamics of springbok antelopes Antidorcas
marsupialis. We propose complementary statistical analysis techniques combined
with machine learning approaches to analyze, across multiple time scales, the
springbok motion recorded in long-term GPS-tracking of collared springboks at a
private wildlife reserve in Namibia. As a new result, we are able to predict
the springbok movement within the next hour with a certainty of about 20\%. The
remaining 80\% are stochastic in nature and are induced by unaccounted factors
in the modeling algorithm and by individual behavioral features of springboks.
We find that directedness of motion contributes approximately 17\% to this
predicted fraction. We find that the measure for directedeness is strongly
dependent on the daily cycle. The previously known daily affinity of springboks
to their water points, as predicted from our machine learning algorithm,
overall accounts for only 3\% of this predicted deterministic component of
springbok motion. Moreover, the resting points are found to affect the motion
of springboks at least as much as the formally studied effects of water points.
The generality of these statements for the motion patterns and their underlying
behavioral reasons for other ruminants can be examined on the basis of our
statistical analysis tools.


  We model the effects of disease and other exogenous damage during human
aging. Even when the exogenous damage is repaired at the end of acute disease,
propagated secondary damage remains. We consider both short-term mortality
effects due to (acute) exogenous damage and long-term mortality effects due to
propagated damage within the context of a generic network model (GNM) of
individual aging that simulates a U.S. population. Across a wide range of
disease durations and severities we find that while excess short-term mortality
is highest for the oldest individuals, the long-term years of life lost are
highest for the youngest individuals. These appear to be universal effects of
human disease. We support this conclusion with a phenomenological model
coupling damage and mortality. Our results are consistent with previous
lifetime mortality studies of atom bomb survivors and post-recovery health
studies of COVID-19. We suggest that short-term health impact studies could
complement lifetime mortality studies to better characterize the lifetime
impacts of disease on both individuals and populations.


  The contact process is an emblematic model of a non-equilibrium system,
containing a phase transition between inactive and active dynamical regimes. In
the epidemiological context, the model is known as the
susceptible-infected-susceptible (SIS) model, and widely used to describe
contagious spreading. In this work, we demonstrate how accurate and efficient
representations of the full probability distribution over all configurations of
the contact process on a one-dimensional chain can be obtained by means of
Matrix Product States (MPS). We modify and adapt MPS methods from many-body
quantum systems to study the classical distributions of the driven contact
process at late times. We give accurate and efficient results for the
distribution of large gaps, and illustrate the advantage of our methods over
Monte Carlo simulations. Furthermore, we study the large deviation statistics
of the dynamical activity, defined as the total number of configuration changes
along a trajectory, and investigate quantum-inspired entropic measures, based
on the second R\'enyi entropy.


  We consider a cell population subject to a parasite infection. Cells divide
at a constant rate and, at division, share the parasites they contain between
their two daughter cells. The sharing may be asymmetrical, and its law may
depend on the quantity of parasites in the mother. Cells die at a rate which
may depend on the quantity of parasites they carry, and are also killed when
this quantity explodes. We study the survival of the cell population as well as
the mean quantity of parasites in the cells, and focus on the role of the
parasites partitioning kernel at division.


  The recognition of the importance of drug-like properties beyond potency to
reduce clinical attrition of biologics has driven significant progress in the
development of in vitro and in silico tools for developability assessment of
antibody sequences. It is now routine to identify and eliminate or optimize
antibody hits with poor developability profiles. To further accelerate
discovery timelines and reduce clinical and non-clinical development attrition
rates, more proactive in silico approaches to design sequence spaces with
favorable developability profiles are required. From pragmatically
front-loading structure based drug design for developability, to combining next
generation sequencing with machine learning to shape screening libraries, to
adapting the use of artificial intelligence and deep learning for
immunoglobulins, we review herein progressively more proactive approaches to
developability by design.


  Deep generative models have recently achieved superior performance in 3D
molecule generation. Most of them first generate atoms and then add chemical
bonds based on the generated atoms in a post-processing manner. However, there
might be no corresponding bond solution for the temporally generated atoms as
their locations are generated without considering potential bonds. We define
this problem as the atom-bond inconsistency problem and claim it is the main
reason for current approaches to generating unrealistic 3D molecules. To
overcome this problem, we propose a new diffusion model called MolDiff which
can generate atoms and bonds simultaneously while still maintaining their
consistency by explicitly modeling the dependence between their relationships.
We evaluated the generation ability of our proposed model and the quality of
the generated molecules using criteria related to both geometry and chemical
properties. The empirical studies showed that our model outperforms previous
approaches, achieving a three-fold improvement in success rate and generating
molecules with significantly better quality.


  We consider the SIRWJS epidemiological model that includes the waning and
boosting of immunity via secondary infections. We carry out combined analytical
and numerical investigations of the dynamics. The formulae describing the
existence and stability of equilibria are derived. Combining this analysis with
numerical continuation techniques, we construct global bifurcation diagrams
with respect to several epidemiological parameters. The bifurcation analysis
reveals a very rich structure of possible global dynamics. We show that
backward bifurcation is possible at the critical value of the basic
reproduction number, $\mathcal{R}_0 = 1$. Furthermore, we find stability
switches and Hopf bifurcations from steady states forming multiple endemic
bubbles, and saddle-node bifurcations of periodic orbits. Regions of
bistability are also found, where either two stable steady states, or a stable
steady state and a stable periodic orbit coexist. This work provides an insight
to the rich and complicated infectious disease dynamics that can emerge from
the waning and boosting of immunity.


  The internal behaviour of a population is an important feature to take
account of when modelling their dynamics. In line with kin selection theory,
many social species tend to cluster into distinct groups in order to enhance
their overall population fitness. Temporal interactions between populations are
often modelled using classical mathematical models, but these sometimes fail to
delve deeper into the, often uncertain, relationships within populations. Here,
we introduce a stochastic framework that aims to capture the interactions of
animal groups and an auxiliary population over time. We demonstrate the model's
capabilities, from a Bayesian perspective, through simulation studies and by
fitting it to predator-prey count time series data. We then derive an
approximation to the group correlation structure within such a population,
while also taking account of the effect of the auxiliary population. We finally
discuss how this approximation can lead to ecologically realistic
interpretations in a predator-prey context. This approximation also serves as
verification to whether the population in question satisfies our various
assumptions. Our modelling approach will be useful for empiricists for
monitoring groups within a conservation framework and also theoreticians
wanting to quantify interactions, to study cooperation and other phenomena
within social populations.


  Computational and machine learning approaches to model the conformational
landscape of macrocyclic peptides have the potential to enable rational design
and optimization. However, accurate, fast, and scalable methods for modeling
macrocycle geometries remain elusive. Recent deep learning approaches have
significantly accelerated protein structure prediction and the generation of
small-molecule conformational ensembles, yet similar progress has not been made
for macrocyclic peptides due to their unique properties. Here, we introduce
CREMP, a resource generated for the rapid development and evaluation of machine
learning models for macrocyclic peptides. CREMP contains 36,198 unique
macrocyclic peptides and their high-quality structural ensembles generated
using the Conformer-Rotamer Ensemble Sampling Tool (CREST). Altogether, this
new dataset contains nearly 31.3 million unique macrocycle geometries, each
annotated with energies derived from semi-empirical extended tight-binding
(xTB) DFT calculations. We anticipate that this dataset will enable the
development of machine learning models that can improve peptide design and
optimization for novel therapeutics.


  Protein-protein interactions (PPIs) are crucial in various biological
processes and their study has significant implications for drug development and
disease diagnosis. Existing deep learning methods suffer from significant
performance degradation under complex real-world scenarios due to various
factors, e.g., label scarcity and domain shift. In this paper, we propose a
self-ensembling multigraph neural network (SemiGNN-PPI) that can effectively
predict PPIs while being both efficient and generalizable. In SemiGNN-PPI, we
not only model the protein correlations but explore the label dependencies by
constructing and processing multiple graphs from the perspectives of both
features and labels in the graph learning process. We further marry GNN with
Mean Teacher to effectively leverage unlabeled graph-structured PPI data for
self-ensemble graph learning. We also design multiple graph consistency
constraints to align the student and teacher graphs in the feature embedding
space, enabling the student model to better learn from the teacher model by
incorporating more relationships. Extensive experiments on PPI datasets of
different scales with different evaluation settings demonstrate that
SemiGNN-PPI outperforms state-of-the-art PPI prediction methods, particularly
in challenging scenarios such as training with limited annotations and testing
on unseen data.


  The genetic information that dictates the structure and function of all life
forms is encoded in the DNA. In 1953, Watson and Crick first presented the
double helical structure of a DNA molecule. Their findings unearthed the desire
to elucidate the exact composition and sequence of DNA molecules. Discoveries
and the subsequent development and optimization of techniques that allowed for
deciphering the DNA sequence has opened new doors in research, biotech, and
healthcare. The application of high-throughput sequencing technologies in these
industries has positively impacted and will continue to contribute to the
betterment of humanity and the global economy. Improvements, such as the use of
radioactive molecules for DNA sequencing to the use of florescent dyes and the
implementation of polymerase chain reaction (PCR) for amplification, led to
sequencing a few hundred base pairs in days, to automation, where sequencing of
thousands of base pairs in hours became possible. Significant advances have
been made, but there is still room for improvement. Here, we look at the
history and the technology of the currently available next-generation
sequencing platforms and the possible applications of such technologies to
biomedical research and beyond.


  Climate changes may affect ecosystems destabilising relationships among
species. We investigate the spatial rock-paper-scissors models with a regional
unevenness that reduces the selection capacity of organisms of one species. Our
results show that the regionally weak species predominates in the local
ecosystem, while spiral patterns appear far from the region, where individuals
of every species play the rock-paper-scissors game with the same strength.
Because the weak species controls all local territory, it is attractive for the
other species to enter the local ecosystem to conquer the territory. However,
our stochastic simulations show that the transitory waves formed when organisms
of the strong species reach the region are quickly destroyed because of local
strength unbalance in the selection game rules. Computing the effect of the
topology on population dynamics, we find that the prevalence of the weak
species becomes more significant if the transition of the selection capacity to
the area of uneven rock-paper-scissors rules is smooth. Finally, our findings
show that the biodiversity loss due to the arising of regional unevenness is
minimised if the transition to the region where the cyclic game is unbalanced
is abrupt. Our results may be helpful to biologists in comprehending the
consequences of changes in the environmental conditions on species coexistence
and spatial patterns in complex systems.


  Deep learning-based approaches, such as AlphaFold2 (AF2), have significantly
advanced protein tertiary structure prediction, achieving results comparable to
real biological experimental methods. While AF2 has shown limitations in
predicting the effects of mutations, its robustness against sequence mutations
remains to be determined. Starting with the wild-type (WT) sequence, we
investigate adversarial sequences generated via an evolutionary approach, which
AF2 predicts to be substantially different from WT. Our experiments on CASP14
reveal that by modifying merely three residues in the protein sequence using a
combination of replacement, deletion, and insertion strategies, the alteration
in AF2's predictions, as measured by the Local Distance Difference Test (lDDT),
reaches 46.61. Moreover, when applied to a specific protein, SPNS2, our
proposed algorithm successfully identifies biologically meaningful residues
critical to protein structure determination and potentially indicates
alternative conformations, thus significantly expediting the experimental
process.


  Many physical and biological systems rely on the progression of material
through multiple independent stages. In viral replication, for example, virions
enter a cell to undergo a complex process comprising several disparate stages
before the eventual accumulation and release of replicated virions. While such
systems may have some control over the internal dynamics that make up this
progression, a challenge for many is to regulate behaviour under what are often
highly variable external environments acting as system inputs. In this work, we
study a simple analogue of this problem through a linear multicompartment model
subject to a stochastic input in the form of a mean-reverting
Ornstein-Uhlenbeck process, a type of Gaussian process. By expressing the
system as a multidimensional Gaussian process, we derive several closed-form
analytical results relating to the covariances and autocorrelations of the
system, quantifying the smoothing effect discrete compartments afford
multicompartment systems. Semi-analytical results demonstrate that feedback and
feedforward loops can enhance system robustness, and simulation results probe
the intractable problem of the first passage time distribution, which has
specific relevance to eventual cell lysis in the viral replication cycle.
Finally, we demonstrate that the smoothing seen in the process is a consequence
of the discreteness of the system, and does not manifest in an equivalent
continuum limit description. While we make progress through analysis of a
simple linear problem, many of our insights are applicable more generally, and
our work enables future analysis into multicompartment processes subject to
stochastic inputs.


  In both natural and artificial studies, evolution is often seen as synonymous
to natural selection. Individuals evolve under pressures set by environments
that are either reset or do not carry over significant changes from previous
generations. Thus, niche construction (NC), the reciprocal process to natural
selection where individuals incur inheritable changes to their environment, is
ignored. Arguably due to this lack of study, the dynamics of NC are today
little understood, especially in real-world settings. In this work, we study NC
in simulation environments that consist of multiple, diverse niches and
populations that evolve their plasticity, evolvability and niche-constructing
behaviors. Our empirical analysis reveals many interesting dynamics, with
populations experiencing mass extinctions, arms races and oscillations. To
understand these behaviors, we analyze the interaction between NC and
adaptability and the effect of NC on the population's genomic diversity and
dispersal, observing that NC diversifies niches. Our study suggests that
complexifying the simulation environments studying NC, by considering multiple
and diverse niches, is necessary for understanding its dynamics and can lend
testable hypotheses to future studies of both natural and artificial systems.


  Antibodies are crucial proteins produced by the immune system in response to
foreign substances or antigens. The specificity of an antibody is determined by
its complementarity-determining regions (CDRs), which are located in the
variable domains of the antibody chains and form the antigen-binding site.
Previous studies have utilized complex techniques to generate CDRs, but they
suffer from inadequate geometric modeling. Moreover, the common iterative
refinement strategies lead to an inefficient inference. In this paper, we
propose a \textit{simple yet effective} model that can co-design 1D sequences
and 3D structures of CDRs in a one-shot manner. To achieve this, we decouple
the antibody CDR design problem into two stages: (i) geometric modeling of
protein complex structures and (ii) sequence-structure co-learning. We develop
a novel macromolecular structure invariant embedding, typically for protein
complexes, that captures both intra- and inter-component interactions among the
backbone atoms, including C$\alpha$, N, C, and O atoms, to achieve
comprehensive geometric modeling. Then, we introduce a simple cross-gate MLP
for sequence-structure co-learning, allowing sequence and structure
representations to implicitly refine each other. This enables our model to
design desired sequences and structures in a one-shot manner. Extensive
experiments are conducted to evaluate our results at both the sequence and
structure levels, which demonstrate that our model achieves superior
performance compared to the state-of-the-art antibody CDR design methods.


  A central task in computational drug discovery is to construct models from
known active molecules to find further promising molecules for subsequent
screening. However, typically only very few active molecules are known.
Therefore, few-shot learning methods have the potential to improve the
effectiveness of this critical phase of the drug discovery process. We
introduce a new method for few-shot drug discovery. Its main idea is to enrich
a molecule representation by knowledge about known context or reference
molecules. Our novel concept for molecule representation enrichment is to
associate molecules from both the support set and the query set with a large
set of reference (context) molecules through a Modern Hopfield Network.
Intuitively, this enrichment step is analogous to a human expert who would
associate a given molecule with familiar molecules whose properties are known.
The enrichment step reinforces and amplifies the covariance structure of the
data, while simultaneously removing spurious correlations arising from the
decoration of molecules. Our approach is compared with other few-shot methods
for drug discovery on the FS-Mol benchmark dataset. On FS-Mol, our approach
outperforms all compared methods and therefore sets a new state-of-the art for
few-shot learning in drug discovery. An ablation study shows that the
enrichment step of our method is the key to improve the predictive quality. In
a domain shift experiment, we further demonstrate the robustness of our method.
Code is available at https://github.com/ml-jku/MHNfs.


  Image-based profiling techniques have become increasingly popular over the
past decade for their applications in target identification,
mechanism-of-action inference, and assay development. These techniques have
generated large datasets of cellular morphologies, which are typically used to
investigate the effects of small molecule perturbagens. In this work, we extend
the impact of such dataset to improving quantitative structure-activity
relationship (QSAR) models by introducing Molecule-Morphology Contrastive
Pretraining (MoCoP), a framework for learning multi-modal representation of
molecular graphs and cellular morphologies. We scale MoCoP to approximately
100K molecules and 600K morphological profiles using data from the JUMP-CP
Consortium and show that MoCoP consistently improves performances of graph
neural networks (GNNs) on molecular property prediction tasks in ChEMBL20
across all dataset sizes. The pretrained GNNs are also evaluated on internal
GSK pharmacokinetic data and show an average improvement of 2.6% and 6.3% in
AUPRC for full and low data regimes, respectively. Our findings suggest that
integrating cellular morphologies with molecular graphs using MoCoP can
significantly improve the performance of QSAR models, ultimately expanding the
deep learning toolbox available for QSAR applications.


  To facilitate a rapid response to pandemic threats, this paper focuses on
developing a mechanistic simulation model for in vitro transcription (IVT)
process, a crucial step in mRNA vaccine manufacturing. To enhance production
and support industry 4.0, this model is proposed to improve the prediction and
analysis of IVT enzymatic reaction network. It incorporates a novel stochastic
molecular reaction queueing network with a regulatory kinetic model
characterizing the effect of bioprocess state variables on reaction rates. The
empirical study demonstrates that the proposed model has a promising
performance under different production conditions and it could offer potential
improvements in mRNA product quality and yield.


  As Nature's version of machine learning, evolution has solved many
extraordinarily complex problems, none perhaps more remarkable than learning to
harness an increase in chemical entropy (disorder) to generate directed
chemical forces (order). Using muscle as a model system, here I unpack the
basic mechanism by which life creates order from disorder. In short, evolution
tuned the physical properties of certain proteins to contain changes in
chemical entropy. As it happens, these are the "sensible" properties Gibbs
postulated were needed to solve his paradox.


  Generative models for molecules based on sequential line notation (e.g.
SMILES) or graph representation have attracted an increasing interest in the
field of structure-based drug design, but they struggle to capture important 3D
spatial interactions and often produce undesirable molecular structures. To
address these challenges, we introduce Lingo3DMol, a pocket-based 3D molecule
generation method that combines language models and geometric deep learning
technology. A new molecular representation, fragment-based SMILES with local
and global coordinates, was developed to assist the model in learning molecular
topologies and atomic spatial positions. Additionally, we trained a separate
noncovalent interaction predictor to provide essential binding pattern
information for the generative model. Lingo3DMol can efficiently traverse
drug-like chemical spaces, preventing the formation of unusual structures. The
Directory of Useful Decoys-Enhanced (DUD-E) dataset was used for evaluation.
Lingo3DMol outperformed state-of-the-art methods in terms of drug-likeness,
synthetic accessibility, pocket binding mode, and molecule generation speed.


  Phylogenetic networks are used to represent the evolutionary history of
species. They are versatile when compared to traditional phylogenetic trees, as
they capture more complex evolutionary events such as hybridization and
horizontal gene transfer. Distance-based methods such as the Neighbor-Net
algorithm are widely used to compute phylogenetic networks from data. However,
the output is necessarily an undirected graph, posing a great challenge to
deduce the direction of genetic flow in order to infer the true evolutionary
history. Recently, Huber et al. investigated two different computational
problems relevant to orienting undirected phylogenetic networks into directed
ones. In this paper, we consider the problem of orienting an undirected binary
network into a tree-child network. We give some necessary conditions for
determining the tree-child orientability, such as a tight upper bound on the
size of tree-child orientable graphs, as well as many interesting examples. In
addition, we introduce new families of undirected phylogenetic networks, the
jellyfish graphs and ladder graphs, that are orientable but not tree-child
orientable. We also prove that any ladder graph can be made tree-child
orientable by adding extra leaves, and describe a simple algorithm for
orienting a ladder graph to a tree-child network with the minimum number of
extra leaves. We pose many open problems as well.


  Due to their high predictive performance and flexibility, machine learning
models are an appropriate and efficient tool for ecologists. However,
implementing a machine learning model is not yet a trivial task and may seem
intimidating to ecologists with no previous experience in this area. Here we
provide a series of tips to help ecologists in implementing machine learning
models. We focus on classification problems as many ecological studies aim to
assign data into predefined classes such as ecological states or biological
entities. Each of the nine tips identifies a common error, trap or challenge in
developing machine learning models and provides recommendations to facilitate
their use in ecological studies.


  Functional connectivity (FC) is one of the most common inputs to fMRI-based
predictive models, due to a combination of its simplicity and robustness.
However, there may be a lack of theoretical models for the generation of FC. In
this work, we present a straightforward decomposition of FC into a set of basis
states of sine waves with an additional jitter component. We show that the
decomposition matches the predictive ability of FC after including 5-10 bases.
We also find that both the decomposition and its residual have approximately
equal predictive value, and when combined into an ensemble, exceed the AUC of
FC-based prediction by up to 5%. Additionally, we find the residual can be used
for subject fingerprinting, with 97.3% same-subject, different-scan
identifiability, compared to 62.5% for FC. Unlike PCA or Factor Analysis
methods, our method does not require knowledge of a population to perform its
decomposition; a single subject is enough. Our decomposition of FC into two
equally-predictive components may lead to a novel appreciation of group
differences in patient populations. Additionally, we generate synthetic patient
FC based on user-specified characteristics such as age, sex, and disease
diagnosis. By creating synthetic datasets or augmentations we may reduce the
high financial burden associated with fMRI data acquisition.


  The graph traversal edit distance (GTED), introduced by Ebrahimpour Boroojeny
et al.~(2018), is an elegant distance measure defined as the minimum edit
distance between strings reconstructed from Eulerian trails in two edge-labeled
graphs. GTED can be used to infer evolutionary relationships between species by
comparing de Bruijn graphs directly without the computationally costly and
error-prone process of genome assembly. Ebrahimpour Boroojeny et al.~(2018)
propose two ILP formulations for GTED and claim that GTED is polynomially
solvable because the linear programming relaxation of one of the ILPs always
yields optimal integer solutions. The claim that GTED is polynomially solvable
is contradictory to the complexity results of existing string-to-graph matching
problems. We resolve this conflict in complexity results by proving that GTED
is NP-complete and showing that the ILPs proposed by Ebrahimpour Boroojeny et
al. do not solve GTED but instead solve for a lower bound of GTED and are not
solvable in polynomial time. In addition, we provide the first two, correct ILP
formulations of GTED and evaluate their empirical efficiency. These results
provide solid algorithmic foundations for comparing genome graphs and point to
the direction of heuristics. The source code to reproduce experimental results
is available at https://github.com/Kingsford-Group/gtednewilp/.


  Designing biological sequences is an important challenge that requires
satisfying complex constraints and thus is a natural problem to address with
deep generative modeling. Diffusion generative models have achieved
considerable success in many applications. Score-based generative stochastic
differential equations (SDE) model is a continuous-time diffusion model
framework that enjoys many benefits, but the originally proposed SDEs are not
naturally designed for modeling discrete data. To develop generative SDE models
for discrete data such as biological sequences, here we introduce a diffusion
process defined in the probability simplex space with stationary distribution
being the Dirichlet distribution. This makes diffusion in continuous space
natural for modeling discrete data. We refer to this approach as Dirchlet
diffusion score model. We demonstrate that this technique can generate samples
that satisfy hard constraints using a Sudoku generation task. This generative
model can also solve Sudoku, including hard puzzles, without additional
training. Finally, we applied this approach to develop the first human promoter
DNA sequence design model and showed that designed sequences share similar
properties with natural promoter sequences.


  This work describes a Tree Ring Detection method for complete Cross-Sections
of trees (CS-TRD). The method is based on the detection, processing, and
connection of edges corresponding to the tree's growth rings. The method
depends on the parameters for the Canny Devernay edge detector ($\sigma$ and
two thresholds), a resize factor, the number of rays, and the pith location.
The first five parameters are fixed by default. The pith location can be marked
manually or using an automatic pith detection algorithm. Besides the pith
localization, the CS-TRD method is fully automated and achieves an F-Score of
89\% in the UruDendro dataset (of Pinus Taeda) with a mean execution time of 17
seconds and of 97\% in the Kennel dataset (of Abies Alba) with an average
execution time 11 seconds.


  Motivation: Studies including more than one type of 'omics data sets are
becoming more prevalent. Integrating these data sets can be a way to solidify
findings and even to make new discoveries. However, integrating multi-omics
data sets is challenging. Typically, data sets are integrated by performing an
all-vs-all correlation analysis, where each feature of the first data set is
correlated to each feature of the second data set. However, all-vs-all
association testing produces unstructured results that are hard to interpret,
and involves potentially unnecessary hypothesis testing that reduces
statistical power due to false discovery rate (FDR) adjustment.
  Implementation: Here, we present the anansi framework, and accompanying R
package, as a way to improve upon all-vs-all association analysis. We take a
knowledge-based approach where external databases like KEGG are used to
constrain the all-vs-all association hypothesis space, only considering
pairwise associations that are a priori known to occur. This produces
structured results that are easier to interpret, and increases statistical
power by skipping unnecessary hypothesis tests. In this paper, we present the
anansi framework and demonstrate its application to learn metabolite-function
interactions in the context of host-microbe interactions. We further extend our
framework beyond pairwise association testing to differential association
testing, and show how anansi can be used to identify associations that differ
in strength or degree based on sample covariates such as case/control status.
  Availability: https://github.com/thomazbastiaanssen/anansi


  The SARS-CoV-2 pandemic has emphasised the importance of developing a
universal vaccine that can protect against current and future variants of the
virus. The present study proposes a novel conditional protein Language Model
architecture, called Vaxformer, which is designed to produce natural-looking
antigenicity-controlled SARS-CoV-2 spike proteins. We evaluate the generated
protein sequences of the Vaxformer model using DDGun protein stability measure,
netMHCpan antigenicity score, and a structure fidelity score with AlphaFold to
gauge its viability for vaccine development. Our results show that Vaxformer
outperforms the existing state-of-the-art Conditional Variational Autoencoder
model to generate antigenicity-controlled SARS-CoV-2 spike proteins. These
findings suggest promising opportunities for conditional Transformer models to
expand our understanding of vaccine design and their role in mitigating global
health challenges. The code used in this study is available at
https://github.com/aryopg/vaxformer .


  A biological transition from a state N to a state T is characterized by a
rearrangement of the gene expression profile in the system, quantitatively
measured through the differential expression of genes. In contrast, changes in
genetic pathways are usually evaluated by means of hypothesis testing schemes.
We introduce a quantitative measure in order to evaluate the average
fold-change of genetic pathways in biological transitions and applied it to
characterize in general grounds the transition from a normal tissue to a tumor.
Additionally, we study the transformation from primary to metastatic melanoma.
Gene expression data from the TCGA portal for 16 tumors and the Reactome
compilation of pathways are used for this purpose.


  When studying social dilemma games, a crucial question arises regarding the
impact of general heterogeneity on cooperation, which has been shown to have
positive effects in numerous studies. Here, we demonstrate that heterogeneity
in the contribution value for the focal Public Goods Game can jeopardize
cooperation. We show that there is an optimal contribution value in the
homogeneous case that most benefits cooperation depending on the lattice. In a
heterogeneous scenario, where strategy and contribution coevolve, cooperators
making contributions higher than the optimal value end up harming those who
contribute lower. This effect is notably detrimental to cooperation in the
square lattice with von Neumann neighborhood, while it can have no impact in
others lattices. Furthermore, in parameter regions where a higher-contributing
cooperator cannot normally survive alone, the exploitation of lower value
contribution cooperators allows their survival, resembling a parasitic
behavior. To obtain these results, we employed various distributions for the
contribution values in the initial condition and conducted Monte Carlo
simulations.


  Identifying molecules that exhibit some pre-specified properties is a
difficult problem to solve. In the last few years, deep generative models have
been used for molecule generation. Deep Graph Variational Autoencoders are
among the most powerful machine learning tools with which it is possible to
address this problem. However, existing methods struggle in capturing the true
data distribution and tend to be computationally expensive. In this work, we
propose RGCVAE, an efficient and effective Graph Variational Autoencoder based
on: (i) an encoding network exploiting a new powerful Relational Graph
Isomorphism Network; (ii) a novel probabilistic decoding component. Compared to
several state-of-the-art VAE methods on two widely adopted datasets, RGCVAE
shows state-of-the-art molecule generation performance while being
significantly faster to train.


  The microbiome-gut-brain-axis field is multidisciplinary, benefiting from the
expertise of microbiology, ecology, psychiatry, computational biology, and
epidemiology amongst other disciplines. As the field matures and moves beyond a
basic demonstration of its relevance, it is critical that study design and
analysis are robust and foster reproducibility.
  In this companion piece to Bugs as Features (part I), we present techniques
from adjacent and disparate fields to enrich and inform the analysis of
microbiome-gut-brain-axis data. Emerging techniques built specifically for the
microbiome-gut-brain axis are also demonstrated. All of these methods are
contextualised to inform several common challenges: how do we establish
causality? How can we integrate data from multiple 'omics techniques? How might
we account for the dynamicism of host-microbiome interactions?
  This perspective is offered to experienced and emerging microbiome scientists
alike, to assist with these questions and others, at the study conception,
design, analysis and interpretation stages of research.


  Finely-tuned enzymatic pathways control cellular processes, and their
dysregulation can lead to disease. Creating predictive and interpretable models
for these pathways is challenging because of the complexity of the pathways and
of the cellular and genomic contexts. Here we introduce Elektrum, a deep
learning framework which addresses these challenges with data-driven and
biophysically interpretable models for determining the kinetics of biochemical
systems. First, it uses in vitro kinetic assays to rapidly hypothesize an
ensemble of high-quality Kinetically Interpretable Neural Networks (KINNs) that
predict reaction rates. It then employs a novel transfer learning step, where
the KINNs are inserted as intermediary layers into deeper convolutional neural
networks, fine-tuning the predictions for reaction-dependent in vivo outcomes.
Elektrum makes effective use of the limited, but clean in vitro data and the
complex, yet plentiful in vivo data that captures cellular context. We apply
Elektrum to predict CRISPR-Cas9 off-target editing probabilities and
demonstrate that Elektrum achieves state-of-the-art performance, regularizes
neural network architectures, and maintains physical interpretability.


  The quality of the inferences we make from pathogen sequence data is
determined by the number and composition of pathogen sequences that make up the
sample used to drive that inference. However, there remains limited guidance on
how to best structure and power studies when the end goal is phylogenetic
inference. One question that we can attempt to answer with molecular data is
whether some people are more likely to transmit a pathogen than others. Here we
present an estimator to quantify differential transmission, as measured by the
ratio of reproductive numbers between people with different characteristics,
using transmission pairs linked by molecular data, along with a sample size
calculation for this estimator. We also provide extensions to our method to
correct for imperfect identification of transmission linked pairs,
overdispersion in the transmission process, and group imbalance. We validate
this method via simulation and provide tools to implement it in an R package,
phylosamp.


  This paper presents a real-time simulation involving ''protozoan-like'' cells
that evolve by natural selection in a physical 2D ecosystem. Selection pressure
is exerted via the requirements to collect mass and energy from the
surroundings in order to reproduce by cell-division. Cells do not have fixed
morphologies from birth; they can use their resources in construction projects
that produce functional nodes on their surfaces such as photoreceptors for
light sensitivity or flagella for motility. Importantly, these nodes act as
modular components that connect to the cell's control system via IO channels,
meaning that the evolutionary process can replace one function with another
while utilising pre-developed control pathways on the other side of the
channel. A notable type of node function is the adhesion receptors that allow
cells to bind together into multicellular structures in which individuals can
share resource and signal to one another. The control system itself is modelled
as an artificial neural network that doubles as a gene regulatory network,
thereby permitting the co-evolution of form and function in a single data
structure and allowing cell specialisation within multicellular groups.


  The competitive exclusion principle (CEP) is a fundamental concept in the
niche theory, which posits that the number of available resources constrains
the coexistence of species. While the CEP offers an intuitive explanation for
coexistence, it has been challenged by counterexamples observed in nature. One
prominent counterexample is the phytoplankton community, known as the paradox
of the plankton. Diverse phytoplankton species coexist in the ocean even though
they demand a limited number of resources. To shed light on this remarkable
biodiversity in large ecosystems quantitatively, we consider self-regulation
into the generalized MacArthur's consumer-resource model and study the relative
diversity, the number ratio between coexisting consumers and resource kinds. By
employing the cavity method and generating functional analysis, we analytically
show how the bounds of the relative diversity can exceed unity and its
dependency on the strength of the self-regulation. We confirm the analysis with
numerical simulations and reveal that the self-regulation suppresses the
emergence of dominant species, thereby fostering high biodiversity.
Furthermore, we study the effect of the self-regulation on different
environments and show that the effect relies on the environmental condition.
Our work presents a comprehensive framework within the niche theory that
encompasses the CEP and its counterexamples by introducing the role of
self-regulation.


  Designing new molecules is essential for drug discovery and material science.
Recently, deep generative models that aim to model molecule distribution have
made promising progress in narrowing down the chemical research space and
generating high-fidelity molecules. However, current generative models only
focus on modeling either 2D bonding graphs or 3D geometries, which are two
complementary descriptors for molecules. The lack of ability to jointly model
both limits the improvement of generation quality and further downstream
applications. In this paper, we propose a new joint 2D and 3D diffusion model
(JODO) that generates complete molecules with atom types, formal charges, bond
information, and 3D coordinates. To capture the correlation between molecular
graphs and geometries in the diffusion process, we develop a Diffusion Graph
Transformer to parameterize the data prediction model that recovers the
original data from noisy data. The Diffusion Graph Transformer interacts node
and edge representations based on our relational attention mechanism, while
simultaneously propagating and updating scalar features and geometric vectors.
Our model can also be extended for inverse molecular design targeting single or
multiple quantum properties. In our comprehensive evaluation pipeline for
unconditional joint generation, the results of the experiment show that JODO
remarkably outperforms the baselines on the QM9 and GEOM-Drugs datasets.
Furthermore, our model excels in few-step fast sampling, as well as in inverse
molecule design and molecular graph generation. Our code is provided in
https://github.com/GRAPH-0/JODO.


  The use of oncolytic viruses as cancer treatment has received considerable
attention in recent years, however the spatial dynamics of this viral infection
is still poorly understood. We present here a stochastic agent-based model
describing infected and uninfected cells for solid tumours, which interact with
viruses in the absence of an immune response. Two kinds of movement, namely
undirected random and pressure-driven movements, are considered: the continuum
limit of the models is derived and a systematic comparison between the systems
of partial differential equations and the individual-based model, in one and
two dimensions, is carried out.
  In the case of undirected movement, a good agreement between agent-based
simulations and the numerical and well-known analytical results for the
continuum model is possible. For pressure-driven motion, instead, we observe a
wide parameter range in which the infection of the agents remains confined to
the center of the tumour, even though the continuum model shows traveling waves
of infection; outcomes appear to be more sensitive to stochasticity and
uninfected regions appear harder to invade, giving rise to irregular,
unpredictable growth patterns.
  Our results show that the presence of spatial constraints in tumours'
microenvironments limiting free expansion has a very significant impact on
virotherapy. Outcomes for these tumours suggest a notable increase in
variability. All these aspects can have important effects when designing
individually tailored therapies where virotherapy is included.


  While visual and auditory information conveyed by wavelength of light and
frequency of sound have been decoded, predicting olfactory information encoded
by the combination of odorants remains challenging due to the unknown and
potentially discontinuous perceptual space of smells and odorants. Herein, we
develop a deep learning model called Mol-PECO (Molecular Representation by
Positional Encoding of Coulomb Matrix) to predict olfactory perception from
molecular structures. Mol-PECO updates the learned atom embedding by
directional graph convolutional networks (GCN), which model the Laplacian
eigenfunctions as positional encoding, and Coulomb matrix, which encodes atomic
coordinates and charges. With a comprehensive dataset of 8,503 molecules,
Mol-PECO directly achieves an area-under-the-receiver-operating-characteristic
(AUROC) of 0.813 in 118 odor descriptors, superior to the machine learning of
molecular fingerprints (AUROC of 0.761) and GCN of adjacency matrix (AUROC of
0.678). The learned embeddings by Mol-PECO also capture a meaningful odor space
with global clustering of descriptors and local retrieval of similar odorants.
Our work may promote the understanding and decoding of the olfactory sense and
mechanisms.


  There is a vast amount of literature concerning the appropriateness of
various perturbation parameters for the standard quasi-steady state
approximation in the Michaelis-Menten reaction mechanism, and also concerning
the relevance of these parameters for the accuracy of the approximation by the
familiar Michaelis-Menten equation. Typically, the arguments in the literature
are based on (heuristic) timescale estimates, from which one cannot obtain
reliable quantitative estimates for the error of the quasi-steady state
approximation. We take a different approach. By combining phase plane analysis
with differential inequalities, we derive sharp explicit upper and lower
estimates for the duration of the initial transient and substrate depletion
during this transitory phase. In addition, we obtain rigorous bounds on the
accuracy of the standard quasi-steady state approximation in the slow dynamics
regime. Notably, under the assumption that the quasi-steady state approximation
is valid over the entire time course of the reaction, our error estimate is of
order one in the Segel-Slemrod parameter.


  Cancer is a systemic heterogeneous disease involving complex molecular
networks. Tumor formation involves epithelial-mesenchymal transition (EMT),
which promotes both metastasis and plasticity of cancer cells. Recent
experiments proposed that cancer cells can be transformed into adipocytes with
combination drugs. However, the underlying mechanisms for how these drugs work
from molecular network perspective remain elusive. To reveal the mechanism of
cancer-adipose conversion (CAC), we adopt a systems biology approach by combing
mathematical modeling and molecular experiments based on the underlying
molecular regulatory network. We identified four types of attractors which
correspond to epithelial (E), mesenchymal (M), adipose (A) and
partial/intermediate EMT (P) cell states on the CAC landscape. Landscape and
transition path results illustrate that the intermediate states play critical
roles in cancer to adipose transition. Through a landscape control strategy, we
identified two new therapeutic strategies for drug combinations to promote CAC.
We further verified these predictions by molecular experiments in different
cell lines. Our combined computational and experimental approach provides a
powerful tool to explore molecular mechanisms for cell fate transitions in
cancer networks. Our results revealed the underlying mechanism for intermediate
cell states governing the CAC, and identified new potential drug combinations
to induce cancer adipogenesis.


  With the recent availability of tissue-specific gene expression data, e.g.,
provided by the GTEx Consortium, there is interest in comparing gene
co-expression patterns across tissues. One promising approach to this problem
is to use a multilayer network analysis framework and perform multilayer
community detection. Communities in gene co-expression networks reveal groups
of genes similarly expressed across individuals, potentially involved in
related biological processes responding to specific environmental stimuli or
sharing common regulatory variations. We construct a multilayer network in
which each of the four layers is an exocrine gland tissue-specific gene
co-expression network. We develop methods for multilayer community detection
with correlation matrix input and an appropriate null model. Our correlation
matrix input method identifies five groups of genes that are similarly
co-expressed in multiple tissues (a community that spans multiple layers, which
we call a generalist community) and two groups of genes that are co-expressed
in just one tissue (a community that lies primarily within just one layer,
which we call a specialist community). We further found gene co-expression
communities where the genes physically cluster across the genome significantly
more than expected by chance (on chromosomes 1 and 11). This clustering hints
at underlying regulatory elements determining similar expression patterns
across individuals and cell types. We suggest that KRTAP3-1, KRTAP3-3, and
KRTAP3-5 share regulatory elements in skin and pancreas. Furthermore, we find
that CELA3A and CELA3B share associated expression quantitative trait loci in
the pancreas. The results indicate that our multilayer community detection
method for correlation matrix input extracts biologically interesting
communities of genes.


  This paper reports on patterns exhibiting self-replication with spontaneous,
inheritable mutations and exponential genetic drift in Neural Cellular
Automata. Despite the models not being explicitly trained for mutation or
inheritability, the descendant patterns exponentially drift away from ancestral
patterns, even when the automaton is deterministic. While this is far from
being the first instance of evolutionary dynamics in a cellular automaton, it
is the first to do so by exploiting the power and convenience of Neural
Cellular Automata, arguably increasing the space of variations and the
opportunity for Open Ended Evolution.


  Searching for potential active compounds in large databases is a necessary
step to reduce time and costs in modern drug discovery pipelines. Such virtual
screening methods seek to provide predictions that allow the search space to be
narrowed down. Although cheminformatics has made great progress in exploiting
the potential of available big data, caution is needed to avoid introducing
bias and provide useful predictions with new compounds. In this work, we
propose the decision-support tool ALMERIA (Advanced Ligand Multiconformational
Exploration with Robust Interpretable Artificial Intelligence) for estimating
compound similarities and activity prediction based on pairwise molecular
contrasts while considering their conformation variability. The methodology
covers the entire pipeline from data preparation to model selection and
hyperparameter optimization. It has been implemented using scalable software
and methods to exploit large volumes of data -- in the order of several
terabytes -- , offering a very quick response even for a large batch of
queries. The implementation and experiments have been performed in a
distributed computer cluster using a benchmark, the public access DUD-E
database. In addition to cross-validation, detailed data split criteria have
been used to evaluate the models on different data partitions to assess their
true generalization ability with new compounds. Experiments show
state-of-the-art performance for molecular activity prediction (ROC AUC:
$0.99$, $0.96$, $0.87$), proving that the chosen data representation and
modeling have good properties to generalize. Molecular conformations --
prediction performance and sensitivity analysis -- have also been evaluated.
Finally, an interpretability analysis has been performed using the SHAP method.


  Generating desirable molecular structures in 3D is a fundamental problem for
drug discovery. Despite the considerable progress we have achieved, existing
methods usually generate molecules in atom resolution and ignore intrinsic
local structures such as rings, which leads to poor quality in generated
structures, especially when generating large molecules. Fragment-based molecule
generation is a promising strategy, however, it is nontrivial to be adapted for
3D non-autoregressive generations because of the combinational optimization
problems. In this paper, we utilize a coarse-to-fine strategy to tackle this
problem, in which a Hierarchical Diffusion-based model (i.e.~HierDiff) is
proposed to preserve the validity of local segments without relying on
autoregressive modeling. Specifically, HierDiff first generates coarse-grained
molecule geometries via an equivariant diffusion process, where each
coarse-grained node reflects a fragment in a molecule. Then the coarse-grained
nodes are decoded into fine-grained fragments by a message-passing process and
a newly designed iterative refined sampling module. Lastly, the fine-grained
fragments are then assembled to derive a complete atomic molecular structure.
Extensive experiments demonstrate that HierDiff consistently improves the
quality of molecule generation over existing methods


  Molecular biologists frequently interpret gene lists derived from
high-throughput experiments and computational analysis. This is typically done
as a statistical enrichment analysis that measures the over- or
under-representation of biological function terms associated with genes or
their properties, based on curated assertions from a knowledge base (KB) such
as the Gene Ontology (GO). Interpreting gene lists can also be framed as a
textual summarization task, enabling the use of Large Language Models (LLMs),
potentially utilizing scientific texts directly and avoiding reliance on a KB.
  We developed SPINDOCTOR (Structured Prompt Interpolation of Natural Language
Descriptions of Controlled Terms for Ontology Reporting), a method that uses
GPT models to perform gene set function summarization as a complement to
standard enrichment analysis. This method can use different sources of gene
functional information: (1) structured text derived from curated ontological KB
annotations, (2) ontology-free narrative gene summaries, or (3) direct model
retrieval.
  We demonstrate that these methods are able to generate plausible and
biologically valid summary GO term lists for gene sets. However, GPT-based
approaches are unable to deliver reliable scores or p-values and often return
terms that are not statistically significant. Crucially, these methods were
rarely able to recapitulate the most precise and informative term from standard
enrichment, likely due to an inability to generalize and reason using an
ontology. Results are highly nondeterministic, with minor variations in prompt
resulting in radically different term lists. Our results show that at this
point, LLM-based methods are unsuitable as a replacement for standard term
enrichment analysis and that manual curation of ontological assertions remains
necessary.


  Eigenvalue disparity, also known as timescale separation, permits the
systematic reduction of deterministic models of enzyme kinetics. Geometric
singular perturbation theory, of which eigenvalue disparity is central,
provides a coordinate-free framework for deriving reduced mass action models in
the deterministic realm. Moreover, homologous deterministic reductions are
often employed in stochastic models to reduce the computational complexity
required to simulate reactions with the Gillespie algorithm. Interestingly,
several detailed studies indicate that timescale separation does not always
guarantee the accuracy of reduced stochastic models. In this work, we examine
the roles of timescale separation and coordinate transformations in the
reduction of the Linear Noise Approximation (LNA) and, unlike previous studies,
we do not require the system to be comprised of distinct fast and slow
variables. Instead, we adopt a coordinate-free approach. We demonstrate that
eigenvalue disparity does not guarantee the accuracy of the reduced LNA, known
as the slow scale LNA (ssLNA). However, the inaccuracy of the ssLNA can often
be eliminated with a proper coordinate transformation. For planar systems in
separated (standard) form, we prove that the error between the variances of the
slow variable generated by the LNA and the ssLNA is $\mathcal{O}(\varepsilon)$.
We also address a nilpotent Jacobian scenario and use the blow-up method to
construct a reduced equation that is accurate near the singular limit in the
deterministic regime. However, this reduction in the stochastic regime is far
less accurate, which illustrates that eigenvalue disparity plays a central role
in stochastic model reduction.


  Intraspecific interactions are key drivers of population dynamics because
they establish relations between individual fitness and population density. The
component Allee effect is defined as a positive correlation between any fitness
component of a focal organism and population density, and it can lead to
positive density dependence in the population per capita growth rate. The
spatial structure is key to determining whether and to which extent a component
Allee effect will manifest at the demographic level because it determines how
individuals interact with one another. However, existing spatial models to
study the Allee effect impose a fixed spatial structure, which limits our
understanding of how a component Allee effect and the spatial dynamics jointly
determine the existence of demographic Allee effects. To fill this gap, we
introduce a spatially-explicit theoretical framework where spatial structure
and population dynamics are emergent properties of the individual-level
demographic rates. Depending on the intensity of the individual processes the
population exhibits a variety of spatial patterns that determine the
demographic-level by-products of an existing individual-level component Allee
effect. We find that aggregation increases population abundance and allows
populations to survive in harsher environments and at lower global population
densities when compared with uniformly distributed organisms. Moreover,
aggregation can prevent the component Allee effect from manifesting at the
population level or restrict it to the level of each independent group. These
results provide a mechanistic understanding of how component Allee effects
operate for different spatial population structures and show at the population
level. Our results contribute to better understanding population dynamics in
the presence of Allee effects and can potentially inform population management
strategies.


  Introduction: The COVID-19 pandemic highlighted the importance of making
epidemiological data and scientific insights easily accessible and explorable
for public health agencies, the general public, and researchers.
State-of-the-art approaches for sharing data and insights included regularly
updated reports and web dashboards. However, they face a trade-off between the
simplicity and flexibility of data exploration. With the capabilities of recent
large language models (LLMs) such as GPT-4, this trade-off can be overcome.
  Results: We developed the chatbot "GenSpectrum Chat"
(https://cov-spectrum.org/chat) which uses GPT-4 as the underlying large
language model (LLM) to explore SARS-CoV-2 genomic sequencing data. Out of 500
inputs from real-world users, the chatbot provided a correct answer for 453
prompts; an incorrect answer for 13 prompts, and no answer although the
question was within scope for 34 prompts. We also tested the chatbot with
inputs from 10 different languages, and despite being provided solely with
English instructions and examples, it successfully processed prompts in all
tested languages.
  Conclusion: LLMs enable new ways of interacting with information systems. In
the field of public health, GenSpectrum Chat can facilitate the analysis of
real-time pathogen genomic data. With our chatbot supporting interactive
exploration in different languages, we envision quick and direct access to the
latest evidence for policymakers around the world.


  Generating molecules with high binding affinities to target proteins (a.k.a.
structure-based drug design) is a fundamental and challenging task in drug
discovery. Recently, deep generative models have achieved remarkable success in
generating 3D molecules conditioned on the protein pocket. However, most
existing methods consider molecular generation for protein pockets
independently while neglecting the underlying connections such as
subpocket-level similarities. Subpockets are the local protein environments of
ligand fragments and pockets with similar subpockets may bind the same
molecular fragment (motif) even though their overall structures are different.
Therefore, the trained models can hardly generalize to unseen protein pockets
in real-world applications. In this paper, we propose a novel method DrugGPS
for generalizable structure-based drug design. With the biochemical priors, we
propose to learn subpocket prototypes and construct a global interaction graph
to model the interactions between subpocket prototypes and molecular motifs.
Moreover, a hierarchical graph transformer encoder and motif-based 3D molecule
generation scheme are used to improve the model's performance. The experimental
results show that our model consistently outperforms baselines in generating
realistic drug candidates with high affinities in challenging
out-of-distribution settings.


  Ordinary differential equation (ODE) models used in mathematical epidemiology
assume explicitly or implicitly large populations. For the study of infections
in a hospital this is an extremely restrictive assumption as typically a
hospital ward has a few dozen, or even fewer, patients. This work reframes a
well-known model used in the study of the spread of antibiotic-resistant
bacteria in hospitals, to consider the pathogen transmission dynamics in small
populations. In this vein, this paper proposes a Markov chain model to describe
the spread of a single bacterial species in a hospital ward where patients may
be free of bacteria or may carry bacterial strains that are either sensitive or
resistant to antimicrobial agents. We determine the probability law of the
\emph{exact} reproduction number ${\cal R}_{exact,0}$, which is here defined as
the random number of secondary infections generated by those patients who are
accommodated in a predetermined bed before a patient who is free of bacteria is
accommodated in this bed for the first time. Specifically, we decompose the
exact reproduction number ${\cal R}_{exact,0}$ into two contributions allowing
us to distinguish between infections due to the sensitive and the resistant
bacterial strains. Our methodology is mainly based on structured Markov chains
and the use of related matrix-analytic methods. This guarantees the
compatibility of the new, finite-population model, with large population models
present in the literature and takes full advantage, in its mathematical
analysis, of the intrinsic stochasticity.


  We introduce a broad class of spatial models to describe how spatially
heterogeneous populations live, die, and reproduce. Individuals are represented
by points of a point measure, whose birth and death rates can depend both on
spatial position and local population density, defined via the convolution of
the point measure with a nonnegative kernel. We pass to three different scaling
limits: an interacting superprocess, a nonlocal partial differential equation
(PDE), and a classical PDE. The classical PDE is obtained both by first scaling
time and population size to pass to the nonlocal PDE, and then scaling the
kernel that determines local population density; and also (when the limit is a
reaction-diffusion equation) by simultaneously scaling the kernel width,
timescale and population size in our individual based model. A novelty of our
model is that we explicitly model a juvenile phase: offspring are thrown off in
a Gaussian distribution around the location of the parent, and reach (instant)
maturity with a probability that can depend on the population density at the
location at which they land. Although we only record mature individuals, a
trace of this two-step description remains in our population models, resulting
in novel limits governed by a nonlinear diffusion. Using a lookdown
representation, we retain information about genealogies and, in the case of
deterministic limiting models, use this to deduce the backwards in time motion
of the ancestral lineage of a sampled individual. We observe that knowing the
history of the population density is not enough to determine the motion of
ancestral lineages in our model. We also investigate the behaviour of lineages
for three different deterministic models of a population expanding its range as
a travelling wave: the Fisher-KPP equation, the Allen-Cahn equation, and a
porous medium equation with logistic growth.


  Computational RNA design has broad applications across synthetic biology and
therapeutic development. Fundamental to the diverse biological functions of RNA
is its conformational flexibility, enabling single sequences to adopt a variety
of distinct 3D states. Currently, computational biomolecule design tasks are
often posed as inverse problems, where sequences are designed based on adopting
a single desired structural conformation. In this work, we propose gRNAde, a
geometric RNA design pipeline that operates on sets of 3D RNA backbone
structures to explicitly account for and reflect RNA conformational diversity
in its designs. We demonstrate the utility of gRNAde for improving native
sequence recovery over single-state approaches on a new large-scale 3D RNA
design dataset, especially for multi-state and structurally diverse RNAs. Our
code is available at https://github.com/chaitjo/geometric-rna-design


  Competing strategies in an evolutionary game model, or species in a
biosystem, can easily form a larger unit which protects them from the invasion
of an external actor. Such a defensive alliance may have two, three, four or
even more members. But how effective can be such formation against an
alternative group composed by other competitors? To address this question we
study a minimal model where a two-member and a four-member alliances fight in a
symmetric and balanced way. By presenting representative phase diagrams, we
systematically explore the whole parameter range which characterizes the inner
dynamics of the alliances and the intensity of their interactions. The group
formed by a pair, who can exchange their neighboring positions, prevail in the
majority of the parameter region. The rival quartet can only win if their inner
cyclic invasion rate is significant while the mixing rate of the pair is
extremely low. At specific parameter values, when neither of the alliances is
strong enough, new four-member solutions emerge where a
rock-paper-scissors-like trio is extended by the other member of the pair.
These new solutions coexist hence all six competitors can survive. The
evolutionary process is accompanied by serious finite-size effects which can be
mitigated by appropriately chosen prepared initial states.


  Phylogenetic networks are a flexible model of evolution that can represent
reticulate evolution and handle complex data. Tree-based networks, which are
phylogenetic networks that have a spanning tree with the same root and leaf-set
as the network itself, have been well studied. However, not all networks are
tree-based. Francis-Semple-Steel (2018) thus introduced several indices to
measure the deviation of rooted binary phylogenetic networks $N$ from being
tree-based, such as the minimum number $\delta^\ast(N)$ of additional leaves
needed to make $N$ tree-based, and the minimum difference $\eta^\ast(N)$
between the number of vertices of $N$ and the number of vertices of a subtree
of $N$ that shares the root and leaf set with $N$. Hayamizu (2021) has
established a canonical decomposition of almost-binary phylogenetic networks of
$N$, called the maximal zig-zag trail decomposition, which has many
implications including a linear time algorithm for computing $\delta^\ast(N)$.
The Maximum Covering Subtree Problem (MCSP) is the problem of computing
$\eta^\ast(N)$, and Davidov et al. (2022) showed that this can be solved in
polynomial time (in cubic time when $N$ is binary) by an algorithm for the
minimum cost flow problem. In this paper, under the assumption that $N$ is
almost-binary (i.e. each internal vertex has in-degree and out-degree at most
two), we show that $\delta^\ast(N)\leq \eta^\ast (N)$ holds, which is tight,
and give a characterisation of such phylogenetic networks $N$ that satisfy
$\delta^\ast(N)=\eta^\ast(N)$. Our approach uses the canonical decomposition of
$N$ and focuses on how the maximal W-fences (i.e. the forbidden subgraphs of
tree-based networks) are connected to maximal M-fences in the network $N$. Our
results introduce a new class of phylogenetic networks for which MCSP can be
solved in linear time, which can be seen as a generalisation of tree-based
networks.


  Recent studies have shown competitive performance in protein design that aims
to find the amino acid sequence folding into the desired structure. However,
most of them disregard the importance of predictive confidence, fail to cover
the vast protein space, and do not incorporate common protein knowledge. After
witnessing the great success of pretrained models on diverse protein-related
tasks and the fact that recovery is highly correlated with confidence, we
wonder whether this knowledge can push the limits of protein design further. As
a solution, we propose a knowledge-aware module that refines low-quality
residues. We also introduce a memory-retrieval mechanism to save more than 50\%
of the training time. We extensively evaluate our proposed method on the CATH,
TS50, and TS500 datasets and our results show that our Knowledge-Design method
outperforms the previous PiFold method by approximately 9\% on the CATH
dataset. Specifically, Knowledge-Design is the first method that achieves 60+\%
recovery on CATH, TS50 and TS500 benchmarks. We also provide additional
analysis to demonstrate the effectiveness of our proposed method. The code will
be publicly available.


  Is there a unified framework for graph-based retrosynthesis prediction?
Through analysis of full-, semi-, and non-template retrosynthesis methods, we
discovered that they strive to strike an optimal balance between combinability
and consistency: \textit{Should atoms be combined as motifs to simplify the
molecular editing process, or should motifs be broken down into atoms to reduce
the vocabulary and improve predictive consistency?}
  Recent works have studied several specific cases, while none of them explores
different combinability-consistency trade-offs. Therefore, we propose
MotifRetro, a dynamic motif editing framework for retrosynthesis prediction
that can explore the entire trade-off space and unify graph-based models.
MotifRetro comprises two components: RetroBPE, which controls the
combinability-consistency trade-off, and a motif editing model, where we
introduce a novel LG-EGAT module to dynamiclly add motifs to the molecule. We
conduct extensive experiments on USPTO-50K to explore how the trade-off affects
the model performance and finally achieve state-of-the-art performance.


  The regulation of various cellular processes heavily relies on the protein
complexes within a living cell, necessitating a comprehensive understanding of
their three-dimensional structures to elucidate the underlying mechanisms.
While neural docking techniques have exhibited promising outcomes in binary
protein docking, the application of advanced neural architectures to multimeric
protein docking remains uncertain. This study introduces SyNDock, an automated
framework that swiftly assembles precise multimeric complexes within seconds,
showcasing performance that can potentially surpass or be on par with recent
advanced approaches. SyNDock possesses several appealing advantages not present
in previous approaches. Firstly, SyNDock formulates multimeric protein docking
as a problem of learning global transformations to holistically depict the
placement of chain units of a complex, enabling a learning-centric solution.
Secondly, SyNDock proposes a trainable two-step SE(3) algorithm, involving
initial pairwise transformation and confidence estimation, followed by global
transformation synchronization. This enables effective learning for assembling
the complex in a globally consistent manner. Lastly, extensive experiments
conducted on our proposed benchmark dataset demonstrate that SyNDock
outperforms existing docking software in crucial performance metrics, including
accuracy and runtime. For instance, it achieves a 4.5% improvement in
performance and a remarkable millionfold acceleration in speed.


  Classically, endemic diseases are expected to display relatively stable,
predictable infection dynamics. Indeed, diseases like influenza show yearly
recurring infection waves that can be anticipated accurately enough to develop
and distribute new vaccines. In contrast, newly-emerging diseases may cause
more complex, unpredictable dynamics, like COVID-19 has demonstrated. Here we
show that complex infection dynamics can also occur in the endemic state of
seasonal diseases when including human behaviour. We implement human behaviour
as a feedback between incidence and disease mitigation and study the system as
an epidemiological oscillator driven by seasonality. When behaviour and
seasonality have a comparable impact, we find a rich structure in parameter and
state space with Arnold tongues, co-existing attractors, and chaos. Moreover,
we demonstrate that if a disease requires active mitigation, balancing costs of
mitigation and infections can lead societies right into this complex regime. We
observe indications of this when comparing past COVID-19 and influenza data to
model simulations. Our results challenge the intuition that endemicity implies
predictability and seasonal waves, and show that complex dynamics can dominate
even in the endemic phase.


  The Kolmogorov-Feller equation for the probability density of a Markov
process on a half-axis, which arises in important problems of biology, is
considered. This process consists of random jumps distributed according to
Laplace's law and a deterministic return to zero. It is shown that Green's
function for such an equation can be found both in the form of a series and in
explicit form for some ratios of the parameters. This allows one to explicitly
find solutions to the Kolmogorov-Feller equation for many initial data.


  In this work we present a systematic mathematical approximation scheme that
exposes the way that information, about the evolutionary forces of selection
and random genetic drift, is encoded in gene-frequency trajectories.
  We determine approximate, time-dependent, gene-frequency trajectory
statistics, assuming additive selection. We use the probability of fixation to
test and illustrate the approximation scheme introduced. For the case where the
strength of selection and the effective population size have constant values,
we show how a standard result for the probability of fixation, under the
diffusion approximation, systematically emerges, when increasing numbers of
approximate trajectory statistics are taken into account. We then provide
examples of how time-dependent parameters influence gene-frequency statistics.


  Colorectal cancer is a public health problem.Approximately 30 to 50 \% of
colorectal tumors are caused by mutations in the KRAS gene.These mutations
induce uncontrolled proliferation.To date,There is no approved effective
treatment for the mutated KRAS oncogene.Farnesyltransferase (FTI) inhibitors
are considered a therapeutic target against the mutated KRAS
oncogene.Tipifarnib is a farnesyltransferase inhibitor that was analyzed in a
Phase II trial.In the present study, the three-dimensional structure of
farnesyltransferase complexed with tipifarnib [1SA4] was used as a basis to
exploit the characteristics of tipifarnib.A pharmacophore model was generated
based on the structure using the Asinex (Gold and Platinum Collections)
database.A total of 299 molecules were obtained after screening.The 299
molecules were anchored to the tipifarnib binding site in the
farnesyltransferase crystal structure for docking analysis.During the molecular
docking process, the pharmacophore that was modeled, and was used as a
constraint to eliminate the molecules that do not satisfy the
pharmacophore.Finally, four Hits identified as farnesyltransferase inhibitors
for biological tests.
  Keywords: colorectal cancer, structure-based pharmacophore, molecular
docking, KRAS, farnesyltransferase inhibitors, Virtual Screening.