The turbulent soundscape of intertidal oyster reefs.

Turbulence and sound are important cues for oyster reef larval recruitment. Numerous studies have found a relationship between turbulence intensity and swimming behaviors of marine larvae, while others have documented the importance of sounds in enhancing larval recruitment to oyster reefs. However, the relationship between turbulence and the reef soundscape is not well understood. In this study we made side-by-side acoustic Doppler velocimeter turbulence measurements and hydrophone soundscape recordings over 2 intertidal oyster reefs (1 natural and 1 restored) and 1 adjacent bare mudflat as a reference. Sound pressure levels (SPL) were similar across all three sites, although SPL > 2000 Hz was highest at the restored reef, likely due to its larger area that contained a greater number of sound-producing organisms. Flow noise (FN), defined as the mean of pressure fluctuations recorded by the hydrophone at f < 100 Hz, was significantly related to mean flow speed, turbulent kinetic energy, and turbulence dissipation rate (ε), agreeing with theoretical calculations for turbulence. Our results also show a similar relationship between ε and FN to what has been previously reported for ε vs. downward larval swimming velocity (w b ), with both FN and w b demonstrating rapid growth at ε > 0.1 cm2 s-3. These results suggest that reef turbulence and sounds may attract oyster larvae in complementary and synergistic ways.

Signal integration and integral feedback control with biochemical reaction networks.

Biochemical reaction networks perform a variety of signal processing functions, one of which is computing the integrals of signal values. This is often used in integral feedback control, where it enables a system's output to respond to changing inputs, but to then return exactly back to some pre-determined setpoint value afterward. To gain a deeper understanding of how biochemical networks are able to both integrate signals and perform integral feedback control, we investigated these abilities for several simple reaction networks. We found imperfect overlap between these categories, with some networks able to perform both tasks, some able to perform integration but not integral feedback control, and some the other way around. Nevertheless, networks that could either integrate or perform integral feedback control shared key elements. In particular, they included a chemical species that was neutrally stable in the open loop system (no feedback), meaning that this species does not have a unique stable steady-state concentration. Neutral stability could arise from zeroth order decay reactions, binding to a partner that was produced at a constant rate (which occurs in antithetic control), or through a long chain of covalent cycles. Mathematically, it arose from rate equations for the reaction network that were underdetermined when evaluated at steady-state.

Design patterns of biological cells.

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

Design Patterns of Biological Cells.

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology in order to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

BioSimulators: a central registry of simulation engines and services for recommending specific tools.

