  We consider the problem of efficiently simulating stochastic models of
chemical kinetics. The Gillespie Stochastic Simulation algorithm (SSA) is often
used to simulate these models, however, in many scenarios of interest, the
computational cost quickly becomes prohibitive. This is further exasperated in
the Bayesian inference context when estimating parameters of chemical models,
as the intractability of the likelihood requires multiple simulations of the
underlying system. To deal with issues of computational complexity in this
paper, we propose a novel hybrid $\tau$-leap algorithm for simulating
well-mixed chemical systems. In particular, the algorithm uses $\tau$-leap when
appropriate (high population densities), and SSA when necessary (low population
densities, when discrete effects become non-negligible). In the intermediate
regime, a combination of the two methods, which leverages the properties of the
underlying Poisson formulation, is employed. As illustrated through a number of
numerical experiments the hybrid $\tau$ offers significant computational
savings when compared to SSA without however sacrificing the overall accuracy.
This feature is particularly welcomed in the Bayesian inference context, as it
allows for parameter estimation of stochastic chemical kinetics at reduced
computational cost.


  The task of deducing three-dimensional molecular configurations from their
two-dimensional graph representations holds paramount importance in the fields
of computational chemistry and pharmaceutical development. The rapid
advancement of machine learning, particularly within the domain of deep
generative networks, has revolutionized the precision of predictive modeling in
this context. Traditional approaches often adopt a two-step strategy: initially
estimating interatomic distances and subsequently refining the spatial
molecular structure by solving a distance geometry problem. However, this
sequential approach occasionally falls short in accurately capturing the
intricacies of local atomic arrangements, thereby compromising the fidelity of
the resulting structural models. Addressing these limitations, this research
introduces a cutting-edge generative framework named \method{}. This framework
is grounded in the principles of diffusion observed in classical
non-equilibrium thermodynamics. \method{} views atoms as discrete entities and
excels in guiding the reversal of diffusion, transforming a distribution of
stochastic noise back into coherent molecular structures through a process akin
to a Markov chain. This transformation commences with the initial
representation of a molecular graph in an abstract latent space, culminating in
the realization of three-dimensional structures via a sophisticated bilevel
optimization scheme meticulously tailored to meet the specific requirements of
the task. One of the formidable challenges in this modeling endeavor involves
preserving roto-translational invariance to ensure that the generated molecular
conformations adhere to the laws of physics. Extensive experimental evaluations
confirm the efficacy of the proposed \method{} in comparison to
state-of-the-art methods.


  The effect of coastal upwelling on the recruitment and connectivity of
coastal marine populations has rarely been characterized to a level of detail
to be included into sound fishery management strategies. The gooseneck barnacle
(Pollicipes pollicipes) fishery at the Cantabrian Coast (Northern Spain) is
located at the fringes of the NW Spanish Upwelling system. This fishery is
being co-managed through a fine-scale, interspersed set of protected rocks
where each rock receives a distinct level of protection. Such interspersion is
potentially beneficial, but the extent to which such spacing is consistent with
mean larval dispersal distances is as yet unknown. We have simulated the spread
of gooseneck barnacle larvae in the Central Cantabrian Coast using a
high-resolution time-series of current profiles measured at a nearshore
location. During a year of high upwelling activity (2009), theoretical
recruitment success was 94% with peak recruitment predicted 56 km west of the
emission point. However, for a year of low upwelling activity (2011)
theoretical recruitment success dropped to 15.4% and peak recruitment was
expected 13 km east of the emission point. This is consistent with a positive
correlation between catch rates and the Integrated Upwelling Index, using a
4-year lag to allow recruits to reach commercial size. Furthermore, a net
long-term westward larval transport was estimated by means of mitochondrial
cytochrome c oxidase subunit I (COI) sequences for five populations in the
Cantabrian Sea. Our results call into question the role of long distance
dispersal, driven by the mesoscale processes in the area, in gooseneck barnacle
populations and point to the prevalent role of small-scale, asymmetric
connectivity more consistent with the typical scale of the co-management
process in this fishery.


  We propose a novel definition of life in terms of which its emergence in the
universe is expected, and its ever-creative open-ended evolution is entailed by
no law. Living organisms are Kantian Wholes that achieve Catalytic Closure,
Constraint Closure, and Spatial Closure. We here unite for the first time two
established mathematical theories, namely Collectively Autocatalytic Sets and
the Theory of the Adjacent Possible. The former establishes that a first-order
phase transition to molecular reproduction is expected in the chemical
evolution of the universe where the diversity and complexity of molecules
increases; the latter posits that, under loose hypotheses, if the system starts
with a small number of beginning molecules, each of which can combine with
copies of itself or other molecules to make new molecules, over time the number
of kinds of molecules increases slowly but then explodes upward hyperbolically.
Together these theories imply that life is expected as a phase transition in
the evolving universe. The familiar distinction between software and hardware
loses its meaning in living cells. We propose new ways to study the phylogeny
of metabolisms, new astronomical ways to search for life on exoplanets, new
experiments to seek the emergence of the most rudimentary life, and the hint of
a coherent testable pathway to prokaryotes with template replication and
coding.


  Correlation is not causation. As simple as this widely agreed-upon statement
may seem, scientifically defining causality and using it to drive our modern
biomedical research is immensely challenging. In this perspective, we attempt
to synergise the partly disparate fields of systems biology, causal reasoning,
and machine learning, to inform future approaches in the field of systems
biology and molecular networks.


  Aotearoa New Zealand implemented a Covid-19 elimination strategy in 2020 and
2021, which enabled a large majority of the population to be vaccinated before
being exposed to the virus. This strategy delivered one of the lowest pandemic
mortality rates in the world. However, quantitative estimates of the
population-level health benefits of vaccination are lacking. Here, we use a
validated mathematical model to investigate counterfactual scenarios with
differing levels of vaccine coverage in different age and ethnicity groups. The
model builds on earlier research by adding age- and time-dependent case
ascertainment, the effect of antiviral medications, improved hospitalisation
rate estimates, and the impact of relaxing control measures. The model was used
for scenario analysis and policy advice for the New Zealand Government in 2022
and 2023. We compare the number of Covid-19 hospitalisations, deaths, and years
of life lost in each counterfactual scenario to a baseline scenario that is
fitted to epidemiological data between January 2022 and June 2023. Our results
estimate that vaccines saved 6650 (95% credible interval [4424, 10180]) lives,
and prevented 74500 [51000, 115400] years of life lost and 45100 [34400, 55600]
hospitalisations during this 18-month period. Making the same comparison before
the benefit of antiviral medications is accounted for, the estimated number of
lives saved by vaccines increases to 7604 [5080, 11942]. Due to inequities in
the vaccine rollout, vaccination rates among M\=aori were lower than in people
of European ethnicity. Our results show that, if vaccination rates had been
equitable, an estimated 11-26% of the 292 M\=aori Covid-19 deaths that were
recorded in this time period could have been prevented. We conclude that
Covid-19 vaccination greatly reduced health burden in New Zealand and that
equity needs to be a key focus of future vaccination programmes.


  A rational design of new therapeutic drugs aims to find a molecular structure
with desired biological functionality, e.g., an ability to activate or suppress
a specific protein via binding to it. Molecular docking is a common technique
for evaluating protein-molecule interactions. Recently, Reinforcement Learning
(RL) has emerged as a promising approach to generating molecules with the
docking score (DS) as a reward. In this work, we reproduce, scrutinize and
improve the recent RL model for molecule generation called FREED
(arXiv:2110.01219). Extensive evaluation of the proposed method reveals several
limitations and challenges despite the outstanding results reported for three
target proteins. Our contributions include fixing numerous implementation bugs
and simplifying the model while increasing its quality, significantly extending
experiments, and conducting an accurate comparison with current
state-of-the-art methods for protein-conditioned molecule generation. We show
that the resulting fixed model is capable of producing molecules with superior
docking scores compared to alternative approaches.


  Biopolymer topology is critical for determining interactions inside cell
environments, exemplified by DNA where its response to mechanical perturbation
is as important as biochemical properties to its cellular roles. The dynamic
structures of chiral biopolymers exhibit complex dependence with extension and
torsion, however the physical mechanisms underpinning the emergence of
structural motifs upon physiological twisting and stretching are poorly
understood due to technological limitations in correlating force, torque and
spatial localization information. We present COMBI-Tweez (Combined Optical and
Magnetic BIomolecule TWEEZers), a transformative tool that overcomes these
challenges by integrating optical trapping, time-resolved electromagnetic
tweezers, and fluorescence microscopy, demonstrated on single DNA molecules,
that can controllably form and visualise higher order structural motifs
including plectonemes. This technology combined with cutting-edge MD
simulations provides quantitative insight into complex dynamic structures
relevant to DNA cellular processes and can be adapted to study a range of
filamentous biopolymers.


  Predicting interactions between proteins is one of the most important yet
challenging problems in structural bioinformatics. Intrinsically, potential
function sites in protein surfaces are determined by both geometric and
chemical features. However, existing works only consider handcrafted or
individually learned chemical features from the atom type and extract geometric
features independently. Here, we identify two key properties of effective
protein surface learning: 1) relationship among atoms: atoms are linked with
each other by covalent bonds to form biomolecules instead of appearing alone,
leading to the significance of modeling the relationship among atoms in
chemical feature learning. 2) hierarchical feature interaction: the neighboring
residue effect validates the significance of hierarchical feature interaction
among atoms and between surface points and atoms (or residues). In this paper,
we present a principled framework based on deep learning techniques, namely
Hierarchical Chemical and Geometric Feature Interaction Network (HCGNet), for
protein surface analysis by bridging chemical and geometric features with
hierarchical interactions. Extensive experiments demonstrate that our method
outperforms the prior state-of-the-art method by 2.3% in site prediction task
and 3.2% in interaction matching task, respectively. Our code is available at
https://github.com/xmed-lab/HCGNet.


  Developing equivariant neural networks for the E(3) group plays an important
