Towards Drugs Targeting Multiple Proteins in a Systems Biology Approach

Systematic Investigation of Dual-Target-Directed Ligands

Multitarget-directed ligands (MTDLs) are compounds rationally designed to affect multiple targets, aiming for a better therapeutic profile. For over 20 years, MTDLs have garnered increasing attention, the idea being that their full potential would have been achieved, thanks to unprecedented target combinations and advanced medicinal chemistry strategies. This study presents a literature mining effort resulting in a data set of dual-target-directed ligands (DTDLs), the fundamental example of MTDLs. We used this data set to evaluate the rationale behind target selection and the chemical novelty of DTDLs targeting specific protein combinations. Our analysis focused on DTDL targets in terms of biological function, disease association, structure, and chemogenomic traits. We also compared DTDLs with single-target compounds. We found that well-known target pathology associations often guide DTDL design, leveraging existing chemical scaffolds and binding pocket similarities. These findings highlight the current state of the field and suggest substantial untapped potential for rational polypharmacology.

Systems biology informed neural networks (SBINN) predict response and novel combinations for PD-1 checkpoint blockade

Anti-PD-1 immunotherapy has recently shown tremendous success for the treatment of several aggressive cancers. However, variability and unpredictability in treatment outcome have been observed, and are thought to be driven by patient-specific biology and interactions of the patient's immune system with the tumor. Here we develop an integrative systems biology and machine learning approach, built around clinical data, to predict patient response to anti-PD-1 immunotherapy and to improve the response rate. Using this approach, we determine biomarkers of patient response and identify potential mechanisms of drug resistance. We develop systems biology informed neural networks (SBINN) to calculate patient-specific kinetic parameter values and to predict clinical outcome. We show how transfer learning can be leveraged with simulated clinical data to significantly improve the response prediction accuracy of the SBINN. Further, we identify novel drug combinations and optimize the treatment protocol for triple combination therapy consisting of IL-6 inhibition, recombinant IL-12, and anti-PD-1 immunotherapy in order to maximize patient response. We also find unexpected differences in protein expression levels between response phenotypes which complement recent clinical findings. Our approach has the potential to aid in the development of targeted experiments for patient drug screening as well as identify novel therapeutic targets.

A mean-field approach for modeling the propagation of perturbations in biochemical reaction networks

Often, the time evolution of a biochemical reaction network is crucial for determining the effects of combining multiple pharmaceuticals. Here we illustrate a mathematical framework for modeling the dominant temporal behaviour of a complicated molecular pathway or biochemical reaction network in response to an arbitrary perturbation, such as resulting from the administration of a therapeutic agent. The method enables the determination of the temporal evolution of a target protein as the perturbation propagates through its regulatory network. The mathematical approach is particularly useful when the experimental data that is available for characterizing or parameterizing the regulatory network is limited or incomplete. To illustrate the method, we consider the examples of the regulatory networks for the target proteins c-Myc and Chop, which play an important role in venetoclax resistance in acute myeloid leukemia. First we show how the networks that regulate each target protein can be reduced to a mean-field model by identifying the distinct effects that groups of proteins in the regulatory network have on the target protein. Then we show how limited protein-level data can be used to further simplify the mean-field model to pinpoint the dominant effects of the network perturbation on the target protein. This enables a further reduction in the number of parameters in the model. The result is an ordinary differential equation model that captures the temporal evolution of the expression of a target protein when one or more proteins in its regulatory network have been perturbed. Finally, we show how the dominant effects predicted by the mathematical model agree with RNA sequencing data for the regulatory proteins comprising the molecular network, despite the model not having a priori knowledge of this data. Thus, while the approach gives a simplified model for the expression of the target protein, it allows for the interpretation of the effects of the perturbation on the regulatory network itself. This method can be easily extended to sets of target proteins to model components of a larger systems biology model, and provides an approach for partially integrating RNA sequencing data and protein expression data. Moreover, it is a general approach that can be used to study drug effects on specific protein(s) in any disease or condition.

Big data, integrative omics and network biology

