Modular protein-oligonucleotide signal exchange.

While many methods are available to measure the concentrations of proteins in solution, the development of a method to quantitatively report both increases and decreases in different protein concentrations in real-time using changes in the concentrations of other molecules, such as DNA outputs, has remained a challenge. Here, we present a biomolecular reaction process that reports the concentration of an input protein in situ as the concentration of an output DNA oligonucleotide strand. This method uses DNA oligonucleotide aptamers that bind either to a specific protein selectively or to a complementary DNA oligonucleotide reversibly using toehold-mediated DNA strand-displacement. It is possible to choose the sequence of output strand almost independent of the sensing protein. Using this strategy, we implemented four different exchange processes to report the concentrations of clinically relevant human α-thrombin and vascular endothelial growth factor using changes in concentrations of DNA oligonucleotide outputs. These exchange processes can operate in tandem such that the same or different output signals can indicate changes in concentration of distinct or identical input proteins. The simplicity of our approach suggests a pathway to build devices that can direct diverse output responses in response to changes in concentrations of specific proteins.

Distinct timescales of RNA regulators enable the construction of a genetic pulse generator.

To build complex genetic networks with predictable behaviors, synthetic biologists use libraries of modular parts that can be characterized in isolation and assembled together to create programmable higher-order functions. Characterization experiments and computational models for gene regulatory parts operating in isolation are routinely used to predict the dynamics of interconnected parts and guide the construction of new synthetic devices. Here, we individually characterize two modes of RNA-based transcriptional regulation, using small transcription activating RNAs (STARs) and clustered regularly interspaced short palindromic repeats interference (CRISPRi), and show how their distinct regulatory timescales can be used to engineer a composed feedforward loop that creates a pulse of gene expression. We use a cell-free transcription-translation system (TXTL) to rapidly characterize the system, and we apply Bayesian inference to extract kinetic parameters for an ordinary differential equation-based mechanistic model. We then demonstrate in simulation and verify with TXTL experiments that the simultaneous regulation of a single gene target with STARs and CRISPRi leads to a pulse of gene expression. Our results suggest the modularity of the two regulators in an integrated genetic circuit, and we anticipate that construction and modeling frameworks that can leverage this modularity will become increasingly important as synthetic circuits increase in complexity.

Quantifiable identification of flow-limited ventilator dyssynchrony with the deformed lung ventilator model.

BACKGROUND: Ventilator dyssynchrony (VD) can worsen lung injury and is challenging to detect and quantify due to the complex variability in the dyssynchronous breaths. While machine learning (ML) approaches are useful for automating VD detection from the ventilator waveform data, scalable severity quantification and its association with pathogenesis and ventilator mechanics remain challenging. OBJECTIVE: We develop a systematic framework to quantify pathophysiological features observed in ventilator waveform signals such that they can be used to create feature-based severity stratification of VD breaths. METHODS: A mathematical model was developed to represent the pressure and volume waveforms of individual breaths in a feature-based parametric form. Model estimates of respiratory effort strength were used to assess the severity of flow-limited (FL)-VD breaths compared to normal breaths. A total of 93,007 breath waveforms from 13 patients were analyzed. RESULTS: A novel model-defined continuous severity marker was developed and used to estimate breath phenotypes of FL-VD breaths. The phenotypes had a predictive accuracy of over 97% with respect to the previously developed ML-VD identification algorithm. To understand the incidence of FL-VD breaths and their association with the patient state, these phenotypes were further successfully correlated with ventilator-measured parameters and electronic health records. CONCLUSION: This work provides a computational pipeline to identify and quantify the severity of FL-VD breaths and paves the way for a large-scale study of VD causes and effects. This approach has direct application to clinical practice and in meaningful knowledge extraction from the ventilator waveform data.

Observation of locked phase dynamics and enhanced frequency stability in synchronized micromechanical oscillators.

Even though synchronization in autonomous systems has been observed for over three centuries, reports of systematic experimental studies on synchronized oscillators are limited. Here, we report on observations of internal synchronization in coupled silicon micromechanical oscillators associated with a reduction in the relative phase random walk that is modulated by the magnitude of the reactive coupling force between the oscillators. Additionally, for the first time, a significant improvement in the frequency stability of synchronized micromechanical oscillators is reported. The concept presented here is scalable and could be suitably engineered to establish the basis for a new class of highly precise miniaturized clocks and frequency references.

A Damaged-Informed Lung Ventilator Model for Ventilator Waveforms.