role in modeling 3D data across real-world applications. Enforcing this
equivariance primarily involves the tensor products of irreducible
representations (irreps). However, the computational complexity of such
operations increases significantly as higher-order tensors are used. In this
work, we propose a systematic approach to substantially accelerate the
computation of the tensor products of irreps. We mathematically connect the
commonly used Clebsch-Gordan coefficients to the Gaunt coefficients, which are
integrals of products of three spherical harmonics. Through Gaunt coefficients,
the tensor product of irreps becomes equivalent to the multiplication between
spherical functions represented by spherical harmonics. This perspective
further allows us to change the basis for the equivariant operations from
spherical harmonics to a 2D Fourier basis. Consequently, the multiplication
between spherical functions represented by a 2D Fourier basis can be
efficiently computed via the convolution theorem and Fast Fourier Transforms.
This transformation reduces the complexity of full tensor products of irreps
from $\mathcal{O}(L^6)$ to $\mathcal{O}(L^3)$, where $L$ is the max degree of
irreps. Leveraging this approach, we introduce the Gaunt Tensor Product, which
serves as a new method to construct efficient equivariant operations across
different model architectures. Our experiments on the Open Catalyst Project and
3BPA datasets demonstrate both the increased efficiency and improved
performance of our approach.


  Boolean networks have been widely used to explore aspects of gene regulation,
traditionally with a single network. A modified form of the model to explore
the effects of increasing the number of gene states has also recently been
introduced. In this paper, these discrete dynamical networks are evolved as
diploids within rugged fitness landscapes to explore their behaviour. Results
suggest the general properties of haploid networks in similar circumstances
remain for diploids. The previously proposed inherent fitness landscape
smoothing properties of eukaryotic sex are shown to be exhibited in these
dynamical systems, as is their propensity to change in size based upon the
characteristics of the network and fitness landscape.


  Ecological connectivity is critical for population dynamics but in many
benthic species it is complicated by a planktonic larval phase, whose dispersal
remains poorly understood. Using a plankton pump, we examine the distribution
of intertidal mussel larvae along three axes: alongshore, cross-shelf and by
depth during a large scale (600 km) cruise over the Agulhas Bank off southern
Africa in August/September 2010. As a general pattern, higher veliger
abundances were found close to the coast. Our analyses of the nearshore flow,
estimated from ADCP data and the vertical distribution of larvae, show that
onshore larval retention may be mediated by active vertical swimming through
the water column guided by light and wind-induced turbulence. A massive
offshore export of larvae off St Francis Bay was, however, observed during an
Agulhas Current meander which influenced inner shelf waters. We hypothesize
that, by increasing and homogenizing flow, the Agulhas Current may erase the
effects of larval vertical positioning on onshore retention and transport
larvae offshore. Our study highlights the need to integrate the effects of
complex, region-specific physical dynamics with the swimming behaviour of
larvae in order to explain their spatial distribution, population connectivity
and the consequences for population dynamics.


  Dengue is a vector-borne disease transmitted to humans by vectors of genus
Aedes and is a global threat with health, social, and economic impact in many
of the tropical countries including Sri Lanka. The virus transmission is
significantly impacted by environmental conditions, with a notable contribution
from elevated per-capita vector density. These conditions are dynamic in nature
and specially having the tropical climate, Sri Lanka experiences seasonal
weather patterns dominated by monsoons. In this work, we investigate the
dynamic influence of environmental conditions on dengue emergence in Colombo
district where dengue is extremely prevalent in Sri Lanka. A novel approach
leveraging the Markov chain Monte Carlo simulations has been employed to
identify seasonal patterns of dengue disease emergence, utilizing the dynamics
of weather patterns governing in the region. The newly developed algorithm
allows us to estimate the timing of dengue outbreaks with uncertainty, enabling
accurate forecasts of upcoming disease emergence patterns for better
preparedness.


  In living cells, intrinsically disordered proteins (IDPs), such as FUS and
DDX4, undergo phase separation, forming biomolecular condensates. Using
molecular dynamics simulations, we investigate their behavior in their
respective homogenous droplets. We find that the proteins exhibit transient
subdiffusion due to the viscoelastic nature and confinement effects in the
droplets. The conformation and the instantaneous diffusivity of the proteins
significantly vary between the interior and the interface of the droplet,
resulting in non-Gaussianity in the displacement distributions. This study
highlights key aspects of IDP behavior in biomolecular condensates.


  Ecosystem models are often used to predict the consequences of management
decisions in applied ecology, including fisheries management and threatened
species conservation. These models are high-dimensional, parameter-rich, and
nonlinear, yet limited data is available to calibrate them, and they are rarely
tested or validated. Consequently, the accuracy of their forecasts, and their
utility as decision-support tools is a matter of debate. In this paper, we
calibrate ecosystem models to time-series data from 110 different experimental
microcosm ecosystems, each containing between three and five interacting
species. We then assess how often these calibrated models offer accurate and
useful predictions about how the ecosystem will respond to a set of standard
management interventions. Our results show that for each timeseries dataset, a
large number of very different parameter sets offer equivalent, good fits.
However, these calibrated ecosystem models have poor predictive accuracy when
forecasting future dynamics and offer ambiguous predictions about how species
in the ecosystem will respond to management interventions. Closer inspection
reveals that the ecosystem models fail because calibration cannot determine the
types of interactions that occur within the ecosystem. Our findings call into
question claims that ecosystem modelling can support applied ecological
decision-making when they are calibrated against real-world datasets.


  In this paper, we present dSASA (differentiable SASA), an exact geometric
method to calculate solvent accessible surface area (SASA) analytically along
with atomic derivatives on GPUs. The atoms in a molecule are first assigned to
tetrahedra in groups of four atoms by Delaunay tetrahedrization adapted for
efficient GPU implementation and the SASA values for atoms and molecules are
calculated based on the tetrahedrization information and inclusion-exclusion
method. The SASA values from the numerical icosahedral-based method can be
reproduced with more than 98% accuracy for both proteins and RNAs. Having been
implemented on GPUs and incorporated into the software Amber, we can apply
dSASA to implicit solvent molecular dynamics simulations with inclusion of this
nonpolar term. The current GPU version of GB/SA simulations has been
accelerated up to nearly 20-fold compared to the CPU version and it outperforms
LCPO as the system size increases. The performance and importance of the
nonpolar part in implicit solvent modeling are demonstrated in GB/SA
simulations of proteins and accurate SASA calculation of nucleic acids.


  With the steady rise of the use of AI in bio-technical applications and the
widespread adoption of genomics sequencing, an increasing amount of AI-based
algorithms and tools is entering the research and production stage affecting
critical decision-making streams like drug discovery and clinical outcomes.
This paper demonstrates the vulnerability of AI models often utilized
downstream tasks on recognized public genomics datasets. We undermine model
robustness by deploying an attack that focuses on input transformation while
mimicking the real data and confusing the model decision-making, ultimately
yielding a pronounced deterioration in model performance. Further, we enhance
our approach by generating poisoned data using a variational autoencoder-based
model. Our empirical findings unequivocally demonstrate a decline in model
performance, underscored by diminished accuracy and an upswing in false
positives and false negatives. Furthermore, we analyze the resulting
adversarial samples via spectral analysis yielding conclusions for
countermeasures against such attacks.


  Protein (receptor)--ligand interaction prediction is a critical component in
computer-aided drug design, significantly influencing molecular docking and
virtual screening processes. Despite the development of numerous scoring
functions in recent years, particularly those employing machine learning,
accurately and efficiently predicting binding affinities for protein--ligand
complexes remains a formidable challenge. Most contemporary methods are
tailored for specific tasks, such as binding affinity prediction, binding pose
prediction, or virtual screening, often failing to encompass all aspects. In
this study, we put forward DeepRLI, a novel protein--ligand interaction
prediction architecture. It encodes each protein--ligand complex into a fully
connected graph, retaining the integrity of the topological and spatial
structure, and leverages the improved graph transformer layers with cosine
envelope as the central module of the neural network, thus exhibiting superior
scoring power. In order to equip the model to generalize to conformations
beyond the confines of crystal structures and to adapt to molecular docking and
virtual screening tasks, we propose a multi-objective strategy, that is, the
model outputs three scores for scoring and ranking, docking, and screening, and
the training process optimizes these three objectives simultaneously. For the
latter two objectives, we augment the dataset through a docking procedure,
incorporate suitable physics-informed blocks and employ an effective
contrastive learning approach. Eventually, our model manifests a balanced
performance across scoring, ranking, docking, and screening, thereby
demonstrating its ability to handle a range of tasks. Overall, this research
contributes a multi-objective framework for universal protein--ligand
interaction prediction, augmenting the landscape of structure-based drug
design.


  Recent strides in the field of neural computation has seen the adoption of
Winner Take All (WTA) circuits to facilitate the unification of hierarchical
Bayesian inference and spiking neural networks as a neurobiologically plausible
model of information processing. Current research commonly validates the
performance of these networks via classification tasks, particularly of the
MNIST dataset. However, researchers have not yet reached consensus about how
best to translate the stochastic responses from these networks into discrete
decisions, a process known as population decoding. Despite being an often
underexamined part of SNNs, in this work we show that population decoding has a
significanct impact on the classification performance of WTA networks. For this
purpose, we apply a WTA network to the problem of cancer subtype diagnosis from
multi omic data, using datasets from The Cancer Genome Atlas (TCGA). In doing
so we utilise a novel implementation of gene similarity networks, a feature
encoding technique based on Kohoens self organising map algorithm. We further
show that the impact of selecting certain population decoding methods is
amplified when facing imbalanced datasets.


  New technologies for acquiring biological information such as eDNA, acoustic
