Modelling the actions of chaperones and their role in ageing

Systems Biology of Ageing

The ageing process is highly complex involving multiple processes operating at different biological levels. Systems Biology presents an approach using integrative computational and laboratory study that allows us to address such complexity. The approach relies on the computational analysis of knowledge and data to generate predictive models that may be validated with further laboratory experimentation. Our understanding of ageing is such that translational opportunities are within reach and systems biology offers a means to ensure that optimal decisions are made. We present an overview of the methods employed from bioinformatics and computational modelling and describe some of the insights into ageing that have been gained.

Protein intake and loss of proteostasis in the eldery

Ageing is a process of declining bodily function and a major risk factor of chronic diseases. The declining bodily function in ageing can cause loss of proteostasis (protein homeostasis), which is a balance between protein synthesis, folding, modification and degradation. For the elderly, adequate protein intake is necessary to prevent sarcopenia, frailty, fracture and osteoporosis as well as reduced resistance to infection. However, increasing the protein intake can enhance the risk of oxidized protein formation, loss of proteostasis and degenerative disorder occurrence. On the other hand, several studies show that protein restriction would increase longevity. The aim of this review was to explain the importance of determining the right amount and composition of protein intake for the elderly. Oxidative stress and molecular mechanism of proteostasis loss in ageing cells as well as its suppression pathway by protein restriction are discussed in this review. Keywords: ageing, dietary proteins, mTOR, oxidative stress, proteostasis loss

Genomic Variation, Evolvability, and the Paradox of Mental Illness

Twentieth-century genetics was hard put to explain the irregular behavior of neuropsychiatric disorders. Autism and schizophrenia defy a principle of natural selection; they are highly heritable but associated with low reproductive success. Nevertheless, they persist. The genetic origins of such conditions are confounded by the problem of variable expression, that is, when a given genetic aberration can lead to any one of several distinct disorders. Also, autism and schizophrenia occur on a spectrum of severity, from mild and subclinical cases to the overt and disabling. Such irregularities reflect the problem of missing heritability; although hundreds of genes may be associated with autism or schizophrenia, together they account for only a small proportion of cases. Techniques for higher resolution, genomewide analysis have begun to illuminate the irregular and unpredictable behavior of the human genome. Thus, the origins of neuropsychiatric disorders in particular and complex disease in general have been illuminated. The human genome is characterized by a high degree of structural and behavioral variability: DNA content variation, epistasis, stochasticity in gene expression, and epigenetic changes. These elements have grown more complex as evolution scaled the phylogenetic tree. They are especially pertinent to brain development and function. Genomic variability is a window on the origins of complex disease, neuropsychiatric disorders, and neurodevelopmental disorders in particular. Genomic variability, as it happens, is also the fuel of evolvability. The genomic events that presided over the evolution of the primate and hominid lineages are over-represented in patients with autism and schizophrenia, as well as intellectual disability and epilepsy. That the special qualities of the human genome that drove evolution might, in some way, contribute to neuropsychiatric disorders is a matter of no little interest.

Low-dimensional projection of stochastic cell-signalling dynamics via a variational approach

Noise and fluctuations play vital roles in signal transduction in cells. Various numerical techniques for its simulation have been proposed, most of which are not efficient in cellular networks with a wide spectrum of timescales. In this paper, based on a recently developed variational technique, low-dimensional structures embedded in complex stochastic reaction dynamics are unfolded which sheds light on new design principles of efficient simulation algorithm for treating noise in the mesoscopic world. This idea is effectively demonstrated in several popular regulation models with an empirical selection of test functions according to their reaction geometry, which not only captures complex distribution profiles of different molecular species but also considerably speeds up the computation.

Proteostasis collapse is a driver of cell aging and death

Cells and organisms grow old and die. We develop a biophysical model of the mechanism. Young cells are kept healthy by the positive processes of protein synthesis, degradation, and chaperoning (the activity of keeping proteins properly folded). But, with age, negative processes increase: Oxidative damage accumulates randomly in the cell’s proteins, healthy synthesis and degradation slow down, and—like overfilled garbage cans—chaperone capacity is exceeded. The chaperones are distracted trying to fold irreversibly damaged proteins, leading to accumulating misfolded and aggregated proteins in the cell. The tipping point to death happens when the negative overwhelms the positive. The model makes several quantitative predictions of the life span of the worm Caenorhabditis elegans.

What molecular processes drive cell aging and death? Here, we model how proteostasis—i.e., the folding, chaperoning, and maintenance of protein function—collapses with age from slowed translation and cumulative oxidative damage. Irreparably damaged proteins accumulate with age, increasingly distracting the chaperones from folding the healthy proteins the cell needs. The tipping point to death occurs when replenishing good proteins no longer keeps up with depletion from misfolding, aggregation, and damage. The model agrees with experiments in the worm Caenorhabditis elegans that show the following: Life span shortens nonlinearly with increased temperature or added oxidant concentration, and life span increases in mutants having more chaperones or proteasomes. It predicts observed increases in cellular oxidative damage with age and provides a mechanism for the Gompertz-like rise in mortality observed in humans and other organisms. Overall, the model shows how the instability of proteins sets the rate at which damage accumulates with age and upends a cell’s normal proteostasis balance.

