Exploring and mapping chemical space with molecular assembly trees

The Chemical Space Spanned by Manually Curated Datasets of Natural and Synthetic Compounds with Activities against SARS-CoV-2

Diseases caused by viruses are challenging to contain, as their outbreak and spread could be very sudden, compounded by rapid mutations, making the development of drugs and vaccines a continued endeavour that requires fast discovery and preparedness. Targeting viral infections with small molecules remains one of the treatment options to reduce transmission and the disease burden. A lesson learned from the recent coronavirus disease (COVID-19) is to collect ready-to-screen small molecule libraries in preparation for the next viral outbreak, and potentially find a clinical candidate before it becomes a pandemic. Public availability of diverse compound libraries, well annotated in terms of chemical structures and scaffolds, modes of action, and bioactivities are, therefore, crucial to ensure the participation of academic laboratories in these screening efforts, especially in resource-limited settings where synthesis, testing and computing capacity are scarce. Here, we demonstrate a low-resource approach to populate the chemical space of naturally occurring and synthetic small molecules that have shown in vitro and/or in vivo activities against the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its target proteins. We have manually curated two datasets of small molecules (naturally occurring and synthetically derived) by reading and collecting (hand-curating) the published literature. Information from the literature reveals that a majority of the reported SARS-CoV-2 compounds act by inhibiting the main protease, while 25% of the compounds currently have no known target. Scaffold analysis and principal component analysis revealed that the most common scaffolds in the datasets are quite distinct. We then expanded the initially manually curated dataset of over 1200 compounds via an ultra-large scale 2D and 3D similarity search, obtaining an expanded collection of over 150 k purchasable compounds. The spanned chemical space significantly extends beyond that of a commercially available coronavirus library of more than 20 k small molecules and constitutes a good starting collection for virtual screening campaigns given its manageable size and proximity to hand-curated compounds.

"Assembly Theory" in life-origin models: A critical review

Any homeostatic protometabolism would have required orchestration of disparate biochemical pathways into integrated circuits. Extraordinarily specific molecular assemblies were also required at the right time and place. Assembly Theory conflated with its cousins-Complexity Theory, Chaos theory, Quantum Mechanics, Irreversible Nonequilibrium Thermodynamics and Molecular Evolution theory- collectively have great naturalistic appeal in hopes of their providing the needed exquisite steering and controls. They collectively offer the best hope of circumventing the need for active selection required to formally orchestrate bona fide formal organization (as opposed to the mere self-ordering of chaos theory) (Abel and Trevors, 2006b). This paper focuses specifically on AT's contribution to naturalistic life-origin models.

Discovery of Highly Bioactive Peptides through Hierarchical Structural Information and Molecular Dynamics Simulations

Peptide drugs play an essential role in modern therapeutics, but the computational design of these molecules is hindered by several challenges. Traditional methods like molecular docking and molecular dynamics (MD) simulation, as well as recent deep learning approaches, often face limitations related to computational resource demands, complex binding affinity assessments, extensive data requirements, and poor model interpretability. Here, we introduce PepHiRe, an innovative methodology that utilizes the hierarchical structural information in peptide sequences and employs a novel strategy called Ladderpath, rooted in algorithmic information theory, to rapidly generate and enhance the efficiency and clarity of novel peptide design. We applied PepHiRe to develop BH3-like peptide inhibitors targeting myeloid cell leukemia-1, a protein associated with various cancers. By analyzing just eight known bioactive BH3 peptide sequences, PepHiRe effectively derived a hierarchy of subsequences used to create new BH3-like peptides. These peptides underwent screening through MD simulations, leading to the selection of five candidates for synthesis and subsequent in vitro testing. Experimental results demonstrated that these five peptides possess high inhibitory activity, with IC50 values ranging from 28.13 ± 7.93 to 167.42 ± 22.15 nM. Our study explores a white-box model driven technique and a structured screening pipeline for identifying and generating novel peptides with potential bioactivity.

Experimentally measured assembly indices are required to determine the threshold for life

