Living systems are intrinsically nonequilibrium: They use metabolically derived chemical energy to power their emergent dynamics and self-organization. A crucial driver of these dynamics is the cellular cytoskeleton, a defining example of an active material where the energy injected by molecular motors cascades across length scales, allowing the material to break the constraints of thermodynamic equilibrium and display emergent nonequilibrium dynamics only possible due to the constant influx of energy. Notwithstanding recent experimental advances in the use of local probes to quantify entropy production and the breaking of detailed balance, little is known about the energetics of active materials or how energy propagates from the molecular to emergent length scales. Here, we use a recently developed picowatt calorimeter to experimentally measure the energetics of an active microtubule gel that displays emergent large-scale flows. We find that only approximately one-billionth of the system's total energy consumption contributes to these emergent flows. We develop a chemical kinetics model that quantitatively captures how the system's total thermal dissipation varies with ATP and microtubule concentrations but that breaks down at high motor concentration, signaling an interference between motors. Finally, we estimate how energy losses accumulate across scales. Taken together, these results highlight energetic efficiency as a key consideration for the engineering of active materials and are a powerful step toward developing a nonequilibrium thermody-namics of living systems.

Life is a nonequilibrium phenomenon: metabolism provides a continuous supply of energy that drives nearly all cellular processes. However, very little is known about how much energy different cellular processes use, i.e. their energetic costs. The most direct experimental measurements of these costs involve modulating the activity of cellular processes and determining the resulting changes in energetic fluxes. In this review, we present a flux balance framework to aid in the design and interpretation of such experiments, and discuss the challenges associated with measuring the relevant metabolic fluxes. We then describe selected techniques that enable measurement of these fluxes. Finally, we review prior experimental and theoretical work that has employed techniques from biochemistry and nonequilibrium physics to determine the energetic costs of cellular processes.

Multiplex assays of variant effect (MAVEs) are a family of methods that includes deep mutational scanning experiments on proteins and massively parallel reporter assays on gene regulatory sequences. Despite their increasing popularity, a general strategy for inferring quantitative models of genotype-phenotype maps from MAVE data is lacking. Here we introduce MAVE-NN, a neural-network-based Python package that implements a broadly applicable information-theoretic framework for learning genotype-phenotype maps-including biophysically interpretable models-from MAVE datasets. We demonstrate MAVE-NN in multiple biological contexts, and highlight the ability of our approach to deconvolve mutational effects from otherwise confounding experimental nonlinearities and noise.

DNA replication in eukaryotic cells initiates from replication origins that bind the Origin Recognition Complex (ORC). Origin establishment requires well-defined DNA sequence motifs in Saccharomyces cerevisiae and some other budding yeasts, but most eukaryotes lack sequence-specific origins. A 3.9 angstrom structure of S. cerevisiae ORC-Cdc6-Cdt1-Mcm2-7 (OCCM) bound to origin DNA revealed that a loop within Orc2 inserts into a DNA minor groove and an alpha-helix within Orc4 inserts into a DNA major groove. Using a massively parallel origin selection assay coupled with a custom mutual-information-based modeling approach, and a separate analysis of whole-genome replication profiling, here we show that the Orc4 alpha-helix contributes to the DNA sequence-specificity of origins in S. cerevisiae and Orc4 alpha-helix mutations change genome-wide origin firing patterns. The DNA sequence specificity of replication origins, mediated by the Orc4 alpha-helix, has co-evolved with the gain of ORC-Sir4-mediated gene silencing and the loss of RNA interference.

Advances in DNA sequencing have revolutionized our ability to read genomes. However, even in the most well-studied of organisms, the bacterium Escherichia coli, for ≈65% of promoters we remain ignorant of their regulation. Until we crack this regulatory Rosetta Stone, efforts to read and write genomes will remain haphazard. We introduce a new method, Reg-Seq, that links massively parallel reporter assays with mass spectrometry to produce a base pair resolution dissection of more than a E. coli promoters in 12 growth conditions. We demonstrate that the method recapitulates known regulatory information. Then, we examine regulatory architectures for more than 80 promoters which previously had no known regulatory information. In many cases, we also identify which transcription factors mediate their regulation. This method clears a path for highly multiplexed investigations of the regulatory genome of model organisms, with the potential of moving to an array of microbes of ecological and medical relevance.