Modelling the role of redox-related mechanisms in musculoskeletal ageing

The decline in the musculoskeletal system with age is driven at the cellular level by random molecular damage. Cells possess mechanisms to repair or remove damage and many of the pathways involved in this response are regulated by redox signals. However, with ageing there is an increase in oxidative stress which can lead to chronic inflammation and disruption of redox signalling pathways. The complexity of the processes involved has led to the use of computational modelling to help increase our understanding of the system, test hypotheses and make testable predictions. This paper will give a brief background of the biological systems that have been modelled, an introduction to computational modelling, a review of models that involve redox-related mechanisms that are applicable to musculoskeletal ageing, and finally a discussion of the future potential for modelling in this field.

Systems Biology in Aging Research

Systems biology is an approach to collect high-dimensional data and analyze in an integrated manner. As aging is a complicated physiological functional decline in biological system, the methods in systems biology could be utilized in aging studies. Here we reviewed recent advances in systems biology in aging research and divide them into two major parts. One is the data resource, which includes omics data from DNA, RNA, proteins, epigenetic changes, metabolisms, and recently single-cell-level variations. The other is the data analysis methods consisting of network and modeling approaches. With all the data and the tools to analyze them, we could further promote our understanding of the systematic aging.

pSSAlib: The partial-propensity stochastic chemical network simulator

Chemical reaction networks are ubiquitous in biology, and their dynamics is fundamentally stochastic. Here, we present the software library pSSAlib, which provides a complete and concise implementation of the most efficient partial-propensity methods for simulating exact stochastic chemical kinetics. pSSAlib can import models encoded in Systems Biology Markup Language, supports time delays in chemical reactions, and stochastic spatiotemporal reaction-diffusion systems. It also provides tools for statistical analysis of simulation results and supports multiple output formats. It has previously been used for studies of biochemical reaction pathways and to benchmark other stochastic simulation methods. Here, we describe pSSAlib in detail and apply it to a new model of the endocytic pathway in eukaryotic cells, leading to the discovery of a stochastic counterpart of the cut-out switch motif underlying early-to-late endosome conversion. pSSAlib is provided as a stand-alone command-line tool and as a developer API. We also provide a plug-in for the SBMLToolbox. The open-source code and pre-packaged installers are freely available from http://mosaic.mpi-cbg.de.

Modelling the molecular mechanisms of aging

The aging process is driven at the cellular level by random molecular damage that slowly accumulates with age. Although cells possess mechanisms to repair or remove damage, they are not 100% efficient and their efficiency declines with age. There are many molecular mechanisms involved and exogenous factors such as stress also contribute to the aging process. The complexity of the aging process has stimulated the use of computational modelling in order to increase our understanding of the system, test hypotheses and make testable predictions. As many different mechanisms are involved, a wide range of models have been developed. This paper gives an overview of the types of models that have been developed, the range of tools used, modelling standards and discusses many specific examples of models that have been grouped according to the main mechanisms that they address. We conclude by discussing the opportunities and challenges for future modelling in this field.

The Impact of Mathematical Modeling in Understanding the Mechanisms Underlying Neurodegeneration: Evolving Dimensions and Future Directions

Neurodegenerative diseases are a heterogeneous group of disorders that are characterized by the progressive dysfunction and loss of neurons. Here, we distil and discuss the current state of modeling in the area of neurodegeneration, and objectively compare the gaps between existing clinical knowledge and the mechanistic understanding of the major pathological processes implicated in neurodegenerative disorders. We also discuss new directions in the field of neurodegeneration that hold potential for furthering therapeutic interventions and strategies.

Formal Models of Biological Systems

Recent biomedical research studies are focused in the mechanisms by which misfolded proteins lead to the generation of oxidative stress in the form of reactive oxygen species (ROS), often implicated in neurodegenerative diseases and aging. Moreover, biological experiments are designed to investigate how proteostasis depends on the balance between the folding capacity of chaperone networks and the continuous flux of potentially nonnative proteins. Nevertheless, biological experimental methods can examine the protein folding quality control mechanisms only in individual cells, but not in a multicellular level. Formal models offer a dynamic form of modelling, which allows to explore various parameter values in an integrated time-dependent system. This paper aims to present a formal approach of a mathematical descriptive model using as example a representation of a known molecular chaperone system and its relation to diseases associated to protein misfolding and neurodegeneration.

The molecular chaperone Hsp70 promotes the proteolytic removal of oxidatively damaged proteins by the proteasome

