User:Johnsy/Lipoprotein Modelling/Model Development: Difference between revisions

From OpenWetWare
Jump to navigationJump to search
No edit summary
No edit summary
Line 1: Line 1:
==A Very Simple Model==
===Introduction===
To start exploring the method which is used to formulate models, we consider a very simple model derived from the biology of LDL particles in Brown & Goldstein, 1979.  We are assuming a cell with a fixed number of LDL particles already bound to the surface of the cell.  The rate of internalization is only dependent upon the number of LDL particles already on the surface.  Futhermore, once the LDL particle is internalized, it is degraded by the lysosome into its constituent components (ie cholesterol esters, phospholipids, apoproteins, etc.). 
We only consider three variables:
#The concentration of surface bound LDL particles
#The concentration of internalized LDL particles
#The concentration of the degradation products of LDL particles (directly proportional to the intracellular cholesterol levels
===Developing the Model===
Equation 1:
<center><math>
\frac{d[SB]}{dt} = -k_{I}[SB]
</math></center>
In the above equation, there are no inputs as we are assuming that the LDL particles are already bound to the surface of the cell.  We have also only one sink which is due to the internalization of LDL particles, determined by the rate of internalization, ''k<sub>I</sub>'', and the concentration of LDL particles bound to the surface of the cell.
Equation 2:
<center><math>
\frac{d[INT]}{dt} = k_{I}[SB] - k_{deg}[INT]
</math></center>
In the above equation 2, we now have a source which is the LDL particles which were internalized from equation 1.  We also have a sink which is the degradation of the internalized LDL particles deteremined by rate ''k<sub>deg</sub>'' and the concentration of internalzed LDL particles.
Equation 3:
<center><math>
\frac{d[DEG]}{dt} = k_{deg}[INT]
</math></center>
In the above equation 3, the degradation products are produced from the LDL particles which were internalized.  We are not considering the fate of these degradation products, but we are assuming that they stay in the cell and accumulate.  The primary degradation products are the cholesterol, fatty acids, lipoproteins, and proteins contained in the LDL particle and these are further processed and recycled by the cell.  The choelsterol is converted to bile acid and cholesterol derivatives such as steriod hormones, while the fatty acids and proteins are broken down or stored for use by other cellular pathways.  The lipoproteins are also recycled or are degraded by the lysosome.
===Parameters and Initial Conditions===
From the Panovska paper discussed below, we obtain the two parameters from literature that govern the above model.
#''k<sub>I</sub>'' - the rate of internalzation of LDL particles = log(2)/22 s<sup>-1</sup>
#''k<sub>deg</sub>'' - the rate of degradation of LDL particles = 1/180 s<sup>-1</sup>
The initial conditions we enter into our model are as follows:
#[SB]<sub>0</sub> - the initial concentration of surface bound LDL particles we set arbitrarily at 45
#[INT]<sub>0</sub> - the initial concentration of internalized LDL particles we initial set to zero (the cells are starved of cholesterol)
#[DEG]<sub>0</sub> - the initial concentration of degradation products we also set to zero since we assume that the cells have not had any cholesterol in them for a lengthy time period.
==Developing the Simple Model==
==HMG-CoA Reductase and Cholesterol==
==HMG-CoA Reductase and Cholesterol==



Revision as of 06:52, 1 November 2007

A Very Simple Model

Introduction

To start exploring the method which is used to formulate models, we consider a very simple model derived from the biology of LDL particles in Brown & Goldstein, 1979. We are assuming a cell with a fixed number of LDL particles already bound to the surface of the cell. The rate of internalization is only dependent upon the number of LDL particles already on the surface. Futhermore, once the LDL particle is internalized, it is degraded by the lysosome into its constituent components (ie cholesterol esters, phospholipids, apoproteins, etc.).

