Microfluidics-based single-cell functional proteomics microchip for portraying protein signal transduction networks within the framework of physicochemical principles, with applications in fundamental and translational cancer research

Single-cell functional proteomics assays can connect genomic information to biological function through quantitative and multiplex protein measurements. Tools for single-cell proteomics have developed rapidly over the past 5 years and are providing unique opportunities. This thesis describes an emerging microfluidics-based toolkit for single cell functional proteomics, focusing on the development of the single cell barcode chips (SCBCs) with applications in fundamental and translational cancer research.

The microchip designed to simultaneously quantify a panel of secreted, cytoplasmic and membrane proteins from single cells will be discussed at the beginning, which is the prototype for subsequent proteomic microchips with more sophisticated design in preclinical cancer research or clinical applications. The SCBCs are a highly versatile and information rich tool for single-cell functional proteomics. They are based upon isolating individual cells, or defined number of cells, within microchambers, each of which is equipped with a large antibody microarray (the barcode), with between a few hundred to ten thousand microchambers included within a single microchip. Functional proteomics assays at single-cell resolution yield unique pieces of information that significantly shape the way of thinking on cancer research. An in-depth discussion about analysis and interpretation of the unique information such as functional protein fluctuations and protein-protein correlative interactions will follow.

The SCBC is a powerful tool to resolve the functional heterogeneity of cancer cells. It has the capacity to extract a comprehensive picture of the signal transduction network from single tumor cells and thus provides insight into the effect of targeted therapies on protein signaling networks. We will demonstrate this point through applying the SCBCs to investigate three isogenic cell lines of glioblastoma multiforme (GBM).

The cancer cell population is highly heterogeneous with high-amplitude fluctuation at the single cell level, which in turn grants the robustness of the entire population. The concept that a stable population existing in the presence of random fluctuations is reminiscent of many physical systems that are successfully understood using statistical physics. Thus, tools derived from that field can probably be applied to using fluctuations to determine the nature of signaling networks. In the second part of the thesis, we will focus on such a case to use thermodynamics-motivated principles to understand cancer cell hypoxia, where single cell proteomics assays coupled with a quantitative version of Le Chatelier's principle derived from statistical mechanics yield detailed and surprising predictions, which were found to be correct in both cell line and primary tumor model.

The third part of the thesis demonstrates the application of this technology in the preclinical cancer research to study the GBM cancer cell resistance to molecular targeted therapy. Physical approaches to anticipate therapy resistance and to identify effective therapy combinations will be discussed in detail. Our approach is based upon elucidating the signaling coordination within the phosphoprotein signaling pathways that are hyperactivated in human GBMs, and interrogating how that coordination responds to the perturbation of targeted inhibitor. Strongly coupled protein-protein interactions constitute most signaling cascades. A physical analogy of such a system is the strongly coupled atom-atom interactions in a crystal lattice. Similar to decomposing the atomic interactions into a series of independent normal vibrational modes, a simplified picture of signaling network coordination can also be achieved by diagonalizing protein-protein correlation or covariance matrices to decompose the pairwise correlative interactions into a set of distinct linear combinations of signaling proteins (i.e. independent signaling modes). By doing so, two independent signaling modes – one associated with mTOR signaling and a second associated with ERK/Src signaling have been resolved, which in turn allow us to anticipate resistance, and to design combination therapies that are effective, as well as identify those therapies and therapy combinations that will be ineffective. We validated our predictions in mouse tumor models and all predictions were borne out.

In the last part, some preliminary results about the clinical translation of single-cell proteomics chips will be presented. The successful demonstration of our work on human-derived xenografts provides the rationale to extend our current work into the clinic. It will enable us to interrogate GBM tumor samples in a way that could potentially yield a straightforward, rapid interpretation so that we can give therapeutic guidance to the attending physicians within a clinical relevant time scale. The technical challenges of the clinical translation will be presented and our solutions to address the challenges will be discussed as well. A clinical case study will then follow, where some preliminary data collected from a pediatric GBM patient bearing an EGFR amplified tumor will be presented to demonstrate the general protocol and the workflow of the proposed clinical studies.

Physical characterization of the rack effect and hydrogen bond networks in blue copper proteins

A summary of previous research is presented that indicates that the purpose of a blue copper protein's fold and hydrogen bond network, aka, the rack effect, enforce a copper(II) geometry around the copper(I) ion in the metal site. In several blue copper proteins, the C-terminal histidine ligand becomes protonated and detaches from the copper in the reduced forms. Mutants of amicyanin from Paracoccus denitrificans were made to alter the hydrogen bond network and quantify the rack effect by pK_a shifts.

