NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

Protein-protein interaction between phototaxis receptor sensory Rhodopsin I and its transducer HtrI

The molecular complex containing the seven transmembrane helix photoreceptor S&barbelow;ensory R&barbelow;hodopsin I&barbelow; (SRI) and transducer protein HtrI (H&barbelow;alobacterial Transducer for SRI&barbelow;) mediates color-sensitive phototaxis responses in the archaeon Halobacterium salinarum. Orange light causes an attractant response by a one-photon reaction and white light (orange + UV light) a repellent response by a two-photon reaction. Three aspects of SRI-HtrI structure/function and the signal transduction pathway were explored. First, the coupling of HtrI to the photoactive site of SRI was analyzed by mutagenesis and kinetic spectroscopy. Second, SRI-HtrI mutations and suppressors were selected and characterized to elucidate the color-sensing mechanism. Third, the signal relay through the transducer-bound histidine kinase was analyzed using an in vitro reconstitution system with known and newly identified taxis components. ^ Twenty-one mutations on HtrI were introduced by site-directed mutagenesis. Several replacements of charged residues perturbed the photochemical kinetics of SRI which led to the finding of a cluster of residues at the membrane/cytoplasm interface in HtrI electrostatically coupled to the photoactive site of SRI. We found by laser-flash kinetic spectroscopy that the transducer and these residues have specific effects on the light-induced proton transfer between the retinal chromophore and the protein. ^ One of the mutations showed an unusual mutant phenotype we called “inverted” signaling, in which the cell produces a repellent response to normally attractant light. Therefore, this mutant (E56Q of HtrI) had lost the color-discrimination by the SRI-HtrI complex. We used suppressor analysis to better understand the phenotype. Certain suppressors resulted in return of attractant responses to orange light but with inversion of the normally repellent response to white light to an attractant response. To explain this and other results, we formulated the Conformational Shuttling model in which the HtrI-SRI complex is poised in a metastable equilibrium of two conformations shifted in opposite directions by orange and white light. We tested this model by behavioral analysis (computerized cell tracking and motion study) of double mutants of inverting and suppressing mutations and the results confirmed the equilibrium-shift explanation. ^ We developed an in vitro system for measuring the effect of purified transducer on the histidine-kinase CheAH that controls the flagellar motor switch. The rate of kinase autophosphorylation was stimulated >2 fold in the reconstitution of the complete signal transduction system from purified components from H. salinarum. The in vitro assay also showed that the kinase activity was reduced in the absence and in the presence of high levels of linker protein CheWH. (Abstract shortened by UMI.) ^

Inhibitor of cytochrome c messenger RNA -protein interaction in stimulated rat skeletal muscle

The purpose of this work was to examine the possible mechanisms for the regulation of cytochrome c gene expression in response to increased contractile activity in rat skeletal muscle. The working hypothesis was that increased contractile activity enhances cytochrome c gene expression through a cis-element. A 110% increase in cytochrome c mRNA concentration was observed in tibialis anterior (TA) muscle after 9 days of chronic stimulation. Similar difference (120%) exists between soleus (SO) muscle of higher contractile activity and white vastus lateralis (WV) muscle of lower contractile activity. These results suggest that the endogenous cytochrome c gene expression is regulated by contractile activity. Cytochrome c-reporter genes were injected into skeletal muscles to identify the cis-element that is responsible for the regulation. Although the data was inconclusive, part of it suggested the importance of the 3$\sp\prime$-untranslated region (3$\sp\prime$-UTR) in mediating the response to increased contractile activity.^ RNA gel mobility shift (GMSA) and ultraviolet (UV) cross-linking assays revealed specific RNA-protein interaction in a 50-nucleotide region of the 3$\sp\prime$-UTR in unstimulated TA muscle. Computer analysis predicted a stem-loop structure of 17 nucleotides, which provides a structural basis for RNA-protein interaction. These 17 nucleotides are 100% conserved among rat, mouse and human cytochrome c genes and their 13 pseudogenes, suggesting a functional role for this region. The RNA-protein interaction was significantly less in highly active SO muscle than in inactive WV muscle and was dramatically decreased in stimulated TA muscle due to a protein inhibitor(s) associated with ribosome. It is possible that cytochrome c mRNAs undergoing translation are subject to a compartmentalized regulatory influence.^ The conclusion from these results is that increases in contractile activity induce or activate a protein inhibitor(s) associated with ribosome in rat skeletal muscle. The inhibitor decreases RNA-protein interaction in the 3$\sp\prime$-UTR of cytochrome c mRNA, which may result in increased mRNA stability and/or translation. ^

