BioCRNpyler: Compiling chemical reaction networks from biomolecular parts in diverse contexts

Impact of Chemical Dynamics of Commercial PURE Systems on Malachite Green Aptamer Fluorescence

The malachite green aptamer (MGapt) is known for its utility in RNA measurement in vivo and in lysate-based cell-free protein systems. However, MGapt fluorescence dynamics do not accurately reflect RNA concentration. Our study finds that MGapt fluorescence is unstable in commercial PURE systems. We discovered that the chemical composition of the cell-free reaction strongly influences MGapt fluorescence, which leads to inaccurate RNA calculations. Specific to the commercial system, we posit that MGapt fluorescence is significantly affected by the system's chemical properties, governed notably by the presence of dithiothreitol (DTT). We propose a model that, on average, accurately predicts MGapt measurement within a 10% margin, leveraging DTT concentration as a critical factor. This model sheds light on the complex dynamics of MGapt in cell-free systems and underscores the importance of considering environmental factors in RNA measurements using aptamers.

Fine-tuned protein-lipid interactions in biological membranes: exploration and implications of the ORMDL-ceramide negative feedback loop in the endoplasmic reticulum

Biological membranes consist of a lipid bilayer in which integral membrane proteins are embedded. Based on the compositional complexity of the lipid species found in membranes, and on their specific and selective interactions with membrane proteins, we recently suggested that membrane bilayers can be best described as "finely-tuned molecular machines." We now discuss one such set of lipid-protein interactions by describing a negative feedback mechanism operating in the de novo sphingolipid biosynthetic pathway, which occurs in the membrane of the endoplasmic reticulum, and describe the atomic interactions between the first enzyme in the pathway, namely serine palmitoyl transferase, and the product of the fourth enzyme in the pathway, ceramide. We explore how hydrogen-bonding and hydrophobic interactions formed between Asn13 and Phe63 in the serine palmitoyl transferase complex and ceramide can influence the ceramide content of the endoplasmic reticulum. This example of finely-tuned biochemical interactions raises intriguing mechanistic questions about how sphingolipids and their biosynthetic enzymes could have evolved, particularly in light of their metabolic co-dependence.

Chemical Reaction Models in Synthetic Promoter Design in Bacteria

We discuss the formalism of chemical reaction networks (CRNs) as a computer-aided design interface for using formal methods in engineering living technologies. We set out by reviewing formal methods within a broader view of synthetic biology. Based on published results, we illustrate, step by step, how mathematical and computational techniques on CRNs can be used to study the structural and dynamic properties of the designed systems. As a case study, we use an E. coli two-component system that relays the external inorganic phosphate concentration signal to genetic components. We show how CRN models can scan and explore phenotypic regimes of synthetic promoters with varying detection thresholds, thereby providing a means for fine-tuning the promoter strength to match the specification.

Catalyst: Fast and flexible modeling of reaction networks

We introduce Catalyst.jl, a flexible and feature-filled Julia library for modeling and high-performance simulation of chemical reaction networks (CRNs). Catalyst supports simulating stochastic chemical kinetics (jump process), chemical Langevin equation (stochastic differential equation), and reaction rate equation (ordinary differential equation) representations for CRNs. Through comprehensive benchmarks, we demonstrate that Catalyst simulation runtimes are often one to two orders of magnitude faster than other popular tools. More broadly, Catalyst acts as both a domain-specific language and an intermediate representation for symbolically encoding CRN models as Julia-native objects. This enables a pipeline of symbolically specifying, analyzing, and modifying CRNs; converting Catalyst models to symbolic representations of concrete mathematical models; and generating compiled code for numerical solvers. Leveraging ModelingToolkit.jl and Symbolics.jl, Catalyst models can be analyzed, simplified, and compiled into optimized representations for use in numerical solvers. Finally, we demonstrate Catalyst's broad extensibility and composability by highlighting how it can compose with a variety of Julia libraries, and how existing open-source biological modeling projects have extended its intermediate representation.

Characterization of Integrase and Excisionase Activity in a Cell-Free Protein Expression System Using a Modeling and Analysis Pipeline

We present a full-stack modeling, analysis, and parameter identification pipeline to guide the modeling and design of biological systems starting from specifications to circuit implementations and parametrizations. We demonstrate this pipeline by characterizing the integrase and excisionase activity in a cell-free protein expression system. We build on existing Python tools─BioCRNpyler, AutoReduce, and Bioscrape─to create this pipeline. For enzyme-mediated DNA recombination in a cell-free system, we create detailed chemical reaction network models from simple high-level descriptions of the biological circuits and their context using BioCRNpyler. We use Bioscrape to show that the output of the detailed model is sensitive to many parameters. However, parameter identification is infeasible for this high-dimensional model; hence, we use AutoReduce to automatically obtain reduced models that have fewer parameters. This results in a hierarchy of reduced models under different assumptions to finally arrive at a minimal ODE model for each circuit. Then, we run sensitivity analysis-guided Bayesian inference using Bioscrape for each circuit to identify the model parameters. This process allows us to quantify integrase and excisionase activity in cell extracts enabling complex-circuit designs that depend on accurate control over protein expression levels through DNA recombination. The automated pipeline presented in this paper opens up a new approach to complex circuit design, modeling, reduction, and parametrization.

Rapid modeling of experimental molecular kinetics with simple electronic circuits instead of with complex differential equations

Kinetic modeling has relied on using a tedious number of mathematical equations to describe molecular kinetics in interacting reactions. The long list of differential equations with associated abstract variables and parameters inevitably hinders readers’ easy understanding of the models. However, the mathematical equations describing the kinetics of biochemical reactions can be exactly mapped to the dynamics of voltages and currents in simple electronic circuits wherein voltages represent molecular concentrations and currents represent molecular fluxes. For example, we theoretically derive and experimentally verify accurate circuit models for Michaelis-Menten kinetics. Then, we show that such circuit models can be scaled via simple wiring among circuit motifs to represent more and arbitrarily complex reactions. Hence, we can directly map reaction networks to equivalent circuit schematics in a rapid, quantitatively accurate, and intuitive fashion without needing mathematical equations. We verify experimentally that these circuit models are quantitatively accurate. Examples include 1) different mechanisms of competitive, noncompetitive, uncompetitive, and mixed enzyme inhibition, important for understanding pharmacokinetics; 2) product-feedback inhibition, common in biochemistry; 3) reversible reactions; 4) multi-substrate enzymatic reactions, both important in many metabolic pathways; and 5) translation and transcription dynamics in a cell-free system, which brings insight into the functioning of all gene-protein networks. We envision that circuit modeling and simulation could become a powerful scientific communication language and tool for quantitative studies of kinetics in biology and related fields.