| Peer-Reviewed

An Accurate and Stable Filtered Explicit Scheme for Biopolymerization Processes in the Presence of Perturbations

Received: 31 August 2021     Accepted: 13 October 2021     Published: 5 November 2021
Views:       Downloads:
Abstract

The focus of this paper is the development, numerical simulation and parameter analysis of a model of the transcription of ribosomal RNA in highly transcribed genes. Inspired by the well-known classic Lighthill-Whitham-Richards (LWR) traffic flow model, a linear advection continuum model is used to describe the DNA transcription process. In this model, elongation velocity is assumed to be essentially constant as RNA polymerases move along the strand through different phases of gene transcription. One advantage of using the linear model is that it allows one to quantify how small perturbations in elongation velocity and inflow parameters affect important biology measures such as Average Transcription Time (ATT) for the gene. The ATT per polymerase is the amount of time an individual RNAP spends traveling through the DNA strand. The numerical treatment for model simulations includes introducing a low complexity and time accurate method by adding a simple linear time filter to the classic upwind scheme. This improved method is modular and requires a minimal modification of adding only one line of code resulting in increased accuracy without increased computational expense. In addition, it removes the overdamping of upwind. A stability condition for the new algorithm is derived, and numerical computations illustrate stability and convergence of the filtered scheme as well as improved ATT estimation.

Published in Applied and Computational Mathematics (Volume 10, Issue 6)
DOI 10.11648/j.acm.20211006.11
Page(s) 121-137
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2021. Published by Science Publishing Group

Keywords

Advection Equation, Lighthill-Whitham-Richards Model, Ribosomal RNA, RNA Polymerases, Time Filter, Traffic Flow Models, Transcription Time, Upwind Scheme