or optical sensors, make it possible to generate spatial community observations
at unprecedented scales. The potential of these novel community data to
standardize community observations at high spatial, temporal, and taxonomic
resolution and at large spatial scale ('many rows and many columns') has been
widely discussed, but so far, there has been little integration of these data
with ecological models and theory. Here, we review these developments and
highlight emerging solutions, focusing on statistical methods for analyzing
novel community data, in particular joint species distribution models; the new
ecological questions that can be answered with these data; and the potential
implications of these developments for policy and conservation.


  In order to efficiently explore the chemical space of all possible small
molecules, a common approach is to compress the dimension of the system to
facilitate downstream machine learning tasks. Towards this end, we present a
data driven approach for clustering potential energy landscapes of molecular
structures by applying recently developed Network Embedding techniques, to
obtain latent variables defined through the embedding function. To scale up the
method, we also incorporate an entropy sensitive adaptive scheme for
hierarchical sampling of the energy landscape, based on Metadynamics and
Transition Path Theory. By taking into account the kinetic information implied
by a system's energy landscape, we are able to interpret dynamical node-node
relationships in reduced dimensions. We demonstrate the framework through
Lennard-Jones (LJ) clusters and a human DNA sequence.


  While the importance of extreme conditions is recognised, patterns in species
abundances are often interpreted through average environmental conditions
within their distributional range. For marine species with pelagic larvae,
temperature and phytoplankton concentration are key variables. Along the south
coast of South Africa, conspicuous spatial patterns in recruitment rates and
the abundances of different mussel species exist, with focal areas
characterized by large populations. We studied 15 years of sea surface
temperature (SST) and chlorophyll-a (chl-a) satellite data, using spectral
analyses to partition their temporal variability over ecologically relevant
time periods, including seasonal (101 to 365 days) and intra-seasonal cycles
(20 to 100 days). Adult cover and mussel recruitment were measured at 10 sites
along the south coast and regression models showed that about 70 percent of the
variability in recruitment and adult cover was explained by seasonal
variability in chl-a, while mean annual chl-a and SST only explained 30 percent
of the recruitment, with no significant effect for adult cover. SST and chl-a
at two upwelling centres showed less predictable seasonal cycles during the
second half of the study period with a significant cooling trend during austral
autumn, coinciding with one of the mussel reproductive peaks. This likely
reflects recent changes in the Agulhas Current, the world largest western
boundary current, which affects coastal ecosystems by driving upwelling.


  Recent advances in protein language models have catalyzed significant
progress in peptide sequence representation. Despite extensive exploration in
this field, pre-trained models tailored for peptide-specific needs remain
largely unaddressed due to the difficulty in capturing the complex and
sometimes unstable structures of peptides. This study introduces a novel
multi-view contrastive learning framework PepHarmony for the sequence-based
peptide encoding task. PepHarmony innovatively combines both sequence- and
structure-level information into a sequence-level encoding module through
contrastive learning. We carefully select datasets from the Protein Data Bank
(PDB) and AlphaFold database to encompass a broad spectrum of peptide sequences
and structures. The experimental data highlights PepHarmony's exceptional
capability in capturing the intricate relationship between peptide sequences
and structures compared with the baseline and fine-tuned models. The robustness
of our model is confirmed through extensive ablation studies, which emphasize
the crucial roles of contrastive loss and strategic data sorting in enhancing
predictive performance. The proposed PepHarmony framework serves as a notable
contribution to peptide representations, and offers valuable insights for
future applications in peptide drug discovery and peptide engineering. We have
made all the source code utilized in this study publicly accessible via GitHub
at https://github.com/zhangruochi/PepHarmony or
http://www.healthinformaticslab.org/supp/.


  Deep learning is now widely used in drug discovery, providing significant
acceleration and cost reduction. As the most fundamental building block,
molecular representation is essential for predicting molecular properties to
enable various downstream applications. Most existing methods attempt to
incorporate more information to learn better representations. However, not all
features are equally important for a specific task. Ignoring this would
potentially compromise the training efficiency and predictive accuracy. To
address this issue, we propose a novel approach, which treats language models
as an agent and molecular pretraining models as a knowledge base. The agent
accentuates task-relevant features in the molecular representation by
understanding the natural language description of the task, just as a tailor
customizes clothes for clients. Thus, we call this approach MolTailor.
Evaluations demonstrate MolTailor's superior performance over baselines,
validating the efficacy of enhancing relevance for molecular representation
learning. This illustrates the potential of language model guided optimization
to better exploit and unleash the capabilities of existing powerful molecular
representation methods. Our codes and appendix are available at
https://github.com/SCIR-HI/MolTailor.


  Proliferation of uninfected as well as infected hepatocytes and recycling of
DNA-containing
  capsids are two major mechanisms playing significant roles in the clearance
of hepatitis B
  virus (HBV) infection. In this study, the temporal dynamics of this infection
are investigated
  through two in silico bio-mathematical models considering both proliferation
of hepatocytes
  and the recycling of capsids. Both models are formulated on the basis of a
key finding in the existing literature: mitosis of infected yields in two
uninfected progenies. In the first model,
  we examine regular proliferation (occurs continuously), while the second
model deals with the
  irregular proliferation (happens when the total number of liver cells
decreases to less than 70%
  of its initial volume). The models are calibrated with the experimental data
obtained from
  an adult chimpanzee. Results of this study suggest that when both hepatocytes
proliferate
  with equal rate, proliferation aids the individual in a rapid recovery from
the acute infection
  whereas in the case of chronic infection, the severity of the infection
increases if the proliferation
  occurs frequently. On the other hand, if the infected cells proliferate at a
slower rate than uninfected cells, the proliferation of uninfected hepatocytes
contributes to increase the infection,
  but the proliferation of infected hepatocytes acts to reduce the infection
from the long-term
  perspective. Furthermore, it is also observed that the differences between
the outcomes of
  regular and irregular proliferations are substantial and noteworthy.


  Although defocus can be used to generate partial phase contrast in
transmission electron microscope images, cryo-electron microscopy (cryo-EM) can
be further improved by the development of phase plates which increase contrast
by applying a phase shift to the unscattered part of the electron beam. Many
approaches have been investigated, including the ponderomotive interaction
between light and electrons. We review the recent successes achieved with this
method in high-resolution, single-particle cryo-EM. We also review the status
of using pulsed or near-field enhanced laser light as alternatives, along with
approaches that use scanning transmission electron microscopy (STEM) with a
segmented detector rather than a phase plate.


  DNA sequence alignment is an important workload in computational genomics.
Reference-guided DNA assembly involves aligning many read sequences against
candidate locations in a long reference genome. To reduce the computational
load of this alignment, candidate locations can be pre-filtered using simpler
alignment algorithms like edit distance. Prior work has explored accelerating
filtering on simulated compute-in-DRAM, due to the massive parallelism of
compute-in-memory architectures. In this paper, we present work-in-progress on
accelerating filtering using a commercial compute-in-SRAM accelerator. We
leverage the recently released Gemini accelerator platform from GSI Technology,
which is the first, to our knowledge, commercial-scale compute-in-SRAM system.
We accelerate the Myers' bit-parallel edit distance algorithm, producing
average speedups of 14.1x over single-core CPU performance. Individual
query/candidate alignments produce speedups of up to 24.1x. These early results
suggest this novel architecture is well-suited to accelerating the filtering
step of sequence-to-sequence DNA alignment.


  Multiplayer games on graphs are at the heart of theoretical descriptions of
key evolutionary processes that govern vital social and natural systems.
However, a comprehensive theoretical framework for solving multiplayer games
with an arbitrary number of strategies on graphs is still missing. Here, we
solve this by drawing an analogy with the Ball-and-Box problem, based on which
we show that the local configuration of multiplayer games on graphs is
equivalent to distributing $k$ identical co-players among $n$ distinct
strategies. We use this to derive the replicator equation for any $n$-strategy
multiplayer game under weak selection, which can be solved in polynomial time.
As an example, we revisit the second-order free-riding problem, where costly
punishment cannot truly resolve social dilemmas in a well-mixed population.
Yet, in structured populations, we derive an accurate threshold for the
punishment strength, beyond which punishment can either lead to the extinction
of defection or transform the system into a rock-paper-scissors-like cycle. The
analytical solution also qualitatively agrees with the phase diagrams that were
previously obtained for non-marginal selection strengths. Our framework thus
allows an exploration of any multi-strategy multiplayer game on regular graphs.


  Oligo pools are array-synthesized, user-defined mixtures of single-stranded
oligonucleotides that can be used as a source of synthetic DNA for library
cloning. While currently offering the most affordable source of synthetic DNA,
oligo pools also come with limitations such as a maximum synthesis length
(approximately 350 bases), a higher error rate compared to alternative
synthesis methods, and the presence of truncated molecules in the pool due to
incomplete synthesis. Here, we provide users with a comprehensive protocol that
details how oligo pools can be used in combination with Golden Gate cloning to
create user-defined protein mutant libraries, as well as single guide RNA
libraries for CRISPR applications. Our methods are optimized to work within the
Yeast Toolkit Golden Gate scheme, but are in principle compatible with any
other Golden Gate-based modular cloning toolkit and extendable to other
restriction enzyme-based cloning methods beyond Golden Gate. Our methods yield
high-quality, affordable, in-house variant libraries.


  This research investigates the impact of dynamic interactions with
