Skip to main content

Life123 : Quantitative Dynamical Modeling of Biological Systems

(UPDATED 8/2022) - Are we ready to embark on a next-generation detailed quantitative modeling of complex biological systems, including whole-cell simulations? 

An anticipated up-jump in computing power may be imminent from Photonics computers (which I discuss here), and GPU's are rapidly gaining power as well...  Are we in ready state to put existing - and upcoming - power to good use?

This is a manifest, and a call to action

What's Life123?

It's about detailed quantitative modeling of biological systems in 1-D, 2-D and full 3-D, as well as a multi-faceted software platform for doing so.

What's (pseudo-)1D?  For now, let's say it's like the inside of a long, thin tube - with no interactions with the tube.  Likewise, (pseudo-)2D can be thought of as a Petri dish, with no interactions with the lid or the bottom.

Website: https://life123.science

A purposeful decision to also utilize 1D and 2D

But why?  Yes, it's in part about "walk before you run"...  but, more specifically, it's about the freedom to sometimes choose to avoid the distraction of higher dimensions, and focus on the essence of the features... and focus on incorporating good habits very early on, before the complexity, the long runs and the difficult visualizations of 3D come into being.

What good habits?  Just to mention a few:

COMPUTING:  GPU-assisted computing, parallelized computations utilizing multiple CPU cores and/or multiple computers.

MODELING:  variable spacial and temporal resolutions.  A modular approach of coarse-to-fine models as needed.  From very early on, the project will model reaction rates, chemical diffusion, membranes and cellular compartments.  Manage both "normal" concentrations and extremely-low ones (such as macro-molecules that are few in numbers.)

INTERACTIVE VISUALIZATION:  plots, graphs, heatmaps, etc, with interactive adjustable controls.  Extremely personalizable to deal with things like membranes and compartments.

UNIT TESTING: "an ounce of prevention is worth a pound of cure"

MODULARITY: being very disciplined in tackling a large software and data science project.

INFRASTRUCTURE / PLATFORM: a tight alliance between the tools of large-scale software engineering (such as IDE's), the tools of data science (such as JupyterLab) and of data engineering (such as Neo4j graph database.)

MULTIPLE AUDIENCES: address, in the platform and in the documentation, the background and needs of different classes of people, such as programmers, data scientists, chemists, biologists.

 

EXAMPLE.  In the modeling category, take variable spacial and temporal resolutions.  Let's say we have several chemical species that are in near-equilibrium, both in terms of diffusion and reactivity...  and then we have other species that - in some locations and some time periods - have highly variable concentration rates, perhaps because they're produced or consumed in reactions...

Do we really want to waste a lot of high resolution and computing power on the species, locations, time periods that are near-equilibrium?  Conversely, do we want to only coarsely simulate the highly dynamic species/locations/time periods?

Biology informs us that there's quite a range of time scales: for example, impulse firing in neurons is of the order of milliseconds...  while DNA replication is of the order of a hour.  Many orders of magnitude!

How do we best model variable spacial and temporal resolutions for some of the chemical species?  Well, Life123 is a great environment to tackle those design decisions, without being immediately slammed with the intricacies of full 3D!

Fundamental Goals 

  1. Detailed, quantitative biological simulations, including whole prokaryotic cells (bacteria), and later eukaryotic cells

  2. Deeper quantitative insight into human tissue/organ/system physiology, for the advancement of medicine

  3. A very integrative approach that is ultimate conducive to body-wide insights, with an eye to Longevity Science

  4. Explore the minimalist essence of life-like dynamical systems, including their evolution under "genetic algorithms" and other machine-learning approaches.  Also, explore chaotic states

  5. Investigate potential paths for the emergence of life on Earth and on Exoplanets

  6.  A community effort bringing together biologists, system biologists, programmers, machine-learning specialists, biochemists, power-computing engineers, doctors, chemical engineers, data scientists, visualization/UX experts, members of the public & institutions willing to share computing resources, etc.

What Life123 is *NOT*

  • A tool primarily for educational purposes
  • A computer game
  • An art project
  • A Reaction-Diffusion exercise as an end to itself
  • Molecular Dynamics (the reactions rates are assumed already known, or at least with interim estimated values)
  • Detailed modeling of the biophysics of membranes, etc.
  • A mathematical exercise in dynamical systems as an end to itself
  • Wolfram-style "cellular automata"
  • An exploration of hypothetical Physics in one or two dimensions

Note that "as an end to itself" is the operative phrase here; some of those categories do overlap with Life123

Broad Strategies ("Guiding Principles")

  • Aim big, for a simulation scale that may be impractical at the present... but attain a ready state to pounce on the latest advances in computing capabilities - in particular GPU computing, and possibly the upcoming Photonics Computing GPU accelerators

  • Not attempting to create a description and simulation of what happens at the molecular, or near-molecular level.  It's going to be much a much-coarser model that's still accurate enough to capture the essence of the cell's (or system's) behavior: for example, how it respond to inputs and what outputs it produces - what it absorbs, what it secretes, what it does with its internal state, including replication, etc.

    Excessive simulation detail with an ultra-fine spatial/temporal grid, could result in a possibly-excessive information content about the cell (or system): "excessive", in the sense that it may be beyond the information content required to describe the system's "computational capability." (broadly defined; see discussion on dynamical systems and theoretical aspects )

  • A mix of "bottom up" and "top down" approaches:

    "bottom up", as in starting with a minimalist scenarios with fictional molecules and rates, even in 2D or 1D, and gradually advancing to full 3D, real molecules and plausible diffusion/reactions rates, concentrations, etc...

    "top down", as in always setting the target on the final goal of realistic biological systems.  That will be the internal compass always guiding this project.