We only consider three variables:

  1. The concentration of surface bound LDL particles
  2. The concentration of internalized LDL particles
  3. The concentration of the degradation products of LDL particles (directly proportional to the intracellular cholesterol levels

Developing the Model

Equation 1:

[math]\displaystyle{ \frac{d[SB]}{dt} = -k_{I}[SB] }[/math]

In the above equation, there are no inputs as we are assuming that the LDL particles are already bound to the surface of the cell. We have also only one sink which is due to the internalization of LDL particles, determined by the rate of internalization, kI, and the concentration of LDL particles bound to the surface of the cell.

Equation 2:

[math]\displaystyle{ \frac{d[INT]}{dt} = k_{I}[SB] - k_{deg}[INT] }[/math]

In the above equation 2, we now have a source which is the LDL particles which were internalized from equation 1. We also have a sink which is the degradation of the internalized LDL particles deteremined by rate kdeg and the concentration of internalzed LDL particles.

Equation 3:

[math]\displaystyle{ \frac{d[DEG]}{dt} = k_{deg}[INT] }[/math]

In the above equation 3, the degradation products are produced from the LDL particles which were internalized. We are not considering the fate of these degradation products, but we are assuming that they stay in the cell and accumulate. The primary degradation products are the cholesterol, fatty acids, lipoproteins, and proteins contained in the LDL particle and these are further processed and recycled by the cell. The choelsterol is converted to bile acid and cholesterol derivatives such as steriod hormones, while the fatty acids and proteins are broken down or stored for use by other cellular pathways. The lipoproteins are also recycled or are degraded by the lysosome.

Parameters and Initial Conditions

From the Panovska paper discussed below, we obtain the two parameters from literature that govern the above model.

  1. kI - the rate of internalzation of LDL particles = log(2)/22 s-1
  2. kdeg - the rate of degradation of LDL particles = 1/180 s-1

The initial conditions we enter into our model are as follows:

  1. [SB]0 - the initial concentration of surface bound LDL particles we set arbitrarily at 45
  2. [INT]0 - the initial concentration of internalized LDL particles we initial set to zero (the cells are starved of cholesterol)
  3. [DEG]0 - the initial concentration of degradation products we also set to zero since we assume that the cells have not had any cholesterol in them for a lengthy time period.

Developing the Simple Model

HMG-CoA Reductase and Cholesterol

Cholesterol Producers (Sources)

  • de novo - Acetyl CoA via HMG CoA and enzyme HMG CoA reductase
  • IDL particles accepted into the cell via LDL receptors
  • LDL particles accepted into the cell via LDL receptors

Cholesterol Fates (Sinks)

  • Bile acids via enzyme cholesterol 7-α-hydroxylase

de novo Pathway

HMG CoA reductase is the rate limiting enzyme in the biosynthesis of cholesterol. Its transcription is upregulated by the sterol regulatory element binding protein (SREBP) which binds to the streol regulatory element (SRE) to transcribe the gene for HMG CoA reductase. SREBP is usually situated on the membrane of the endoplasmic reticulum or the nuclear membrane, but when bound by cholesterol, the protein is released via proteolysis and migrates to the nucleus where it binds to the SRE to initiate transcription.

Statins limit the action of HMG CoA reductase by acting as a competitive inhibitor since it resembles the HMG CoA molecule.

To simplify the model, we can simply say that both the concentration of cholesterol and statins will negatively affect the rate of the enzyme HMG CoA reductase and will limit the production of mevalonate (another key intermediary in cholesterol biosynthesis).

Adjusting the Model Developed by Barahona, et al.

If we have a look again at the equation which governed the intracellular concentration of cholesterol, we see that we must target this equation if we want to take into account the de novo pathway for cholesterol biosynthesis.

[math]\displaystyle{ \frac{d[IC]}{dt} = (\chi_{I}d_{I}[IDL] + \chi_{L}d_{L}[LDL])\phi_{LR} + \chi_{L}d[LDL] - d_{IC}[IC] }[/math]

From the biology, we know that the 'de novo' pathway will be a source of intracellular cholesterol that is dependent upon the concentration of cholesterol there is in the cell. Thus, we obtain the straightforward addition below, where n denotes the rate constant for the biosynthesis of cholesterol determined by the maximum rate of HMG-CoA reductase, the limiting step in the biosynthesis pathway.

[math]\displaystyle{ \frac{d[IC]}{dt} = n[IC] + (\chi_{I}d_{I}[IDL] + \chi_{L}d_{L}[LDL])\phi_{LR} + \chi_{L}d[LDL] - d_{IC}[IC] }[/math]