  The epidermis is a specialized epithelium that constitutes the outermost
layer of the skin, and it provides a protective barrier against environmental
assaults. Primarily consisting of multilayered keratinocytes, the epidermis is
continuously renewed by proliferation of stem cells and the differentiation of
their progeny, which undergo terminal differentiation as they leave the basal
layer and move upward toward the surface, where they die and slough off. Basal
keratinocytes rest on a basement membrane at the dermal-epidermal junction that
is composed of specific extracellular matrix proteins organized into
interactive and mechanically supportive networks. Firm attachment of basal
keratinocytes, and their dynamic regulation via focal adhesions and
hemidesmosomes, is essential for maintaining major skin processes, such as
self-renewal, barrier function, and resistance to physical and chemical
stresses. The adhesive integrin receptors expressed by epidermal cells serve
structural, signaling, and mechanosensory roles that are critical for epidermal
cell anchorage and tissue homeostasis. More specifically, the basement membrane
components play key roles in preserving the stem cell pool, and establishing
cell polarity cues enabling asymmetric cell divisions, which result in the
transition from a proliferative basal cell layer to suprabasal cells committed
to terminal differentiation. Finally, through a well-regulated sequence of
synthesis and remodeling, the components of the dermal-epidermal junction play
an essential role in regeneration of the epidermis during skin healing. Here
too, they provide biological and mechanical signals that are essential to the
restoration of barrier function.


  We propose a novel approach to detect the binding between proteins making use
of the anomalous diffraction of natively present heavy elements inside the
molecule 3D structure. In particular, we suggest considering sulfur atoms
contained in protein structures at lower percentages than the other atomic
species. Here, we run an extensive preliminary investigation to probe both the
feasibility and the range of usage of the proposed protocol. In particular, we
(i) analytically and numerically show that the diffraction patterns produced by
the anomalous scattering of the sulfur atoms in a given direction depend
additively on the relative distances between all couples of sulfur atoms. Thus
the differences in the patterns produced by bound proteins with respect to
their non-bonded states can be exploited to rapidly assess protein complex
formation. Next, we (ii) carried out analyses on the abundances of sulfurs in
the different proteomes and molecular dynamics simulations on a representative
set of protein structures to probe the typical motion of sulfur atoms. Finally,
we (iii) suggest a possible experimental procedure to detect protein-protein
binding. Overall, the completely label-free and rapid method we propose may be
readily extended to probe interactions on a large scale even between other
biological molecules, thus paving the way to the development of a novel field
of research based on a synchrotron light source.


  In this paper, we introduce and study an age-structured epidemiological
compartment model and its respective problem, applied but not limited to the
COVID-19 pandemic, in order to investigate the role of the age of the
individuals in the evolution of epidemiological phenomena. We investigate the
well-posedness of the model, as well as the global dynamics of it in the sense
of basic reproduction number via constructing Lyapunov functions.


  Electron cryomicroscopy (cryo-EM) is an imaging technique widely used in
structural biology to determine the three-dimensional structure of biological
molecules from noisy two-dimensional projections with unknown orientations. As
the typical pipeline involves processing large amounts of data, efficient
algorithms are crucial for fast and reliable results. The stochastic gradient
descent (SGD) algorithm has been used to improve the speed of ab initio
reconstruction, which results in a first, low-resolution estimation of the
volume representing the molecule of interest, but has yet to be applied
successfully in the high-resolution regime, where expectation-maximization
algorithms achieve state-of-the-art results, at a high computational cost. In
this article, we investigate the conditioning of the optimization problem and
show that the large condition number prevents the successful application of
gradient descent-based methods at high resolution. Our results include a
theoretical analysis of the condition number of the optimization problem in a
simplified setting where the individual projection directions are known, an
algorithm based on computing a diagonal preconditioner using Hutchinson's
diagonal estimator, and numerical experiments showing the improvement in the
convergence speed when using the estimated preconditioner with SGD. The
preconditioned SGD approach can potentially enable a simple and unified
approach to ab initio reconstruction and high-resolution refinement with faster
convergence speed and higher flexibility, and our results are a promising step
in this direction.


  Therapeutic antibodies are an essential and rapidly expanding drug modality.
The binding specificity between antibodies and antigens is decided by
complementarity-determining regions (CDRs) at the tips of these Y-shaped
proteins. In this paper, we propose a hierarchical training paradigm (HTP) for
the antibody sequence-structure co-design. HTP consists of four levels of
training stages, each corresponding to a specific protein modality within a
particular protein domain. Through carefully crafted tasks in different stages,
HTP seamlessly and effectively integrates geometric graph neural networks
(GNNs) with large-scale protein language models to excavate evolutionary
information from not only geometric structures but also vast antibody and
non-antibody sequence databases, which determines ligand binding pose and
strength. Empirical experiments show that HTP sets the new state-of-the-art
performance in the co-design problem as well as the fix-backbone design. Our
research offers a hopeful path to unleash the potential of deep generative
architectures and seeks to illuminate the way forward for the antibody sequence
and structure co-design challenge.


  RNA protein Interactions (RPIs) play an important role in biological systems.
Recently, we have enumerated the RPIs at the residue level and have elucidated
the minimum structural unit (MSU) in these interactions to be a stretch of five
residues (Nucleotides/amino acids). Pseudouridine is the most frequent
modification in RNA. The conversion of uridine to pseudouridine involves
interactions between pseudouridine synthase and RNA. The existing models to
predict the pseudouridine sites in a given RNA sequence mainly depend on
user-defined features such as mono and dinucleotide composition/propensities of
RNA sequences. Predicting pseudouridine sites is a non-linear classification
problem with limited data points. Deep Learning models are efficient
discriminators when the data set size is reasonably large and fail when there
is a paucity of data ($<1000$ samples). To mitigate this problem, we propose a
Support Vector Machine (SVM) Kernel based on utility theory from Economics, and
using data-driven parameters (i.e. MSU) as features. For this purpose, we have
used position-specific tri/quad/pentanucleotide composition/propensity
(PSPC/PSPP) besides nucleotide and dineculeotide composition as features. SVMs
are known to work well in small data regimes and kernels in SVM are designed to
classify non-linear data. The proposed model outperforms the existing
state-of-the-art models significantly (10%-15% on average).


  The binding between proteins and ligands plays a crucial role in the realm of
drug discovery. Previous deep learning approaches have shown promising results
over traditional computationally intensive methods, but resulting in poor
generalization due to limited supervised data. In this paper, we propose to
learn protein-ligand binding representation in a self-supervised learning
manner. Different from existing pre-training approaches which treat proteins
and ligands individually, we emphasize to discern the intricate binding
patterns from fine-grained interactions. Specifically, this self-supervised
learning problem is formulated as a prediction of the conclusive binding
complex structure given a pocket and ligand with a Transformer based
interaction module, which naturally emulates the binding process. To ensure the
representation of rich binding information, we introduce two pre-training
tasks, i.e.~atomic pairwise distance map prediction and mask ligand
reconstruction, which comprehensively model the fine-grained interactions from
both structure and feature space. Extensive experiments have demonstrated the
superiority of our method across various binding tasks, including
protein-ligand affinity prediction, virtual screening and protein-ligand
docking.


  The package mvlearnR and accompanying Shiny App is intended for integrating
data from multiple sources or views or modalities (e.g. genomics, proteomics,
clinical and demographic data). Most existing software packages for multiview
learning are decentralized and offer limited capabilities, making it difficult
for users to perform comprehensive integrative analysis. The new package wraps
statistical and machine learning methods and graphical tools, providing a
convenient and easy data integration workflow. For users with limited
programming language, we provide a Shiny Application to facilitate data
integration anywhere and on any device. The methods have potential to offer
deeper insights into complex disease mechanisms.
  Availability and Implementation: mvlearnR is available from the following
GitHub repository: https://github.com/lasandrall/mvlearnR. The web application
is hosted on shinyapps.io and available at:
https://multi-viewlearn.shinyapps.io/MultiView_Modeling/


  We present Symphony, an $E(3)$-equivariant autoregressive generative model
for 3D molecular geometries that iteratively builds a molecule from molecular
fragments. Existing autoregressive models such as G-SchNet and G-SphereNet for
molecules utilize rotationally invariant features to respect the 3D symmetries
of molecules. In contrast, Symphony uses message-passing with higher-degree
$E(3)$-equivariant features. This allows a novel representation of probability
distributions via spherical harmonic signals to efficiently model the 3D
geometry of molecules. We show that Symphony is able to accurately generate
small molecules from the QM9 dataset, outperforming existing autoregressive
models and approaching the performance of diffusion models.


  The rapid evolution of artificial intelligence in drug discovery encounters
challenges with generalization and extensive training, yet Large Language
Models (LLMs) offer promise in reshaping interactions with complex molecular
data. Our novel contribution, InstructMol, a multi-modal LLM, effectively
aligns molecular structures with natural language via an instruction-tuning
approach, utilizing a two-stage training strategy that adeptly combines limited
domain-specific data with molecular and textual information. InstructMol
showcases substantial performance improvements in drug discovery-related
molecular tasks, surpassing leading LLMs and significantly reducing the gap
with specialized models, thereby establishing a robust foundation for a
versatile and dependable drug discovery assistant.


  In a previous paper, we showed that a compartmental stochastic process model
of SARS-CoV-2 transmission could be fit to time series data and then
reinterpreted as a collection of interacting branching processes drawn from a
dynamic degree distribution. We called this reinterpretation a transmission
forest. This paper builds on that idea. Specifically, leveraging generating
function methods from analytic combinatorics, we develop a theory describing
the transmission forest's properties, allowing us to show for example that
transmission tree interactions fade with increasing disease prevalence. We then
validate the theory by computing forest statistics, like the tree survival
function, which we compare to estimates based on the sampling method developed
previously. The accuracy and flexibility of the analytic approach is clear, and
it allows us to comment on multi-scale features of more general transmission
processes.


  The quest for accurate prediction of drug molecule properties poses a