References
[1] T. W. Turowski and E. Petfalski, B. D. Goddard, S. L. French, A. Helwak and D. Tollervey (2020) “Nascent Transcript Folding Plays a Major Role in Determining RNA Polymerase Elongation Rates,” Molecular Cell, 79: 3, pp. 488-503.
[2] A. Klindziuk and A. Kolomeisky (2021) “Long-Range Supercoiling-Mediated RNA Polymerase Cooperation in Transcription,” J. Phys. Chem. B 2021, 125: 18, 4692- 4700.
[3] J. Saba, X. Y. Chua, T. V. Mishanina, D. Nayak, T. A. Windgassen, R. A. Mooney and R. Landick (2019) The elemental mechanism of transcriptional pausing. Elife. 2019 Jan 8; 8: e40981. doi: 10.7554/eLife.40981.
[4] J. Y. Kang, T. V. Mishanina, R. Landick, S. A. Darst (2019) Mechanisms of Transcriptional Pausing in Bacteria. J Mol Biol. 2019; 431 (20): 4007-4029. doi: 10.1016/j.jmb.2019.07.017.
[5] C. T. MacDonald, J. H. Gibbs, J. H. and A. C. Pipkin, A. C., (1968) “Kinetics of biopolymerization on nucleic acid templates”, Biopolymers, Vol. 6, p. 1-25.
[6] C. T. MacDonald and J. H. Gibbs (1969) “Concerning the kinetics of polypeptide synthesis on polyribosomes”, Biopolymers, Vol. 7, p. 707-725.
[7] F. Spitzer, Frank (1970) “Interaction of Markov processes”, Advances in Math., Vol. 5, pp. 246 - 290.
[8] R. Heinrich and T. A. Rapoport (1980) “Mathematical modelling of translation of mRNA in Eucaryotes: Steady states, time-dependent processes and application to reticulocytes”, J. Theor. Biol., Vol. 86, pp. 279-313.
[9] B. Derrida, M. Evans, V. Hakim and V. Pasquier (1993) “Exact solution of a 1D asymmetric exclusion model using a matrix formalism,” J. Phys. A: Math. Gen., Vol. 26, pp. 1493-1517.
[10] L. B. Shaw, R. K. Z. Zia and K. H. Lee, K. H. (2003) “Totally asymmetric exclusion process with extended objects: A model for protein synthesis,” Phys. Rev. E, Vol. 68, 021910.
[11] A. Kolomeisky (1998) “Asymetric simple exclusion model with local inhomogeneity”, J. Phys. A, Vol. 31, pp. 1153-1164.
[12] S. Klumpp and T. Hwa (2008) “Stochasticity and traffic jams in the transcription of ribosomal RNA: Intriguing role of termination and antitermination,” PNAS, Vol. 105, No. 47, pp. 18159-164.
[13] S. Klumpp (2011) “Pausing and backtracking in transcription under dense traffic conditions,” J. Stat. Phys., Vol. 142, pp. 1252 - 1267.
[14] C. A. Brackley, M. C. Romano, M. C., C. Grebogi and M. Thiel, M. (2010) “Limited resources in a driven diffusion process,”, Phys. Rev. Let., Vol. 105.
[15] L. Ciandrini, L., I. Stansfield and M. C. Romano, M. C. (2010) “Role of the particle’s stepping cycle in an asymmetric exclusion process: model of mRNA translation,” Phys. Rev. E, Vol. 81, No. 2.
[16] J. Wang, J., B. Pfeuty, Q. Thommen, M. C. Romano, M. C. and M. Lefranc, M. (2014) “Minimal model of transcriptional elongation processes with pauses,” Phys. Rev. E, Vol. 90.
[17] L. Mier-y-Teran-R, M. Silber and V. Hatzimanakatis (2010) “The origins of time-delay in template bio- polymerization processes,” PloS Comput Biology, Vol. 6, No. 4, pp. e1000726.
[18] R. K. P. Zia, J. J. Dong, and B Schmittmann, B. (2011) “Modeling translation in protein synthesis with TASEP: a tutorial and recent developments,” J. Stat. Phys., Vol. 144, No. 2, pp. 405 - 428.
[19] L. Davis, T. Gedeon, J. Gedeon, and J. Thorenson (2013) “A Traffic Flow Model for Bio-polymerization Processes,” J. Math. Bio., pp. 1416-1432.
[20] L. Davis, T. Gedeon and J. Thorenson (2014) “Discontinuous Galerkin Calculations for a Nonlinear PDE Model of DNA Transcription with Short, Transient and Frequent Pausing,” J. Comp. Math. Vol. 32, No. 6, pp. 601-629.
[21] L. Davis, T. Gedeon and J. Thorenson (2014) “Parameter Studies for a Nonlinear Continuum Model of Transcription,” Proceedings of 21st International Symposium on the Mathematical Theory of Networks and Systems.
[22] C. F. Daganzo (1995) “Requiem for second-order fluid approximations of traffic flow,” Transportation Research Part B: Methodological, Vol. 29, No. 4, pp. 277-286.
[23] A. Aw and M. Rascle (2000) “Resurrection of “second order” models of traffic flow,” SIAM J. Appl. Math., Vol. 60, No. 3, pp. 916–938.
[24] H. M. Zhang (2002) “A non-equilibrium traffic model devoid of gas-like behavior,” Transportation Research Part B: Methodological, 36: 3 275-290.
[25] R. Mohan and G. Ramadurai (2015) “Submission to the DTA2012 Special Issue: A Case for Higher-Order Traffic Flow Models in DTA,” Networks and spatial economics, 15: 3 765-790.
[26] R. Haberman (1998) Mathematical models: mechanical vibrations, population dynamics, and traffic flow: an introduction to applied mathematics, in series Classics in Applied Mathematics, SIAM, Philadelphia.
[27] M. D. Rosini (2013) Macroscopic Models for Vehicular Flows and Crowd Dynamics: Theory and Applications Classical and Non-Classical Advanced Mathematics for Real Life Applications in series Understanding Complex Systems, Springer.
[28] R. J. LeVeque (2002) “Finite volume methods for hyperbolic problems,” in series Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge.
[29] K. W. Morton and D. F. Mayers (2005) “Numerical Solution of Partial Differential Equations: An Introduction,” Ed. 2, Cambridge University Press, Cambridge.
[30] J. W. Thomas (2013) Numerical Partial Differential Equations: Finite Difference Methods in series Texts in Applied Mathematics, Vol. 22, Springer.
[31] J. W. Thomas (2013) Numerical Partial Differential Equations: Conservation Laws and Elliptic Equations, in series Texts in Applied Mathematics, Vol. 33, Springer, New York.
[32] A. Guzel and W. Layton (2018) “Time filters increase accuracy of the fully implicit method,” BIT Numerical Mathematics, 58: 2 301-315.
[33] R. E. Schlesinger, L. W. Uccellini, Louis W and D. R. Johnson (1983) “The Effects of the Asselin Time Filter on Numerical Solutions to the Linearized Shallow-Water Wave Equations,” Monthly Weather Review, 111: 3 455- 467.
[34] P. P. Dennis, M. Ehrenberg, D. Fange and H. Bremer (2009) “Varying Rate of RNA Chain Elongation during rrn Transcription in Escherichia coli,” Journal of Bacteriology, 191: 11 3740-3746.34
Cite This Article
  • APA Style

    Lisa Davis, Faranak Pahlevani, Timmy Susai Rajan. (2021). An Accurate and Stable Filtered Explicit Scheme for Biopolymerization Processes in the Presence of Perturbations. Applied and Computational Mathematics, 10(6), 121-137. https://doi.org/10.11648/j.acm.20211006.11

    Copy | Download

    ACS Style

    Lisa Davis; Faranak Pahlevani; Timmy Susai Rajan. An Accurate and Stable Filtered Explicit Scheme for Biopolymerization Processes in the Presence of Perturbations. Appl. Comput. Math. 2021, 10(6), 121-137. doi: 10.11648/j.acm.20211006.11

    Copy | Download

    AMA Style

    Lisa Davis, Faranak Pahlevani, Timmy Susai Rajan. An Accurate and Stable Filtered Explicit Scheme for Biopolymerization Processes in the Presence of Perturbations. Appl Comput Math. 2021;10(6):121-137. doi: 10.11648/j.acm.20211006.11

    Copy | Download

  • @article{10.11648/j.acm.20211006.11,
      author = {Lisa Davis and Faranak Pahlevani and Timmy Susai Rajan},
      title = {An Accurate and Stable Filtered Explicit Scheme for Biopolymerization Processes in the Presence of Perturbations},
      journal = {Applied and Computational Mathematics},
      volume = {10},
      number = {6},
      pages = {121-137},
      doi = {10.11648/j.acm.20211006.11},
      url = {https://doi.org/10.11648/j.acm.20211006.11},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.acm.20211006.11},
      abstract = {The focus of this paper is the development, numerical simulation and parameter analysis of a model of the transcription of ribosomal RNA in highly transcribed genes. Inspired by the well-known classic Lighthill-Whitham-Richards (LWR) traffic flow model, a linear advection continuum model is used to describe the DNA transcription process. In this model, elongation velocity is assumed to be essentially constant as RNA polymerases move along the strand through different phases of gene transcription. One advantage of using the linear model is that it allows one to quantify how small perturbations in elongation velocity and inflow parameters affect important biology measures such as Average Transcription Time (ATT) for the gene. The ATT per polymerase is the amount of time an individual RNAP spends traveling through the DNA strand. The numerical treatment for model simulations includes introducing a low complexity and time accurate method by adding a simple linear time filter to the classic upwind scheme. This improved method is modular and requires a minimal modification of adding only one line of code resulting in increased accuracy without increased computational expense. In addition, it removes the overdamping of upwind. A stability condition for the new algorithm is derived, and numerical computations illustrate stability and convergence of the filtered scheme as well as improved ATT estimation.},
     year = {2021}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - An Accurate and Stable Filtered Explicit Scheme for Biopolymerization Processes in the Presence of Perturbations
    AU  - Lisa Davis
    AU  - Faranak Pahlevani
    AU  - Timmy Susai Rajan
    Y1  - 2021/11/05
    PY  - 2021
    N1  - https://doi.org/10.11648/j.acm.20211006.11
    DO  - 10.11648/j.acm.20211006.11
    T2  - Applied and Computational Mathematics
    JF  - Applied and Computational Mathematics
    JO  - Applied and Computational Mathematics
    SP  - 121
    EP  - 137
    PB  - Science Publishing Group
    SN  - 2328-5613
    UR  - https://doi.org/10.11648/j.acm.20211006.11
    AB  - The focus of this paper is the development, numerical simulation and parameter analysis of a model of the transcription of ribosomal RNA in highly transcribed genes. Inspired by the well-known classic Lighthill-Whitham-Richards (LWR) traffic flow model, a linear advection continuum model is used to describe the DNA transcription process. In this model, elongation velocity is assumed to be essentially constant as RNA polymerases move along the strand through different phases of gene transcription. One advantage of using the linear model is that it allows one to quantify how small perturbations in elongation velocity and inflow parameters affect important biology measures such as Average Transcription Time (ATT) for the gene. The ATT per polymerase is the amount of time an individual RNAP spends traveling through the DNA strand. The numerical treatment for model simulations includes introducing a low complexity and time accurate method by adding a simple linear time filter to the classic upwind scheme. This improved method is modular and requires a minimal modification of adding only one line of code resulting in increased accuracy without increased computational expense. In addition, it removes the overdamping of upwind. A stability condition for the new algorithm is derived, and numerical computations illustrate stability and convergence of the filtered scheme as well as improved ATT estimation.
    VL  - 10
    IS  - 6
    ER  - 

    Copy | Download

Author Information
  • Department of Mathematical Sciences, Montana State University, Bozeman, United States

  • Division of Science & Engineering, Penn State University-Abington, Abington, United States

  • Division of Science & Engineering, Penn State University-Abington, Abington, United States

  • Sections