One hallmark of aging is the accumulation of protein aggregates, promoted by the unfolding of oxidized proteins. Unraveling the mechanism by which oxidized proteins are degraded may provide a basis to delay the early onset of features, such as protein aggregate formation, that contribute to the aging phenotype. In order to prevent aggregation of oxidized proteins, cells recur to the 20S proteasome, an efficient turnover proteolysis complex. It has previously been shown that upon oxidative stress the 26S proteasome, another form, dissociates into the 20S form. A critical player implicated in its dissociation is the Heat Shock Protein 70 (Hsp70), which promotes an increase in free 20S proteasome and, therefore, an increased capability to degrade oxidized proteins. The aim of this study was to test whether or not Hsp70 is involved in cooperating with the 20S proteasome for a selective degradation of oxidatively damaged proteins. Our results demonstrate that Hsp70 expression is induced in HT22 cells as a result of mild oxidative stress conditions. Furthermore, Hsp70 prevents the accumulation of oxidized proteins and directly promotes their degradation by the 20S proteasome. In contrast the expression of the Heat shock cognate protein 70 (Hsc70) was not changed in recovery after oxidative stress and Hsc70 has no influence on the removal of oxidatively damaged proteins. We were able to demonstrate in HT22 cells, in brain homogenates from 129/SV mice and in vitro, that there is an increased interaction of Hsp70 with oxidized proteins, but also with the 20S proteasome, indicating a role of Hsp70 in mediating the interaction of oxidized proteins with the 20S proteasome. Thus, our data clearly implicate an involvement of Hsp70 oxidatively damaged protein degradation by the 20S proteasome.

Aging and computational systems biology

Aging research is undergoing a paradigm shift, which has led to new and innovative methods of exploring this complex phenomenon. The systems biology approach endeavors to understand biological systems in a holistic manner, by taking account of intrinsic interactions, while also attempting to account for the impact of external inputs, such as diet. A key technique employed in systems biology is computational modeling, which involves mathematically describing and simulating the dynamics of biological systems. Although a large number of computational models have been developed in recent years, these models have focused on various discrete components of the aging process, and to date no model has succeeded in completely representing the full scope of aging. Combining existing models or developing new models may help to address this need and in so doing could help achieve an improved understanding of the intrinsic mechanisms which underpin aging.

Integration of 'omics' data in aging research: from biomarkers to systems biology

Age is the strongest risk factor for many diseases including neurodegenerative disorders, coronary heart disease, type 2 diabetes and cancer. Due to increasing life expectancy and low birth rates, the incidence of age-related diseases is increasing in industrialized countries. Therefore, understanding the relationship between diseases and aging and facilitating healthy aging are major goals in medical research. In the last decades, the dimension of biological data has drastically increased with high-throughput technologies now measuring thousands of (epi) genetic, expression and metabolic variables. The most common and so far successful approach to the analysis of these data is the so-called reductionist approach. It consists of separately testing each variable for association with the phenotype of interest such as age or age-related disease. However, a large portion of the observed phenotypic variance remains unexplained and a comprehensive understanding of most complex phenotypes is lacking. Systems biology aims to integrate data from different experiments to gain an understanding of the system as a whole rather than focusing on individual factors. It thus allows deeper insights into the mechanisms of complex traits, which are caused by the joint influence of several, interacting changes in the biological system. In this review, we look at the current progress of applying omics technologies to identify biomarkers of aging. We then survey existing systems biology approaches that allow for an integration of different types of data and highlight the need for further developments in this area to improve epidemiologic investigations.

Adaptive deployment of model reductions for tau-leaping simulation

Multiple time scales in cellular chemical reaction systems often render the tau-leaping algorithm inefficient. Various model reductions have been proposed to accelerate tau-leaping simulations. However, these are often identified and deployed manually, requiring expert knowledge. This is time-consuming and prone to error. In previous work, we proposed a methodology for automatic identification and validation of model reduction opportunities for tau-leaping simulation. Here, we show how the model reductions can be automatically and adaptively deployed during the time course of a simulation. For multiscale systems, this can result in substantial speedups.

BioModels: Content, Features, Functionality, and Use

BioModels is a reference repository hosting mathematical models that describe the dynamic interactions of biological components at various scales. The resource provides access to over 1,200 models described in literature and over 140,000 models automatically generated from pathway resources. Most model components are cross-linked with external resources to facilitate interoperability. A large proportion of models are manually curated to ensure reproducibility of simulation results. This tutorial presents BioModels' content, features, functionality, and usage.

The proteome profiles of the olfactory bulb of juvenile, adult and aged rats - an ontogenetic study

Background

In this study, we searched for proteins that change their expression in the olfactory bulb (oB) of rats during ontogenesis. Up to now, protein expression differences in the developing animal are not fully understood. Our investigation focused on the question whether specific proteins exist which are only expressed during different development stages. This might lead to a better characterization of the microenvironment and to a better determination of factors and candidates that influence the differentiation of neuronal progenitor cells.

Results

After analyzing the samples by two-dimensional polyacrylamide gel electrophoresis (2DE) and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS), it could be shown that the number of expressed proteins differs depending on the developmental stages. Especially members of the functional classes, like proteins of biosynthesis, regulatory proteins and structural proteins, show the highest differential expression in the stages of development analyzed.

Conclusion

In this study, quantitative changes in the expression of proteins in the oB at different developmental stages (postnatal days (P) 7, 90 and 637) could be observed. Furthermore, the expression of many proteins was found at specific developmental stages. It was possible to identify these proteins which are involved in processes like support of cell migration and differentiation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12953-014-0058-x) contains supplementary material, which is available to authorized users.

Understanding the role of heat shock protein isoforms in male fertility, aging and apoptosis