fundamental challenge in the realm of Artificial Intelligence Drug Discovery
(AIDD). An effective representation of drug molecules emerges as a pivotal
component in this pursuit. Contemporary leading-edge research predominantly
resorts to self-supervised learning (SSL) techniques to extract meaningful
structural representations from large-scale, unlabeled molecular data,
subsequently fine-tuning these representations for an array of downstream
tasks. However, an inherent shortcoming of these studies lies in their singular
reliance on one modality of molecular information, such as molecule image or
SMILES representations, thus neglecting the potential complementarity of
various molecular modalities. In response to this limitation, we propose MolIG,
a novel MultiModaL molecular pre-training framework for predicting molecular
properties based on Image and Graph structures. MolIG model innovatively
leverages the coherence and correlation between molecule graph and molecule
image to execute self-supervised tasks, effectively amalgamating the strengths
of both molecular representation forms. This holistic approach allows for the
capture of pivotal molecular structural characteristics and high-level semantic
information. Upon completion of pre-training, Graph Neural Network (GNN)
Encoder is used for the prediction of downstream tasks. In comparison to
advanced baseline models, MolIG exhibits enhanced performance in downstream
tasks pertaining to molecular property prediction within benchmark groups such
as MoleculeNet Benchmark Group and ADMET Benchmark Group.


  The basic reproduction number R0 is a concept which originated in population
dynamics, mathematical epidemiology, and ecology and is closely related to the
mean number of children in branching processes.We offer below three new
contributions to the literature: 1) We order a universal algorithmic definition
of a (F, V) gradient decomposition (and hence of the resulting R0), which
requires a minimal input from the user, namely the specification of an
admissible set of disease/infection variables. We also present examples where
other choices may be more reasonable, with more terms in F, or more terms in V .
  2) We glean out from the works of Bacaer a fixed point equation (8) for the
extinction probabilities of a stochastic model associated to a deterministic
ODE model, which may be expressed in terms of the (F, V ) decomposition. The
fact that both R0 and the extinction probabilities are functions of (F, V )
underlines the centrality of this pair, which may be viewed as more fundamental
than the famous next generation matrix FV^{-1}.
  3) We suggest introducing a new concept of sufficient/minimal
disease/infection set (sufficient for determining R0). More precisely, our
universal recipe of choosing "new infections" once the "infections" are
specified suggests focusing on the choice of the latter, which is also not
unique. The maximal choice of choosing all compartments which become 0 at the
given boundary point seems to always work, but is the least useful for analytic
computations, therefore we propose to investigate the minimal one. As a bonus,
this idea seems useful for understanding the Jacobian factorization approach
for computing R0 . Last but not least, we offer Mathematica scripts and
implement them for a large variety of examples, which illustrate that our
recipe others always reasonable results, but that sometimes other reasonable
(F, V ) decompositions are available as well.


  A closed formula multiallelic Walsh (or Hadamard) transform is introduced.
Basic results are derived, and a statistical interpretation of some of the
resulting linear forms is discussed.


  Molecular dynamics (MD) simulations and computer-aided drug design (CADD)
have advanced substantially over the past two decades, thanks to continuous
computer hardware and software improvements. Given these advancements, MD
simulations are poised to become even more powerful tools for investigating the
dynamic interactions between potential small-molecule drugs and their target
proteins, with significant implications for pharmacological research.


  Single-cell multi-view clustering enables the exploration of cellular
heterogeneity within the same cell from different views. Despite the
development of several multi-view clustering methods, two primary challenges
persist. Firstly, most existing methods treat the information from both
single-cell RNA (scRNA) and single-cell Assay of Transposase Accessible
Chromatin (scATAC) views as equally significant, overlooking the substantial
disparity in data richness between the two views. This oversight frequently
leads to a degradation in overall performance. Additionally, the majority of
clustering methods necessitate manual specification of the number of clusters
by users. However, for biologists dealing with cell data, precisely determining
the number of distinct cell types poses a formidable challenge. To this end, we
introduce scUNC, an innovative multi-view clustering approach tailored for
single-cell data, which seamlessly integrates information from different views
without the need for a predefined number of clusters. The scUNC method
comprises several steps: initially, it employs a cross-view fusion network to
create an effective embedding, which is then utilized to generate initial
clusters via community detection. Subsequently, the clusters are automatically
merged and optimized until no further clusters can be merged. We conducted a
comprehensive evaluation of scUNC using three distinct single-cell datasets.
The results underscored that scUNC outperforms the other baseline methods.


  In recent years, the field of single-cell data analysis has seen a marked
advancement in the development of clustering methods. Despite advancements,
most of these algorithms still concentrate on analyzing the provided
single-cell matrix data. However, in medical applications, single-cell data
often involves a wealth of exogenous information, including gene networks.
Overlooking this aspect could lead to information loss and clustering results
devoid of significant clinical relevance. An innovative single-cell deep
clustering method, incorporating exogenous gene information, has been proposed
to overcome this limitation. This model leverages exogenous gene network
information to facilitate the clustering process, generating discriminative
representations. Specifically, we have developed an attention-enhanced graph
autoencoder, which is designed to efficiently capture the topological features
between cells. Concurrently, we conducted a random walk on an exogenous
Protein-Protein Interaction (PPI) network, thereby acquiring the gene's
topological features. Ultimately, during the clustering process, we integrated
both sets of information and reconstructed the features of both cells and genes
to generate a discriminative representation. Extensive experiments have
validated the effectiveness of our proposed method. This research offers
enhanced insights into the characteristics and distribution of cells, thereby
laying the groundwork for early diagnosis and treatment of diseases.


  Scope: Biomarkers of food intake (BFIs) are useful tools for objective
assessment of food intake and compliance. The aim of this study was to discover
and identify urinary BFIs for onion. Methods and results: In a randomized
controlled cross-over trial, 6 overweight participants (age 24-62 years)
consumed meals with 20 g/d onion powder or no onion for 2 weeks. Untargeted
UPLC-qTOF-MS metabolic profiling analysis was performed on urine samples and
the profiles were analysed by multilevel-PLSDA, modified PLS, and nearest
shrunken centroid to select features associated with onion intake. Eight
biomarkers were tentatively identified; six of them originated from
S-substituted cysteine derivatives such as isoalliin and propiin, which are
considered the most specific for onion intake. Most of the biomarkers were
completely excreted within 24 hours and no accumulation was observed during 2
weeks indicating their ability to reflect only recent intake of onions.
Receiver-operator curves were made to evaluate the performance of individual
biomarkers for predicting onion intake. The area under the curve values for
these biomarkers ranged from 0.81 to 1. Conclusion: Promising biomarkers of
recent onion intake have been identified in human urine. Further studies with
complex diets are needed to validate the robustness of these biomarkers.


  As the impending consequences of climate change loom over the Earth, it has
become vital for researchers to understand the role microorganisms play in this
process. In this paper, we examine how environmental factors, including
moisture levels and temperature, affect the expression of certain fungal
characteristics on a microscale, and how these in turn affect fungal
biodiversity and ecosystem decomposition rates over time.
  We first present a differential equation model to understand how the
distribution of different fungal isolates depends on regional moisture levels.
We introduce both slow and sudden variations into the environment in order to
represent the various ways climate change will impact fungal ecosystems. This
model demonstrates that increased variability in moisture (both short-term and
long-term) increases biodiversity and that fungal populations will shift
towards more stress-tolerant fungi as aridity increases. The model further
suggests the lack of any direct link between biodiversity and decomposition
rates.
  To better describe fungal competition with respect to space, we develop a
local agent-based model (ABM). Unlike the previous model, our ABM focuses on
individuals, tracking each fungus and the result of its interactions. Our ABM
also features a more accurate spatial combat system, allowing us to precisely
discern the influence of fungal interactions on the environment. This model
corroborates the results of the differential equation model and further
suggests that moisture, through its link with temperature and effects on fungal
population, also plays a strong role in determining fungal decomposition rates.
  Together, these models suggest that climate change, which portends increasing
variability in regional conditions and higher average temperatures worldwide,
will lead to an increase in both wood decomposition rates and, independently,
fungal biodiversity.


  Biological evolution can be conceptualized as a search process in the space
of gene sequences guided by the fitness landscape, a mapping that assigns a
measure of reproductive value to each genotype. Here we discuss probabilistic
models of fitness landscapes with a focus on their evolutionary accessibility,
where a path in a fitness landscape is said to be accessible if the fitness
values encountered along the path increase monotonically. For uncorrelated
(random) landscapes with independent and identically distributed fitness
values, the probability of existence of accessible paths between genotypes at a
distance linear in the sequence length $L$ becomes nonzero at a nontrivial
threshold value of the fitness difference between the initial and final
genotype, which can be explicitly computed for large classes of genotype
graphs. The behaviour in uncorrelated random landscapes is contrasted with
landscape models that display additional, biologically motivated structural
features. In particular, landscapes defined by a tradeoff between adaptation to
environmental extremes have been found to display a combinatorially large
number of accessible paths to all local fitness maxima. We show that this
property is characteristic of a broad class of models that satisfy a certain
global constraint, and provide further examples from this class.


  The translation of comparative genomics into clinical decision support tools
