Systems Biology Brings New Dimensions for Structure-Based Drug Design

Diverse ways of perturbing the human arachidonic acid metabolic network to control inflammation

Inflammation and other common disorders including diabetes, cardiovascular disease, and cancer are often the result of several molecular abnormalities and are not likely to be resolved by a traditional single-target drug discovery approach. Though inflammation is a normal bodily reaction, uncontrolled and misdirected inflammation can cause inflammatory diseases such as rheumatoid arthritis and asthma. Nonsteroidal anti-inflammatory drugs including aspirin, ibuprofen, naproxen, or celecoxib are commonly used to relieve aches and pains, but often these drugs have undesirable and sometimes even fatal side effects. To facilitate safer and more effective anti-inflammatory drug discovery, a balanced treatment strategy should be developed at the biological network level. In this Account, we focus on our recent progress in modeling the inflammation-related arachidonic acid (AA) metabolic network and subsequent multiple drug design. We first constructed a mathematical model of inflammation based on experimental data and then applied the model to simulate the effects of commonly used anti-inflammatory drugs. Our results indicated that the model correctly reproduced the established bleeding and cardiovascular side effects. Multitarget optimal intervention (MTOI), a Monte Carlo simulated annealing based computational scheme, was then developed to identify key targets and optimal solutions for controlling inflammation. A number of optimal multitarget strategies were discovered that were both effective and safe and had minimal associated side effects. Experimental studies were performed to evaluate these multitarget control solutions further using different combinations of inhibitors to perturb the network. Consequently, simultaneous control of cyclooxygenase-1 and -2 and leukotriene A4 hydrolase, as well as 5-lipoxygenase and prostaglandin E2 synthase were found to be among the best solutions. A single compound that can bind multiple targets presents advantages including low risk of drug-drug interactions and robustness regarding concentration fluctuations. Thus, we developed strategies for multiple-target drug design and successfully discovered several series of multiple-target inhibitors. Optimal solutions for a disease network often involve mild but simultaneous interventions of multiple targets, which is in accord with the philosophy of traditional Chinese medicine (TCM). To this end, our AA network model can aptly explain TCM anti-inflammatory herbs and formulas at the molecular level. We also aimed to identify activators for several enzymes that appeared to have increased activity based on MTOI outcomes. Strategies were then developed to predict potential allosteric sites and to discover enzyme activators based on our hypothesis that combined treatment with the projected activators and inhibitors could balance different AA network pathways, control inflammation, and reduce associated adverse effects. Our work demonstrates that the integration of network modeling and drug discovery can provide novel solutions for disease control, which also calls for new developments in drug design concepts and methodologies. With the rapid accumulation of quantitative data and knowledge of the molecular networks of disease, we can expect an increase in the development and use of quantitative disease models to facilitate efficient and safe drug discovery.

Molecular Determinants Underlying Binding Specificities of the ABL Kinase Inhibitors: Combining Alanine Scanning of Binding Hot Spots with Network Analysis of Residue Interactions and Coevolution

Quantifying binding specificity and drug resistance of protein kinase inhibitors is of fundamental importance and remains highly challenging due to complex interplay of structural and thermodynamic factors. In this work, molecular simulations and computational alanine scanning are combined with the network-based approaches to characterize molecular determinants underlying binding specificities of the ABL kinase inhibitors. The proposed theoretical framework unveiled a relationship between ligand binding and inhibitor-mediated changes in the residue interaction networks. By using topological parameters, we have described the organization of the residue interaction networks and networks of coevolving residues in the ABL kinase structures. This analysis has shown that functionally critical regulatory residues can simultaneously embody strong coevolutionary signal and high network centrality with a propensity to be energetic hot spots for drug binding. We have found that selective (Nilotinib) and promiscuous (Bosutinib, Dasatinib) kinase inhibitors can use their energetic hot spots to differentially modulate stability of the residue interaction networks, thus inhibiting or promoting conformational equilibrium between inactive and active states. According to our results, Nilotinib binding may induce a significant network-bridging effect and enhance centrality of the hot spot residues that stabilize structural environment favored by the specific kinase form. In contrast, Bosutinib and Dasatinib can incur modest changes in the residue interaction network in which ligand binding is primarily coupled only with the identity of the gate-keeper residue. These factors may promote structural adaptability of the active kinase states in binding with these promiscuous inhibitors. Our results have related ligand-induced changes in the residue interaction networks with drug resistance effects, showing that network robustness may be compromised by targeted mutations of key mediating residues. This study has outlined mechanisms by which inhibitor binding could modulate resilience and efficiency of allosteric interactions in the kinase structures, while preserving structural topology required for catalytic activity and regulation.

Maximizing computational tools for successful drug discovery

Drug discovery is an iterative cycle of identifying promising hits followed by lead optimization via bioisosteric replacements. In the search for compounds affording good bioactivity, equal importance should also be placed on achieving those with favorable pharmacokinetic properties. Thus, the balance and realization of both key properties is an intricate problem that requires great caution. In this editorial, the authors explore the available computational tools in the context of the extant of big data that has borne out via advents of the Omics revolution. As such, the selection of appropriate computational tools for analyzing the vast number of chemical libraries, target proteins and interactomes is the first step toward maximizing the chance for success. However, in order to realize this, it is also necessary to have a solid foundation on the big concepts of drug discovery as well as knowing which tools are available in order to give drug discovery scientists the best opportunity.