Physical and functional interaction between the dopamine transporter and the synaptic vesicle protein synaptogyrin-3.

Uptake through the dopamine transporter (DAT) represents the primary mechanism used to terminate dopaminergic transmission in brain. Although it is well known that dopamine (DA) taken up by the transporter is used to replenish synaptic vesicle stores for subsequent release, the molecular details of this mechanism are not completely understood. Here, we identified the synaptic vesicle protein synaptogyrin-3 as a DAT interacting protein using the split ubiquitin system. This interaction was confirmed through coimmunoprecipitation experiments using heterologous cell lines and mouse brain. DAT and synaptogyrin-3 colocalized at presynaptic terminals from mouse striatum. Using fluorescence resonance energy transfer microscopy, we show that both proteins interact in live neurons. Pull-down assays with GST (glutathione S-transferase) proteins revealed that the cytoplasmic N termini of both DAT and synaptogyrin-3 are sufficient for this interaction. Furthermore, the N terminus of DAT is capable of binding purified synaptic vesicles from brain tissue. Functional assays revealed that synaptogyrin-3 expression correlated with DAT activity in PC12 and MN9D cells, but not in the non-neuronal HEK-293 cells. These changes were not attributed to changes in transporter cell surface levels or to direct effect of the protein-protein interaction. Instead, the synaptogyrin-3 effect on DAT activity was abolished in the presence of the vesicular monoamine transporter-2 (VMAT2) inhibitor reserpine, suggesting a dependence on the vesicular DA storage system. Finally, we provide evidence for a biochemical complex involving DAT, synaptogyrin-3, and VMAT2. Collectively, our data identify a novel interaction between DAT and synaptogyrin-3 and suggest a physical and functional link between DAT and the vesicular DA system.

The roles of electrostatic charge pairing in the interaction between cytochrome P-450 and NADPH-cytochrome P-450 reductase

The cytochrome P450 enzyme catalysis requires two electrons transferred from NADPH-cytochrome P450 reductase (reductase) to P450. Electrostatic charge-pairing has been proposed to be one of the major forces in the interaction between P450 and reductase. In order to obtain further insight into the molecular basis for the protein interaction, I used two methods, chemical modification and specific anti-peptide antibodies, to study the involvement and importance of charged amino acid residues. Acetylation of lysine residues of P450c and P450b by acetic anhydride dramatically inhibited the reductase-supported P450c-dependent ethoxycoumarin hydroxylation activity, but P450 activity supported by cumene hydroperoxide is relatively unchanged. The modification of lysine residues of P450c and P450b did not grossly disturb the protein conformation as revealed by several spectral studies. This differential effect of lysine modification on the P450 activity in the system reconstituted with reductase versus the system supported by cumene hydroperoxide suggested an important role for P450 lysine residues in the interaction with reductase. Using $\rm\sp{14}C$-acetic anhydride, P450 lysine residues were labelled and further identified on P450c and P450b. Those lysine residues are at position 97, 271, 279, and 407 for P450c, and 251, 384, 422, 433, and 473 for P450b. Alignment of those identified lysine residues on P450c and P450b with amino acid residues identified in other studies indicated those residues reside in three major sequence areas. Modification of arginine residues of P450b by phenylglyoxal and 2, 3-butanedione have no significant effect on P450 activity either supported by NADPH and reductase or supported by cumene hydroperoxide. Further studies using $\rm\sp{14}C$-phenylglyoxal reveals that no incorporation of phenylglyoxal into P450b was found. These results demonstrated a predominant role of lysine residues of P450 in the electrostatic interaction with reductase. To understand the protein binding sites on each of P450 and reductase, I generated three anti-peptide antibodies against regions on reductase and five anti-peptide antibodies against five putative reductase binding sites on P450c. These anti-peptide antibodies were affinity purified and characterized on ELISA and by Western blot analysis. Inhibition experiments using these antibodies demonstrated that regions 109-120 and 204-220 of reductase are probably the two major binding sites for P450. The association of reductase with cytochromes P450 and cytochrome c may rely on different mechanisms. The data from experiments using anti-peptide (P450c) antibodies supports the important role of P450c lysine residues 271/279 and 458/460 in the interaction with reductase. ^

Molecular mechanisms of prostacyclin receptor signaling: The intracellular domains in G protein coupling