often depends on the quality of sequence alignments. However, currently used
methods of multiple sequence alignments suffer from significant biases and
problems with aligning diverged sequences. The objective of this study was to
develop and test a new multiple sequence alignment (MSA) algorithm suitable for
the high-throughput comparative analysis of different microbial genomes. This
algorithm employs an innovative tensor indexing method for partitioning the
dynamic programming hyper-cube space for parallel processing. We have used the
clinically relevant task of identifying regions that determine resistance to
antibiotics to test the new algorithm and to compare its performance with
existing MSA methods. The new method "mmDst" performed better than existing MSA
algorithms for more divergent sequences because it employs a simultaneous
alignment scoring recurrence, which effectively approximated the score for edge
missing cell scores that fall outside the scoring region.


  In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss
data requirements for deposition and validation of cryoEM structures, with a
focus on single-particle analysis. The meeting was attended by 45 experts in
data processing, model building and refinement, validation, and archiving of
such structures. This report describes the workshop's motivation and history,
the topics discussed, and consensus recommendations resulting from the
workshop. Some challenges for future methods-development efforts in this area
are also highlighted, as is the implementation to date of some of the
recommendations.


  Escherichia coli bacterium is a rod-shaped organism composed of a complex
double membrane structure. Knowledge of electric field driven ion transport
through both membranes and the evolution of their induced permeabilization has
important applications in biomedical engineering, delivery of genes and
antibacterial agents. However, few studies have been conducted on Gram-negative
bacteria in this regard considering the contribution of all ion types. To
address this gap in knowledge, we have developed a deterministic and stochastic
Brownian dynamics model to simulate in 3D space the motion of ions through
pores formed in the plasma membranes of E. coli cells during electroporation.
The diffusion coefficient, mobility, and translation time of Ca$^{2+}$,
Mg$^{2+}$, Na$^+$, K$^+$, and Cl$^-$ ions within the pore region are estimated
from the numerical model. Calculations of pore's conductance have been
validated with experiments conducted at Gustave Roussy. From the simulations,
it was found that the main driving force of ionic uptake during the pulse is
the one due to the externally applied electric field. The results from this
work provide a better understanding of ion transport during electroporation,
aiding in the design of electrical pulses for maximizing ion throughput,
primarily for application in cancer treatment.


  Small molecule drugs have an important role to play in combating viral
infections, and biophysics support has been central for contributing to the
discovery and design of direct acting antivirals. Perhaps one of the most
successful biophysical tools for this purpose is NMR spectroscopy when utilized
strategically and pragmatically within team workflows and timelines. This
report describes some clear examples of how NMR applications contributed to the
design of antivirals when combined with medicinal chemistry, biochemistry,
X-ray crystallography and computational chemistry. Overall, these
multidisciplinary approaches allowed teams to reveal and expose compound
physical properties from which design ideas were spawned and tested to achieve
the desired successes. Examples are discussed for the discovery of antivirals
that target HCV, HIV and SARS-CoV-2.


  The aim of this study was to develop a method that would identify the cluster
centroids and the optimal number of clusters for a given sensitivity level and
could work equally well for the different sequence datasets. A novel method
that combines the linear mapping hash function and multiple sequence alignment
(MSA) was developed. This method takes advantage of the already sorted by
similarity sequences from the MSA output, and identifies the optimal number of
clusters, clusters cut-offs, and clusters centroids that can represent
reference gene vouchers for the different species. The linear mapping hash
function can map an already ordered by similarity distance matrix to indices to
reveal gaps in the values around which the optimal cut-offs of the different
clusters can be identified. The method was evaluated using sets of closely
related (16S rRNA gene sequences of Nocardia species) and highly variable (VP1
genomic region of Enterovirus 71) sequences and outperformed existing
unsupervised machine learning clustering methods and dimensionality reduction
methods. This method does not require prior knowledge of the number of clusters
or the distance between clusters, handles clusters of different sizes and
shapes, and scales linearly with the dataset. The combination of MSA with the
linear mapping hash function is a computationally efficient way of gene
sequence clustering and can be a valuable tool for the assessment of
similarity, clustering of different microbial genomes, identifying reference
sequences, and for the study of evolution of bacteria and viruses.


  The intra- and inter-species genetic diversity of bacteria and the absence of
'reference', or the most representative, sequences of individual species
present a significant challenge for sequence-based identification. The aims of
this study were to determine the utility, and compare the performance of
several clustering and classification algorithms to identify the species of 364
sequences of 16S rRNA gene with a defined species in GenBank, and 110 sequences
of 16S rRNA gene with no defined species, all within the genus Nocardia. A
total of 364 16S rRNA gene sequences of Nocardia species were studied. In
addition, 110 16S rRNA gene sequences assigned only to the Nocardia genus level
at the time of submission to GenBank were used for machine learning
classification experiments. Different clustering algorithms were compared with
a novel algorithm or the linear mapping (LM) of the distance matrix. Principal
Components Analysis was used for the dimensionality reduction and
visualization. Results: The LM algorithm achieved the highest performance and
classified the set of 364 16S rRNA sequences into 80 clusters, the majority of
which (83.52%) corresponded with the original species. The most representative
16S rRNA sequences for individual Nocardia species have been identified as
'centroids' in respective clusters from which the distances to all other
sequences were minimized; 110 16S rRNA gene sequences with identifications
recorded only at the genus level were classified using machine learning
methods. Simple kNN machine learning demonstrated the highest performance and
classified Nocardia species sequences with an accuracy of 92.7% and a mean
frequency of 0.578.


  Transcription factors (TFs) play a vital role in the regulation of gene
expression thereby making them critical to many cellular processes. In this
study, we used graph machine learning methods to create a compendium of TF
cascades using data extracted from the STRING database. A TF cascade is a
sequence of TFs that regulate each other, forming a directed path in the TF
network. We constructed a knowledge graph of 81,488 unique TF cascades, with
the longest cascade consisting of 62 TFs. Our results highlight the complex and
intricate nature of TF interactions, where multiple TFs work together to
regulate gene expression. We also identified 10 TFs with the highest regulatory
influence based on centrality measurements, providing valuable information for
researchers interested in studying specific TFs. Furthermore, our pathway
enrichment analysis revealed significant enrichment of various pathways and
functional categories, including those involved in cancer and other diseases,
as well as those involved in development, differentiation, and cell signaling.
The enriched pathways identified in this study may have potential as targets
for therapeutic intervention in diseases associated with dysregulation of
transcription factors. We have released the dataset, knowledge graph, and
graphML methods for the TF cascades, and created a website to display the
results, which can be accessed by researchers interested in using this dataset.
Our study provides a valuable resource for understanding the complex network of
interactions between TFs and their regulatory roles in cellular processes.


  Biomolecules composed of a limited set of chemical building blocks can
co-localize into distinct, spatially segregated compartments known as
biomolecular condensates. Although recent studies of intracellular condensates
have shown that coexisting, immiscible condensates can form spontaneously via
phase separation, it has remained unclear how coexisting and multiphase
condensates assemble from chemical building blocks with limited specificity.
Here we establish a connection between the interdependencies among biomolecular
interactions and the thermodynamic stability of multiphase condensates. We then
introduce an inverse design approach for computing the minimum interaction
specificity required to assemble condensates with prescribed molecular
compositions in a multicomponent biomolecular mixture. As a proof of principle,
we apply our theory to design mixtures of model heteropolymers using a minimal
number of distinct monomer types, and we use molecular simulations to verify
that our designs produce coexisting condensates with the target molecular
compositions. Our theoretical approach explains how multiphase condensates
arise in naturally occurring biomolecular mixtures and provides a rational
algorithm for engineering complex artificial condensates from simple chemical
building blocks.


  Phylogenetic network is an evolutionary model that uses a rooted directed
acyclic graph (instead of a tree) to model an evolutionary history of species
in which reticulate events (e.g., hybrid speciation or horizontal gene
transfer) occurred. Tree-child network is a kind of phylogenetic network with
structural constraints. Existing approaches for tree-child network
reconstruction can be slow for large data. In this paper, we present several
computational approaches for bounding from below the number of reticulations in
a tree-child network that displays a given set of rooted binary phylogenetic
trees. In addition, we also present some theoretical results on bounding from
above the number of reticulations. Through simulation, we demonstrate that the
new lower bounds on the reticulation number for tree-child networks can
practically be computed for large tree data. The bounds can provide estimates
of reticulation for relatively large data.


  The advent of highly accurate protein structure prediction methods has fueled
an exponential expansion of the protein structure database. Consequently, there
is a rising demand for rapid and precise structural homolog search. Traditional
alignment-based methods are dedicated to precise comparisons between pairs,
exhibiting high accuracy. However, their sluggish processing speed is no longer
adequate for managing the current massive volume of data. In response to this
challenge, we propose a novel deep-learning approach FoldExplorer. It harnesses
the powerful capabilities of graph attention neural networks and protein large
language models for protein structures and sequences data processing to
generate embeddings for protein structures. The structural embeddings can be
used for fast and accurate protein search. The embeddings also provide insights
into the protein space. FoldExplorer demonstrates a substantial performance
improvement of 5% to 8% over the current state-of-the-art algorithm on the
benchmark datasets. Meanwhile, FoldExplorer does not compromise on search speed
and excels particularly in searching on a large-scale dataset.


  AI algorithms proved excellent predictors of protein structure, but whether
their exceptional accuracy is merely due to megascale regression or these
algorithms learn the underlying physics remains an open question. Here, we
perform a stringent test for the existence of such learning in the Alphafold2
(AF) algorithm: We use AF to predict the subtle structural deformation induced
by single mutations, quantified by strain, and compare with experimental
datasets of corresponding perturbations in folding free energy $\Delta\Delta
G$. Unexpectedly, we find that physical strain alone -- without any additional
data or computation -- correlates almost as well with $\Delta \Delta G$ as
state-of-the-art energy-based and machine-learning predictors.This indicates
that the AF-predicted structures alone encode fine details about the energy
landscape. In particular, the structures encode significant information on
stability, enough to estimate (de-)stabilizing effects of mutations, thus
paving the way for the development of novel, structure-based stability
predictors for protein design and evolution.


  Machine learning prediction of organic materials properties is an efficient