A cell integrates various signals through a network of biomolecules that crosstalk to synergistically regulate the replication, transcription, translation and other metabolic activities of a cell. These networks regulate signal perception and processing that drives biological functions. The biological complexity cannot be fully captured by a single -omics discipline. The holistic study of an organism-in health, perturbation, exposure to environment and disease, is studied under systems biology. The bottom-up molecular approaches (genes, mRNA, protein, metabolite, etc.) have laid the foundation of current biological knowledge covering the horizon from viruses, bacteria, fungi, plants and animals. Yet, these techniques provide a rather myopic view of biology at the molecular level. To understand how the interconnected molecular components are formed and rewired in disease or exposure to environmental stimuli is the holy grail of modern biology. The omics era was heralded by the genomics revolution but advanced sequencing techniques are now also ubiquitous in transcriptomics, proteomics, metabolomics and lipidomics. Multi-omics data analysis and integration techniques are driving the quest for deeper insights into how the different layers of biomolecules talk to each other in diverse contexts.

Biophysical Insight into the Interaction of Human Lysozyme with Anticancer Drug Anastrozole: A Multitechnique Approach

In the present study, we employ fluorescence spectroscopy, dynamic light scattering, and molecular docking methods. Binding of anticancer drug anastrozole with human lysozyme (HL) is studied. Binding of anastrozole to HL is moderate but spontaneous. There is anastrozole persuaded hydrodynamic change in HL, leading to molecular compaction. Binding of anastrozole to HL also decreased in vitro lytic activity of HL. Molecular docking results suggest the electrostatic interactions and van der Waals forces played key role in binding interaction of anastrozole near the catalytic site. Binding interaction of anastrozole to proteins other than major transport proteins in blood can significantly affect pharmacokinetics of this molecule. Hence, rationalizing drug dosage is important. This study also points to unrelated effects that small molecules bring in the body that are considerable and need thorough investigation.

Molecular docking as a popular tool in drug design, an in silico travel

New molecular modeling approaches, driven by rapidly improving computational platforms, have allowed many success stories for the use of computer-assisted drug design in the discovery of new mechanism-or structure-based drugs. In this overview, we highlight three aspects of the use of molecular docking. First, we discuss the combination of molecular and quantum mechanics to investigate an unusual enzymatic mechanism of a flavoprotein. Second, we present recent advances in anti-infectious agents' synthesis driven by structural insights. At the end, we focus on larger biological complexes made by protein-protein interactions and discuss their relevance in drug design. This review provides information on how these large systems, even in the presence of the solvent, can be investigated with the outlook of drug discovery.

"Omics" in pharmaceutical research: overview, applications, challenges, and future perspectives

In the post-genomic era, biological studies are characterized by the rapid development and wide application of a series of "omics" technologies, including genomics, proteomics, metabolomics, transcriptomics, lipidomics, cytomics, metallomics, ionomics, interactomics, and phenomics. These "omics" are often based on global analyses of biological samples using high through-put analytical approaches and bioinformatics and may provide new insights into biological phenomena. In this paper, the development and advances in these omics made in the past decades are reviewed, especially genomics, transcriptomics, proteomics and metabolomics; the applications of omics technologies in pharmaceutical research are then summarized in the fields of drug target discovery, toxicity evaluation, personalized medicine, and traditional Chinese medicine; and finally, the limitations of omics are discussed, along with the future challenges associated with the multi-omics data processing, dynamics omics analysis, and analytical approaches, as well as amenable solutions and future prospects.

BalestraWeb: efficient online evaluation of drug-target interactions

SUMMARY

BalestraWeb is an online server that allows users to instantly make predictions about the potential occurrence of interactions between any given drug-target pair, or predict the most likely interaction partners of any drug or target listed in the DrugBank. It also permits users to identify most similar drugs or most similar targets based on their interaction patterns. Outputs help to develop hypotheses about drug repurposing as well as potential side effects.

AVAILABILITY AND IMPLEMENTATION

BalestraWeb is accessible at http://balestra.csb.pitt.edu/. The tool is built using a probabilistic matrix factorization method and DrugBank v3, and the latent variable models are trained using the GraphLab collaborative filtering toolkit. The server is implemented using Python, Flask, NumPy and SciPy.