To better understand the mechanisms of how the human prostacyclin receptor (1P) mediates vasodilation and platelet anti-aggregation through Gs protein coupling, a strategy integrating multiple approaches including high resolution NMR experiments, synthetic peptide, fluorescence spectroscopy, molecular modeling, and recombinant protein was developed and used to characterize the structure/function relationship of important segments and residues of the IP receptor and the α-subunit of the Gs protein (Gαs). The first (iLP1) and third (iLP3) intracellular loops of the IP receptor, as well as the Gαs C-terminal domain, relevant to the Gs-mediated IP receptor signaling, were first identified by observation of the effects of the mini gene-expressed corresponding protein segments in HEK293 cells which co-expressed the receptor and Gαs. Evidence of the IP iLP1 domain interacted with the Gαs C-terminal domain was observed by fluorescence and NMR spectroscopic studies using a constrained synthetic peptide, which mimicked the IP iLP1 domain, and the synthetic peptide, which mimicked Gαs C-terminal domain. The solution structural models and the peptide-peptide interaction of the two synthetic protein segments were determined by high resolution NMR spectroscopy. The important residues in the corresponding domains of the IP receptor and the Gαs predicted by NMR chemical shift mapping were used to guide the identification of their protein-protein interaction in cells. A profile of the residues Arg42 - Ala48 of the IP iLP1 domain and the three residues Glu392 ∼ Leu394 of the Gαs C-terminal domain involved in the IP/Gs protein coupling were confirmed by recombinant proteins. The data revealed an intriguing speculation on the mechanisms of how the signal of the ligand-activated IP receptor is transmitted to the Gs protein in regulating vascular functions and homeostasis, and also provided substantial insights into other prostanoid receptor signaling. ^

Macromolecular interaction studies on FF1-3 of human CA150 and STAT1-importin alpha5 using biophysical and biochemical methods

Macromolecular interactions, such as protein-protein interactions and protein-DNA interactions, play important roles in executing biological functions in cells. However the complexity of such interactions often makes it very challenging to elucidate the structural details of these subjects. In this thesis, two different research strategies were applied on two different two macromolecular systems: X-ray crystallography on three tandem FF domains of transcription regulator CA150 and electron microscopy on STAT1-importin α5 complex. The results from these studies provide novel insights into the function-structure relationships of transcription coupled RNA splicing mediated by CA150 and the nuclear import process of the JAK-STAT signaling pathway. ^ The first project aimed at the protein-protein interaction module FF domain, which often occurs as tandem repeats. Crystallographic structure of the first three FF domains of human CA150 was determined to 2.7 Å resolution. This is the only crystal structure of an FF domain and the only structure on tandem FF domains to date. It revealed a striking connectivity between an FF domain and the next. Peptide binding assay with the potential binding ligand of FF domains was performed using fluorescence polarization. Furthermore, for the first time, FF domains were found to potentially interact with DNA. DNA binding assays were also performed and the results were supportive to this newly proposed functionality of an FF domain. ^ The second project aimed at understanding the molecular mechanism of the nuclear import process of transcription factor STAT1. The first structural model of pSTAT1-importin α5 complex in solution was built from the images of negative staining electron microscopy. Two STAT1 molecules were observed to interact with one molecule of importin α5 in an asymmetric manner. This seems to imply that STAT1 interacts with importin α5 with a novel mechanism that is different from canonical importin α-cargo interactions. Further in vitro binding assays were performed to obtain more details on the pSTAT1-importin α5 interaction. ^

Active water in protein-protein communication within the membrane: the case of SRII-HtrII signal relay.

We detect internal water molecules in a membrane-embedded receptor-transducer complex and demonstrate water structure changes during formation of the signaling state. Time-resolved FTIR spectroscopy reveals stimulus-induced repositioning of one or more structurally active water molecules to a significantly more hydrophobic environment in the signaling state of the sensory rhodopsin II (SRII)-transducer (HtrII) complex. These waters, distinct from bound water molecules within the SRII receptor, appear to be in the middle of the transmembrane interface region near the Tyr199(SRII)-Asn74(HtrII) hydrogen bond. We conclude that water potentially plays an important role in the SRII --> HtrII signal transfer mechanism in the membrane's hydrophobic core.

CCTeta: A novel soluble guanylyl cyclase-interacting protein