Heat shock proteins (HSPs) play a role in the homeostasis, apoptosis regulation and the maintenance of the various other physiological processes. Aging is accompanied by a decrease in the resistance to environmental stress, while mitochondria are primary targets in the process of aging, their expression decreasing with age. Mitochondrion also plays a significant role in the process of spermatogenesis. HSPs have been shown to be involved in apoptosis with some of acting as apoptotic inhibitors and are involved in cytoprotection. In this review we discuss the roles of Hsp 27, 60, 70, and 90 in aging and male infertility and have concluded that these particular HSPs can be used as a molecular markers for mitochondrially- mediated apoptosis, aging and male infertility.

Proteome analysis in the assessment of ageing

Based on demographic trends, the societies in many developed countries are facing an increasing number and proportion of people over the age of 65. The raise in elderly populations along with improved health-care will be concomitant with an increased prevalence of ageing-associated chronic conditions like cardiovascular, renal, and respiratory diseases, arthritis, dementia, and diabetes mellitus. This is expected to pose unprecedented challenges both for individuals and societies and their health care systems. An ultimate goal of ageing research is therefore the understanding of physiological ageing and the achievement of 'healthy' ageing by decreasing age-related pathologies. However, on a molecular level, ageing is a complex multi-mechanistic process whose contributing factors may vary individually, partly overlap with pathological alterations, and are often poorly understood. Proteome analysis potentially allows modelling of these multifactorial processes. This review summarises recent proteomic research on age-related changes identified in animal models and human studies. We combined this information with pathway analysis to identify molecular mechanisms associated with ageing. We identified some molecular pathways that are affected in most or even all organs and others that are organ-specific. However, appropriately powered studies are needed to confirm these findings based in in silico evaluation.

Age-correlated stress resistance improves fitness of yeast: support from agent-based simulations

BACKGROUND

Resistance to stress is often heterogeneous among individuals within a population, which helps protect against intermittent stress (bet hedging). This is also the case for heat shock resistance in the budding yeast Saccharomyces cerevisiae. Interestingly, the resistance appears to be continuously distributed (vs. binary, switch-like) and correlated with replicative age (vs. random). Older, slower-growing cells are more resistant than younger, faster-growing ones. Is there a fitness benefit to age-correlated stress resistance?

RESULTS

Here this hypothesis is explored using a simple agent-based model, which simulates a population of individual cells that grow and replicate. Cells age by accumulating damage, which lowers their growth rate. They synthesize trehalose at a metabolic cost, which helps protect against heat shock. Proteins Tsl1 and Tps3 (trehalose synthase complex regulatory subunit TSL1 and TPS3) represent the trehalose synthesis complex and they are expressed using constant, age-dependent and stochastic terms. The model was constrained by calibration and comparison to data from the literature, including individual-based observations obtained using high-throughput microscopy and flow cytometry. A heterogeneity network was developed, which highlights the predominant sources and pathways of resistance heterogeneity. To determine the best trehalose synthesis strategy, model strains with different Tsl1/Tps3 expression parameters were placed in competition in an environment with intermittent heat shocks.

CONCLUSIONS

For high severities and low frequencies of heat shock, the winning strain used an age-dependent bet hedging strategy, which shows that there can be a benefit to age-correlated stress resistance. The study also illustrates the utility of combining individual-based observations and modeling to understand mechanisms underlying population heterogeneity, and the effect on fitness.

Stochastic simulation of chemically reacting systems using multi-core processors

In recent years, computer simulations have become increasingly useful when trying to understand the complex dynamics of biochemical networks, particularly in stochastic systems. In such situations stochastic simulation is vital in gaining an understanding of the inherent stochasticity present, as these models are rarely analytically tractable. However, a stochastic approach can be computationally prohibitive for many models. A number of approximations have been proposed that aim to speed up stochastic simulations. However, the majority of these approaches are fundamentally serial in terms of central processing unit (CPU) usage. In this paper, we propose a novel simulation algorithm that utilises the potential of multi-core machines. This algorithm partitions the model into smaller sub-models. These sub-models are then simulated, in parallel, on separate CPUs. We demonstrate that this method is accurate and can speed-up the simulation by a factor proportional to the number of processors available.

Comparative physiology and transcriptional networks underlying the heat shock response in Populus trichocarpa, Arabidopsis thaliana and Glycine max

The heat shock response continues to be layered with additional complexity as interactions and crosstalk among heat shock proteins (HSPs), the reactive oxygen network and hormonal signalling are discovered. However, comparative analyses exploring variation in each of these processes among species remain relatively unexplored. In controlled environment experiments, photosynthetic response curves were conducted from 22 to 42 °C and indicated that temperature optimum of light-saturated photosynthesis was greater for Glycine max relative to Arabidopsis thaliana or Populus trichocarpa. Transcript profiles were taken at defined states along the temperature response curves, and inferred pathway analysis revealed species-specific variation in the abiotic stress and the minor carbohydrate raffinose/galactinol pathways. A weighted gene co-expression network approach was used to group individual genes into network modules linking biochemical measures of the antioxidant system to leaf-level photosynthesis among P. trichocarpa, G. max and A. thaliana. Network-enabled results revealed an expansion in the G. max HSP17 protein family and divergence in the regulation of the antioxidant and heat shock modules relative to P. trichocarpa and A. thaliana. These results indicate that although the heat shock response is highly conserved, there is considerable species-specific variation in its regulation.

Modelling the role of the Hsp70/Hsp90 system in the maintenance of protein homeostasis

