skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An Accurate and Stable Filtered Explicit Scheme for Biopolymerization Processes in the Presence of Perturbations
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.  more » « less
Award ID(s):
1951510
PAR ID:
10339533
Author(s) / Creator(s):
Date Published:
Journal Name:
Applied and computational mathematics
Volume:
10
Issue:
6
ISSN:
2328-5605
Page Range / eLocation ID:
121-137
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Daròs, José-Antonio (Ed.)
    Viroids, a fascinating group of plant pathogens, are subviral agents composed of single-stranded circular noncoding RNAs. It is well-known that nuclear-replicating viroids exploit host DNA-dependent RNA polymerase II (Pol II) activity for transcription from circular RNA genome to minus-strand intermediates, a classic example illustrating the intrinsic RNA-dependent RNA polymerase activity of Pol II. The mechanism for Pol II to accept single-stranded RNAs as templates remains poorly understood. Here, we reconstituted a robust in vitro transcription system and demonstrated that Pol II also accepts minus-strand viroid RNA template to generate plus-strand RNAs. Further, we purified the Pol II complex on RNA templates for nano-liquid chromatography-tandem mass spectrometry analysis and identified a remodeled Pol II missing Rpb4, Rpb5, Rpb6, Rpb7, and Rpb9, contrasting to the canonical 12-subunit Pol II or the 10-subunit Pol II core on DNA templates. Interestingly, the absence of Rpb9, which is responsible for Pol II fidelity, explains the higher mutation rate of viroids in comparison to cellular transcripts. This remodeled Pol II is active for transcription with the aid of TFIIIA-7ZF and appears not to require other canonical general transcription factors (such as TFIIA, TFIIB, TFIID, TFIIE, TFIIF, TFIIH, and TFIIS), suggesting a distinct mechanism/machinery for viroid RNA-templated transcription. Transcription elongation factors, such as FACT complex, PAF1 complex, and SPT6, were also absent in the reconstituted transcription complex. Further analyses of the critical zinc finger domains in TFIIIA-7ZF revealed the first three zinc finger domains pivotal for RNA template binding. Collectively, our data illustrated a distinct organization of Pol II complex on viroid RNA templates, providing new insights into viroid replication, the evolution of transcription machinery, as well as the mechanism of RNA-templated transcription. 
    more » « less
  2. Kursula, Petri (Ed.)
    Short segments of RNA displace one strand of a DNA duplex during diverse processes including transcription and CRISPR-mediated immunity and genome editing. These strand exchange events involve the intersection of two geometrically distinct helix types—an RNA:DNA hybrid (A-form) and a DNA:DNA homoduplex (B-form). Although previous evidence suggests that these two helices can stack on each other, it is unknown what local geometric adjustments could enable A-on-B stacking. Here we report the X-ray crystal structure of an RNA-5′/DNA-3′ strand exchange junction at an anisotropic resolution of 1.6 to 2.2 Å. The structure reveals that the A-to-B helical transition involves a combination of helical axis misalignment, helical axis tilting and compression of the DNA strand within the RNA:DNA helix, where nucleotides exhibit a mixture of A- and B-form geometry. These structural principles explain previous observations of conformational stability in RNA/DNA exchange junctions, enabling a nucleic acid architecture that is repeatedly populated during biological strand exchange events. 
    more » « less
  3. Abstract Transcription by RNA polymerases is an exquisitely regulated step of the central dogma. Transcription is the primary determinant of cell-state, and most cellular perturbations impact transcription by altering polymerase activity. Thus, detecting changes in polymerase activity yields insight into most cellular processes. Nascent run-on sequencing provides a direct readout of polymerase activity, but no tools exist to model all aspects of this activity at genes. We focus on RNA polymerase II—responsible for transcribing protein-coding genes. We present the first model to capture the complete process of gene transcription. For individual genes, this model parameterizes each distinct stage of transcription—loading, initiation, elongation, and termination, hence LIET—in a biologically interpretable Bayesian mixture, which is applied to nascent run-on data. Our improved modeling of loading/initiation demonstrates these stages are characteristically different between sense and antisense strands. Applying LIET to 24 human cell-types, our analysis indicates the position of dissociation (the last step of termination) appears to be highly consistent, indicative of a tightly regulated process. Furthermore, by applying LIET to perturbation experiments, we demonstrate its ability to detect specific changes in pausing (5′ end), strand-bias, and dissociation location (3′ end)—opening the door to differential assessment of transcription at individual stages of individual genes. 
    more » « less
  4. null (Ed.)
    The sequence of the DNA template has long been thought to influence the rate of transcription by DNA-dependent RNA polymerases, but the influence of DNA sequence on transcription elongation properties of eukaryotic RNA polymerase I (Pol I) from Saccharomyces cerevisiae has not been defined. In this study, we observe changes in dinucleotide production, transcription elongation complex stability, and Pol I pausing in vitro in response to downstream DNA. In vitro studies demonstrate that AT-rich downstream DNA enhances pausing by Pol I and inhibits Pol I nucleolytic cleavage activity. Analysis of Pol I native elongating transcript sequencing data in Saccharomyces cerevisiae suggests that these downstream sequence elements influence Pol I in vivo . Native elongating transcript sequencing studies reveal that Pol I occupancy increases as downstream AT content increases and decreases as downstream GC content increases. Collectively, these data demonstrate that the downstream DNA sequence directly impacts the kinetics of transcription elongation prior to the sequence entering the active site of Pol I both in vivo and in vitro . 
    more » « less
  5. Abstract Transcription is the critical first step in expressing a gene, during which an RNA polymerase (RNAP) synthesizes an RNA copy of one strand of the DNA that encodes a gene. Here we describe a laboratory experiment that uses a single assay to probe two important steps in transcription: (1) RNAP binding to DNA, and (2) the transcriptional activity of the polymerase. Students probe both these steps in a single experiment using a fluorescence‐based electrophoretic mobility shift assay (EMSA) and commercially availableEscherichia coliRNAP. As an inquiry‐driven component, students add the transcriptional inhibitor rifampicin to reactions and draw conclusions about its mechanism of inhibition by determining whether it blocks polymerase binding to DNA or transcriptional activity. Depending on the curriculum and learning goals of individual courses, this experimental module could be easily expanded to include additional experimentation that mimics a research environment more closely. After completing the experiment students understand basic principles of transcription, mechanisms of inhibition, and the use of EMSAs to probe protein/DNA interactions. 
    more » « less