Search

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

Metabolic heterogeneity is widely recognised as the next challenge in our understanding of non-genetic variation. A growing body of evidence suggests that metabolic heterogeneity may result from the inherent stochasticity of intracellular events. However, metabolism has been traditionally viewed as a purely deterministic process, on the basis that highly abundant metabolites tend to filter out stochastic phenomena. Here we bridge this gap with a general method for prediction of metabolite distributions across single cells. By exploiting the separation of time scales between enzyme expression and enzyme kinetics, our method produces estimates for metabolite distributions without the lengthy stochastic simulations that would be typically required for large metabolic models. The metabolite distributions take the form of Gaussian mixture models that are directly computable from single-cell expression data and standard deterministic models for metabolic pathways. The proposed mixture models provide a systematic method to predict the impact of biochemical parameters on metabolite distributions. Our method lays the groundwork for identifying the molecular processes that shape metabolic heterogeneity and its functional implications in disease., Comment: 5 Figures, 3 Tables

Metabolism plays a central role in cell physiology because it provides the molecular machinery for growth. At the genome-scale, metabolism is made up of thousands of reactions interacting with one another. Untangling this complexity is key to understand how cells respond to genetic, environmental, or therapeutic perturbations. Here we discuss the roles of two complementary strategies for the analysis of genome-scale metabolic models: Flux Balance Analysis (FBA) and network science. While FBA estimates metabolic flux on the basis of an optimisation principle, network approaches reveal emergent properties of the global metabolic connectivity. We highlight how the integration of both approaches promises to deliver insights on the structure and function of metabolic systems with wide-ranging implications in discovery science, precision medicine and industrial biotechnology.

Heterologous gene expression draws resources from host cells. These resources include vital components to sustain growth and replication, and the resulting cellular burden is a widely recognised bottleneck in the design of robust circuits. In this tutorial we discuss the use of computational models that integrate gene circuits and the physiology of host cells. Through various use cases, we illustrate the power of host-circuit models to predict the impact of design parameters on both burden and circuit functionality. Our approach relies on a new generation of computational models for microbial growth that can flexibly accommodate resource bottlenecks encountered in gene circuit design. Adoption of this modelling paradigm can facilitate fast and robust design cycles in synthetic biology.

Phenotypic variation is a hallmark of cellular physiology. Metabolic heterogeneity, in particular, underpins single-cell phenomena such as microbial drug tolerance and growth variability. Much research has focussed on transcriptomic and proteomic heterogeneity, yet it remains unclear if such variation permeates to the metabolic state of a cell. Here we propose a stochastic model to show that complex forms of metabolic heterogeneity emerge from fluctuations in enzyme expression and catalysis. The analysis predicts clonal populations to split into two or more metabolically distinct subpopulations. We reveal mechanisms not seen in deterministic models, in which enzymes with unimodal expression distributions lead to metabolites with a bimodal or multimodal distribution across the population. Based on published data, the results suggest that metabolite heterogeneity may be more pervasive than previously thought. Our work casts light on links between gene expression and metabolism, and provides a theory to probe the sources of metabolite heterogeneity.

We present a simple, top-down approach for the calculation of minimum energy consumption of electrosorptive ion separation using variational form of the (Gibbs) free energy. We focus and expand on the case of electrostatic capacitive deionization (CDI), and the theoretical framework is independent of details of the double-layer charge distribution and is applicable to any thermodynamically consistent model, such as the Gouy-Chapman-Stern (GCS) and modified Donnan (mD) models. We demonstrate that, under certain assumptions, the minimum required electric work energy is indeed equivalent to the free energy of separation. Using the theory, we define the thermodynamic efficiency of CDI. We explore the thermodynamic efficiency of current experimental CDI systems and show that these are currently very low, less than 1% for most existing systems. We applied this knowledge and constructed and operated a CDI cell to show that judicious selection of the materials, geometry, and process parameters can be used to achieve a 9% thermodynamic efficiency (4.6 kT energy per removed ion). This relatively high value is, to our knowledge, by far the highest thermodynamic efficiency ever demonstrated for CDI. We hypothesize that efficiency can be further improved by further reduction of CDI cell series resistances and optimization of operational parameters.