Assembly theory (AT) aims to distinguish living from non-living systems by explaining and quantifying selection and evolution. The theory proposes that the degree of assembly depends on the number of complex objects, with complexity measured using a combination of the object's assembly index (AI) and its abundance. We previously provided experimental evidence supporting AT's predictive power, finding that abiotic systems do not randomly produce organic molecules with an AI greater than approximately 15 in detectable amounts. Hazen et al. (Hazen et al. 2024 J. R. Soc. Interface 21, 20230632. (doi:10.1098/rsif.2023.0632)) proposed inorganic molecules that theoretically have AIs greater than 15, suggesting similar complexity to biological molecules. However, our AIs are experimentally measured for organic, covalently bonded molecules, whereas Hazen's are theoretical, derived from crystal structures of charged units that are not isolable in solution. This distinction underscores the challenge in experimentally validating theoretical AIs.

Reply to 'Experimental measurement of assembly indices are required to determine the threshold for life'

We clarify misunderstandings of Walker et al. (Walker et al. 2024 J. R. Soc. Interface 21, 20240367 (doi:10.1098/rsif.2024.0367)) related to studies of the assembly pathways of molecular subunits in minerals. The finding that these subunits have calculated assembly pathways less than approximately 25 informs a central premise of Assembly Theory-that only life can produce numerous copies of molecules with assembly indices above a threshold value. What that threshold value might be, and whether the same value applies to chemical systems as different as organic and inorganic molecules, are questions deserving of additional study.

Autocatalytic Sets and Assembly Theory: A Toy Model Perspective

Assembly Theory provides a promising framework to explain the complexity of systems such as molecular structures and the origins of life, with broad applicability across various disciplines. In this study, we explore and consolidate different aspects of Assembly Theory by introducing a simplified Toy Model to simulate the autocatalytic formation of complex structures. This model abstracts the molecular formation process, focusing on the probabilistic control of catalysis rather than the intricate interactions found in organic chemistry. We establish a connection between probabilistic catalysis events and key principles of Assembly Theory, particularly the probability of a possible construction path in the formation of a complex object, and examine how the assembly of complex objects is impacted by the presence of autocatalysis. Our findings suggest that this presence of autocatalysis tends to favor longer consecutive construction sequences in environments with a low probability of catalysis, while this bias diminishes in environments with higher catalysis probabilities, highlighting the significant influence of environmental factors on the assembly of complex structures.

Self-organization in computation and chemistry: Return to AlChemy

How do complex adaptive systems, such as life, emerge from simple constituent parts? In the 1990s, Walter Fontana and Leo Buss proposed a novel modeling approach to this question, based on a formal model of computation known as the λ calculus. The model demonstrated how simple rules, embedded in a combinatorially large space of possibilities, could yield complex, dynamically stable organizations, reminiscent of biochemical reaction networks. Here, we revisit this classic model, called AlChemy, which has been understudied over the past 30 years. We reproduce the original results and study the robustness of those results using the greater computing resources available today. Our analysis reveals several unanticipated features of the system, demonstrating a surprising mix of dynamical robustness and fragility. Specifically, we find that complex, stable organizations emerge more frequently than previously expected, that these organizations are robust against collapse into trivial fixed points, but that these stable organizations cannot be easily combined into higher order entities. We also study the role played by the random generators used in the model, characterizing the initial distribution of objects produced by two random expression generators, and their consequences on the results. Finally, we provide a constructive proof that shows how an extension of the model, based on the typed λ calculus, could simulate transitions between arbitrary states in any possible chemical reaction network, thus indicating a concrete connection between AlChemy and chemical reaction networks. We conclude with a discussion of possible applications of AlChemy to self-organization in modern programming languages and quantitative approaches to the origin of life.

Single-Molecule Mixture: A Concept in Polymer Science

In theory, two extreme forms of substances exist: the pure form and the single-molecule mixture form. The latter contains a mixture of molecules with molecularly different structures. Inspired by the "chemical space" concept, in this paper, I report a study of the single-molecule mixture state that combines model construction and mathematical analysis, obtaining some interesting results. These results provide theoretical evidence that the single-molecule mixture state may indeed exist in realistic synthetic or natural polymer systems.

Scientific Discovery Framework Accelerating Advanced Polymeric Materials Design

