ISC Keynote: The Algorithms of Life – Scientific Computing for Systems Biology

By John Russell

June 19, 2019

Systems biology has existed loosely under many definitions for a couple of decades. It’s the notion of describing living systems using first-principle physics and mathematics to capture life in equations that are both descriptive and predictive – and let’s add productive by which we mean being able to deliver therapies (drugs et. al) to enhance health and fight disease.

Doing that has proven difficult at best and disappointing at worst as even a cursory glance at the state of healthcare reveals; that’s notwithstanding many marvelous breakthroughs such as sequencing the human genome and the steady chipping away at functional genomics (and other ‘omics) to understand better how DNA informs what we become.

Ivo Sbalzarini

With apologies to ISC organizers I’ve stolen the name of the opening keynote by Ivo Sbalzarini –  The Algorithms of Life – Scientific Computing for Systems Biology – for the headline of this article in an attempt to capture his expansive presentation. Thanks also to Sbalzarini for providing a few of his slides.

Given all we know today and the steady gush of experimental data from modern instruments, what we are missing, said Sbalzarini, are the algorithms to make sense of it all. Having poked away at this problem for nearly as long as it has been around, Sbalzarini presented a sweeping approach to digging out those algorithms by capitalizing on recent advances in imaging technology, immersive virtual/augmented reality, a sophisticated analysis approach that leverages particle-mesh mathematics and which has been built into a software platform (OpenFPM), and lastly, no surprise, the steadily growing power of HPC.

 

As in many important life sciences advances the ‘lowly’ fruit fly took center stage. In this instance the analysis was to investigate a dysregulation in embryogenesis – specifically the failure of tissue to fold properly. In the end, the researchers identified the DNA influence, the chemical environment influence, and the mechanical environment influence, and delivered a predictive understanding of the embryo’s tissue response. Lest you think this is old work, it was presented last week at the New York Scientific Data Summit.

Getting from Sbalzarini’s early nascent research 15 years ago to the impressive results (and tool suite) presented is a long journey. We’ll summarize as practical but the ISC is likely to archive its keynote; for biologists it is well worth watching.

Advanced imaging, such as light sheet microscopy, now makes it possible to observe life science phenomena in 3D and great detail at the cellular and intracellular level.

“We can image an embryo from the time it is a fertilized to the time it moves out of the microscope field by itself and continues its life. When we image the fruit fly embryo over the 72 hours of development, we gather 180 TB of image data. If you would like to visualize that in real-time. That means a rendering performance or a rendering throughput in real-time of about 1.8 Gigapixels per second,” said Sbalzarini[i]. A key advantage here is the animal stays alive unlike older approaches requiring stains and fixing.

Hardly just pretty pictures, the extensive image data captured (and visualizations possible) are the raw input for building hypotheses and predictive models. The other primary driver is Sbalzarini clever adaption of particle-mesh technology to convert the data into actionable, in silico simulation. Underlying HPC infrastructure, of course, is the engine without which the whole process would grind to a halt.

“The numerical methods are particle methods or hybrid particle mesh methods. They comprise an interesting class of numerical methods. They discretize the system by particles, so if you have a complex geometry, you don’t need to generate the mesh for the simulation, but you simply fill the geometry with particles that store the variables; there can be a mesh in addition in order to do far field equations in order to compute for example forces for far field equations, for example,” he said.

“This is a classic framework of particle-mesh methods to solve partial differential equations, but particle methods as an algorithm are much more general than that. I would define everything as a particle method that is composed of dots of zero dimension elements that are characterized by a position in some space and some properties that they carry. Such an algorithm can be used to solve partial differential equations where the particles are the colocation points of your various discretization and they store the values of the field at that position.”

He adds quickly, “There is nothing that limits us to having particles interacting in a deterministic fashion and this then also allows us to solve stochastic different equations, numerically or to perform agent based simulation or agent based modeling.”