virtual screening method ahead of more expensive screening methods. However,
this approach has suffered from insufficient labeled data on organic materials
to train state-of-the-art machine learning models. In this study, we
demonstrate that drug-like small molecule and chemical reaction databases can
be used to pretrain the BERT model for the virtual screening of organic
materials. Among the BERT models fine-tuned by five virtual screening tasks on
organic materials, the USPTO-SMILES pretrained BERT model had R2 > 0.90 for two
tasks and R2 > 0.82 for one, which was generally superior to the same models
pretrained by the small molecule or organic materials databases, as well as to
the other three traditional machine learning models trained directly on the
virtual screening task data. The superior performance of the USPTO-SMILES
pretrained BERT model is due to the greater variety of organic building blocks
in the USPTO database and the broader coverage of the chemical space. The even
better performance of the BERT model pretrained externally from a chemical
reaction database with additional sources of chemical reactions strengthens our
proof of concept that transfer learning across different chemical domains is
practical for the virtual screening of organic materials.


  Using asymmetric flow field flow fractionation, the polymer mass (Mw),
gyration radius (Rw) and the polydispersity index (PI) of glutenin polymers
(GPs) have been previously studied. Here, using the same multi-location trials
(4 years, 11 locations and 192 cultivars), the factors associated with the
conformation (Conf) of the polymers were reported. Conf, which is the slope of
Log (Rw) = f [Log (Mw)], varied between 0.285 and 0.740. Conf had low
broad-sense heritability (H 2 = 16.8) and was significantly influenced by
temperature occurring over the last month of grain filling. Higher temperatures
were found to increase Rw and the compactness and sphericity of GPs. The HMW-GS
and LMW-GS alleles had a significant influence on the Conf value. Assuming a
Gaussian distribution for Mw, the number of polymers present in wheat grain was
computed for different kernel weights and protein concentrations: it exceeded
10 12 GPs per grain. Using atomic force microscopy and cryo-microscopy, photos
of GPs were revealed for the first time. Under high temperatures, GPs became
larger, less reticulated, more spherical and consequently less prone to rapid
hydrolysis. Some orientations aimed at reducing the possible impact of the
numerous GPs for people suffering from non-celiac gluten sensitivity (NCGS)
were proposed.


  Molecular docking is a key computational tool utilized to predict the binding
conformations of small molecules to protein targets, which is fundamental in
the design of novel drugs. Despite recent advancements in geometric deep
learning-based approaches leading to improvements in blind docking efficiency,
these methods have encountered notable challenges, such as limited
generalization performance on unseen proteins, the inability to concurrently
address the settings of blind docking and site-specific docking, and the
frequent occurrence of physical implausibilities such as inter-molecular steric
clash. In this study, we introduce DeltaDock, a robust and versatile framework
designed for efficient molecular docking to overcome these challenges.
DeltaDock operates in a two-step process: rapid initial complex structures
sampling followed by multi-scale iterative refinement of the initial
structures. In the initial stage, to sample accurate structures with high
efficiency, we develop a ligand-dependent binding site prediction model founded
on large protein models and graph neural networks. This model is then paired
with GPU-accelerated sampling algorithms. The sampled structures are updated
using a multi-scale iterative refinement module that captures both
protein-ligand atom-atom interactions and residue-atom interactions in the
following stage. Distinct from previous geometric deep learning methods that
are conditioned on the blind docking setting, DeltaDock demonstrates superior
performance in both blind docking and site-specific docking settings.
Comprehensive experimental results reveal that DeltaDock consistently surpasses
baseline methods in terms of docking accuracy. Furthermore, it displays
remarkable generalization capabilities and proficiency for predicting
physically valid structures, thereby attesting to its robustness and
reliability in various scenarios.


  Spatio-temporal pathogen spread is often partially observed at the
metapopulation scale. Available data correspond to proxies and are incomplete,
censored and heterogeneous. Moreover, representing such biological systems
often leads to complex stochastic models. Such complexity together with data
characteristics make the analysis of these systems a challenge. Our objective
was to develop a new inference procedure to estimate key parameters of
stochastic metapopulation models of animal disease spread from longitudinal and
spatial datasets, while accurately accounting for characteristics of census
data. We applied our procedure to provide new knowledge on the regional spread
of \emph{Mycobacterium avium} subsp. \emph{paratuberculosis} (\emph{Map}),
which causes bovine paratuberculosis, a worldwide endemic disease. \emph{Map}
spread between herds through trade movements was modeled with a stochastic
mechanistic model. Comprehensive data from 2005 to 2013 on cattle movements in
12,857 dairy herds in Brittany (western France) and partial data on animal
infection status in 2,278 herds sampled from 2007 to 2013 were used. Inference
was performed using a new criterion based on a Monte-Carlo approximation of a
composite likelihood, coupled to a numerical optimization algorithm
(Nelder-Mead Simplex-like). Our criterion showed a clear superiority to
alternative ones in identifying the right parameter values, as assessed by an
empirical identifiability on simulated data. Point estimates and profile
likelihoods allowed us to establish the initial state of the system, identify
the risk of pathogen introduction from outside the metapopulation, and confirm
the assumption of the low sensitivity of the diagnostic test. Our inference
procedure could easily be applied to other spatio-temporal infection dynamics,
especially when ABC-like methods face challenges in defining relevant summary
statistics.


  Dispersal is a well recognized driver of ecological and evolutionary
dynamics, and simultaneously an evolving trait. Dispersal evolution has
traditionally been studied in single-species metapopulations so that it remains
unclear how dispersal evolves in spatially structured communities and food
webs. Since most natural systems are biodiverse and spatially structured, and
thus affected by dispersal and its evolution, this knowledge gap should be
bridged.
  Here we discuss whether knowledge established in single-species systems holds
in spatially structured multispecies systems and highlight generally valid and
fundamental principles. Most biotic interactions form the ecological theatre
for the evolutionary dispersal play because interactions mediate patterns of
fitness expectations in space and time. While this allows for a simple
transposition of certain known drivers to a multispecies context, other drivers
may require more complex transpositions, or might not be transferred. We
discuss an important quantitative modulator of dispersal evolution in the
increased trait dimensionality of biodiverse meta-systems and an additional
driver in co-dispersal.
  We speculate that scale and selection pressure mismatches due to
co-dispersal, together with increased trait dimensionality may lead to slower
and more "diffuse" evolution in biodiverse meta-systems. Open questions and
potential consequences in both ecological and evolutionary terms call for more
investigation.


  Modeling the interactions between drugs, targets, and diseases is paramount
in drug discovery and has significant implications for precision medicine and
personalized treatments. Current approaches frequently consider drug-target or
drug-disease interactions individually, ignoring the interdependencies among
all three entities. Within human metabolic systems, drugs interact with protein
targets in cells, influencing target activities and subsequently impacting
biological pathways to promote healthy functions and treat diseases. Moving
beyond binary relationships and exploring tighter triple relationships is
essential to understanding drugs' mechanism of action (MoAs). Moreover,
identifying the heterogeneity of drugs, targets, and diseases, along with their
distinct characteristics, is critical to model these complex interactions
appropriately. To address these challenges, we effectively model the
interconnectedness of all entities in a heterogeneous graph and develop a novel
Heterogeneous Graph Triplet Attention Network (\texttt{HeTriNet}).
\texttt{HeTriNet} introduces a novel triplet attention mechanism within this
heterogeneous graph structure. Beyond pairwise attention as the importance of
an entity for the other one, we define triplet attention to model the
importance of pairs for entities in the drug-target-disease triplet prediction
problem. Experimental results on real-world datasets show that
\texttt{HeTriNet} outperforms several baselines, demonstrating its remarkable
proficiency in uncovering novel drug-target-disease relationships.


  Deep learning promises to dramatically improve scoring functions for
molecular docking, leading to substantial advances in binding pose prediction
and virtual screening. To train scoring functions-and to perform molecular
docking-one must generate a set of candidate ligand binding poses.
Unfortunately, the sampling protocols currently used to generate candidate
poses frequently fail to produce any poses close to the correct, experimentally
determined pose, unless information about the correct pose is provided. This
limits the accuracy of learned scoring functions and molecular docking. Here,
we describe two improved protocols for pose sampling: GLOW (auGmented sampLing
with sOftened vdW potential) and a novel technique named IVES (IteratiVe
Ensemble Sampling). Our benchmarking results demonstrate the effectiveness of
our methods in improving the likelihood of sampling accurate poses, especially
for binding pockets whose shape changes substantially when different ligands
bind. This improvement is observed across both experimentally determined and
AlphaFold-generated protein structures. Additionally, we present datasets of
candidate ligand poses generated using our methods for each of around 5,000
protein-ligand cross-docking pairs, for training and testing scoring functions.
To benefit the research community, we provide these cross-docking datasets and
an open-source Python implementation of GLOW and IVES at
https://github.com/drorlab/GLOW_IVES .


  Accurate epidemic forecasting is a critical task in controlling disease
transmission. Many deep learning-based models focus only on static or dynamic
graphs when constructing spatial information, ignoring their relationship.
Additionally, these models often rely on recurrent structures, which can lead
to error accumulation and computational time consumption. To address the
aforementioned problems, we propose a novel model called Backbone-based Dynamic
Graph Spatio-Temporal Network (BDGSTN). Intuitively, the continuous and smooth
changes in graph structure, make adjacent graph structures share a basic
pattern. To capture this property, we use adaptive methods to generate static
backbone graphs containing the primary information and temporal models to
generate dynamic temporal graphs of epidemic data, fusing them to generate a
backbone-based dynamic graph. To overcome potential limitations associated with
recurrent structures, we introduce a linear model DLinear to handle temporal
dependencies and combine it with dynamic graph convolution for epidemic
forecasting. Extensive experiments on two datasets demonstrate that BDGSTN
outperforms baseline models and ablation comparison further verifies the
effectiveness of model components. Furthermore, we analyze and measure the
significance of backbone and temporal graphs by using information metrics from
different aspects. Finally, we compare model parameter volume and training time
to confirm the superior complexity and efficiency of BDGSTN.


  High-throughput drug screening -- using cell imaging or gene expression