Organic polymer materials, as the most abundantly produced materials, possess a flammable nature, making them potential hazards to human casualties and property losses. Target polymer design is still hindered due to the lack of a scientific foundation. Herein, we present a robust, generalizable, yet intelligent polymer discovery framework, which synergizes diverse capabilities, including the in situ burning analyzer, virtual reaction generator, and material genomic model, to achieve results that surpass the sum of individual parts. Notably, the high-throughput analyzer created for the first time, grounded in multiple spectroscopic principles, enables in situ capturing of massive combustion intermediates; then, the created realistic apparatus transforming to the virtual reaction generator acquires exponentially more intermediate information; further, the proposed feature engineering tool, which embedded both polymer hierarchical structures and massive intermediate data, develops the generalizable genomic model with excellent universality (adapting over 20 kinds of polymers) and high accuracy (88.8%), succeeding discovering series of novel polymers. This emerging approach addresses the target polymer design for flame-retardant application and underscores a pivotal role in accelerating polymeric materials discovery.

Pan-Evo: The Evolution of Information and Biology's Part in This

Many people wonder whether biology, including humans, will benefit or experience harm from new developments in information such as artificial intelligence (AI). Here, it is proposed that biological and non-biological information might be components of a unified process, 'Panevolution' or 'Pan-Evo', based on four basic operations-innovation, transmission, adaptation, and movement. Pan-Evo contains many types of variable objects, from molecules to ecosystems. Biological innovation includes mutations and behavioural changes; non-biological innovation includes naturally occurring physical innovations and innovation in software. Replication is commonplace in and outside biology, including autocatalytic chemicals and autonomous software replication. Adaptation includes biological selection, autocatalytic chemicals, and 'evolutionary programming', which is used in AI. The extension of biological speciation to non-biological information creates a concept called 'Panspeciation'. Panevolution might benefit or harm biology, but the harm might be minimal if AI and humans behave intelligently because humans and the machines in which an AI resides might split into vastly different environments that suit them. That is a possible example of Panspeciation and would be the first speciation event involving humans for thousands of years. This event will not be particularly hostile to humans if humans learn to evaluate information and cooperate better to minimise both human stupidity and artificial simulated stupidity (ASS-a failure of AI).

In silico discovery of potential PPI inhibitors for anti-lung cancer activity by targeting the CCND1-CDK4 complex via the P21 inhibition mechanism

Non-Small Cell Lung Cancer (NSCLC) is a prevalent and deadly form of lung cancer worldwide with a low 5-year survival rate. Current treatments have limitations, particularly for advanced-stage patients. P21, a protein that inhibits the CCND1-CDK4 complex, plays a crucial role in cell proliferation. Computer-Aided Drug Design (CADD) based on pharmacophores can screen and design PPI inhibitors targeting the CCND1-CDK4 complex. By analyzing known inhibitors, key pharmacophores are identified, and computational methods are used to screen potential PPI inhibitors. Molecular docking, pharmacophore matching, and structure-activity relationship studies optimize the inhibitors. This approach accelerates the discovery of CCND1-CDK4 PPI inhibitors for NSCLC treatment. Molecular dynamics simulations of CCND1-CDK4-P21 and CCND1-CDK4 complexes showed stable behavior, comprehensive sampling, and P21's impact on complex stability and hydrogen bond formation. A pharmacophore model facilitated virtual screening, identifying compounds with favorable binding affinities. Further simulations confirmed the stability and interactions of selected compounds, including 513457. This study demonstrates the potential of CADD in optimizing PPI inhibitors targeting the CCND1-CDK4 complex for NSCLC treatment. Extended simulations and experimental validations are necessary to assess their efficacy and safety.

Assembly Theory: What It Does and What It Does Not Do

A recent publication in Nature has generated much heated discussion about evolution, its tendency towards increasing diversity and complexity, and its potential status above and beyond the known laws of fundamental physics. The argument at the heart of this controversy concerns assembly theory, a method to detect and quantify the influence of higher-level emergent causal constraints in computational worlds made of basic objects and their combinations. In this short essay, I briefly review the theory, its basic principles and potential applications. I then go on to critically examine its authors' assertions, concluding that assembly theory has merit but is not nearly as novel or revolutionary as claimed. It certainly does not provide any new explanation of biological evolution or natural selection, or a new grounding of biology in physics. In this regard, the presentation of the paper is starkly distorted by hype, which may explain some of the outrage it created.