Homopharma: a new concept for exploring the molecular binding mechanisms and drug repurposing

Background

Drugs that simultaneously target multiple proteins often improve efficacy, particularly in the treatment of complex diseases such as cancers and central nervous system disorders. Many approaches have been proposed to identify the potential targets of a drug. Recently, we have introduced Space-Related Pharmamotif (SRPmotif) method to recognize the proteins that share similar binding environments. In addition, compounds with similar topology may bind to similar proteins and have similar protein-compound interactions. However, few studies have focused on exploring the relationships between binding environments and protein-compound interactions, which is important for understanding molecular binding mechanisms and helpful to be used in discovering drug repurposing.

Results

In this study, we propose a new concept of "Homopharma", combining similar binding environments and protein-compound interaction profiles, to explore the molecular binding mechanisms and drug repurposing. A Homopharma consists of a set of proteins which have the conserved binding environment and a set of compounds that share similar structures and functional groups. These proteins and compounds present conserved interactions and similar physicochemical properties. Therefore, these compounds are often able to inhibit the proteins in a Homopharma. Our experimental results show that the proteins and compounds in a Homopharma often have similar protein-compound interactions, comprising conserved specific residues and functional sites. Based on the Homopharma concept, we selected four flavonoid derivatives and 32 human protein kinases for enzymatic profiling. Among these 128 bioassays, the IC₅₀ of 56 and 25 flavonoid-kinase inhibitions are less than 10 μM and 1 μM, respectively. Furthermore, these experimental results suggest that these flavonoids can be used as anticancer compounds, such as oral and colorectal cancer drugs.

Conclusions

The experimental results show that the Homopharma is useful for identifying key binding environments of proteins and compounds and discovering new inhibitory effects. We believe that the Homopharma concept can have the potential for understanding molecular binding mechanisms and providing new clues for drug development.

Discovery of novel drugs for promising targets

BACKGROUND

Once a promising drug target is identified, the steps to actually discover and optimize a drug are diverse and challenging.

OBJECTIVE

The goal of this study was to provide a road map to navigate drug discovery.

METHODS

Review general steps for drug discovery and provide illustrating references.

RESULTS

A number of approaches are available to enhance and accelerate target identification and validation. Consideration of a variety of potential mechanisms of action of potential drugs can guide discovery efforts. The hit to lead stage may involve techniques such as high-throughput screening, fragment-based screening, and structure-based design, with informatics playing an ever-increasing role. Biologically relevant screening models are discussed, including cell lines, 3-dimensional culture, and in vivo screening. The process of enabling human studies for an investigational drug is also discussed.

CONCLUSIONS

Drug discovery is a complex process that has significantly evolved in recent years.

Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review

Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.

Computational prediction of protein hot spot residues

Most biological processes involve multiple proteins interacting with each other. It has been recently discovered that certain residues in these protein-protein interactions, which are called hot spots, contribute more significantly to binding affinity than others. Hot spot residues have unique and diverse energetic properties that make them challenging yet important targets in the modulation of protein-protein complexes. Design of therapeutic agents that interact with hot spot residues has proven to be a valid methodology in disrupting unwanted protein-protein interactions. Using biological methods to determine which residues are hot spots can be costly and time consuming. Recent advances in computational approaches to predict hot spots have incorporated a myriad of features, and have shown increasing predictive successes. Here we review the state of knowledge around protein-protein interactions, hot spots, and give an overview of multiple in silico prediction techniques of hot spot residues.

Proline rich motifs as drug targets in immune mediated disorders

The current version of the human immunome network consists of nearly 1400 interactions involving approximately 600 proteins. Intermolecular interactions mediated by proline-rich motifs (PRMs) are observed in many facets of the immune response. The proline-rich regions are known to preferentially adopt a polyproline type II helical conformation, an extended structure that facilitates transient intermolecular interactions such as signal transduction, antigen recognition, cell-cell communication and cytoskeletal organization. The propensity of both the side chain and the backbone carbonyls of the polyproline type II helix to participate in the interface interaction makes it an excellent recognition motif. An advantage of such distinct chemical features is that the interactions can be discriminatory even in the absence of high affinities. Indeed, the immune response is mediated by well-orchestrated low-affinity short-duration intermolecular interactions. The proline-rich regions are predominantly localized in the solvent-exposed regions such as the loops, intrinsically disordered regions, or between domains that constitute the intermolecular interface. Peptide mimics of the PRM have been suggested as potential antagonists of intermolecular interactions. In this paper, we discuss novel PRM-mediated interactions in the human immunome that potentially serve as attractive targets for immunomodulation and drug development for inflammatory and autoimmune pathologies.

