Open Access

Systems pharmacology and genome medicine: a future perspective

Genome Medicine20091:11

DOI: 10.1186/gm11

Published: 22 January 2009


Genome medicine uses genomic information in the diagnosis of disease and in prescribing treatment. This transdisciplinary field brings together knowledge on the relationships between genetics, pathophysiology and pharmacology. Systems pharmacology aims to understand the actions and adverse effects of drugs by considering targets in the context of the biological networks in which they exist. Genome medicine forms the base on which systems pharmacology can develop. Experimental and computational approaches enable systems pharmacology to obtain holistic, mechanistic information on disease networks and drug responses, and to identify new drug targets and specific drug combinations. Network analyses of interactions involved in pathophysiology and drug response across various scales of organization, from molecular to organismal, will allow the integration of the systems-level understanding of drug action with genome medicine. The interface of the two fields will enable drug discovery for personalized medicine. Here we provide a perspective on the questions and approaches that drive the development of these new interrelated fields.


Our knowledge of the mechanisms by which drugs act physiologically advanced radically during the twentieth century. With the advent of biochemistry and molecular biology, the targets of drugs became increasingly well characterized. The development of receptor theory by Clark [1] and Black [2, 3], followed by analyses that distinguished between competitive and non-competitive inhibition, began to shed light on the mechanisms by which drugs worked at the molecular level [4]. The influence and relevance of receptor theory in modern pharmacology is derived from the large number of drugs that target membrane receptors, the majority of which are G protein-coupled receptors (GPCRs). The theory of enzyme kinetics led to substrate-based inhibitor design of drugs. These theoretical underpinnings, the size of the market for specific classes of drugs and the ease of drug design for a proven target have resulted in many similar drugs that can target a single protein. ACE inhibitors that are used to treat hypertension are good examples of this approach. The drug pipeline has evolved, with the appearance of targeted therapies and biological therapeutics, such as monoclonal antibody therapies. Many diseases, such as hypertension, ulcers and several types of cancer, that could not be treated two generations ago, can now successfully be managed, if not cured. Yet the 'drugome' (the proteins and genes that are targeted by drugs approved by national regulators such as the US Food and Drug Administration, FDA) covers only a small fraction of the proteome or the 'diseaseome' (genes that have been linked with disease), and many drugs are focused in just a few areas (Figure 1) [5, 6]. This disparity reflects the current relationship between basic biological science and its use for therapeutic purposes. There are substantial opportunities to use the accumulated knowledge of biological processes for drug discovery and clinical applications. If we are to take advantage of such opportunities, genome medicine and systems pharmacology need to be well integrated.
Figure 1

Relationships between the genome, proteome, diseaseome and drugome. The number of distinct protein species (about 400,000) comprising the proteome (green circle, scaled down by 25% relative to the other circles), is estimated by taking the approximately 25,000 currently annotated genes (yellow circle) and assuming about four splice variants per gene and about four post-translationally modified proteins per splice variant. The genome, diseaseome and drugome form a Venn diagram. The red circle represents the approximately 1,800 genes known to be involved in various diseases (the diseaseome). Of these, a small fraction (the drugome) is targeted by FDA-approved drugs. Not all drug targets have been characterized as disease genes. In total, proteins encoded by approximately 400 genes (0.1% of the proteome) are targeted by about 1,200 FDA-approved drugs. There are more drugs than protein targets because more than one drug can target the same protein.

As the systems-level understanding of biological processes expands, it is becoming a crucial driver of pharmacology that is anchored in the human genome and personalized medicine. The path from laboratory research to clinical application is becoming short as translational research grows, facilitating collaborations between basic researchers and clinicians. Genomic and proteomic technologies drive discovery of biomarker sets for the classification of diseases and the stages of their progression, as exemplified by microarray-based marker sets that have been developed to identify stages of cancer progression [7, 8]. Although more of these approaches need to be discovered and then standardized before they are routinely used in clinical practice, the importance of using systems-type methodologies to characterize therapeutic interventions, to delineate the pathways (or more often networks) involved in disease, and to identify the mechanisms of action and off-target effects of current drugs is becoming clearer. A multi-faceted understanding of therapeutic intervention is necessary, given the complexity of human physiology and the increasing availability of numerous clinical parameters and analyses.

Here, we explain the reasoning underlying the assertion that systems-level knowledge of pharmacology and pathophysiology, rooted in genomic information, will increase the efficacy of existing drugs by aiding in the development of personalized medicine and will facilitate the rational discovery of new drugs using a much wider target base. For pharmacology to be understood at a systems level, it is necessary to use a genome-based approach to systems-level studies of physiology and pathophysiology.

Genome medicine

An operational definition of genome medicine is the way in which genomic information from a patient helps in the diagnosis or treatment of disease. This includes areas such as the genetics-based diagnosis of the origin of diseases, their progression and the response to drugs. In the area of drug response, genome medicine overlaps with pharmacogenomics, a field that studies how genome variation affects drug response. It is well established that some diseases arise from mutations in single genes and can be treated using the wild-type product from that particular gene. Fabry's disease is an example of such a monogenic disease [9, 10]. However, it is now known that many complex diseases arise from interactions and changes in multiple genes. The two-hit model for cancer [11] was an early example of this recognition. Mapping of genetic variations, such as single-nucleotide polymorphisms (SNPs), in the human genome has given rise to the idea that combinations of variant genes can alter susceptibility to various diseases. Genome-wide association studies have become popular over the past few years [12], although their ability to predict disease susceptibility is only beginning to be determined. Nevertheless, it is clear there is sufficient genomic variation between individuals to affect the origin and progression of diseases, as well as drug response. For example, genetic testing for cancer drugs, such as imatinib [13, 14], trastuzumab [15], gefitinib [16] and bucindolol [17], provides information about whether the form of the protein target that is expressed in each individual will be responsive to the drug. Testing for the gene BRCA1 is used to determine the feasibility of preventative measures for breast cancer [18]. Other types of genetic testing are more pharmacogenomic in nature, such as those for warfarin [19, 20] or tamoxifen [21]. These tests screen for polymorphisms in distinct cytochrome P450 isoforms, which are responsible for differences in how the drug is metabolized and thus affect the therapeutic index, the ratio between the toxicity and effectiveness of the drug. Pharmacogenomics can thus be considered to be a very important part of genome medicine. Recent views on pharmacogenomics describe in greater depth the relationship between genomics and drug response [22, 23].