measurements as readouts of drug effect -- is a critical tool in biotechnology
to assess and understand the relationship between the chemical structure and
biological activity of a drug. Since large-scale screens have to be divided
into multiple experiments, a key difficulty is dealing with batch effects,
which can introduce systematic errors and non-biological associations in the
data. We propose InfoCORE, an Information maximization approach for COnfounder
REmoval, to effectively deal with batch effects and obtain refined molecular
representations. InfoCORE establishes a variational lower bound on the
conditional mutual information of the latent representations given a batch
identifier. It adaptively reweighs samples to equalize their implied batch
distribution. Extensive experiments on drug screening data reveal InfoCORE's
superior performance in a multitude of tasks including molecular property
prediction and molecule-phenotype retrieval. Additionally, we show results for
how InfoCORE offers a versatile framework and resolves general distribution
shifts and issues of data fairness by minimizing correlation with spurious
features or removing sensitive attributes. The code is available at
https://github.com/uhlerlab/InfoCORE.


  We generalize the Bush--Mosteller learning, the Roth--Erev learning, and the
social learning to include mistakes such that the nonlinear replicator-mutator
equation with either additive or multiplicative mutation is generated in an
asymptotic limit. Subsequently, we exhaustively investigate the ubiquitous
Rock-Paper-Scissors game for some analytically tractable motifs of mutation
pattern. We consider both symmetric and asymmetric game interactions, and
reveal that mistakes can some-times help the players learn. While the
replicator-mutator flow exhibits rich dynamics that include limit cycles and
chaotic orbits, it can control chaos as well to lead to rational Nash
equilibrium outcome. Moreover, we also report an instance of hitherto unknown
Hamiltonian structure of the replicator-mutator equation.


  Determining the physicochemical properties of a protein can reveal important
insights in their structure, biological functions, stability, and interactions
with other molecules. Although tools for computing properties of proteins
already existed, we could not find a comprehensive tool that enables the
calculations of multiple properties for multiple input proteins on the proteome
level at once. Facing this limitation, we have developed Multiple Protein
Profiler (MPP) 1.0 as an integrated tool that allows the profiling of 12
individual properties of multiple proteins in a significant manner. MPP
provides a tabular and graphic visualization of properties of multiple
proteins. The tool is freely accessible at
https://mproteinprofiler.microbiologyandimmunology.dal.ca/


  The recent COVID-19 pandemic underscores the significance of early-stage
non-pharmacological intervention strategies. The widespread use of masks and
the systematic implementation of contact tracing strategies provide a
potentially equally effective and socially less impactful alternative to more
conventional approaches, such as large-scale mobility restrictions. However,
manual contact tracing faces strong limitations in accessing the network of
contacts, and the scalability of currently implemented protocols for
smartphone-based digital contact tracing becomes impractical during the rapid
expansion phases of the outbreaks, due to the surge in exposure notifications
and associated tests. A substantial improvement in digital contact tracing can
be obtained through the integration of probabilistic techniques for risk
assessment that can more effectively guide the allocation of new diagnostic
tests. In this study, we first quantitatively analyze the diagnostic and social
costs associated with these containment measures based on contact tracing,
employing three state-of-the-art models of SARS-CoV-2 spreading. Our results
suggest that probabilistic techniques allow for more effective mitigation at a
lower cost. Secondly, our findings reveal a remarkable efficacy of
probabilistic contact-tracing techniques in capturing backward propagations and
super-spreading events, relevant features of the diffusion of many pathogens,
including SARS-CoV-2.


  Notwithstanding its long history and compelling motivation, research seeking
to explicate the emergence life (EoL) has throughout been a cacophony of
unresolved speculation and dispute; absent still any clear convergence or other
inarguable evidence of progress. This notwithstanding that it has also produced
a rich and varied supply of putatively promising technical advances. Not
surprising then the effort being advanced by some to establish a shared basis
in fundamental assumptions upon which a more productive community research
effort might arise. In this essay however, I press a case in opposition. First,
that a chasm divides the rich fauna of contesting EoL models into two
conceptually incommensurate classes; here named "chemistry models" (to which
class belongs nearly all thinking and work in the field, past and present) and
"engine models" (advanced in various more-or-less partial forms by a marginal
minority of voices dating from Boltzmann forward). Second, that contemporary
non-equilibrium thermodynamics dictates that 'engine-less' (i.e. 'chemistry')
models cannot in principle generate non-equilibrium, organized states of matter
and are in consequence inherently incapable of prizing life out of inanimate
matter.


  Knowledge amount is an integral indicator of the development of society.
Humanity produces knowledge in response to challenges from nature and society.
Knowledge production depends on population size and human productivity.
Productivity is a function of knowledge amount. The purpose of this study is to
find this function and verify it on empirical material, including global
demographic and information data. The productivity function is a basic element
of the theory that results in the dynamic equations of knowledge production and
population growth. A separate problem is the quantitative assessment of
knowledge. To solve it, we consider knowledge representations in the form of
patents, articles and books. Knowledge is stored in various types of devices,
which together form a global informational storage. Storage capacity is
increasing rapidly as digital technology advances. We compare storage capacity
with the memory occupied by the forms of knowledge representation. The results
obtained in this study contribute to the theory of knowledge production and
related demographic dynamics and allow us to deepen our understanding of
civilization development.


  Predicting the structure of interacting chains is crucial for understanding
biological systems and developing new drugs. Large-scale pre-trained Protein
Language Models (PLMs), such as ESM2, have shown impressive abilities in
extracting biologically meaningful representations for protein structure
prediction. In this paper, we show that ESMFold, which has been successful in
computing accurate atomic structures for single-chain proteins, can be adapted
to predict the heterodimer structures in a lightweight manner. We propose
Linker-tuning, which learns a continuous prompt to connect the two chains in a
dimer before running it as a single sequence in ESMFold. Experiment results
show that our method successfully predicts 56.98% of interfaces on the i.i.d.
heterodimer test set, with an absolute improvement of +12.79% over the
ESMFold-Linker baseline. Furthermore, our model can generalize well to the
out-of-distribution (OOD) test set HeteroTest2 and two antibody test sets Fab
and Fv while being $9\times$ faster than AF-Multimer.


  Deep reinforcement learning has successfully been applied for molecular
discovery as shown by the Molecule Deep Q-network (MolDQN) algorithm. This
algorithm has challenges when applied to optimizing new molecules: training
such a model is limited in terms of scalability to larger datasets and the
trained model cannot be generalized to different molecules in the same dataset.
In this paper, a distributed reinforcement learning algorithm for antioxidants,
called DA-MolDQN is proposed to address these problems. State-of-the-art bond
dissociation energy (BDE) and ionization potential (IP) predictors are
integrated into DA-MolDQN, which are critical chemical properties while
optimizing antioxidants. Training time is reduced by algorithmic improvements
for molecular modifications. The algorithm is distributed, scalable for up to
512 molecules, and generalizes the model to a diverse set of molecules. The
proposed models are trained with a proprietary antioxidant dataset. The results
have been reproduced with both proprietary and public datasets. The proposed
molecules have been validated with DFT simulations and a subset of them
confirmed in public "unseen" datasets. In summary, DA-MolDQN is up to 100x
faster than previous algorithms and can discover new optimized molecules from
proprietary and public antioxidants.


  The escalating drug addiction crisis in the United States underscores the
urgent need for innovative therapeutic strategies. This study embarked on an
innovative and rigorous strategy to unearth potential drug repurposing
candidates for opioid and cocaine addiction treatment, bridging the gap between
transcriptomic data analysis and drug discovery. We initiated our approach by
conducting differential gene expression analysis on addiction-related
transcriptomic data to identify key genes. We propose a novel topological
differentiation to identify key genes from a protein-protein interaction (PPI)
network derived from DEGs. This method utilizes persistent Laplacians to
accurately single out pivotal nodes within the network, conducting this
analysis in a multiscale manner to ensure high reliability. Through rigorous
literature validation, pathway analysis, and data-availability scrutiny, we
identified three pivotal molecular targets, mTOR, mGluR5, and NMDAR, for drug
repurposing from DrugBank. We crafted machine learning models employing two
natural language processing (NLP)-based embeddings and a traditional 2D
fingerprint, which demonstrated robust predictive ability in gauging binding
affinities of DrugBank compounds to selected targets. Furthermore, we
elucidated the interactions of promising drugs with the targets and evaluated
their drug-likeness. This study delineates a multi-faceted and comprehensive
analytical framework, amalgamating bioinformatics, topological data analysis
and machine learning, for drug repurposing in addiction treatment, setting the
stage for subsequent experimental validation. The versatility of the methods we
developed allows for applications across a range of diseases and transcriptomic
datasets.


  In the domain of network biology, the interactions among heterogeneous
genomic and molecular entities are represented through networks. Link
prediction (LP) methodologies are instrumental in inferring missing or
prospective associations within these biological networks. In this review, we
systematically dissect the attributes of local, centrality, and embedding-based
LP approaches, applied to static and dynamic biological networks. We undertake
an examination of the current applications of LP metrics for predicting links
between diseases, genes, proteins, RNA, microbiomes, drugs, and neurons. We
carry out comprehensive performance evaluations on established biological
network datasets to show the practical applications of standard LP models.
Moreover, we compare the similarity in prediction trends among the models and
the specific network attributes that contribute to effective link prediction,
before underscoring the role of LP in addressing the formidable challenges
prevalent in biological systems, ranging from noise, bias, and data sparseness
to interpretability. We conclude the review with an exploration of the
essential characteristics expected from future LP models, poised to advance our
comprehension of the intricate interactions governing biological systems.


  During an epidemic, such as the COVID-19 pandemic, policy-makers are dealing