Neurodegeneration is an age-related disorder which is characterised by the accumulation of aggregated protein and neuronal cell death. There are many different neurodegenerative diseases which are classified according to the specific proteins involved and the regions of the brain which are affected. Despite individual differences, there are common mechanisms at the sub-cellular level leading to loss of protein homeostasis. The two central systems in protein homeostasis are the chaperone system, which promotes correct protein folding, and the cellular proteolytic system, which degrades misfolded or damaged proteins. Since these systems and their interactions are very complex, we use mathematical modelling to aid understanding of the processes involved. The model developed in this study focuses on the role of Hsp70 (IPR00103) and Hsp90 (IPR001404) chaperones in preventing both protein aggregation and cell death. Simulations were performed under three different conditions: no stress; transient stress due to an increase in reactive oxygen species; and high stress due to sustained increases in reactive oxygen species. The model predicts that protein homeostasis can be maintained during short periods of stress. However, under long periods of stress, the chaperone system becomes overwhelmed and the probability of cell death pathways being activated increases. Simulations were also run in which cell death mediated by the JNK (P45983) and p38 (Q16539) pathways was inhibited. The model predicts that inhibiting either or both of these pathways may delay cell death but does not stop the aggregation process and that eventually cells die due to aggregated protein inhibiting proteasomal function. This problem can be overcome if the sequestration of aggregated protein into inclusion bodies is enhanced. This model predicts responses to reactive oxygen species-mediated stress that are consistent with currently available experimental data. The model can be used to assess specific interventions to reduce cell death due to impaired protein homeostasis.

Systems biology of ageing and longevity

Ageing is intrinsically complex, being driven by multiple causal mechanisms. Each mechanism tends to be partially supported by data indicating that it has a role in the overall cellular and molecular pathways underlying the ageing process. However, the magnitude of this role is usually modest. The systems biology approach combines (i) data-driven modelling, often using the large volumes of data generated by functional genomics technologies, and (ii) hypothesis-driven experimental studies to investigate causal pathways and identify their parameter values in an unusually quantitative manner, which enables the contributions of individual mechanisms and their interactions to be better understood, and allows for the design of experiments explicitly to test the complex predictions arising from such models. A clear example of the success of the systems biology approach in unravelling the complexity of ageing can be seen in recent studies on cell replicative senescence, revealing interactions between mitochondrial dysfunction, telomere erosion and DNA damage. An important challenge also exists in connecting the network of (random) damage-driven proximate mechanisms of ageing with the higher level (genetically specified) signalling pathways that influence longevity. This connection is informed by actions of natural selection on the determinants of ageing and longevity.

Stiffness detection and reduction in discrete stochastic simulation of biochemical systems

Typical multiscale biochemical models contain fast-scale and slow-scale reactions, where "fast" reactions fire much more frequently than "slow" ones. This feature often causes stiffness in discrete stochastic simulation methods such as Gillespie's algorithm and the Tau-Leaping method leading to inefficient simulation. This paper proposes a new strategy to automatically detect stiffness and identify species that cause stiffness for the Tau-Leaping method, as well as two stiffness reduction methods. Numerical results on a stiff decaying dimerization model and a heat shock protein regulation model demonstrate the efficiency and accuracy of the proposed methods for multiscale biochemical systems.

Efficient formulations for exact stochastic simulation of chemical systems

One can generate trajectories to simulate a system of chemical reactions using either Gillespie's direct method or Gibson and Bruck's next reaction method. Because one usually needs many trajectories to understand the dynamics of a system, performance is important. In this paper, we present new formulations of these methods that improve the computational complexity of the algorithms. We present optimized implementations, available from http://cain.sourceforge.net/, that offer better performance than previous work. There is no single method that is best for all problems. Simple formulations often work best for systems with a small number of reactions, while some sophisticated methods offer the best performance for large problems and scale well asymptotically. We investigate the performance of each formulation on simple biological systems using a wide range of problem sizes. We also consider the numerical accuracy of the direct and the next reaction method. We have found that special precautions must be taken in order to ensure that randomness is not discarded during the course of a simulation.

Computational systems biology of aging

Computational systems biology is expected to make major contributions to unravel the complex molecular mechanisms underlying the progression of aging in cells, tissues, and organisms. The development of computational approaches is, however, challenged by a wide spectrum of aging mechanisms participating on different levels of biological organization. The tight connectivity between the molecular constituents, functions, and cell states requires frameworks and strategies that extend beyond current practice to model, simulate, and predict the progression of aging and the emerging aging phenotypes. We provide a general overview of the specific computational tasks and opportunities in aging research, and discuss some illustrative systems level concepts in more detail. One example provided here is the assembly of a conceptual whole cell model that considers the temporal dynamics of the aging process grounded on molecular mechanisms. Another application is the assembly of interactomes, such as protein networks that allow us to analyze changes in network topology and interaction of proteins that have been implicated in aging with other cellular constituents and processes. We introduce the necessary key steps to build these applications and discuss their merits and future extensions for aging research. WIREs Syst Biol Med 2011 3 414-428 DOI: 10.1002/wsbm.126

Global aging and the brain