Nitric oxide (NO) transduces most of its biological effects through activation of the heterodimeric enzyme, soluble guanylyl cyclase (sGC). Activation of sGC results in the production of 3′,5 ′-cyclic guanosine monophosphate (cGMP) from 5′ -guanosine triphosphate (GTP). In this thesis, we demonstrate a novel protein interaction between CCT (chaperonin containing t-complex polypeptide) subunit η and the α1β1 isoform of sGC. Using the yeast-two-hybrid system, CCTη was found to interact with the N-terminal portion of β1 subunit of sGC. This interaction was then confirmed in vitro with a co-immunoprecipitation from mouse brain. The interaction between these two proteins was further supported by a co-localization of the proteins within rat brain. Using the yeast-two-hybrid system, CCTη was found to bind to the N-terminal portion of sGC. In vitro assays with purified CCTη and Sf9 lysate expressing sGC resulted in a 33% inhibition of sodium nitroprusside (SNP)-stimulated sGC activity. The same assays were then performed using BAY41-2272, an NO-independent allosteric sGC activator, and CCTη had no effect on this activity. Furthermore, CCTη had no effect on the activity of αβCys105 sGC a constitutively active mutant that lacks a heme group. Of note is the fact that the full-length CCTη-expressing bacterial lysate inhibited the activity of sGC-expressing Sf9 lysate by 48% compared with GST alone. This indicates that the amino terminal 94 amino acids of CCTη are important to the inhibition of sGC activity. Lastly, a 45% inhibition of sGC activity by CCTη was seen in vivo in BE2 cells stably transfected with CCTη and treated with SNP. The fact that the inhibition of sGC was more pronounced with bacterial lysate expressing CCTη versus the purified CCTη implies that some factor in the bacterial lysate enhances the inhibitory effect of CCTη. Because the level of inhibition seen in bacterial lysate and in vivo experiments is similar, might imply that the factor that aids in CCTη effect on sGC is conserved. Together, these data suggest that CCTη is a novel type of sGC inhibitor that inhibits sGC by modifying the binding of NO to the heme group or the subsequent conformational changes induced by NO binding. ^

Tackling ErbB2: ErbB2-targeting peptide therapy and combination therapy with trastuzumab plus PI3K inhibitors

ErbB2 is an excellent target for cancer therapies because its overexpression was found in about 30% of breast cancers and correlated with poor prognosis of the patients. Unfortunately, current therapies for ErbB2-positive breast cancers remain unsatisfying due to side effects and resistance, and new therapies for ErbB2 overexpressing breast cancers are needed. Peptide/protein therapy using cell-penetrating peptides (CPPs) as carriers is promising because the internalization is highly efficient and the cargos can be bioactive. The major obstacle in using CPPs for therapy is their lack of specificity. We sought to develop a peptide carrier specifically introducing therapeutics to ErbB2-overexpressing breast cancer cells. By modifying the TAT-derived CPP, and attaching anti-HER2/neu peptide mimetic (AHNP), we developed the peptide carrier (P3-AHNP) specifically targeted ErbB2-overexpressing breast cancers in vitro and in vivo. A STAT3 SH2 domain-binding peptide conjugated to this peptide carrier (P3-AHNP-STAT3BP) was delivered preferentially into ErbB2-overexpressing breast cancer cells in vitro and in vivo. P3-AHNP-STAT3BP inhibited growth and induced apoptosis in vitro, with ErbB2-overexpressing 435.eB cells being more sensitive than the ErbB2-lowexpressing MDA-MB-435 cells. P3-AHNP-STAT3BP preferentially accumulated and inhibited growth in 435.eB xenografts, comparing with MDA-MB-435 xenografts or normal tissues with low levels of ErbB2. This ErbB2-targeting peptide delivery system provided the basis for future development of novel cancer target-specific treatments with low toxicity to normal cells. ^ Another urgent issue in treating ErbB2-positive breast cancers is trastuzumab resistance. Trastuzumab is the only FDA-approved ErbB2-targeting antibody for treatment of metastatic breast cancers overexpressing ErbB2, and has remarkable therapeutic efficacy in certain patients. The overall trastuzumab response rate, however, is limited, and understanding the mechanisms of trastuzumab resistance is needed to overcome this problem. We report that PTEN activation contributes to trastuzumab's anti-tumor activity. Trastuzumab treatment quickly inactivated Src, which reduced PTEN tyrosine phosphorylation, increased PTEN membrane localization and its phosphatase activity in cancer cells. Reducing PTEN expression in breast cancer cells by antisense oligonucleotides conferred trastuzumab resistance in vitro and in vivo. Importantly, PI3K inhibitors sensitized PTEN-deficient breast cancers to the growth inhibition by trastuzumab in vitro and in vivo, suggesting that combination therapies with PI3K inhibitors plus trastuzumab could overcome trastuzumab resistance. ^