Building the computational tools to deliver these models has been a challenging and lengthy task for which Sbalzarini is well-qualified. He is the chair of scientific computing for systems biology on the faculty of computer science of TU Dresden, as well as the faculty of mathematics, and director of the TUD-Department in the Center for Systems Biology Dresden. He also is a permanent Senior Research Group Leader with the Max Planck Institute of Molecular Cell Biology and Genetics in Dresden.

Leaving out many details and with regrets for over-simplification, Sbalzarini and colleagues imaged the fruit fly embryo; used machine learning to identify ‘algorithms’, converted the data and algorithms into models based on particle-mesh approaches using their home-developed platform; ran computational experiments to test their hypotheses; and used immersive visualization technology as a step to allow researchers to see the real process and simulations unfold. “It is possible to walk around inside the simulation,” he said. Informed by what they saw and their knowledge, researchers tweaked parameters and hypotheses, iteratively converging on a solution.

“To me it is a very nice example of how HPC and these numerically intricate simulations that we can do with these machines allow us to bridge really from the molecular scale to the tissue scaler in order to explain how things work and in order to propose remedies,” said Sbalzarini.

Sbalzarini reminded the audience living systems are computing machines themselves, “[A fruit fly embryo] is a massively parallel and fully self-organized system in which we can view every single cell as a processing element that executes programs. [It’s a] highly interconnected computer and able to solve NP hard problems with billions or hundreds of billions of processing elements. We know a lot about the hardware of this computer – the proteins, the molecules, the lipids, the fats out of which this computer is made – and thanks to sequencing technology, [we’re] able to read the source code of this computer, which is the genomic sequence. However we have no idea what algorithms this source code implements on his hardware.”

Now, advanced imaging and machine learning capabilities are catalyzing researchers’ ability to identify ‘mechanistic’ guidelines and incorporate traditional formulations (ODEs/PDEs) of physics laws and mathematics into the life sciences tool box. Chemical diffusion. Fluid dynamics. EMI influences. Activation energy thresholds. These are the kinds of attributes that can be captured in particle-mesh models.

When Sbalzarini began his studies in earnest, he used an NEC SX-5 with 512 processors housed at CSCS (Swiss Supercomputer Center). In 2005 that became a Cray XT-3 with 1664 processors. A lot has changed since. The first iteration of the system biology software platform his team developed was Parallel Particle Mesh Library (PPM) written in Fortran 90 many years ago.  It served as layer between MPI and Client Applications for simulations of physical systems using Particle-Mesh methods. The PPM library runs on single and multi-processor architectures, and handles 2D and 3D problems.

“The PPM library had two parts, what we call the PPM core, which is implemented in all the communication primitives, the load balancing, the file IO, [and] the distributed data structures. And the PPM numerics using frequently used numerical solvers; it does this in part by using the abstractions from the core and in part by renting third party libraries such as PETSc or FFTW. On top of PPM there is a domain specific programming language called PPM Language which provides a reasonably simple way of coding PPM but you could also directly interface with the Fortran API,”

PPML used overloading and generic interfaces and provided for the limitations of the important routines for different hardware platforms such as vector processors, like the NEC system, shared memory, distributed memory, even single processor systems, said Sbalzarini.

It was a beast to maintain. “Because of overloading the amount of source code in the PPM library was huge, several millions of lines of code that needed to be maintained here and ported. What we liked about PPML was the abstraction on which it is based. It’s a set of abstract data types and abstract operators for computing that are in our opinion the most coarse-grained abstractions possible that still cleanly separate computation from communication. So in PPM an abstraction would either only compute but not incur any communication overhead or it would only communicate but not do any computation,” he said.

Five years ago the platform was upgraded, “We decided to keep the abstractions, to keep the definitions of the data types and the operators, but now implement a C++ library which is called OpenFPM (Open Framework for Particle Method Library) and make use of template metaprogrammingin C++ for compiled time code generation. OpenFPM can do much more than PPM, for example it can do simulations in arbitrary dimensional spaces where PPM is limited to 2D and 3D. OpenFPM allowed particle properties to be objects of any C++ that the user can define and all the communication and file IO will work for it,” he said.