with the decision of effective, but socioeconomically costly interference, such
as horizontal lockdowns, including school and workplace closure, physical
distancing e.t.c. Investigating the role of the age of the individuals in the
evolution of epidemiological phenomena, we propose a scheme for the comparison
of certain epidemiological strategies. Then, we put the proposed scheme to the
test by employing an age-based epidemiological compartment model introduced in
a previous work of the authors, coupled with data from the literature, in order
to compare the effectiveness of age-based and horizontal interventions. In
general, our results suggest that these two are comparable mainly at low or
medium level of austerity.


  Longer timelines and lower success rates of drug candidates limit the
productivity of clinical trials in the pharmaceutical industry. Promising de
novo drug design techniques help solve this by exploring a broader chemical
space, efficiently generating new molecules, and providing improved therapies.
However, optimizing for molecular characteristics found in approved oral drugs
remains a challenge, limiting de novo usage. In this work, we propose NovoMol,
a novel de novo method using recurrent neural networks to mass-generate drug
molecules with high oral bioavailability, increasing clinical trial time
efficiency. Molecules were optimized for desirable traits and ranked using the
quantitative estimate of drug-likeness (QED). Generated molecules meeting QED's
oral bioavailability threshold were used to retrain the neural network, and,
after five training cycles, 76% of generated molecules passed this strict
threshold and 96% passed the traditionally used Lipinski's Rule of Five. The
trained model was then used to generate specific drug candidates for the
cancer-related PDGFR{\alpha} receptor and 44% of generated candidates had
better binding affinity than the current state-of-the-art drug, Imatinib (with
a receptor binding affinity of -9.4 kcal/mol), and the best-generated candidate
at -12.9 kcal/mol. NovoMol provides a time/cost-efficient AI-based de novo
method offering promising drug candidates for clinical trials.


  In a state-space framework, temporal variations in fishery-dependent
processes, such as selectivity and catchability, can be modeled as random
effects. This makes state-space models (SSMs) powerful tools for data-limited
assessments, especially when conventional CPUE standardization is inapplicable.
However, the flexibility of this modeling approach can lead to challenges such
as overfitting and parameter non-identifiability. To demonstrate and address
these challenges, we developed a state-space length-based age-structured model,
which we applied to the Korea chub mackerel ($\textit{Scomber japonicus}$)
stock as a case study. The model underwent rigorous scrutiny using various
model checking methods to detect potential model mis-specification and
non-identifiability under diverse scenarios. Our results demonstrated that
incorporating temporal variations in fishery-dependent processes through random
effects resolved model mis-specification, but excessive inclusion of random
effects rendered the model sensitive to a small number of observations, even
when the model was identifiable. For the non-identifiability issue, we employed
a non-degenerate estimator, using a gamma distribution as a penalty for the
standard deviation (SD) parameters of observation errors. This approach made
the SD parameters identifiable and facilitated the simultaneous estimation of
both process and observation error variances with minimal bias, known to be a
challenging task in SSMs. These findings underscore the importance of model
checking in SSMs and emphasize the need for careful consideration of
overfitting and non-identifiability when developing such models for
data-limited assessments. Additionally, novel assessment results for the
mackerel stock were presented, and implications for future stock assessment and
management were discussed.


  Biological cells utilize membranes and liquid-like droplets, known as
biomolecular condensates, to structure their interior. The interaction of
droplets and membranes, despite being involved in several key biological
processes, is so far little-understood. Here, we present a first numerical
method to simulate the continuum dynamics of droplets interacting with
deformable membranes via wetting. The method combines the advantages of the
phase field method for multi-phase flow simulation and the arbitrary
Lagrangian-Eulerian (ALE) method for an explicit description of the elastic
surface. The model is thermodynamically consistent, coupling bulk hydrodynamics
with capillary forces, as well as bending, tension, and stretching of a thin
membrane. The method is validated by comparing simulations for single droplets
to theoretical results of shape equations, and its capabilities are illustrated
in 2D and 3D axisymmetric scenarios.


  The integration of multi-omics data has emerged as a promising approach for
gaining comprehensive insights into complex diseases such as cancer. This paper
proposes a novel approach to identify cancer subtypes through the integration
of multi-omics data for clustering. The proposed method, named LIDAF utilises
affinity matrices based on linear relationships between and within different
omics datasets (Linear Inter and Intra Dataset Affinity Fusion (LIDAF)).
Canonical Correlation Analysis is in this paper employed to create distance
matrices based on Euclidean distances between canonical variates. The distance
matrices are converted to affinity matrices and those are fused in a three-step
process. The proposed LIDAF addresses the limitations of the existing method
resulting in improvement of clustering performance as measured by the Adjusted
Rand Index and the Normalized Mutual Information score. Moreover, our proposed
LIDAF approach demonstrates a notable enhancement in 50% of the log10 rank
p-values obtained from Cox survival analysis, surpassing the performance of the
best reported method, highlighting its potential of identifying distinct cancer
subtypes.


  Generating protein sequences that fold into a intended 3D structure is a
fundamental step in de novo protein design. De facto methods utilize
autoregressive generation, but this eschews higher order interactions that
could be exploited to improve inference speed. We describe a non-autoregressive
alternative that performs inference using a constant number of calls resulting
in a 23 times speed up without a loss in performance on the CATH benchmark.
Conditioned on the 3D structure, we fine-tune ProteinMPNN to perform discrete
diffusion with a purity prior over the index sampling order. Our approach gives
the flexibility in trading off inference speed and accuracy by modulating the
diffusion speed. Code: https://github.com/johnyang101/pmpnndiff


  Horizontal gene transfer inference approaches are usually based on gene
sequences: parametric methods search for patterns that deviate from a
particular genomic signature, while phylogenetic methods use sequences to
reconstruct the gene and species trees. However, it is well-known that
sequences have difficulty identifying ancient transfers since mutations have
enough time to erase all evidence of such events. In this work, we ask whether
character-based methods can predict gene transfers. Their advantage over
sequences is that homologous genes can have low DNA similarity, but still have
retained enough important common motifs that allow them to have common
character traits, for instance the same functional or expression profile. A
phylogeny that has two separate clades that acquired the same character
independently might indicate the presence of a transfer even in the absence of
sequence similarity. We introduce perfect transfer networks, which are
phylogenetic networks that can explain the character diversity of a set of taxa
under the assumption that characters have unique births, and that once a
character is gained it is rarely lost. Examples of such traits include
transposable elements, biochemical markers and emergence of organelles, just to
name a few. We study the differences between our model and two similar models:
perfect phylogenetic networks and ancestral recombination networks. Our goals
are to initiate a study on the structural and algorithmic properties of perfect
transfer networks. We then show that in polynomial time, one can decide whether
a given network is a valid explanation for a set of taxa, and show how, for a
given tree, one can add transfer edges to it so that it explains a set of taxa.
We finally provide lower and upper bounds on the number of transfers required
to explain a set of taxa, in the worst case.


  Advances in gene sequencing have enabled in silico analyses of microbial
genomes and have led to the revision of concepts of microbial taxonomy and
evolution. We explore deficiencies in existing multiple sequence global
alignment algorithms and introduce a new indexing scheme to partition the
dynamic programming algorithm hypercube scoring tensor over processors based on
the dependency between partitions to be scored in parallel. The performance of
algorithms is compared in the study of rpoB gene sequences of Mycoplasma
species.


  The spotted lanternfly (Lycorma delicatula) has recently spread from its
native range to several other countries and forecasts predict that it may
become a global invasive pest. In particular, since its confirmed presence in
the United States in 2014 it has established itself as a major invasive pest in
the Mid-Atlantic region where it is damaging both naturally occurring and
commercially important farmed plants. Quarantine zones have been introduced to
contain the infestation, but the spread to new areas continues. At present the
pathways and drivers of spread are not well-understood. In particular, several
human activity related factors have been proposed to contribute to the spread;
however, which features of the current spread can be attributed to these
factors remains unclear. Here we collect county level data on infestation
status and four human activity related factors and use statistical methods to
determine whether there is evidence for an association between the factors and
infestation. Then we construct a mechanistic network model based on the factors
found to be associated with infestation and use it to simulate local spread. We
find that the model reproduces key features of the spread 2014 to 2021. In
particular, the growth of the main infestation region and the opening of spread
corridors in the westward and southwestern directions is consistent with data
and the model accurately forecasts the correct infestation status at the county
level in 2021 with $81\%$ accuracy. We then use the model to forecast the
spread up to 2025 in a larger region. Given that this model is based on a few
human activity related factors that can be targeted, it may prove useful in
informing management and further modeling efforts related to the current
spotted lanternfly infestation in the US and potentially for current and future
invasions elsewhere globally.


  Recently, artificial intelligence for drug discovery has raised increasing
interest in both machine learning and chemistry domains. The fundamental
building block for drug discovery is molecule geometry and thus, the molecule's
geometrical representation is the main bottleneck to better utilize machine
learning techniques for drug discovery. In this work, we propose a pretraining
method for molecule joint auto-encoding (MoleculeJAE). MoleculeJAE can learn
both the 2D bond (topology) and 3D conformation (geometry) information, and a
diffusion process model is applied to mimic the augmented trajectories of such
two modalities, based on which, MoleculeJAE will learn the inherent chemical
structure in a self-supervised manner. Thus, the pretrained geometrical
representation in MoleculeJAE is expected to benefit downstream
geometry-related tasks. Empirically, MoleculeJAE proves its effectiveness by
reaching state-of-the-art performance on 15 out of 20 tasks by comparing it
with 12 competitive baselines.


  Online information ecosystems are now central to our everyday social