What does pseudo-1D/-2D mean?

PSEUDO is the operative word!  A true 1D or 2D world would be a speculative exercise in the Foundations of Physics...  In 1D, it'd be molecules on a string...  Would they be able to "pass" each other?  Fermions that can't occupy the same space?  Quantum tunneling to "slide over each other"?  How would force transmit across 1D?  (In 3D, force fields dilute with the square of distance... in 2D, with the distance... but in 1D?)
 
NO!  Life123 isn't about any of that!  If any physicist wants to chat on the side on fascinating concepts of actual 1D or 2D, by all means - but that's completely outside the scope and goals of this project.

Let's think of 1D more like a very thin, long tube of aqueous solution - minus the tube!  In particular, for example, no capillary effects!  Likewise, in 2D, think of the Petri-dish minus the actual dish!

How could 1D or 2D simulations ever be realistic?

They can't - and aren't meant to.  Remember, simulations in 1D and 2D is all about setting up good practices in computing, data modeling, chemical and biological modeling, as a crutch to fully 3D- simulations.

Computing 

  • Python
  • NumPy
  • JupyterLab
  • Custom visualizations with D3.js and Vue.js (a great alliance I discuss here)
  • Network visualizations with cytoscape.js
  • GPU acceleration
  • Multi-core computing, perhaps making use of Dask
  • Distributed computing
  • Neo4j Graph databases (a lot of full-stack infrastructure for that is available thru the open-source sister project Brain Annex)
  • TensorFlow machine learning
  • Unit testing with PyTest

Chemistry 

  • Chemical diffusion
  • Chemical reaction rates
  • Coupled reactions, together with diffusion 
  • Diffusion across membranes (passive and active transport)
  • Temperature effects

Biology

  • Cellular compartments
  • Macro molecules
  • Transcription/Translation
  • Replication

Isn't this overly ambitious?

To paraphrase former president J. F. Kennedy's famous speech about going to the Moon:

We choose to purse quantitative dynamical modeling of complex biological systems not because they are easy, but because they are hard;  because that challenge is one that we are willing to accept, one we are unwilling to postpone, and one we intend to win!


Where does the project currently stand?

[8/2022 update in box below!]  

