Timing molecular motion and production with a synthetic transcriptional clock

The realization of artificial biochemical reaction networks with unique functionality is one of the main challenges for the development of synthetic biology. Due to the reduced number of components, biochemical circuits constructed in vitro promise to be more amenable to systematic design and quantitative assessment than circuits embedded within living organisms. To make good on that promise, effective methods for composing subsystems into larger systems are needed. Here we used an artificial biochemical oscillator based on in vitro transcription and RNA degradation reactions to drive a variety of “load” processes such as the operation of a DNA-based nanomechanical device (“DNA tweezers”) or the production of a functional RNA molecule (an aptamer for malachite green). We implemented several mechanisms for coupling the load processes to the oscillator circuit and compared them based on how much the load affected the frequency and amplitude of the core oscillator, and how much of the load was effectively driven. Based on heuristic insights and computational modeling, an “insulator circuit” was developed, which strongly reduced the detrimental influence of the load on the oscillator circuit. Understanding how to design effective insulation between biochemical subsystems will be critical for the synthesis of larger and more complex systems

Transform-limited pulses are not optimal for resonant multiphoton transitions

Maximizing nonlinear light-matter interactions is a primary motive for compressing laser pulses to achieve ultrashort transform limited pulses. Here we show how, by appropriately shaping the pulses, resonant multiphoton transitions can be enhanced significantly beyond the level achieved by maximizing the pulse's peak intensity. We demonstrate the counterintuitive nature of this effect with an experiment in a resonant two-photon absorption, in which, by selectively removing certain spectral bands, the peak intensity of the pulse is reduced by a factor of 40, yet the absorption rate is doubled. Furthermore, by suitably designing the spectral phase of the pulse, we increase the absorption rate by a factor of 7.Comment: 4 pages, 3 figure

Rule-based modeling of biochemical systems with BioNetGen

Totowa, NJ. Please cite this article when referencing BioNetGen in future publications. Rule-based modeling involves the representation of molecules as structured objects and molecular interactions as rules for transforming the attributes of these objects. The approach is notable in that it allows one to systematically incorporate site-specific details about proteinprotein interactions into a model for the dynamics of a signal-transduction system, but the method has other applications as well, such as following the fates of individual carbon atoms in metabolic reactions. The consequences of protein-protein interactions are difficult to specify and track with a conventional modeling approach because of the large number of protein phosphoforms and protein complexes that these interactions potentially generate. Here, we focus on how a rule-based model is specified in the BioNetGen language (BNGL) and how a model specification is analyzed using the BioNetGen software tool. We also discuss new developments in rule-based modeling that should enable the construction and analyses of comprehensive models for signal transduction pathways and similarly large-scale models for other biochemical systems. Key Words: Computational systems biology; mathematical modeling; combinatorial complexity; software; formal languages; stochastic simulation; ordinary differential equations; protein-protein interactions; signal transduction; metabolic networks. 1

Two-Dimensional Electronic Spectroscopy of Chlorophyll a: Solvent Dependent Spectral Evolution

The interaction of the monomeric chlorophyll Q-band electronic transition with solvents of differing physical-chemical properties is investigated through two-dimensional electronic spectroscopy (2DES). Chlorophyll constitutes the key chromophore molecule in light harvesting complexes. It is well-known that the surrounding protein in the light harvesting complex fine-tunes chlorophyll electronic transitions to optimize energy transfer. Therefore, an understanding of the influence of the environment on the monomeric chlorophyll electronic transitions is important. The Q-band 2DES is inhomogeneous at early times, particularly in hydrogen bonding polar solvents, but also in nonpolar solvents like cyclohexane. Interestingly this inhomogeneity persists for long times, even up to the nanosecond time scale in some solvents. The reshaping of the 2DES occurs over multiple time scales and was assigned mainly to spectral diffusion. At early times the reshaping is Gaussian-like, hinting at a strong solvent reorganization effect. The temporal evolution of the 2DES response was analyzed in terms of a Brownian oscillator model. The spectral densities underpinning the Brownian oscillator fitting were recovered for the different solvents. The absorption spectra and Stokes shift were also properly described by this model. The extent and nature of inhomogeneous broadening was a strong function of solvent, being larger in H-bonding and viscous media and smaller in nonpolar solvents. The fastest spectral reshaping components were assigned to solvent dynamics, modified by interactions with the solute