Expanding the conformational selection paradigm in protein-ligand docking

Conformational selection emerges as a theme in macromolecular interactions. Data validate it as a prevailing mechanism in protein-protein, protein-DNA, protein-RNA, and protein-small molecule drug recognition. This raises the question of whether this fundamental biomolecular binding mechanism can be used to improve drug docking and discovery. Actually, in practice this has already been taking place for some years in increasing numbers. Essentially, it argues for using not a single conformer, but an ensemble. The paradigm of conformational selection holds that because the ensemble is heterogeneous, within it there will be states whose conformation matches that of the ligand. Even if the population of this state is low, since it is favorable for binding the ligand, it will bind to it with a subsequent population shift toward this conformer. Here we suggest expanding it by first modeling all protein interactions in the cell by using Prism, an efficient motif-based protein-protein interaction modeling strategy, followed by ensemble generation. Such a strategy could be particularly useful for signaling proteins, which are major targets in drug discovery and bind multiple partners through a shared binding site, each with some-minor or major-conformational change.

Curcumin analogues as possible anti-proliferative & anti-inflammatory agents

A series of novel curcumin analogues has been designed, synthesized and tested in vitro/in vivo as potential multi-target agents. Their anti-proliferative and anti-inflammatory activities were studied. Compounds 1b and 2b were stronger inhibitors of soybean lipoxygenase (LOX) than curcumin. Analogue 1b was also the most potent aldose reductase (ALR2) inhibitor. Two compounds, (1a and 1f) exhibited in vivo anti-inflammatory activity comparable to that of indomethacin, whereas derivative 1i exhibited even higher activity. The derivatives were also tested for their anti-proliferative activity using three different human cancer cell lines. Compounds 1a, 1b, 1d and 2b exhibited significant growth inhibitory activity as compared to curcumin, against all three cancer cell lines. Lipophilicity was determined as R(M) values using RPTLC and theoretically. The results are discussed in terms of the structural characteristics of the compounds. Docking simulations were performed on LOX and ALR2 inhibitor 1b and curcumin. Compound 1b is well fitted in the active site of ALR2, binding to the ALR2 enzyme in a similar way to curcumin. Allosteric interactions may govern the LOX-inhibitor binding.

The potassium channel KcsA: a model protein in studying membrane protein oligomerization and stability of oligomeric assembly?

Many membrane proteins are functional as stable oligomers. An understanding of the conditions that elicit and enhance oligomerization is important in many therapeutics. In this regard, protein-protein and protein-lipid interactions play crucial roles in the assembly and stability of oligomeric complexes. Recent years have seen a rapid increase in the mechanistic information on the importance of cytoplasmic termini in determining subunit assembly and stability of oligomeric complexes. In addition, the role of specific protein-lipid interaction between anionic phospholipids and "hot spots" on the protein surface has also become evident in stabilizing oligomeric assemblies. This review focuses on several contemporary developments of membrane proteins that stabilize oligomers by taking the potassium channel KcsA as an exemplary ion channel.

Disruption of cell wall fatty acid biosynthesis in Mycobacterium tuberculosis using a graph theoretic approach

Fatty acid biosynthesis of Mycobacterium tuberculosis was analyzed using graph theory and influential (impacting) proteins were identified. The graphs (digraphs) representing this biological network provide information concerning the connectivity of each protein or metabolite in a given pathway, providing an insight into the importance of various components in the pathway, and this can be quantitatively analyzed. Using a graph theoretic algorithm, the most influential set of proteins (sets of {1, 2, 3}, etc.), which when eliminated could cause a significant impact on the biosynthetic pathway, were identified. This set of proteins could serve as drug targets. In the present study, the metabolic network of Mycobacterium tuberculosis was constructed and the fatty acid biosynthesis pathway was analyzed for potential drug targeting. The metabolic network was constructed using the KEGG LIGAND database and subjected to graph theoretical analysis. The nearness index of a protein was used to determine the influence of the said protein on other components in the network, allowing the proteins in a pathway to be ordered according to their nearness indices. A method for identifying the most strategic nodes to target for disrupting the metabolic networks is proposed, aiding the development of new drugs to combat this deadly disease.

