An adaptive stepsize method for the chemical Langevin equation

The Journal of Chemical Physics
Silvana Ilie, Alexandra Teslya

Abstract

Mathematical and computational modeling are key tools in analyzing important biological processes in cells and living organisms. In particular, stochastic models are essential to accurately describe the cellular dynamics, when the assumption of the thermodynamic limit can no longer be applied. However, stochastic models are computationally much more challenging than the traditional deterministic models. Moreover, many biochemical systems arising in applications have multiple time-scales, which lead to mathematical stiffness. In this paper we investigate the numerical solution of a stochastic continuous model of well-stirred biochemical systems, the chemical Langevin equation. The chemical Langevin equation is a stochastic differential equation with multiplicative, non-commutative noise. We propose an adaptive stepsize algorithm for approximating the solution of models of biochemical systems in the Langevin regime, with small noise, based on estimates of the local error. The underlying numerical method is the Milstein scheme. The proposed adaptive method is tested on several examples arising in applications and it is shown to have improved efficiency and accuracy compared to the existing fixed stepsize schemes.

References

Jul 5, 2001·Proceedings of the National Academy of Sciences of the United States of America·M Thattai, A van Oudenaarden
Mar 2, 2002·Science·Hiroaki Kitano
Aug 17, 2002·Science·Michael B ElowitzPeter S Swain
Nov 15, 2002·Nature·Hiroaki Kitano
Jul 21, 2004·Computational Biology and Chemistry·T E TurnerK Burrage
Nov 20, 2004·The Journal of Chemical Physics·Tianhai Tian, Kevin Burrage
Jan 11, 2005·The Journal of Chemical Physics·Abhijit ChatterjeeMarkos A Katsoulakis
Feb 8, 2006·The Journal of Chemical Physics·Yang CaoLinda R Petzold
May 19, 2006·Journal of Computational Biology : a Journal of Computational Molecular Cell Biology·Andrew Golightly, Darren J Wilkinson
Mar 9, 2007·The Journal of Chemical Physics·Michel F Pettigrew, Haluk Resat
Jun 22, 2007·The Journal of Chemical Physics·Yang CaoLinda R Petzold
Feb 13, 2008·The Journal of Chemical Physics·David F Anderson
Sep 15, 2011·The Journal of Chemical Physics·Simon L CotterRadek Erban

❮ Previous
Next ❯

Citations

Jul 2, 2019·Wiley Interdisciplinary Reviews. Systems Biology and Medicine·Giulia SimoniLuca Marchetti

❮ Previous
Next ❯

Related Concepts

Trending Feeds

COVID-19

Coronaviruses encompass a large family of viruses that cause the common cold as well as more serious diseases, such as the ongoing outbreak of coronavirus disease 2019 (COVID-19; formally known as 2019-nCoV). Coronaviruses can spread from animals to humans; symptoms include fever, cough, shortness of breath, and breathing difficulties; in more severe cases, infection can lead to death. This feed covers recent research on COVID-19.

Blastomycosis

Blastomycosis fungal infections spread through inhaling Blastomyces dermatitidis spores. Discover the latest research on blastomycosis fungal infections here.

Nuclear Pore Complex in ALS/FTD

Alterations in nucleocytoplasmic transport, controlled by the nuclear pore complex, may be involved in the pathomechanism underlying multiple neurodegenerative diseases including Amyotrophic Lateral Sclerosis and Frontotemporal Dementia. Here is the latest research on the nuclear pore complex in ALS and FTD.

Applications of Molecular Barcoding

The concept of molecular barcoding is that each original DNA or RNA molecule is attached to a unique sequence barcode. Sequence reads having different barcodes represent different original molecules, while sequence reads having the same barcode are results of PCR duplication from one original molecule. Discover the latest research on molecular barcoding here.

Chronic Fatigue Syndrome

Chronic fatigue syndrome is a disease characterized by unexplained disabling fatigue; the pathology of which is incompletely understood. Discover the latest research on chronic fatigue syndrome here.

Evolution of Pluripotency

Pluripotency refers to the ability of a cell to develop into three primary germ cell layers of the embryo. This feed focuses on the mechanisms that underlie the evolution of pluripotency. Here is the latest research.

Position Effect Variegation

Position Effect Variagation occurs when a gene is inactivated due to its positioning near heterochromatic regions within a chromosome. Discover the latest research on Position Effect Variagation here.

STING Receptor Agonists

Stimulator of IFN genes (STING) are a group of transmembrane proteins that are involved in the induction of type I interferon that is important in the innate immune response. The stimulation of STING has been an active area of research in the treatment of cancer and infectious diseases. Here is the latest research on STING receptor agonists.

Microbicide

Microbicides are products that can be applied to vaginal or rectal mucosal surfaces with the goal of preventing, or at least significantly reducing, the transmission of sexually transmitted infections. Here is the latest research on microbicides.