Anil Korkut

Perturbation biology: inferring signaling networks in cellular systems.

We present a powerful experimental-computational technology for inferring network models that predict the response of cells to perturbations, and that may be useful in the design of combinatorial therapy against cancer. The experiments are systematic series of perturbations of cancer cell lines by targeted drugs, singly or in combination. The response to perturbation is quantified in terms of relative changes in the measured levels of proteins, phospho-proteins and cellular phenotypes such as viability. Computational network models are derived de novo, i.e., without prior knowledge of signaling pathways, and are based on simple non-linear differential equations. The prohibitively large solution space of all possible network models is explored efficiently using a probabilistic algorithm, Belief Propagation (BP), which is three orders of magnitude faster than standard Monte Carlo methods. Explicit executable models are derived for a set of perturbation experiments in SKMEL-133 melanoma cell lines, which are resistant to the therapeutically important inhibitor of RAF kinase. The resulting network models reproduce and extend known pathway biology. They empower potential discoveries of new molecular interactions and predict efficacious novel drug perturbations, such as the inhibition of PLK1, which is verified experimentally. This technology is suitable for application to larger systems in diverse areas of molecular biology.

Drug Synergy Screen and Network Modeling in Dedifferentiated Liposarcoma Identifies CDK4 and IGF1R as Synergistic Drug Targets

Drug screening and computational modeling of oncogenic signaling pathways identifies synergistic drug pairs for liposarcoma. Predicting Synergistic Therapies Identifying oncogenic targets has improved therapeutic outcomes for cancer patients, but cancers notoriously show primary or acquired resistance to single-agent therapies. Using computational modeling derived from cell viability and high-throughput proteomics data, Miller et al. identified several synergistic pairs of targets in dedifferentiated liposarcoma (DDLS). In two patient-derived DDLS cell lines, combined inhibition of CDK4 (cyclin-dependent kinase 4) and IGF1R (insulin-like growth factor 1 receptor) induced a synergistic decrease in cell proliferation by repressing two pathways: that of retinoblastoma by CDK4 inhibitors and that of AKT and mTOR (mammalian target of rapamycin) by IGF1R inhibitors. The findings suggest that dual inhibition of CDK4 and IGF1R may be a treatment strategy for DDLS and that computational modeling may be applied to various cancers to predict improved combination therapies. Dedifferentiated liposarcoma (DDLS) is a rare but aggressive cancer with high recurrence and low response rates to targeted therapies. Increasing treatment efficacy may require combinations of targeted agents that counteract the effects of multiple abnormalities. To identify a possible multicomponent therapy, we performed a combinatorial drug screen in a DDLS-derived cell line and identified cyclin-dependent kinase 4 (CDK4) and insulin-like growth factor 1 receptor (IGF1R) as synergistic drug targets. We measured the phosphorylation of multiple proteins and cell viability in response to systematic drug combinations and derived computational models of the signaling network. These models predict that the observed synergy in reducing cell viability with CDK4 and IGF1R inhibitors depends on the activity of the AKT pathway. Experiments confirmed that combined inhibition of CDK4 and IGF1R cooperatively suppresses the activation of proteins within the AKT pathway. Consistent with these findings, synergistic reductions in cell viability were also found when combining CDK4 inhibition with inhibition of either AKT or epidermal growth factor receptor (EGFR), another receptor similar to IGF1R that activates AKT. Thus, network models derived from context-specific proteomic measurements of systematically perturbed cancer cells may reveal cancer-specific signaling mechanisms and aid in the design of effective combination therapies.

A force field for virtual atom molecular mechanics of proteins

Activities of many biological macromolecules involve large conformational transitions for which crystallography can specify atomic details of alternative end states, but the course of transitions is often beyond the reach of computations based on full-atomic potential functions. We have developed a coarse-grained force field for molecular mechanics calculations based on the virtual interactions of Cα atoms in protein molecules. This force field is parameterized based on the statistical distribution of the energy terms extracted from crystallographic data, and it is formulated to capture features dependent on secondary structure and on residue-specific contact information. The resulting force field is applied to energy minimization and normal mode analysis of several proteins. We find robust convergence in minimizations to low energies and energy gradients with low degrees of structural distortion, and atomic fluctuations calculated from the normal mode analyses correlate well with the experimental B-factors obtained from high-resolution crystal structures. These findings suggest that the virtual atom force field is a suitable tool for various molecular mechanics applications on large macromolecular systems undergoing large conformational changes.