Adopting template metaprogramming reduced the amount of code needed to “about a factor of ten less complexity than the PPM.”

Sbalzarini presented many more details in his rich talk. It will be interesting to watch how widely OpenFPM is used and if it gains tractions in other domains. Ease of use is a key question for many biomedical researchers and clinicians. Sbalzarini said, “This hopefully makes HPC so easy to use that every science-based application in biology, in computational biology, and also in other fields can benefit.”

That said computer expertise, particularly HPC expertise, has historically been lacking in life science although that is changing and fairly quickly.

The main motivation is to understand biology and to understand how cells form tissues, and eventually to be able to provide novel explanations for disease phenotypes and maybe therapies for disease, said Sbalzarini. Nevertheless, “For us as computer scientists it’s also just a lot of fun because what we do combines several technologies that we think are fun to work with, technologies like virtual reality, HPC, massively scalable software systems, building microscopes and playing with optics, or using and developing artificial intelligence and learning algorithms to interface with the living things in the microscope.”

[i]Some quotes have been very lightly edited to improve readability.

 

Subscribe to HPCwire's Weekly Update!

Be the most informed person in the room! Stay ahead of the tech trends with industy updates delivered to you every week!

The Case for an Edge-Driven Future for Supercomputing

September 24, 2021

“Exascale only becomes valuable when it’s creating and using data that we care about,” said Pete Beckman, co-director of the Northwestern-Argonne Institute of Science and Engineering (NAISE), at the most recent HPC Read more…

Three Universities Team for NSF-Funded ‘ACES’ Reconfigurable Supercomputer Prototype

September 23, 2021

As Moore’s law slows, HPC developers are increasingly looking for speed gains in specialized code and specialized hardware – but this specialization, in turn, can make testing and deploying code trickier than ever. Now, researchers from Texas A&M University, the University of Illinois at Urbana... Read more…

Qubit Stream: Monte Carlo Advance, Infosys Joins the Fray, D-Wave Meeting Plans, and More

September 23, 2021

It seems the stream of quantum computing reports never ceases. This week – IonQ and Goldman Sachs tackle Monte Carlo on quantum hardware, Cambridge Quantum pushes chemistry calculations forward, D-Wave prepares for its Read more…

Asetek Announces It Is Exiting HPC to Protect Future Profitability

September 22, 2021

Liquid cooling specialist Asetek, well-known in HPC circles for its direct-to-chip cooling technology that is inside some of the fastest supercomputers in the world, announced today that it is exiting the HPC space amid multiple supply chain issues related to the pandemic. Although pandemic supply chain... Read more…

TACC Supercomputer Delves Into Protein Interactions

September 22, 2021

Adenosine triphosphate (ATP) is a compound used to funnel energy from mitochondria to other parts of the cell, enabling energy-driven functions like muscle contractions. For ATP to flow, though, the interaction between the hexokinase-II (HKII) enzyme and the proteins found in a specific channel on the mitochondria’s outer membrane. Now, simulations conducted on supercomputers at the Texas Advanced Computing Center (TACC) have simulated... Read more…

AWS Solution Channel

Introducing AWS ParallelCluster 3

Running HPC workloads, like computational fluid dynamics (CFD), molecular dynamics, or weather forecasting typically involves a lot of moving parts. You need a hundreds or thousands of compute cores, a job scheduler for keeping them fed, a shared file system that’s tuned for throughput or IOPS (or both), loads of libraries, a fast network, and a head node to make sense of all this. Read more…

The Latest MLPerf Inference Results: Nvidia GPUs Hold Sway but Here Come CPUs and Intel

September 22, 2021

The latest round of MLPerf inference benchmark (v 1.1) results was released today and Nvidia again dominated, sweeping the top spots in the closed (apples-to-apples) datacenter and edge categories. Perhaps more interesti Read more…

The Case for an Edge-Driven Future for Supercomputing

September 24, 2021