Our understanding of a patient's genome will increasingly drive medical practice in the 21st century. The development of technologies for rapid sequencing of whole genomes is an important factor in this changing approach to medical practice. The growing number of treatment options for pathophysiology, and the accumulation of knowledge about the risk-to-benefit ratio of prescribing one therapy over another when anchored in individual genomic information, should allow the individualized tailoring of therapies. Such tailoring would be based on an individual definition of the therapeutic index of a drug for a specific individual. Currently, many clinical practices are based on empirical trial-and-error approaches, and drug usage typically follows a 'one size fits all' approach.

Even for diseases for which therapies have been very successful, we generally do not understand why one type of therapy works for one individual and not for another. For example, there are four general types of therapy to treat hypertension (thiazide-based diuretics, angiotension converting enzyme inhibitors or angiotension II receptor blockers, calcium channel blockers and β-adrenergic receptor blockers or beta-blockers) [24]. There is relatively little predictability as to which hypertension treatment will be more effective for any given patient, and there is a large patient-to-patient variability in response to each therapy and required dosage. Current data suggest that 50% of patients who do not respond to one type of therapy will respond to another and that 70-80% of patients will respond if switched for a second time to yet another type of antihypertensive drug [25]. The characterization of all of the hypertension-related genes in the patient's genome may facilitate the construction and systems-level analyses of the regulatory networks, using those hypertension-related genes as seed nodes. From the physiological functions of these networks, we may be able to identify the pathways involved in disease of an individual to predict the action of a drug. Such personalized disease networks for patients would allow clinical practice to move away from a trial-and-error approach to prescribing drugs and to advance to a more genome-informed disease management. This type of clinical practice could easily be called 'genome medicine'.

Although the state of the genome will be a major determinant of how diseases originate and progress, it is unlikely to be the only cause in many cases. Environmental factors also have an important role. In many cases these two factors are related, whereby a particular genotype changes the risk of an environmentally induced disease. This relationship can be understood by considering the biochemistry and physiology of the system. Genes are most often transcribed and translated into components of biochemical networks that underlie cellular processes. Often there is more than one type of change in a particular gene; some of these types of change lead to a change in the function of a key cellular component and ultimately to a disease state, whereas others are inconsequential. Sometimes, even if a mutant gene product behaves aberrantly at a biochemical level, this behavior can be compensated for and the physiology remains normal. Alternatively, the aberrant behavior of diseased cells might not be only genetic in origin but also a result of changes in normal cell signaling processes that regulate transcription, translation and effector protein function. The hypothesis that inflammation underlies the origins of several diseases, for example, is based on the assumption that the normal process has gone awry because of environmental signals, resulting in sustained activation of intracellular signaling networks that results in pathophysiologies. There is accumulating evidence to support this hypothesis [26, 27].

The effects of post-transcriptional and translational regulation of signaling can increase the difficulty of developing methods for individualized mechanism-based therapeutic interventions. A systems understanding of disease and drug action at the level of cellular biochemical reactions and physiological function would be useful in framing the effect of genomic variations in the context of environmental cues. This is where systems pharmacology comes in.

Systems pharmacology

Systems pharmacology seeks to develop a global understanding of the interactions between pathophysiology and drug action. To develop such an understanding it is necessary to analyze interactions across and between various scales of organization. The representations of different scales are illustrated in Figure 2. The biological insights gained from multi-scale analyses of physiological processes have been noted previously [28, 29]. Analysis of such multi-scale systems requires one to 'zoom' in and out depending on the type of analysis being conducted. During the process of multi-scale analysis (zooming), it is essential that we develop a mechanistic understanding of the relationships across the various levels. Simply correlating structural information or molecular interactions with clinical phenotypes is a good starting point, but it will not yield the ability to predict disease progression or drug treatment outcomes. The type of 'system' to be analyzed can vary depending on the zoom level (Figure 2) of information desired. The system can be considered at the organismal, organ, tissue, cellular or molecular levels. The effects of a drug on pathophysiology that are seen at the organismal level, that is, symptoms or clinical measurements, are zoomed-out observations. These observations usually consist of clinical data, ranging from blood chemistry to measurements reflecting organismal function, such as blood pressure and stress tests, all of which are documented in the electronic medical records of the patient that will aid future computational analyses across scales.
Figure 2

Multi-scale analyses in systems pharmacology. The top half of the figure is a schematic representation of different scales of organization involved in human pathophysiology and systems pharmacology. Clinical indicators and analyses (left) indicate measurements of various types of blood concentrations, blood pressure, stress and so on; these parameters are available in the electronic medical records of patients. From left to right, the scale becomes smaller, or 'zoomed in'. The human body (or organism) can be analyzed at the levels of organs, tissues, cells (represented here together with tissues) or molecules. Drugs are prescribed and taken at the organismal level but exert their effects by interacting with their target at the molecular level (red arrow). The gradient from white to blue corresponds to the various levels of interaction systems: white represents a clinical setting; blue represents a laboratory setting. Studies in systems pharmacology fully span all levels shown here.

The physiological responses to a drug or disease are manifestations of events that can be studied at zoomed-in perspectives. Emerging imaging technologies are complementing pathology to observe disease progression and treatment outcome. The entire network of events that are involved in therapeutic intervention can be analyzed in proteomic or genomic studies, or, at a more detailed level, specific protein-protein interactions or enzymatic reactions can be examined. Much of the zoomed-in systems pharmacology understanding comes from basic laboratory research. The methodologies available to address the questions across zoom levels include biochemical experiments, microarray and other high-throughput methodologies, animal models and computational modeling and analysis of molecular and cellular functions. Generally, a systems approach that integrates knowledge from analyses across multiple zoom levels is likely to uncover a new target in a pathway of interest, or the basis for an observed adverse effect, that would probably not have been discovered using traditional techniques. In order for a discovery made at the molecular or cellular level to be applicable clinically, the effect must be demonstrated at the organ and organismal levels and established to be applicable to a population. Thus, systems pharmacology uses a range of approaches that spans multiple scales of organization, as is shown in Figure 2. The different areas of study in systems pharmacology and their relationship to genome medicine are briefly discussed below.