time-varying topologies on the evolution of cooperative behaviours in social
dilemmas. Traditional research has focused on deterministic rules governing
pairwise interactions, yet the impact of interaction frequency and
synchronicity on cooperation remains underexplored. Addressing this gap, our
work introduces two temporal interaction mechanisms to model the stochastic or
periodic participation of individuals in these games, acknowledging real-life
variances due to exogenous temporal factors and geographical time differences.
We consider that the interaction state significantly influences both game
payoff calculations and the strategy updating process, offering new insights
into the emergence and sustainability of cooperation. Our results indicate that
maximum game participation frequency is suboptimal under a stochastic
interaction mechanism. Instead, an intermediate region of activation
probability yields the highest cooperation level, especially under strong
dilemma conditions. This suggests that a balance between inactivity security
and interaction frequency is crucial. Furthermore, local synchronization of
interactions within specific areas is shown to be beneficial, as time
differences hinder the spread of cross-structures but promote the formation of
dense cooperative clusters with smoother boundaries. Our findings provide an
intuitive understanding of node-based temporality and probabilistic
interactions, contributing to the broader discourse on resolving social
dilemmas.


  Given non-sequential snapshots from instances of a dynamical system, we
design a compressed sensing based algorithm that reconstructs the dynamical
system. We formally prove that successful reconstruction is possible under the
assumption that we can construct an approximate clock from a subset of the
coordinates of the underlying system.
  As an application, we argue that our assumption is likely true for genomic
datasets, and we recover the underlying nuclear receptor networks and predict
pathways, as opposed to genes, that may differentiate phenotypes in some
publicly available datasets.


  Motivation: Assessing the match between two biomolecular structures is at the
heart of structural analyses such as superposition, alignment and docking.
These tasks are typically solved with specialized structure-matching techniques
implemented in software for protein structural alignment, rigid-body docking,
or rigid fitting into cryo-EM maps. Results: We present a unifying framework to
compare biomolecular structures by applying ideas from computer vision. The
structures are represented as three-dimensional point clouds and compared by
quantifying their overlap. We use the kernel correlation to measure point cloud
overlap, and discuss local and global optimization strategies for maximizing
the kernel correlation over the space of rigid transformations. We derive a
majorization-minimization procedure that can be used to register two point
clouds without establishing a point-to-point correspondence. We demonstrate
that the majorization-minimization algorithms outperform the commonly used
Iterative Closest Point registration algorithm. Furthermore, we discuss and
benchmark a randomization strategy for globally optimizing the kernel
correlation. We illustrate the approach on various 3D fitting problems such as
the comparison of circularly permuted structures and rigid fitting of cryo-EM
maps or bead models from small-angle scattering.


  The delimitation of biological species, i.e., deciding which individuals
belong to the same species and whether and how many different species are
represented in a data set, is key to the conservation of biodiversity. Much
existing work uses only genetic data for species delimitation, often employing
some kind of cluster analysis. This can be misleading, because geographically
distant groups of individuals can be genetically quite different even if they
belong to the same species. This paper investigates the problem of testing
whether two potentially separated groups of individuals can belong to a single
species or not based on genetic and spatial data. Various approaches are
compared (some of which already exist in the literature) based on simulated
metapopulations generated with SLiM and GSpace - two software packages that can
simulate spatially-explicit genetic data at an individual level. Approaches
involve partial Mantel testing, maximum likelihood mixed-effects models with a
population effect, and jackknife-based homogeneity tests. A key challenge is
that most tests perform on genetic and geographical distance data, violating
standard independence assumptions. Simulations showed that partial Mantel tests
and mixed-effects models have larger power than jackknife-based methods, but
tend to display type-I-error rates slightly above the significance level.
Moreover, a multiple regression model neglecting the dependence in the
dissimilarities did not show inflated type-I-error rate. An application on
brassy ringlets concludes the paper.


  The pharmaceutical Research and development (R&D) process is lengthy and
costly, taking nearly a decade to bring a new drug to the market. However,
advancements in biotechnology, computational methods, and machine learning
algorithms have the potential to revolutionize drug discovery, speeding up the
process and improving patient outcomes. The COVID-19 pandemic has further
accelerated and deepened the recognition of the potential of these techniques,
especially in the areas of drug repurposing and efficacy predictions.
Meanwhile, non-small molecule therapeutic modalities such as cell therapies,
monoclonal antibodies, and RNA interference (RNAi) technology have gained
importance due to their ability to target specific disease pathways and/or
patient populations. In the field of RNAi, many experiments have been carried
out to design and select highly efficient siRNAs. However, the established
patterns for efficient siRNAs are sometimes contradictory and unable to
consistently determine the most potent siRNA molecules against a target mRNA.
Thus, this paper focuses on developing machine learning models based on the
cheminformatics representation of the nucleotide composition (i.e. AUTGC) of
siRNA to predict their potency and aid the selection of the most efficient
siRNAs for further development. The PLS (Partial Least Square) and SVR (Support
Vector Regression) machine learning models built in this work outperformed
previously published models. These models can help in predicting siRNA potency
and aid in selecting the best siRNA molecules for experimental validation and
further clinical development. The study has demonstrated the potential of
AI/machine learning models to help expedite siRNA-based drug discovery
including the discovery of potent siRNAs against SARS-CoV-2.


  While quite a few recent papers have explored game-resource feedback using
the framework of evolutionary game theory, almost all the studies are confined
to using time-continuous dynamical equations. Moreover, in such literature, the
effect of ubiquitous chaos in the resulting eco-evolutionary dynamics is rather
missing. Here, we present a deterministic eco-evolutionary discrete-time
dynamics in generation-wise non-overlapping population of two types of
harvesters, one harvesting at a faster rate than the other, consuming a
self-renewing resource capable of showing chaotic dynamics. In the light of our
finding that sometimes chaos is confined exclusively to either the dynamics of
the resource or that of the consumer fractions, an interesting scenario is
realized: The resource state can keep oscillating chaotically, and hence, it
does not vanish to result in the tragedy of the commons, extinction of the
resource due to selfish indiscriminate exploitation, and yet the consumer
population, whose dynamics depends directly on the state of the resource, may
end up being composed exclusively of defectors, i.e., high harvesters. This
appears non-intuitive because it is well known that prevention of tragedy of
the commons usually requires substantial cooperation to be present.


  China is one of the countries with a high burden of tuberculosis, and
although the number of new cases of tuberculosis has been decreasing year by
year, the number of new infections per year has remained high and the diagnosis
rate of tuberculosis-infected patients has remained low. Based on the analysis
of TB infection data, we develop a model of TB transmission dynamics that
include potentially infected individuals and BCG vaccination, fit the model
parameters to the data on new TB cases, calculate the basic reproduction number
\mathcal{R}_v= 0.4442. A parametric sensitivity analysis of \mathcal{R}_v is
performed, and we obtained the correlation coefficients of BCG vaccination rate
and effectiveness rate with \mathcal{R}_v as -0.810, -0.825. According to the
model, we estimate that there are 614,186 (95% CI [562,631,665,741])
potentially infected TB cases in China, accounting for about 39.5% of the total
number of TB cases. We assess the feasibility of achieving the goals of the WHO
strategy to end tuberculosis in China and find that reducing the number of new
cases by 90 per cent by 2035 is very difficult with the current tuberculosis
control measures. However, with an effective combination of control measures
such as increased detection of potentially infected persons, improved drug
treatment, and reduction of overall exposure to tuberculosis patients, it is
feasible to reach the WHO strategic goal of ending tuberculosis by 2035.


  The concept of control is central to understanding and applications of
biological network models. Some of their key structural features relate to
control functions, through gene regulation, signaling, or metabolic mechanisms,
and computational models need to encode these. Applications of models often
focus on model-based control, such as in biomedicine or metabolic engineering.
This paper presents an approach to model-based control that exploits two common
features of biological networks, namely their modular structure and canalizing
features of their regulatory mechanisms. The paper focuses on intracellular
regulatory networks, represented by Boolean network models. A main result of
this paper is that control strategies can be identified by focusing on one
module at a time. This paper also presents a criterion based on canalizing
features of the regulatory rules to identify modules that do not contribute to
network control and can be excluded. For even moderately sized networks,
finding global control inputs is computationally very challenging. The modular
approach presented here leads to a highly efficient approach to solving this
problem. This approach is applied to a published Boolean network model of blood
cancer large granular lymphocyte (T-LGL) leukemia to identify a minimal control
set that achieves a desired control objective.


  Noise plays a crucial role in the regulation of cellular and organismal
function and behavior.
  Exploring noise's impact is key to understanding fundamental biological
processes, such as gene expression, signal transduction, and the mechanisms of
development and evolution.
  Currently, a comprehensive method to quantify dynamical behavior of cellular
noise within these biochemical systems is lacking.
  In this study, we introduce an optically-controlled perturbation system
utilizing the light-sensitive Phytochrome B (PhyB) from \textit{Arabidopsis
thaliana}, which enables precise noise modulation with high spatial-temporal
resolution.
  Our system exhibits exceptional sensitivity to light, reacting consistently
to pulsed light signals, distinguishing it from other photoreceptor-based
promoter systems that respond to a single light wavelength.
  To characterize our system, we developed a stochastic model for phytochromes
that accounts for photoactivation/deactivation, thermal reversion, and the
dynamics of the light-activated gene promoter system.
  To precisely control our system, we determined the rate constants for this
model using an omniscient deep neural network that can directly map rate
constant combinations to time-dependent state joint distributions.
  By adjusting the activation rates through light intensity and degradation
rates via N-terminal mutagenesis, we illustrate that out optical-controlled
perturbation can effectively modulate molecular expression level as well as
noise.
  Our results highlight the potential of employing an optically-controlled gene
perturbation system as a noise-controlled stimulus source.
  This approach, when combined with the analytical capabilities of a
sophisticated deep neural network, enables the accurate estimation of rate
constants from observational data in a broad range of biochemical reaction
networks.


  We introduce a biologically detailed, stochastic model of gene expression