Xp95 is phosphorylated within the proline rich domain during Xenopus oocyte maturation

Xp95 is the Xenopus ortholog of a conserved family of scaffold proteins that have in common an N-terminal Bro1 domain and a C-terminal proline rich domain (PRD). The regulation of this protein family is poorly understood. We previously showed that Xp95 undergoes a phosphorylation-dependant gel mobility shift during meiotic maturation of Xenopus oocytes, the only natural biological system in which post-translational modifications of this family has been demonstrated. Here we characterized Xp95 phosphorylation via two approaches. First, we tested a series of Xp95 fragments for the ability to gel-shift during oocyte maturation, and found that a fragment containing amino acids 705-786 is sufficient to cause a gel-shift. This fragment is within the N-terminal region of Xp95's PRD (N-PRD). Second, we purified phosphorylated Xp95 and by mass spectrometry found that a 5080 Da peptide which maps to N-PRD (amino acids 706-756) contains two phosphorylation sites, one of which is T745, within the conserved CIN85 binding motif. By in vitro protein interaction assays, we that T745 is critical for CIN85/Xp95 interaction, and that Xp95 phosphorylation correlates with loss of binding to CIN85. We also show that an Alix fragment (amino acids 604-789) also undergoes a gel-shift during oocyte maturation and during colcemid-induced mitotic arrest of HeLa cells. These findings indicate that Xp95/Alix is phosphorylated on the PRD during M phase induction and that the PRD phosphorylation regulates partner protein interaction. ^

A mechanism of insulin-like growth factor binding protein 2-induced cell mobility in glioma

Insulin-like growth factor binding protein 2 (IGFBP2) is a protein known to be overexpressed in a majority of glioblastoma multiforme (GBM) tumors. While it is known the IGFBP2 is involved in promoting GBM tumor cell invasion, no mechanism exists for how the protein is involved in signal transduction pathways leading to enhanced cell invasion. ^ We follow up on preliminary microarray data on IGFBP2-overexpressing GBM cells and protein sequence analysis of IGFBP2 in generating the hypothesis that IGFBP2 interacts with integnn α5 in regulating cell mobility. Microarray data showing upregulation of integrin α5 by IGFBP2 is validated and evidence of protein-protein interaction between IGFBP2 and integrin α5 is found. The exact binding domain on IGFBP2 responsible for its interaction with integrin α5 is also determined, confirming our initial findings and reaffirming that the IGFBP2/integrin α5 interaction is specific. Disruption of this interaction resulted in attenuation of IGFBP2-enhanced cell mobility. Further, we found that cell mobility is only enhanced when IGFBP2 and integrin α5 are both overexpressed and able to interact with each other. ^ We also determined fibronectin to be a critical player in the activation of the IGFBP2/integrin α5 pathway. The activation of this pathway appears to be progressive and initiates once GBM cells have sufficiently established anchorage. ^

The role of AKT-mediated phosphorylation in suppression of MAD1 and the development of new approach in protein complex detection

MAX dimerization protein 1 (MAD1) is a basic-helix-loop-helix transcription factors that recruits transcription repressor such as HDAC to suppress target genes transcription. It antagonizes to MYC because the promoter binding sites for MYC are usually also serve as the binding sites for MAD1 so they compete for it. However, the mechanism of the switch between MYC and MAD1 in turning on and off of genes' transcription is obscure. In this study, we demonstrated that AKT-mediated MAD1 phosphorylation inhibits MAD1 transcription repression function. The association between MAD1 and its target genes' promoter is reduced after been phosphorylated by AKT; therefore, consequently, allows MYC to occupy the binding site and activates transcription. Mutation of such phosphorylation site abrogates the inhibition from AKT. In addition, functional assays demonstrated that AKT suppressed MAD1-mediated transcription repression of its target genes hTERT and ODC. Cell cycle and cell growth were also been released from inhibition by MAD1 in the presents of AKT. Taken together, our study suggests that MAD1 is a novel substrate of AKT and AKT-mediated MAD1 phosphorylation inhibits MAD1function; therefore, activates MAD1 target genes expression. ^ Furthermore, analysis of protein-protein interaction is indispensable for current molecular biology research, but multiplex protein dynamics in cells is too complicated to be analyzed by using existing biochemical methods. To overcome the disadvantage, we have developed a single molecule level detection system with nanofluidic chip. Single molecule was analyzed based on their fluorescent profile and their profiles were plotted into 2 dimensional time co-incident photon burst diagram (2DTP). From this 2DTP, protein complexes were characterized. These results demonstrate that the nanochannel protein detection system is a promising tool for future molecular biology. ^