Research in systems pharmacology integrated into genome medicine tends to follow two directions: the 'physiology to therapeutic intervention' direction (bottom-up) comprises approaches that facilitate new drug or drug-target discovery; and the 'therapeutic intervention to physiology' direction (top-down) focuses on the characterization of current drugs across scales of interactions in terms of their mechanism of action, off-target effects or similarity to other drugs. Studies in the bottom-up direction include network analysis, which provides a foundation for systems pathophysiology. Network analysis can be used to identify new targets for therapeutic intervention and to understand adverse events and drug resistance from genomic information. Bottom-up analyses are based on the integration between systems pharmacology and genome medicine (Figure 3). The bottom-up approach can be complemented by top-down studies and include analyses of networks at various levels to study the effect of a drug on a system, commonalities among drugs (global analyses of the 'drugome'), drug resistance and effects of drug combinations.
Figure 3

The relationship between genome medicine and systems pharmacology. The diagram summarizes various aspects of genome medicine (in blue) and systems pharmacology (in yellow). Overlapping aspects of analyses and practice are in green (intersection of circles). The positioning of the circles indicates the operational classification of 'genome medicine to systems pharmacology' as top-down and 'systems pharmacology to genome medicine' as bottom-up. The key analyses and practices are in the circle for the field that uses them. Approaches and practices that are used in both fields are in the overlapping region. Genome medicine starts with genetic and genomic testing. Experimental data are computationally processed using statistical genetics tools to yield information that is used in personalized medicine for therapeutic-index targeting (such as dosage of warfarin) and combination therapy. Network analysis is a common approach that integrates genome medicine and systems pharmacology. Systems pharmacology starts from cataloguing the characteristics of individual drugs and targets from biochemistry and cell-physiology experiments. Computational methods and genomic and proteomic data together enable the use of this catalog of information to make predictions regarding drug discovery, drug action and adverse events. Such predictions can be experimentally and clinically tested. Approaches common to both genome medicine and systems pharmacology are based on network analyses that underlie systems pathophysiology, whereby the origins of disease are understood in the context of multi-scale systems. Such understanding enables network-based drug screening and whole genome-based predictions of adverse events and drug resistance. Thus, ultimately, therapeutics intervention will be guided by integrating genome medicine and systems pharmacology.

The role of systems pathophysiology in bottom-up intervention

Many efforts to systematically understand the cellular processes involved in disease lead to sets of proteins that are mutated to become under- or overactive, deleted, over-expressed or aberrantly post-translationally modified in the disease state. Genomic and proteomic studies can be designed to globally screen for cellular components that stand out in the disease model. Subsequent small-scale studies often focus on verifying one or several of the predictions from the large-scale screen. In order to obtain a multi-faceted analysis, techniques such as RNA interference (RNAi) screens, gene-expression profiling, chromatin immunoprecipitation (ChIP) analyses or protein micro-arrays are combined. Analyses on patient-derived samples allow data from humans to constrain and validate animal or cell models of disease. Approaches such as these have been applied to cancer [3032], malaria [33], heart failure [34] and HIV [35, 36].

These global analyses to identify genes, proteins and other cellular components related to the origin or progression of a specific disease yield lists of 'seed nodes' that can be used to computationally construct disease networks by integrating knowledge about protein interactions reported in biomedical literature. With additional experiments or analyses, it becomes possible to place seemingly unrelated proteins in the context of a pathophysiology and thus to explore whether they are potential therapeutic targets or biomarkers. The use of a network of known pathways surrounding genes with annotated cardiovascular involvement, for example, enabled the identification of a series of mass spectrometry biomarkers for major adverse cardiac events [37]. Similarly, a network-based approach was able to identify sets of interactions between nuclear receptors and insulin signaling-pathway proteins related to type 2 diabetes [38]. Identification of previously unrelated relationships between disease pathways and cellular functions is a first step in identifying drug targets.

The accumulation of published systems-level datasets related to drug action calls for new ways to compile the data into a computable format. An example of this type of effort is the Connectivity Map [39], which finds correlations between gene-expression signatures and the sets of proteins involved with the action of a class of drugs or with a particular disease. Although the correlations found through the Connectivity Map approach are an excellent starting point to gain insights into how the gene-expression signatures indicate aberrant cellular or physiological behavior, more detailed mechanistic studies are necessary to develop predictive capabilities.

Even at the level of a single protein, a systems-level understanding becomes useful. A mutant protein implicated in disease can have multiple changes in its behavior in the network. These changes must be put into the context of the disease 'system' in order to determine which property has the highest impact. Different mutations in the small G protein Ras, which are found in many cancers, increase Ras activity by several different mechanisms [40]. To understand which mutation has a more important role in oncogenesis, the effect of two mechanistically distinct Ras mutations on the entire downstream signaling network was analyzed computationally and the results were confirmed experimentally [41]. The computational model predicted that a drug that specifically binds GTP-bound Ras will have a more specific effect on the mutant than it will on the wild-type system, and thus also on cancerous than on normal tissue. Such a prediction could not have been made by simply correlating mutations with disease phenotype. Thus, a systems understanding of a single protein function in a disease signaling network can predict new targets for therapeutic intervention in a mechanism-based manner.

Top down approaches: global drug analyses

Direct relationships between diseases, drugs and proteins enable statistical inferences about less obvious relationships between them. Global analysis of FDA-approved drugs has found that, unlike ubiquitously expressed essential genes, drug targets tend to be expressed in specific tissues yet tend to interact with many other proteins in the cellular network while being independently regulated [5, 6]. A bipartite graph of 1,052 FDA-approved drugs interacting with 485 targets contained 179 'islands' (sets of nodes that are not connected to any other node in the graph). Most of these islands are made up of 10-30 interacting cellular components. A single large island of 481 components consisted of drugs that target GPCRs [5]. This large island exists in part because of the many physiological processes, as varied as cardiac contractility, acid secretion and airway constriction, that are regulated by GPCRs, and in part because the extracellular ligand-binding domains of these receptors make ideal drug targets.