Measurements in biology have reached a level of precision that demands quantitative modeling. This is particularly true in the field of gene regulation, where concepts from physics such as thermodynamics have allowed for accurate models to be made. Many issues remain. DNA sequencing is routine enough to sequence new genomes in days and cheap enough to use deep sequencing to perform precision measurements, but our ability to interpret the wealth of genomic data is lagging behind, especially in the realm of gene regulation. The primary reason is that we lack any information what so ever as to the basic regulatory details of approximately 65 percent of operons even in E. coli, the best understood organism in biology. As a result we cannot use our hard won modeling efforts to understand any of these operons. This work takes steps to address these issues. First we use 30 LacI mutants as a test case to prove that we can make quantitatively accurate models of gene expression and sequence-dependent binding energies of transcription factors and RNA polymerase. Next we note that much of the quantitative insight available on transcriptional regulation relies on work on only a few model regulatory systems such as LacI as was considered above. We develop an approach, through a combination of massively parallel reporter assays, mass spectrometry, and information-theoretic modeling that can be used to dissect bacterial promoters in a systematic and scalable way. We demonstrate that we can uncover a qualitative list of transcription factor binding sites as well as their associated quantitative details from both well-studied and previously uncharacterized promoters in E. coli. Finally we extend the above method to over 100 E. coli promoters using over 12 growth conditions. We show the method recapitulates known regulatory information. Then, we examine regulatory architectures for more than 80 promoters which previously had no known regulation. In many cases, we identify which transcription factors mediate their regulation. The method introduced clears a path for fully characterizing the regulatory genome of E. coli and advances towards the goal of using this method on a wide variety of other organisms including other prokaryotes and eukaryotes such as Drosophila melanogaster.

Despite the central importance of transcriptional regulation in biology, it has proven difficult to determine the regulatory mechanisms of individual genes, let alone entire gene networks. It is particularly difficult to decipher the biophysical mechanisms of transcriptional regulation in living cells and determine the energetic properties of binding sites for transcription factors and RNA polymerase. In this work, we present a strategy for dissecting transcriptional regulatory sequences using in vivo methods (massively parallel reporter assays) to formulate quantitative models that map a transcription factor binding site's DNA sequence to transcription factor-DNA binding energy. We use these models to predict the binding energies of transcription factor binding sites to within 1 k(B)T of their measured values. We further explore how such a sequence-energy mapping relates to the mechanisms of trancriptional regulation in various promoter contexts. Specifically, we show that our models can be used to design specific induction responses, analyze the effects of amino acid mutations on DNA sequence preference, and determine how regulatory context affects a transcription factor's sequence specificity.

Gene regulation is one of the most ubiquitous processes in biology. But while the catalog of bacterial genomes continues to expand rapidly, we remain ignorant about how almost all of the genes in these genomes are regulated. At present, characterizing the molecular mechanisms by which individual regulatory sequences operate requires focused efforts using low-throughput methods. Here we show how a combination of massively parallel reporter assays, mass spectrometry, and information-theoretic modeling can be used to dissect bacterial promoters in a systematic and scalable way. We demonstrate this method on both well-studied and previously uncharacterized promoters in the enteric bacterium Escherichia coli. In all cases we recover nucleotide-resolution models of promoter mechanism. For some promoters, including previously unannotated ones, the approach allowed us to further extract quantitative biophysical models describing input-output relationships. This method opens up the possibility of exhaustively dissecting the mechanisms of promoter function in E. coli and a wide range of other bacteria.

Massively parallel assays (MPAs) are being rapidly adopted for studying a wide range of DNA, RNA, and protein sequence-function relationships. However, the software available for quantitatively modeling these relationships is severely limited. Here we describe MPAthic, a software package that enables the rapid inference of such models from a variety of MPA datasets. Using both simulated and previously published data, we show that the modeling capabilities of MPAthic greatly improve on those of existing software. In particular, only MPAthic can accurately quantify the strength of epistatic interactions. These capabilities address a major need in the analysis of MPA data.