interactions. Of the many opportunities and challenges this presents, the
capacity for artificial agents to shape individual and collective human
decision-making in such environments is of particular importance. In order to
assess and manage the impact of artificial agents on human well-being, we must
consider not only the technical capabilities of such agents, but the impact
they have on human social dynamics at the individual and population level. We
approach this problem by modelling the potential for artificial agents to
"nudge" attitudes to fairness and cooperation in populations of human agents,
who update their behavior according to a process of social learning. We show
that the presence of artificial agents in a population playing the ultimatum
game generates highly divergent, multi-stable outcomes in the learning dynamics
of human agents' behaviour. These outcomes correspond to universal fairness
(successful nudging), universal selfishness (failed nudging), and a strategy of
fairness towards artificial agents and selfishness towards other human agents
(unintended consequences of nudging). We then consider the consequences of
human agents shifting their behavior when they are aware that they are
interacting with an artificial agent. We show that under a wide range of
circumstances artificial agents can achieve optimal outcomes in their
interactions with human agents while avoiding deception. However we also find
that, in the donation game, deception tends to make nudging easier to achieve.


  The identification of essential proteins in protein-protein interaction
networks (PINs) can help to discover drug targets and prevent disease. In order
to improve the accuracy of the identification of essential proteins,
researchers attempted to obtain a refined PIN by combining multiple biological
information to filter out some unreliable interactions in the PIN.
Unfortunately, such approaches drastically reduce the number of nodes in the
PIN after multiple refinements and result in a sparser PIN. It makes a
considerable portion of essential proteins unidentifiable. In this paper, we
propose a multi-layer refined network (MR-PIN) that addresses this problem.
Firstly, four refined networks are constructed by respectively integrating
different biological information into the static PIN to form a multi-layer
heterogeneous network. Then scores of proteins in each network layer are
calculated by the existing node ranking method, and the importance score of a
protein in the MR-PIN is evaluated in terms of the geometric mean of its scores
in all layers. Finally, all nodes are sorted by their importance scores to
determine their essentiality. To evaluate the effectiveness of the multi-layer
refined network model, we apply 16 node ranking methods on the MR-PIN, and
compare the results with those on the SPIN, DPIN and RDPIN. Then the predictive
performances of these ranking methods are validated in terms of the
identification number of essential protein at top100 - top600, sensitivity,
specificity, positive predictive value, negative predictive value, F-measure,
accuracy, Jackknife, ROCAUC and PRAUC. The experimental results show that the
MR-PIN is superior to the existing refined PINs in the identification accuracy
of essential proteins.


  The rapid expansion of genomic sequence data calls for new methods to achieve
robust sequence representations. Existing techniques often neglect intricate
structural details, emphasizing mainly contextual information. To address this,
we developed k-mer embeddings that merge contextual and structural string
information by enhancing De Bruijn graphs with structural similarity
connections. Subsequently, we crafted a self-supervised method based on
Contrastive Learning that employs a heterogeneous Graph Convolutional Network
encoder and constructs positive pairs based on node similarities. Our
embeddings consistently outperform prior techniques for Edit Distance
Approximation and Closest String Retrieval tasks.


  Predicting protein stability changes induced by single-point mutations has
been a persistent challenge over the years, attracting immense interest from
numerous researchers. The ability to precisely predict protein thermostability
is pivotal for various subfields and applications in biochemistry, including
drug development, protein evolution analysis, and enzyme synthesis. Despite the
proposition of multiple methodologies aimed at addressing this issue, few
approaches have successfully achieved optimal performance coupled with high
computational efficiency. Two principal hurdles contribute to the existing
challenges in this domain. The first is the complexity of extracting and
aggregating sufficiently representative features from proteins. The second
refers to the limited availability of experimental data for protein mutation
analysis, further complicating the comprehensive evaluation of model
performance on unseen data samples. With the advent of Large Language
Models(LLM), such as the ESM models in protein research, profound
interpretation of protein features is now accessibly aided by enormous training
data. Therefore, LLMs are indeed to facilitate a wide range of protein
research. In our study, we introduce an ESM-assisted efficient approach that
integrates protein sequence and structural features to predict the
thermostability changes in protein upon single-point mutations. Furthermore, we
have curated a dataset meticulously designed to preclude data leakage,
corresponding to two extensively employed test datasets, to facilitate a more
equitable model comparison.


  Molecular docking is critical to structure-based virtual screening, yet the
throughput of such workflows is limited by the expensive optimization of
scoring functions involved in most docking algorithms. We explore how machine
learning can accelerate this process by learning a scoring function with a
functional form that allows for more rapid optimization. Specifically, we
define the scoring function to be the cross-correlation of multi-channel ligand
and protein scalar fields parameterized by equivariant graph neural networks,
enabling rapid optimization over rigid-body degrees of freedom with fast
Fourier transforms. The runtime of our approach can be amortized at several
levels of abstraction, and is particularly favorable for virtual screening
settings with a common binding pocket. We benchmark our scoring functions on
two simplified docking-related tasks: decoy pose scoring and rigid conformer
docking. Our method attains similar but faster performance on crystal
structures compared to the widely-used Vina and Gnina scoring functions, and is
more robust on computationally predicted structures. Code is available at
https://github.com/bjing2016/scalar-fields.


  Cooperation in human society is sustained by reputation. In general, the
reputation of an individual is determined by others who observe his behavior,
but this rarely happens in private situations. This may cause people to behave
inconsistently, cooperating in public and not cooperating in private. A
previous experiment showed that people gave a lower reputation to an individual
who cooperated in public but defected in private rather than a consistently
uncooperative individual regardless of public and private situations. However,
the reason behind this is unclear. Here, we study how cooperation and the
reputational mechanism co-evolve on the condition that two types of interaction
(public and private) exist. The simulation results show that the evolved social
norm is characterized by at least one of the following: preference for
consistent or aversion of inconsistent behavior in both interactions when the
risk that behaviors in private interactions are observed exceeds a certain
threshold. We also find that such social norms promote cooperation in private
situations as well as in public ones.


  Fractional order compartment models derived from an underlying physical
stochastic process introduce a more nuanced and phenomenological approach to
system modelling than compartment models obtained by replacing integer order
derivatives with fractional order derivatives. In these models, inherent
stochasticity becomes important, particularly when dealing with the dynamics of
small populations; far from the continuum limit of large particle numbers. The
need for stochastic simulations is motivated by such deviations from mean
states of governing equations. This work presents an exact stochastic
simulation algorithm designed for fractional order compartment models, based on
a semi-Markov process. We have considered a fractional re-susceptibility
epidemic model as an illustrative example, highlighting pronounced disparities
between deterministic and stochastic dynamics when the total population is
small. Beyond this singular example, the algorithm presented is a general tool
for modelling fractional order differential equations via Monte Carlo
simulations.


  The success of therapeutic antibodies relies on their ability to selectively
bind antigens. AI-based antibody design protocols have shown promise in
generating epitope-specific designs. Many of these protocols use an inverse
folding step to generate diverse sequences given a backbone structure. Due to
prohibitive screening costs, it is key to identify candidate sequences likely
to bind in vitro. Here, we compare the efficacy of 8 common scoring paradigms
based on open-source models to classify antibody designs as binders or
non-binders. We evaluate these approaches on a novel surface plasmon resonance
(SPR) dataset, spanning 5 antigens. Our results show that existing methods
struggle to detect binders, and performance is highly variable across antigens.
We find that metrics computed on flexibly docked antibody-antigen complexes are
more robust, and ensembles scores are more consistent than individual metrics.
We provide experimental insight to analyze current scoring techniques,
highlighting that the development of robust, zero-shot filters is an important
research gap.


  Five distinct full-length mineralocorticoid receptor (MR) genes have been
identified in humans. These human MRs can be distinguished by the presence or
absence of an in-frame insertion of 12 base pairs coding for Lys, Cys, Ser, Trp
(KCSW) in their DNA-binding domain (DBD) and the presence of two amino acid
mutations in their amino terminal domain (NTD). Two human MRs with the KCSW
insertion (MR-KCSW) and three human MRs without KCSW in the DBD have been
identified. The three human MRs without KCSW contain either (Ile-180, Ala-241)
or (Val-180, Val-241) or (Ile-180, Val-241) in their NTD. The two human MRs
with KCSW contain either (Val-180, Val-241) or (Ile-180, Val-241) in their NTD.
Human MR-KCSW with (Ile-180, Ala-241) has not been cloned. In contrast,
chimpanzees contain two MRs with KCSW and two MRs without KCSW in their DBD and
both contain only Ile180, Val-241 in their NTDs. Each pair of chimpanzee MRs
differ at another amino acid in the NTD. A chimpanzee MR with either Val-180,
Val-241 or Ile-180, Ala-241 in the NTD has not been cloned. Gorillas and
orangutans each contain one MR with KCSW in the DBD and one MR without KCSW.
Both gorilla and orangutan MRs contain I-180, Val-241 in their NTD. Neither
Val-180, Val-241 nor Ile-180, Ala-241 are found in the NTD in either a gorilla
MR or an orangutan MR. These data suggest that human MRs with Val-180, Val-241
or Ile-180, Ala-241 in the NTD evolved after humans and chimpanzees diverged
from their common ancestor. These unique human MRs may have had a role in the
divergent evolution of humans from chimpanzees. Studies are underway to
characterize transcriptional activation of the five human MRs by aldosterone,
cortisol, and other corticosteroids for comparison with each other to elucidate
the roles of these MRs in human physiology.


  Epidemiological compartmental models are useful for understanding infectious