Many drugs target GPCRs and, despite the physiological diversity of these proteins, all GPCR-mediated signaling is coupled through heterotrimeric G proteins. Commonalities of signaling mechanism among existing drug targets, such as this, create a basis for identifying what network properties make a protein a potentially good drug target. Statistical metrics from network analyses, such as a centrality measure, which quantifies the relative importance of a protein in communicating between different modules within a network, have been suggested for identifying nodes (proteins in a network) that have attractive properties as potential drug targets [42]. Recent work has generated a network connecting drugs on the basis of their structural similarity and similarity of side-effect profiles [43]. This method has been proven effective at identifying groups of drugs that share common targets. When drugs with predicted shared targets did not in fact have known targets in common, the predictions provided a framework for testing these drugs for binding against each other's targets. Hence, this approach led to the identification of new targets for existing drugs.

Integrating the bottom-up and top-down approaches: network analysis

Studies on cellular signaling pathways and networks drive systems pharmacology because they can lead to discoveries that enable new therapeutic interventions. These networks are often based on genomic information. Network analysis integrates genome medicine and systems pharmacology, as is shown by its central location in the overlapping part of the Venn diagram in Figure 3. Proteomic and genomic studies have increased dramatically in the past five years, and datasets describing global behavior, such as genomic or proteomic analyses of many cellular pathways and processes, are now available [44]. These data expand our knowledge of signaling pathways by implicating many more signaling molecules, which may serve as new drug targets or have implications for drug resistance or off-target effects.

Certain proteins or protein families are targeted more frequently in therapeutic interventions than others, because of their involvement in disease and the straightforwardness of designing drugs to target the proteins [45]. Systems pharmacology studies on these systems are usually more common. For example, there are at least eight different approved or pipeline drugs that target one or several members of the ErbB family of receptor tyrosine kinases, which are modified in many kinds of cancer and other diseases [46]. There are also many studies that computationally model [4752] or experimentally analyze [51, 5356] the dynamics of the ErbB signaling network. For example, the ErbB network has been analyzed by examining the interactions between the known phosphotyrosine interaction domains (PTB and SH2 domains) and phosphopeptides representing the phosphorylated forms of the four ErbB family members using microarrays [57]. This study [57] uncovered behavior of the ErbB network that is pertinent to its role in oncogenesis and the implications of therapeutically modulating this family of receptors.

Network analysis has facilitated the development of computational prediction algorithms that predict all possible molecules affected by specified perturbations of upstream targets by hypothetical drugs [58]. These algorithms can also solve the 'minimum knockout' problem by yielding predictions of the smallest number of drug targets needed to fully block a cellular process [58]. A similar computational framework for identifying targeted disruption of signaling networks allows identification of all sets of potential drug targets to block a process, while not affecting other processes [59]. These types of network analysis suggest that computational analyses of cell signaling networks are likely to be an important aspect of systems pharmacology and future drug discovery [60].

Adverse events and drug resistance: relationship to genome characteristics and systems pathophysiology

Most studies that focus on understanding resistance to therapies are in oncology. Many cancer therapies, such as microtubule stabilizers [61], tamoxifen or endocrine therapy [62] and drugs targeted at epithelial growth factor receptors (EGFRs) [63], are effective for a limited period and/or in a limited population and then resistance develops. Screens using RNAi or with genomics and proteomics identify proteins that are up- or down-regulated or are necessary for drug action in resistant cell populations. The receptor tyrosine kinase c-MET has recently been implicated in resistance to EGFR-targeted therapies and, on the basis of these discoveries, therapies directed against c-MET are in development [64]. In order to understand better the convergence of signals from these two receptors, the tyrosine signaling networks of several cancer cell lines that overexpress c-MET or EGFR or express a mutant form of EGFR were examined to identify a core network of 50 proteins mediating drug response [65]. This type of network analysis can form the basis for the selection of drugs that target the proteins common to both pathways and thus overcome drug-induced resistance.

Genome medicine and systems pharmacology need to be integrated for defining the genes and proteins involved in drug treatment or drug resistance. Such integrated analyses can lead to identification of targets that are likely to synergize with or add to the effect of the drug or therapy. For example, studies focused on defining the genes involved in radiation therapy aim to discover new targets that would increase the effectiveness of radiation therapy [66, 67]. In fact, many drugs are prescribed in combination with other drugs, because the effectiveness of both drugs is increased when they are combined. Predicting which combinations will show this effect, and the doses of such combinations, is not entirely intuitive or simple. Sometimes, serendipitous results from genomic or proteomic screens identify targets that might have synergistic effects in conjunction with commonly used therapies, such as ceramide transport protein in taxane-based therapy [68]. In proactive approaches, network analysis [42] or computational modeling [69] can provide information on the effect of intracellular signaling on inhibition of the two nodes, predicting how dosing with pairs or groups of drugs that target different proteins would work. Such a model was used to predict combination therapies in the EGFR pathway [70]. These types of study can also be used to look for ways to lower dosages while sustaining effectiveness, avoiding unnecessary drug toxicity.

Adverse events from drugs are a major concern in the development and prescription of pharmaceuticals. The susceptibility and severity of adverse events can vary tremendously between patients as a result of, among other confounders, genomic factors. At the simplest level, adverse events can be due to dosage effects of the drug. Pharmacogenomic factors affecting drug metabolism can lead to increased levels of active drug in the body. For drugs such as warfarin, which acts as an anticoagulant, too high a blood level owing to reduced metabolism (resulting from variant CYPs) can lead to uncontrolled bleeding and cerebral hemorrhage [20]. For patients with CYPs that have lower warfarin metabolizing capabilities, the dosage of warfarin needs to be reduced so as to obtain its therapeutic benefits without increasing the risk of adverse events.