describing the multiple rate-limiting steps of transcription, nuclear pre-mRNA
processing, nuclear mRNA export, cytoplasmic mRNA degradation and translation
of mRNA into protein. The processes in sub-cellular compartments are described
by an arbitrary number of processing stages, thus accounting for a
significantly finer molecular description of gene expression than conventional
models such as the telegraph, two-stage and three-stage models of gene
expression. We use two distinct tools, queueing theory and model reduction
using the slow-scale linear-noise approximation, to derive exact or approximate
analytic expressions for the moments or distributions of nuclear mRNA,
cytoplasmic mRNA and protein fluctuations, as well as lower bounds for their
Fano factors in steady-state conditions. We use these to study the phase
diagram of the stochastic model; in particular we derive parametric conditions
determining three types of transitions in the properties of mRNA fluctuations:
from sub-Poissonian to super-Poissonian noise, from high noise in the nucleus
to high noise in the cytoplasm, and from a monotonic increase to a monotonic
decrease of the Fano factor with the number of processing stages. In contrast,
protein fluctuations are always super-Poissonian and show weak dependence on
the number of mRNA processing stages. Our results delineate the region of
parameter space where conventional models give qualitatively incorrect results
and provide insight into how the number of processing stages, e.g. the number
of rate-limiting steps in initiation, splicing and mRNA degradation, shape
stochastic gene expression by modulation of molecular memory.


  How cooperation evolves and particularly maintains at a large scale remains
an open problem for improving humanity across domains ranging from climate
change to pandemic response. To shed light on how behavioral norms can resolve
the social dilemma of cooperation, here we present a formal mathematical model
of individuals' decision making under general social norms, encompassing a
variety of concerns and motivations an individual may have beyond simply
maximizing their own payoffs. Using the canonical game of the Prisoner's
Dilemma, we compare four different norms: compassion, universalizability,
reciprocity, and equity, to determine which social forces can facilitate the
evolution of cooperation, if any. We analyze our model through a variety of
limiting cases, including weak selection, low mutation, and large population
sizes. This is complemented by computer simulations of population dynamics via
a Fisher process, which confirm our theoretical results. We find that the first
two norms lead to the emergence of cooperation in a wide range of games, but
the latter two do not on their own. Due to its generality, our framework can be
used to investigate many more norms, as well as how norms themselves emerge and
evolve. Our work complements recent work on fair-minded learning dynamics and
provides a useful bottom-up perspective into understanding the impact of
top-down social norms on collective cooperative intelligence.


  A species' genetic code or genome encodes valuable evolutionary, biological,
and phylogenetic information that aids in species recognition, taxonomic
classification, and understanding genetic predispositions like drug resistance
and virulence. However, the vast number of potential species poses significant
challenges in developing a general-purpose whole genome classification tool.
Traditional bioinformatics tools have made notable progress but lack
scalability and are computationally expensive. Machine learning-based
frameworks show promise but must address the issue of large classification
vocabularies with long-tail distributions. In this study, we propose addressing
this problem through zero-shot learning using TEPI, Taxonomy-aware Embedding
and Pseudo-Imaging. We represent each genome as pseudo-images and map them to a
taxonomy-aware embedding space for reasoning and classification. This embedding
space captures compositional and phylogenetic relationships of species,
enabling predictions in extensive search spaces. We evaluate TEPI using two
rigorous zero-shot settings and demonstrate its generalization capabilities
qualitatively on curated, large-scale, publicly sourced data.


  Although continuous advances in theoretical modelling of Molecular
Communications (MC) are observed, there is still an insuperable gap between
theory and experimental testbeds, especially at the microscale. In this paper,
the development of the first testbed incorporating engineered yeast cells is
reported. Different from the existing literature, eukaryotic yeast cells are
considered for both the sender and the receiver, with {\alpha}-factor molecules
facilitating the information transfer. The use of such cells is motivated
mainly by the well understood biological mechanism of yeast mating, together
with their genetic amenability. In addition, recent advances in yeast
biosensing establish yeast as a suitable detector and a neat interface to
in-body sensor networks. The system under consideration is presented first, and
the mathematical models of the underlying biological processes leading to an
end-to-end (E2E) system are given. The experimental setup is then described and
used to obtain experimental results which validate the developed mathematical
models. Beyond that, the ability of the system to effectively generate output
pulses in response to repeated stimuli is demonstrated, reporting one event per
two hours. However, fast RNA fluctuations indicate cell responses in less than
three minutes, demonstrating the potential for much higher rates in the future.


  Inverse molecular design with diffusion models holds great potential for
advancements in material and drug discovery. Despite success in unconditional
molecule generation, integrating multiple properties such as synthetic score
and gas permeability as condition constraints into diffusion models remains
unexplored. We introduce multi-conditional diffusion guidance. The proposed
Transformer-based denoising model has a condition encoder that learns the
representations of numerical and categorical conditions. The denoising model,
consisting of a structure encoder-decoder, is trained for denoising under the
representation of conditions. The diffusion process becomes graph-dependent to
accurately estimate graph-related noise in molecules, unlike the previous
models that focus solely on the marginal distributions of atoms or bonds. We
extensively validate our model for multi-conditional polymer and small molecule
generation. Results demonstrate our superiority across metrics from
distribution learning to condition control for molecular properties. An inverse
polymer design task for gas separation with feedback from domain experts
further demonstrates its practical utility.


  Language Models (LMs) have greatly influenced diverse domains. However, their
inherent limitation in comprehending 3D molecular structures has considerably
constrained their potential in the biomolecular domain. To bridge this gap, we
focus on 3D molecule-text interpretation, and propose 3D-MoLM: 3D-Molecular
Language Modeling. Specifically, 3D-MoLM enables an LM to interpret and analyze
3D molecules by equipping the LM with a 3D molecular encoder. This integration
is achieved by a 3D molecule-text projector, bridging the 3D molecular
encoder's representation space and the LM's input space. Moreover, to enhance
3D-MoLM's ability of cross-modal molecular understanding and instruction
following, we meticulously curated a 3D molecule-centric instruction tuning
dataset -- 3D-MoIT. Through 3D molecule-text alignment and 3D molecule-centric
instruction tuning, 3D-MoLM establishes an integration of 3D molecular encoder
and LM. It significantly surpasses existing baselines on downstream tasks,
including molecule-text retrieval, molecule captioning, and more challenging
open-text molecular QA tasks, especially focusing on 3D-dependent properties.


  Nearly all circadian clocks maintain a period that is insensitive to
temperature changes, a phenomenon known as temperature compensation (TC). Yet,
it is unclear whether there is any common feature among different systems that
exhibit TC. From a general timescale invariance, we show that TC relies on
existence of certain period-lengthening reactions wherein the period of the
system increases strongly with the rates in these reactions. By studying
several generic oscillator models, we show that this counter-intuitive
dependence is nonetheless a common feature of oscillators in the nonlinear
(far-from-onset) regime where the oscillation can be separated into fast and
slow phases. The increase of the period with the period-lengthening reaction
rates occurs when the amplitude of the slow phase in the oscillation increases
with these rates while the progression-speed in the slow phase is controlled by
other rates of the system. The positive dependence of the period on the
period-lengthening rates balances its inverse dependence on other kinetic rates
in the system, which gives rise to robust TC in a wide range of parameters. We
demonstrate the existence of such period-lengthening reactions and their
relevance for TC in all four model systems we considered. Theoretical results
for a model of the Kai system are supported by experimental data. A study of
the energy dissipation also shows that better TC performance requires higher
energy consumption. Our study unveils a general mechanism by which a
biochemical oscillator achieves TC by operating at regimes far from the onset
where period-lengthening reactions exist.


  Recent studies in DNA sequence classification have leveraged sophisticated
machine learning techniques, achieving notable accuracy in categorizing complex
genomic data. Among these, methods such as k-mer counting have proven effective
in distinguishing sequences from varied species like chimpanzees, dogs, and
humans, becoming a staple in contemporary genomic research. However, these
approaches often demand extensive computational resources, posing a challenge
in terms of scalability and efficiency. Addressing this issue, our study
introduces a novel adaptation of Jiang et al.'s compressor-based,
parameter-free classification method, specifically tailored for DNA sequence
analysis. This innovative approach utilizes a variety of compression
algorithms, such as Gzip, Brotli, and LZMA, to efficiently process and classify
genomic sequences. Not only does this method align with the current
state-of-the-art in terms of accuracy, but it also offers a more
resource-efficient alternative to traditional machine learning methods. Our
comprehensive evaluation demonstrates the proposed method's effectiveness in
accurately classifying DNA sequences from multiple species. We present a
detailed analysis of the performance of each algorithm used, highlighting the
strengths and limitations of our approach in various genomic contexts.
Furthermore, we discuss the broader implications of our findings for
bioinformatics, particularly in genomic data processing and analysis. The
results of our study pave the way for more efficient and scalable DNA sequence
classification methods, offering significant potential for advancements in
genomic research and applications.


  Exosomes are significant facilitators of inter-cellular communication that
can unveil cell-cell interactions, signaling pathways, regulatory mechanisms
and disease diagnostics. Nonetheless, current analysis required large amount of
data for exosome identification that it hampers efficient and timely mechanism
study and diagnostics. Here, we used a machine-learning assisted
Surface-enhanced Raman spectroscopy (SERS) method to detect exosomes derived
from six distinct cell lines (HepG2, Hela, 143B, LO-2, BMSC, and H8) with small
amount of data. By employing sodium borohydride-reduced silver nanoparticles
and sodium borohydride solution as an aggregating agent, 100 SERS spectra of
the each types of exosomes were collected and then subjected to multivariate
and machine learning analysis. By integrating Principal Component Analysis with
Support Vector Machine (PCA-SVM) models, our analysis achieved a high accuracy
rate of 94.4% in predicting exosomes originating from various cellular sources.
In comparison to other machine learning analysis, our method used small amount
of SERS data to allow a simple and rapid exosome detection, which enables a
timely subsequent study of cell-cell interactions, communication mechanisms,
and disease mechanisms in life sciences.


  Motivated by the evolution of complex bird songs, an abstract imitation game