As of today, with the propitious date (in U.S. format) of  3.14, I am officially releasing to open source the early Beta version of the Life123 platform: GitHub repository

Accompanying the software at this relatively early stage is:

  • This project's "manifesto"
  • A "call to arms" to the community, to be possibly followed by a Discord channel to coordinate. For now, ALL DISCUSSION will take place here
  • A new website: https://life123.science
This project originates from years of discussions I've had with colleagues in the pharmaceutical industry, researchers, and academicians...  as well as papers such as A Whole-Cell Computational Model Predicts Phenotype from Genotype...  and from my own professional background in computing, machine learning, molecular biology, bioinformatics and systems biology.  
 
The current platform stems from early research I carried out in 2020, and active development of the alpha version in late 2021 thru the present.

The underlying philosophical framework is detailed in this entry I wrote in 2019.

August 2022 UPDATE: the platform (repo) is shaping up nicely, and entering a late Beta stage.

In particular, the infrastructure to create interactive visualizations, possibly custom ones (with plotly, Vue.js and D3.js) is getting more polished and streamlined. 
The cytoscape.js library has been brought in for network visualizations.

The website is much more substantial, and existing notebooks and their visualizations are just a click away.
 
A live demo of all JupyterLab notebooks is also one click away!  
 
Ways for you to get involved are getting spelled out in more detail

 
Micro-blogging, ready to later turn into an early discussion platform, was added.

With the recently-released Beta 12 version (change-log), Life123 has by now established a substantial amount of the infrastructures for data science, software management and visualization - for diffusion and (single binding site) reactions.

The really interesting stuff is now just starting!  Upcoming releases will be introducing elements that have been on the table from the very inception of the project:

- membranes
- macro-molecules (with multiple binding sites)
- passive/active transport across membranes

How to scale up?

This is a community project meant to bring together a wide variety of skill sets

In the early fall,  as the foundations get more solid, I'll start publicizing the project, and actively seeking collaborators.

Down the line, some non-profits and companies may opt to get involved.  Perhaps philanthropists and/or investors as well.

A Call to the Community

This is an open call to researchers, academicians, computer scientists, students, colleagues in the pharmaceutical industry, doctors, philanthropists, members of the public.

In the Longevity-Science community, which I've been active in - and at times working in - for a number of years, I hear a lot of "how can I help?"  Well, here's a way!

Do you have skills in:

  • Biology
  • Python programming
  • CUDA programming, or other ways to utilize GPU computing
  • TensorFlow and/or other machine learning 
  • Plotly or D3.js visualization
  • UX 
  • Chemical engineering 
  • Biochemistry
  • Biophysics
  • Systems biology
  • Bioinformatics
  • Medicine
  • Web design
  • Technical writing
  • QA / DevOps
Or do you have access to:
  • Computing resources - such as your own gaming PC, or your company's computing (if you're authorized)
  • Funds for research projects, philanthropy or investments
  • Ways to spread the word (social media, etc)
We need you in all of the above scenarios, and more!  Help us design, implement, test, refine and run the simulations and the platform to scale.
 
I expect to later start a Discord channel.  For now, all discussions (as well as micro-blogging) will take place here.  Also, you're welcome to reach me on LinkedIn or on my professional Facebook account.

The website has a page detailing areas of potential involvement for people with particular skillsets.

Comments

Popular posts from this blog

Discussing Neuroscience with ChatGPT

UPDATED Apr. 2023 - I'm excited by ChatGPT 's possibilities in terms of facilitating advanced learning .  For example, I got enlightening answers to questions that I had confronted when I first studied neuroscience.  The examples below are taken from a very recent session I had with ChatGPT (mid Jan. 2023.) Source: https://neurosciencestuff.tumblr.com In case you're not familiar with ChatGPT, it's a very sophisticated "chatbot" - though, if you call it that way, it'll correct you!  'I am not a "chatbot", I am a language model, a sophisticated type of AI algorithm trained on vast amounts of text data to generate human-like text'. For a high-level explanation of how ChatGPT actually works - which also gives immense insight into its weaknesses, there's an excellent late Jan. 2023 talk by Stephen Wolfram, the brilliant author of the Mathematica software and of Wolfram Alpha , a product that could be combined with ChatGPT to imp...