Computational models have great potential to accelerate bioscience, bioengineering, and medicine. However, it remains challenging to reproduce and reuse simulations, in part, because the numerous formats and methods for simulating various subsystems and scales remain siloed by different software tools. For example, each tool must be executed through a distinct interface. To help investigators find and use simulation tools, we developed BioSimulators (https://biosimulators.org), a central registry of the capabilities of simulation tools and consistent Python, command-line and containerized interfaces to each version of each tool. The foundation of BioSimulators is standards, such as CellML, SBML, SED-ML and the COMBINE archive format, and validation tools for simulation projects and simulation tools that ensure these standards are used consistently. To help modelers find tools for particular projects, we have also used the registry to develop recommendation services. We anticipate that BioSimulators will help modelers exchange, reproduce, and combine simulations.

Dynamics and Sensitivity of Signaling Pathways.

Purpose of Review: Signaling pathways serve to communicate information about extracellular conditions into the cell, to both the nucleus and cytoplasmic processes to control cell responses. Genetic mutations in signaling network components are frequently associated with cancer and can result in cells acquiring an ability to divide and grow uncontrollably. Because signaling pathways play such a significant role in cancer initiation and advancement, their constituent proteins are attractive therapeutic targets. In this review, we discuss how signaling pathway modeling can assist with identifying effective drugs for treating diseases, such as cancer. An achievement that would facilitate the use of such models is their ability to identify controlling biochemical parameters in signaling pathways, such as molecular abundances and chemical reaction rates, because this would help determine effective points of attack by therapeutics. Recent Findings: We summarize the current state of understanding the sensitivity of phosphorylation cycles with and without sequestration. We also describe some basic properties of regulatory motifs including feedback and feedforward regulation. Summary: Although much recent work has focused on understanding the dynamics and particularly the sensitivity of signaling networks in eukaryotic systems, there is still an urgent need to build more scalable models of signaling networks that can appropriately represent their complexity across different cell types and tumors.

Python interfaces for the Smoldyn simulator.

MOTIVATION: Smoldyn is a particle-based biochemical simulator that is frequently used for systems biology and biophysics research. Previously, users could only define models using text-based input or a C/C++ application programming interface (API), which were convenient, but limited extensibility. RESULTS: We added a Python API to Smoldyn to improve integration with other software tools, such as Jupyter notebooks, other Python code libraries and other simulators. It includes low-level functions that closely mimic the existing C/C++ API and higher-level functions that are more convenient to use. These latter functions follow modern object-oriented Python conventions. AVAILABILITY AND IMPLEMENTATION: Smoldyn is open source and free, available at http://www.smoldyn.org and can be installed with the Python package manager pip. It runs on Mac, Windows and Linux.Documentation is available at http://www.smoldyn.org/SmoldynManual.pdf and https://smoldyn.readthedocs.io/en/latest/python/api.html.

Effects of surfaces and macromolecular crowding on bimolecular reaction rates.

Biological cells are complex environments that are densely packed with macromolecules and subdivided by membranes, both of which affect the rates of chemical reactions. It is well known that crowding reduces the volume available to reactants, which increases reaction rates, and also inhibits reactant diffusion, which decreases reaction rates. This work investigates these effects quantitatively using analytical theory and particle-based simulations. A reaction rate equation based on only these two processes turned out to be inconsistent with simulation results. However, accounting for diffusion inhibition by the surfaces of nearby obstacles, which affects access to reactants, it led to perfect agreement for reactions near impermeable planar membranes and improved agreement for reactions in crowded spaces. A separate model that quantified reactant occlusion by crowders, and extensions to a thermodynamic 'cavity' model proposed by Berezhkovskii and Szabo [25], were comparably successful. These results help elucidate reaction dynamics in confined spaces and improve prediction of in vivo reaction rates from in vitro ones.

Rule-Based Modeling Using Wildcards in the Smoldyn Simulator.

Many biological molecules exist in multiple variants, such as proteins with different posttranslational modifications, DNAs with different sequences, and phospholipids with different chain lengths. Representing these variants as distinct species, as most biochemical simulators do, leads to the problem that the number of species, and chemical reactions that interconvert them, typically increase combinatorially with the number of ways that the molecules can vary. This can be alleviated by "rule-based modeling methods," in which software generates the chemical reaction network from relatively simple "rules." This chapter presents a new approach to rule-based modeling. It is based on wildcards that match to species names, much as wildcards can match to file names in computer operating systems. It is much simpler to use than the formal rule-based modeling approaches developed previously but can lead to unintended consequences if not used carefully. This chapter demonstrates rule-based modeling with wildcards through examples for signaling systems, protein complexation, polymerization, nucleic acid sequence copying and mutation, the "SMILES" chemical notation, and others. The method is implemented in Smoldyn, a spatial and stochastic biochemical simulator, for both generate-first and on-the-fly expansion, meaning whether the reaction network is generated before or during the simulation.

Transferring information without distortion.

Despite employing diverse molecular mechanisms, many different cell signaling systems avoid losing information by transmitting it in a linear manner.

Smoldyn: particle-based simulation with rule-based modeling, improved molecular interaction and a library interface.

Motivation: Smoldyn is a spatial and stochastic biochemical simulator. It treats each molecule of interest as an individual particle in continuous space, simulating molecular diffusion, molecule-membrane interactions and chemical reactions, all with good accuracy. This article presents several new features. Results: Smoldyn now supports two types of rule-based modeling. These are a wildcard method, which is very convenient, and the BioNetGen package with extensions for spatial simulation, which is better for complicated models. Smoldyn also includes new algorithms for simulating the diffusion of surface-bound molecules and molecules with excluded volume. Both are exact in the limit of short time steps and reasonably good with longer steps. In addition, Smoldyn supports single-molecule tracking simulations. Finally, the Smoldyn source code can be accessed through a C/C ++ language library interface. Availability and Implementation: Smoldyn software, documentation, code, and examples are at http://www.smoldyn.org . Contact: steven.s.andrews@gmail.com.

Push-Pull and Feedback Mechanisms Can Align Signaling System Outputs with Inputs.

Many cell signaling systems, including the yeast pheromone response system, exhibit "dose-response alignment" (DoRA), in which output of one or more downstream steps closely matches the fraction of occupied receptors. DoRA can improve the fidelity of transmitted dose information. Here, we searched systematically for biochemical network topologies that produced DoRA. Most networks, including many containing feedback and feedforward loops, could not produce DoRA. However, networks including "push-pull" mechanisms, in which the active form of a signaling species stimulates downstream activity and the nominally inactive form reduces downstream activity, enabled perfect DoRA. Networks containing feedbacks enabled DoRA, but only if they also compared feedback to input and adjusted output to match. Our results establish push-pull as a non-feedback mechanism to align output with variable input and maximize information transfer in signaling systems. They also suggest genetic approaches to determine whether particular signaling systems use feedback or push-pull control.

Crosstalk between the lipopolysaccharide and phospholipid pathways during outer membrane biogenesis in Escherichia coli.

The outer membrane of gram-negative bacteria is composed of phospholipids in the inner leaflet and lipopolysaccharides (LPS) in the outer leaflet. LPS is an endotoxin that elicits a strong immune response from humans, and its biosynthesis is in part regulated via degradation of LpxC (EC 3.5.1.108) and WaaA (EC 2.4.99.12/13) enzymes by the protease FtsH (EC 3.4.24.-). Because the synthetic pathways for both molecules are complex, in addition to being produced in strict ratios, we developed a computational model to interrogate the regulatory mechanisms involved. Our model findings indicate that the catalytic activity of LpxK (EC 2.7.1.130) appears to be dependent on the concentration of unsaturated fatty acids. This is biologically important because it assists in maintaining LPS/phospholipids homeostasis. Further crosstalk between the phospholipid and LPS biosynthetic pathways was revealed by experimental observations that LpxC is additionally regulated by an unidentified protease whose activity is independent of lipid A disaccharide concentration (the feedback source for FtsH-mediated LpxC regulation) but could be induced in vitro by palmitic acid. Further experimental analysis provided evidence on the rationale for WaaA regulation. Overexpression of waaA resulted in increased levels of 3-deoxy-d-manno-oct-2-ulosonic acid (Kdo) sugar in membrane extracts, whereas Kdo and heptose levels were not elevated in LPS. This implies that uncontrolled production of WaaA does not increase the LPS production rate but rather reglycosylates lipid A precursors. Overall, the findings of this work provide previously unidentified insights into the complex biogenesis of the Escherichia coli outer membrane.

A Method and On-Line Tool for Maximum Likelihood Calibration of Immunoblots and Other Measurements That Are Quantified in Batches.

Experimental measurements require calibration to transform measured signals into physically meaningful values. The conventional approach has two steps: the experimenter deduces a conversion function using measurements on standards and then calibrates (or normalizes) measurements on unknown samples with this function. The deduction of the conversion function from only the standard measurements causes the results to be quite sensitive to experimental noise. It also implies that any data collected without reliable standards must be discarded. Here we show that a "1-step calibration method" reduces these problems for the common situation in which samples are measured in batches, where a batch could be an immunoblot (Western blot), an enzyme-linked immunosorbent assay (ELISA), a sequence of spectra, or a microarray, provided that some sample measurements are replicated across multiple batches. The 1-step method computes all calibration results iteratively from all measurements. It returns the most probable values for the sample compositions under the assumptions of a statistical model, making them the maximum likelihood predictors. It is less sensitive to measurement error on standards and enables use of some batches that do not include standards. In direct comparison of both real and simulated immunoblot data, the 1-step method consistently exhibited smaller errors than the conventional "2-step" method. These results suggest that the 1-step method is likely to be most useful for cases where experimenters want to analyze existing data that are missing some standard measurements and where experimenters want to extract the best results possible from their data. Open source software for both methods is available for download or on-line use.

Multiscale reaction-diffusion simulations with Smoldyn.

UNLABELLED: Smoldyn is a software package for stochastic modelling of spatial biochemical networks and intracellular systems. It was originally developed with an accurate off-lattice particle-based model at its core. This has recently been enhanced with the addition of a computationally efficient on-lattice model, which can be run stand-alone or coupled together for multiscale simulations using both models in regions where they are most required, increasing the applicability of Smoldyn to larger molecule numbers and spatial domains. Simulations can switch between models with only small additions to their configuration file, enabling users with existing Smoldyn configuration files to run the new on-lattice model with any reaction, species or surface descriptions they might already have. AVAILABILITY AND IMPLEMENTATION: Source code and binaries freely available for download at www.smoldyn.org, implemented in C/C++ and supported on Linux, Mac OSX and MS Windows.

An integrated model of transcription factor diffusion shows the importance of intersegmental transfer and quaternary protein structure for target site finding.

We present a computational model of transcription factor motion that explains both the observed rapid target finding of transcription factors, and how this motion influences protein and genome structure. Using the Smoldyn software, we modelled transcription factor motion arising from a combination of unrestricted 3D diffusion in the nucleoplasm, sliding along the DNA filament, and transferring directly between filament sections by intersegmental transfer. This presents a fine-grain picture of the way in which transcription factors find their targets two orders of magnitude faster than 3D diffusion alone allows. Eukaryotic genomes contain sections of nucleosome free regions (NFRs) around the promoters; our model shows that the presence and size of these NFRs can be explained as their acting as antennas on which transcription factors slide to reach their targets. Additionally, our model shows that intersegmental transfer may have shaped the quaternary structure of transcription factors: sequence specific DNA binding proteins are unusually enriched in dimers and tetramers, perhaps because these allow intersegmental transfer, which accelerates target site finding. Finally, our model shows that a 'hopping' motion can emerge from 3D diffusion on small scales. This explains the apparently long sliding lengths that have been observed for some DNA binding proteins observed in vitro. Together, these results suggest that transcription factor diffusion dynamics help drive the evolution of protein and genome structure.

Methods for modeling cytoskeletal and DNA filaments.

This review summarizes the models that researchers use to represent the conformations and dynamics of cytoskeletal and DNA filaments. It focuses on models that address individual filaments in continuous space. Conformation models include the freely jointed, Gaussian, angle-biased chain (ABC), and wormlike chain (WLC) models, of which the first three bend at discrete joints and the last bends continuously. Predictions from the WLC model generally agree well with experiment. Dynamics models include the Rouse, Zimm, stiff rod, dynamic WLC, and reptation models, of which the first four apply to isolated filaments and the last to entangled filaments. Experiments show that the dynamic WLC and reptation models are most accurate. They also show that biological filaments typically experience strong hydrodynamic coupling and/or constrained motion. Computer simulation methods that address filament dynamics typically compute filament segment velocities from local forces using the Langevin equation and then integrate these velocities with explicit or implicit methods; the former are more versatile and the latter are more efficient. Much remains to be discovered in biological filament modeling. In particular, filament dynamics in living cells are not well understood, and current computational methods are too slow and not sufficiently versatile. Although primarily a review, this paper also presents new statistical calculations for the ABC and WLC models. Additionally, it corrects several discrepancies in the literature about bending and torsional persistence length definitions, and their relations to flexural and torsional rigidities.

A Complete Pathway Model for Lipid A Biosynthesis in Escherichia coli.

Lipid A is a highly conserved component of lipopolysaccharide (LPS), itself a major component of the outer membrane of Gram-negative bacteria. Lipid A is essential to cells and elicits a strong immune response from humans and other animals. We developed a quantitative model of the nine enzyme-catalyzed steps of Escherichia coli lipid A biosynthesis, drawing parameters from the experimental literature. This model accounts for biosynthesis regulation, which occurs through regulated degradation of the LpxC and WaaA (also called KdtA) enzymes. The LpxC degradation signal appears to arise from the lipid A disaccharide concentration, which we deduced from prior results, model results, and new LpxK overexpression results. The model agrees reasonably well with many experimental findings, including the lipid A production rate, the behaviors of mutants with defective LpxA enzymes, correlations between LpxC half-lives and cell generation times, and the effects of LpxK overexpression on LpxC concentrations. Its predictions also differ from some experimental results, which suggest modifications to the current understanding of the lipid A pathway, such as the possibility that LpxD can replace LpxA and that there may be metabolic channeling between LpxH and LpxB. The model shows that WaaA regulation may serve to regulate the lipid A production rate when the 3-deoxy-D-manno-oct-2-ulosonic acid (KDO) concentration is low and/or to control the number of KDO residues that get attached to lipid A. Computation of flux control coefficients showed that LpxC is the rate-limiting enzyme if pathway regulation is ignored, but that LpxK is the rate-limiting enzyme if pathway regulation is present, as it is in real cells. Control also shifts to other enzymes if the pathway substrate concentrations are not in excess. Based on these results, we suggest that LpxK may be a much better drug target than LpxC, which has been pursued most often.

The autoregulation of a eukaryotic DNA transposon.

How do DNA transposons live in harmony with their hosts? Bacteria provide the only documented mechanisms for autoregulation, but these are incompatible with eukaryotic cell biology. Here we show that autoregulation of Hsmar1 operates during assembly of the transpososome and arises from the multimeric state of the transposase, mediated by a competition for binding sites. We explore the dynamics of a genomic invasion using a computer model, supported by in vitro and in vivo experiments, and show that amplification accelerates at first but then achieves a constant rate. The rate is proportional to the genome size and inversely proportional to transposase expression and its affinity for the transposon ends. Mariner transposons may therefore resist post-transcriptional silencing. Because regulation is an emergent property of the reaction it is resistant to selfish exploitation. The behavior of distantly related eukaryotic transposons is consistent with the same mechanism, which may therefore be widely applicable. DOI:http://dx.doi.org/10.7554/eLife.00668.001.

Spatial and stochastic cellular modeling with the Smoldyn simulator.

This chapter describes how to use Smoldyn, which is a computer program for modeling cellular systems with spatial and stochastic detail. Smoldyn represents each molecule of interest as an individual point-like particle. These simulated molecules diffuse, interact with surfaces (e.g., biological membranes), and undergo chemical reactions much as they would in real biochemical systems. Smoldyn has been used to model signal transduction within bacterial cells, pheromone signaling between yeast cells, bacterial carboxysome function, diffusion in crowded spaces, and many other systems. A new "rule-based modeling" feature automatically generates chemical species and reactions as they arise in simulations due to protein modifications and complexation. Smoldyn is easy to use, quantitatively accurate, and computationally efficient. It is generally best for systems with length scales between nanometers and several microns, time scales from tens of nanoseconds to tens of minutes, and up to about 10(5) individual molecules. Smoldyn runs on Macintosh, Linux, or Windows systems, is open source, and can be downloaded from http://www.smoldyn.org.