Signaling pathway cloud regulation for in silico screening and ranking of the potential geroprotective drugs
- 1Department of Biological and Medical Physics, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
- 2First Oncology Research and Advisory Center, Moscow, Russia
- 3The Biogerontology Research Foundation, London, UK
- 4Department of Experimental and Molecular Medicine, D. Rogachyov Federal Research Center of Pediatric Hematology, Oncology and Immunology, Moscow, Russia
- 5Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Moscow, Russia
- 6Burnasyan Federal Medical Biophysical Center, Moscow, Russia
- 7Department of Ecology, Syktyvkar State University, Syktyvkar, Russia
- 8Laboratory of Molecular Radiobiology and Gerontology, Institute of Biology of Komi Science Center of Ural Branch of Russian Academy of Sciences, Syktyvkar, Russia
The major challenges of aging research include absence of the comprehensive set of aging biomarkers, the time it takes to evaluate the effects of various interventions on longevity in humans and the difficulty extrapolating the results from model organisms to humans. To address these challenges we propose the in silico method for screening and ranking the possible geroprotectors followed by the high-throughput in vivo and in vitro validation. The proposed method evaluates the changes in the collection of activated or suppressed signaling pathways involved in aging and longevity, termed signaling pathway cloud, constructed using the gene expression data and epigenetic profiles of young and old patients' tissues. The possible interventions are selected and rated according to their ability to regulate age-related changes and minimize differences in the signaling pathway cloud. While many algorithmic solutions to simulating the induction of the old into young metabolic profiles in silico are possible, this flexible and scalable approach may potentially be used to predict the efficacy of the many drugs that may extend human longevity before conducting pre-clinical work and expensive clinical trials.
The increasing burden of the aging on the economies of the developed countries is turning the quest to increase healthy life spans from an altruistic cause into a pressing economic priority required to maintain the current standards of living and facilitate economic growth (Zhavoronkov and Litovchenko, 2013). There is an urgent need to develop and validate interventions with geroprotective properties to increase the productive health spans of the working population and maintaining performance and avoiding loss of function (Kennedy, 2012).
While no doubt exists that aging is a complex multifactorial process with no single cause or treatment (Zhavoronkov and Cantor, 2011; Trindade et al., 2013), the issue whether aging can be classified as the disease is widely debated (Rattan, 2013). However, many strategies for extending organismal life spans have been proposed including replacing cells (Rodgerson and Harris, 2011) and organs, comprehensive strategies for repairing the accumulated damage, using hormetins to activate endogenous repair processes (Gems and Partridge, 2008; Gaman et al., 2011), modulating the aging processes through specific mutations, gene therapy (Bernardes De Jesus et al., 2012) and small molecule drugs (Kennedy and Pennypacker, 2013). An animal's survival strongly depends on its ability to maintain homeostasis and achieved through intracellular and intercellular communication within and among different tissues (Alcedo et al., 2013). Many strategies for the development and validation of drugs with geroprotective properties have been proposed to help maintain the homeostasis including drugs that act on specific targets or combinations of molecular pathways (Moskalev and Shaposhnikov, 2010, 2011; Zhavoronkov et al., 2012; Danilov et al., 2013) and epigenetic drugs (Vaiserman and Pasyukova, 2012). However, none of the proposed strategies for aging-suppressive drug development provide a roadmap for rapid screening, validation, and clinical deployment. No methods currently exist to predict the effects of currently available drugs on human longevity and health span in a timely manner. This is partly due to the absence of the clear panel of human processes involved in aging to effectively run clinical trials.
Many processes are involved in the aging of cells and organisms including telomere length (Lehmann et al., 2013), intracellular and extracellular aggregates, racemization of the amino acids and genetic instability. Both gene expression (Wolters and Schumacher, 2013) and DNA methylation profiles (Horvath et al., 2012; Horvath, 2013; Mendelsohn and Larrick, 2013) change during aging and may be used as biomarkers of aging. Many studies analyzing transcriptomes of biopsies in a variety of diseases indicated that age and sex of the patient had significant effects on gene expression (Chowers et al., 2003) and that there are noticeable changes in gene expression with age in mice (Weindruch et al., 2002; Park et al., 2009) resulting in development of mouse aging gene expression databases (Zahn et al., 2007) and in humans (Blalock et al., 2003; Welle et al., 2003; Park and Prolla, 2005; Hong et al., 2008; De Magalhaes et al., 2009).
Combination of protein-protein interaction and gene expression in both flies and humans demonstrated that aging is mainly associated with a small number of biological processes, might preferentially attack key regulatory nodes that are important for network stability (Xue et al., 2007).
Our prior work with gene expression and epigenetics of various solid tumors (Kuzmin et al., 2010; Mityaev et al., 2010; Zabolotneva et al., 2012a,b) using the OncoFinder system (www.oncofinder.com), provided clues that transcription profiles of cancer cells mapped onto the signaling pathways may be used to screen for and rate the targeted drugs that regulate pathways directly and indirectly related to aging and longevity. Instead of focusing on individual network elements, this approach involves creating the signaling pathway cloud, a collection of signaling pathways involved in aging and longevity each comprised of multiple network elements and evaluating the individual pathway activation strength. Despite significant advances in aging research, the knowledge of the aging processes is still poor, and combining all available factors involved in cellular aging, aging of the organisms, age-related diseases, stress-resistance, and stress-response along the many other factors into a comprehensive signaling pathway cloud may be more beneficial than focusing on the narrow collection of elements. The creation of the pathway cloud may allow for the annotated databases of molecules and other factors to be screened for effectiveness of individual compounds in replicating the “young” signaling activation profiles in silico.
Several new methods evaluating the robustness and response ability of the gene regulatory network have been developed and applied to gene expression data sets from young and old patients (Tu and Chen, 2013). Prior studies suggested that a combination of pathways, termed pathway cloud, instead of one element of the pathway or the whole pathway might be responsible for pathological changes in the cell (Voronkov and Krauss, 2013). Long-lived species like the sea urchin (Strongylocentrotus franciscanus) and naked mole rat (Heterocephalus glaber) that senesce at a slower rate than members of the same order show less transcriptome changes with age (Kim et al., 2011; Loram and Bodnar, 2012). Gene network analysis using gene expression data was effective in identifying the possible drug targets (Imoto et al., 2003; Savoie et al., 2003). In silico drug discovery algorithms that attempt to transform the metabolism to the healthy state have been proposed and validated (Yizhak et al., 2013). We theorize that in order to be effective, the geroprotector or a combination of aging-suppressive drugs must regulate the pathway cloud in a way that minimizes the difference in the net differences in pathway cloud activation or downregulation between samples of young and old patients. Small molecules and other factors that may influence gene expression may be ranked by their ability to minimize the net difference between the pathway activation profiles of young and old cells. The algorithms for calculating the ability of the potential geroprotector to minimize signaling disturbance may be parametric and account for the effects on specific targets within signaling pathways or machine learned.
Despite the differences in life span and aging phenotypes, many molecular mechanisms of aging are common in all eukaryotes. Pathway analysis revealed that there are many common age-related transcriptomic changes between different species, including yeasts, worms, flies, rodents, and human (Murphy et al., 2006). Hypothetically, the human orthologs of aging-related genes of model organisms are also involved in aging process. To select longevity-associated pathways for future analysis we perform the following procedure. Using (Tacutu et al., 2013) database we selected genes, where knockout, loss-of-function mutation, deletion or RNA interference significantly extended lifespan in several model organisms (yeasts Schizosaccharomyces pombe and Saccharomyces cerevisiae, nematode Caenorhabditis elegans, fruitfly Drosophila melanogaster and mouse Mus musculus) from 10 to 200%. We converted the obtained gene lists from different models to the general list of human orthologs where it is possible. 226 genes of 315 from our set were subjected to over-representation pathway analysis in (Kamburov et al., 2011). P-value and corrected by FDR P-value are calculated according to the hypergeometric test based on the number of physical entities present in both the predefined set of genes to each (Kanehisa et al., 2004) pathway and our list of aging-associated genes (Table 1). We established limits of 2 genes of minimum overlap with input list and 0.01 p-value cut off threshold.
As a result we revealed overrepresented cell signaling pathways (mTOR, insulin/IGF-1, PI3K-Akt, PPAR, HIF-1, TGF-beta, chemokines, adipocytokine, prolactin, estrogen), general metabolism (TCA cycle, ribosome, oxidative phosphorylation), RNA transport, cell cycle and meiosis, gap junction, peroxisome, cyrcadian rhythm, different synapse types (dopaminergic, glutamatergic, cholinergic, serotonergic, GABAergic), gastric acid secretion as well as age-related diseases pathways (Parkinson's disease, type II diabetes mellitus, Huntington's disease, long-term depression, amyotrophic lateral sclerosis, Alzheimer's disease), Hepatitis B, HTLV-I infection and cancer pathways (prostate cancer, colorectal cancer, glioma, pancreatic cancer, chronic myeloid leukemia, proteoglycans in cancer). We considered obtained such a way pathways as probably associated with the human longevity. Human genes known as key activators/repressors of these pathways may be used in provided further mathematical model.
The methods that may be applied for the possible analysis of geroprotector efficiency by pathways regulation have been arisen from our research experience of cell signaling pathways. As far as we have seen before (Kiyatkin et al., 2006; Borisov et al., 2009; Kuzmina and Borisov, 2011), most signal transduction proteins are essentially far from saturation even at the peak concentrations of the activated form in comparison with the total protein abundances. Thus, we can consider that all activator/repressor genes/proteins have equal importance for the pathway activation/downregulation, and then arrive at the following assessment function for the overall signal pathway cloud disturbance outcome (SPCD) is proportional to the following estimator function,
Here the multiplication is done over all possible activator and repressor proteins in the pathway, and [AGEL]i and [RGEL]j are gene expression levels of an activator i and repressor j, respectively. To obtain an additive rather than multiplicative value, it is enough just turn from the absolute values of the expression levels to their logarithms, arriving at the pathway activation strength (PAS) value for each pathway (see Figure 1). To obtain the values of Old (case)-to-Young ratio, YORn, one just has to divide the expression levels for a gene n in the sample taken for the senescent person by the same average value for the normalized young group. The discrete value of ARR (activator/repressor role) equals to the following numbers:
Figure 1. Gene expression-based approach to in silico screening for drugs with geroprotective properties and estimating the predicted efficacy. (A) Using signaling pathway cloud regulation for theoretical in silico aging-suppressive drug identification and ranking. The proposed method for identifying and ranking of geroprotective drugs by evaluating the net effect on the many elements of signaling pathway cloud that brings the “old” metabolic state closer to the “young.” (B) An example of how multiple pathways are activated and down-regulated during aging. (C) Pathway Activation Strength (PAS) is the logarithmic additive value that characterizes the up-/downregulation of signaling pathways. (D) Function for the overall signal pathway cloud disturbance outcome (SPCD).
−1, when the gene/protein n is a repressor of pathway excitation;
1, if the gene/protein n is an activator of pathway excitation;
0, when the gene/protein n can be both an activator and a repressor of signal transduction;
0.5 and −0.5, respectively, if the gene/protein n is more an activator or repressor of the signaling pathway p.
The information about the activator/repressor role of a particular gene/protein may be obtained from the analysis of open-access or customized pathway databases and from the literature.
The Boolean flag of BTIF (beyond tolerance interval flag) equals to zero when the OYR value lies within the tolerance limit, and to one when otherwise. During the current study, we have admitted that the OYR lies beyond the tolerance limit if it satisfies simultaneously the two criteria. First, it either higher than 3/2 or lower than 2/3, and, second, the expression level for a corresponding gene from an old patient of an individual patient differs by more than two standard deviations from the average expression level for the same gene from a set of analogous young tissue/organ samples.
We propose a new computational approach for identifying and rating the variety of factors including small molecules, peptides, stress factors and conditions with the known effects on the transcriptomes at different ages of one or more cell or tissue types or known targets (Figure 1A). The approach may be used for general geroprotector screening, but after the validation of the algorithms in vivo and in vitro may be expanded to identify and predict the efficacy of personalized aging-suppressive intervention regimens for individual patients based on the transcriptome information from various tissue biopsies and blood samples.
The generic geroprotector rating approach involves collecting the transcriptome data sets from young and old patients and normalizing the data for each cell and tissue type, evaluating the pathway activation strength (PAS) for each individual pathway (Figure 1B) and constructing the pathway cloud (PC, Figure 1C) and screen for drugs or combinations that minimize the signaling pathway cloud disturbance (SPCD, Figure 1D) by acting on one or multiple elements of the pathway cloud. Drugs and combinations may be rated by their ability to compensate the changes in signaling pathway activation patterns that are related to aging, thus bringing all the PAS values for the general set of the pathways as close to zero as possible. Since many of the drugs approved for use in humans have known molecular targets and some have been screened for the impact on longevity in model organisms (Ye et al., 2013), the predictions may be then tested both in vitro and in vivo on human cells and on model organisms such as rodents, nematodes and flies to validate the screening and rating algorithms.
Longevity studies of aging-suppressive drug efficiency in higher mammals take several years and decades and may cost millions of dollars. An intelligent process for predicting the activity and ranking the geroprotective activity of various factors and strengthening the prediction in rapid and cost-effective studies on cell cultures and model organisms may help increase the longevity dividend of these studies. In this paper we propose a method for in silico screening and ranking of drugs and other factors that act on many signaling pathways implicated in aging processes by calculating their ability to minimize the difference between signaling pathway activation patterns in cells of young and old patients and confirming the results using in vivo and in vitro studies.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank the reviewers for the many useful comments, corrections, and recommendations. We would like to thank Charles R. Cantor, Brian Kennedy, Robert Shmookler-Reis and many others for helpful comments and discussions on strategies for geroprotector screening. We thank UMA Foundation for their help in preparation of this manuscript and Dr. Kristen Swithers from Yale University for her assistance with editing the manuscript. We would like to thank Alex Kim and ASUSTek for equipment and support of this research.
Bernardes De Jesus, B., Vera, E., Schneeberger, K., Tejera, A. M., Ayuso, E., Bosch, F., et al. (2012). Telomerase gene therapy in adult and old mice delays aging and increases longevity without increasing cancer. EMBO Mol. Med. 4, 691–704. doi: 10.1002/emmm.201200245
Blalock, E. M., Chen, K. C., Sharrow, K., Herman, J. P., Porter, N. M., Foster, T. C., et al. (2003). Gene microarrays in hippocampal aging: statistical profiling identifies novel processes correlated with cognitive impairment. J. Neurosci. 23, 3807–3819.
Borisov, N., Aksamitiene, E., Kiyatkin, A., Legewie, S., Berkhout, J., Maiwald, T., et al. (2009). Systems-level interactions between insulin-EGF networks amplify mitogenic signaling. Mol. Syst. Biol. 5, 256. doi: 10.1038/msb.2009.19
Chowers, I., Liu, D., Farkas, R. H., Gunatilaka, T. L., Hackam, A. S., Bernstein, S. L., et al. (2003). Gene expression variation in the adult human retina. Hum. Mol. Genet. 12, 2881–2893. doi: 10.1093/hmg/ddg326
De Magalhaes, J. P., Curado, J., and Church, G. M. (2009). Meta-analysis of age-related gene expression profiles identifies common signatures of aging. Bioinformatics 25, 875–881. doi: 10.1093/bioinformatics/btp073
Horvath, S., Zhang, Y., Langfelder, P., Kahn, R. S., Boks, M. P., Van Eijk, K., et al. (2012). Aging effects on DNA methylation modules in human brain and blood tissue. Genome Biol. 13:R97. doi: 10.1186/gb-2012-13-10-r97
Imoto, S., Savoie, C. J., Aburatani, S., Kim, S., Tashiro, K., Kuhara, S., et al. (2003). Use of gene networks for identifying and validating drug targets. J. Bioinform. Comput. Biol. 1, 459–474. doi: 10.1142/S0219720003000290
Kamburov, A., Pentchev, K., Galicka, H., Wierling, C., Lehrach, H., and Herwig, R. (2011). ConsensusPathDB: toward a more complete picture of cell biology. Nucleic Acids Res. 39, D712–D717. doi: 10.1093/nar/gkq1156
Kim, E. B., Fang, X., Fushan, A. A., Huang, Z., Lobanov, A. V., Han, L., et al. (2011). Genome sequencing reveals insights into physiology and longevity of the naked mole rat. Nature 479, 223–227. doi: 10.1038/nature10533
Kiyatkin, A., Aksamitiene, E., Markevich, N. I., Borisov, N. M., Hoek, J. B., and Kholodenko, B. N. (2006). Scaffolding protein Grb2-associated binder 1 sustains epidermal growth factor-induced mitogenic and survival signaling by multiple positive feedback loops. J. Biol. Chem. 281, 19925–19938. doi: 10.1074/jbc.M600482200
Kuzmin, D., Gogvadze, E., Kholodenko, R., Grzela, D. P., Mityaev, M., Vinogradova, T., et al. (2010). Novel strong tissue specific promoter for gene expression in human germ cells. BMC Biotechnol. 10:58. doi: 10.1186/1472-6750-10-58
Kuzmina, N. B., and Borisov, N. M. (2011). Handling complex rule-based models of mitogenic cell signaling (on the example of ERK Activation upon EGF Stimulation). Int. Proc. Chem. Biol. Environ. Eng. 5, 67–82.
Mityaev, M. V., Kopantzev, E. P., Buzdin, A. A., Vinogradova, T. V., and Sverdlov, E. D. (2010). Enhancer element potentially involved in human survivin gene promoter regulation in lung cancer cell lines. Biochemistry (Mosc.) 75, 182–191. doi: 10.1134/S0006297910020082
Moskalev, A. A., and Shaposhnikov, M. V. (2010). Pharmacological inhibition of phosphoinositide 3 and TOR kinases improves survival of Drosophila melanogaster. Rejuvenation Res. 13, 246–247. doi: 10.1089/rej.2009.0903
Murphy, G. G., Shah, V., Hell, J. W., and Silva, A. J. (2006). Investigation of age-related cognitive decline using mice as a model system: neurophysiological correlates. Am. J. Geriatr. Psychiatry 14, 1012–1021. doi: 10.1097/01.JGP.0000209404.54310.b3
Park, S. K., Kim, K., Page, G. P., Allison, D. B., Weindruch, R., and Prolla, T. A. (2009). Gene expression profiling of aging in multiple mouse strains: identification of aging biomarkers and impact of dietary antioxidants. Aging Cell 8, 484–495. doi: 10.1111/j.1474-9726.2009.00496.x
Savoie, C. J., Aburatani, S., Watanabe, S., Eguchi, Y., Muta, S., Imoto, S., et al. (2003). Use of gene networks from full genome microarray libraries to identify functionally relevant drug-affected genes and gene regulation cascades. DNA Res. 10, 19–25. doi: 10.1093/dnares/10.1.19
Tacutu, R., Craig, T., Budovsky, A., Wuttke, D., Lehmann, G., Taranukha, D., et al. (2013). Human ageing genomic resources: integrated databases and tools for the biology and genetics of ageing. Nucleic Acids Res. 41, D1027–D1033. doi: 10.1093/nar/gks1155
Trindade, L. S., Aigaki, T., Peixoto, A. A., Balduino, A., Manica Da Cruz, I. B., and Heddle, J. G. (2013). A novel classification system for evolutionary aging theories. Front. Genet. 4:25. doi: 10.3389/fgene.2013.00025
Welle, S., Brooks, A. I., Delehanty, J. M., Needler, N., and Thornton, C. A. (2003). Gene expression profile of aging in human muscle. Physiol. Genomics 14, 149–159. doi: 10.1152/physiolgenomics.00049.2003
Ye, X., Linton, J. M., Schork, N. J., Buck, L. B., and Petrascheck, M. (2013). A pharmacological network for lifespan extension in Caenorhabditis elegans. Aging Cell. doi: 10.1111/acel.12163. [Epub ahead of print].
Yizhak, K., Gabay, O., Cohen, H., and Ruppin, E. (2013). Model-based identification of drug targets that revert disrupted metabolism and its application to ageing. Nat. Commun. 4:2632. doi: 10.1038/ncomms3632.
Zabolotneva, A. A., Bantysh, O., Suntsova, M. V., Efimova, N., Malakhova, G. V., Schumann, G. G., et al. (2012a). Transcriptional regulation of human-specific SVAF(1) retrotransposons by cis-regulatory MAST2 sequences. Gene 505, 128–136. doi: 10.1016/j.gene.2012.05.016
Zabolotneva, A. A., Zhavoronkov, A., Garazha, A. V., Roumiantsev, S. A., and Buzdin, A. A. (2012b). Characteristic patterns of microRNA expression in human bladder cancer. Front. Genet. 3:310. doi: 10.3389/fgene.2012.00310
Zahn, J. M., Poosala, S., Owen, A. B., Ingram, D. K., Lustig, A., Carter, A., et al. (2007). AGEMAP: a gene expression database for aging in mice. PLoS Genet. 3:e201. doi: 10.1371/journal.pgen.0030201
Zhavoronkov, A., and Litovchenko, M. (2013). Biomedical progress rates as new parameters for models of economic growth in developed countries. Int. J. Environ. Res. Public Health 10, 5936–5952. doi: 10.3390/ijerph10115936
Zhavoronkov, A., Smit-Mcbride, Z., Guinan, K. J., Litovchenko, M., and Moskalev, A. (2012). Potential therapeutic approaches for modulating expression and accumulation of defective lamin A in laminopathies and age-related diseases. J. Mol. Med. (Berl.) 90, 1361–1389. doi: 10.1007/s00109-012-0962-4
Keywords: geroprotector, aging-suppressive drug, signaling pathway cloud, validation of drugs, transcriptomics
Citation: Zhavoronkov A, Buzdin AA, Garazha AV, Borissov NM and Moskalev AA (2014) Signaling pathway cloud regulation for in silico screening and ranking of the potential geroprotective drugs. Front. Genet. 5:49. doi: 10.3389/fgene.2014.00049
Received: 09 December 2013; Accepted: 16 February 2014;
Published online: 03 March 2014.
Edited by:Elena G. Pasyukova, Institute of Molecular Genetics of Russian Academy of Sciences, Russia
Reviewed by:Hao Mei, University of Mississippi Medical Center, USA
Taru Deva, Biosymfonix, Denmark
Douglas Gray, Ottawa Hospital Research Institute, Canada
Copyright © 2014 Zhavoronkov, Buzdin, Garazha, Borissov and Moskalev. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alex Zhavoronkov, The Biogerontology Research Foundation, 4 Hill Street, London W1J 5NE, UK e-mail: firstname.lastname@example.org;
Alexey A. Moskalev, Komi Science Center of Russian Academy of Sciences, Institute of Biology, Kommunisticheskaja St. 28, Syktyvkar 167982, Russia e-mail: email@example.com