Molecular assembly indices of mineral heteropolyanions: some abiotic molecules are as complex as large biomolecules

Molecular assembly indices, which measure the number of unique sequential steps theoretically required to construct a three-dimensional molecule from its constituent atomic bonds, have been proposed as potential biosignatures. A central hypothesis of assembly theory is that any molecule with an assembly index ≥15 found in significant local concentrations represents an unambiguous sign of life. We show that abiotic molecule-like heteropolyanions, which assemble in aqueous solution as precursors to some mineral crystals, range in molecular assembly indices from 2 for H2CO3 or Si(OH)4 groups to as large as 21 for the most complex known molecule-like subunits in the rare minerals ewingite and ilmajokite. Therefore, values of molecular assembly indices ≥15 do not represent unambiguous biosignatures.

False Positives and the Challenge of Testing the Alien Hypothesis

The origin of life and the detection of alien life have historically been treated as separate scientific research problems. However, they are not strictly independent. Here, we discuss the need for a better integration of the sciences of life detection and origins of life. Framing these dual problems within the formalism of Bayesian hypothesis testing, we demonstrate via simple examples how high confidence in life detection claims require either (1) a strong prior hypothesis about the existence of life in a particular alien environment, or conversely, (2) signatures of life that are not susceptible to false positives. As a case study, we discuss the role of priors and hypothesis testing in recent results reporting potential detection of life in the venusian atmosphere and in the icy plumes of Enceladus. While many current leading biosignature candidates are subject to false positives because they are not definitive of life, our analyses demonstrate why it is necessary to shift focus to candidate signatures that are definitive. This indicates a necessity to develop methods that lack substantial false positives, by using observables for life that rely on prior hypotheses with strong theoretical and empirical support in identifying defining features of life. Abstract theories developed in pursuit of understanding universal features of life are more likely to be definitive and to apply to life-as-we-don't-know-it. We discuss Molecular Assembly theory as an example of such an observable which is applicable to life detection within the solar system. In the absence of alien examples these are best validated in origin of life experiments, substantiating the need for better integration between origins of life and biosignature science research communities. This leads to a conclusion that extraordinary claims in astrobiology (e.g., definitive detection of alien life) require extraordinary explanations, whereas the evidence itself could be quite ordinary.

Assembly theory explains and quantifies selection and evolution

Scientists have grappled with reconciling biological evolution1,2 with the immutable laws of the Universe defined by physics. These laws underpin life's origin, evolution and the development of human culture and technology, yet they do not predict the emergence of these phenomena. Evolutionary theory explains why some things exist and others do not through the lens of selection. To comprehend how diverse, open-ended forms can emerge from physics without an inherent design blueprint, a new approach to understanding and quantifying selection is necessary3-5. We present assembly theory (AT) as a framework that does not alter the laws of physics, but redefines the concept of an 'object' on which these laws act. AT conceptualizes objects not as point particles, but as entities defined by their possible formation histories. This allows objects to show evidence of selection, within well-defined boundaries of individuals or selected units. We introduce a measure called assembly (A), capturing the degree of causation required to produce a given ensemble of objects. This approach enables us to incorporate novelty generation and selection into the physics of complex objects. It explains how these objects can be characterized through a forward dynamical process considering their assembly. By reimagining the concept of matter within assembly spaces, AT provides a powerful interface between physics and biology. It discloses a new aspect of physics emerging at the chemical scale, whereby history and causal contingency influence what exists.

Preparation of 3,5-Methanobenzo[b]azepines: An sp3-Rich Quinolone Isostere

The replacement of the aromatic ring in bioactive compounds with saturated bioisosteres has become a popular tactic to obtain novel structures with improved physicochemical profiles. In this paper, we describe an efficient synthesis of 3,5-methanobenzo[b]azepine analogues and suggest them as isosteres of quinolones. Quinolones are heteroaromatic, flat rings and considered as privileged scaffolds. An isosteric version of this scaffold with more 3D character would offer new options to expand their use.

