Registry of Standard Biological Models/CellML Practical/Practical 1

From OpenWetWare
Jump to navigationJump to search

Goal: Describe a simple genetic assembly (promoter + RBS + protein coding region + stop codon) from the description of each sub-component.

Example of such a part: <bbpart>BBa_I7101</bbpart>

Simple Modeling

Assumptions:
> Central Dogma: Gene --> mRNA --> Protein
> unlimited resources to produce mRNA and proteins (i.e. excess of Polymerases, Ribosomes, tRNA ...)
> No stochastic effect taken into account, continous approach

<amsmath>\frac{d[mRNA]}{dt} = k_M*[Gene] - \gamma_M*[mRNA]</amsmath>

<amsmath>\frac{d[Protein]}{dt} = k_P*[mRNA] - \gamma_P*[Protein]</amsmath>

Where,

<amsmath>\emph [Gene]</amsmath>= number of gene copies per cell, [Gene]
<amsmath>\emph [mRNA]</amsmath>= mRNA concentration per cell, [mRNA]
<amsmath>\emph [Protein]</amsmath>= Protein concentration per cell, [Protein]
<amsmath>\emph k_{M}</amsmath>= mRNA production rate (PoPS), unit:[mRNA]/[DNA]/s
<amsmath>\emph k_{P}</amsmath>= Protein production rate (RiPS), unit:[Protein]/[mRNA]/s
<amsmath>\gamma_M</amsmath>= mRNA degradation rate, 1/s
<amsmath>\gamma_P</amsmath>= Protein degradation rate., 1/s

BioBrick matching

Simple BioBrick: BBa_I7101
Component Function Related parameters Comments
Promoter tanscription <amsmath>\emph [Gene], \emph k_{M} </amsmath> the promoter can be characterized with a given rate of transcription rate (PoPs) which depends on its intrinsic nature and the number of gene copies available in the system
RBS translation <amsmath>\emph k_{P}</amsmath> the RBS defines the rate of translation per mRNA
mRNA degradation <amsmath>\gamma_M</amsmath> the mRNA component is not a BioBrick, however it has a specific degradation rate
Protein degradation <amsmath>\gamma_P</amsmath> without considering its function, the protein should be characterized by its half-life.

CellML matching

Conversion of BBa_I into a CellML formalism
Detail view
Abstraction)

>>>>Comments

  • About the Promoter:
    • Characterized by nb gene copies + transcription rate
    • PoPs is defined as a MathML expression. Can easily define repression, activation or constitutive expression behaviour.
    • Inputs: none. BC: Should numbers of copies be here? VR: true.
    • Output: PoPs (Polymerases per second). It is used to feed a Protein component. BC: Although it feeds into a protein component here, more generally it could feed into other components, such as tRNA. VR: Not too sure to understand how you define PoPs for tRNAs. Can you explain ? BC:I mean that PoPS should be independent of whatever piece of DNA you are sitting on when you measure it. So the polymerases could be transcribing an mRNA or a tRNA or whatever. Since we can easily imagine wanting to transcribe non-protein components such as tRNA, we need to keep the output of the promoter general in my opinion. VR: I am a bit confused. To me polymerases are only binding DNA (chromosomal or plasmid). tRNA are compounds formed by an anti-codon and amino-acid site, to me they onlyinteract with the ribosome during translation of the protein. Am I missing something ?
  • About the RBS:
    • Characterized by translation rate
    • RiPS is defined by a MathML expression. Could describe limited resource at this point.
    • Inputs: needs [mRNA]
    • Ouputs: RiPs (Ribosomes per second). It is used to feed a Protein component.


  • About the Protein:
  • I felt the need to encapsulate in the same component the concept of mRNA and protein. mRNA is not a BioBrick anyway. This component might benefit being implemented as a group, or maybe simply an import of a 'protein_molecule' component and a 'mRNA_molecule' component.
  • We should think about riboswitch parts which are affecting mRNA behaviour
  • Characterized by degradation rate of mRNA and Protein
  • Inputs: need PoPs from a promoter + RiPs from a RBS
    • BC:Do we really need PoPS as an input to a protein component. I feel that the PoPS defines the numbers of mRNA and after that, it doesn't really affect the protein component. What do you think? VR: I have to admit that I haven't found yet a design with which I am happy. My problem comes from the fact that the mRNA component is not a BioBrick, so I would like to hide it some how. I have also tried to encapsulate it with the RBS component here.
  • Outputs: [mRNA] + [Protein]
  • The overall part:
    • Inputs: none BC:Number of DNA copies of the part? VR: true.
    • Outputs: Protein and mRNA expression levelss
    • Import all the other components: Promoter + RBS + Protein