-
SBML2Julia: interfacing SBML with efficient nonlinear Julia modelling and solution tools for parameter optimization
Authors:
Paul F. Lang,
Sungho Shin,
Victor M. Zavala
Abstract:
Motivation: Estimating model parameters from experimental observations is one of the key challenges in systems biology and can be computationally very expensive. While the Julia programming language was recently developed as a high-level and high-performance language for scientific computing, systems biologists have only started to realise its potential. For instance, we have recently used Julia t…
▽ More
Motivation: Estimating model parameters from experimental observations is one of the key challenges in systems biology and can be computationally very expensive. While the Julia programming language was recently developed as a high-level and high-performance language for scientific computing, systems biologists have only started to realise its potential. For instance, we have recently used Julia to cut down the optimization time of a microbial community model by a factor of 140. To facilitate access of the systems biology community to the efficient nonlinear solvers used for this optimisation, we developed SBML2Julia. SBML2Julia translates optimisation problems specified in SBML and TSV files (PEtab format) into Julia for Mathematical Programming (JuMP), executes the optimization and returns the results in tabular format.
Availability and implementation: SBML2Julia is freely available under the MIT license. It comes with a command line interface and Python API. Internally, SBML2Julia calls the Julia LTS release v1.0.5 for optimisation. All necessary dependencies can be pulled from Docker Hub (https://hub.docker.com/repository/docker/paulflang/sbml2julia). Source code and documentation are available at https://github.com/paulflang/SBML2Julia.
△ Less
Submitted 4 November, 2020;
originally announced November 2020.
-
Water Dynamics Around Proteins: T- and R-States of Hemoglobin and Melittin
Authors:
Marco Pezzella,
Krystel El Hage,
Michiel J. M. Niesen,
Sucheol Shin,
Adam P. Willard,
Markus Meuwly,
Martin Karplus
Abstract:
The water dynamics, as characterized by the local hydrophobicity (LH), is investigated for tetrameric hemoglobin and dimeric melittin. For the T0 to R0 transition in Hb it is found that LH provides additional molecular-level insight into the Perutz mechanism, i.e., the breaking and formation of salt bridges at the alpha1 / beta2 and alpha2 / beta1 interface is accompanied by changes in LH. For Hb…
▽ More
The water dynamics, as characterized by the local hydrophobicity (LH), is investigated for tetrameric hemoglobin and dimeric melittin. For the T0 to R0 transition in Hb it is found that LH provides additional molecular-level insight into the Perutz mechanism, i.e., the breaking and formation of salt bridges at the alpha1 / beta2 and alpha2 / beta1 interface is accompanied by changes in LH. For Hb in cubic water boxes with 90 Aengstroem and 120 Aengstroem edge length it is observed that following a decrease in LH as a consequence of reduced water density or change of water orientation at the protein/water interface the alpha / beta interfaces are destabilized; this is a hallmark of the Perutz stereochemical model for the T to R transition in Hb. The present work thus provides a dynamical view of the classical structural model relevant to the molecular foundations of Hb function. For dimeric melittin, earlier results by Cheng and Rossky (Nature, 1998, 392, 696-699) are confirmed and interpreted on the basis of LH from simulations in which the protein structure is frozen. For the flexible melittin dimer the changes in the local hydration can be as much as 30 % than for the rigid dimer, reflecting the fact that protein and water dynamics are coupled.
△ Less
Submitted 26 June, 2020;
originally announced June 2020.
-
A domain-level DNA strand displacement reaction enumerator allowing arbitrary non-pseudoknotted secondary structures
Authors:
Casey Grun,
Karthik Sarma,
Brian Wolfe,
Seung Woo Shin,
Erik Winfree
Abstract:
DNA strand displacement systems have proven themselves to be fertile substrates for the design of programmable molecular machinery and circuitry. Domain-level reaction enumerators provide the foundations for molecular programming languages by formalizing DNA strand displacement mechanisms and modeling interactions at the "domain" level - one level of abstraction above models that explicitly descri…
▽ More
DNA strand displacement systems have proven themselves to be fertile substrates for the design of programmable molecular machinery and circuitry. Domain-level reaction enumerators provide the foundations for molecular programming languages by formalizing DNA strand displacement mechanisms and modeling interactions at the "domain" level - one level of abstraction above models that explicitly describe DNA strand sequences. Unfortunately, the most-developed models currently only treat pseudo-linear DNA structures, while many systems being experimentally and theoretically pursued exploit a much broader range of secondary structure configurations. Here, we describe a new domain-level reaction enumerator that can handle arbitrary non-pseudoknotted secondary structures and reaction mechanisms including association and dissociation, 3-way and 4-way branch migration, and direct as well as remote toehold activation. To avoid polymerization that is inherent when considering general structures, we employ a time-scale separation technique that holds in the limit of low concentrations. This also allows us to "condense" the detailed reactions by eliminating fast transients, with provable guarantees of correctness for the set of reactions and their kinetics. We hope that the new reaction enumerator will be used in new molecular programming languages, compilers, and tools for analysis and verification that treat a wider variety of mechanisms of interest to experimental and theoretical work. We have implemented this enumerator in Python, and it is included in the DyNAMiC Workbench Integrated Development Environment.
△ Less
Submitted 11 May, 2015;
originally announced May 2015.
-
Verifying Chemical Reaction Network Implementations: A Pathway Decomposition Approach
Authors:
Seung Woo Shin,
Chris Thachuk,
Erik Winfree
Abstract:
Here we focus on the challenge of verifying the correctness of molecular implementations of abstract chemical reaction networks, where operation in a well-mixed "soup" of molecules is stochastic, asynchronous, concurrent, and often involves multiple intermediate steps in the implementation, parallel pathways, and side reactions. This problem relates to the verification of Petri nets, but existing…
▽ More
Here we focus on the challenge of verifying the correctness of molecular implementations of abstract chemical reaction networks, where operation in a well-mixed "soup" of molecules is stochastic, asynchronous, concurrent, and often involves multiple intermediate steps in the implementation, parallel pathways, and side reactions. This problem relates to the verification of Petri nets, but existing approaches are not sufficient for providing a single guarantee covering an infinite set of possible initial states (molecule counts) and an infinite state space potentially explored by the system given any initial state. We address these issues by formulating a new theory of pathway decomposition that provides an elegant formal basis for comparing chemical reaction network implementations, and we present an algorithm that computes this basis. Our theory naturally handles certain situations that commonly arise in molecular implementations, such as what we call "delayed choice," that are not easily accommodated by other approaches. We further show how pathway decomposition can be combined with weak bisimulation to handle a wider class that includes most currently known enzyme-free DNA implementation techniques. We anticipate that our notion of logical equivalence between chemical reaction network implementations will be valuable for other molecular implementations such as biochemical enzyme systems, and perhaps even more broadly in concurrency theory.
△ Less
Submitted 10 August, 2017; v1 submitted 3 November, 2014;
originally announced November 2014.