Motivated by a desire to understand pulmonary physiology, scientists have developed physiological lung models of varying complexity. However, pathophysiology and interactions between human lungs and ventilators, e.g., ventilator-induced lung injury (VILI), present challenges for modeling efforts. This is because the real-world pressure and volume signals may be too complex for simple models to capture, and while complex models tend not to be estimable with clinical data, limiting clinical utility. To address this gap, in this manuscript we developed a new damaged-informed lung ventilator (DILV) model. This approach relies on mathematizing ventilator pressure and volume waveforms, including lung physiology, mechanical ventilation, and their interaction. The model begins with nominal waveforms and adds limited, clinically relevant, hypothesis-driven features to the waveform corresponding to pulmonary pathophysiology, patient-ventilator interaction, and ventilator settings. The DILV model parameters uniquely and reliably recapitulate these features while having enough flexibility to reproduce commonly observed variability in clinical (human) and laboratory (mouse) waveform data. We evaluate the proof-in-principle capabilities of our modeling approach by estimating 399 breaths collected for differently damaged lungs for tightly controlled measurements in mice and uncontrolled human intensive care unit data in the absence and presence of ventilator dyssynchrony. The cumulative value of mean squares error for the DILV model is, on average, ≈12 times less than the single compartment lung model for all the waveforms considered. Moreover, changes in the estimated parameters correctly correlate with known measures of lung physiology, including lung compliance as a baseline evaluation. Our long-term goal is to use the DILV model for clinical monitoring and research studies by providing high fidelity estimates of lung state and sources of VILI with an end goal of improving management of VILI and acute respiratory distress syndrome.

Reconfiguring DNA Nanotube Architectures via Selective Regulation of Terminating Structures.

Molecular assemblies inside cells often undergo structural reconfiguration in response to stimuli to alter their function. Adaptive reconfiguration of cytoskeletal networks, for example, enables cellular shape change, movement, and cargo transport and plays a key role in driving complex processes such as division and differentiation. The cellular cytoskeleton is a self-assembling polymer network composed of simple filaments, so reconfiguration often occurs through the rearrangement of its component filaments' connectivities. DNA nanotubes have emerged as promising building blocks for constructing programmable synthetic analogs of cytoskeletal networks. Nucleating seeds can control when and where nanotubes grow, and capping structures can bind nanotube ends to stop growth. Such seeding and capping structures, collectively called termini, can organize nanotubes into larger architectures. However, these structures cannot be selectively activated or inactivated in response to specific stimuli to rearrange nanotube architectures, a key property of cytoskeletal networks. Here, we demonstrate how selective regulation of the binding affinity of DNA nanotube termini for DNA nanotube monomers or nanotube ends can direct the reconfiguration of nanotube architectures. Using DNA hybridization and strand displacement reactions that specifically activate or inactivate four orthogonal nanotube termini, we demonstrate that nanotube architectures can be reconfigured by selective addition or removal of distinct termini. Finally, we show how terminus activation could be a sensitive detector and amplifier of a DNA sequence signal. These results could enable the development of adaptive and multifunctional materials or diagnostic tools.

Mathematical Models of Protease-Based Enzymatic Biosensors.

An important goal of synthetic biology is to build biosensors and circuits with well-defined input-output relationships that operate at speeds found in natural biological systems. However, for molecular computation, most commonly used genetic circuit elements typically involve several steps from input detection to output signal production: transcription, translation, and post-translational modifications. These multiple steps together require up to several hours to respond to a single stimulus, and this limits the overall speed and complexity of genetic circuits. To address this gap, molecular frameworks that rely exclusively on post-translational steps to realize reaction networks that can process inputs at a time scale of seconds to minutes have been proposed. Here, we build mathematical models of fast biosensors capable of producing Boolean logic functionality. We employ protease-based chemical and light-induced switches, investigate their operation, and provide selection guidelines for their use as on-off switches. As a proof of concept, we implement a rapamycin-induced switch in vitro and demonstrate that its response qualitatively agrees with the predictions from our models. We then use these switches as elementary blocks, developing models for biosensors that can perform OR and XOR Boolean logic computation while using reaction conditions as tuning parameters. We use sensitivity analysis to determine the time-dependent sensitivity of the output to proteolytic and protein-protein binding reaction parameters. These fast protease-based biosensors can be used to implement complex molecular circuits with a capability of processing multiple inputs controllably and algorithmically. Our framework for evaluating and optimizing circuit performance can be applied to other molecular logic circuits.

In vitro implementation of robust gene regulation in a synthetic biomolecular integral controller.

Feedback mechanisms play a critical role in the maintenance of cell homeostasis in the presence of disturbances and uncertainties. Motivated by the need to tune the dynamics and improve the robustness of gene circuits, biological engineers have proposed various designs that mimic natural molecular feedback control mechanisms. However, practical and predictable implementations have proved challenging because of the complexity of synthesis and analysis of complex biomolecular networks. Here, we analyze and experimentally validate a synthetic biomolecular controller executed in vitro. The controller ensures that gene expression rate tracks an externally imposed reference level, and achieves this goal even in the presence of certain kinds of disturbances. Our design relies upon an analog of the well-known principle of integral feedback in control theory. We implement the controller in an Escherichia coli cell-free transcription-translation system, which allows rapid prototyping and implementation. Modeling and theory guide experimental implementation with well-defined operational predictability.

Mathematical Modeling of RNA-Based Architectures for Closed Loop Control of Gene Expression.

