JuliaCon 2023

Julia Systems Biology
07-28, 10:30–11:30 (US/Eastern), 32-G449 (Kiva)

Julia has had the most developed ecosystem for differential equation modeling in simulation through the SciML organization for a while. Here we present a collection of talks by computational systems biologists in the community. The focus of the symposium will be to look at how SciML tools are being used in systems biology, how they can improve, and how we can take steps to increase collaboration throughout industry and academia.


Julia SciML is a collection of hundreds of scientific computing packages for the Julia programming language. Historically, SciML began as a group of fast ordinary differential equation (ODE) solvers. However, in the last half a decade, SciML has grown to include an extensive collection of packages geared toward modern scientific computing, with its own Computer Algebra System (CAS) and Domain Specific Language (DSL), the ability for heterogeneous computing (e.g., running code on GPUs), and seamless connection to machine learning packages.

Julia SysBio community has significantly benefited from the SciML ecosystem. One of the first efforts to make SciML usable to the SysBio community was to provide importers and bridging packages (e.g., CellMLToolkit.jl, SBMLToolkit.jl) that allow for standard and well-curated models in open-source repositories (CellML, SBML) to be imported and used by the SciML packages. CellML models are generally in the form of ODE and are usually solved with the help of ODE and partial differential equation (PDE) solvers. On the other hand, SBML models are usually presented as reaction networks and are better suited to Catalyst.jl.

With the maturation of the SciML infrastructure and the availability of bridging packages between the biological models and the SciML ecosystem, it is time for applications! In this mini-symposium, we showcase some recent applications of SciML in SysBio in academia and industry.

In addition to presenting various applications, our plan is to have a forum to provide feedbacks from the community to the core SciML developers in order to guide the future development of SciML and related packages to make Julia and SciML the best environment for SysBio applications.

Presenters:

  • Alex Cohen will be presenting on techniques for learning physically-constrained models in mode space to characterize the locomotion and behavioral states of animals that use undulatory locomotion.

  • Andrew Stine (United Therapeutics) will be presenting how Julia and SciML tools are used in QSP modeling teams at United Therapeutics.

  • Torsten Schenkel/Harry Saxton will be presenting lumped parameter models We introduce CirculatorySystemModels.jl (CSM), an acausal model library built on ModelingToolkit.jl, containing a wide range of cardiovascular elements for creating complex models. CSM provides a common API, with automatic generation of ODEs meaning every 0D model can be solved with the same workflow, allowing for full integration with the SciML framework.

  • Katy Norman (Sanofi) will be presenting how Julia and SciML tools are used in QSP modeling teams at Sanofi.

  • Otto Ritter (Merck) will be presenting how Julia, Applied Category Theory, and SciML tools are used in QSP modeling teams at Merck with ReactiveDynamics.jl.

  • Sam Isaacson will be discussing high level progression of the systems biology ecosystem in Julia, providing important context on future steps and ways that the community can improve and collaborate.

  • Don Elbert will be presenting on variable volume models and his experience as a SciML user. This will provide important context on pain points for adoption and what we can do to make Julia the primary choice for modeling and simulation of biology.

  • Wiktor Phillips will be presenting his acausal neuronal modeling and simulation package, Conductor.jl.

  • Shahriar Iravanian will be presenting simulating cardiac electrophysiological ionic models as 2D/3D PDEs with the help of Julia and SciML and showcase the power of composition in Julia by using dual complex numbers in CUDA kernels to solve a biological problem (calculating the membrane impedance of cardiac cells).

  • Torkel Loman will be presenting Catalyst.jl and how chemical reaction networks are simulated in Julia.

  • Sebastian Micluța-Câmpeanu/Paul Lang/Elisabeth Roesch (JuliaHub) will explore strategies to handle parameter unidentifiability and integrate open standards and neural model autocomplete in QSP workflows.

mind body

This speaker also appears in:

Paul Lang works in developing comprehensive models of biochemical reaction systems and
supporting software tools. He obtained a BSc degree in Molecular Biology at the University of
Graz (Austria) and an MSc degree in Molecular Health Sciences at ETH Zurich (Switzerland). He
then switched from experimental to computational research. During his PhD at the University of
Oxford (UK), Paul developed a rule-based cell cycle model that explains location and dynamics
of 16 observables in RPE1 cells. As visiting scholar at the Icahn School of Medicine at Mount
Sinai in New York (USA), he co-developed BpForms, a toolkit for concretely describing non-
canonical polymers to enable the construction of whole-cell models. Paul also co-developed SBMLToolkit.jl, which imports SBML models into the SciML ecosystem. Currently, Paul works for JuliaHub on developing a parameter optimization tool for quantitative systems pharmacology. He also helps clients to translate their models to ModelingToolkit.jl ODESystems.

This speaker also appears in:

Torsten Schenkel is a theoretical engineer. He is Associate Professor for Continuum Mechanics at Sheffield Hallam University.

With a background in aerospace engineering, Torsten's research is focused on biomedical problems, mostly concerning numerical modelling of the cardio-vascular system. Models range from 3-dimensional heart models to time-resolved 3d-flows in atherosclerotic vessels to 0D-lumped models of the whole circulatory system. Having started scientific programming with Fortran77, he now feels like he has come home with Julia.

This speaker also appears in:

Harry Saxton is a 2nd Year PhD student at Sheffield Hallam University. With a background in mathematics, Harry is interested in Lumped parameter modelling and inverse problems. Harry's main interests lie in the areas of Sensitivity, Identifiability and Bayesian analysis. However, any biological problem in which mathematics is applied is of interest. Harry has been a part of developing the Julia package CirculatorySystemModels.jl which is a package which allows simple, modular lumped parameter modelling. Harry has been using Julia a year and feels like it is the best programming community that is out there.

This speaker also appears in: