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.

CONFORMATIONAL DYNAMICS OF K-RAS AND H-RAS PROTEINS: IS THERE FUNCTIONAL SPECIFICITY AT THE CATALYTIC DOMAIN?

Ras proteins serve as crucial signaling modulators in cell proliferation through their ability to hydrolyze GTP and exist in a GTP “on” state and GTP “off” state. There are three different human Ras isoforms: H-ras, N-ras and K-ras (4A and 4B). Although their sequence identity is very high at the catalytic domain, these isoforms differ in their ability to activate different effectors and hence different signaling pathways. Much of the previous work on this topic has attributed this difference to the hyper variable region of Ras proteins, which contains most of the sequence variance among the isoforms and encodes specificity for differential distribution in the membrane. However, we hypothesize that sequence variation on lobe II of Ras catalytic domain alters dynamics and leads to differential preference for different effectors or modulators. In this work, we used all atom molecular dynamics to analyze the dynamics in the catalytic domain of H-ras and K-ras. We have also analyzed the dynamics of a transforming mutant of H-ras and K-ras and further studied the dynamics of an effectorselective mutant of H-ras. Collectively we have determined that wild type K-ras is more dynamic than H-ras and that the structure of the effector binding loop more closely resembles that of the T35S Raf-selective mutant, possibly giving us a new view and insight into the v mode of effector specificity. Furthermore we have determined that specific mutations at the same location perturb the conformational equilibrium differently in H-ras and K-ras and that an enhanced oncogenic potential may arise from different structural perturbations for each point mutation of a specific isoform.