Systems biology approaches and tools for analysis of interactomes and multi-target drugs

Systems biology is essentially a proteomic and epigenetic exercise because the relatively condensed information of genomes unfolds on the level of proteins. The flexibility of cellular architectures is not only mediated by a dazzling number of proteinaceous species but moreover by the kinetics of their molecular changes: The time scales of posttranslational modifications range from milliseconds to years. The genetic framework of an organism only provides the blue print of protein embodiments which are constantly shaped by external input. Indeed, posttranslational modifications of proteins represent the scope and velocity of these inputs and fulfil the requirements of integration of external spatiotemporal signal transduction inside an organism. The optimization of biochemical networks for this type of information processing and storage results in chemically extremely fine tuned molecular entities. The huge dynamic range of concentrations, the chemical diversity and the necessity of synchronisation of complex protein expression patterns pose the major challenge of systemic analysis of biological models. One further message is that many of the key reactions in living systems are essentially based on interactions of moderate affinities and moderate selectivities. This principle is responsible for the enormous flexibility and redundancy of cellular circuitries. In complex disorders such as cancer or neurodegenerative diseases, which initially appear to be rooted in relatively subtle dysfunctions of multimodal physiologic pathways, drug discovery programs based on the concept of high affinity/high specificity compounds ("one-target, one-disease"), which has been dominating the pharmaceutical industry for a long time, increasingly turn out to be unsuccessful. Despite improvements in rational drug design and high throughput screening methods, the number of novel, single-target drugs fell much behind expectations during the past decade, and the treatment of "complex diseases" remains a most pressing medical need. Currently, a change of paradigm can be observed with regard to a new interest in agents that modulate multiple targets simultaneously, essentially "dirty drugs." Targeting cellular function as a system rather than on the level of the single target, significantly increases the size of the drugable proteome and is expected to introduce novel classes of multi-target drugs with fewer adverse effects and toxicity. Multiple target approaches have recently been used to design medications against atherosclerosis, cancer, depression, psychosis and neurodegenerative diseases. A focussed approach towards "systemic" drugs will certainly require the development of novel computational and mathematical concepts for appropriate modelling of complex data. But the key is the extraction of relevant molecular information from biological systems by implementing rigid statistical procedures to differential proteomic analytics.

Intrinsically disordered proteins in a physics-based world

Intrinsically disordered proteins (IDPs) are a newly recognized class of functional proteins that rely on a lack of stable structure for function. They are highly prevalent in biology, play fundamental roles, and are extensively involved in human diseases. For signaling and regulation, IDPs often fold into stable structures upon binding to specific targets. The mechanisms of these coupled binding and folding processes are of significant importance because they underlie the organization of regulatory networks that dictate various aspects of cellular decision-making. This review first discusses the challenge in detailed experimental characterization of these heterogeneous and dynamics proteins and the unique and exciting opportunity for physics-based modeling to make crucial contributions, and then summarizes key lessons from recent de novo simulations of the structure and interactions of several regulatory IDPs.

Systems biology and modeling in neuroblastoma: practicalities and perspectives