Genetic factors can affect how patients respond to therapeutic doses of a drug. Patients with a deficiency in glucose 6-phosphatase dehydrogenase (G6PD) cannot respond sufficiently to multiple forms of oxidative stress and, thus, many drugs, including many anti-malarials, analgesics and antibiotics, can cause these patients to develop hemolytic anemia [71]. Adverse events can also be modulated by the immune system, whereby the patient's immune system responds to the drug, leading to damage of otherwise healthy cells. For example, this can lead to drug-induced neutropenia [72]. Several of this class of adverse event have been associated genetically with specific HLA alleles, the genes that control presentation of immunogenic epitopes [73].

Drug-induced cardiac arrhythmias comprise another important class of side-effects. One example, the long QT syndrome, has been extensively studied because it can in many cases lead to fatal arrhythmias. Increased understanding of a congenital form of the syndrome, the frequency of the side-effect and the variability in its severity has led researchers to suggest that individual variability in cardiac repolarization reserve is among the major risk modifiers [74]. Cardiac repolarization reserve is the capability of cells to compensate for changes in ion channel function that underlie the normal myocyte action potential. Thus, understanding how drugs regulate specific ion channels to modulate electrical interactions between channels to produce a physiological event is critical. This type of adverse event, which reflects an interaction between a drug and a complex biological system, is where systems pharmacology will prove most effective. Systems pharmacology projects have begun to identify targets that are important for causing the side-effect by grouping drugs according to their side-effect profiles [43]. In addition, large-scale studies have begun to identify cellular modules that are important for causing the side-effects. For example, a large-scale study of mitochondrial function identified mitochondrial expression profiles affected by drugs that can cause drug-induced myopathies [75].


Both genome medicine and systems pharmacology, as interdisciplinary research fields, are in their infancy. The convergent goals of both fields are to treat disease in each patient on the basis of the patient's genome and their unique environmental interactions. Thus, they can be considered to be the twin pillars supporting the gateway to personalized medicine. In both fields, there are considerable opportunities for new conceptual and technological developments. Success in personalized medicine will require advances in technology and concepts in both fields to be well integrated. The most influential driver for this integration is a mechanism-based understanding of disease and therapy across scales of biological organization. As a multi-scale understanding of the human body develops, it is likely to influence not only the treatment and prevention of disease, but also the economics and social aspects of health care.



angiotensin converting enzyme


breast cancer 1, early onset (mutations in this gene greatly increase predisposition to breast cancer)


chromatin immunoprecipitation


cytochrome P450 (drug metabolizing enzymes)


epithelial growth factor receptor


US Food and Drug Administration


glucose 6-phosphatase dehydrogenase


G-protein coupled receptor


RNA interference


single nucleotide polymorphism.



This work was supported by the NIH grants GM54508 and New York Systems Biology Center grant P50-GM071558. SB is supported by pre-doctoral training grant in Pharmacological Sciences GM-062754. We thank Avi Ma'ayan and Emmanuel Landau for comments.

Authors’ Affiliations