disease propagation and directing public health policy decisions. Calibration
of these models is an important step in offering accurate forecasts of disease
dynamics and the effectiveness of interventions. In this study, we present an
overview of calibrating strategies that can be employed, including several
optimization methods and reinforcement learning (RL). We discuss the benefits
and drawbacks of these methods and highlight relevant practical conclusions
from our experiments. Optimization methods iteratively adjust the parameters of
the model until the model output matches the available data, whereas RL uses
trial and error to learn the optimal set of parameters by maximizing a reward
signal. Finally, we discuss how the calibration of parameters of
epidemiological compartmental models is an emerging field that has the
potential to improve the accuracy of disease modeling and public health
decision-making. Further research is needed to validate the effectiveness and
scalability of these approaches in different epidemiological contexts. All
codes and resources are available on
\url{https://github.com/Nikunj-Gupta/On-the-Calibration-of-Compartmental-Epidemiological-Models}.
We hope this work can facilitate related research.


  Contact tracing is one of the most important tools for preventing the spread
of infectious diseases, but as the experience of COVID-19 showed, it is also
next-to-impossible to implement when the disease is spreading rapidly. We show
how to substantially improve the efficiency of contact tracing by combining
standard microeconomic tools that measure heterogeneity in how infectious a
sick person is with ideas from machine learning about sequential optimization.
Our contributions are twofold. First, we incorporate heterogeneity in
individual infectiousness in a multi-armed bandit to establish optimal
algorithms. At the heart of this strategy is a focus on learning. In the
typical conceptualization of contact tracing, contacts of an infected person
are tested to find more infections. Under a learning-first framework, however,
contacts of infected persons are tested to ascertain whether the infected
person is likely to be a "high infector" and to find additional infections only
if it is likely to be highly fruitful. Second, we demonstrate using three
administrative contact tracing datasets from India and Pakistan during COVID-19
that this strategy improves efficiency. Using our algorithm, we find 80% of
infections with just 40% of contacts while current approaches test twice as
many contacts to identify the same number of infections. We further show that a
simple strategy that can be easily implemented in the field performs at nearly
optimal levels, allowing for, what we call, feasible contact tracing. These
results are immediately transferable to contact tracing in any epidemic.


  The epidemiological behavior of Plasmodium vivax malaria occurs across
spatial scales including within-host, population, and metapopulation levels. On
the within-host scale, P. vivax sporozoites inoculated in a host may form
latent hypnozoites, the activation of which drives secondary infections and
accounts for a large proportion of P. vivax illness; on the metapopulation
level, the coupled human-vector dynamics characteristic of the population level
are further complicated by the migration of human populations across patches
with different malaria forces of (re-)infection. To explore the interplay of
all three scales in a single two-patch model of Plasmodium vivax dynamics, we
construct and study a system of eight integro-differential equations with
periodic forcing (arising from the single-frequency sinusoidal movement of a
human sub-population). Under the numerically-informed ansatz that the limiting
solutions to the system are closely bounded by sinusoidal ones for certain
regions of parameter space, we derive a single nonlinear equation from which
all approximate limiting solutions may be drawn, and devise necessary and
sufficient conditions for the equation to have only a disease-free solution.
Our results illustrate the impact of movement on P. vivax transmission and
suggest a need to focus vector control efforts on forest mosquito populations.
The three-scale model introduced here provides a more comprehensive framework
for studying the clinical, behavioral, and geographical factors underlying P.
vivax malaria endemicity.


  A central feature of vertebrate immune response is affinity maturation,
wherein antibody-producing B cells undergo evolutionary selection in
substructures of lymph nodes, called germinal centers. While it has been shown
that the median B cell affinity dependably increases over the course of
maturation, the exact logic behind this evolution remains vague. Three
potential selection methods include encouraging the reproduction of high
affinity cells (``birth/positive selection''), encouraging cell death in low
affinity cells (``death/negative selection''), and adjusting the mutation rate
based on cell affinity (``mutational selection''). While all three forms of
selection would lead to a net increase in affinity, different selection methods
may lead to distinct statistical dynamics. We present a tractable model of
selection, and analyze proposed signatures of negative selection. We find that
net selection strength and mutational distribution dynamics can distinguish
between birth and death selection. However, preferential ancestry ratios,
terminal node counts, and the mutational distribution shape are all ill-suited
for detecting selection method. We ultimately encourage thinking about
evolutionary fitness as more than a singular number.


  Mutated genes may lead to cancer development in numerous tissues. While more
than 600 cancer-causing genes are known today, some of the most widespread
mutations are connected to the RAS gene: RAS mutations are found in
approximately 25% of all human tumors. Specifically, KRAS mutations are
involved in the three most lethal cancers in U.S.: pancreatic ductal
adenocarcinoma, colorectal adenocarcinoma, and lung adenocarcinoma. These
cancers are among the most difficult to treat, and they are frequently excluded
from chemotherapeutic attacks as hopeless cases. The mutated KRAS proteins have
specific 3-dimensional conformations, which perturb functional interaction with
the GAP protein on the GAP:RAS complex surface leading to a signaling cascade
and uncontrolled cell growth. Here we describe a gluing docking method for
finding small molecules that bind to both the GAP and the mutated KRAS
molecules. These small molecules glue together the GAP and the mutated KRAS
molecules and may serve as new cancer drugs for the most lethal, most
difficult-to-treat carcinomas. As a proof of concept, we identify two new,
drug-like small molecules with the new method: these compounds specifically
inhibit the growth of PANC-1 cell line with KRAS mutation G12D in vitro and in
vivo. Importantly, the two new compounds show significantly lower IC50 and
higher specificity against the G12D KRAS mutant as compared to the recently
described MRTX-1133 inhibitor against the G12D KRAS mutant.


  Some proteins can find their targets on DNA faster than by pure diffusion in
the 3-dimensional cytoplasm, through the process of facilitated diffusion: They
can loosely bind to DNA and temporarily slide along it, thus being guided by
the DNA molecule itself to the target. This chapter examines this process in
mathematical detail with a focus on including the effect of DNA coiling on the
search process.


  In a developing embryo, information about the position of cells is encoded in
the concentrations of "morphogen" molecules. In the fruit fly, the local
concentrations of just a handful of proteins encoded by the gap genes are
sufficient to specify position with a precision comparable to the spacing
between cells along the anterior--posterior axis. This matches the precision of
downstream events such as the striped patterns of expression in the pair-rule
genes, but is not quite sufficient to define unique identities for individual
cells. We demonstrate theoretically that this information gap can be bridged if
positional errors are spatially correlated, with relatively long correlation
lengths. We then show experimentally that these correlations are present, with
the required strength, in the fluctuating positions of the pair-rule stripes,
and this can be traced back to the gap genes. Taking account of these
correlations, the available information matches the information needed for
unique cellular specification, within error bars of ~2%. These observation
support a precisionist view of information flow through the underlying genetic
networks, in which accurate signals are available from the start and preserved
as they are transformed into the final spatial patterns.


  Exponential-family random graph models (ERGMs) are a family of network models
originating in social network analysis, which have also been applied to
biological networks. Advances in estimation algorithms have increased the
practical scope of these models to larger networks, however it is still not
always possible to estimate a model without encountering problems of model
near-degeneracy, particularly if it is desired to use only simple model
parameters, rather than more complex parameters designed to overcome the
problem of near-degeneracy. Two new network models related to the ERGM, the
Tapered ERGM, and the latent order logistic (LOLOG) model, have recently been
proposed to overcome this problem. In this work I illustrate the application of
the Tapered ERGM and the LOLOG to a set of biological networks, including
protein-protein interaction (PPI) networks, gene regulatory networks, and
neural networks. I find that the Tapered ERGM and the LOLOG are able to
estimate models for networks for which it was not possible to estimate a
conventional ERGM, and are able to do so using only simple model parameters. In
the case of two neural networks where data on the spatial position of neurons
is available, this allows the estimation of models including terms for spatial
distance and triangle structures, allowing triangle motif statistical
significance to be estimated while accounting for the effect of spatial
proximity on connection probability. For some larger networks, however, Tapered
ERGM and LOLOG estimation was not possible in practical time, while
conventional ERGM models were able to be estimated only by using the
Equilibrium Expectation (EE) algorithm.


  The diversity and heterogeneity of biomarkers has made the development of
general methods for single-step quantification of analytes difficult. For
individual biomarkers, electrochemical methods that detect a conformational
change in an affinity binder upon analyte binding have shown promise. However,
because the conformational change must operate within a nanometer-scale working
distance, an entirely new sensor, with a unique conformational change, must be
developed for each analyte. Here, we demonstrate a modular electrochemical
biosensor, built from DNA origami, which is easily adapted to diverse molecules
by merely replacing its analyte binding domains. Instead of relying on a unique
nanometer-scale movement of a single redox reporter, all sensor variants rely
on the same 100-nanometer scale conformational change, which brings dozens of
reporters close enough to a gold electrode surface that a signal can be
measured via square wave voltammetry, a standard electrochemical technique. To
validate our sensor's mechanism, we used single-stranded DNA as an analyte, and
optimized the number of redox reporters and various linker lengths. Adaptation
of the sensor to streptavidin and PDGF-BB analytes was achieved by simply
adding biotin or anti-PDGF aptamers to appropriate DNA linkers.
Geometrically-optimized streptavidin sensors exhibited signal gain and limit of
detection markedly better than comparable reagentless electrochemical sensors.
After use, the same sensors could be regenerated under mild conditions:
performance was largely maintained over four cycles of DNA strand displacement
and rehybridization. By leveraging the modularity of DNA nanostructures, our
work provides a straightforward route to the single-step quantification of
arbitrary nucleic acids and proteins.