-
Distinguishing Feedback Mechanisms in Clock Models
Authors:
Alexander D. Golden,
Joris Paijmans,
David K. Lubensky
Abstract:
Biological oscillators are very diverse but can be classified based on dynamical motifs such as the types of feedback loops present. The S. Elongatus circadian clock is a remarkable phosphorylation-based oscillator that can be reconstituted in vitro with only 3 different purified proteins: the clock proteins KaiA, KaiB, and KaiC. Despite a growing body of knowledge about the biochemistry of the Ka…
▽ More
Biological oscillators are very diverse but can be classified based on dynamical motifs such as the types of feedback loops present. The S. Elongatus circadian clock is a remarkable phosphorylation-based oscillator that can be reconstituted in vitro with only 3 different purified proteins: the clock proteins KaiA, KaiB, and KaiC. Despite a growing body of knowledge about the biochemistry of the Kai proteins, basic questions about how their interactions lead to sustained oscillations remain unanswered. Here, we compare models of this system that make opposing assumptions about whether KaiA sequestration introduces a positive or a negative feedback loop. We find that the two different feedback mechanisms can be distinguished experimentally by the introduction of a protein that binds competitively with KaiA. Understanding the dynamical mechanism responsible for oscillations in the Kai system may shed light on the broader question of what clock architectures have been selected by evolution and why.
△ Less
Submitted 4 June, 2018; v1 submitted 28 May, 2018;
originally announced May 2018.
-
Sequencing single-stranded libraries on the Illumina NextSeq 500 platform
Authors:
Johanna L. A. Paijmans,
Sina Baleka,
Kirstin Henneberger,
Ulrike H. Taron,
Alexandra Trinks,
Michael V. Westbury,
Axel Barlow
Abstract:
Single-stranded libraries generated from ancient DNA extracts have specific sequencing requirements. The short library molecules typically associated with ancient DNA templates are expected to behave differently during the annealing to the flow cell and subsequent cluster generation (bridge PCR), requiring optimisation of loading amounts to obtain optimal and consistent cluster densities. For the…
▽ More
Single-stranded libraries generated from ancient DNA extracts have specific sequencing requirements. The short library molecules typically associated with ancient DNA templates are expected to behave differently during the annealing to the flow cell and subsequent cluster generation (bridge PCR), requiring optimisation of loading amounts to obtain optimal and consistent cluster densities. For the past 3 years, we have carried out sequencing of single-stranded libraries on the Illumina NextSeq 500 sequencing platform at the Institute for Biochemistry and Biology, University of Potsdam. We report our optimisations here. This document may be useful for other researchers wishing to sequence single-stranded libraries on the NextSeq 500 platform. It does not replace the excellent documentation provided by Illumina (links provided below), but rather serves as additional information specific to single-stranded libraries. The original papers describing the library protocol should also be studied in detail, and complement the information presented here.
△ Less
Submitted 29 November, 2017;
originally announced November 2017.
-
eGFRD in all dimensions
Authors:
Thomas R. Sokolowski,
Joris Paijmans,
Laurens Bossen,
Martijn Wehrens,
Thomas Miedema,
Nils B. Becker,
Kazunari Kaizu,
Koichi Takahashi,
Marlieen Dogterom,
Pieter Rein ten Wolde
Abstract:
Biochemical reactions typically occur at low copy numbers, but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. However, while particle-based models…
▽ More
Biochemical reactions typically occur at low copy numbers, but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. However, while particle-based models guarantee the level of detail necessary to accurately describe the microscopic dynamics at very low copy numbers, the algorithms used to simulate them oftentimes imply trade-offs between computational efficiency and accuracy. eGFRD (enhanced Green's Function Reaction Dynamics) is an exact algorithm that evades such trade-offs by partitioning the N-particle system into M<N analytically tractable one- and two-particle systems; the analytical solutions (Green's functions) then are used to implement an event-driven particle-based scheme that allows particles to make large jumps in time and space while retaining access to their state variables at any moment. Here we present "eGFRD2", a new eGFRD version that implements the principle of eGFRD in all dimensions, enabling efficient simulation of biochemical reaction-diffusion processes in the 3D cytoplasm, on 2D planes representing membranes, and on 1D elongated cylinders representative of, e.g., cytoskeletal tracks or DNA; in 1D, it also incorporates convective motion used to model active transport. We find that, for low particle densities, eGFRD2 is up to 3 orders of magnitude faster than optimized Brownian Dynamics. We exemplify the capabilities of eGFRD2 by simulating an idealized model of Pom1 gradient formation, which involves 3D diffusion, active transport on microtubules, and autophosphorylation on the membrane, confirming recent results on this system and demonstrating that it can efficiently operate under genuinely stochastic conditions.
△ Less
Submitted 30 August, 2017;
originally announced August 2017.
-
Period Robustness and Entrainability of the Kai System to Changing Nucleotide Concentrations
Authors:
Joris Paijmans,
David K Lubensky,
Pieter Rein ten Wolde
Abstract:
Circadian clocks must be able to entrain to time-varying signals to keep their oscillations in phase with the day-night rhythm. On the other hand, they must also exhibit input compensation: their period must remain about one day in different constant environments. The post-translational oscillator of the Kai system can be entrained by transient or oscillatory changes in the ATP fraction, yet is in…
▽ More
Circadian clocks must be able to entrain to time-varying signals to keep their oscillations in phase with the day-night rhythm. On the other hand, they must also exhibit input compensation: their period must remain about one day in different constant environments. The post-translational oscillator of the Kai system can be entrained by transient or oscillatory changes in the ATP fraction, yet is insensitive to constant changes in this fraction. We study in three different models of this system how these two seemingly conflicting criteria are met: the Van Zon model (Van Zon et al., PNAS, 2007), the Rust model (Phong et al., PNAS, 2013), and our new model (Paijmans et al., arXiv:1612.02715). We find that the new model exhibits the best trade-off between input compensation and entrainability: on the footing of equal phase-response curves, it exhibits the strongest input compensation. Performing stochastic simulations at the level of individual hexamers allows us to identify a new mechanism, which is employed by the new model to achieve input compensation: At lower ATP fraction, the individual hexamers make a shorter cycle in the phosphorylation state space, which compensates for the slower pace at which they traverse the cycle.
△ Less
Submitted 22 June, 2017; v1 submitted 25 December, 2016;
originally announced December 2016.
-
Robustness of synthetic oscillators in growing and dividing cells
Authors:
Joris Paijmans,
David K Lubensky,
Pieter Rein ten Wolde
Abstract:
Synthetic biology sets out to implement new functions in cells, and to develop a deeper understanding of biological design principles. In 2000, Elowitz and Leibler showed that by rational design of the reaction network, and using existing biological components, they could create a network that exhibits periodic gene expression, dubbed the repressilator (Elowitz and Leibler, Nature, 2000). More rec…
▽ More
Synthetic biology sets out to implement new functions in cells, and to develop a deeper understanding of biological design principles. In 2000, Elowitz and Leibler showed that by rational design of the reaction network, and using existing biological components, they could create a network that exhibits periodic gene expression, dubbed the repressilator (Elowitz and Leibler, Nature, 2000). More recently, Stricker et al. presented another synthetic oscillator, called the dual-feedback oscillator (Stricker et al., 2008), which is more stable. How the stability of these oscillators is affected by the intrinsic noise of the interactions between the components and the stochastic expression of their genes, has been studied in considerable detail. However, as all biological oscillators reside in growing and dividing cells, an important question is how these oscillators are perturbed by the cell cycle. In previous work we showed that the periodic doubling of the gene copy numbers due to DNA replication can couple not only natural, circadian oscillators to the cell cycle (Paijmans et al., PNAS, \textbf{113}, 4063, (2016)), but also these synthetic oscillators. Here we expand this study. We find that the strength of the locking between oscillators depends not only on the positions of the genes on the chromosome, but also on the noise in the timing of gene replication: noise tends to weaken the coupling. Yet, even in the limit of high levels of noise in the replication times of the genes, both synthetic oscillators show clear signatures of locking to the cell cycle. This work enhances our understanding of the design of robust biological oscillators inside growing and diving cells.
△ Less
Submitted 23 December, 2016;
originally announced December 2016.
-
A Thermodynamically consistent model of the post-translational Kai circadian clock
Authors:
Joris Paijmans,
David K Lubensky,
Pieter Rein ten Wolde
Abstract:
The principal pacemaker of the circadian clock of the cyanobacterium S. elongatus is a protein phosphorylation cycle consisting of three proteins, KaiA, KaiB and KaiC. KaiC forms a homohexamer, with each monomer consisting of two domains, CI and CII. Both domains can bind and hydrolyze ATP, but only the CII domain can be phosphorylated, at two residues, in a well-defined sequence. While this syste…
▽ More
The principal pacemaker of the circadian clock of the cyanobacterium S. elongatus is a protein phosphorylation cycle consisting of three proteins, KaiA, KaiB and KaiC. KaiC forms a homohexamer, with each monomer consisting of two domains, CI and CII. Both domains can bind and hydrolyze ATP, but only the CII domain can be phosphorylated, at two residues, in a well-defined sequence. While this system has been studied extensively, how the clock is driven thermodynamically has remained elusive. Inspired by recent experimental observations and building on ideas from previous mathematical models, we present a new, thermodynamically consistent, statistical-mechanical model of the clock. At its heart are two main ideas: i) ATP hydrolysis in the CI domain provides the thermodynamic driving force for the clock, switching KaiC between an active conformational state in which its phosphorylation level tends to rise and an inactive one in which it tends to fall; ii) phosphorylation of the CII domain provides the timer for the hydrolysis in the CI domain. The model also naturally explains how KaiA, by acting as a nucleotide exchange factor, can stimulate phosphorylation of KaiC, and how the differential affinity of KaiA for the different KaiC phosphoforms generates the characteristic temporal order of KaiC phosphorylation. As the phosphorylation level in the CII domain rises, the release of ADP from CI slows down, making the inactive conformational state of KaiC more stable. In the inactive state, KaiC binds KaiB, which not only stabilizes this state further, but also leads to the sequestration of KaiA, and hence to KaiC dephosphorylation. Using a dedicated kinetic Monte Carlo algorithm, we show that the model can describe a wealth of experimental data.
△ Less
Submitted 8 December, 2016;
originally announced December 2016.
-
Discrete gene replication events drive coupling between the cell cycle and circadian clocks
Authors:
Joris Paijmans,
Mark Bosman,
Pieter Rein ten Wolde,
David K. Lubensky
Abstract:
Many organisms possess both a cell cycle to control DNA replication and a circadian clock to anticipate changes between day and night. In some cases, these two rhythmic systems are known to be coupled by specific, cross-regulatory interactions. Here, we use mathematical modeling to show that, additionally, the cell cycle generically influences circadian clocks in a non-specific fashion: The regula…
▽ More
Many organisms possess both a cell cycle to control DNA replication and a circadian clock to anticipate changes between day and night. In some cases, these two rhythmic systems are known to be coupled by specific, cross-regulatory interactions. Here, we use mathematical modeling to show that, additionally, the cell cycle generically influences circadian clocks in a non-specific fashion: The regular, discrete jumps in gene-copy number arising from DNA replication during the cell cycle cause a periodic driving of the circadian clock, which can dramatically alter its behavior and impair its function. A clock built on negative transcriptional feedback either phase locks to the cell cycle, so that the clock period tracks the cell division time, or exhibits erratic behavior. We argue that the cyanobacterium Synechococcus elongatus has evolved two features that protect its clock from such disturbances, both of which are needed to fully insulate it from the cell cycle and give it its observed robustness: a phosphorylation-based protein modification oscillator, together with its accompanying push-pull read-out circuit that responds primarily to the ratios of the different phosphoforms, makes the clock less susceptible to perturbations in protein synthesis; and the presence of multiple, asynchronously replicating copies of the same chromosome diminishes the effect of replicating any single copy of a gene.
△ Less
Submitted 20 October, 2015;
originally announced October 2015.
-
Analysis of whole mitogenomes from ancient samples
Authors:
Gloria G. Fortes,
Johanna L. A. Paijmans
Abstract:
Ancient mitochondrial DNA has been used in a wide variety of palaeontological and archaeological studies, ranging from population dynamics of extinct species to patterns of domestication. Most of these studies have traditionally been based on the analysis of short fragments from the mitochondrial control region, analysed using PCR coupled with Sanger sequencing. With the introduction of high-throu…
▽ More
Ancient mitochondrial DNA has been used in a wide variety of palaeontological and archaeological studies, ranging from population dynamics of extinct species to patterns of domestication. Most of these studies have traditionally been based on the analysis of short fragments from the mitochondrial control region, analysed using PCR coupled with Sanger sequencing. With the introduction of high-throughput sequencing, as well as new enrichment technologies, the recovery of full mitochondrial genomes (mitogenomes) from ancient specimens has become significantly less complicated. Here we present a protocol to build ancient extracts into Illumina high-throughput sequencing libraries, and subsequent Agilent array-based capture to enrich for the desired mitogenome. Both are based on previously published protocols, with the introduction of several improvements aimed to increase the recovery of short DNA fragments, while keeping the cost and effort requirements low. This protocol was designed for enrichment of mitochondrial DNA in ancient or degraded samples. However, the protocols can be easily adapted for using for building libraries for shotgun-sequencing of whole genomes, or enrichment of other genomic regions.
△ Less
Submitted 25 April, 2016; v1 submitted 17 March, 2015;
originally announced March 2015.
-
The lower bound on the precision of transcriptional regulation
Authors:
Joris Paijmans,
Pieter Rein ten Wolde
Abstract:
The diffusive arrival of transcription factors at the promoter sites on the DNA sets a lower bound on how accurately a cell can regulate its protein levels. Using results from the literature on diffusion-influenced reactions, we derive an analytical expression for the lower bound on the precision of transcriptional regulation. In our theory, transcription factors can perform multiple rounds o…
▽ More
The diffusive arrival of transcription factors at the promoter sites on the DNA sets a lower bound on how accurately a cell can regulate its protein levels. Using results from the literature on diffusion-influenced reactions, we derive an analytical expression for the lower bound on the precision of transcriptional regulation. In our theory, transcription factors can perform multiple rounds of 1D diffusion along the DNA and 3D diffusion in the cytoplasm before binding to the promoter. Comparing our expression for the lower bound on the precision against results from Green's Function Reaction Dynamics simulations shows that the theory is highly accurate under biologically relevant conditions. Our results demonstrate that, to an excellent approximation, the promoter switches between the transcription-factor bound and unbound state in a Markovian fashion. This remains true even in the presence of sliding, i.e. with 1D diffusion along the DNA. This has two important implications: (1) minimizing the noise in the promoter state is equivalent to minimizing the search time of transcription factors for their promoters; (2) the complicated dynamics of 3D diffusion in the cytoplasm and 1D diffusion along the DNA can be captured in a well-stirred model by renormalizing the promoter association and dissociation rates, making it possible to efficiently simulate the promoter dynamics using Gillespie simulations. Based on the recent experimental observation that sliding can speed up the promoter search by a factor of 4, our theory predicts that sliding can enhance the precision of transcriptional regulation by a factor of 2.
△ Less
Submitted 3 July, 2014;
originally announced July 2014.