Pan-Cancer Analysis of Mutation Hotspots in Protein Domains

In cancer genomics, recurrence of mutations in independent tumor samples is a strong indicator of functional impact. However, rare functional mutations can escape detection by recurrence analysis owing to lack of statistical power. We enhance statistical power by extending the notion of recurrence of mutations from single genes to gene families that share homologous protein domains. Domain mutation analysis also sharpens the functional interpretation of the impact of mutations, as domains more succinctly embody function than entire genes. By mapping mutations in 22 different tumor types to equivalent positions in multiple sequence alignments of domains, we confirm well-known functional mutation hotspots, identify uncharacterized rare variants in one gene that are equivalent to well-characterized mutations in another gene, detect previously unknown mutation hotspots, and provide hypotheses about molecular mechanisms and downstream effects of domain mutations. With the rapid expansion of cancer genomics projects, protein domain hotspot analysis will likely provide many more leads linking mutations in proteins to the cancer phenotype.

Computation of conformational transitions in proteins by virtual atom molecular mechanics as validated in application to adenylate kinase

Many proteins function through conformational transitions between structurally disparate states, and there is a need to explore transition pathways between experimentally accessible states by computation. The sizes of systems of interest and the scale of conformational changes are often beyond the scope of full atomic models, but appropriate coarse-grained approaches can capture significant features. We have designed a comprehensive knowledge-based potential function based on a Cα representation for proteins that we call the virtual atom molecular mechanics (VAMM) force field. Here, we describe an algorithm for using the VAMM potential to describe conformational transitions, and we validate this algorithm in application to a transition between open and closed states of adenylate kinase (ADK). The VAMM algorithm computes normal modes for each state and iteratively moves each structure toward the other through a series of intermediates. The move from each side at each step is taken along that normal mode showing greatest engagement with the other state. The process continues to convergence of terminal intermediates to within a defined limit—here, a root-mean-square deviation of 1 Å. Validations show that the VAMM algorithm is highly effective, and the transition pathways examined for ADK are compatible with other structural and biophysical information. We expect that the VAMM algorithm can address many biological systems.

Perturbation biology nominates upstream–downstream drug combinations in RAF inhibitor resistant melanoma cells

Resistance to targeted cancer therapies is an important clinical problem. The discovery of anti-resistance drug combinations is challenging as resistance can arise by diverse escape mechanisms. To address this challenge, we improved and applied the experimental-computational perturbation biology method. Using statistical inference, we build network models from high-throughput measurements of molecular and phenotypic responses to combinatorial targeted perturbations. The models are computationally executed to predict the effects of thousands of untested perturbations. In RAF-inhibitor resistant melanoma cells, we measured 143 proteomic/phenotypic entities under 89 perturbation conditions and predicted c-Myc as an effective therapeutic co-target with BRAF or MEK. Experiments using the BET bromodomain inhibitor JQ1 affecting the level of c-Myc protein and protein kinase inhibitors targeting the ERK pathway confirmed the prediction. In conclusion, we propose an anti-cancer strategy of co-targeting a specific upstream alteration and a general downstream point of vulnerability to prevent or overcome resistance to targeted drugs. DOI: http://dx.doi.org/10.7554/eLife.04640.001

Structural Plasticity and Conformational Transitions of HIV Envelope Glycoprotein gp120

HIV envelope glycoproteins undergo large-scale conformational changes as they interact with cellular receptors to cause the fusion of viral and cellular membranes that permits viral entry to infect targeted cells. Conformational dynamics in HIV gp120 are also important in masking conserved receptor epitopes from being detected for effective neutralization by the human immune system. Crystal structures of HIV gp120 and its complexes with receptors and antibody fragments provide high-resolution pictures of selected conformational states accessible to gp120. Here we describe systematic computational analyses of HIV gp120 plasticity in such complexes with CD4 binding fragments, CD4 mimetic proteins, and various antibody fragments. We used three computational approaches: an isotropic elastic network analysis of conformational plasticity, a full atomic normal mode analysis, and simulation of conformational transitions with our coarse-grained virtual atom molecular mechanics (VAMM) potential function. We observe collective sub-domain motions about hinge points that coordinate those motions, correlated local fluctuations at the interfacial cavity formed when gp120 binds to CD4, and concerted changes in structural elements that form at the CD4 interface during large-scale conformational transitions to the CD4-bound state from the deformed states of gp120 in certain antibody complexes.