“Exascale only becomes valuable when it’s creating and using data that we care about,” said Pete Beckman, co-director of the Northwestern-Argonne Institut Read more…

Three Universities Team for NSF-Funded ‘ACES’ Reconfigurable Supercomputer Prototype

September 23, 2021

As Moore’s law slows, HPC developers are increasingly looking for speed gains in specialized code and specialized hardware – but this specialization, in turn, can make testing and deploying code trickier than ever. Now, researchers from Texas A&M University, the University of Illinois at Urbana... Read more…

Qubit Stream: Monte Carlo Advance, Infosys Joins the Fray, D-Wave Meeting Plans, and More

September 23, 2021

It seems the stream of quantum computing reports never ceases. This week – IonQ and Goldman Sachs tackle Monte Carlo on quantum hardware, Cambridge Quantum pu Read more…

Asetek Announces It Is Exiting HPC to Protect Future Profitability

September 22, 2021

Liquid cooling specialist Asetek, well-known in HPC circles for its direct-to-chip cooling technology that is inside some of the fastest supercomputers in the world, announced today that it is exiting the HPC space amid multiple supply chain issues related to the pandemic. Although pandemic supply chain... Read more…

TACC Supercomputer Delves Into Protein Interactions

September 22, 2021

Adenosine triphosphate (ATP) is a compound used to funnel energy from mitochondria to other parts of the cell, enabling energy-driven functions like muscle contractions. For ATP to flow, though, the interaction between the hexokinase-II (HKII) enzyme and the proteins found in a specific channel on the mitochondria’s outer membrane. Now, simulations conducted on supercomputers at the Texas Advanced Computing Center (TACC) have simulated... Read more…

The Latest MLPerf Inference Results: Nvidia GPUs Hold Sway but Here Come CPUs and Intel

September 22, 2021