Neuroblastoma (NB) is a common pediatric malignancy characterized by clinical and biological heterogeneity. A host of prognostic markers are available, contributing to accurate risk stratification and appropriate treatment allocation. Unfortunately, outcome is still poor for many patients, indicating the need for a new approach with enhanced utilization of the available biological data. Systems biology is a holistic approach in which all components of a biological system carry equal importance. Systems biology uses mathematical modeling and simulation to investigate dynamic interactions between system components, as a means of explaining overall system behavior. Systems biology can benefit the biomedical sciences by providing a more complete understanding of human disease, enhancing the development of targeted therapeutics. Systems biology is largely contiguous with current approaches in NB, which already employ an integrative and pseudo-holistic approach to disease management. Systems modeling of NB offers an optimal method for continuing progression in this field, and conferring additional benefit to current risk stratification and management. Likewise, NB provides an opportunity for systems biology to prove its utility in the context of human disease, since the biology of NB is comprehensively characterized and, therefore, suited to modeling. The purpose of this review is to outline the benefits, challenges and fundamental workings of systems modeling in human disease, using a specific example of bottom-up modeling in NB. The intention is to demonstrate practical requirements to begin bridging the gap between biological research and applied mathematical approaches for the mutual gain of both fields, and with additional benefits for clinical management.

Towards inferring time dimensionality in protein-protein interaction networks by integrating structures: the p53 example

Structural data, efficient structural comparison algorithms and appropriate datasets and filters can assist in getting an insight into time dimensionality in interaction networks; in predicting which interactions can and cannot co-exist; and in obtaining concrete predictions consistent with experiment.

Inspection of protein–protein interaction maps illustrates that a hub protein can interact with a very large number of proteins, reaching tens and even hundreds. Since a single protein cannot interact with such a large number of partners at the same time, this presents a challenge: can we figure out which interactions can occur simultaneously and which are mutually excluded? Addressing this question adds a fourth dimension into interaction maps: that of time. Including the time dimension in structural networks is an immense asset; time dimensionality transforms network node-and-edge maps into cellular processes, assisting in the comprehension of cellular pathways and their regulation. While the time dimensionality can be further enhanced by linking protein complexes to time series of mRNA expression data, current robust, network experimental data are lacking. Here we outline how, using structural data, efficient structural comparison algorithms and appropriate datasets and filters can assist in getting an insight into time dimensionality in interaction networks; in predicting which interactions can and cannot co-exist; and in obtaining concrete predictions consistent with experiment. As an example, we present p53-linked processes.

A hybrid clustering of protein binding sites

The Protein Data Bank contains the description of approximately 27 000 protein-ligand binding sites. Most of the ligands at these sites are biologically active small molecules, affecting the biological function of the protein. The classification of their binding sites may lead to relevant results in drug discovery and design. Clusters of similar binding sites were created here by a hybrid, sequence and spatial structure-based approach, using the OPTICS clustering algorithm. A dissimilarity measure was defined: a distance function on the amino acid sequences of the binding sites. All the binding sites were clustered in the Protein Data Bank according to this distance function, and it was found that the clusters characterized well the Enzyme Commission numbers of the entries. The results, carefully color coded by the Enzyme Commission numbers of the proteins, containing the 20 967 binding sites clustered, are available as html files in three parts at http://pitgroup.org/seqclust/.

Protein-protein interaction networks: how can a hub protein bind so many different partners?

How can a single hub protein bind so many different partners? Numerous studies have sought differences between hubs and non-hubs to explain what makes a protein a hub and how a shared hub-binding site can be promiscuous, yet at the same time be specific. Here, we suggest that the problem is largely non-existent and resides in the popular representation of protein interaction networks: protein products derived from a single gene, even if different, are clustered in maps into a single node. This leads to the impression that a single protein binds to a very large number of partners. In reality, it does not; rather, protein networks reflect the combination of multiple proteins, each with a distinct conformation.

Evaluating docking methods for prediction of binding affinities of small molecules to the G protein betagamma subunits

Several studies have suggested that disrupting interactions of the G protein betagamma subunits with downstream binding partners might be a valuable study for pharmaceutical development. Recently, small molecules have been found which bind to Gbetagamma with high apparent affinity in an enzyme-linked immunosorbent assay (ELISA), have demonstrated selective inhibition of interactions of Gbetagamma with downstream signaling partners, and have been shown to increase antinociceptive effects of morphine and inhibit inflammation in vivo. In this paper we examine several docking and scoring protocols for estimating binding affinities for a set of 830 ligands from the NCI diversity set to the Gbeta1gamma2 subunit and compared these with IC50s measured in a competition ELISA with a high-affinity peptidic ligand. The best-performing docking protocol used a consensus score and ensemble docking and resulted in a 6-fold enrichment of high-affinity compounds in the top-ranked 5% of the ligand data set.