Spatial Normalization of Reverse Phase Protein Array Data

Reverse phase protein arrays (RPPA) are an efficient, high-throughput, cost-effective method for the quantification of specific proteins in complex biological samples. The quality of RPPA data may be affected by various sources of error. One of these, spatial variation, is caused by uneven exposure of different parts of an RPPA slide to the reagents used in protein detection. We present a method for the determination and correction of systematic spatial variation in RPPA slides using positive control spots printed on each slide. The method uses a simple bi-linear interpolation technique to obtain a surface representing the spatial variation occurring across the dimensions of a slide. This surface is used to calculate correction factors that can normalize the relative protein concentrations of the samples on each slide. The adoption of the method results in increased agreement between technical and biological replicates of various tumor and cell-line derived samples. Further, in data from a study of the melanoma cell-line SKMEL-133, several slides that had previously been rejected because they had a coefficient of variation (CV) greater than 15%, are rescued by reduction of CV below this threshold in each case. The method is implemented in the R statistical programing language. It is compatible with MicroVigene and SuperCurve, packages commonly used in RPPA data analysis. The method is made available, along with suggestions for implementation, at http://bitbucket.org/rppa_preprocess/rppa_preprocess/src.

Causal interactions from proteomic profiles: molecular data meets pathway knowledge

Measurement of changes in protein levels and in post-translational modifications, such as phosphorylation, can be highly informative about the phenotypic consequences of genetic differences or about the dynamics of cellular processes. Typically, such proteomic profiles are interpreted intuitively or by simple correlation analysis. Here, we present a computational method to generate causal explanations for proteomic profiles using prior mechanistic knowledge in the literature, as recorded in cellular pathway maps. To demonstrate its potential, we use this method to analyze the cascading events after EGF stimulation of a cell line, to discover new pathways in platelet activation, to identify influential regulators of oncoproteins in breast cancer, to describe signaling characteristics in predefined subtypes of ovarian and breast cancers, and to highlight which pathway relations are most frequently activated across 32 cancer types. Causal pathway analysis, that combines molecular profiles with prior biological knowledge captured in computational form, may become a powerful discovery tool as the amount and quality of cellular profiling rapidly expands. The method is freely available at http://causalpath.org.

The RNA-binding Protein MEX3B Mediates Resistance to Cancer Immunotherapy by Downregulating HLA-A Expression

Purpose: Cancer immunotherapy has shown promising clinical outcomes in many patients. However, some patients still fail to respond, and new strategies are needed to overcome resistance. The purpose of this study was to identify novel genes and understand the mechanisms that confer resistance to cancer immunotherapy. Experimental Design: To identify genes mediating resistance to T-cell killing, we performed an open reading frame (ORF) screen of a kinome library to study whether overexpression of a gene in patient-derived melanoma cells could inhibit their susceptibility to killing by autologous tumor-infiltrating lymphocytes (TIL). Results: The RNA-binding protein MEX3B was identified as a top candidate that decreased the susceptibility of melanoma cells to killing by TILs. Further analyses of anti–PD-1–treated melanoma patient tumor samples suggested that higher MEX3B expression is associated with resistance to PD-1 blockade. In addition, significantly decreased levels of IFNγ were secreted from TILs incubated with MEX3B-overexpressing tumor cells. Interestingly, this phenotype was rescued upon overexpression of exogenous HLA-A2. Consistent with this, we observed decreased HLA-A expression in MEX3B-overexpressing tumor cells. Finally, luciferase reporter assays and RNA-binding protein immunoprecipitation assays suggest that this is due to MEX3B binding to the 3′ untranslated region (UTR) of HLA-A to destabilize the mRNA. Conclusions: MEX3B mediates resistance to cancer immunotherapy by binding to the 3′ UTR of HLA-A to destabilize the HLA-A mRNA and thus downregulate HLA-A expression on the surface of tumor cells, thereby making the tumor cells unable to be recognized and killed by T cells. Clin Cancer Res; 24(14); 3366–76. ©2018 AACR. See related commentary by Kalbasi and Ribas, p. 3239