Intrinsic aging of the brain and its relationship to age-related neurodegenerative diseases need to be understood as part of global aging, which results from gradual accumulation of a variety of kinds of cellular and molecular damage. While certain kinds of molecular lesions are particularly associated with specific diseases, there is substantial overlap between the pathways causing disease and those contributing to global aging. Growing evidence indicates that the global aging process is more malleable than used to be thought. This needs to be taken into account in efforts to improve health and retention of mental capital across the life course.

Rule-based cell systems model of aging using feedback loop motifs mediated by stress responses

Investigating the complex systems dynamics of the aging process requires integration of a broad range of cellular processes describing damage and functional decline co-existing with adaptive and protective regulatory mechanisms. We evolve an integrated generic cell network to represent the connectivity of key cellular mechanisms structured into positive and negative feedback loop motifs centrally important for aging. The conceptual network is casted into a fuzzy-logic, hybrid-intelligent framework based on interaction rules assembled from a priori knowledge. Based upon a classical homeostatic representation of cellular energy metabolism, we first demonstrate how positive-feedback loops accelerate damage and decline consistent with a vicious cycle. This model is iteratively extended towards an adaptive response model by incorporating protective negative-feedback loop circuits. Time-lapse simulations of the adaptive response model uncover how transcriptional and translational changes, mediated by stress sensors NF-κB and mTOR, counteract accumulating damage and dysfunction by modulating mitochondrial respiration, metabolic fluxes, biosynthesis, and autophagy, crucial for cellular survival. The model allows consideration of lifespan optimization scenarios with respect to fitness criteria using a sensitivity analysis. Our work establishes a novel extendable and scalable computational approach capable to connect tractable molecular mechanisms with cellular network dynamics underlying the emerging aging phenotype.

The global process of aging disturbs a broad range of cellular mechanisms in a complex fashion and is not well understood. One important goal of computational approaches in aging is to develop integrated models in terms of a unifying aging theory, predicting progression of aging phenotypes grounded on molecular mechanisms. However, current experimental data incoherently reflects many isolated processes from a large diversity of approaches, biological model systems, and species, which makes such integration a challenging task. In an attempt to close this gap, we iteratively develop a fuzzy-logic cell systems model considering the interplay of damage, metabolism, and signaling by positive and negative feedback-loop motifs using relationships drawn from literature data. Because cellular biodynamics may be considered a complex control system, this approach seems particularly suitable. Here, we demonstrate that rule-based fuzzy-logic models provide semi-quantitative predictions that enhance our understanding of complex and interlocked molecular mechanisms and their implications on the aging physiome.

Specific regulation of noncanonical p38alpha activation by Hsp90-Cdc37 chaperone complex in cardiomyocyte

RATIONALE

p38 is an important stress activated protein kinase involved in gene regulation, proliferation, differentiation, and cell death regulation in heart. p38 kinase activity can be induced through canonical pathway via upstream kinases or by noncanonical autophosphorylation. The intracellular p38 kinase activity is tightly regulated and maintained at low level under basal condition. The underlying regulatory mechanism for canonical p38 kinase activation is well-studied, but the regulation of noncanonical p38 autophosphorylation remains poorly understood.

OBJECTIVE

We investigated the molecular basis for the regulation of noncanonical p38 autophosphorylation and its potential functional impact in cardiomyocytes.

METHODS AND RESULTS

Using both proteomic and biochemical tools, we established that heat shock protein (Hsp)90-Cdc37 chaperones are part of the p38alpha signaling complex in mammalian cells both in vitro and in vivo. The Hsp90-Cdc37 chaperone complex interacts with p38 via direct binding between p38 and Cdc37. Cdc37 expression is both sufficient and necessary to suppress noncanonical p38 activation via autophosphorylation at either basal state or under TAB1 (TAK1 binding protein-1) induction. In contrast, Cdc37 expression has no impact on p38 activation by canonical upstream kinase MKK3 or oxidative stress. Furthermore, Hsp90 inhibition results in p38 activation via autophosphorylation, and p38 activity contribute to apoptotic cell death induced by Hsp90 inhibition.

CONCLUSION

Our study has revealed a so far uncharacterized function of Hsp90-Cdc37 as an endogenous regulator of noncanonical p38 activity.

Protein aggregation as a paradigm of aging

The process of physiological decline leading to death of the individual is driven by the deteriorating capacity to withstand extrinsic and intrinsic hazards, resulting in damage accumulation with age. The dynamic changes with time of the network governing the outcome of misfolded proteins, exemplifying as intrinsic hazards, is considered here as a paradigm of aging. The main features of the network, namely, the non-linear increase of damage and the presence of amplifying feedback loops within the system are presented through a survey of the different components of the network and related cellular processes in aging and disease.

HSP 70 and atherosclerosis--protector or activator?

BACKGROUND

Atherosclerosis and its complications represent the leading cause of morbidity and mortality. Heat shock protein 70 (HSP70) protects cellular elements from injury by reducing oxidation, inflammation and apoptosis and by refolding damaged proteins. HSP70 improves viability of stressed vascular smooth muscle cells, possibly via its chaperone functions. It has been proposed that the response mounted against bacterial HSPs results in an autoimmune reaction, which has the potential to cause complement-mediated endothelial injury, and hence accelerate atherogenesis.

OBJECTIVE