Hippocampal and Medial Prefrontal Cortex Fractal Spiking Patterns Encode Episodes and Rules

A central question in neuroscience is how the brain represents and processes information to guide behavior. The principles that organize brain computations are not fully known, and could include scale-free, or fractal patterns of neuronal activity. Scale-free brain activity may be a natural consequence of the relatively small subsets of neuronal populations that respond to task features, i.e., sparse coding. The size of the active subsets constrains the possible sequences of inter-spike intervals (ISI), and selecting from this limited set may produce firing patterns across wide-ranging timescales that form fractal spiking patterns. To investigate the extent to which fractal spiking patterns corresponded with task features, we analyzed ISIs in simultaneously recorded populations of CA1 and medial prefrontal cortical (mPFC) neurons in rats performing a spatial memory task that required both structures. CA1 and mPFC ISI sequences formed fractal patterns that predicted memory performance. CA1 pattern duration, but not length or content, varied with learning speed and memory performance whereas mPFC patterns did not. The most common CA1 and mPFC patterns corresponded with each region's cognitive function: CA1 patterns encoded behavioral episodes which linked the start, choice, and goal of paths through the maze whereas mPFC patterns encoded behavioral "rules" which guided goal selection. mPFC patterns predicted changing CA1 spike patterns only as animals learned new rules. Together, the results suggest that CA1 and mPFC population activity may predict choice outcomes by using fractal ISI patterns to compute task features.

Dynamic Environmental Conditions Affect the Composition of a Model Prebiotic Reaction Network

Prebiotic environments are dynamic, containing a range of periodic and aperiodic variations in reaction conditions. However, the impact of the temporal dynamics of environmental conditions upon prebiotic chemical reaction networks has not been investigated. Here, we demonstrate how the magnitude and rate of temporal fluctuations of the catalysts Ca²⁺ and hydroxide control the product distributions of the formose reaction. Surprisingly, the product compositions of the formose reaction under dynamic conditions deviate significantly from those under steady state conditions. We attribute these compositional changes to the non-uniform propagation of fluctuations through the network, thereby shaping reaction outcomes. An examination of temporal concentration patterns showed that collections of compounds responded collectively to perturbations, indicating that key gating reactions branching from the Breslow cycle may be important responsive features of the formose reaction. Our findings show how the compositions of prebiotic reaction networks were shaped by sequential environmental events, illustrating the necessity for considering the temporal traits of prebiotic environments that supported the origin of life.

Ladderpath Approach: How Tinkering and Reuse Increase Complexity and Information

The notion of information and complexity are important concepts in many scientific fields such as molecular biology, evolutionary theory and exobiology. Many measures of these quantities are either difficult to compute, rely on the statistical notion of information, or can only be applied to strings. Based on assembly theory, we propose the notion of a ladderpath, which describes how an object can be decomposed into hierarchical structures using repetitive elements. From the ladderpath, two measures naturally emerge: the ladderpath-index and the order-index, which represent two axes of complexity. We show how the ladderpath approach can be applied to both strings and spatial patterns and argue that all systems that undergo evolution can be described as ladderpaths. Further, we discuss possible applications to human language and the origin of life. The ladderpath approach provides an alternative characterization of the information that is contained in a single object (or a system) and could aid in our understanding of evolving systems and the origin of life in particular.

Automated Exploration of Prebiotic Chemical Reaction Space: Progress and Perspectives

Prebiotic chemistry often involves the study of complex systems of chemical reactions that form large networks with a large number of diverse species. Such complex systems may have given rise to emergent phenomena that ultimately led to the origin of life on Earth. The environmental conditions and processes involved in this emergence may not be fully recapitulable, making it difficult for experimentalists to study prebiotic systems in laboratory simulations. Computational chemistry offers efficient ways to study such chemical systems and identify the ones most likely to display complex properties associated with life. Here, we review tools and techniques for modelling prebiotic chemical reaction networks and outline possible ways to identify self-replicating features that are central to many origin-of-life models.