Using Schema in Graph Databases such as Neo4j

UPDATED Feb. 2024 - Graph databases have an easygoing laissez-faire attitude: "express yourself (almost) however you want"... By contrast, relational databases come across with an attitude like a micro-manager:  "my way or the highway"... Is there a way to take the best of both worlds and distance oneself from their respective excesses, as best suited for one's needs?  A way to marry the flexibility of Graph Databases and the discipline of Relational Databases? This article is part 5 of a growing,  ongoing  series  on Graph Databases and Neo4j Let's Get Concrete Consider a simple scenario with scientific data such as the Sample, Experiment, Study, Run Result , where Samples are used in Experiments, and where Experiments are part of Studies and produce Run Results.  That’s all very easy and intuitive to represent and store in a Labeled Graph Database such as Neo4j .   For example, a rough draft might go like this:   The “labels” (b...

Graph Databases (Neo4j) - a revolution in modeling the real world!

UPDATED Oct. 2023 - I was "married" to Relational Databases for many years... and it was a good "relationship" full of love and productivity - but SOMETHING WAS MISSING! Let me backtrack.   In college, I got a hint of the "pre-relational database" days...  Mercifully, that was largely before my time, but  - primarily through a class - I got a taste of what the world was like before relational databases.  It's an understatement to say: YUCK! Gratitude for the power and convenience of Relational Databases and SQL - and relief at having narrowly averted life before it! - made me an instant mega-fan of that technology.  And for many years I held various jobs that, directly or indirectly, made use of MySQL and other relational databases - whether as a Database Administrator, Full-Stack Developer, Data Scientist, CTO or various other roles. UPDATE: This article is now part 1 of a growing, ongoing series on Graph Databases and Neo4j But ther...

Anti-Aging Research: Science, not Hype

Last updated May 2023 Q: "How is aging a disease?" A: It's a dynamic system that veers away from its homeostasis (normal equilibrium point): hence a form of slow-progressing illness. Labeling it as 'natural' is a surrender to our traditional state of ignorance and powerlessness, which fortunately is beginning to be changed! Aging is "normal" only from the point of view of the "selfish gene", for whom the body is a disposable carrier. Individuals organisms - for whom self-preservation has a different meaning than for genes - have received scant help from evolution... with rare exceptions such as the T. dohrnii jellyfish (which I discuss here )... but now the time has finally arrived for our rational design to remedy some of the cellular flaws that evolution never bothered to correct!   The above is my standard answer to an oft-asked question. The science of aging is by all evidence very misunderstood by the general public.  Hype,...

What are Graph Databases - and Why Should I Care?? : "Graph Databases for Poets"

  This is a very gentle introduction to the subject.  The subtitle is inspired by university courses such as "Physics for Poets"!  (if you're technically inclined, there's an alternate article for you.) It has been said that "The language of physics (or of God) is math".  On a similar note, it could be said that: The language of the biological world - or of any subject or endeavor involving complexity - is networks ('meshes') What is a network?  Think of  it as the familiar 'friends of friends' diagram from social media. Everywhere one turns in biology, there's a network – at the cellular level, tissue level, organ level, ecosystem level.  The weather and other earth systems are networks.  Human societal organization is a network.  Electrical circuits, the Internet, our own brains...  Networks are everywhere! What can we do with networks, to better understand the world around us, or to create something that we need? Broadly ...

Using Neo4j with Python : the Open-Source Library "NeoAccess"

So, you want to build a python app or Jupyter notebook to utilize Neo4j, but aren't too keen on coding a lot of string manipulation to programmatic create ad-hoc Cypher queries?   You're in the right place: the NeoAccess library can do take care of all that, sparing you from lengthy, error-prone development that requires substantial graph-database and software-development expertise! This article is part 4 of a growing,  ongoing  series  on Graph Databases and Neo4j   "NeoAccess" is the bottom layer of the technology stack provided by the BrainAnnex open-source project .  All layers are very modular, and the NeoAccess library may also be used by itself , entirely separately from the rest of the technology stack.  (A diagram of the full stack is shown later in this article.) NeoAccess interacts with the Neo4j Python driver , which is provided by the Neo4j company, to access the database from Python; the API to access that driver is very p...

Visualization of Graph Databases Using Cytoscape.js

(UPDATED APR. 2024)   I have ample evidence from multiple sources that there are strong unmet needs in the area of visualization of graph databases. And whenever there's a vacuum, vendors circle like vultures - with incomplete, non-customizable, and at times ridiculously expensive, closed-box proprietary solutions.   Fortunately, coming to the rescue is the awesome open-source cytoscape.js library ,  an offshoot of the "Cytoscape" project of the  Institute for Systems Biology , a project with a long history that goes back to 2002. One can do amazing custom solutions, relatively easily, when one combines this Cytoscape library with:   1) a front-end framework such as Vue.js   2) backend libraries (for example in python) to prepare and serve the data   For example, a while back I created a visualizer for networks of chemical reactions, for another open-source project I lead ( life123.science )   This visualizer will look and feel generally ...

Neo4j Sandbox Tutorial : try Neo4j and learn Cypher - free and easy!

So, you have an itch to test-drive Neo4j and its Cypher query language.  Maybe you want to learn it, or evaluate it, or introduce colleagues/clients to it.  And you wish for: fast, simple and free! Well, good news: the Neo4j company kindly provides a free, short-term hosted solution called "the Neo4j sandbox" .  Extremely easy to set up and use! This article is part 2 of a growing, ongoing series on Graph Databases and Neo4j Register (free) for the Neo4j "Sandbox" Go to sandbox.neo4j.com , and register with a working email and a password.  That's it! Note that this same email/password will also let you into the Neo4j Community Forums and Support ; the same login for all: very convenient! Launch your instance - blank or pre-populated After registering, go to  sandbox.neo4j.com  , and follow the steps in the diagram below (the choices might differ, but the "Blank Sandbox" should always be there): Too good to be true?  Is there...

Full-Text Search with the Neo4j Graph Database

(UPDATED May 2024)   Now that we have discussed a full technology stack based on Neo4j (or other graph databases), and that we a design and implementation available from the open-source project BrainAnnex.org  , what next?  What shall we build on top? Well, how about  Full-Text Search ?  This article is part of a growing, ongoing series on Graph Databases and Neo4j Full-Text Searching/Indexing The Brain Annex open-source project includes an implementation of a design that uses the convenient services of its Schema Layer , to provide indexing of word-based documents using Neo4j. The python class FullTextIndexing ( source code ) provides the necessary methods, and it can parse both plain-text and HTML documents (for example, used in "formatted notes"); parsing of PDF files and other formats will be added at a later date. No grammatical analysis ( stemming or lemmatizing ) is done on the text.  However, a long list of common word ("stop ...

Neo4j & Cypher Tutorial : Getting Started with a Graph Database and its Query Language

You have a general idea of what Graph Databases - and Neo4j in particular - are...  But how to get started?  Read on! This article is part 3 of a growing,  ongoing  series  on Graph Databases and Neo4j   If you're new to graph databases, please check out part 1 for an intro and motivation about them.  There, we discussed an example about an extremely simple database involving actors, movies and directors...  and saw how easy the Cypher query language makes it to answer questions such as "which directors have worked with Tom Hanks in 2016" - questions that, when done with relational databases and SQL, turn into a monster of a query and an overly-complicated data model involving a whopping 5 tables! In this tutorial, we will actually carry out that query - and get acquainted with Cypher and the Neo4j browser interface in the process.  This is the dataset we'll be constructing: Get the database in place If you don't already have a datab...