to examine the roles of HSPs in atherosclerosis.

METHODS

A literature review.

RESULTS/CONCLUSIONS

The role of HSPs in atherosclerosis is controversial. HSP60 probably acts as an autoantigen, and may trigger both cell- and antibody-mediated immune responses, while HSP70 is likely to be involved in cytoprotection. The significance of this inverse relation between HSP70 and atherosclerosis has not yet been elucidated. Whether HSPs will become therapeutic targets remains to be established.

Stochastic modelling for quantitative description of heterogeneous biological systems

Two related developments are currently changing traditional approaches to computational systems biology modelling. First, stochastic models are being used increasingly in preference to deterministic models to describe biochemical network dynamics at the single-cell level. Second, sophisticated statistical methods and algorithms are being used to fit both deterministic and stochastic models to time course and other experimental data. Both frameworks are needed to adequately describe observed noise, variability and heterogeneity of biological systems over a range of scales of biological organization.

Differential proteomics in the aging Noble rat ventral prostate

Incidence of prostatic diseases increases dramatically with age which may be related to a decline in androgen support. However, the key mechanisms underlying prostate aging remain unclear. In the present study, we investigated the aging process in the ventral prostate (VP) of Noble rats by identifying differentially expressed prostate proteins between 3- and 16-month-old animals using ICAT and MS. In total, 472 proteins were identified with less than a 1% false positive rate, among which 34 were determined to have a greater than two-fold increase or 1.7-fold decrease in expression in the aged VPs versus their younger counterparts. The majority of the differentially expressed proteins identified have not been previously reported to be associated with prostate aging, and they fall into specific functional categories, including oxidative stress/detoxification, chaperones, protein biosynthesis, vesicle transport, and intracellular trafficking. The expression of GST, ferritin, clusterin, kininogen, oxygen regulated protein 150, spermidine synthase, ADP ribosylation factor, and cyclophilin B was verified by Western blot analyses on samples used for the ICAT study, as well as on those obtained from an independent group of animals comprised of three age groups. To the best of our knowledge, this is the first study on the proteome of the aging rat prostate.

Gene-expression profiles and age of donor kidney biopsies obtained before transplantation distinguish medium term graft function

BACKGROUND

Donor factors such as age profoundly influence long-term graft function after cadaveric renal transplantation, but the molecular signature of these aspects in the allograft remains unknown.

METHODS

We analyzed the genome-wide gene expression signature of donor kidney biopsies of different ages obtained before transplantation. Subsequent analysis compared expression profiles from allografts with excellent function versus impaired function at 1 yr after engraftment. Differential expression profiles were analyzed on the level of molecular function and biologic role, as well as by analysis of co-regulation through transcription factors, regulatory networks, and protein-protein interaction data utilizing extended bioinformatics.

RESULTS

The 15 subjects with excellent transplant function defined as calculated GFR>or=45 mL/min/1.73 m2 at 1 yr exhibited a distinctly different gene expression profile than the matched 16 subjects with impaired function defined as calculated GFR<45 mL/min/1.73 m2. Donor kidneys from recipients with impaired allograft function showed activation of genes mainly belonging to the functional classes of immunity, signal transduction, and oxidative stress response. Two-thirds of these genes exhibited at least one protein interacting partner, suggesting choreographed intracellular events differentiating the two recipient groups. However, donor age may have confounded some of the associations found between gene profiles and graft function.

CONCLUSION

In summary, a distinctive gene expression profile in the donor kidney at transplantation together with donor age predicts medium term allograft function in recipients of cadaveric allografts.

Stress chaperones, mortalin, and pex19p mediate 5-aza-2' deoxycytidine-induced senescence of cancer cells by DNA methylation-independent pathway

DNA demethylating agents are used to reverse epigenetic silencing of tumor suppressors in cancer therapeutics. Understanding of the molecular and cellular factors involved in DNA demethylation-induced gene desilencing and senescence is still limited. We have tested the involvement of two stress chaperones, Pex19p and mortalin, in 5-Aza-2' deoxycytidine (5AZA-dC; DNA demethylating agent)-induced senescence. We found that the cells overexpressing these chaperones were highly sensitive to 5AZA-dC, and their partial silencing eliminated 5AZA-dC-induced senescence in human osteosarcoma cells. We demonstrate that these chaperones modulate the demethylation and chromatin remodeling-dependent (as accessed by p16(INK4A) expression) and remodeling-independent (such as activation of tumor suppressor p53 pathway) senescence response of cells. Furthermore, we found the direct interactions of 5AZA-dC with these chaperones that may alter their functions. We conclude that both mortalin and Pex19p are important mediators, prognostic indicators, and tailoring tools for 5AZA-dC-induced senescence in cancer cells.

An in silico model of the ubiquitin-proteasome system that incorporates normal homeostasis and age-related decline

Background

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.

Results

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.

Conclusion

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.

On why decreasing protein synthesis can increase lifespan

