- Research article
- Open Access
An in silico model of the ubiquitin-proteasome system that incorporates normal homeostasis and age-related decline
© Proctor et al; licensee BioMed Central Ltd. 2007
Received: 03 November 2006
Accepted: 21 March 2007
Published: 21 March 2007
The ubiquitin-proteasome system is responsible for homeostatic degradation of intact protein substrates as well as the elimination of damaged or misfolded proteins that might otherwise aggregate. During ageing there is a decline in proteasome activity and an increase in aggregated proteins. Many neurodegenerative diseases are characterised by the presence of distinctive ubiquitin-positive inclusion bodies in affected regions of the brain. These inclusions consist of insoluble, unfolded, ubiquitinated polypeptides that fail to be targeted and degraded by the proteasome. We are using a systems biology approach to try and determine the primary event in the decline in proteolytic capacity with age and whether there is in fact a vicious cycle of inhibition, with accumulating aggregates further inhibiting proteolysis, prompting accumulation of aggregates and so on. A stochastic model of the ubiquitin-proteasome system has been developed using the Systems Biology Mark-up Language (SBML). Simulations are carried out on the BASIS (Biology of Ageing e-Science Integration and Simulation) system and the model output is compared to experimental data wherein levels of ubiquitin and ubiquitinated substrates are monitored in cultured cells under various conditions. The model can be used to predict the effects of different experimental procedures such as inhibition of the proteasome or shutting down the enzyme cascade responsible for ubiquitin conjugation.
The model output shows good agreement with experimental data under a number of different conditions. However, our model predicts that monomeric ubiquitin pools are always depleted under conditions of proteasome inhibition, whereas experimental data show that monomeric pools were depleted in IMR-90 cells but not in ts20 cells, suggesting that cell lines vary in their ability to replenish ubiquitin pools and there is the need to incorporate ubiquitin turnover into the model. Sensitivity analysis of the model revealed which parameters have an important effect on protein turnover and aggregation kinetics.
We have developed a model of the ubiquitin-proteasome system using an iterative approach of model building and validation against experimental data. Using SBML to encode the model ensures that it can be easily modified and extended as more data become available. Important aspects to be included in subsequent models are details of ubiquitin turnover, models of autophagy, the inclusion of a pool of short-lived proteins and further details of the aggregation process.
Ageing is due to the gradual accumulation of unrepaired random molecular faults, which leads to an increased fraction of damaged cells and eventually to the functional impairment of older tissues and organs [reviewed in ]. Cellular components are continuously renewed by the processes of catabolism and resynthesis as they wear out or become damaged. However the mechanisms of renewal are not perfectly efficient and over time there may be an increase in damaged molecular structures. For example, there is an accumulation of incompletely degraded intralysosomal waste material, known as lipofuscin, in postmitotic cells . In 2001, Terman proposed the "garbage catastrophe theory" of ageing  which states that the process of ageing may derive from imperfect clearance of oxidatively damaged, relatively indigestible material, the accumulation of which further hinders cellular catabolic and anabolic functions. Terman points out that this theory applies mainly to postmitotic cells and not to proliferating cells. This is supported by the evidence that tissues composed of postmitotic cells such as the heart and brain show the most pronounced age-related changes [4, 5], whereas tissues with constantly proliferating cell populations such as the intestinal epithelium show only minor changes with age . A possible reason for this difference is that during cell growth and division, most cellular structures are renewed and there is a dilution of any undigested material.
Terman mainly applied his theory to the lysosomal pathway of degradation. He suggested that lysosomes gradually accumulate lipofuscin and that the lipofuscin granules attract lysosomal enzymes which, however, fail to degrade the pigment. This leads to a shortage in lysosomal enzymes, resulting in a decrease in degradation by this pathway. However, the major pathway for the rapid degradation of proteins is the ubiquitin-proteasome pathway. Proteins need to be degraded either because they are short-lived e.g. key regulatory proteins or because they have become damaged and partially unfolded. Partially unfolded proteins have exposed hydrophobic surfaces and are liable to form aggregates. Molecular chaperones bind to the exposed surfaces and an attempt is made to refold the protein. If refolding is unsuccessful then chaperones assist in the removal of these proteins and so prevent their aggregation .
Interestingly, proteasomes themselves are degraded by the lysosomal system and there is evidence for cross-talk between the two degradation pathways [8, 9]. It has been observed that there is a progressive decline in the overall proteolytic capacity of the cell with age [10, 11] resulting in an accumulation of oxidized and cross-linked proteins [12, 13]. These observations have been linked to an increase in free radicals with age . However, there is also evidence that free radicals do not increase with age [14, 15], suggesting that the observed increase in damaged proteins may be a result of a decline in removal of damage rather than an increase in the rate of damage itself. Proteins that have been damaged by free radicals appear resistant to proteolytic degradation and can act as inhibitors of lysosomal hydrolases and proteasome activity . The small aggregates which precede the formation of mature amyloid fibres are particularly pathogenic in many amyloid diseases . Sitte et al found a dramatic decline of proteasome activity, but not of proteasome enzyme levels in nondividing human BJ fibroblasts during hyperoxic ageing. Such a decline is not a cell culture artefact because it has been observed in many tissues including cardiac and skeletal muscle, skin, lens and brain (reviewed in ). Sitte et al did not observe a decrease in lysosomal activity with time in nondividing BJ fibroblasts. However, both proteasome and lysosomal activity decreased in proliferating BJ fibroblasts .
It has been reported that aggregated protein directly impairs the function of the ubiquitin-proteasome system [19, 20]. Another set of data indicates that aggregation occurs when the capacity of the proteasome degradation pathway is exceeded but no distinction could be drawn between the possibility of increased substrate expression or a decline in proteasome activity . Therefore, it is not clear whether aggregated protein directly impairs proteolysis or whether its presence is a symptom of impaired proteolysis. The issues of cause versus effect require further investigation particularly in the role of ageing.
Many neurodegenerative diseases are characterised by the presence of distinctive ubiquitin-positive inclusion bodies in affected regions of the brain. These inclusions consist of insoluble, unfolded, ubiquitinated polypeptides that fail to be targeted and degraded by the proteasome . It has been suggested that cells possess a protective mechanism to sequester and isolate toxic misfolded cytoplasmic intermediates by delivering them to ubiquitin-rich structures, called aggresomes at the microtubule-organising centre (MTOC) . Nuclear inclusions may serve a similar function, sequestering aggregation-prone polypeptides into relatively inert higher-order complexes. It has been demonstrated that misfolded proteins are less toxic to cells capable of forming nuclear inclusions .
In order to be recognised and degraded by the 26S proteasome, proteins must first be tagged with chains of four or more ubiquitin molecules. Ubiquitin is a small protein found in all eukaryotic organisms. Ubiquitination occurs as a result of the sequential activity of three classes of enzymes, E1 (ubiquitin activating enzyme), E2 (ubiquitin conjugating enzyme), and E3 (ubiquitin protein ligase). In the first step, ubiquitin is activated through ATP hydrolysis mediated by an E1 enzyme. The activated ubiquitin is then transferred to an E2 (ubiquitin conjugate) enzyme. There are approximately 50 different E2 proteins in mammalian cells  that have limited substrate specificities. E2 bound to ubiquitin then forms a complex with an E3 (ubiquitin ligase) enzyme. There are hundreds of E3 proteins  specific for different substrates. Three classes of E3 proteins have been identified: the HECT (homologous to E6-AP C-terminus), the RING (really interesting new gene) finger, and the U-box domain types (see review by Robinson and Ardley ). The following describes the action of a HECT E3. As well as binding to the E2-ubiquitin complex, E3 binds to the protein which is to be degraded. The activated ubiquitin is then transferred to the substrate with the release of E2 and E3. At this stage the protein is mono-ubiquitinated. Further ubiquitin is then attached to form a chain by the further action of E1 and E2 enzymes. A ubiquitin chain of at least four ubiquitin molecules has physical affinity for the proteasome and delivers the substrate for degradation in an ATP-dependent manner . During degradation, the ubiquitin is released for recycling and the protein is cleaved into short peptides, which are then reduced to amino-acids by cytosolic proteases. Note that for simplicity our model ignores alternative mechanisms of chain assembly (conjugation of preassembled chains onto substrates etc., as discussed in , and does not require the activity of E4 enzymes, which for some substrates may promote chain elongation .
A large number of de-ubiquitinating enzymes (DUBs) are also found within all eukaryotic cells . In addition to processing the primary pro-protein translation products of the ubiquitin genes DUBs may edit polyubiquitin chains on substrates and thereby prolong the half-life of such substrates by preventing the formation of threshold length chains required for proteasome binding.
A model of protein turnover and the role of the molecular chaperone Hsp90 in maintaining protein homeostasis has previously been developed . This model contains a pool of proteins in its correctly folded native state (represented by a species NatP). At any point in time, a native protein can become misfolded with the rate of this reaction depending on the level of reactive oxygen species (ROS) within the cell (higher levels of ROS lead to an increase in the rate of misfolding). The pool of misfolded proteins is represented by a species MisP. There are three possible outcomes for a misfolded protein. The first possibility is that a misfolded protein binds to Hsp90 and once bound, it can either be refolded and return to the pool of native protein or it may dissociate from Hsp90 and so remain in its misfolded state. Secondly, a misfolded protein may be degraded and so removed. Thirdly, a misfolded protein may bind to another misfolded protein to form a small aggregate (or bind to a previously formed aggregate). The model contains details of the regulation of Hsp90 through its interaction with Heat Shock Factor-1 (HSF1). However, no detail is included of either the degradation pathway or the process of aggregation.
List of species with database terms and initial conditions
Ubiquitin activating enzyme
Ubiquitin conjugating enzyme
Ubiquitin ligase (HECT)
Reactive oxygen species
Misfolded protein bound by E3
E1 bound by Ub
E2 bound by Ub
Monoubiquitinated misfolded protein
Polyubiquitinated misfolded protein
MisP_Ub2, ... MisP_Ub8
Polyubiquitinated protein bound to proteasome
MisP_Ub4_Proteasome ... MisP_Ub8_Proteasome
Sequestered aggregated protein
Aggregated protein bound to proteasome
Dummy species to track chain length of substrate at degradation
degUb4, ..., degUb8
List of reactions and parameter values for the ubiquitin-proteasome system
Kinetic rate law
0.01 mol s-1
2.0E-6 mol-1 s-1
1.7E-5 mol-1 s-1
E1, Ub, ATP
2.0E-4 mol-1 s-1
1.0E-3 mol-1 s-1
E2, E3, MisP_Ub
1.0E-3 mol-1 s-1
1.0E-2 mol-1 s-1
1.0E-2 mol-1 s-1
1.0E-2 mol-1 s-1
DUB, MisP_Ub7, Ub
1.0E-5 mol-1 s-1
DUB, MisP, Ub
1.0E-5 mol-1 s-1
1.0E-5 mol-1 s-1
1.0E-5 mol-1 s-1
DUB, MisP_Ub3, Ub, Proteasome
1.0E-5 mol-1 s-1
DUB, Ub, MisP_Ub4_Proteasome
1.0E-5 mol-1 s-1
DUB, Ub, MisP_Ub7_Proteasome
1.0E-5 mol-1 s-1
ATP, MisP_Ub4_ Proteasome
ADP, Proteasome, 4Ub,
ATP, MisP_Ub8_ Proteasome
ADP, Proteasome, 8Ub,
We carried out a sensitivity analysis to see how changing each of the model parameters affects the simulation output. The only parameters which have a significant effect on the results when increased or decreased by an order of magnitude were k1, k2, k67, k68, and k69. Varying the parameter for protein synthesis, k1, affects the levels of native protein but only an increase in k1 has any affect on any of the other species. Increasing protein synthesis leads to a depletion in monomeric ubiquitin pools due to the increase pool of protein requiring degradation. An increase in the parameter for misfolding, k2, leads to an increase in misfolded protein, a decrease in native protein, all the ubiquitin ending up in conjugates and an increase in degradation. Conversely, a decrease in k2, leads to higher levels of native protein, more monomeric ubiquitin, less ubiquitin conjugates and less degradation. Increasing the parameter for binding of polyubiquitinated misfolded protein to the proteasome, k67, has no effect on the model output, but a decrease in k67 leads to a depletion in monomeric ubiquitin pools, an increase in unbound polyubiquitinated misfolded protein and less degradation. Increasing the parameter for the chain shortening of polyubiquitinated misfolded protein, k68, has the same effect as decreasing k67. However, decreasing k68 has little effect on model output except that nearly all the degradation reactions are via the substrates with the longest ubiquitin chain. Increasing the parameter for proteasome activity, k69, leads to more degradation and less ubiquitin conjugates, whereas decreasing k69 has the opposite effect and also results in depletion of monomeric ubiquitin due to less recycling of ubiquitin.
We also varied each of the parameters in turn to see which parameters affect the kinetics of aggregation. We found that apart from the parameters actually involved in the aggregation steps (k71, k72, k73, k74), the only parameters which affect aggregation are k1, k2, k3 and k61, Either an increase in protein synthesis (k1) or in the misfolding rate (k2) leads to an increase in aggregated protein which in turn leads to both an increase in sequestered aggregated protein (SeqAggP) and binding of aggregated protein to the proteasome (AggP_Proteasome). This is due to the increase level of misfolded protein which can not be degraded. On the other hand an increase in refolding (k3) leads to a reduction in aggregation. Increased refolding could occur if there was an upregulation in molecular chaperones. Increasing the rate of binding of E3 to misfolded protein (k61) also reduces the level of aggregation suggesting that, if a particular protein was involved in aggregation, then overexpressing the relevant E3 for a misfolded protein could help to prevent the accumulation of aggregates. However, aggregates are usually composed of a variety of proteins, so this approach may not be feasible. As expected, increasing the rate of aggregation by either increasing k71 or k72 leads to more aggregation, and varying the parameters k73 or k74 changes the ratio of SeqAggP to AggP_Proteasome.
Parameters which affect protein turnover
Shutting down E1 activity
We have built a computer model of the ubiquitin-proteasome system using an iterative process of model-building and validation of the model against experimental data. The initial conditions and parameters for the model were chosen so that the model output fitted the experimental data under normal conditions. The model was then validated by mimicking the experimental conditions such as proteasome inhibition and shutting down E1 activity. For example, with our original choice of parameters our model predicted that if E1 activity is shut down, then there is an almost immediate conversion of ubiquitin conjugates to monomeric ubiquitin. However, experimental data show that this conversion takes about two hours (Fig 5a). Therefore, we adjusted the parameters to reduce the rates of the reactions involved in ubiquitination, de-ubiquitination and degradation until the model output fitted the experimental data.
With this final set of parameters, the model output also agrees with experimental data under normal conditions and when the proteasome is inhibited. However, under conditions of proteasome inhibition our model only agrees with data for IMR-90 cells and not ts20 cells with respect to the pool of monomeric ubiquitin. IMR-90 cells are diploid fibroblast-like cells isolated from human fetal lung ; they are nontransformed and with prolonged culture will reach the Hayflick limit and undergo senescence. The ts20 cell line is an immortalized line derived from the E36 Chinese hamster ovary line . Whereas cell lines were found to be variable with respect to ubiquitin content (Fig 3b) there did not seem to be a correlation between ubiquitin content and immortalization or oncogenic transformation (comparing the normal diploid IMR-90 to the human glioblastoma line U87MG, for example). It may be that the robust growth of transformed cells places demands on ubiquitin pools and that such cells have devised mechanisms to maintain or replenish their pools; if so the IMR-90 results may be a more accurate representation of normal ubiquitin homeostasis. In order to be able to use our model for different cell types (including transformed cells) it will be necessary to add detail of ubiquitin turnover instead of assuming a constant pool.
We carried out a sensitivity analysis to see which parameters affect the model output, not only under normal conditions but also under experimental conditions such as proteasome inhibition and shutting down E1 activity. Under normal conditions, pools of native protein were affected by changing the protein synthesis or misfolding rate, and pools of monomeric ubiquitin were affected by changing the rate of polyubiquitinated misfolded binding to the proteasome, the rate of chain shortening of polyubiquitinated misfolded protein bound to proteasomes, or the rate of proteasome activity. Under conditions of proteasome inhibition, increasing the rate of either protein synthesis or protein misfolding, or decreasing the rate of either refolding or binding of misfolded proteins to E3 led to an increase in protein aggregation (and vice versa). These results suggest that the availability of chaperones and E3 enzymes are important in maintaining protein homeostasis.
We also examined which parameters affect protein half-life. The most important parameter in this respect was the parameter for misfolding with a ten-fold increase/decrease leading to a ten-fold decrease/increase in protein half-life. Since the misfolding reaction can be driven by the level of ROS in the cell, our model would predict a similar outcome by changing ROS levels, with an increase in ROS leading to higher misfolding and an increase in degradation, consistent with current thought linking protein oxidation and ageing . Our model predicts that a ten-fold decrease in the rate of proteasome activity leads to a four-fold increase in protein half-life but that an increase in proteasome activity did not decrease protein half-life. It has been suggested that since proteasome levels can only increase slowly due to their slow assembly time, that they must have residual capacity for degrading proteins when cells are stressed . Our model allows for this residual capacity and explains why a ten-fold increase in proteasome activity has no effect on protein half-life.
There are also many other modifications that we could make to the model. For example, we have assumed that E3 binds to misfolded protein in one simple step; however it is likely that when a native protein misfolds, the first step would be for a molecular chaperone (e.g. Hsp90) to bind to the exposed hydrophobic surface to prevent inappropriate interactions. We could easily add this detail by including the following steps:
MisP + Hsp90 ↔ Hsp90_MisP (reversible reaction)
Hsp90_MisP + E3 → E3_MisP + Hsp90
Similarly, we could include the possibility that chaperones aid in delivery of ubiquitinated proteins to the proteasome.
Our model includes a HECT E3, but it would be easy to adapt the model for other types of E3 enzymes. For example if E3 is a RING ligase, then it has to be in complex with E2 before it can bind to the misfolded protein. This could be represented by the following reactions:
E3 + E2_Ub ↔ E2_E3_Ub (reversible reaction)
MisP + E2_E3_Ub → E2_E3_MisP_Ub
E2_E3_MisP_Ub → MisP_Ub + E2 + E3
Our model only includes one pool of long-lived proteins and so far ignores the fact that the majority of proteasome activity is involved in the turnover of short-lived regulatory proteins and also the elimination of incorrect newly-synthesised proteins. It would make the model more realistic to include an additional pool of short-lived proteins and this work is currently in progress.
We have so far restricted our models of protein turnover to the ubiquitin-proteasome system; however work is currently in progress to develop other models of protein turnover. It will be important to incorporate models of lysosomal pathways since studies suggest that protein aggregates such as mutant huntingtin are removed by autophagy [reviewed in ].
We have included the possibility of aggregated protein inhibiting proteasomes. Another possible cause of proteasome inhibition is via unaggregated damaged protein which could also sequester proteasomes. It would be fairly straight forward to add this extra detail to our current model. By running simulations over long time periods, the models can then be used to see how proteasome activity is affected by damaged and/or aggregated protein. This would be more realistic, than simply setting the parameter for proteasome activity to zero.
Our model predicted that levels of native protein did not change when misfolded protein levels rise, but in reality we would not expect total protein levels to continually rise. This prediction was due to our assumption of a constant protein synthesis rate. We could later include detail of feedback mechanisms such as the unfolded protein response which downregulates translation when unfolded protein accumulates in the endoplasmic reticulum.
We have encoded our model using SBML which not only allows for easy portability but also enables the model to be modified and extended in a very straight-forward way. The model is available on the BASIS website and can be freely accessed from the public repository, so that the interested reader can run simulations, make their own modifications, or download the SBML code. The BASIS website has a stochastic simulator and a database for storing models and simulation results.
Computer models of other molecular mechanisms involved in ageing are being developed at Newcastle University, for example the role of mitochondria, the antioxidant system, and DNA damage signalling pathways. Methods to link models in an automated way are also under development which will greatly facilitate the development of integrated models. For example, our model of the ubiquitin-proteasome system could be linked to a model of the mitochondria. This would be very valuable since neurodegeneration is not only associated with an increase in aggregated protein but also an increase in damaged mitochondrial DNA .
We have developed a model of the ubiquitin-proteasome system using an iterative approach of model building and validation against experimental data. We have used SBML to encode the model to ensure that it can be easily modified and extended as more data become available. Important aspects to be included in subsequent models are details of ubiquitin turnover, models of autophagy, the inclusion of a pool of short-lived proteins and further details of the aggregation process. The model and its extensions will be an invaluable aid to further our understanding of the cellular mechanisms involved in maintaining protein homeostasis and how it is disturbed during ageing.
U87MG human glioblastoma cells, NIH-3T3 mouse fibroblasts and IMR90 human fibroblasts were maintained in 5% CO2 at 37°C in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% of a 2:1 mixture of fetal bovine and donor bovine serum. ts20 cells harbouring a temperature sensitive E1 ubiquitin activating enzyme (a gift from Dr. Sergio Grinstein, University of Toronto) were maintained in 5% CO2 at 34°C in DMEM as above. Inactivation of E1 was achieved by incubating the cells at 42°C for the indicated time. When proteasome inhibitor (Calbiochem, San Diego, CA, USA) was used, it was added to the culture medium 3 or 4 hours prior to harvesting the cells for western analysis at a final concentration of 50 μM.
Preparation of cell extracts for western blot analysis
Cells were harvested in protein lysis buffer (20 mM Tris-HCl pH 7.5, 150 mM NaCl, 1% Nonidet P-40 (Sigma, St. Louis, Missouri, USA), 0.5 mM EDTA, 5% glycerol) containing the following protease and phosphates inhibitors: phenylmethylsulfonyl fluoride, leupeptin (Sigma, St. Louis, Missouri, USA), aprotinin (Sigma, St. Louis, Missouri, USA), sodium fluoride (NaF), sodium pyrophosphate (NaPPi) and N-ethylmaleimide. The cell extracts were sonicated and centrifuged at 14,000 rpm for 20 minutes at 4°C to pellet cellular debris. The soluble fractions were recovered and the protein concentration was determined using the Bradford protein assay (Bio-Rad Laboratories Inc., Mississauga, Ontario, Canada). 40 μg of cytoplasmic extracts were then resolved on a two-phase SDS-polyacrylamide gel (15 and 8%) and electroblotted onto a hybond C nitrocellulose membrane (Amersham Pharmacia Biotech, Baie D'Urfé, Québec, Canada). The membranes were autoclaved on a liquid cycle for 45 minutes to enhance the detection of poly-ubiquitinated proteins, stained with Ponceau S (Sigma, St. Louis, Missouri, USA) and analyzed by western blotting with the indicated antibody. The proteins were visualized by a horseradish peroxidase method using the ECL kit from Kirkegaard and Perry Laboratories Inc., Gaithersburg, Maryland, USA.
The rabbit polyclonal antibody used to detect ubiquitin was purchased from Dako. The mouse monoclonal antibodies recognizing actin and GAPDH were from Sigma (St. Louis, Missouri, USA) and Stressgen Bioreagents (Victoria, British Columbia) respectively.
The model was built by first drawing a diagram of the system (see Figure 1) and listing all the elements of the model and the interactions between them. We then specified the initial amount of each element, the rate laws and the values of the parameters. With this information, the model can be translated into computer code and simulated. To encode the model we used SBML shorthand  which can then be converted into the Systems Biology Markup Language (SBML), a computer-readable format for representing biochemical models . The model was imported into the BASIS (Biology of Ageing e-Science Integration and Simulation) system [39, 40] and simulations were carried out using a stochastic simulator based on the Gillespie algorithm . The SBML code is provided as supplementary material (see Additional file 1). The model is also in the public repository of BASIS (urn:basis.ncl:model:2215) and will soon be available from the Biomodels database [42, 43].
Elements of the model
It is first necessary to specify all the elements of the model and the interactions between the elements. Using the terminology of SBML, we refer to elements as "species" and interactions as "reactions". The species needed for the model are a pool of native proteins, misfolded proteins, ubiquitin, ubiquitin activating enzyme (E1), ubiquitin conjugating enzyme (E2), ubiquitin ligase (E3), de-ubiquitinating enzymes (DUB), proteasomes, ATP, ADP, AMP, reactive oxygen species (ROS), and complexes to represent binding of elements e.g. ubiquitinated protein. A full list of the species with their names, database terms and their initial amounts is shown in Table 1.
Interactions between elements of the model
In our previous model of the chaperone system , we included a single step for the degradation of misfolded proteins:
where MisP represents the pool of misfolded protein.
We now replace this single step with a more detailed representation of the ubiquitin-proteasome system.
Binding of E3 to misfolded protein
The first step in the degradation pathway is that the substrate (in this case MisP) is bound to its specific ubiquitin ligase (E3). The binding of E3 to MisP produces a complex E3_MisP. This reaction is reversible and so we have two reactions as follows:
We use mass action stochastic kinetics for the rate laws . The binding reaction is a second-order reaction since there are two reactants and is given by k61<#E3><#MisP>, where # denotes the number of molecules. The disassociation reaction is a first-order reaction and is given by k61r<#E3_MisP>. The values for k61 and k61rcan be estimated from knowledge of the protein half-life and the steady state level of protein in its misfolded state.
Action of ubiquitin enzymes, E1 and E2 leading to first ubiquitination event
Before a substrate can be ubiquitinated, ubiquitin must be activated by the ubiquitin-activating enzyme (E1). This reaction requires one molecule of ATP. The binding of E1 to ubiquitin (Ub) produces a complex which we denote by E1_Ub:
The activated ubiquitin is then transferred to the ubiquitin-conjugating enzyme (E2) which releases E1 making it available for further ubiquitin-activating reactions:
Although it is possible that this reaction is reversible, it seems unlikely that E2 would return the Ub to E1. So we assume that the reverse reaction is negligible and do not include it in the model. Ubiquitin bound to E2 is then transferred to the misfolded protein bound to E3, releasing E2 and E3:
Extending the ubiquitin chain
Further ubiquitin is then attached to the monoubiquitinated protein via the action of E2. Since each ubiquitin molecule must be activated by E1 before being transferred to E2, each step in the formation of polyubiquitin chain uses one molecule of ATP. The extension of the ubiquitin chain is modelled by the following set of reactions:
Each of the above reactions is reversible, with the reversible reactions requiring the activity of a deubiquitinating enzyme (DUB). We assume that the shortening of ubiquitin chains by DUBs is a step-wise process, although we could later modify the model to allow whole chains to be removed if it becomes of interest to model this. The step-wise removal of ubiquitin is modelled by the following reactions:
Binding to the proteasome and degradation
We assume that ubiquitinated protein first binds to the proteasome and waits for degradation and that it can remain bound so long as at least 4 ubiquitin sub-units are attached. This can be modelled by the following reactions:
where MisP_Ub4_Proteasome etc. represents an ubiquitinated protein bound to the proteasome.
Bound ubiquitinated protein can also be de-ubiquitinated in a step-wise process, so that proteins with longer chains stay longer at the proteasome. This is modelled by the following set of reactions:
When a bound protein has a chain shorter than four, then it dissociates from the proteasome. This is shown in the following reaction:
Any protein which has a chain of four or more ubiquitin molecules may be degraded by the proteasome in an ATP dependent reaction:
Protein synthesis and refolding
This model can either be combined with the chaperone model  or can be explored as a separate module. In order to investigate how the various parameters and species amounts affect the model predictions for the rate of protein turnover, we first carried out simulations on the separate module. In order to do this we added the reactions of protein synthesis, protein misfolding and protein refolding to the model:
The rate of misfolding depends on the level of reactive oxygen species (ROS) within the cell. We have omitted details of the chaperone activity, where a misfolded protein would bind to Hsp90 before refolding can take place. However, when we combine the models, this simple reaction will be replaced with the more detailed reactions given in Proctor et al., . In this model, we assume that ROS takes a constant value throughout the simulation but we can look at the effect of varying the amount of ROS. It would also be simple to add reactions for the generation and removal of ROS which would then make it possible to examine the effects of increasing ROS over time.
If a misfolded protein is not removed immediately by refolding or degradation, then there is a chance that it will interact with another misfolded protein to form a small aggregate, or it may interact with an existing aggregate to form a larger aggregate. (Here we will not be concerned with the size of aggregates – a more detailed model is in preparation). An aggregate may be sequestered to prevent it interfering with the cellular machinery or it may bind to the proteasome. We include enough detail to be able to examine whether an increase in misfolding (for example by an increase in levels of ROS) leads to an increase in aggregation and inhibition of the proteasome which in turn leads to an even greater level of aggregated protein.
The following set of reactions show how protein aggregation can be simply modelled:
where AggP represents the pool of aggregates. The aggregates may be sequestered (SeqAggP) which keeps them out of harms way, or they may bind to the proteasome (AggP_Proteasome) and so inhibit its function. Ubiquitinated misfolded protein may also form aggregates but this will not normally occur unless the proteasome is inhibited. As there are many ways in which misfolded proteins with ubiquitin chains of different lengths can interact, we only list a subset of the reactions:
Setting the initial amounts of the species and parameter values
Before the model can be simulated, it is necessary to specify the initial amounts of each species and the parameter values. There is experimental data on levels of ubiquitin and proteasomes in human cells. For example, Haas et al.  measured levels of monomeric ubiquitin and conjugated ubiquitin in IMR-90 cells giving values of 71.5 pmol/106 cells and 61.1 pmol/106 cells respectively. From this we can calculate that there are approximately 107 ubiquitin molecules per cell (estimated at between 1 and 10% of total protein molecules). Since we are only modelling part of the total cellular protein, we have scaled down the level of ubiquitin to 500 molecules per cell. The total level of proteasomes per cell has been estimated to be 8 × 105 per cell in L929 cells  which gives a ratio of about one proteasome to every 10 ubiquitin molecules. Since ubiquitin has many other functions apart from the ubiquitin-proteasome pathway, the ratio of proteasomes to ubiquitin must be less than 1:10 and so we assumed a ratio of 1:5.
The level of E1 has been measured in IMR90 cells by western analysis of cell lysate compared to recombinant protein (Tsirigotis and Gray, unpublished data) and is estimated to be 1 million molecules per cell. This is one order less than the level of ubiquitin and so we chose an initial amount of 100 for E1. We do not have data for E2 or E3 enzymes but for the sake of building the model assume that for specific substrates these would be present at roughly the same abundance as E1. We do not have any experimental data for the abundance of DUBs within cells but assume that these are less abundant than ubiquitin and set the initial amount to 200.
The level of native protein, NatP, was set at 500 and at this level the ratio of native protein to ubiquitin and proteasomes enabled efficient degradation to take place. If we increased the level of NatP to 5000, then it was also necessary to increase the pool of ubiquitin and proteasomes to prevent depletion of monomeric ubiquitin under normal conditions. The levels of ATP, ADP and AMP were set at a fairly arbitrary level of 10000, 1000 and 1000 respectively and the levels were kept constant by imposing a boundary condition to these species. We chose to do this instead of allowing the levels to fluctuate as cells possess mechanisms to maintain ATP levels at fairly constant levels and it is unlikely that the ubiquitin-proteasome system would cause major changes to levels of these molecules. ROS levels were set at an arbitrary value of 10 and the level remained constant as there are no mechanisms in this model that would cause ROS levels to change. The initial amounts of all the species are shown in Table 1.
List of reactions and parameter values for the aggregation process
Kinetic rate law
1.0E-8 mol-1 s-1
1.0E-8 mol-1 s-1
1.0E-8 mol-1 s-1
1.0E-8 mol-1 s-1
1.0E-8 mol-1 s-1
1.0E-8 mol-1 s-1
We now describe each of the parameters in turn and give an interpretation of their values in terms of reaction kinetics. The parameters k1, k2 and k3 were set so that the half-life of the pool of protein equals about 10 hours, and under normal conditions, there is only a very low proportion of misfolded protein . The parameter k61 for E3/MisP binding is set so that the average time for one E3 to bind to MisP is about 10 minutes. The parameter k61ris set so that the reverse reaction (dissociation of E3 from MisP) happens once every 100 minutes. So the forward reaction is 10 times stronger than the reverse reaction. The parameter k62 for E1/Ub binding is set so that about 40 ubiquitin molecules are activated every minute and the parameter k63 for E2/Ub binding is set to give about 6 reactions per minute. The parameter k64 for monoubiquitination is set so that once MisP is bound by E3 it takes about 10 seconds to receive the first Ub molecule. The parameter k65 for polyubiquitination is set so that it took about one second for each additional Ub molecule to be added to the chain. So each polyubiquitination step is about 10 times faster than the first ubiquitination reaction. The parameter k66 for the de-ubiquitination reactions is set so that it takes about 8 minutes for each reaction. It is necessary that this reaction is very slow compared to the ubiquitination reactions so that the model output fits the experimental data when E1 activity is shut down (see section 4.4). The parameter k67 for the binding of ubiquitinated proteins (with chains of 4 or more ubiquitin molecules) to the proteasome is set so that under normal conditions when there are about 200 ubiquitin conjugates present in the cell, on average about 5 reactions per second are occurring. The parameter k68 for the de-ubiquitination of proteasome-bound ubiquitinated proteins is set to be equal to k66. We assume that this reaction is slow so that proteins with long chains are highly likely to reside long enough at the proteasome for degradation to take place, whereas proteins with short chains may be shortened to below the required threshold and would then dissociate from the proteasome. The parameter k69 for proteasome activity is set so that it takes about 15 seconds to degrade each proteasome-bound ubiquitinated protein. We initially set k71 and k72 to be equal to 10-7 as in the model of Proctor et al. . However, with these values, aggregates sometimes formed even under normal conditions (data not shown). Therefore we lowered the values until we obtained results in which the appearance of aggregates was very unlikely under normal conditions. This was satisfied with k71 and k72 set equal to 10-8, The parameter k73 for sequestering aggregates is set so that if there is a low level of aggregates in the cell (1–10), then on average there are 0.005 reactions per second, so each reaction takes about 3 minutes for each reaction.
The parameter k74 is set so that there is an equal probability of an aggregated protein being sequestered or binding to the proteasome when all the proteasomes are available for binding.
Adapting the model to mimic different experimental procedures
We carried out an in silico experiment of inhibiting the proteasome by setting the parameter k69 = 0. We also varied each of the other parameters in turn to see which parameters affect the aggregation kinetics.
Blocking protein synthesis
By setting the parameter for protein synthesis, k1, to zero we can use our model to check the half-life of the protein pool. We also varied each of the other parameters in turn to see which parameters affect the half-life. For each parameter we increased and decreased its value ten-fold, re-ran the simulations and plotted the model output.
Shutting down E1 activity
Event structure for shutting down E1 activity in the computer model
t >= 1800
k62 = 0
Statistical analysis of repeat runs
Statistical analysis of repeated simulations was carried using the R programming language.
The work was supported by a pilot project award from the Institute of Aging, Canadian Institutes of Health Research. CJP is funded by BBSRC and EPSRC. We thank the BASIS team for their technical help. BASIS was funded by BBSRC, MRC, DTI and Unilever plc. We also thank Darren Wilkinson for reading the manuscript and for his helpful suggestions.
- Kirkwood TBL: Evolution of ageing. Nature. 1977, 270: 301-304. 10.1038/270301a0PubMedView ArticleGoogle Scholar
- Brunk UT, Terman A: The mitochondrial-lysosomal axis theory of aging - Accumulation of damaged mitochondria as a result of imperfect autophagocytosis. Eur J Biochem. 2002, 269 (8): 1996-2002. 10.1046/j.1432-1033.2002.02869.xPubMedView ArticleGoogle Scholar
- Terman A: Garbage catastrophe theory of aging: imperfect removal of oxidative damage?. Redox Rep. 2001, 6 (1): 15-26. 10.1179/135100001101535996PubMedView ArticleGoogle Scholar
- Smith CD, Carney JM, Starke-Reed PE, Oliver CN, Stadtman ER, Floyd RA, Markesbery WR: Excess Brain Protein Oxidation and Enzyme Dysfunction in Normal Aging and in Alzheimer Disease. PNAS. 1991, 88 (23): 10540-10543. 10.1073/pnas.88.23.10540PubMed CentralPubMedView ArticleGoogle Scholar
- Terman A, Brunk UT: Lipofuscin: Mechanisms of formation and increase with age. Apmis. 1998, 106 (2): 265-276.PubMedView ArticleGoogle Scholar
- Martin K, Kirkwood TBL, Potten CS: Age changes in stem cells of murine small intestinal crypts. Exp Cell Res. 1998, 241 (2): 316-323. 10.1006/excr.1998.4001PubMedView ArticleGoogle Scholar
- Wickner S, Maurizi MR, Gottesman S: Posttranslational quality control: Folding, refolding, and degrading proteins. Science. 1999, 286 (5446): 1888-1893. 10.1126/science.286.5446.1888PubMedView ArticleGoogle Scholar
- Martinez-Vicente M, Sovak G, Cuervo AM: Protein degradation and aging. Experimental Gerontology. 2005, 40 (8-9): 622-633. 10.1016/j.exger.2005.07.005PubMedView ArticleGoogle Scholar
- Bjorkoy G, Lamark T, Brech A, Outzen H, Perander M, Overvatn A, Stenmark H, Johansen T: p62/SQSTM1 forms protein aggregates degraded by autophagy and has a protective effect on huntingtin-induced cell death. J Cell Biol. 2005, jcb.200507002-Google Scholar
- Friguet B, Bulteau AL, Chondrogianni N, Conconi M, : Protein degradation by the proteasome and its implications in aging. The Annals of the New York Academy of Sciences. 2000, 908: 143-154.View ArticleGoogle Scholar
- Szweda PA, Friguet B, Szweda LI: Proteolysis, free radicals, and aging. Free Radic Biol Med. 2002, 33 (1): 29-36. 10.1016/S0891-5849(02)00837-7PubMedView ArticleGoogle Scholar
- Sitte N, Merker K, Von Zglinicki T, Davies KJA, Grune T: Protein oxidation and degradation during cellular senescence of human BJ fibroblasts: part II - aging of nondividing cells. Faseb J. 2000, 14 (15): 2503-2510. 10.1096/fj.00-0210comPubMedView ArticleGoogle Scholar
- Sitte N, Merker K, Von Zglinicki T, Grune T, Davies KJA: Protein oxidation and degradation during cellular senescence of human BJ fibroblasts: part I - effects of proliferative senescence. Faseb J. 2000, 14 (15): 2495-2502. 10.1096/fj.00-0209comPubMedView ArticleGoogle Scholar
- Gredilla R, Sanz A, Lopez-Torres M, Barja G: Caloric restriction decreases mitochondrial free radical generation at complex I and lowers oxidative damage to mitochondrial DNA in the rat heart. Faseb J. 2001, 15 (9): 1589-1591.PubMedGoogle Scholar
- Miwa S, Riyahi K, Partridge L, Brand MD: Strategies for Engineered Negligible Senescence: Why Genuine Control of Aging May Be Foreseeable. 2004, 1019: 388-391. Lack of correlation between mitochondrial reactive oxygen species production and life span in Drosophila, Annals of the New York Academy of SciencesGoogle Scholar
- Sitte N, Huber M, Grune T, Ladhoff A, Doecke WD, Von Zglinicki T, Davies KJA: Proteasome inhibition by lipofuscin/ceroid during postmitotic aging of fibroblasts. Faseb J. 2000, 14 (11): 1490-1498. 10.1096/fj.14.11.1490PubMedView ArticleGoogle Scholar
- Chiti F, Calamai M, Taddei N, Stefani M, Ramponi G, Dobson CM: Studies of the aggregation of mutant proteins in vitro provide insights into the genetics of amyloid diseases. Proc Natl Acad Sci U S A. 2002, 99: 16419-16426. 10.1073/pnas.212527999PubMed CentralPubMedView ArticleGoogle Scholar
- Chondrogianni N, Gonos ES: Proteasome dysfunction in mammalian aging: Steps and factors involved. Experimental Gerontology. 2005, 40 (12): 931-938. 10.1016/j.exger.2005.09.004PubMedView ArticleGoogle Scholar
- Bence NF, Sampat RM, Kopito RR: Impairment of the ubiquitin-proteasome system by protein aggregation. Science. 2001, 292 (5521): 1552-1555. 10.1126/science.292.5521.1552PubMedView ArticleGoogle Scholar
- Bennett EJ, Bence NF, Jayakumar R, Kopito RR: Global impairment of the ubiquitin-proteasome system by nuclear or cytoplasmic protein aggregates precedes inclusion body formation. Mol Cell. 2005, 17 (3): 351-365. 10.1016/j.molcel.2004.12.021PubMedView ArticleGoogle Scholar
- Johnston JA, Ward CL, Kopito RR: Aggresomes: A cellular response to misfolded proteins. J Cell Biol. 1998, 143 (7): 1883-1898. 10.1083/jcb.143.7.1883PubMed CentralPubMedView ArticleGoogle Scholar
- Ardley HC, Hung CC, Robinson PA: The aggravating role of the ubiquitin-proteasome system in neurodegeneration. FEBS Lett. 2005, 579 (3): 571-576. 10.1016/j.febslet.2004.12.058PubMedView ArticleGoogle Scholar
- Arrasate M, Mitra S, Schweitzer ES, Segal MR, Finkbeiner S: Inclusion body formation reduces levels of mutant huntingtin and the risk of neuronal death. Nature. 2004, 431 (7010): 805-810. 10.1038/nature02998PubMedView ArticleGoogle Scholar
- Semple CAM, , : The Comparative Proteomics of Ubiquitination in Mouse. Genome Res. 2003, 13 (6b): 1389-1394. 10.1101/gr.980303PubMed CentralPubMedView ArticleGoogle Scholar
- Robinson PA, Ardley HC: Ubiquitin-protein ligases. Journal of Cell Science. 2004, 117: 5191-5194. 10.1242/jcs.01539PubMedView ArticleGoogle Scholar
- Pickart CM, Fushman D: Polyubiquitin chains: polymeric protein signals. Current Opinion in Chemical Biology. 2004, 8 (6): 610-616. 10.1016/j.cbpa.2004.09.009PubMedView ArticleGoogle Scholar
- Hochstrasser M: Lingering Mysteries of Ubiquitin-Chain Assembly. Cell. 2006, 124 (1): 27-34. 10.1016/j.cell.2005.12.025PubMedView ArticleGoogle Scholar
- Hoppe T: Multiubiquitylation by E4 enzymes: `one size' doesn't fit all. Trends in Biochemical Sciences. 2005, 30 (4): 183-187. 10.1016/j.tibs.2005.02.004PubMedView ArticleGoogle Scholar
- Wilkinson KD: Regulation of ubiquitin-dependent processes by deubiquitinating enzymes. Faseb J. 1997, 11 (14): 1245-1256.PubMedGoogle Scholar
- Proctor CJ, Söti C, Boys RJ, Gillespie CS, Shanley DP, Wilkinson DJ, Kirkwood TBL: Modelling the actions of chaperones and their role in ageing. Mechanisms of Ageing and Development. 2005, 126: 119-131. 10.1016/j.mad.2004.09.031.PubMedView ArticleGoogle Scholar
- Kulka RG, Raboy B, Schuster R, Parag HA, Diamond G, Ciechanover A, Marcus M: A Chinese hamster cell cycle mutant arrested at G2 phase has a temperature-sensitive ubiquitin-activating enzyme, E1. J Biol Chem. 1988, 263 (30): 15726-15731.PubMedGoogle Scholar
- Nichols WW, Murphy DG, Cristofalo VJ, Toji LH, Greene AE, Dwight SA: Characterization of a new human diploid cell strain, IMR-90. 1977, 196 (4285): 60-63.Google Scholar
- Sohal RS: Role of oxidative stress and protein oxidation in the aging process. Free Radic Biol Med. 2002, 33 (1): 37-44. 10.1016/S0891-5849(02)00856-0PubMedView ArticleGoogle Scholar
- Princiotta MF, Finzi D, Qian SB, Gibbs J, Schuchmann S, Buttgereit F, Bennink JR, Yewdell JW: Quantitating Protein Synthesis, Degradation, and Endogenous Antigen Processing. Immunity. 2003, 18 (3): 343-354. 10.1016/S1074-7613(03)00051-7PubMedView ArticleGoogle Scholar
- Ravikumar B, Rubinsztein DC: Role of autophagy in the clearance of mutant huntingtin: A step towards therapy?. Molecular Aspects of Medicine. 2006, 27 (5-6): 520-527. 10.1016/j.mam.2006.08.008PubMedView ArticleGoogle Scholar
- Bender A, Krishnan KJ, Morris CM, Taylor GA, Reeve AK, Perry RH, Jaros E, Hersheson JS, Betts J, Klopstock et : High levels of mitochondrial DNA deletions in substantia nigra neurons in aging and Parkinson disease. Nature Genetics. 2006, 38 (5): 515-517. 10.1038/ng1769PubMedView ArticleGoogle Scholar
- Wilkinson DJ: Stochastic Modelling for Systems Biology. 2006, Chapman & Hall/CRC PressGoogle Scholar
- Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, , Arkin AP, Bornstein BJ, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, Juty NS, Kasberger JL, Kremling A, Kummer U, Le Novere N, Loew LM, Lucio D, Mendes P, Minch E, Mjolsness ED, Nakayama Y, Nelson MR, Nielsen PF, Sakurada T, Schaff JC, Shapiro BE, Shimizu TS, Spence HD, Stelling J, Takahashi K, Tomita M, Wagner J, Wang J: The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics. 2003, 19 (4): 524-531. 10.1093/bioinformatics/btg015PubMedView ArticleGoogle Scholar
- Basis. http://www.basis.ncl.ac.uk
- Kirkwood TBL, Boys RJ, Gillespie CS, Proctor CJ, Shanley DP, Wilkinson DJ: Towards an e-biology of ageing: integrating theory and data. Nat Rev Mol Cell Biol. 2003, 4 (3): 243-249. 10.1038/nrm1051PubMedView ArticleGoogle Scholar
- Gillespie DT: Exact stochastic simulation of coupled chemical reactions. The Journal of Physical Chemistry. 1977, 31: 2340-2361. 10.1021/j100540a008.View ArticleGoogle Scholar
- Biomodels. , http://www.ebi.ac.uk/biomodels/
- Le Novere N, Bornstein B, Broicher A, Courtot M, Donizelli M, Dharuri H, Li L, Sauro H, Schilstra M, Shapiro B, Snoep JL, Hucka M: BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems. Nucleic Acids Research. 2006, 34 (suppl_1): D689-691. 10.1093/nar/gkj092.PubMed CentralPubMedView ArticleGoogle Scholar
- Haas AL, Bright PM: The dynamics of ubiquitin pools within cultured human lung fibroblasts. J Biol Chem. 1987, 262 (1): 345-351.PubMedGoogle Scholar
- Lodish H, Berk A, Zipursky SL, Matsudaira P, Baltimore D, Darnell J: Molecular Cell Biology. 2000, New York , W.H. Freeman and Company, 4thGoogle Scholar
- Interpro. http://www.ebi.ac.uk/interpro/
- Gene Ontology. http://www.geneontology.org
- Chemical Entities of Biological Interest database. http://www.ebi.ac.uk/chebi
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.