Stem Cell Reprogramming
|
|
- Charlotte Warren
- 6 years ago
- Views:
Transcription
1 Stem Cell Reprogramming Colin McDonnell 1 Introduction Since the demonstration of adult cell reprogramming by Yamanaka in 2006, there has been immense research interest in modelling and understanding the complicated un differentiation process. The addition of four transcription factors, or nucleic acid encoding the same, can fully un differentiate a mature adult somatic cell in a matter of days. These un differentiated cells can then be reprogrammed to a new state, providing a source of individualized pluripotent stem cells for use in autologous cell replacement therapies, which have the potential to treat or cure Alzheimer's disease, Parkinson's disease, cardiovascular disease, and diabetes, among others. Here, we study a degenerate case involving just two proteins X 1 and X 2, representing an arbitrary pair of the four Yamanaka factors. Each positively autoregulates itself and negatively regulates the other, as shown in Figure 1. Figure 1: Each protein activates itself and represses the other. Without positive autoregulation, this circuit reduces to a standard bi stable toggle switch. We develop a model that demonstrates both tri stability within certain parameter ranges and the existence of a weakly stable un differentiated state.
2 2 Tri stable ODE Model The behavior of the system is modelled using standard ODEs as follows. 2.1 Parameters We see that X 1 has some basal expression a 01, a positive autoregulation term with coefficient a 1, a X 2 dependent repression term with coefficient b 1, and a degradation term with constant δ 1. Similar terms exist for the ODE governing X 2 expression. The autoregulation and repression terms follow a Hill model with dissociation constants of K 1 and K 2 and cooperativity constants of n and m. These are all the parameters of the system. 2.2 Conditions for tri stability We now examine the qualitative parameter ranges over which this system demonstrates tristable behavior. For now, we assume zero basal expression ( a 01 = a 02 = 0 ) and a fully symmetric system, that is a 1 = a 2, b 1 = b 2, etc. We therefore eliminate the subscripts for simplicity. a The positive autoregulation constant a is related to the x 3 behavior of the system required for 5 equilibria to emerge. As a approaches zero, the system degenerates to the standard bi stable genetic toggle switch, which is incapable of tristable behavior. Thus a should be the same order of magnitude as b to avoid bi stability. b K n = m Increasing b actually increases the stability of the un differentiated state. As it gets much larger than a however, there is a tipping point where instantaneously goes unstable and bi stable behavior emerges. The dissociation constants have the opposite effect of a. Increasing K acts as a damper on the x 3 behavior, thus it should be small relative to a. Though any increase in K can be counteracted by an appropriate increase in a, the absolute values of the equilibria in terms of X 1 and concentrations shift outwards. Decreasing the cooperativity exponent damps x 3 behavior and, if sufficiently large, results in a non linear change from tri stability to mono stability. Increasing cooperativity also results in X 2
3 nonlinear behavior, even resulting in tetra stability for sufficiently large values. There is a happy medium around 2 where tristable behavior can exist. δ The degradation coefficient is inherently a damping term. At large values the system is mono stable. As the value decreases, the system becomes briefly bi stable before the onset of tri stability. After this onset, increasing δ increases the stability of. 2.3 Parameter selections After some trial and error, we settled on a set of parameter values that resulted in tri stable behavior while maximizing ease of implementation and analysis. The results of nullcline stability analysis are shown in Figure 2. a = 2 b =. 5 K = 1 n = m = 2 δ = 1 Figure 2: Nullclines for given parameter values. Note the 5 equilibria, three of which are stable. 3 Noise mediated Differentiation un differentiated cells frequently differentiate due to stochastic noise in the system during the creation of induced PSCs. Thus, to make this model realistic, the un differentiated state (the equilibrium point that falls on the line y = x above) should be less stable than the differentiated states S 1 and S 2. The parameter set above does not fulfill this criterion, as a transient perturbation in X 1 of over 1 nm is required to differentiate from to S 1, though a perturbation of approximately. 4nM is sufficient to revert back to. 3.1 Decreasing Stability To rectify this, we define a new set of parameters with a shallower potential well for. We changed only a.
4 a = 1.6 b =. 5 K = 1 n = m = 2 δ = 1 Figure 3: Nullclines for new parameter values. Note that the perturbation 5 equilibria, three of which are stable. Note that the necessary perturbation to differentiate is now much smaller, and the stability of differentiated states is far greater, requiring a perturbation of up to. 8nM. The state occurs at concentrations X 1 = X 2 = 1.1, S 1 occurs at [ X 1 = 1.65, X 2 = 0.2], and S 2 occurs at [ X 1 = 0.2, X 2 = 1.65]. By simulating the system with these parameters and introducing Gaussian white noise, we were able to see spontaneous differentiation occur. Here are plots resulting from the simulation with a intensity I 0 = 1 and standard deviation σ = 2. Figure 4: Time series data of and according to the model with the described X 1 X 2 parameters and the introduction of randomly generated Gaussian noise.
5 We see that, from an initial condition of, the system spontaneously differentiates into S 1 or S 2 within tens of hours, after which is reliably stays at the differentiated equilibrium for over 100 hours. 4 Reprogramming Strategies Now that the parameterized model is sufficiently reflective of reality, we begin looking into strategies for reverting differentiated cells back to reliably and easily. We propose two strategies. 4.1 Tuning cooperativity It is possible to consistents un differentiate a cell in state or by transiently decreasing the cooperativity of the positive autoregulation. These exponents are proportional to the number of ligand binding sites on the X i promoter available for transcription factor. Thus, this value can be changed by designing a transcription factor N that binds to the domain but activates protein expression to a lesser degree than a bound X 1 protein. In our lumped model, this would decrease the net cooperativity, as there are fewer available sites for X 1 to bind to and activate its own promoter. By delivering a transient pulse of N that degrades over time, we can temporarily change the cooperativity then have it return to it s default value of 2. We show the results of a nullcline stability analysis of the system for different values of cooperativity. S 1 S 2 Figure 5: Nullclines for cooperativity of 1.67 (red) and 1.87 (green). Notice how the red nullcline is more linear than the green, eliminating the possibility of tri stability. At 1.87, we see the green nullcline is tangent to itself: the onset of tri stability in the system.
6 As the system moves from n = 1.67 through n = 1.86 and back to the natural value n = 2, the operating point of the system returns back to before noise induced fluctuations again cause differentiation to occur. The behavior of the system as it experienced this transient change in cooperativity is shown below. Figure 6: The time series data showing the concentrations of X 1 and X 2. The system was initialized to S 1 and experienced a transient pulse of neutral transcription factor N that decreased the effective cooperativity n eff to 1.5 at t = 100 hrs. As N was degraded and diluted, the deviation n died off exponentially with half life of As the system underwent a transient decrease in n eff 7 hours. n i eff, note how the system was briefly restored to for approximately 24 hours before redifferentiating. In the first two trials, the system returned to its original equilibrium, and in the third trial the system switched states. 4.2 Transient overexpression The above methodology required protein engineering, which is notoriously difficult. To make the process of reprogramming more easily implemented, we describe an alternative method involving transient overexpression of X 1 and X 2. This can be implemented by introducing mrna that encodes these proteins. Over time the mrna will be degraded and diluted, but not before being many rounds of translation. This results in an exponentially decaying transient pulse of protein to the system. In our model, this is represented as a transient increase in our basal production rates a 0 = a 01 = a 02.
7 We have conducted nullcline analysis on the system at different levels of basal expression. a 0 Figure 7: Nullclines of the system at different basal production rates. At a 0 =. 1, the system is monostable, thus, for sufficiently large basal expression, the operating point of the system will migrate back towards an un differentiated state. At a 0 =. 05, we see the onset of tri stability, as before in the cooperativity case. At a basal rate of zero, we return to the original nullclines of our system. We consider a transient pulse of mrna delivered at time t = 100 that translates into a spike in a 0 of. 1 that dies off with half life of 34 hours. Below are the results of a simulation without noise. Figure 8: un differentiation through transient overexpression of X 1 and X 2. The system becomes temporarily monostable. We also simulated the system with Gaussian noise, again with intensity standard deviation σ = 2. I 0 = 1 and
8 Figure 9: un differentiation through transient overexpression of X 1 and X 2. The system becomes temporarily monostable, resulting in un differentiation in spite of the strong stochastic noise. This demonstrates the feasibility of reverting a cell back to through transient overexpression that can be implemented via any standard transfection mechanism. 5 Stability with Feedback One problem remains to be solved. Though we have successfully reverted the cells to an un differentiated state, they will still spontaneously re differentiate from stochastic fluctuations in concentration. An additional mechanism must be put in place to enhance the stability properties of relative to and. There are many mechanisms that can accomplish this without dramatically changing the absolute values of. We consider two: genome editing with Crispr and high gain feedback. 5.1 Genome editing There are myriad ways to use genome editing with Crispr to change the parameters such that enjoys stronger stability properties. As we mentioned in Section 2.2, increasing or increases the stability of the un differentiated state. We can use a b Crispr to identify the and promoters, after which we can excise and insert X 1 X 2 sections of DNA. We can modify the and parameters by inserting stronger a b ligand binding sites to increase the strength and duration of the positive autoregulation and repression effects. S 1 S 2
9 Alternatively, we can insert a strand of DNA into the genome that constitutively produces X 1 and X 2. We conducted nullcline analysis for a basal expression rate of 0.4nM/hr. Recall from the previous section that. 5nM/hr is the onset of tri stability, and as the basal expression rate decreases, gets progressively less stable. Figure 10: Increasing the stability of through Crispr mediated constitutive basal expression of X 1 and X 2. Note that large perturbations are required to knock the system from to a differentiated state. All these genomic changes can be reversed or mitigated through another application of Crispr designed to excise the inserted DNA and return the genome to its unmodified state. However, this is difficult to implement with high fidelity because Crispr has poor success rates and may have unanticipated effects on the cells viability or ability to re differentiate. 5.2 High gain Feedback We can implement a high gain feedback loop to amplify the stability of. Using an understood mechanism such as phosphorylation as a source of gain in the system, we can implement a feedback loop as shown in Figure 11. Figure 11: High gain feedback loop. G represents a gain in the system, implemented transcriptional activation or phosphorylation processes.
10 The output y of this system is the concentration of X i. By using a sufficiently high 1 gain, we can have the output concentration approach K and experience arbitrarily small impact from noise and other effects. By tuning K, we can implement the system such that is the feedback loop s desired value. Any perturbation from this value will result in a large corrective effect, conferring strong stability properties to the un differentiated cells. 6 Conclusion We have described here a realistic two dimensional model describing the process of cell differentiation and re programming. Based on a Hill model of repression and positive autoregulation, we derived a system of ODEs that could be tuned to demonstrate spontaneous differentiation, as seen in true mammalian cells. We then describe two strategies for reverting these cells to an un differentiated state: by decreasing cooperativity and by transiently overexpressing the two proteins X 1 and X 2. We proceed to describe strategies of stabilizing this un differentiated state through genome editing and through the introduction of a high gain negative feedback loop. We hope this model can be generalized to the four dimensions indicated by the Yamanaka factors and shed light on the complicated processes at play during the creation of induced pluripotent stem cells.
7.32/7.81J/8.591J: Systems Biology. Fall Exam #1
7.32/7.81J/8.591J: Systems Biology Fall 2013 Exam #1 Instructions 1) Please do not open exam until instructed to do so. 2) This exam is closed- book and closed- notes. 3) Please do all problems. 4) Use
More informationPrinciples of Synthetic Biology: Midterm Exam
Principles of Synthetic Biology: Midterm Exam October 28, 2010 1 Conceptual Simple Circuits 1.1 Consider the plots in figure 1. Identify all critical points with an x. Put a circle around the x for each
More informationLecture 7: Simple genetic circuits I
Lecture 7: Simple genetic circuits I Paul C Bressloff (Fall 2018) 7.1 Transcription and translation In Fig. 20 we show the two main stages in the expression of a single gene according to the central dogma.
More informationBi 8 Lecture 11. Quantitative aspects of transcription factor binding and gene regulatory circuit design. Ellen Rothenberg 9 February 2016
Bi 8 Lecture 11 Quantitative aspects of transcription factor binding and gene regulatory circuit design Ellen Rothenberg 9 February 2016 Major take-home messages from λ phage system that apply to many
More informationBasic Synthetic Biology circuits
Basic Synthetic Biology circuits Note: these practices were obtained from the Computer Modelling Practicals lecture by Vincent Rouilly and Geoff Baldwin at Imperial College s course of Introduction to
More informationL3.1: Circuits: Introduction to Transcription Networks. Cellular Design Principles Prof. Jenna Rickus
L3.1: Circuits: Introduction to Transcription Networks Cellular Design Principles Prof. Jenna Rickus In this lecture Cognitive problem of the Cell Introduce transcription networks Key processing network
More informationCS-E5880 Modeling biological networks Gene regulatory networks
CS-E5880 Modeling biological networks Gene regulatory networks Jukka Intosalmi (based on slides by Harri Lähdesmäki) Department of Computer Science Aalto University January 12, 2018 Outline Modeling gene
More informationGenetic transcription and regulation
Genetic transcription and regulation Central dogma of biology DNA codes for DNA DNA codes for RNA RNA codes for proteins not surprisingly, many points for regulation of the process DNA codes for DNA replication
More information56:198:582 Biological Networks Lecture 9
56:198:582 Biological Networks Lecture 9 The Feed-Forward Loop Network Motif Subgraphs in random networks We have discussed the simplest network motif, self-regulation, a pattern with one node We now consider
More informationGenetic transcription and regulation
Genetic transcription and regulation Central dogma of biology DNA codes for DNA DNA codes for RNA RNA codes for proteins not surprisingly, many points for regulation of the process https://www.youtube.com/
More information56:198:582 Biological Networks Lecture 8
56:198:582 Biological Networks Lecture 8 Course organization Two complementary approaches to modeling and understanding biological networks Constraint-based modeling (Palsson) System-wide Metabolism Steady-state
More informationBiomolecular Feedback Systems
Biomolecular Feedback Systems Domitilla Del Vecchio MIT Richard M. Murray Caltech Version 1.0b, September 14, 2014 c 2014 by Princeton University Press All rights reserved. This is the electronic edition
More informationWritten Exam 15 December Course name: Introduction to Systems Biology Course no
Technical University of Denmark Written Exam 15 December 2008 Course name: Introduction to Systems Biology Course no. 27041 Aids allowed: Open book exam Provide your answers and calculations on separate
More informationIs chaos possible in 1d? - yes - no - I don t know. What is the long term behavior for the following system if x(0) = π/2?
Is chaos possible in 1d? - yes - no - I don t know What is the long term behavior for the following system if x(0) = π/2? In the insect outbreak problem, what kind of bifurcation occurs at fixed value
More informationBME 5742 Biosystems Modeling and Control
BME 5742 Biosystems Modeling and Control Lecture 24 Unregulated Gene Expression Model Dr. Zvi Roth (FAU) 1 The genetic material inside a cell, encoded in its DNA, governs the response of a cell to various
More informationBioinformatics: Network Analysis
Bioinformatics: Network Analysis Kinetics of Gene Regulation COMP 572 (BIOS 572 / BIOE 564) - Fall 2013 Luay Nakhleh, Rice University 1 The simplest model of gene expression involves only two steps: the
More informationA synthetic oscillatory network of transcriptional regulators
A synthetic oscillatory network of transcriptional regulators Michael B. Elowitz & Stanislas Leibler, Nature, 403, 2000 igem Team Heidelberg 2008 Journal Club Andreas Kühne Introduction Networks of interacting
More informationCombined Model of Intrinsic and Extrinsic Variability for Computational Network Design with Application to Synthetic Biology
Combined Model of Intrinsic and Extrinsic Variability for Computational Network Design with Application to Synthetic Biology The MIT Faculty has made this article openly available. Please share how this
More informationSparse, stable gene regulatory network recovery via convex optimization
Sparse, stable gene regulatory network recovery via convex optimization Arwen Meister June, 11 Gene regulatory networks Gene expression regulation allows cells to control protein levels in order to live
More informationSupplementary Materials for
www.sciencesignaling.org/cgi/content/full/7/345/ra91/dc1 Supplementary Materials for TGF-β induced epithelial-to-mesenchymal transition proceeds through stepwise activation of multiple feedback loops Jingyu
More informationLipniacki 2004 Ground Truth
Abstract Lipniacki 2004 Ground Truth The two-feedback-loop regulatory module of nuclear factor kb (NF-kB) signaling pathway is modeled by means of ordinary differential equations. signaling pathway: https://en.wikipedia.org/wiki/signaling_pathway
More informationBE 150 Problem Set #2 Issued: 16 Jan 2013 Due: 23 Jan 2013
M. Elowitz and R. M. Murray Winter 2013 CALIFORNIA INSTITUTE OF TECHNOLOGY Biology and Biological Engineering (BBE) BE 150 Problem Set #2 Issued: 16 Jan 2013 Due: 23 Jan 2013 1. (Shaping pulses; based
More informationSUPPLEMENTARY MATERIAL
Activity [RLU] Activity [RLU] RAB:TALA TALA:RAB TALB:RAB RAB:TALB RAB:TALB TALB:RAB TALA:RAB RAB:TALA Activity [nrlu] Activity [nrlu] SUPPLEMENTARY MATERIAL a b 100 100 50 50 0 0 5 25 50 5 25 50 0 50 50
More informationNetworks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationREVIEW SESSION. Wednesday, September 15 5:30 PM SHANTZ 242 E
REVIEW SESSION Wednesday, September 15 5:30 PM SHANTZ 242 E Gene Regulation Gene Regulation Gene expression can be turned on, turned off, turned up or turned down! For example, as test time approaches,
More informationMeasuring TF-DNA interactions
Measuring TF-DNA interactions How is Biological Complexity Achieved? Mediated by Transcription Factors (TFs) 2 Regulation of Gene Expression by Transcription Factors TF trans-acting factors TF TF TF TF
More informationProblem Set 2. 1 Competitive and uncompetitive inhibition (12 points) Systems Biology (7.32/7.81J/8.591J)
Problem Set 2 1 Competitive and uncompetitive inhibition (12 points) a. Reversible enzyme inhibitors can bind enzymes reversibly, and slowing down or halting enzymatic reactions. If an inhibitor occupies
More informationBE 150 Problem Set #4 Issued: 24 April 2014 Due: May
M. Elowitz and R. M. Murray Spring 2014 CALIFORNIA INSTITUTE OF TECHNOLOGY Biology and Biological Engineering (BBE) BE 150 Problem Set #4 Issued: 24 April 2014 Due: May 1 2014 1. (Repressilator modeling
More informationIntroduction. Dagmar Iber Jörg Stelling. CSB Deterministic, SS 2015, 1.
Introduction Dagmar Iber Jörg Stelling joerg.stelling@bsse.ethz.ch CSB Deterministic, SS 2015, 1 Origins of Systems Biology On this assumption of the passage of blood, made as a basis for argument, and
More informationLecture 4: Transcription networks basic concepts
Lecture 4: Transcription networks basic concepts - Activators and repressors - Input functions; Logic input functions; Multidimensional input functions - Dynamics and response time 2.1 Introduction The
More informationModeling Multiple Steady States in Genetic Regulatory Networks. Khang Tran. problem.
Modeling Multiple Steady States in Genetic Regulatory Networks Khang Tran From networks of simple regulatory elements, scientists have shown some simple circuits such as the circadian oscillator 1 or the
More informationComputational Cell Biology Lecture 4
Computational Cell Biology Lecture 4 Case Study: Basic Modeling in Gene Expression Yang Cao Department of Computer Science DNA Structure and Base Pair Gene Expression Gene is just a small part of DNA.
More information(Refer Slide Time: 1:42)
Control Engineering Prof. Madan Gopal Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture - 21 Basic Principles of Feedback Control (Contd..) Friends, let me get started
More informationBistability in ODE and Boolean network models
Bistability in ODE and Boolean network models Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson)
More informationWelcome to Class 21!
Welcome to Class 21! Introductory Biochemistry! Lecture 21: Outline and Objectives l Regulation of Gene Expression in Prokaryotes! l transcriptional regulation! l principles! l lac operon! l trp attenuation!
More informationIntroduction. ECE/CS/BioEn 6760 Modeling and Analysis of Biological Networks. Adventures in Synthetic Biology. Synthetic Biology.
Introduction ECE/CS/BioEn 6760 Modeling and Analysis of Biological Networks Chris J. Myers Lecture 17: Genetic Circuit Design Electrical engineering principles can be applied to understand the behavior
More informationAP3162D: Lecture 2 - Gene-regulatory circuits: Deterministic dynamics
AP3162D: Lecture 2 - Gene-regulatory circuits: Deterministic dynamics Hyun Youk Delft University of Technology (Dated: February 22, 2018) In this lecture, we will discuss how to model standard gene-regulatory
More informationStochastic simulations
Stochastic simulations Application to molecular networks Literature overview Noise in genetic networks Origins How to measure and distinguish between the two types of noise (intrinsic vs extrinsic)? What
More information1 Periodic stimulations of the Incoherent Feedforward Loop network
1 Periodic stimulations of the Incoherent Feedforward Loop network In this Additional file, we give more details about the mathematical analysis of the periodic activation of the IFFL network by a train
More informationMulti-modality in gene regulatory networks with slow promoter kinetics. Abstract. Author summary
Multi-modality in gene regulatory networks with slow promoter kinetics M. Ali Al-Radhawi 1, D. Del Vecchio 2, E. D. Sontag 3*, 1 Department of Electrical and Computer Engineering, Northeastern University,
More informationStochastic and delayed stochastic models of gene expression. and regulation
*Manuscript Stochastic and delayed stochastic models of gene expression and regulation Andre S. Ribeiro a,b a Computational Systems Biology Research Group, Dept. of Signal Processing, Tampere University
More informationRegulation and signaling. Overview. Control of gene expression. Cells need to regulate the amounts of different proteins they express, depending on
Regulation and signaling Overview Cells need to regulate the amounts of different proteins they express, depending on cell development (skin vs liver cell) cell stage environmental conditions (food, temperature,
More informationThe Effect of Stochasticity on the Lac Operon: An Evolutionary Perspective
The Effect of Stochasticity on the Lac Operon: An Evolutionary Perspective Milan van Hoek *, Paulien Hogeweg Theoretical Biology/Bioinformatics Group, Utrecht University, Utrecht, The Netherlands The role
More informationDynamical-Systems Perspective to Stem-cell Biology: Relevance of oscillatory gene expression dynamics and cell-cell interaction
Dynamical-Systems Perspective to Stem-cell Biology: Relevance of oscillatory gene expression dynamics and cell-cell interaction Kunihiko Kaneko Universal Biology Inst., Center for Complex-Systems Biology,
More informationChapter 15 Active Reading Guide Regulation of Gene Expression
Name: AP Biology Mr. Croft Chapter 15 Active Reading Guide Regulation of Gene Expression The overview for Chapter 15 introduces the idea that while all cells of an organism have all genes in the genome,
More informationCONTROL OF STOCHASTIC GENE EXPRESSION BY A NONLINEAR BIOLOGICAL OSCILLATOR
PHYSCON 27, Florence, Italy, 7 9 July, 27 CONTROL OF STOCHASTIC GENE EXPRESSION BY A NONLINEAR BIOLOGICAL OSCILLATOR Federica Colombo Division of Genetics and Cell Biology San Raffaele Scientific Institute
More informationProkaryotic Regulation
Prokaryotic Regulation Control of transcription initiation can be: Positive control increases transcription when activators bind DNA Negative control reduces transcription when repressors bind to DNA regulatory
More informationNewly made RNA is called primary transcript and is modified in three ways before leaving the nucleus:
m Eukaryotic mrna processing Newly made RNA is called primary transcript and is modified in three ways before leaving the nucleus: Cap structure a modified guanine base is added to the 5 end. Poly-A tail
More informationLecture 6: The feed-forward loop (FFL) network motif
Lecture 6: The feed-forward loop (FFL) network motif Chapter 4 of Alon x 4. Introduction x z y z y Feed-forward loop (FFL) a= 3-node feedback loop (3Loop) a=3 Fig 4.a The feed-forward loop (FFL) and the
More informationThe role of mrna and protein stability in the function of coupled positive and negative feedback systems in eukaryotic cells.
The role of mrna and protein stability in the function of coupled positive and negative feedback systems in eukaryotic cells. Kristian Moss Bendtsen 1, Mogens H. Jensen 1a, Sandeep Krishna 1,2, Szabolcs
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature13238 Additional Experimental Results Degradation tag experiments In addition to exploring the effect of variable-length linker (TS repeats) on the phase-shift in module degradation (Supplementary
More informationthe noisy gene Biology of the Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB)
Biology of the the noisy gene Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB) day III: noisy bacteria - Regulation of noise (B. subtilis) - Intrinsic/Extrinsic
More informationSimulation of Gene Regulatory Networks
Simulation of Gene Regulatory Networks Overview I have been assisting Professor Jacques Cohen at Brandeis University to explore and compare the the many available representations and interpretations of
More informationGene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back
Gene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back June 19, 2007 Motivation & Basics A Stochastic Approach to Gene Expression Application to Experimental
More informationRegulation of Gene Expression
Chapter 18 Regulation of Gene Expression Edited by Shawn Lester PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley
More informationGene Autoregulation via Intronic micrornas and its Functions
Gene Autoregulation via Intronic micrornas and its Functions Carla Bosia Department of Theoretical Physics University of Torino and INFN, Italy cbosia@to.infn.it Annecy-le-Vieux 20-22/10/2010 OUTLINE microrna
More informationBistability and a differential equation model of the lac operon
Bistability and a differential equation model of the lac operon Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Fall 2016 M. Macauley
More informationModeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool
Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht
More informationBioinformatics 3. V18 Kinetic Motifs. Fri, Jan 8, 2016
Bioinformatics 3 V18 Kinetic Motifs Fri, Jan 8, 2016 Modelling of Signalling Pathways Curr. Op. Cell Biol. 15 (2003) 221 1) How do the magnitudes of signal output and signal duration depend on the kinetic
More informationA systems approach to biology
A systems approach to biology SB200 Lecture 5 30 September 2008 Jeremy Gunawardena jeremy@hms.harvard.edu Recap of Lecture 4 matrix exponential exp(a) = 1 + A + A2/2 +... + Ak/k! +... dx/dt = Ax matrices
More informationBioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014"
Bioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014" Modelling of Signalling Pathways" Curr. Op. Cell Biol. 15 (2003) 221" 1) How do the magnitudes of signal output and signal duration depend on the
More informationTesting the transition state theory in stochastic dynamics of a. genetic switch
Testing the transition state theory in stochastic dynamics of a genetic switch Tomohiro Ushikubo 1, Wataru Inoue, Mitsumasa Yoda 1 1,, 3, and Masaki Sasai 1 Department of Computational Science and Engineering,
More informationLaboratory Exercise 1 DC servo
Laboratory Exercise DC servo Per-Olof Källén ø 0,8 POWER SAT. OVL.RESET POS.RESET Moment Reference ø 0,5 ø 0,5 ø 0,5 ø 0,65 ø 0,65 Int ø 0,8 ø 0,8 Σ k Js + d ø 0,8 s ø 0 8 Off Off ø 0,8 Ext. Int. + x0,
More informationRegulated gene expression is the process through which cells
Noise-based switches and amplifiers for gene expression Jeff Hasty*, Joel Pradines*, Milos Dolnik*, and J. J. Collins* *Center for BioDynamics and Department of Biomedical Engineering, Boston University,
More informationA complementation test would be done by crossing the haploid strains and scoring the phenotype in the diploids.
Problem set H answers 1. To study DNA repair mechanisms, geneticists isolated yeast mutants that were sensitive to various types of radiation; for example, mutants that were more sensitive to UV light.
More informationIntroduction. Gene expression is the combined process of :
1 To know and explain: Regulation of Bacterial Gene Expression Constitutive ( house keeping) vs. Controllable genes OPERON structure and its role in gene regulation Regulation of Eukaryotic Gene Expression
More informationDynamics of the Mixed Feedback Loop Integrated with MicroRNA
The Second International Symposium on Optimization and Systems Biology (OSB 08) Lijiang, China, October 31 November 3, 2008 Copyright 2008 ORSC & APORC, pp. 174 181 Dynamics of the Mixed Feedback Loop
More informationHonors Biology Reading Guide Chapter 11
Honors Biology Reading Guide Chapter 11 v Promoter a specific nucleotide sequence in DNA located near the start of a gene that is the binding site for RNA polymerase and the place where transcription begins
More informationBasic Principles of Kinetics and Thermodynamics Arthur L. Haas, Ph.D. Department of Biochemistry and Molecular Biology
Basic Principles of Kinetics and Thermodynamics Arthur L. Haas, Ph.D. Department of Biochemistry and Molecular Biology Review: Garrett and Grisham, Enzyme Kinetics (Chapt. 13), in Biochemistry, Third Edition,
More informationTopic 4 - #14 The Lactose Operon
Topic 4 - #14 The Lactose Operon The Lactose Operon The lactose operon is an operon which is responsible for the transport and metabolism of the sugar lactose in E. coli. - Lactose is one of many organic
More informationMolecular Biology (9)
Molecular Biology (9) Translation Mamoun Ahram, PhD Second semester, 2017-2018 1 Resources This lecture Cooper, Ch. 8 (297-319) 2 General information Protein synthesis involves interactions between three
More informationComputational Models of the Notch Network Elucidate Mechanisms of Context-dependent Signaling
Computational Models of the Notch Network Elucidate Mechanisms of Context-dependent Signaling Smita Agrawal 1,2,3 a, Colin Archer 1,2,3 b, David V. Schaffer 1,2,3 * 1 Department of Chemical Engineering,
More informationProtein production. What can we do with ODEs
Protein production What can we do with ODEs 2018 Modeling protein production with ODE We have seen various examples of wiring networks. Once you know how a regulatory network is wired in a particular way,
More informationOld FINAL EXAM BIO409/509 NAME. Please number your answers and write them on the attached, lined paper.
Old FINAL EXAM BIO409/509 NAME Please number your answers and write them on the attached, lined paper. Gene expression can be regulated at several steps. Describe one example for each of the following:
More informationNature Genetics: doi: /ng Supplementary Figure 1. Plasmid diagrams of primary strains used in this work.
Supplementary Figure 1 Plasmid diagrams of primary strains used in this work. Refer to Supplementary Table 1 for additional information and full list of strains. Diagrams are not to scale. (a) Strain used
More information2 Dilution of Proteins Due to Cell Growth
Problem Set 1 1 Transcription and Translation Consider the following set of reactions describing the process of maing a protein out of a gene: G β g G + M M α m M β m M + X X + S 1+ 1 X S 2+ X S X S 2
More informationTopic 4: Equilibrium binding and chemical kinetics
Topic 4: Equilibrium binding and chemical kinetics Outline: Applications, applications, applications use Boltzmann to look at receptor-ligand binding use Boltzmann to look at PolII-DNA binding and gene
More informationNoisy Attractors and Ergodic Sets in Models. of Genetic Regulatory Networks
Noisy Attractors and Ergodic Sets in Models of Genetic Regulatory Networks Andre S. Ribeiro Institute for Biocomplexity and Informatics, Univ. of Calgary, Canada Department of Physics and Astronomy, Univ.
More informationStochastic Model for Adaptation Using Basin Hopping Dynamics
Stochastic Model for Adaptation Using Basin Hopping Dynamics Peter Davis NTT Communication Science Laboratories 2-4 Hikaridai, Keihanna Science City, Kyoto, Japan 619-0237 davis@cslab.kecl.ntt.co.jp Abstract
More informationarxiv:physics/ v1 [physics.bio-ph] 30 Mar 2000
arxiv:physics/0003105v1 [physics.bio-ph] 30 Mar 2000 Noise-Based Switches and Amplifiers for Gene Expression Jeff Hasty 1, Joel Pradines 1, Milos Dolnik 1,2 and J.J. Collins 1 September 23, 1999 1 Center
More informationProkaryotic Gene Expression (Learning Objectives)
Prokaryotic Gene Expression (Learning Objectives) 1. Learn how bacteria respond to changes of metabolites in their environment: short-term and longer-term. 2. Compare and contrast transcriptional control
More informationStochastic Process Algebra models of a Circadian Clock
Stochastic Process Algebra models of a Circadian Clock Jeremy T. Bradley Thomas Thorne Department of Computing, Imperial College London 180 Queen s Gate, London SW7 2BZ, United Kingdom Email: jb@doc.ic.ac.uk
More informationEffective interaction graphs arising from resource limitations in gene networks
Effective interaction graphs arising from resource limitations in gene networks The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
More informationLecture 4: Importance of Noise and Fluctuations
Lecture 4: Importance of Noise and Fluctuations Jordi Soriano Fradera Dept. Física de la Matèria Condensada, Universitat de Barcelona UB Institute of Complex Systems September 2016 1. Noise in biological
More informationComputational Systems Biology Exam
Computational Systems Biology Exam Dr. Jürgen Pahle Aleksandr Andreychenko, M.Sc. 31 July, 2012 Name Matriculation Number Do not open this exam booklet before we ask you to. Do read this page carefully.
More informationComplete all warm up questions Focus on operon functioning we will be creating operon models on Monday
Complete all warm up questions Focus on operon functioning we will be creating operon models on Monday 1. What is the Central Dogma? 2. How does prokaryotic DNA compare to eukaryotic DNA? 3. How is DNA
More information10-810: Advanced Algorithms and Models for Computational Biology. microrna and Whole Genome Comparison
10-810: Advanced Algorithms and Models for Computational Biology microrna and Whole Genome Comparison Central Dogma: 90s Transcription factors DNA transcription mrna translation Proteins Central Dogma:
More informationComputer Simulation and Data Analysis in Molecular Biology and Biophysics
Victor Bloomfield Computer Simulation and Data Analysis in Molecular Biology and Biophysics An Introduction Using R Springer Contents Part I The Basics of R 1 Calculating with R 3 1.1 Installing R 3 1.1.1
More informationBasic modeling approaches for biological systems. Mahesh Bule
Basic modeling approaches for biological systems Mahesh Bule The hierarchy of life from atoms to living organisms Modeling biological processes often requires accounting for action and feedback involving
More informationBi 1x Spring 2014: LacI Titration
Bi 1x Spring 2014: LacI Titration 1 Overview In this experiment, you will measure the effect of various mutated LacI repressor ribosome binding sites in an E. coli cell by measuring the expression of a
More information13.4 Gene Regulation and Expression
13.4 Gene Regulation and Expression Lesson Objectives Describe gene regulation in prokaryotes. Explain how most eukaryotic genes are regulated. Relate gene regulation to development in multicellular organisms.
More informationDynamics and Control of Synthetic Bio-molecular Networks
Dynamics and Control of Synthetic io-molecular Networks D. Del Vecchio and E. D. Sontag bstract This tutorial paper presents an introduction to systems and synthetic molecular biology. It provides an introduction
More informationTuning an Activator-Repressor Clock Employing Retroactivity
Tuning an Activator-Repressor Clock Employing Retroactivity Alexander Rosenberg, Shridhar Jayanthi and Domitilla Del Vecchio Electrical Engineering Dept., University of Washington, Seattle WA 9895 Electrical
More informationCybergenetics: Control theory for living cells
Department of Biosystems Science and Engineering, ETH-Zürich Cybergenetics: Control theory for living cells Corentin Briat Joint work with Ankit Gupta and Mustafa Khammash Introduction Overview Cybergenetics:
More informationarxiv: v2 [q-bio.qm] 12 Jan 2017
Approximation and inference methods for stochastic biochemical kinetics - a tutorial review arxiv:1608.06582v2 [q-bio.qm] 12 Jan 2017 David Schnoerr 1,2,3, Guido Sanguinetti 2,3, and Ramon Grima 1,3,*
More informationBrief contents. Chapter 1 Virus Dynamics 33. Chapter 2 Physics and Biology 52. Randomness in Biology. Chapter 3 Discrete Randomness 59
Brief contents I First Steps Chapter 1 Virus Dynamics 33 Chapter 2 Physics and Biology 52 II Randomness in Biology Chapter 3 Discrete Randomness 59 Chapter 4 Some Useful Discrete Distributions 96 Chapter
More information2. Mathematical descriptions. (i) the master equation (ii) Langevin theory. 3. Single cell measurements
1. Why stochastic?. Mathematical descriptions (i) the master equation (ii) Langevin theory 3. Single cell measurements 4. Consequences Any chemical reaction is stochastic. k P d φ dp dt = k d P deterministic
More informationChapter 18: Control of Gene Expression
Chapter 18: Control of Gene Expression 海洋生物研究所 曾令銘 海事大樓 426 室分機 : 5326 Differential Expression of Genes Prokaryotes and eukaryotes precisely regulate gene expression in response to environmental conditions
More informationA(γ A D A + γ R D R + γ C R + δ A )
Title: Mechanisms of noise-resistance in genetic oscillators Authors: José M. G. Vilar 1,2, Hao Yuan Kueh 1, Naama Barkai 3, and Stanislas Leibler 1,2 1 Howard Hughes Medical Institute, Departments of
More informationProblem Set 5. 1 Waiting times for chemical reactions (8 points)
Problem Set 5 1 Waiting times for chemical reactions (8 points) In the previous assignment, we saw that for a chemical reaction occurring at rate r, the distribution of waiting times τ between reaction
More information