An explanation is offered for the increased lifespan of Caenorhabditis elegans when mRNA translation is inhibited due to loss of the initiation factor IFE-2 [Hansen, M., Taubert, T., Crawford, D., Libina, N., Lee, S.-J., Kenyon, C., 2007. Lifespan extension by conditions that inhibit translation in Caenorhabditis elegans. Ageing Cell 6, 95-110; Pan, K.Z., Palter, J.E., Rogers, A.N., Olsen, A., Chen, D., Lithgow, G.J., Kapahi, P., 2007. Inhibition of mRNA translation extends lifespan in Caenorhabditis elegans. Ageing Cell 6, 111-119; Syntichaki, P., Troulinaki, K., Tavernarakis, N., 2007. eIF4E function in somatic cells modulates ageing in Caenorhabditis elegans. Nature 445, 922-926]. It is suggested that the general reduction of protein synthesis, due to the decreased frequency of mRNA translation, also lowers the cellular load of erroneously synthesized polypeptides which the constitutive protein homeostatic apparatus (proteases and chaperones proteins) normally eliminates. This situation results in "spare" proteolytic and chaperone function which can then deal with those proteins modified post-synthetically, e.g. by oxidation and/or glycation, which are thought to contribute to the senescent phenotype. This increased availability of proteolytic and chaperone functions may thereby contribute to the observed increase in organism stress resistance and lifespan.

The use of conformation-specific ligands and assays to dissect the molecular mechanisms of neurodegenerative diseases

The use of conformation-specific ligands has been closely linked to progress in the molecular characterization of neurodegenerative diseases. Deposition of misfolded or misprocessed proteins is now recognized as a hallmark of all neurodegenerative diseases. Initially, dyes like Congo red and thioflavin T were used as crudely conformation-specific ligands for staining the beta-sheeted protein components of amyloid deposits in neurodegenerative diseases such as Alzheimer disease (AD) and prion disease, the two diseases in which protein conformations were distinguished early on. This conformational characterization of extracellular protein deposits with dyes ultimately led to the identification of key players in the disease processes. The recent discovery of intermediate conformational species, i.e., soluble oligomers for AD and PK-sensitive PrP(Sc) for prion disease, whose conformation and assembly are thought to be distinct from both the physiological and the fibrillar conformational states, replaced the former notion that the microscopic protein deposits themselves caused disease. This insight and the generation of conformation-specific monoclonal antibodies to these conformers further advanced diagnosis and the understanding of molecular mechanisms of AD and are likely to do so in other neurodegenerative diseases. Here we review how conformer distinction performed by a variety of different techniques, including biophysical, biochemical, and antibody-based methods, led to the current molecular concepts of AD and the prion diseases. We provide an outlook on the application of these techniques in advancing the understanding of molecular mechanisms of other neurodegenerative diseases or degenerative brain conditions.

Transcriptional profiling of Alzheimer blood mononuclear cells by microarray

We evaluated pathomechanisms and systemic manifestations of Alzheimer disease (AD), an aging-related dementing neurodegenerative disorder, by expression profiling. Blood mononuclear cell (BMC) transcriptomes of sporadic AD subjects and aged-matched normal elderly controls (NEC) were compared using the human NIA microarray. Relative to the NEC samples, the Alzheimer BMC exhibited a significant decline in the expression of genes concerned with cytoskeletal maintenance, cellular trafficking, cellular stress response, redox homeostasis, transcription and DNA repair. We observed decreased expression of several genes which may impact amyloid-beta production and the processing of the microtubule-associated protein tau. The microarray results were validated by quantitative real time PCR and revealed gender differences in the levels of altered gene expression. Our findings attest to the systemic nature of gene dys-regulation in sporadic AD, implicate disruption of cytoskeletal integrity, DNA repair mechanisms and cellular defenses in this condition, and suggest novel pathways of beta-amyloid deposition in this disease. BMC are highly accessible and may reflect molecular events germane to the neuropathophysiology of AD.

PML clastosomes prevent nuclear accumulation of mutant ataxin-7 and other polyglutamine proteins

The pathogenesis of spinocerebellar ataxia type 7 and other neurodegenerative polyglutamine (polyQ) disorders correlates with the aberrant accumulation of toxic polyQ-expanded proteins in the nucleus. Promyelocytic leukemia protein (PML) nuclear bodies are often present in polyQ aggregates, but their relation to pathogenesis is unclear. We show that expression of PML isoform IV leads to the formation of distinct nuclear bodies enriched in components of the ubiquitin-proteasome system. These bodies recruit soluble mutant ataxin-7 and promote its degradation by proteasome-dependent proteolysis, thus preventing the aggregate formation. Inversely, disruption of the endogenous nuclear bodies with cadmium increases the nuclear accumulation and aggregation of mutant ataxin-7, demonstrating their role in ataxin-7 turnover. Interestingly, β-interferon treatment, which induces the expression of endogenous PML IV, prevents the accumulation of transiently expressed mutant ataxin-7 without affecting the level of the endogenous wild-type protein. Therefore, clastosomes represent a potential therapeutic target for preventing polyQ disorders.

Tools for the SBML Community

MOTIVATION

SBML is quickly becoming the standard format to exchange biochemical models. The tools presented in this paper are loosely-coupled, and are intended to be incorporated into SBML aware applications. The rationale for this is to reduce the amount of repeated work carried out within the community and to create tools that offer a greater number of features to the end-user.

AVAILABILITY

All tools described are available from http://www.basis.ncl.ac.uk/software and are licensed under GNU General Public License.