leloup2

reaction_0

Mb synthesized

∅ > species_0

reaction_1

Mb translated into protein

∅ > species_1

reaction_10

Mp translated into protein

∅ > species_8

reaction_11

Pcp specific degradation

species_9 > ∅

reaction_12

Pc phosphorylation

species_8 > species_9

reaction_13

Cc and Pc produce complex Per_Cry

species_4 + species_8 = species_10

reaction_14

PCc phosphorylation

species_10 > species_11

reaction_15

PCcp specific degradation

species_11 > ∅

reaction_16

PCc transfered into nuclear

species_10 = species_12

reaction_17

PCnp nonspecific degradation

species_14 > ∅

reaction_18

Bcp nonspecific degradation

species_2 > ∅

reaction_19

Bnp nonspecific degradation

species_13 > ∅

reaction_2

Mb nonspecific degradation

species_0 > ∅

reaction_20

Mc synthesis

∅ > species_5

reaction_21

PCn phosphorylation

species_12 > species_14

reaction_22

Mp nonspecific degradation

species_7 > ∅

reaction_23

Per_Cry and Clock_Bmal form inactive complex

species_12 + species_3 = species_15

reaction_24

Mb specific degradation

species_0 > ∅

reaction_25

Mc specific degradation

species_5 > ∅

reaction_26

Mp specific degradation

species_7 > ∅

reaction_27

Pc nonspecific degradation

species_8 > ∅

reaction_28

Cc nonspecific degradation

species_4 > ∅

reaction_29

Pcp nonspecific degradation

species_9 > ∅

reaction_3

Bc phosphorylation

species_1 > species_2

reaction_30

Ccp nonspecific degradation

species_6 > ∅

reaction_31

PCcp nonspecific degradation

species_11 > ∅

reaction_32

PCc nonspecific degradation

species_10 > ∅

reaction_33

PCnp specific degradation

species_14 > ∅

reaction_34

Bc nonspecific degradation

species_1 > ∅

reaction_35

Bcp specific degradation

species_2 > ∅

reaction_36

Bn phosphorylation

species_3 > species_13

reaction_37

Bnp specific degradation

species_13 > ∅

reaction_38

In nonspecific degration

species_15 > ∅

reaction_39

In specific degradation

species_15 > ∅

reaction_4

Bc transfered from cytosolic to nuclear

species_1 = species_3

reaction_40

Bn nonspecific degradation

species_3 > ∅

reaction_41

Bcp dephosphorylation

species_2 > species_1

reaction_42

Bnp dephosphorylation

species_13 > species_3

reaction_43

Ccp dephosphorylation

species_6 > species_4

reaction_44

Pcp dephosphorylation

species_9 > species_8

reaction_45

PCnp dephosphorylation

species_14 > species_12

reaction_46

PCn nonspecific degradation

species_12 > ∅

reaction_47

PCcp dephosphorylation

species_11 > species_10

reaction_5

Mc translated into protein

∅ > species_4

reaction_6

Mc nonspecific degradation

species_5 > ∅

reaction_7

Cc phosphorylation

species_4 > species_6

reaction_8

Ccp specific degradation

species_6 > ∅

reaction_9

Mp synthesis

∅ > species_7

Global parameters
reaction_0
reaction_1
reaction_10
reaction_11
reaction_12
reaction_13
reaction_14
reaction_15
reaction_16
reaction_17
reaction_18
reaction_19
reaction_2
reaction_20
reaction_21
reaction_22
reaction_23
reaction_24
reaction_25
reaction_26
reaction_27
reaction_28
reaction_29
reaction_3
reaction_30
reaction_31
reaction_32
reaction_33
reaction_34
reaction_35
reaction_36
reaction_37
reaction_38
reaction_39
reaction_4
reaction_40
reaction_41
reaction_42
reaction_43
reaction_44
reaction_45
reaction_46
reaction_47
reaction_5
reaction_6
reaction_7
reaction_8
reaction_9

Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.


Species:

Reactions:


Middle-click: pin/unpin nodes
Shift-click: pool/unpool species
Right-click: context menu

Apply alternate model layout to overlapping elements in current model:

log scales

y-axis min/max

x-axis min/max

Toward a detailed computational model for the mammalian circadian clock.

  • Jean-Christophe Leloup
  • Albert Goldbeter
Proc. Natl. Acad. Sci. U.S.A. 2003; 100 (12): 7051-7056
Abstract
We present a computational model for the mammalian circadian clock based on the intertwined positive and negative regulatory loops involving the Per, Cry, Bmal1, Clock, and Rev-Erb alpha genes. In agreement with experimental observations, the model can give rise to sustained circadian oscillations in continuous darkness, characterized by an antiphase relationship between Per/Cry/Rev-Erbalpha and Bmal1 mRNAs. Sustained oscillations correspond to the rhythms autonomously generated by suprachiasmatic nuclei. For other parameter values, damped oscillations can also be obtained in the model. These oscillations, which transform into sustained oscillations when coupled to a periodic signal, correspond to rhythms produced by peripheral tissues. When incorporating the light-induced expression of the Per gene, the model accounts for entrainment of the oscillations by light-dark cycles. Simulations show that the phase of the oscillations can then vary by several hours with relatively minor changes in parameter values. Such a lability of the phase could account for physiological disorders related to circadian rhythms in humans, such as advanced or delayed sleep phase syndrome, whereas the lack of entrainment by light-dark cycles can be related to the non-24h sleep-wake syndrome. The model uncovers the possible existence of multiple sources of oscillatory behavior. Thus, in conditions where the indirect negative autoregulation of Per and Cry expression is inoperative, the model indicates the possibility that sustained oscillations might still arise from the negative autoregulation of Bmal1 expression.
The SBML for this model was obtained from the BioModels database (BioModels ID: BIOMD0000000073) Biomodels notes: This model is according to the paper Toward a detailed computational model for the mammalian circadian clock . In this model only interlocked negative and positive regulation of Per, Cry, Bmal gene are involved. Some initial values were not provided, therefore they were chosen to fit the curve from the paper. JWS Online curation: This model was curated by reproducing the figures as described in the BioModels Notes. No additional changes were made.