Annotations for Rule-Based Models

The chapter reviews the syntax to store machine-readable annotations and describes the mapping between rule-based modelling entities (e.g., agents and rules) and these annotations. In particular, we review an annotation framework and the associated guidelines for annotating rule-based models of molecular interactions, encoded in the commonly used Kappa and BioNetGen languages, and present prototypes that can be used to extract and query the annotations. An ontology is used to annotate models and facilitate their description

Machines vs. Ensembles: Effective MAPK Signaling through Heterogeneous Sets of Protein Complexes

A grant from the One-University Open Access Fund at the University of Kansas was used to defray the author’s publication fees in this Open Access journal. The Open Access Fund, administered by librarians from the KU, KU Law, and KUMC libraries, is made possible by contributions from the offices of KU Provost, KU Vice Chancellor for Research & Graduate Studies, and KUMC Vice Chancellor for Research. For more information about the Open Access Fund, please see http://library.kumc.edu/authors-fund.xml.Despite the importance of intracellular signaling networks, there is currently no consensus regarding the fundamental nature of the protein complexes such networks employ. One prominent view involves stable signaling machines with well-defined quaternary structures. The combinatorial complexity of signaling networks has led to an opposing perspective, namely that signaling proceeds via heterogeneous pleiomorphic ensembles of transient complexes. Since many hypotheses regarding network function rely on how we conceptualize signaling complexes, resolving this issue is a central problem in systems biology. Unfortunately, direct experimental characterization of these complexes has proven technologically difficult, while combinatorial complexity has prevented traditional modeling methods from approaching this question. Here we employ rule-based modeling, a technique that overcomes these limitations, to construct a model of the yeast pheromone signaling network. We found that this model exhibits significant ensemble character while generating reliable responses that match experimental observations. To contrast the ensemble behavior, we constructed a model that employs hierarchical assembly pathways to produce scaffold-based signaling machines. We found that this machine model could not replicate the experimentally observed combinatorial inhibition that arises when the scaffold is overexpressed. This finding provides evidence against the hierarchical assembly of machines in the pheromone signaling network and suggests that machines and ensembles may serve distinct purposes in vivo. In some cases, e.g. core enzymatic activities like protein synthesis and degradation, machines assembled via hierarchical energy landscapes may provide functional stability for the cell. In other cases, such as signaling, ensembles may represent a form of weak linkage, facilitating variation and plasticity in network evolution. The capacity of ensembles to signal effectively will ultimately shape how we conceptualize the function, evolution and engineering of signaling networks

SBML Level 3: an extensible format for the exchange and reuse of biological models

Abstract Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction‐based models and packages that extend the core with features suited to other model types including constraint‐based models, reaction‐diffusion models, logical network models, and rule‐based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single‐cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution

A computational analysis of the dynamic roles of talin, Dok1, and PIPKI for integrin activation

Integrin signaling regulates cell migration and plays a pivotal role in developmental processes and cancer metastasis. Integrin signaling has been studied extensively and much data is available on pathway components and interactions. Yet the data is fragmented and an integrated model is missing. We use a rule-based modeling approach to integrate available data and test biological hypotheses regarding the role of talin, Dok1 and PIPKI in integrin activation. The detailed biochemical characterization of integrin signaling provides us with measured values for most of the kinetics parameters. However, measurements are not fully accurate and the cellular concentrations of signaling proteins are largely unknown and expected to vary substantially across different cellular conditions. By sampling model behaviors over the physiologically realistic parameter range we find that the model exhibits only two different qualitative behaviours and these depend mainly on the relative protein concentrations, which offers a powerful point of control to the cell. Our study highlights the necessity to characterize model behavior not for a single parameter optimum, but to identify parameter sets that characterize different signaling modes