Identification of differentially expressed subnetworks based on multivariate ANOVA

Background

Since high-throughput protein-protein interaction (PPI) data has recently become available for humans, there has been a growing interest in combining PPI data with other genome-wide data. In particular, the identification of phenotype-related PPI subnetworks using gene expression data has been of great concern. Successful integration for the identification of significant subnetworks requires the use of a search algorithm with a proper scoring method. Here we propose a multivariate analysis of variance (MANOVA)-based scoring method with a greedy search for identifying differentially expressed PPI subnetworks.

Results

Given the MANOVA-based scoring method, we performed a greedy search to identify the subnetworks with the maximum scores in the PPI network. Our approach was successfully applied to human microarray datasets. Each identified subnetwork was annotated with the Gene Ontology (GO) term, resulting in the phenotype-related functional pathway or complex. We also compared these results with those of other scoring methods such as t statistic- and mutual information-based scoring methods. The MANOVA-based method produced subnetworks with a larger number of proteins than the other methods. Furthermore, the subnetworks identified by the MANOVA-based method tended to consist of highly correlated proteins.

Conclusion

This article proposes a MANOVA-based scoring method to combine PPI data with expression data using a greedy search. This method is recommended for the highly sensitive detection of large subnetworks.

Network pharmacology: the next paradigm in drug discovery

The dominant paradigm in drug discovery is the concept of designing maximally selective ligands to act on individual drug targets. However, many effective drugs act via modulation of multiple proteins rather than single targets. Advances in systems biology are revealing a phenotypic robustness and a network structure that strongly suggests that exquisitely selective compounds, compared with multitarget drugs, may exhibit lower than desired clinical efficacy. This new appreciation of the role of polypharmacology has significant implications for tackling the two major sources of attrition in drug development--efficacy and toxicity. Integrating network biology and polypharmacology holds the promise of expanding the current opportunity space for druggable targets. However, the rational design of polypharmacology faces considerable challenges in the need for new methods to validate target combinations and optimize multiple structure-activity relationships while maintaining drug-like properties. Advances in these areas are creating the foundation of the next paradigm in drug discovery: network pharmacology.

Knowledgebase for addiction-related genes: is it possible an extrapolation to rational multi-target drug design?

In recent years the single-probe-single-target approach in drug design has started to be smoothly replaced by the single-probe-multiple-target (or multi-target) one, where a single drug is able to tackle different, but disease-related targets in a selective manner. However, the design of multi-target drugs has been hindered by a lack of a systematic network of disease-related common pathways. The recent development of the knowledgebase of addiction-related genes (KARG) has provided important hints on how to rationally design multi-target probes by connecting experimental techniques with available network models. In this perspective, the use of KARG as a template to build knowledgebases of disease-related genes for the rational multi-target drug design is suggested. Moreover, it is proposed that building knowledgebases of disease-related genes will become a necessary and ubiquitous tool in the rational design of multi-target drugs.

Interactive proteomics: what lies ahead?

Interactive proteomics addresses the physical associations among proteins and establishes global, disease-, and pathway-specific protein interaction networks. The inherent chemical and structural diversity of proteins, their different expression levels, and their distinct subcellular localizations pose unique challenges for the exploration of these networks, necessitating the use of a variety of innovative and ingenious approaches. Consequently, recent years have seen exciting developments in protein interaction mapping and the establishment of very large interaction networks, especially in model organisms. In the near future, attention will shift to the establishment of interaction networks in humans and their application in drug discovery and understanding of diseases. In this review, we present an impressive toolbox of different technologies that we expect to be crucial for interactive proteomics in the coming years.

Computer-aided drug design: the next 20 years

This perspectives article has been taken from a talk the author gave at the symposium in honor of Yvonne C. Martin's retirement, held at the American Chemical Society spring meeting in Chicago on March 25, 2007. The talk was intended as a somewhat lighthearted attempt to gaze into the future; inevitably, in print, things will come across more seriously than was intended. As we all know-the past is rarely predictive of the future.