Nitrate is an important pollutant in drinking water worldwide, and a number of methods exist for the removal of nitrate from water including ion exchange and reverse osmosis. However, these approaches suffer from a variety of disadvantages including requirements for supply and disposal of brine used for regeneration in ion exchange and low water recovery ratio for reverse osmosis. Here, we demonstrate the use of high surface area activated carbon electrodes functionalized with moieties having high affinity for adsorption of nitrate from aqueous solution, such as those used in ion exchange. Adsorption of surfactant molecules having a quaternary amine ionic group to the activated carbon surfaces provides functionalization of the surfaces without complex chemistries. The functionalized electrodes have adsorption capacities of about 80 mg NaNO$_3$ per gram of activated carbon material. Unlike a traditional ion exchanger, the functionalized surfaces can be repeatedly regenerated by the application of an electrostatic potential which displaces the bound NO$_3^-$ while leaving an excess of electronic charge on the electrode. The cell is completed by an inert counter electrode where Faradaic reactions occur during regeneration. Up to approximately 40% of the initial capacity of the electrode is accessible following electrical regeneration., Comment: 14 pages, 9 figures, 2 tables

Cells adapt their metabolic fluxes in response to changes in the environment. We present a framework for the systematic construction of flux-based graphs derived from organism-wide metabolic networks. Our graphs encode the directionality of metabolic fluxes via edges that represent the flow of metabolites from source to target reactions. The methodology can be applied in the absence of a specific biological context by modelling fluxes probabilistically, or can be tailored to different environmental conditions by incorporating flux distributions computed through constraint-based approaches such as Flux Balance Analysis. We illustrate our approach on the central carbon metabolism of Escherichia coli and on a metabolic model of human hepatocytes. The flux-dependent graphs under various environmental conditions and genetic perturbations exhibit systemic changes in their topological and community structure, which capture the re-routing of metabolic fluxes and the varying importance of specific reactions and pathways. By integrating constraint-based models and tools from network science, our framework allows the study of context-specific metabolic responses at a system level beyond standard pathway descriptions., Comment: 26 Pages, 9 Figures

In this paper we study how to shape temporal pulses to switch a bistable system between its stable steady states. Our motivation for pulse-based control comes from applications in synthetic biology, where it is generally difficult to implement real-time feedback control systems due to technical limitations in sensors and actuators. We show that for monotone bistable systems, the estimation of the set of all pulses that switch the system reduces to the computation of one non-increasing curve. We provide an efficient algorithm to compute this curve and illustrate the results with a genetic bistable system commonly used in synthetic biology. We also extend these results to models with parametric uncertainty and provide a number of examples and counterexamples that demonstrate the power and limitations of the current theory. In order to show the full potential of the framework, we consider the problem of inducing oscillations in a monotone biochemical system using a combination of temporal pulses and event-based control. Our results provide an insight into the dynamics of bistable systems under external inputs and open up numerous directions for future investigation., Comment: 14 pages, contains material from the paper in Proc Amer Control Conf 2015, (pp. 3138-3143) and "Shaping pulses to control bistable systems analysis, computation and counterexamples", which is due to appear in Automatica

The regulation of metabolic activity by tuning enzyme expression levels is crucial to sustain cellular growth in changing environments. Metabolic networks are often studied at steady state using constraint-based models and optimization techniques. However, metabolic adaptations driven by changes in gene expression cannot be analyzed by steady state models, as these do not account for temporal changes in biomass composition. Here we present a dynamic optimization framework that integrates the metabolic network with the dynamics of biomass production and composition, explicitly taking into account enzyme production costs and enzymatic capacity. In contrast to the established dynamic flux balance analysis, our approach allows predicting dynamic changes in both the metabolic fluxes and the biomass composition during metabolic adaptations. We applied our algorithm in two case studies: a minimal nutrient uptake network, and an abstraction of core metabolic processes in bacteria. In the minimal model, we show that the optimized uptake rates reproduce the empirical Monod growth for bacterial cultures. For the network of core metabolic processes, the dynamic optimization algorithm predicted commonly observed metabolic adaptations, such as a diauxic switch with a preference ranking for different nutrients, re-utilization of waste products after depletion of the original substrate, and metabolic adaptation to an impending nutrient depletion. These examples illustrate how dynamic adaptations of enzyme expression can be predicted solely from an optimization principle.

This paper deals with temporal enzyme distribution in the activation of biochemical pathways. Pathway activation arises when production of a certain biomolecule is required due to changing environmental conditions. Under the premise that biological systems have been optimized through evolutionary processes, a biologically meaningful optimal control problem is posed. In this setup, the enzyme concentrations are assumed to be time dependent and constrained by a limited overall enzyme production capacity, while the optimization criterion accounts for both time and resource usage. Using geometric arguments we establish the bang-bang nature of the solution and reveal that each reaction must be sequentially activated in the same order as they appear in the pathway. The results hold for a broad range of enzyme dynamics which includes, but is not limited to, Mass Action, Michaelis-Menten and Hill Equation kinetics., Comment: 14 pages, 3 figures. Paper to be presented at the 17th IFAC World Congress, Seoul, Korea, July 2008