The latest round of MLPerf inference benchmark (v 1.1) results was released today and Nvidia again dominated, sweeping the top spots in the closed (apples-to-ap Read more…

Why HPC Storage Matters More Now Than Ever: Analyst Q&A

September 17, 2021

With soaring data volumes and insatiable computing driving nearly every facet of economic, social and scientific progress, data storage is seizing the spotlight. Hyperion Research analyst and noted storage expert Mark Nossokoff looks at key storage trends in the context of the evolving HPC (and AI) landscape... Read more…

GigaIO Gets $14.7M in Series B Funding to Expand Its Composable Fabric Technology to Customers

September 16, 2021

Just before the COVID-19 pandemic began in March 2020, GigaIO introduced its Universal Composable Fabric technology, which allows enterprises to bring together Read more…

Ahead of ‘Dojo,’ Tesla Reveals Its Massive Precursor Supercomputer

June 22, 2021

In spring 2019, Tesla made cryptic reference to a project called Dojo, a “super-powerful training computer” for video data processing. Then, in summer 2020, Tesla CEO Elon Musk tweeted: “Tesla is developing a [neural network] training computer called Dojo to process truly vast amounts of video data. It’s a beast! … A truly useful exaflop at de facto FP32.” Read more…

Enter Dojo: Tesla Reveals Design for Modular Supercomputer & D1 Chip

August 20, 2021

Two months ago, Tesla revealed a massive GPU cluster that it said was “roughly the number five supercomputer in the world,” and which was just a precursor to Tesla’s real supercomputing moonshot: the long-rumored, little-detailed Dojo system. “We’ve been scaling our neural network training compute dramatically over the last few years,” said Milan Kovac, Tesla’s director of autopilot engineering. Read more…

Esperanto, Silicon in Hand, Champions the Efficiency of Its 1,092-Core RISC-V Chip

August 27, 2021

Esperanto Technologies made waves last December when it announced ET-SoC-1, a new RISC-V-based chip aimed at machine learning that packed nearly 1,100 cores onto a package small enough to fit six times over on a single PCIe card. Now, Esperanto is back, silicon in-hand and taking aim... Read more…

CentOS Replacement Rocky Linux Is Now in GA and Under Independent Control

June 21, 2021

The Rocky Enterprise Software Foundation (RESF) is announcing the general availability of Rocky Linux, release 8.4, designed as a drop-in replacement for the soon-to-be discontinued CentOS. The GA release is launching six-and-a-half months after Red Hat deprecated its support for the widely popular, free CentOS server operating system. The Rocky Linux development effort... Read more…

Intel Completes LLVM Adoption; Will End Updates to Classic C/C++ Compilers in Future

August 10, 2021

Intel reported in a blog this week that its adoption of the open source LLVM architecture for Intel’s C/C++ compiler is complete. The transition is part of In Read more…

Hot Chips: Here Come the DPUs and IPUs from Arm, Nvidia and Intel

August 25, 2021

The emergence of data processing units (DPU) and infrastructure processing units (IPU) as potentially important pieces in cloud and datacenter architectures was Read more…

AMD-Xilinx Deal Gains UK, EU Approvals — China’s Decision Still Pending

July 1, 2021

AMD’s planned acquisition of FPGA maker Xilinx is now in the hands of Chinese regulators after needed antitrust approvals for the $35 billion deal were receiv Read more…

Google Launches TPU v4 AI Chips

May 20, 2021

Google CEO Sundar Pichai spoke for only one minute and 42 seconds about the company’s latest TPU v4 Tensor Processing Units during his keynote at the Google I Read more…

Leading Solution Providers

Contributors

HPE Wins $2B GreenLake HPC-as-a-Service Deal with NSA

September 1, 2021

In the heated, oft-contentious, government IT space, HPE has won a massive $2 billion contract to provide HPC and AI services to the United States’ National Security Agency (NSA). Following on the heels of the now-canceled $10 billion JEDI contract (reissued as JWCC) and a $10 billion... Read more…

10nm, 7nm, 5nm…. Should the Chip Nanometer Metric Be Replaced?

June 1, 2020

The biggest cool factor in server chips is the nanometer. AMD beating Intel to a CPU built on a 7nm process node* – with 5nm and 3nm on the way – has been i Read more…

Julia Update: Adoption Keeps Climbing; Is It a Python Challenger?

January 13, 2021

The rapid adoption of Julia, the open source, high level programing language with roots at MIT, shows no sign of slowing according to data from Julialang.org. I Read more…

Quantum Roundup: IBM, Rigetti, Phasecraft, Oxford QC, China, and More

July 13, 2021

IBM yesterday announced a proof for a quantum ML algorithm. A week ago, it unveiled a new topology for its quantum processors. Last Friday, the Technical Univer Read more…

Intel Launches 10nm ‘Ice Lake’ Datacenter CPU with Up to 40 Cores

April 6, 2021

The wait is over. Today Intel officially launched its 10nm datacenter CPU, the third-generation Intel Xeon Scalable processor, codenamed Ice Lake. With up to 40 Read more…

Frontier to Meet 20MW Exascale Power Target Set by DARPA in 2008

July 14, 2021

After more than a decade of planning, the United States’ first exascale computer, Frontier, is set to arrive at Oak Ridge National Laboratory (ORNL) later this year. Crossing this “1,000x” horizon required overcoming four major challenges: power demand, reliability, extreme parallelism and data movement. Read more…

Intel Unveils New Node Names; Sapphire Rapids Is Now an ‘Intel 7’ CPU

July 27, 2021

What's a preeminent chip company to do when its process node technology lags the competition by (roughly) one generation, but outmoded naming conventions make it seem like it's two nodes behind? For Intel, the response was to change how it refers to its nodes with the aim of better reflecting its positioning within the leadership semiconductor manufacturing space. Intel revealed its new node nomenclature, and... Read more…

Latest MLPerf Results: Nvidia Shines but Intel, Graphcore, Google Increase Their Presence

June 30, 2021

While Nvidia (again) dominated the latest round of MLPerf training benchmark results, the range of participants expanded. Notably, Google’s forthcoming TPU v4 Read more…

  • arrow
  • Click Here for More Headlines
  • arrow
HPCwire