is proposed to study the increase of dynamical complexity: Artificial "birds"
display a "song" time series to each other, and those that imitate the other's
song better win the game. With the introduction of population dynamics
according to the score of the game and the mutation of parameters for the song
dynamics, the dynamics is found to evolve towards the borderline between chaos
and a periodic window, after punctuated equilibria. The importance of edge of
chaos with topological chaos for complexity is stressed.


  A spacially extended model of the collective behavior of a large number of
locally acting organisms is proposed in which organisms move probabilistically
between local cells in space, but with weights dependent on local morphogenetic
substances, or morphogens. The morphogens are in turn are effected by the
passage of an organism. The evolution of the morphogens, and the corresponding
flow of the organisms constitutes the collective behavior of the group. Such
models have various types of phase transitions and self-organizing properties
controlled both by the level of the noise, and other parameters.
  The model is then applied to the specific case of ants moving on a lattice.
The local behavior of the ants is inspired by the actual behavior observed in
the laboratory, and analytic results for the collective behavior are compared
to the corresponding laboratory results.
  It is hoped that the present model might serve as a paradigmatic example of a
complex cooperative system in nature. In particular swarm models can be used to
explore the relation of nonequilibrium phase transitions to at least three
important issues encountered in artificial life. Firstly, that of emergence as
complex adaptive behavior. Secondly, as an exploration of continuous phase
transitions in biological systems. Lastly, to derive behavioral criteria for
the evolution of collective behavior in social organisms.


  Using an artificial system of self-replicating strings, we show a correlation
between the age of a genotype and its abundance that reflects a punctuated
rather than gradual picture of evolution, as suggested long ago by Willis. In
support of this correlation, we measure genotype abundance distributions and
find universal coefficients. Finally, we propose a simple stochastic model
which describes the dynamics of equilibrium periods and which correctly
predicts most of the observed distributions.


  This work approaches human chromosome mapping by developing algorithms for
ordering markers associated with radiation hybrid data. Motivated by recent
work of Boehnke et al. [1], we formulate the ordering problem by developing
stochastic spin models to search for minimum-break marker configurations. As a
particular application, the methods developed are applied to 14 human
chromosome-21 markers tested by Cox et al. [2]. The methods generate
configurations consistent with the best found by others. Additionally, we find
that the set of low-lying configurations is described by a Markov-like ordering
probability distribution. The distribution displays cluster correlations
reflecting closely linked loci.


  We observe the propagation of information in a system of self-replicating
strings of code (``Artificial Life'') as a function of fitness and mutation
rate. Comparison with theoretical predictions based on the reaction-diffusion
equation shows that the response of the artificial system to fluctuations (\eg
velocity of the information wave as a function of relative fitness) closely
follows that of natural systems. We find that the relaxation time of the system
depends on the speed of propagation of information and the size of the system.
This analysis offers the possibility of determining the minimal system size for
observation of non-equilibrium effects at fixed mutation rate.


  A practical measure for the complexity of sequences of symbols (``strings'')
is introduced that is rooted in automata theory but avoids the problems of
Kolmogorov-Chaitin complexity. This physical complexity can be estimated for
ensembles of sequences, for which it reverts to the difference between the
maximal entropy of the ensemble and the actual entropy given the specific
environment within which the sequence is to be interpreted. Thus, the physical
complexity measures the amount of information about the environment that is
coded in the sequence, and is conditional on such an environment. In practice,
an estimate of the complexity of a string can be obtained by counting the
number of loci per string that are fixed in the ensemble, while the volatile
positions represent, again with respect to the environment, randomness. We
apply this measure to tRNA sequence data.


  Altruistic behaviour is disadvantageous for the individual while is
advantageous for its group. If the target of the selection is the individual,
one would expect the selection process to lead to populations formed by wholly
homogeneous groups, made up of either altruistic or egoistic individuals, where
the winning choice depends on the balance beetwen group advantage and
individual disadvantage. We show in a simple model that populations formed by
inhomogeneous groups can be stabilized in some circumstances. We argue that
this condition is realized when there is a relative advantage conferred by the
presence of a few altruists to all the members of the group.


  The sustainable biodiversity associated with a specific ecological niche as a
function of land area is analysed computationally by considering the
interaction of ant societies over a collection of islands. A power law
relationship between sustainable species and land area is observed. We will
further consider the effect a perturbative inflow of ants has upon the model.


  The study of correlation structure in the primary sequences of DNA is
reviewed. The issues reviewed include: symmetries among 16 base-base
correlation functions, accurate estimation of correlation measures, the
relationship between $1/f$ and Lorentzian spectra, heterogeneity in DNA
sequences, different modeling strategies of the correlation structure of DNA
sequences, the difference of correlation structure between coding and
non-coding regions (besides the period-3 pattern), and source of broad
distribution of domain sizes. Although some of the results remain
controversial, a body of work on this topic constitutes a good starting point
for future studies.


  Extending the classic works of Berg and Purcell on the biophysics of
bacterial chemotaxis, we find the optimal chemotactic strategy for the
peritrichous bacterium E. Coli in the high and low signal to noise ratio
limits. The optimal strategy depends on properties of the environment and
properties of the individual bacterium and is therefore highly adaptive. We
review experiments relevant to testing both the form of the proposed strategy
and its adaptability and propose extensions of them which could test the limits
of the adaptability in this simplest sensory processing system.


  We present a model for biological aging that considers the number of
individuals whose (inherited) genetic charge determines the maximum age for
death: each individual may die before that age due to some external factor, but
never after that limit. The genetic charge of the offspring is inherited from
the parent with some mutations, described by a transition matrix. The model can
describe different strategies of reproduction and it is exactly soluble. We
applied our method to the bit-string model for aging and the results are in
perfect agreement with numerical simulations.


  In this work we check the occurrence of the Azbel assumption of mortality
within the framework of a bit string model for biological ageing. We reproduced
the observed feature of linear correspondence between the fitting parameters of
the death rate as obtained by Azbel with demographic data.


  This commentary discusses a recently proposed measure of heterogeneity of DNA
sequences and compares with the measures of complexity.


  We analyze the geometry of the species- and genotype-size distribution in
evolving and adapting populations of single-stranded self-replicating genomes:
here programs in the Avida world. We find that a scale-free distribution (power
law) emerges in complex landscapes that achieve a separation of two fundamental
time scales: the relaxation time (time for population to return to equilibrium
after a perturbation) and the time between mutations that produce fitter
genotypes. The latter can be dialed by changing the mutation rate. In the
scaling regime, we determine the critical exponent of the distribution of sizes
and strengths of avalanches in a system without coevolution, described by
first-order phase transitions in single finite niches.


  We present a model of decentralized growth for Artificial Neural Networks
(ANNs) inspired by the development and the physiology of real nervous systems.
In this model, each individual artificial neuron is an autonomous unit whose
behavior is determined only by the genetic information it harbors and local
concentrations of substrates modeled by a simple artificial chemistry. Gene
expression is manifested as axon and dendrite growth, cell division and
differentiation, substrate production and cell stimulation. We demonstrate the
model's power with a hand-written genome that leads to the growth of a simple
network which performs classical conditioning. To evolve more complex
structures, we implemented a platform-independent, asynchronous, distributed
Genetic Algorithm (GA) that allows users to participate in evolutionary
experiments via the World Wide Web.


  Epochal dynamics, in which long periods of stasis in population fitness are
punctuated by sudden innovations, is a common behavior in both natural and
artificial evolutionary processes. We use a recent quantitative mathematical
analysis of epochal evolution to estimate, as a function of population size and
mutation rate, the average number of fitness function evaluations to reach the
global optimum. This is then used to derive estimates of and bounds on
evolutionary parameters that minimize search effort.


  Epochal dynamics, in which long periods of stasis in an evolving population
are punctuated by a sudden burst of change, is a common behavior in both
natural and artificial evolutionary processes. We analyze the population
dynamics for a class of fitness functions that exhibit epochal behavior using a
mathematical framework developed recently. In the latter the approximations
employed led to a population-size independent theory that allowed us to
determine optimal mutation rates. Here we extend this approach to include the
destabilization of epochs due to finite-population fluctuations and show that
this dynamical behavior often occurs around the optimal parameter settings for
efficient search. The resulting, more accurate theory predicts the total number
of fitness function evaluations to reach the global optimum as a function of
mutation rate, population size, and the parameters specifying the fitness
function. We further identify a generalized error threshold, smoothly bounding
the two-dimensional regime of mutation rates and population sizes for which
evolutionary search operates efficiently.


  Based on the law of mass action (and its microscopic foundation) and mass
conservation, we present here a method to derive consistent dynamic models for
the time evolution of systems with an arbitrary number of species. Equations
are derived through a mechanistic description, ensuring that all parameters
have ecological meaning. After discussing the biological mechanisms associated
to the logistic and Lotka-Volterra equations, we show how to derive general
models for trophic chains, including the effects of internal states at fast
time scales. We show that conformity with the mass action law leads to
different functional forms for the Lotka-Volterra and trophic chain models. We
use mass conservation to recover the concept of carrying capacity for an
arbitrary food chain.


  We analyze the population dynamics of a broad class of fitness functions that
exhibit epochal evolution---a dynamical behavior, commonly observed in both
natural and artificial evolutionary processes, in which long periods of stasis
in an evolving population are punctuated by sudden bursts of change. Our
approach---statistical dynamics---combines methods from both statistical
mechanics and dynamical systems theory in a way that offers an alternative to
current ``landscape'' models of evolutionary optimization. We describe the
population dynamics on the macroscopic level of fitness classes or phenotype
subbasins, while averaging out the genotypic variation that is consistent with
a macroscopic state. Metastability in epochal evolution occurs solely at the
macroscopic level of the fitness distribution. While a balance between
selection and mutation maintains a quasistationary distribution of fitness,
individuals diffuse randomly through selectively neutral subbasins in genotype
space. Sudden innovations occur when, through this diffusion, a genotypic
portal is discovered that connects to a new subbasin of higher fitness
genotypes. In this way, we identify innovations with the unfolding and
stabilization of a new dimension in the macroscopic state space. The
architectural view of subbasins and portals in genotype space clarifies how
frozen accidents and the resulting phenotypic constraints guide the evolution
to higher complexity.


  The aims of this paper are to propose, construct and analyse microscopic
