A hybrid algorithm for coupling partial differential equation and compartment-based dynamics

Stochastic simulation methods can be applied successfully to model exact spatio-temporally resolved reaction–diffusion systems. However, in many cases, these methods can quickly become extremely computationally intensive with increasing particle numbers. An alternative description of many of these s...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Journal of the Royal Society interface 2016-09, Vol.13 (122), p.20160335
Hauptverfasser: Harrison, Jonathan U., Yates, Christian A.
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
container_end_page
container_issue 122
container_start_page 20160335
container_title Journal of the Royal Society interface
container_volume 13
creator Harrison, Jonathan U.
Yates, Christian A.
description Stochastic simulation methods can be applied successfully to model exact spatio-temporally resolved reaction–diffusion systems. However, in many cases, these methods can quickly become extremely computationally intensive with increasing particle numbers. An alternative description of many of these systems can be derived in the diffusive limit as a deterministic, continuum system of partial differential equations (PDEs). Although the numerical solution of such PDEs is, in general, much more efficient than the full stochastic simulation, the deterministic continuum description is generally not valid when copy numbers are low and stochastic effects dominate. Therefore, to take advantage of the benefits of both of these types of models, each of which may be appropriate in different parts of a spatial domain, we have developed an algorithm that can be used to couple these two types of model together. This hybrid coupling algorithm uses an overlap region between the two modelling regimes. By coupling fluxes at one end of the interface and using a concentration-matching condition at the other end, we ensure that mass is appropriately transferred between PDE- and compartment-based regimes. Our methodology gives notable reductions in simulation time in comparison with using a fully stochastic model, while maintaining the important stochastic features of the system and providing detail in appropriate areas of the domain. We test our hybrid methodology robustly by applying it to several biologically motivated problems including diffusion and morphogen gradient formation. Our analysis shows that the resulting error is small, unbiased and does not grow over time.
doi_str_mv 10.1098/rsif.2016.0335
format Article
fullrecord <record><control><sourceid>proquest_pubme</sourceid><recordid>TN_cdi_pubmed_primary_27628171</recordid><sourceformat>XML</sourceformat><sourcesystem>PC</sourcesystem><sourcerecordid>1820595435</sourcerecordid><originalsourceid>FETCH-LOGICAL-c364t-59b3d62c9abfe8cff34f3e862ad058116620bb1fbc86680c384942960e42c2243</originalsourceid><addsrcrecordid>eNp1kc1LwzAYh4MofkyvHqVHL535XnoRhvgFAy96E0KaJltG28ykFfrfm7o59OApeXmf_PImDwCXCE4RLMRNiM5OMUR8CglhB-AUzSjOGef4cL8XxQk4i3ENIZkRxo7BCZ5xLNAMnYL3ebYayuCqTNVLH1y3ajLrQ6Z9v6ldu8w2KnRO1VnlrDXBtN-F-ehV53ybqbZKaDNCTerlpYqmyqqhVY3T8RwcWVVHc7FbJ-Dt4f717ilfvDw-380XuSacdjkrSlJxrAtVWiO0tYRaYgTHqoJMIJQeA8sS2VILzgXURNCC4oJDQ7HGmJIJuN3mbvqyMZVOkwRVy01wjQqD9MrJv53WreTSf0oGKS8oSQHXu4DgP3oTO9m4qE1dq9b4PkokMGQFo4QldLpFdfAxBmP31yAoRyVyVCJHJXJUkg5c_R5uj_84SADZAsEP6Ze8dqYb5Nr3oU3lf7FfyTCbpA</addsrcrecordid><sourcetype>Open Access Repository</sourcetype><iscdi>true</iscdi><recordtype>article</recordtype><pqid>1820595435</pqid></control><display><type>article</type><title>A hybrid algorithm for coupling partial differential equation and compartment-based dynamics</title><source>MEDLINE</source><source>PubMed Central</source><creator>Harrison, Jonathan U. ; Yates, Christian A.</creator><creatorcontrib>Harrison, Jonathan U. ; Yates, Christian A.</creatorcontrib><description>Stochastic simulation methods can be applied successfully to model exact spatio-temporally resolved reaction–diffusion systems. However, in many cases, these methods can quickly become extremely computationally intensive with increasing particle numbers. An alternative description of many of these systems can be derived in the diffusive limit as a deterministic, continuum system of partial differential equations (PDEs). Although the numerical solution of such PDEs is, in general, much more efficient than the full stochastic simulation, the deterministic continuum description is generally not valid when copy numbers are low and stochastic effects dominate. Therefore, to take advantage of the benefits of both of these types of models, each of which may be appropriate in different parts of a spatial domain, we have developed an algorithm that can be used to couple these two types of model together. This hybrid coupling algorithm uses an overlap region between the two modelling regimes. By coupling fluxes at one end of the interface and using a concentration-matching condition at the other end, we ensure that mass is appropriately transferred between PDE- and compartment-based regimes. Our methodology gives notable reductions in simulation time in comparison with using a fully stochastic model, while maintaining the important stochastic features of the system and providing detail in appropriate areas of the domain. We test our hybrid methodology robustly by applying it to several biologically motivated problems including diffusion and morphogen gradient formation. Our analysis shows that the resulting error is small, unbiased and does not grow over time.</description><identifier>ISSN: 1742-5689</identifier><identifier>EISSN: 1742-5662</identifier><identifier>DOI: 10.1098/rsif.2016.0335</identifier><identifier>PMID: 27628171</identifier><language>eng</language><publisher>England: The Royal Society</publisher><subject>Algorithms ; Deterministic ; Hybrid Algorithms ; Life Sciences–Mathematics interface ; Models, Theoretical ; Multiscale ; Reaction–diffusion ; Stochastic</subject><ispartof>Journal of the Royal Society interface, 2016-09, Vol.13 (122), p.20160335</ispartof><rights>2016 The Authors.</rights><rights>2016 The Authors. 2016</rights><lds50>peer_reviewed</lds50><oa>free_for_read</oa><woscitedreferencessubscribed>false</woscitedreferencessubscribed><citedby>FETCH-LOGICAL-c364t-59b3d62c9abfe8cff34f3e862ad058116620bb1fbc86680c384942960e42c2243</citedby><cites>FETCH-LOGICAL-c364t-59b3d62c9abfe8cff34f3e862ad058116620bb1fbc86680c384942960e42c2243</cites><orcidid>0000-0002-2748-9921</orcidid></display><links><openurl>$$Topenurl_article</openurl><openurlfulltext>$$Topenurlfull_article</openurlfulltext><thumbnail>$$Tsyndetics_thumb_exl</thumbnail><linktopdf>$$Uhttps://www.ncbi.nlm.nih.gov/pmc/articles/PMC5046943/pdf/$$EPDF$$P50$$Gpubmedcentral$$Hfree_for_read</linktopdf><linktohtml>$$Uhttps://www.ncbi.nlm.nih.gov/pmc/articles/PMC5046943/$$EHTML$$P50$$Gpubmedcentral$$Hfree_for_read</linktohtml><link.rule.ids>230,315,728,781,785,886,27926,27927,53793,53795</link.rule.ids><backlink>$$Uhttps://www.ncbi.nlm.nih.gov/pubmed/27628171$$D View this record in MEDLINE/PubMed$$Hfree_for_read</backlink></links><search><creatorcontrib>Harrison, Jonathan U.</creatorcontrib><creatorcontrib>Yates, Christian A.</creatorcontrib><title>A hybrid algorithm for coupling partial differential equation and compartment-based dynamics</title><title>Journal of the Royal Society interface</title><addtitle>J. R. Soc. Interface</addtitle><addtitle>J R Soc Interface</addtitle><description>Stochastic simulation methods can be applied successfully to model exact spatio-temporally resolved reaction–diffusion systems. However, in many cases, these methods can quickly become extremely computationally intensive with increasing particle numbers. An alternative description of many of these systems can be derived in the diffusive limit as a deterministic, continuum system of partial differential equations (PDEs). Although the numerical solution of such PDEs is, in general, much more efficient than the full stochastic simulation, the deterministic continuum description is generally not valid when copy numbers are low and stochastic effects dominate. Therefore, to take advantage of the benefits of both of these types of models, each of which may be appropriate in different parts of a spatial domain, we have developed an algorithm that can be used to couple these two types of model together. This hybrid coupling algorithm uses an overlap region between the two modelling regimes. By coupling fluxes at one end of the interface and using a concentration-matching condition at the other end, we ensure that mass is appropriately transferred between PDE- and compartment-based regimes. Our methodology gives notable reductions in simulation time in comparison with using a fully stochastic model, while maintaining the important stochastic features of the system and providing detail in appropriate areas of the domain. We test our hybrid methodology robustly by applying it to several biologically motivated problems including diffusion and morphogen gradient formation. Our analysis shows that the resulting error is small, unbiased and does not grow over time.</description><subject>Algorithms</subject><subject>Deterministic</subject><subject>Hybrid Algorithms</subject><subject>Life Sciences–Mathematics interface</subject><subject>Models, Theoretical</subject><subject>Multiscale</subject><subject>Reaction–diffusion</subject><subject>Stochastic</subject><issn>1742-5689</issn><issn>1742-5662</issn><fulltext>true</fulltext><rsrctype>article</rsrctype><creationdate>2016</creationdate><recordtype>article</recordtype><sourceid>EIF</sourceid><recordid>eNp1kc1LwzAYh4MofkyvHqVHL535XnoRhvgFAy96E0KaJltG28ykFfrfm7o59OApeXmf_PImDwCXCE4RLMRNiM5OMUR8CglhB-AUzSjOGef4cL8XxQk4i3ENIZkRxo7BCZ5xLNAMnYL3ebYayuCqTNVLH1y3ajLrQ6Z9v6ldu8w2KnRO1VnlrDXBtN-F-ehV53ybqbZKaDNCTerlpYqmyqqhVY3T8RwcWVVHc7FbJ-Dt4f717ilfvDw-380XuSacdjkrSlJxrAtVWiO0tYRaYgTHqoJMIJQeA8sS2VILzgXURNCC4oJDQ7HGmJIJuN3mbvqyMZVOkwRVy01wjQqD9MrJv53WreTSf0oGKS8oSQHXu4DgP3oTO9m4qE1dq9b4PkokMGQFo4QldLpFdfAxBmP31yAoRyVyVCJHJXJUkg5c_R5uj_84SADZAsEP6Ze8dqYb5Nr3oU3lf7FfyTCbpA</recordid><startdate>20160901</startdate><enddate>20160901</enddate><creator>Harrison, Jonathan U.</creator><creator>Yates, Christian A.</creator><general>The Royal Society</general><scope>CGR</scope><scope>CUY</scope><scope>CVF</scope><scope>ECM</scope><scope>EIF</scope><scope>NPM</scope><scope>AAYXX</scope><scope>CITATION</scope><scope>7X8</scope><scope>5PM</scope><orcidid>https://orcid.org/0000-0002-2748-9921</orcidid></search><sort><creationdate>20160901</creationdate><title>A hybrid algorithm for coupling partial differential equation and compartment-based dynamics</title><author>Harrison, Jonathan U. ; Yates, Christian A.</author></sort><facets><frbrtype>5</frbrtype><frbrgroupid>cdi_FETCH-LOGICAL-c364t-59b3d62c9abfe8cff34f3e862ad058116620bb1fbc86680c384942960e42c2243</frbrgroupid><rsrctype>articles</rsrctype><prefilter>articles</prefilter><language>eng</language><creationdate>2016</creationdate><topic>Algorithms</topic><topic>Deterministic</topic><topic>Hybrid Algorithms</topic><topic>Life Sciences–Mathematics interface</topic><topic>Models, Theoretical</topic><topic>Multiscale</topic><topic>Reaction–diffusion</topic><topic>Stochastic</topic><toplevel>peer_reviewed</toplevel><toplevel>online_resources</toplevel><creatorcontrib>Harrison, Jonathan U.</creatorcontrib><creatorcontrib>Yates, Christian A.</creatorcontrib><collection>Medline</collection><collection>MEDLINE</collection><collection>MEDLINE (Ovid)</collection><collection>MEDLINE</collection><collection>MEDLINE</collection><collection>PubMed</collection><collection>CrossRef</collection><collection>MEDLINE - Academic</collection><collection>PubMed Central (Full Participant titles)</collection><jtitle>Journal of the Royal Society interface</jtitle></facets><delivery><delcategory>Remote Search Resource</delcategory><fulltext>fulltext</fulltext></delivery><addata><au>Harrison, Jonathan U.</au><au>Yates, Christian A.</au><format>journal</format><genre>article</genre><ristype>JOUR</ristype><atitle>A hybrid algorithm for coupling partial differential equation and compartment-based dynamics</atitle><jtitle>Journal of the Royal Society interface</jtitle><stitle>J. R. Soc. Interface</stitle><addtitle>J R Soc Interface</addtitle><date>2016-09-01</date><risdate>2016</risdate><volume>13</volume><issue>122</issue><spage>20160335</spage><pages>20160335-</pages><issn>1742-5689</issn><eissn>1742-5662</eissn><abstract>Stochastic simulation methods can be applied successfully to model exact spatio-temporally resolved reaction–diffusion systems. However, in many cases, these methods can quickly become extremely computationally intensive with increasing particle numbers. An alternative description of many of these systems can be derived in the diffusive limit as a deterministic, continuum system of partial differential equations (PDEs). Although the numerical solution of such PDEs is, in general, much more efficient than the full stochastic simulation, the deterministic continuum description is generally not valid when copy numbers are low and stochastic effects dominate. Therefore, to take advantage of the benefits of both of these types of models, each of which may be appropriate in different parts of a spatial domain, we have developed an algorithm that can be used to couple these two types of model together. This hybrid coupling algorithm uses an overlap region between the two modelling regimes. By coupling fluxes at one end of the interface and using a concentration-matching condition at the other end, we ensure that mass is appropriately transferred between PDE- and compartment-based regimes. Our methodology gives notable reductions in simulation time in comparison with using a fully stochastic model, while maintaining the important stochastic features of the system and providing detail in appropriate areas of the domain. We test our hybrid methodology robustly by applying it to several biologically motivated problems including diffusion and morphogen gradient formation. Our analysis shows that the resulting error is small, unbiased and does not grow over time.</abstract><cop>England</cop><pub>The Royal Society</pub><pmid>27628171</pmid><doi>10.1098/rsif.2016.0335</doi><orcidid>https://orcid.org/0000-0002-2748-9921</orcidid><oa>free_for_read</oa></addata></record>
fulltext fulltext
identifier ISSN: 1742-5689
ispartof Journal of the Royal Society interface, 2016-09, Vol.13 (122), p.20160335
issn 1742-5689
1742-5662
language eng
recordid cdi_pubmed_primary_27628171
source MEDLINE; PubMed Central
subjects Algorithms
Deterministic
Hybrid Algorithms
Life Sciences–Mathematics interface
Models, Theoretical
Multiscale
Reaction–diffusion
Stochastic
title A hybrid algorithm for coupling partial differential equation and compartment-based dynamics
url https://sfx.bib-bvb.de/sfx_tum?ctx_ver=Z39.88-2004&ctx_enc=info:ofi/enc:UTF-8&ctx_tim=2024-12-17T17%3A14%3A42IST&url_ver=Z39.88-2004&url_ctx_fmt=infofi/fmt:kev:mtx:ctx&rfr_id=info:sid/primo.exlibrisgroup.com:primo3-Article-proquest_pubme&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.genre=article&rft.atitle=A%20hybrid%20algorithm%20for%20coupling%20partial%20differential%20equation%20and%20compartment-based%20dynamics&rft.jtitle=Journal%20of%20the%20Royal%20Society%20interface&rft.au=Harrison,%20Jonathan%20U.&rft.date=2016-09-01&rft.volume=13&rft.issue=122&rft.spage=20160335&rft.pages=20160335-&rft.issn=1742-5689&rft.eissn=1742-5662&rft_id=info:doi/10.1098/rsif.2016.0335&rft_dat=%3Cproquest_pubme%3E1820595435%3C/proquest_pubme%3E%3Curl%3E%3C/url%3E&disable_directlink=true&sfx.directlink=off&sfx.report_link=0&rft_id=info:oai/&rft_pqid=1820595435&rft_id=info:pmid/27628171&rfr_iscdi=true