Feedback allows biological systems to control gene expression precisely and reliably, even in the presence of uncertainty, by sensing and processing environmental changes. Taking inspiration from natural architectures, synthetic biologists have engineered feedback loops to tune the dynamics and improve the robustness and predictability of gene expression. However, experimental implementations of biomolecular control systems are still far from satisfying performance specifications typically achieved by electrical or mechanical control systems. To address this gap, we present mathematical models of biomolecular controllers that enable reference tracking, disturbance rejection, and tuning of the temporal response of gene expression. These controllers employ RNA transcriptional regulators to achieve closed loop control where feedback is introduced via molecular sequestration. Sensitivity analysis of the models allows us to identify which parameters influence the transient and steady state response of a target gene expression process, as well as which biologically plausible parameter values enable perfect reference tracking. We quantify performance using typical control theory metrics to characterize response properties and provide clear selection guidelines for practical applications. Our results indicate that RNA regulators are well-suited for building robust and precise feedback controllers for gene expression. Additionally, our approach illustrates several quantitative methods useful for assessing the performance of biomolecular feedback control systems.

Self-Assembly of Hierarchical DNA Nanotube Architectures with Well-Defined Geometries.

An essential motif for the assembly of biological materials such as actin at the scale of hundreds of nanometers and beyond is a network of one-dimensional fibers with well-defined geometry. Here, we demonstrate the programmed organization of DNA filaments into micron-scale architectures where component filaments are oriented at preprogrammed angles. We assemble L-, T-, and Y-shaped DNA origami junctions that nucleate two or three micron length DNA nanotubes at high yields. The angles between the nanotubes mirror the angles between the templates on the junctions, demonstrating that nanoscale structures can control precisely how micron-scale architectures form. The ability to precisely program filament orientation could allow the assembly of complex filament architectures in two and three dimensions, including circuit structures, bundles, and extended materials.

Terminating DNA Tile Assembly with Nanostructured Caps.

Precise control over the nucleation, growth, and termination of self-assembly processes is a fundamental tool for controlling product yield and assembly dynamics. Mechanisms for altering these processes programmatically could allow the use of simple components to self-assemble complex final products or to design processes allowing for dynamic assembly or reconfiguration. Here we use DNA tile self-assembly to develop general design principles for building complexes that can bind to a growing biomolecular assembly and terminate its growth by systematically characterizing how different DNA origami nanostructures interact with the growing ends of DNA tile nanotubes. We find that nanostructures that present binding interfaces for all of the binding sites on a growing facet can bind selectively to growing ends and stop growth when these interfaces are presented on either a rigid or floppy scaffold. In contrast, nucleation of nanotubes requires the presentation of binding sites in an arrangement that matches the shape of the structure's facet. As a result, it is possible to build nanostructures that can terminate the growth of existing nanotubes but cannot nucleate a new structure. The resulting design principles for constructing structures that direct nucleation and termination of the growth of one-dimensional nanostructures can also serve as a starting point for programmatically directing two- and three-dimensional crystallization processes using nanostructure design.

A self-regulating biomolecular comparator for processing oscillatory signals.

While many cellular processes are driven by biomolecular oscillators, precise control of a downstream on/off process by a biochemical oscillator signal can be difficult: over an oscillator's period, its output signal varies continuously between its amplitude limits and spends a significant fraction of the time at intermediate values between these limits. Further, the oscillator's output is often noisy, with particularly large variations in the amplitude. In electronic systems, an oscillating signal is generally processed by a downstream device such as a comparator that converts a potentially noisy oscillatory input into a square wave output that is predominantly in one of two well-defined on and off states. The comparator's output then controls downstream processes. We describe a method for constructing a synthetic biochemical device that likewise produces a square-wave-type biomolecular output for a variety of oscillatory inputs. The method relies on a separation of time scales between the slow rate of production of an oscillatory signal molecule and the fast rates of intermolecular binding and conformational changes. We show how to control the characteristics of the output by varying the concentrations of the species and the reaction rates. We then use this control to show how our approach could be applied to process different in vitro and in vivo biomolecular oscillators, including the p53-Mdm2 transcriptional oscillator and two types of in vitro transcriptional oscillators. These results demonstrate how modular biomolecular circuits could, in principle, be combined to build complex dynamical systems. The simplicity of our approach also suggests that natural molecular circuits may process some biomolecular oscillator outputs before they are applied downstream.

Modeling nonlinearities in MEMS oscillators.

We present a mathematical model of a microelectromechanical system (MEMS) oscillator that integrates the nonlinearities of the MEMS resonator and the oscillator circuitry in a single numerical modeling environment. This is achieved by transforming the conventional nonlinear mechanical model into the electrical domain while simultaneously considering the prominent nonlinearities of the resonator. The proposed nonlinear electrical model is validated by comparing the simulated amplitude-frequency response with measurements on an open-loop electrically addressed flexural silicon MEMS resonator driven to large motional amplitudes. Next, the essential nonlinearities in the oscillator circuit are investigated and a mathematical model of a MEMS oscillator is proposed that integrates the nonlinearities of the resonator. The concept is illustrated for MEMS transimpedance-amplifier- based square-wave and sine-wave oscillators. Closed-form expressions of steady-state output power and output frequency are derived for both oscillator models and compared with experimental and simulation results, with a good match in the predicted trends in all three cases.