Department of Pharmacology and Systems Therapeutics, Mount Sinai School of Medicine


  1. Kenakin T: Principles: receptor theory in pharmacology. Trends Pharmacol Sci. 2004, 25: 186-192. 10.1016/ Article
  2. Black J, Leff P: Operational models of pharmacological agonist. Proc R Soc Lond B Biol Sci. 1983, 220: 141-162. 10.1098/rspb.1983.0093.PubMedView Article
  3. Maehle AH, Prull CR, Halliwell RF: The emergence of the drug receptor theory. Nat Rev Drug Discov. 2002, 1: 637-641. 10.1038/nrd875.PubMedView Article
  4. Colquhoun D: The quantitative analysis of drug-receptor interactions: a short history. Trends Pharmacol Sci. 2006, 27: 149-157. 10.1016/ Article
  5. Ma'ayan A, Jenkins SL, Goldfarb J, Iyengar R: Network analysis of FDA approved drugs and their targets. Mt Sinai J Med. 2007, 74: 27-32. 10.1002/msj.20002.PubMedPubMed CentralView Article
  6. Yildirim MA, Goh KI, Cusick ME, Barabasi AL, Vidal M: Drug-target network. Nat Biotechnol. 2007, 25: 1119-1126. 10.1038/nbt1338.PubMedView Article
  7. Mook S, Schmidt MK, Viale G, Pruneri G, Eekhout I, Floore A, Glas AM, Bogaerts J, Cardoso F, Piccart-Gebhart MJ, Rutgers ET, Van't Veer LJ, on behalf of the TRANSBIG consortium: The 70-gene prognosis-signature predicts disease outcome in breast cancer patients with 1-3 positive lymph nodes in an independent validation study. Breast Cancer Res Treat. 2008, doi: 101007/s10549-008-0130-2.
  8. Kroon BK, Leijte JA, van Boven H, Wessels LF, Velds A, Horenblas S, van't Veer LJ: Microarray gene-expression profiling to predict lymph node metastasis in penile carcinoma. BJU Int. 2008, 102: 510-515. 10.1111/j.1464-410X.2008.07697.x.PubMedView Article
  9. Desnick RJ: Enzyme replacement therapy for Fabry disease: lessons from two alpha-galactosidase A orphan products and one FDA approval. Expert Opin Biol Ther. 2004, 4: 1167-1176. 10.1517/14712598.4.7.1167.PubMedView Article
  10. Desnick RJ, Schuchman EH: Enzyme replacement and enhancement therapies: lessons from lysosomal disorders. Nat Rev Genet. 2002, 3: 954-966. 10.1038/nrg963.PubMedView Article
  11. Land H, Parada LF, Weinberg RA: Tumorigenic conversion of primary embryo fibroblasts requires at least two cooperating oncogenes. Nature. 1983, 304: 596-602. 10.1038/304596a0.PubMedView Article
  12. Manolio TA, Brooks LD, Collins FS: A HapMap harvest of insights into the genetics of common disease. J Clin Invest. 2008, 118: 1590-1605. 10.1172/JCI34772.PubMedPubMed CentralView Article
  13. DeMatteo RP: The GIST of targeted cancer therapy: a tumor (gastrointestinal stromal tumor), a mutated gene (c-kit), and a molecular inhibitor (STI571). Ann Surg Oncol. 2002, 9: 831-839. 10.1007/BF02557518.PubMedView Article
  14. Schiffer CA, Hehlmann R, Larson R: Perspectives on the treatment of chronic phase and advanced phase CML and Philadelphia chromosome positive ALL. Leukemia. 2003, 17: 691-699. 10.1038/sj.leu.2402879.PubMedView Article
  15. Harris L, Fritsche H, Mennel R, Norton L, Ravdin P, Taube S, Somerfield MR, Hayes DF, Bast RC: American Society of Clinical Oncology 2007 update of recommendations for the use of tumor markers in breast cancer. J Clin Oncol. 2007, 25: 5287-5312. 10.1200/JCO.2007.14.2364.PubMedView Article
  16. Sequist LV, Martins RG, Spigel D, Grunberg SM, Spira A, Jänne PA, Joshi VA, McCollum D, Evans TL, Muzikansky A, Kuhlmann GL, Han M, Goldberg JS, Settleman J, Iafrate AJ, Engelman JA, Haber DA, Johnson BE, Lynch TJ: First-line gefitinib in patients with advanced non-small-cell lung cancer harboring somatic egfr mutations. J Clin Oncol. 2008, 26: 2442-2449. 10.1200/JCO.2007.14.8494.PubMedView Article
  17. Liggett SB, Mialet-Perez J, Thaneemit-Chen S, Weber SA, Greene SM, Hodne D, Nelson B, Morrison J, Domanski MJ, Wagoner LE, Abraham WT, Anderson JL, Carlquist JF, Krause-Steinrauf HJ, Lazzeroni LC, Port JD, Lavori PW, Bristow MR: A polymorphism within a conserved β1-adrenergic receptor motif alters cardiac function and β-blocker response in human heart failure. Proc Natl Acad Sci USA. 2006, 103: 11288-11293. 10.1073/pnas.0509937103.PubMedPubMed CentralView Article
  18. Nusbaum R, Isaacs C: Management updates for women with a BRCA1 or BRCA2 mutation. Mol Diagn Ther. 2007, 11: 133-144.PubMedView Article
  19. Schwarz UI, Ritchie MD, Bradford Y, Li C, Dudek SM, Frye-Anderson A, Kim RB, Roden DM, Stein CM: Genetic determinants of response to warfarin during initial anticoagulation. N Engl J Med. 2008, 358: 999-1008. 10.1056/NEJMoa0708078.PubMedPubMed CentralView Article
  20. Au N, Rettie AE: Pharmacogenomics of 4-hydroxycoumarin anticoagulants. Drug Metab Rev. 2008, 40: 355-375. 10.1080/03602530801952187.PubMedView Article
  21. Goetz MP, Kamal A, Ames MM: Tamoxifen pharmacogenomics: the role of CYP2D6 as a predictor of drug response. Clin Pharmacol Ther. 2008, 83: 160-166. 10.1038/sj.clpt.6100367.PubMedPubMed CentralView Article
  22. Evans WE, Relling MV: Moving towards individualized medicine with pharmacogenomics. Nature. 2004, 429: 464-468. 10.1038/nature02626.PubMedView Article
  23. Court MH: A pharmacogenomics primer. J Clin Pharmacol. 2007, 47: 1087-1103. 10.1177/0091270007303768.PubMedView Article
  24. Frank J: Managing hypertension using combination therapy. Am Fam Physician. 2008, 77: 1279-1286.PubMed
  25. Materson BJ, Reda DJ, Preston RA, Cushman WC, Massie BM, Freis ED, Kochar MS, Hamburger RJ, Fye C, Lakshman R, et al: Response to a second single antihypertensive agent used as monotherapy for hypertension after failure of the initial drug. Department of Veterans Affairs Cooperative Study Group on Antihypertensive Agents. Arch Intern Med. 1995, 155: 1757-1762. 10.1001/archinte.155.16.1757.PubMedView Article
  26. Libby P: Inflammatory mechanisms: the molecular basis of inflammation and disease. Nutr Rev. 2007, 65: S140-146. 10.1111/j.1753-4887.2007.tb00352.x.PubMedView Article
  27. Lucas SM, Rothwell NJ, Gibson RM: The role of inflammation in CNS injury and disease. Br J Pharmacol. 2006, 147 (Suppl 1): S232-S240. 10.1038/sj.bjp.0706400.PubMedPubMed Central
  28. Hunter PJ, Borg TK: Integration from proteins to organs: the Physiome Project. Nat Rev Mol Cell Biol. 2003, 4: 237-243. 10.1038/nrm1054.PubMedView Article
  29. Hunter PJ, Crampin EJ, Nielsen PM: Bioinformatics, multiscale modeling and the IUPS Physiome Project. Brief Bioinform. 2008, 9: 333-343. 10.1093/bib/bbn024.PubMedView Article
  30. Shaffer AL, Emre NC, Lamy L, Ngo VN, Wright G, Xiao W, Powell J, Dave S, Yu X, Zhao H, Zeng Y, Chen B, Epstein J, Staudt LM: IRF4 addiction in multiple myeloma. Nature. 2008, 454: 226-231. 10.1038/nature07064.PubMedPubMed CentralView Article
  31. Wang Z, Shen D, Parsons DW, Bardelli A, Sager J, Szabo S, Ptak J, Silliman N, Peters BA, Heijden van der MS, Parmigiani G, Yan H, Wang TL, Riggins G, Powell SM, Willson JK, Markowitz S, Kinzler KW, Vogelstein B, Velculescu VE: Mutational analysis of the tyrosine phosphatome in colorectal cancers. Science. 2004, 304: 1164-1166. 10.1126/science.1096096.PubMedView Article
  32. Ruiz-Vela A, Aggarwal M, de la Cueva P, Treda C, Herreros B, Martin-Perez D, Dominguez O, Piris MA: Lentiviral (HIV)-based RNA interference screen in human B-cell receptor regulatory networks reveals MCL1-induced oncogenic pathways. Blood. 2008, 111: 1665-1676. 10.1182/blood-2007-09-110601.PubMedView Article
  33. Tarun AS, Peng X, Dumpit RF, Ogata Y, Silva-Rivera H, Camargo N, Daly TM, Bergman LW, Kappe SHI: A combined transcriptome and proteome survey of malaria parasite liver stages. Proc Natl Acad Sci USA. 2008, 105: 305-310. 10.1073/pnas.0710780104.PubMedPubMed CentralView Article
  34. Frost RJ, Engelhardt S: A secretion trap screen in yeast identifies protease inhibitor 16 as a novel antihypertrophic protein secreted from the heart. Circulation. 2007, 116: 1768-1775. 10.1161/CIRCULATIONAHA.107.696468.PubMedView Article
  35. Brass AL, Dykxhoorn DM, Benita Y, Yan N, Engelman A, Xavier RJ, Lieberman J, Elledge SJ: Identification of host proteins required for HIV infection through a functional genomic screen. Science. 2008, 319: 921-926. 10.1126/science.1152725.PubMedView Article
  36. Lu TC, Wang Z, Feng X, Chuang P, Fang W, Chen Y, Neves S, Maayan A, Xiong H, Liu Y, Iyengar R, Klotman PE, He JC: Retinoic acid utilizes CREB and USF1 in a transcriptional feed-forward loop in order to stimulate MKP1 expression in human immunodeficiency virus-infected podocytes. Mol Cell Biol. 2008, 28: 5785-5794. 10.1128/MCB.00245-08.PubMedPubMed CentralView Article
  37. Jin G, Zhou X, Wang H, Zhao H, Cui K, Zhang XS, Chen L, Hazen SL, Li K, Wong ST: The knowledge-integrated network biomarkers discovery for major adverse cardiac events. J Proteome Res. 2008, 7: 4013-4021. 10.1021/pr8002886.PubMedPubMed CentralView Article
  38. Liu M, Liberzon A, Kong SW, Lai WR, Park PJ, Kohane IS, Kasif S: Network-based analysis of affected biological processes in type 2 diabetes models. PLoS Genet. 2007, 3: e96-10.1371/journal.pgen.0030096.PubMedPubMed CentralView Article
  39. Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN, Reich M, Hieronymus H, Wei G, Armstrong SA, Haggarty SJ, Clemons PA, Wei R, Carr SA, Lander ES, Golub TR: The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006, 313: 1929-1935. 10.1126/science.1132939.PubMedView Article
  40. Repasky GA, Chenette EJ, Der CJ: Renewing the conspiracy theory debate: does Raf function alone to mediate Ras oncogenesis?. Trends Cell Biol. 2004, 14: 639-647. 10.1016/j.tcb.2004.09.014.PubMedView Article
  41. Stites EC, Trampont PC, Ma Z, Ravichandran KS: Network analysis of oncogenic Ras activation in cancer. Science. 2007, 318: 463-467. 10.1126/science.1144642.PubMedView Article
  42. Hwang WC, Zhang A, Ramanathan M: Identification of information flow-modulating drug targets: a novel bridging paradigm for drug discovery. Clin Pharmacol Ther. 2008, doi: 101038/clpt.2008.129.
  43. Campillos M, Kuhn M, Gavin AC, Jensen LJ, Bork P: Drug target identification using side-effect similarity. Science. 2008, 321: 263-266. 10.1126/science.1158140.PubMedView Article
  44. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. 2007, D760-D765. 10.1093/nar/gkl887. 35 Database
  45. Lauss M, Kriegner A, Vierlinger K, Noehammer C: Characterization of the drugged human genome. Pharmacogenomics. 2007, 8: 1063-1073. 10.2217/14622416.8.8.1063.PubMedView Article
  46. Johnston JB, Navaratnam S, Pitz MW, Maniate JM, Wiechec E, Baust H, Gingerich J, Skliris GP, Murphy LC, Los M: Targeting the EGFR pathway for cancer therapy. Curr Med Chem. 2006, 13: 3483-3492. 10.2174/092986706779026174.PubMedView Article
  47. Hendriks BS, Orr G, Wells A, Wiley HS, Lauffenburger DA: Parsing ERK activation reveals quantitatively equivalent contributions from epidermal growth factor receptor and HER2 in human mammary epithelial cells. J Biol Chem. 2005, 280: 6157-6169. 10.1074/jbc.M410491200.PubMedView Article
  48. Kholodenko BN, Demin OV, Moehren G, Hoek JB: Quantification of short term signaling by the epidermal growth factor receptor. J Biol Chem. 1999, 274: 30169-30181. 10.1074/jbc.274.42.30169.PubMedView Article
  49. Nakakuki T, Yumoto N, Naka T, Shirouzu M, Yokoyama S, Hatakeyama M: Topological analysis of MAPK cascade for kinetic ErbB signaling. PLoS ONE. 2008, 3: e1782-10.1371/journal.pone.0001782.PubMedPubMed CentralView Article
  50. Schoeberl B, Eichler-Jonsson C, Gilles ED, Muller G: Computational modeling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors. Nat Biotechnol. 2002, 20: 370-375. 10.1038/nbt0402-370.PubMedView Article
  51. Citri A, Yarden Y: EGF-ERBB signalling: towards the systems level. Nat Rev Mol Cell Biol. 2006, 7: 505-516. 10.1038/nrm1962.PubMedView Article
  52. Wiley HS, Shvartsman SY, Lauffenburger DA: Computational modeling of the EGF-receptor system: a paradigm for systems biology. Trends Cell Biol. 2003, 13: 43-50. 10.1016/S0962-8924(02)00009-0.PubMedView Article
  53. Bose R, Molina H, Patterson AS, Bitok JK, Periaswamy B, Bader JS, Pandey A, Cole PA: Phosphoproteomic analysis of Her2/neu signaling and inhibition. Proc Natl Acad Sci USA. 2006, 103: 9773-9778. 10.1073/pnas.0603948103.PubMedPubMed CentralView Article
  54. Sevecka M, MacBeath G: State-based discovery: a multidimensional screen for small-molecule modulators of EGF signaling. Nat Methods. 2006, 3: 825-831. 10.1038/nmeth931.PubMedPubMed CentralView Article
  55. Wolf-Yadlin A, Kumar N, Zhang Y, Hautaniemi S, Zaman M, Kim HD, Grantcharova V, Lauffenburger DA, White FM: Effects of HER2 overexpression on cell signaling networks governing proliferation and migration. Mol Syst Biol. 2006, 2: 54-10.1038/msb4100094.PubMedPubMed CentralView Article
  56. Hatakeyama M: System properties of ErbB receptor signaling for the understanding of cancer progression. Mol Biosyst. 2007, 3: 111-116. 10.1039/b612800a.PubMedView Article
  57. Jones RB, Gordus A, Krall JA, MacBeath G: A quantitative protein interaction network for the ErbB receptors using protein micro-arrays. Nature. 2006, 439: 168-174. 10.1038/nature04177.PubMedView Article
  58. Ruths D, Muller M, Tseng JT, Nakhleh L, Ram PT: The signaling petri net-based simulator: a non-parametric strategy for characterizing the dynamics of cell-specific signaling networks. PLoS Comput Biol. 2008, 4: e1000005-10.1371/journal.pcbi.1000005.PubMedPubMed CentralView Article
  59. Dasika MS, Burgard A, Maranas CD: A computational framework for the topological analysis and targeted disruption of signal transduction networks. Biophys J. 2006, 91: 382-398. 10.1529/biophysj.105.069724.PubMedPubMed CentralView Article
  60. Kumar N, Hendriks BS, Janes KA, de Graaf D, Lauffenburger DA: Applying computational modeling to drug discovery and development. Drug Discov Today. 2006, 11: 806-811. 10.1016/j.drudis.2006.07.010.PubMedView Article
  61. Fojo T, Menefee M: Mechanisms of multidrug resistance: the potential role of microtubule-stabilizing agents. Ann Oncol. 2007, 18 (Suppl 5): v3-8. 10.1093/annonc/mdm172.PubMedView Article
  62. Lord CJ, Iorns E, Ashworth A: Dissecting resistance to endocrine therapy in breast cancer. Cell Cycle. 2008, 7: 1895-1898.PubMedView Article
  63. Ciardiello F, Tortora G: EGFR antagonists in cancer treatment. N Engl J Med. 2008, 358: 1160-1174. 10.1056/NEJMra0707704.PubMedView Article
  64. Jin H, Yang R, Zheng Z, Romero M, Ross J, Bou-Reslan H, Carano RA, Kasman I, Mai E, Young J, Zha J, Zhang Z, Ross S, Schwall R, Colbern G, Merchant M: MetMAb, the one-armed 5D5 anti-c-Met antibody, inhibits orthotopic pancreatic tumor growth and improves survival. Cancer Res. 2008, 68: 4360-4368. 10.1158/0008-5472.CAN-07-5960.PubMedView Article
  65. Guo A, Villen J, Kornhauser J, Lee KA, Stokes MP, Rikova K, Possemato A, Nardone J, Innocenti G, Wetzel R, Wang Y, MacNeill J, Mitchell J, Gygi SP, Rush J, Polakiewicz RD, Comb MJ: Signaling networks assembled by oncogenic EGFR and c-Met. Proc Natl Acad Sci USA. 2008, 105: 692-697. 10.1073/pnas.0707270105.PubMedPubMed CentralView Article
  66. Sudo H, Tsuji AB, Sugyo A, Imai T, Saga T, Harada YN: A loss of function screen identifies nine new radiation susceptibility genes. Biochem Biophys Res Commun. 2007, 364: 695-701. 10.1016/j.bbrc.2007.10.074.PubMedView Article
  67. Amundson SA, Do KT, Vinikoor LC, Lee RA, Koch-Paiz CA, Ahn J, Reimers M, Chen Y, Scudiero DA, Weinstein JN, Trent JM, Bittner ML, Meltzer PS, Fornace AJ: Integrating global gene expression and radiation survival parameters across the 60 cell lines of the National Cancer Institute Anticancer Drug Screen. Cancer Res. 2008, 68: 415-424. 10.1158/0008-5472.CAN-07-2120.PubMedView Article
  68. Swanton C, Marani M, Pardo O, Warne PH, Kelly G, Sahai E, Elustondo F, Chang J, Temple J, Ahmed AA, Brenton JD, Downward J, Nicke B: Regulators of mitotic arrest and ceramide metabolism are determinants of sensitivity to paclitaxel and other chemotherapeutic drugs. Cancer Cell. 2007, 11: 498-512. 10.1016/j.ccr.2007.04.011.PubMedView Article
  69. Fitzgerald JB, Schoeberl B, Nielsen UB, Sorger PK: Systems biology and combination therapy in the quest for clinical efficacy. Nat Chem Biol. 2006, 2: 458-466. 10.1038/nchembio817.PubMedView Article
  70. Araujo RP, Petricoin EF, Liotta LA: A mathematical model of combination therapy using the EGFR signaling network. Biosystems. 2005, 80: 57-69. 10.1016/j.biosystems.2004.10.002.PubMedView Article
  71. Cappellini MD, Fiorelli G: Glucose-6-phosphate dehydrogenase deficiency. Lancet. 2008, 371: 64-74. 10.1016/S0140-6736(08)60073-2.PubMedView Article
  72. Bux J: Molecular nature of antigens implicated in immune neutropenias. Int J Hematol. 2002, 76 (Suppl 1): 399-403.PubMedView Article
  73. Opgen-Rhein C, Dettling M: Clozapine-induced agranulocytosis and its genetic determinants. Pharmacogenomics. 2008, 9: 1101-1111. 10.2217/14622416.9.8.1101.PubMedView Article
  74. Roden DM: Cellular basis of drug-induced torsades de pointes. Br J Pharmacol. 2008, 154: 1502-1507. 10.1038/bjp.2008.238.PubMedPubMed CentralView Article
  75. Wagner BK, Kitami T, Gilbert TJ, Peck D, Ramanathan A, Schreiber SL, Golub TR, Mootha VK: Large-scale chemical dissection of mitochondrial function. Nat Biotechnol. 2008, 26: 343-351. 10.1038/nbt1387.PubMedPubMed CentralView Article


© BioMed Central Ltd 2009