Identification of a conserved cluster in the RH domain of GRK critical for activation by GPCRs

One of the most critical aspects of G Protein Coupled Receptors (GPCRs) regulation is their rapid and acute desensitization following agonist stimulation. Phosphorylation of these receptors by GPCR kinases (GRK) is a major mechanism of desensitization. Considerable evidence from studies of rhodopsin kinase and GRK2 suggests there is an allosteric docking site for the receptor distinct from the GRK catalytic site. While the agonist-activated GPCR appears crucial for GRK activation, the molecular details of this interaction remain unclear. Recent studies suggested an important role for the N- and C-termini and domains in the small lobe of the kinase domain in allosteric activation; however, neither the mechanism of action of that site nor the RH domain contributions have been elucidated. To search for the allosteric site, we first indentified evolutionarily conserved sites within the RH and kinase domains presumably deterministic of protein function employing evolutionary trace (ET) methodology and crystal structures of GRK6. Focusing on a conserved cluster centered on helices 3, 9, and 10 in the RH domain, key residues of GRK5 and 6 were targeted for mutagenesis and functional assays. We found that a number of double mutations within helices 3, 9, and 10 and the N-terminus markedly reduced (50–90%) the constitutive phosphorylation of the β-2 Adrenergic Receptor (β2AR) in intact cells and phosphorylation of light-activated rhodopsin (Rho*) in vitro as compared to wild type (WT) GRK5 or 6. Based on these results, we designed peptide mimetics of GRK5 helix 9 both computationally and through chemical modifications with the goal of both confirming the importance of helix 9 and developing a useful inhibitor to disrupt the GPCR-GRK interaction. Several peptides were found to block Rho* phosphorylation by GRK5 including the native helix 9 sequence, Peptide Builder designed-peptide preserving only the key ET residues, and chemically locked helices. Most peptidomimetics showed inhibition of GRK5 activity greater than 80 % with an IC50 of ∼ 30 µM. Alanine scanning of helix 9 has further revealed both essential and non-essential residues for inhibition. Importantly, substitution of Arg 169 by an alanine in the native helix 9-based peptide gave an almost complete inhibition at 30 µM with an IC50 of ∼ 10 µM. In summary we report a previously unrecognized crucial role for the RH domain of GRK5 and 6, and the subsequent identification of a lead peptide inhibitor of protein-protein interaction with potential for specific blockade of GPCR desensitization. ^

Targeting TRAF6 for Cancer Therapeutical Development

Tumor necrosis factor (TNF)-Receptor Associated Factors (TRAFs) are a family of signal transducer proteins. TRAF6 is a unique member of this family in that it is involved in not only the TNF superfamily, but the toll-like receptor (TLR)/IL-1R (TIR) superfamily. The formation of the complex consisting of Receptor Activator of Nuclear Factor κ B (RANK), with its ligand (RANKL) results in the recruitment of TRAF6, which activates NF-κB, JNK and MAP kinase pathways. TRAF6 is critical in signaling with leading to release of various growth factors in bone, and promotes osteoclastogenesis. TRAF6 has also been implicated as an oncogene in lung cancer and as a target in multiple myeloma. In the hopes of developing small molecule inhibitors of the TRAF6-RANK interaction, multiple steps were carried out. Computational prediction of hot spot residues on the protein-protein interaction of TRAF6 and RANK were examined. Three methods were used: Robetta, KFC2, and HotPoint, each of which uses a different methodology to determine if a residue is a hot spot. These hot spot predictions were considered the basis for resolving the binding site for in silico high-throughput screening using GOLD and the MyriaScreen database of drug/lead-like compounds. Computationally intensive molecular dynamics simulations highlighted the binding mechanism and TRAF6 structural changes upon hit binding. Compounds identified as hits were verified using a GST-pull down assay, comparing inhibition to a RANK decoy peptide. Since many drugs fail due to lack of efficacy and toxicity, predictive models for the evaluation of the LD50 and bioavailability of our TRAF6 hits, and these models can be used towards other drugs and small molecule therapeutics as well. Datasets of compounds and their corresponding bioavailability and LD50 values were curated based, and QSAR models were built using molecular descriptors of these compounds using the k-nearest neighbor (k-NN) method, and quality of these models were cross-validated.