kinetic models for the emergence of long chains of RNA from monomeric
beta-D-ribonucleotide precursors in prebiotic circumstances. Our theory starts
out from similar but more general chemical assumptions to those of Eigen,
namely that catalytic replication can lead to a large population of long
chains. In particular, our models incorporate the possibility of (i) direct
chain growth, (ii) template-assisted synthesis and (iii) catalysis by RNA
replicase ribozymes, all with varying degrees of efficiency. However, in our
models the reaction mechanisms are kept `open'; we do not assume the existence
of closed hypercycles which sustain a population of long chains. Rather it is
the feasibility of the initial emergence of a self-sustaining set of RNA chains
from monomeric nucleotides which is our prime concern. We confront directly the
central nonlinear features of the problem, which have often been overlooked in
previous studies. Our detailed microscopic kinetic models lead to kinetic
equations which are generalisations of the Becker-Doring system (BD) for the
step-wise growth of clusters or polymer chains; they lie within a general
theoretical framework which has recently been successfully applied to a wide
range of complex chemical problems. In fact, the most accurate model we
consider has BD aggregation terms, together with a general Smoluchowski
fragmentation term to model the competing hydrolysis of RNA polymer chains. We
conclude that, starting from reasonable initial conditions of monomeric
nucleotide concentrations within a prebiotic soup and in an acceptable
timescale, it is possible for a self-replicating subset of polyribonucleotide
chains to be selected, while less efficient replicators become extinct.


  We examine the evolution of expression patterns and the organization of
genetic information in populations of self-replicating digital organisms.
Seeding the experiments with a linearly expressed ancestor, we witness the
development of complex, parallel secondary expression patterns. Using
principles from information theory, we demonstrate an evolutionary pressure
towards overlapping expressions causing variation (and hence further evolution)
to sharply drop. Finally, we compare the overlapping sections of dominant
genomes to those portions which are singly expressed and observe a significant
difference in the entropy of their encoding.


  We introduce and analyze a general model of a population evolving over a
network of selectively neutral genotypes. We show that the population's limit
distribution on the neutral network is solely determined by the network
topology and given by the principal eigenvector of the network's adjacency
matrix. Moreover, the average number of neutral mutant neighbors per individual
is given by the matrix spectral radius. This quantifies the extent to which
populations evolve mutational robustness: the insensitivity of the phenotype to
mutations. Since the average neutrality is independent of evolutionary
parameters---such as, mutation rate, population size, and selective
advantage---one can infer global statistics of neutral network topology using
simple population data available from {\it in vitro} or {\it in vivo}
evolution. Populations evolving on neutral networks of RNA secondary structures
show excellent agreement with our theoretical predictions.


  We consider swarms formed by populations of self-propelled particles with
attractive long-range interactions. These swarms represent multistable
dynamical systems and can be found either in coherent traveling states or in an
incoherent oscillatory state where translational motion of the entire swarm is
absent. Under increasing the noise intensity, the coherent traveling state of
the swarms is destroyed and an abrupt transition to the oscillatory state takes
place.


  We analytically study the dynamics of evolving populations that exhibit
metastability on the level of phenotype or fitness. In constant selective
environments, such metastable behavior is caused by two qualitatively different
mechanisms. One the one hand, populations may become pinned at a local fitness
optimum, being separated from higher-fitness genotypes by a {\em fitness
barrier} of low-fitness genotypes. On the other hand, the population may only
be metastable on the level of phenotype or fitness while, at the same time,
diffusing over {\em neutral networks} of selectively neutral genotypes.
Metastability occurs in this case because the population is separated from
higher-fitness genotypes by an {\em entropy barrier}: The population must
explore large portions of these neutral networks before it discovers a rare
connection to fitter phenotypes.
  We derive analytical expressions for the barrier crossing times in both the
fitness barrier and entropy barrier regime. In contrast with ``landscape''
evolutionary models, we show that the waiting times to reach higher fitness
depend strongly on the width of a fitness barrier and much less on its height.
The analysis further shows that crossing entropy barriers is faster by orders
of magnitude than fitness barrier crossing. Thus, when populations are trapped
in a metastable phenotypic state, they are most likely to escape by crossing an
entropy barrier, along a neutral path in genotype space. If no such escape
route along a neutral path exists, a population is most likely to cross a
fitness barrier where the barrier is {\em narrowest}, rather than where the
barrier is shallowest.


  There are 512 two-locus, two-allele, two-phenotype, fully-penetrant disease
models. Using the permutation between two alleles, between two loci, and
between being affected and unaffected, one model can be considered to be
equivalent to another model under the corresponding permutation. These
permutations greatly reduce the number of two-locus models in the analysis of
complex diseases. This paper determines the number of non-redundant two-locus
models (which can be 102, 100, 96, 51, 50, or 48, depending on which
permutations are used, and depending on whether zero-locus and single-locus
models are excluded). Whenever possible, these non-redundant two-locus models
are classified by their property. Besides the familiar features of
multiplicative models (logical AND), heterogeneity models (logical OR), and
threshold models, new classifications are added or expanded: modifying-effect
models, logical XOR models, interference and negative interference models
(neither dominant nor recessive), conditionally dominant/recessive models,
missing lethal genotype models, and highly symmetric models. The following
aspects of two-locus models are studied: the marginal penetrance tables at both
loci, the expected joint identity-by-descent probabilities, and the correlation
between marginal identity-by-descent probabilities at the two loci. These
studies are useful for linkage analyses using single-locus models while the
underlying disease model is two-locus, and for correlation analyses using the
linkage signals at different locations obtained by a single-locus model.


  In the spirit of the many recent simple models of evolution inspired by
statistical physics, we put forward a simple model of the evolution of such
models. Like its objects of study, it is (one supposes) in principle testable
and capable of making predictions, and gives qualitative insights into a
hitherto mysterious process.


  We examine the dynamics of an age-structured population model in which the
life expectancy of an offspring may be mutated with respect to that of the
parent. While the total population of the system always reaches a steady state,
the fitness and age characteristics exhibit counter-intuitive behavior as a
function of the mutational bias. By analytical and numerical study of the
underlying rate equations, we show that if deleterious mutations are favored,
the average fitness of the population reaches a steady state, while the average
population age is a decreasing function of the overall fitness. When
advantageous mutations are favored, the average population fitness grows
linearly with time t, while the average age is independent of fitness. For no
mutational bias, the average fitness grows as t^{2/3}.


  (Abridged) We show that sizeable fractions of incident ionizing radiation
from stochastic astrophysical sources can be redistributed to biologically and
chemically important UV wavelengths, a significant fraction of which can reach
the surface. This redistribution is mediated by secondary electrons, resulting
from Compton scattering and X-ray photoabsorption, with energies low enough to
excite atmospheric molecules and atoms, resulting in a rich aurora-like
spectrum. We calculate the fraction of energy redistributed into biologically
and chemically important wavelength regions for spectra characteristic of
stellar flares and supernovae using a Monte-Carlo transport code written for
this problem and then estimate the fraction of this energy that is transmitted
from the atmospheric altitudes of redistribution to the surface for a few
illustrative cases. Redistributed fractions are found to be of order 1%, even
in the presence of an ozone shield. This result implies that planetary
organisms will be subject to mutationally significant, if intermittent,
fluences of UV-B and harder radiation even in the presence of a narrow-band UV
shield like ozone. We also calculate the surficial transmitted fraction of
ionizing radiation and redistributed ultraviolet radiation for two illustrative
evolving Mars atmospheres whose initial surface pressures were 1 bar. Our
results suggest that coding organisms on planets orbiting low-mass stars (and
on the early Earth) may evolve very differently than on contemporary Earth,
with diversity and evolutionary rate controlled by a stochastically varying
mutation rate and frequent hypermutation episodes.


  (Abridged) The propagation of ionizing radiation through model atmospheres of
terrestrial-like exoplanets is studied for a large range of column densities
and incident photon energies using a Monte Carlo code we have developed to
treat Compton scattering and photoabsorption. Incident spectra from parent star
flares, supernovae, and gamma-ray bursts are modeled and compared to energetic
particles in importance. We find that terrestrial-like exoplanets with
atmospheres thinner than about 100 g cm^-2 transmit and reprocess a significant
fraction of incident gamma-rays, producing a characteristic, flat surficial
spectrum. Thick atmospheres (>~ 100 g cm^-2) efficiently block even gamma-rays,
but nearly all incident energy is redistributed into diffuse UV and visible
aurora-like emission, increasing the effective atmospheric transmission by many
orders of magnitude. Depending on the presence of molecular UV absorbers and
atmospheric thickness, up to 10% of the incident energy can reach the surface
as UV reemission. For the Earth, between 2 x 10^-3 and 4 x 10^-2 of the
incident flux reaches the ground in the biologically effective 200--320 nm
range, depending on O_2/O_3 shielding. Finally, we suggest that transient
atmospheric ionization layers can be frequently created at low altitudes. We
conclude that these events can produce frequent fluctuations in atmospheric
ionization levels and surficial UV fluxes on terrestrial-like planets.


  Gamma-ray bursts (hereafter GRB) produce a flux of radiation detectable
across the observable Universe, and at least some of them are associated with
galaxies. A GRB within our own Ggalaxy could do considerable damage to the
Earth's biosphere; rate estimates suggest that a dangerously near GRB should
occur on average two or more times per billion years. At least five times in
the history of life, the Earth experienced mass extinctions that eliminated a
large percentage of the biota. Many possible causes have been documented, and
GRB may also have contributed. The late Ordovician mass extinction
approximately 440 million years ago may be at least partly the result of a GRB.
A special feature of GRB in terms of terrestrial effects is a nearly impulsive
energy input of order 10 s. Due to expected severe depletion of the ozone
layer, intense solar ultraviolet radiation would result from a nearby GRB, and
some of the patterns of extinction and survivorship at this time may be
attributable to elevated levels of UV radiation reaching the Earth. In addition
a GRB could trigger the global cooling which occurs at the end of the
Ordovician period that follows an interval of relatively warm climate. Intense
rapid cooling and glaciation at that time, previously identified as the
probable cause of this mass extinction, may have resulted from a GRB.


  Based on cosmological rates, it is probable that at least once in the last Gy
the Earth has been irradiated by a gamma-ray burst in our Galaxy from within 2
kpc. Using a two-dimensional atmospheric model we have performed the first
computation of the effects upon the Earth's atmosphere of one such impulsive
event. A ten second burst delivering 100 kJ/m^2 to the Earth penetrates to the
stratosphere and results in globally averaged ozone depletion of 35%, with
depletion reaching 55% at some latitudes. Significant global depletion persists
for over 5 years after the burst. This depletion would have dramatic
implications for life since a 50% decrease in ozone column density results in
approximately three times the normal UVB flux. Widespread extinctions are
likely, based on extrapolation from UVB sensitivity of modern organisms.
Additional effects include a shot of nitrate fertilizer and NO2 opacity in the
visible providing a cooling perturbation to the climate over a similar
timescale. These results lend support to the hypothesis that a GRB may have
initiated the late Ordovician mass extinction (Melott et al. 2004).


  It is likely that one or more gamma-ray bursts within our galaxy have
strongly irradiated the Earth in the last Gy. This produces significant
atmospheric ionization and dissociation, resulting in ozone depletion and
DNA-damaging ultraviolet solar flux reaching the surface for up to a decade.
Here we show the first detailed computation of two other significant effects.
Visible opacity of NO2 is sufficient to reduce solar energy at the surface up
to a few percent, with the greatest effect at the poles, which may be
sufficient to initiate glaciation. Rainout of dilute nitric acid is could have
been important for a burst nearer than our conservative nearest burst. These
results support the hypothesis that the characteristics of the late Ordovician
mass extinction are consistent with GRB initiation.


  The development of prebiotic homochirality on early-Earth or another
planetary platform may be viewed as a critical phenomenon. It is shown, in the
context of spatio-temporal polymerization reaction networks, that environmental
effects -- be them temperature surges or other external disruptions -- may
destroy any net chirality previously produced. In order to understand the
emergence of prebiotic homochirality it is important to model the coupling of
polymerization reaction networks to different planetary environments.


  We describe results of modeling the effects on Earth-like planets of
long-duration gamma-ray bursts (GRBs) within a few kiloparsecs. A primary
effect is generation of nitrogen oxide compounds which deplete ozone. Ozone
depletion leads to an increase in solar UVB radiation at the surface, enhancing
DNA damage, particularly in marine microorganisms such as phytoplankton. In
addition, we expect increased atmospheric opacity due to buildup of nitrogen
dioxide produced by the burst and enhanced precipitation of nitric acid. We
review here previous work on this subject and discuss recent developments,
including further discussion of our estimates of the rates of impacting GRBs
and the possible role of short-duration bursts.


  Recent work has revealed a 62 (+/-) 3-million-year cycle in the fossil
diversity in the past 542 My, however no plausible mechanism has been found. We
propose that the cycle may be caused by modulation of cosmic ray (CR) flux by
the Solar system vertical oscillation (64 My period) in the galaxy, the
galactic north-south anisotropy of CR production in the galactic
halo/wind/termination shock (due to the galactic motion toward the Virgo
cluster), and the shielding by galactic magnetic fields. We revisit the
mechanism of CR propagation and show that CR flux can vary by a factor of about
4.6 and reach a maximum at north-most displacement of the Sun. The very high
statistical significance of (i) the phase agreement between Solar north-ward
excursions and the diversity minima and (ii) the correlation of the magnitude
of diversity drops with CR amplitudes through all cycles provide solid support
for our model. Various observational predictions which can be used to confirm
or falsify our hypothesis are presented.


  Stanek et al. (astro-ph/0604113) have noted that the four low-redshift
long-duration gamma-ray bursts (LSB) observed to date all occurred in faint,
metal-poor galaxies. Given this selection, they argue that it is improbable
that there has been a substantial population of Milky Way galaxy bursts
sufficiently recently to affect life on Earth. This argument ignores the
heterogeneity of stellar populations in the Milky Way, with evidence for
continuing mergers with low-metallicity dwarf galaxies; observational analysis
that points to LSBs being hosted by such galaxies undergoing interaction; and
the existence of a likely recent GRB remnant in our galaxy.


  Gamma-Ray Bursts (GRBs) directed at Earth from within a few kpc may have
damaged the biosphere, primarily though changes in atmospheric chemistry which
admit greatly increased Solar UV. However, GRBs are highly variable in spectrum
and duration. Recent observations indicate that short (~0.1 s) burst GRBs,
which have harder spectra, may be sufficiently abundant at low redshift that
they may offer an additional significant effect. A much longer timescale is
associated with shock breakout luminosity observed in the soft X-ray (~10^3 s)
and UV (~10^5 s) emission, and radioactive decay gamma-ray line radiation
emitted during the light curve phase of supernovae (~10^7 s). Here we
generalize our atmospheric computations to include a broad range of peak photon
energies and investigate the effect of burst duration while holding total
fluence and other parameters constant. The results can be used to estimate the
probable impact of various kinds of ionizing events (such as short GRBs, X-ray
flashes, supernovae) upon the terrestrial atmosphere. We find that the ultimate
intensity of atmospheric effects varies only slightly with burst duration from
10^-1 s to 10^8 s. Therefore, the effect of many astrophysical events causing
atmospheric ionization can be approximated without including time development.
Detailed modeling requires specification of the season and latitude of the
event. Harder photon spectra produce greater atmospheric effects for spectra
with peaks up to about 20 MeV, because of greater penetration into the
stratosphere.


  We have examined the Fourier power spectrum as well as the Hurst exponent of
extinction, origination, and total biodiversity in the marine fossil record,
using a recently improved geologic timescale. We find all of them strongly
inconsistent with past claims of self-similarity as well as inconsistent with
random walk behavior. Instead, they are dominated by low-frequency power, with
approximate f^-2 power over one decade in frequency. The spectrum turns over at
about 10^8 y, lending plausibility to connections with galactic dynamics. Even
in the background of this low-frequency dominance, a previously noted 62 My
biodiversity cycle stands out with better than 99% confidence above the noise
level, accounting for about 35% of the total variance in the fossil
biodiversity record.


  The role of asymmetry on the evolution of prebiotic homochirality is
investigated in the context of autocatalytic polymerization reaction networks.
A model featuring enantiometric cross-inhibition and chiral bias is used to
study the diffusion equations controlling the spatiotemporal development of
left and right-handed domains. Bounds on the chiral bias are obtained
consistent with present-day constraints on the emergence of life on early
Earth. The viability of biasing mechanisms such as weak neutral currents and
circularly polarized UV light is discussed. The results can be applied to any
hypothetical planetary platform.


  In this paper, we review the literature on statistical long-range correlation
in DNA sequences. We examine the current evidence for these correlations, and
conclude that a mixture of many length scales (including some relatively long
ones) in DNA sequences is responsible for the observed 1/f-like spectral
component. We note the complexity of the correlation structure in DNA
sequences. The observed complexity often makes it hard, or impossible, to
decompose the sequence into a few statistically stationary regions. We suggest
that, based on the complexity of DNA sequences, a fruitful approach to
understand long-range correlation is to model duplication, and other
rearrangement processes, in DNA sequences. One model, called
``expansion-modification system", contains only point duplication and point
mutation. Though simplistic, this model is able to generate sequences with 1/f
spectra. We emphasize the importance of DNA duplication in its contribution to
the observed long-range correlation in DNA sequences.


  We study the spatial patterns formed by interacting populations or reacting
chemicals under the influence of chaotic flows. In particular, we have
considered a three-component model of plankton dynamics advected by a
meandering jet. We report general results, stressing the existence of a
smooth-filamental transition in the concentration patterns depending on the
relative strength of the stirring by the chaotic flow and the relaxation
properties of planktonic dynamical system. Patterns obtained in open and closed
flows are compared.


  Integral membrane proteins deform the surrounding bilayer creating
long-ranged forces that influence distant proteins. These forces can be
attractive or repulsive, depending on the proteins' shape, height, contact
angle with the bilayer, as well as the local membrane curvature. Although
interaction energies are not pairwise additive, for sufficiently low protein
density, thermodynamic properties depend only upon pair interactions. Here, we
compute pair interaction potentials and entropic contributions to the
two-dimensional osmotic pressure of a collection of noncircular proteins. In
contrast to direct short-ranged interactions such as van der Waal's,
hydrophobic, or electrostatic interactions, both local membrane Gaussian
curvature and protein ellipticity can induce attractions between two proteins
at distances of up to ten times their typical radii. For flat membranes,
bending rigidities of $\sim 30k_{B}T$, and moderate ellipticities, we find
thermally averaged attractive interactions of order $\sim 2k_{B}T$. These
interactions may play an important role in the intermediate stages of protein
aggregation.


  We study the population profile in a simple discrete time model of population
dynamics. Our model, which is closely related to certain ``bit-string'' models
of evolution, incorporates competition for resources via a population dependent
death probability, as well as a variable reproduction probability for each
individual as a function of age. We first solve for the steady-state of the
model in mean field theory, before developing analytic techniques to compute
Gaussian fluctuation corrections around the mean field fixed point. Our
computations are found to be in good agreement with Monte-Carlo simulations.
Finally we discuss how similar methods may be applied to fluctuations in
continuous time population models.