The pK_a's of mutant amicyanins have been measured by pH-dependent electrochemistry. P94F and P94A mutations loosen the Northern loop, allowing the reduced copper to adopt a relaxed conformation: the ability to relax drives the reduction potentials up. The measured potentials are 265 (wild type), 380 (P94A), and 415 (P94F) mV vs. NHE. The measured pK_a's are 7.0 (wild type), 6.3 (P94A), and 5.0 (P94F). The additional hydrogen bond to the thiolate in the mutants is indicated by a red-shift in the blue copper absorption and an increase in the parallel hyperfine splitting in the EPR spectrum. This hydrogen bond is invoked as the cause for the increased stability of the C-terminal imidazole.

Melting curves give a measure of the thermal stability of the protein. A thermodynamic intermediate with pH-dependent reversibility is revealed. Comparisons with the electrochemistry and apoamicyanin suggest that the intermediate involves the region of the protein near the metal site. This region is destabilized in the P94F mutant; coupled with the evidence that the imidazole is stabilized under the same conditions confirms an original concept of the rack effect: a high energy configuration is stabilized at a cost to the rest of the protein.

Biochemical studies of postsynaptic density signaling proteins with a focus on synaptic GTPase activating protein and PDZ domains

Memory storage in the brain involves adjustment of the strength of existing synapses and formation of new neural networks. A key process underlying memory formation is synaptic plasticity, the ability of excitatory synapses to strengthen or weaken their connections in response to patterns of activity between their connected neurons. Synaptic plasticity is governed by the precise pattern of Ca²⁺ influx through postsynaptic N-methyl-D-aspartate-type glutamate receptors (NMDARs), which can lead to the activation of the small GTPases Ras and Rap. Differential activation of Ras and Rap acts to modulate synaptic strength by promoting the insertion or removal of 2-amino-3-(3-hydroxy-5-methyl-isoxazol-4-yl)propanoic acid receptors (AMPARs) from the synapse. Synaptic GTPase activating protein (synGAP) regulates AMPAR levels by catalyzing the inactivation of GTP-bound (active) Ras or Rap. synGAP is positioned in close proximity to the cytoplasmic tail regions of the NMDAR through its association with the PDZ domains of PSD-95. SynGAP’s activity is regulated by the prominent postsynaptic protein kinase, Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) and cyclin-dependent kinase 5 (CDK5), a known binding partner of CaMKII. Modulation of synGAP’s activity by phosphorylation may alter the ratio of active Ras to Rap in spines, thus pushing the spine towards the insertion or removal of AMPARs, subsequently strengthening or weakening the synapse. To date, all biochemical studies of the regulation of synGAP activity by protein kinases have utilized impure preparations of membrane bound synGAP. Here we have clarified the effects of phosphorylation of synGAP on its Ras and Rap GAP activities by preparing and utilizing purified, soluble recombinant synGAP, Ras, Rap, CaMKII, CDK5, PLK2, and CaM. Using mass spectrometry, we have confirmed the presence of previously identified CaMKII and CDK5 sites in synGAP, and have identified novel sites of phosphorylation by CaMKII, CDK5, and PLK2. We have shown that the net effect of phosphorylation of synGAP by CaMKII, CDK5, and PLK2 is an increase in its GAP activity toward HRas and Rap1. In contrast, there is no effect on its GAP activity toward Rap2. Additionally, by assaying the GAP activity of phosphomimetic synGAP mutants, we have been able to hypothesize the effects of CDK5 phosphorylation at specific sites in synGAP. In the course of this work, we also found, unexpectedly, that synGAP is itself a Ca²⁺/CaM binding protein. While Ca²⁺/CaM binding does not directly affect synGAP activity, it causes a conformational change in synGAP that increases the rate of its phosphorylation and exposes additional phosphorylation sites that are inaccessible in the absence of Ca²⁺/CaM.

The postsynaptic density (PSD) is an electron-dense region in excitatory postsynaptic neurons that contains a high concentration of glutamate receptors, cytoskeletal proteins, and associated signaling enzymes. Within the PSD, three major classes of scaffolding molecules function to organize signaling enzymes and glutamate receptors. PDZ domains present in the Shank and PSD-95 scaffolds families serve to physically link AMPARs and NMDARs to signaling molecules in the PSD. Because of the specificity and high affinity of PDZ domains for their ligands, I reasoned that these interacting pairs could provide the core components of an affinity chromatography system, including affinity resins, affinity tags, and elution agents. I show that affinity columns containing the PDZ domains of PSD-95 can be used to purify active PDZ domain-binding proteins to very high purity in a single step. Five heterologously expressed neuronal proteins containing endogenous PDZ domain ligands (NMDAR GluN2B subunit Tail, synGAP, neuronal nitric oxide synthase PDZ domain, cysteine rich interactor of PDZ three and cypin) were purified using PDZ domain resin, with synthetic peptides having the sequences of cognate PDZ domain ligands used as elution agents. I also show that conjugation of PDZ domain-related affinity tags to Proteins Of Interest (POIs) that do not contain endogenous PDZ domains or ligands does not alter protein activity and enables purification of the POIs on PDZ domain-related affinity resins.